Genetic diversity and population structure of Phlebotomus argentipes: Vector of Leishmania donovani in Sri Lanka

Comparative Analysis of Phlebotomus argentipes Vector of Leishmaniasis in India and Sri Lanka

Phlebotomus argentipes is the predominant sandfly vector of leishmaniasis in the Indian subcontinent. India and Sri Lanka primarily report visceral and cutaneous leishmaniasis caused by Leishmania donovani. We compared Ph. argentipes from two locations, focusing on its morphological, molecular, and salivary protein characteristics. Sandflies were captured using CDC light traps and cattle-baited net traps. Species identification and morphological comparisons were carried out using standard taxonomic keys. DNA extracted from 12 Sri Lankan sandfly samples was PCR-amplified and sequenced for the variable region of Cytochrome oxidase subunit I. Existing DNA sequences of India from GenBank were utilized for a phylogenetic analysis between Sri Lanka and India. Salivary protein profiles were studied using SDS-PAGE, Western blot, and electrospray ionization/LC/MS/MS. The morphological similarities observed between female Ph. argentipes from India and Sri Lanka suggest the presence of Ph. argentipes var. glaucus. A phylogenetic analysis showed genetic divergence between Ph. argentipes populations, but both shared a similar salivary protein profile. A common, strong 30 kDa immunogenic band comprised PagSP05, PagSP06, and PagSP17 proteins of Ph. argentipes. The similarity between the immunogenic salivary proteins suggests their potential use as common markers for vector exposure or immune response stimulants across regions. The use of multiple samples for each category of serum would improve the comprehensiveness of the immunogenic profiles obtained.

Population genetics and genomics of Triatoma brasiliensis (Hemiptera, Reduviidae) in an area of high pressure of domiciliary infestation in Northeastern Brazil

Understanding the population dynamics of vectors is crucial for effective control of vector-borne diseases. In the Northeastern Brazilian semi-arid region, Triatoma brasiliensis persists as the most significant Chagas disease vector, frequently displaying recurrent domiciliary infestations. This situation raises relevant public health concerns in the municipality of Currais Novos in the state of Rio Grande do Norte. This area has experienced a high prevalence of peridomiciliary re-infestations by T. brasiliensis, coupled with elevated rates of Trypanosoma cruzi infection. Therefore, we assessed the distribution of genetic variation via mitochondrial Cytochrome b gene (MT-CYB) sequencing (n = 109) and single nucleotide polymorphisms (SNPs, n = 86) to assess the gene flow among distinct populations distributed in varied geographic spots and environments, mainly sylvatic and peridomiciliary. Insects were collected from rural communities at Currais Novos, enclosed within a 16 km radius. Sampling included 13 populations: one intradomiciliary, eight peridomiciliary, and four sylvatic. Furthermore, an external population located 220 km from Currais Novos was also included in the study. The method employed to obtain SNP information relied on ddRAD-seq genotyping-by-sequencing (GBS), enabling a genome-wide analysis to infer genetic variation. Through AMOVA analysis of MT-CYB gene variation, we identified four distinct population groups with statistical significance (FCT= 0.42; p<0.05). We identified a total of 3,013 SNPs through GBS, with 11 loci showing putative signs of being under selection. The variation based on 3,002 neutral loci evidenced low genetic structuration based on low FST values (p>0.05), indicating local panmixia. However, resampling algorithms pointed out that three samples from the external population were assigned (>98 %) in a cluster contrasting from the ones putatively under local panmixia - validating the newly applied genome-wide marker for studies on the population genetics at finer-scale resolution for T. brasiliensis. The presence of population structuring in some of the sampled points, as suggested by the mitochondrial marker, leads us to assume that infestations were probably initiated by small populations of females - demographic event poses a risk for rapid re-infestations. The local panmictic pattern revealed by the GBS marker poses a challenge for vector control measures, as re-infestation foci may be distributed over a wide geographical and ecological range. In such instances, vectors exhibit reduced susceptibility to conventional insecticide spraying operations since sylvatic populations are beyond the reach of these interventions. The pattern of infestation exhibited by T. brasiliensis necessitates integrating innovative strategies into the existing control framework, holding the potential to create a more resilient and adaptive vector control program. In our dataset, the results demonstrated that the genetic signals from both markers were complementary. Therefore, it is essential to consider the nature and inheritance pattern of each marker when inferring the pattern of re-infestations.

Genetic diversity and phylogeography of Phlebotomus argentipes (Diptera: Psychodidae, Phlebotominae), using COI and ND4 mitochondrial gene sequences

BACKGROUND

