Molecular Identification of Host Blood Meals and Detection of Blood Parasites in Culicoides Latreille (Diptera: Ceratopogonidae) Collected from Phatthalung Province, Southern Thailand

Seasonal dynamics, Leishmania diversity, and nanopore-based metabarcoding of blood meal origins in Culicoides spp. in the newly emerging focus of leishmaniasis in Northern Thailand

BACKGROUND

Clinical cases of leishmaniasis caused by Leishmania (Mundinia) parasites have been increasingly reported in Southeast Asia, particularly Thailand. Recent evidence has shown that Leishmania (Mundinia) parasites successfully developed into infective metacyclic promastigotes in Culicoides biting midges, strongly supporting their putative role in disease transmission. However, Culicoides diversity, host preference, and Leishmania prevalence in endemic areas remain largely unknown.

METHODS

We investigated the seasonal dynamics, infection prevalence, and blood meal identification of Culicoides collected from the emerging focus of visceral leishmaniasis in Lampang Province, Northern Thailand, during 2021-2023. Midge samples were molecularly screened for Leishmania using SSU rRNA-qPCR and ITS1-PCR, followed by Sanger plasmid sequencing, and parasite haplotype diversity was analyzed. Host blood meal origins were comparatively identified using host-specific Cytb-PCRs and a nanopore-based metabarcoding approach.

RESULTS

A total of 501 parous and gravid females and 46 blood-engorged ones belonging to at least 17 species of five subgenera (Remmia, Trithecoides, Avaritia, Hoffmania, and Meijerehelea) and two species groups (Shortti and Calvipalpis) were collected with temporal differences in abundance. Leishmania was detected by SSU rRNA-qPCR in 31 samples of at least 11 midge species, consisting of Culicoides oxystoma, C. guttifer, C. orientalis, C. mahasarakhamense, C (Trithecoides) spp., C. innoxius, C. shortti, C. arakawae, C. sumatrae, C. actoni, and C. fulvus, with the overall infection prevalence of 5.7%. The latter six species represent the new records as putative leishmaniasis vectors in Northern Thailand. The ITS1-PCR and plasmid sequencing revealed that Leishmania martiniquensis was predominantly identified in all qPCR-positive species, whereas L. orientalis was identified only in three C. oxystoma samples. The most dominant haplotype of L. martiniquensis in Thailand was genetically intermixed with those from other geographical regions, confirming its globalization. Neutrality test statistics were also significantly negative on regional and country-wide scales, suggesting rapid population expansion or selective sweeps. Nanopore-based blood meal analysis revealed that most Culicoides species are mammalophilic, with peridomestic and wild mammals (cow, pig, deer, and goat-like species) and humans as hosts, while C. guttifer and C. mahasarakhamense fed preferentially on chickens.

CONCLUSIONS

This study revealed seasonal dynamics and sympatric circulation of L. martiniquensis and L. orientalis in different species of Culicoides. Evidence of human blood feeding was also demonstrated, implicating Culicoides as putative vectors of human leishmaniasis in endemic areas. Further research is therefore urgently needed to develop vector control strategies and assess the infection status of their reservoir hosts to effectively minimize disease transmission.

Seasonal Abundance and Diversity of Culicoides Biting Midges in Livestock Sheds in Kanchanaburi Province, Thailand

Culicoides biting midges were collected using UV light traps from nine livestock farms in Kanchanaburi Province, Thailand. Collections were made one night per month from June 2020 to May 2021 to determine the seasonal changes and spatial distribution of the Culicoides assemblage. The influence of four environmental factors (temperature, rainfall, humidity, and wind speed), types of animals in each shed (cattle, pigs, and avians), and neighboring animals (those in the other sheds and their proximity) were assessed. A subsample of 130,670 out of a total of 224,153 specimens collected was identified and counted. The collections were predominantly female (76.9%), though males were also commonly collected (23.1%). The dominant species included C. peregrinus (97,098 individuals; 43.0%), C. oxystoma (55,579; 24.6%), C. arakawae (45,996; 20.4%), and C. imicola (15,703; 7.0%), while other species accounted for the remaining 9777 individuals (4.4%). Trap captures were strongly influenced by seasons and animal sheds. Cattle exhibited the greatest effect on the community, likely due to their large biomass. Humidity and temperature had a positive effect on trap captures, whereas wind speed exhibited a negative effect. Cattle positively influenced all major species, except for C. arakawae, which showed a positive association with avians. A "neighboring effect" was also observed. Additionally, the lowest dispersal ability of C. arakawae was suggested.

