Nanopore-Based Surveillance of Leishmania Parasites in Culicoides Latrielle (Diptera: Ceratopogonidae) Caught from the Affected Community and Tham Phra Cave in Chiang Rai Province, the Endemic Area of Leishmaniasis in Northern Thailand.

Simple Summary: The number of leishmaniasis cases caused by autochthonous Leishmania (Mundinia) martiniquensis and Leishmania (Mundinia) orientalis has been continuously growing in Southeast Asia over the decades, particularly in Thailand. Recent evidence has suggested that Culicoides biting midges are the most likely natural vectors of autochthonous leishmaniasis in Thailand. Nevertheless, the epidemiology of vector infection and the genetic variation of Leishmania parasites in endemic areas have remained largely unknown. In this study, we report the high infection prevalence, sympatric circulation, and regional genetic diversity of two Leishmania (Mundinia) parasites in Culicoides spp. caught from the affected community and Tham Phra cave in Chiang Rai Province, Northern Thailand. These results provide us with a more complete understanding of the complexity of leishmaniasis transmission, which will be useful for establishing effective management and control measures for this neglected disease, especially in endemic areas of Northern Thailand. In this research, we elucidated the species composition of Culicoides biting midges, infection prevalence, and genetic diversity of Leishmania parasites circulating in the affected community in Chiang Rai Province, being the most endemic area in Northern Thailand. A total of 146 parous and gravid females, belonging to at least twelve Culicoides species in five subgenera and one species group, were trapped from three collection sites with an overall Leishmania prevalence of 26.7% (39/146). Leishmania was detected, using ITS1-PCR, in C. mahasarakamense (15), C. guttifer (11), C. (Trithecoides) spp. (8), C. jacobsoni (2), C. oxystoma (2), and C. orientalis (1). The evidence of Leishmania infection in these last five species represents new records in Northern Thailand. Given a high infection rate in cavernicolous specimens, this indicates an increased risk of parasite exposure when visiting the cave. Using the nanopore amplicon sequencing, L. martiniquensis was ubiquitously identified in all positives, and more than half of these were also co-infected with L. orientalis. The genetic diversity analysis revealed 13 and 17 unique haplotypes for L. martiniquensis and L. orientalis, respectively. Higher haplotype diversity and relatively low nucleotide diversity were observed in both parasite populations, suggesting recent population divergence. Neutrality tests (Tajima's D and Fu and Li's D) showed to be significantly negative, indicating rapid population growth or a selective sweep. Moreover, dominant haplotypes of both Leishmania species were 100% identical to those in all leishmaniasis patients previously reported from Northern Thailand, strongly supporting the imperative role of Culicoides spp. in disease transmission. Essentially, this research provides the first entomological surveillance data representing the sympatric existence, transmission dynamics, and genetic complexity of two autochthonous Leishmania (Mundinia) parasites in several Culicoides species in the endemic area of Northern Thailand. This would contribute to a more complete understanding of the epidemiology of vector infection and facilitate the development of vector control programs to effectively reduce the transmission of this neglected tropical disease in endemic areas of Northern Thailand. [ABSTRACT FROM AUTHOR]