Phlebotomus argentipes complex is the primary vector for cutaneous leishmaniasis, a burgeoning health concern in contemporary Sri Lanka, where effective vector control is important for proper disease management. Understanding the genetic diversity of the P. argentipes population in Sri Lanka is vital before implementing a successful vector control program. Various studies have indicated that genetic divergence, caused by genetic drift or selection, can significantly influence the vector capacity of arthropod species. To devise innovative control strategies for P. argentipes, exploring genetic diversity and phylogeography can offer valuable insights into vector competence, key genetic trait transfer, and impact on disease epidemiology. The primary objective is to analyze the genetic diversity and phylogeography of the P. argentipes complex in Sri Lanka, based on two mitochondrial genomic regions in modern representatives of P. argentipes populations.

METHODOLOGY

A total of 159 P. argentipes specimens were collected from five endemic areas of cutaneous leishmaniasis and identified morphologically. Two mitochondrial regions (Cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 4 (ND4) were amplified using the total DNA and subsequently sequenced. Partial sequences of those mitochondrial genes were utilized to analyze genetic diversity indices and to explore phylogenetic and phylogeographic relationships.

PRINCIPAL FINDINGS

Among five sampling locations, the highest genetic diversity for COI and ND4 was observed in Hambantota (Hd-0.749, π-0.00417) and Medirigiriya (Hd-0.977, π-0.01055), respectively. Phylogeographic analyses conducted using COI sequences and GenBank retrieved sequences demonstrated a significant divergence of P. argentipes haplotypes found in Sri Lanka. Results revealed that they have evolved from the Indian ancestral haplotype due to historical- geographical connections of the Indian subcontinent with Sri Lanka.

CONCLUSIONS

Utilizing high-mutation-rate mitochondrial genes, such as ND4, can enhance the accuracy of genetic variability analysis in P. argentipes populations in Sri Lanka. The phylogeographical analysis of COI gene markers in this study provides insights into the historical geographical relationship between India and P. argentipes in Sri Lanka. Both COI and ND4 genes exhibited consistent genetic homogeneity in P. argentipes in Sri Lanka, suggesting minimal impact on gene flow. This homogeneity also implies the potential for horizontal gene transfer across populations, facilitating the transmission of genes associated with traits like insecticide resistance. This dynamic undermines disease control efforts reliant on vector control strategies.

Impact of climatic factors on temporal variability of sand fly abundance in Sri Lanka: Longitudinal study (2018 to 2020) with two-stage hierarchical analysis

Background

Phlebotomine sand flies serve as vectors for leishmaniasis, a major health concern, but a neglected tropical disease. The risk of vector activity is governed by climatic factors that vary in different geographic zones in the country. Thus, we aimed to quantify the effect of climatic variables on sand fly vector activity in ten sentinel sites across Sri Lanka.

Methods

Mean rainfall, ambient temperature, relative humidity, wind speed, soil temperature, evaporation, sunshine hours, and vector densities were recorded at monthly intervals in each location from March 2018 to February 2020. The association between weather variables and sand fly densities was analysed using a two-staged hierarchical procedure; Distributed Lag Non-Linear (DLNM) modelling framework and the DLNM method implemented in the R package dlnm (version number 2.4.6).

Results

Moderate rainfall values up to 120 mm per month and increasing RH up to 82 at lag of 0 months along with increasing soil temperature and evaporation rate at lag of 2 months were associated with statistically significant increase in the sand fly activity. These associations were heterogeneous across study settings. Whereas increasing ambient and soil temperature, sunshine hours, evaporation rate appeared to reduce the sand fly activity homogeneously at lag of 0 month in all the study settings.

Conclusions

The abundance of sand fly vectors varied in relation to selected climatic variables, either in real-time or with a time lag. This information can be utilized for predicting sand fly densities and for the development of effective strategies to prevent leishmaniasis transmission in specific settings.

Molecular taxonomy of phlebotomine sand flies (Diptera, Psychodidae) with emphasis on DNA barcoding: A review

The taxonomy and systematics of sand flies (Diptera, Psychodidae, Phlebotominae) are one of the pillars of research aimed to identifying vector populations and the agents transmitted by these insects. Traditionally, the use of morphological traits has been the main line of evidence for the definition of species, but the use of DNA sequences is useful as an integrative approach for their delimitation. Here, we discuss the current status of the molecular taxonomy of sand flies, including their most sequenced molecular markers and the main results. Only about 37% of all sand fly species have been processed for any molecular marker and are publicly available in the NCBI GenBank or BOLD Systems databases. The genera Phlebotomus, Nyssomyia, Psathyromyia and Psychodopygus are well-sampled, accounting for more than 56% of their sequenced species. However, less than 34% of the species of Sergentomyia, Lutzomyia, Trichopygomyia and Trichophoromyia have been sampled, representing a major gap in the knowledge of these groups. The most sequenced molecular markers are those within mtDNA, especially the DNA barcoding fragment of the cytochrome c oxidase subunit I (coi) gene, which has shown promising results in detecting cryptic diversity within species. Few sequences of conserved genes have been generated, which hampers higher-level phylogenetic inferences. We argue that sand fly species should be sequenced for at least the coi DNA barcoding marker, but multiple markers with different mutation rates should be assessed, whenever possible, to generate multilocus analysis.