Isolation and comparative analysis of culturable bacterial communities associated with life stages, breeding and rearing substrates of Culicoides oxystoma Kieffer (Diptera: Ceratopogonidae) vector of bluetongue virus

Culicoides oxystoma Kieffer (Diptera: Ceratopogonidae) has been vectoring several arboviruses, protozoa and nematodes, leading to mortality and morbidity of livestock and wild ruminants in the tropics and subtropics. Insight into the bacterial communities associated with the vector species must be worked out. This work tries to inventorize the bacterial communities associated with this important vector species. Acquisition of gut microbiota may be the parental origin, while some are obtained through feeding during larval stages. Culicoides oxystoma possesses semi-aquatic life cycle strategies for egg-laying and larval survival. The bacteria associated with C. oxystoma were compared throughout (i) life stages: egg, larval instars, pupa, adult: male and female obtained from laboratory colony; (ii) field-collected adult: male and age-graded females; and (iii) natural breeding substrate and artificial rearing substrate. The culture-dependent bacteria were identified by Sanger sequencing of 16S rRNA, and haemolytic bacteria were screened on blood agar. Results show that Firmicutes and Proteobacteria are the predominant Phyla, of which Bacillus spp. was the most abundant across the life stages. Across the life history, Bacillus cereus, Bacillus pumilus, Bacillus tropicus, Lysinibacillus sp. and Paenibacillus sp. were retrieved routinely. Bacillus cereus and Alcaligenes faecalis were detected in the lab-reared specimens and shared between the natural breeding site and rearing medium. From the adults trapped across two locations, B. cereus, Bacillus flexus, A. faecalis, Enterococcus faecium and Pseudomonas sp. were isolated. The bacterial species associated with this vector may influence various physiological traits, such as vectorial capacity, digestion and larval development, which need further investigation.

Molecular detection of blood protozoa and identification of black flies of the Simulium varicorne species group (Diptera: Simuliidae) in Thailand

Species of the Simulium varicorne group in Thailand have veterinary significance as vectors of haemosporidian parasites. Accurate identification is, therefore, critical to the study of vectors and parasites. We used morphology and molecular markers to investigate cryptic genetic lineages in samples identified as Simulium chumpornense Takaoka & Kuvangkadilok, 2000. We also tested the efficiency of the nuclear internal transcribed spacer 2 (ITS2) marker for the identification of species in this group. Morphological examinations revealed that S. chumpornense lineage A is most similar to S. khelangense Takaoka, Srisuka & Saeung, 2022, with minor morphological differences. They are also genetically similar based on mitochondrial cytochrome c oxidase I (COI) sequences. Geographically, the sampling site where paratypes of S. khelangense were originally collected is <50 km from where S. chumpornense lineage A was collected. We concluded that cryptic lineage A of S. chumpornense is actually S. khelangense. COI sequences could not differentiate S. kuvangkadilokae Pramual and Tangkawanit, 2008 from S. chumpornense and S. khelangense. In contrast, ITS2 sequences provided perfect accuracy in the identification of these species. Molecular analyses of the blood protozoa Leucocytozoon and Trypanosoma demonstrated that S. khelangense carries L. shoutedeni, Leucocytozoon sp., and Trypanosoma avium. The Leucocytozoon sp. in S. khelangense differs genetically from that in S. asakoae Takaoka & Davies, 1995, signaling the possibility of vector-parasite specificity.

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

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.

Prevalence and Diversity of Blood Parasites (Plasmodium, Leucocytozoon and Trypanosoma) in Backyard Chickens (Gallus gallus domesticus) Raised in Southern Thailand

Avian malaria and leucocytozoonosis can cause fatal diseases, whereas avian trypanosomiasis is reported to be harmless in chickens. Backyard chickens can be infected by several pathogens, including blood parasites, that may shed to industrial poultry production, with a consequently higher economic impact. This study aimed to investigate the presence of several blood parasites (Plasmodium, Leucocytozoon and Trypanosoma) in backyard chickens raised in Southern Thailand, using PCR-based detection and microscopic methods. From June 2021 to June 2022, 57 backyard chickens were sampled. Fresh thin blood smears were prepared from 11 individuals, and buffy coat smears were prepared from 55 of them. Both thin blood smears and buffy coat smears were used for microscopic analysis. Two nested PCR protocols that amplify a fragment of cytochrome b (cytb) and small subunit rRNA (SSU rRNA) genes were used to identify Haemosporida and Trypanosoma parasites, respectively. The number of positive samples was higher with the application of nested PCR than when buffy coat smears were used. Three new Plasmodium lineages (GALLUS47-49) and thirteen Leucocytozoon lineages (GALLUS50-62) were found. Trophozoites, meronts and gametocytes of Plasmodium gallinaceum (GALLUS01) were present in one thin blood smear. All thin blood smears revealed Leucocytozoon infections, but only three samples were a single infection. These three samples revealed the presence of fusiform host cell-parasite complexes, of which the morphological features resembled those of Leucocytozoon macleani (possible synonym is Leucocytozoon sabrazesi), while the cytb showed that this parasite is closely related to the lineage GALLUS06-07, described as Leucocytozoon schouteni. The Trypanosoma prevalence was 33.33%; it was present in only one of the thin blood smears, and it resembles Trypanosoma calmettei. This study showed the prevalence of a high diversity of Plasmodium (64.91%) and Leucocytozoon (89.47%) in Thai chickens. Both nested-PCR and buffy coat smear can be used as the diagnostic tool for the testing of Plasmodium, Leucocytozoon and Trypanosoma for parasitic control in backyard chickens and poultry farms. The information on the parasite species that can be found in chickens raised in Southern Thailand was also considered as the baseline information for further study.

Natural infection with Leishmania (Mundinia) martiniquensis supports Culicoides peregrinus (Diptera: Ceratopogonidae) as a potential vector of leishmaniasis and characterization of a Crithidia sp. isolated from the midges

The prevalence of autochthonous leishmaniasis in Thailand is increasing but the natural vectors that are responsible for transmission remain unknown. Experimental in vivo infections in Culicoides spp. with Leishmania (Mundinia) martiniquensis and Leishmania (Mundinia) orientalis, the major causative pathogens in Thailand, have demonstrated that biting midges can act as competent vectors. Therefore, the isolation and detection of Leishmania and other trypanosomatids were performed in biting midges collected at a field site in an endemic area of leishmaniasis in Tha Ruea and a mixed farm of chickens, goats, and cattle in Khuan Phang, Nakhon Si Thammarat province, southern Thailand. Results showed that Culicoides peregrinus was the abundant species (>84%) found in both locations and only cow blood DNA was detected in engorged females. Microscopic examination revealed various forms of Leishmania promastigotes in the foregut of several C. peregrinus in the absence of bloodmeal remnants, indicating established infections. Molecular identification using ITS1 and 3'UTR HSP70 type I markers showed that the Leishmania parasites found in the midges were L. martiniquensis. The infection rate of L. martiniquensis in the collected flies was 2% in Tha Ruea and 6% in Khuan Phang, but no L. orientalis DNA or parasites were found. Additionally, organisms from two different clades of Crithidia, both possibly new species, were identified using SSU rRNA and gGAPDH genes. Choanomastigotes and promastigotes of both Crithidia spp. were observed in the hindgut of the dissected C. peregrinus. Interestingly, midges infected with both L. martiniquensis and Crithidia were found. Moreover, four strains of Crithidia from one of the clades were successfully isolated into culture. These parasites could grow at 37°C in the culture and infect BALB/c mice macrophages but no multiplication was observed, suggesting they are thermotolerant monoxenous trypanosomatids similar to Cr. thermophila. These findings provide the first evidence of natural infection of L. martiniquensis in C. peregrinus supporting it as a potential vector of L. martiniquensis.

Host Associations of Culicoides Biting Midges in Northeastern Kansas, USA

Culicoides biting midges (Diptera: Ceratopogonidae) are hematophagous flies that transmit several viruses of veterinary concern to livestock. Understanding blood feeding behaviors is integral towards identification of putative vector species and preventing the transmission of these pathogens. PCR-based blood meal analysis was conducted on 440 blood-engorged Culicoides midges collected in northeastern Kansas, with 316 (71.8%) returning non-human vertebrate identifications at the ≥95% identity match level. Broadly, Culicoides sonorensis, Culicoides stellifer, and Culicoides variipennis were found to feed heavily on mammalian hosts, while Culicoides crepuscularis and Culicoides haematopotus fed on avian hosts. The blood meals in all specimens were graded prior to DNA extraction to determine whether blood meal size or digestion status significantly impacted the likelihood of a quality host match. Size had a significant impact on the likelihood of a quality match at grades 3-5, whereas digestion only significantly impacted outcomes at the most extreme grade. These vector-host dynamics have not previously been studied in Culicoides collected in Kansas, which represents a unique tallgrass prairie biome within the United States that is heavily interspersed with livestock operations. Based on these data, the highly abundant species C. crepuscularis and C. haematopotus are unlikely to be major vectors of mammalian viruses.

Diversity, Abundance and Host Blood Meal Analysis of Culicoides Latreille (Diptera: Ceratopogonidae) from Cattle Pens in Different Land Use Types from Thailand

Biting midges of the genus Culicoides Latreille are significant pests and vectors that transmit pathogens to humans and other animals. Cattle are among the important livestock that can potentially be severely affected by Culicoides. In this study, we examined the species diversity, abundance, and host blood meal identification of biting midges in cattle pens located in three different land use types: villages, agricultural areas, and the forest edge. A total of 12,916 biting midges were collected, and most of these were from cattle pens located in villages (34%) and agricultural land (52%). Morphological identification revealed 29 Culicoides species. The most common species were C. oxystoma, C. mahasarakhamense, C. peregrinus, and C. shortti; taken together, these species represented >80% of all specimens collected. Despite midges being less numerous (14% of the total collection), cattle pens located near the forest showed greater diversity (23) than those from villages and agricultural areas. More diverse immature habitats and host blood sources from wildlife in nearby forests possibly explain the greater diversity in the cattle pens near the forest edge. Host blood meal analysis revealed that most (65%) biting midges had fed on buffalo despite the fact that this animal was much less numerous than cows or chickens. Relatively larger size and black-colored skin could be factors that make buffalo more attractive to biting midges than other host species. In this study, we also provided 67 DNA barcoding sequences of 13 species, three of which (C. flaviscutatus, C. geminus, and C. suzukii) were first reported from Thai specimens. DNA barcode analysis indicated cryptic diversity within C. hegneri and C. flavescens in Thailand, and thus, further investigation is required to resolve their species status.

Molecular Identification of Culicoides Species and Host Preference Blood Meal in the African Horse Sickness Outbreak-Affected Area in Hua Hin District, Prachuap Khiri Khan Province, Thailand

African horse sickness (AHS) was reported as an outbreak in Thailand in 2020. Hematophagous insects from the genus Culicoides are the suspected vector responsible for AHS transmission. Horses in Hua Hin district, Prachuab Khiri Khan province, Thailand, were affected and died from AHS in 2020. However, the potential Culicoides species and its host preference blood meal in the affected areas are unknown. To investigate the potential vectors of AHS, Culicoides were collected using ultraviolet light traps placed near horse stables. Six horse farms, including five farms with AHS history and one farm without AHS history, were included in this study. Morphological and molecular identification of the Culicoides species was performed. Polymerase chain reaction (PCR) targeting the cytochrome b oxidase I (COXI) gene for confirmation of the Culicoides species, identification of the prepronociceptin (PNOC) gene for host preference blood meal, and bidirectional sequencing were conducted. Consequently, 1008 female Culicoides were collected, consisting of 708 and 300 samples captured at positions A and B at a distance of <2 and >5 m from the horse, respectively. Twelve Culicoides species identified by morphology were noted, including C. oxystoma (71.92%), C. imicola (20.44%), C. actoni (2.28%), C. flavipunctatus (1.98%), C. asiana (0.99%), C. peregrinus (0.60%), C. huffi (0.60%), C. brevitarsis (0.40%), C. innoxius (0.30%), C. histrio (0.30%), C. minimus (0.10%), and C. geminus (0.10%). The PCR detection of the Culicoides COXI gene confirmed Culicoides species in 23 DNA samples. PCR targeting the PNOC gene revealed that the Culicoides collected in this study fed on Equus caballus (86.25%), Canis lupus familiaris (6.25%), Sus scrofa (3.75%), and Homo sapiens (3.75%) for their blood meal. Human blood was identified from two samples of C. oxystoma and a sample of C. imicola. Three dominant species including C. oxystoma, C. imicola, and C. actoni that were reported in the Hua Hin area prefer to feed on horse blood. Moreover, C. oxystoma, C. imicola, and C. bravatarsis also feed on canine blood. This study revealed the species of Culicoides in Hua Hin district, Thailand, after the AHS outbreak.

First Evidence of Co-Circulation of Emerging Leishmania martiniquensis, Leishmania orientalis, and Crithidia sp. in Culicoides Biting Midges (Diptera: Ceratopogonidae), the Putative Vectors for Autochthonous Transmission in Southern Thailand

Since 1996, autochthonous cases of emerging leishmaniasis caused by Leishmania (Mundinia) martiniquensis and Leishmania (Mundinia) orientalis have been more frequently reported, especially in the northern and southern parts of Thailand. However, the accurate identification of their natural vectors and reservoirs remains unconfirmed. Previous studies have suggested that these emerging parasites might be transmitted by other non-phlebotomine vectors. Herein, we speculated that Culicoides biting midges might act as the competent vectors responsible for autochthonous leishmaniasis in southern Thailand. In this research, 187 non-engorged, parous and gravid Culicoides females and 47 blood-engorged ones were trapped from the residences of two recently diagnosed visceral leishmaniasis patients in Sadao District and the unaffected site in Rattaphum District, Songkhla Province, southern Thailand. Species diversity and abundance of biting midges varied among the trapping sites. Using ITS1-PCR and BLASTn analysis, L. martiniquensis was predominantly detected in several Culicoides species, including C. peregrinus, C. oxystoma, C. mahasarakhamense, and C. huffi from the vicinity of patients’ houses; and in C. fordae and C. fulvus from the unaffected site. L. orientalis was also co-circulated in C. peregrinus and C. oxystoma caught near the second patient’s house. Additionally, Crithidia sp. were also detected using SSU rRNA-PCR across Culicoides spp. Host blood meal analysis of eight different Culicoides species from the unaffected site also revealed that all trapped Culicoides had fed on cows and goats, indicating the possible role of these mammalian species as reservoir hosts. Essentially, this study is the first entomological investigation, revealing the co-circulation of emerging trypanosomatids among several species of Culicoides biting midges and strongly supporting the potential role of this insect group as the main vectors responsible for the epidemiology of autochthonous leishmaniasis in southern Thailand.