Phylogenetic analysis of Phlebotomus species belonging to the subgenus Larroussius (Diptera, psychodidae) by ITS2 rDNA sequences

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.

Phlebotomine sand flies (Diptera, Psychodidae) from Spain: an updated checklist and extended distributions

Phlebotomine sand flies (Diptera: Psychodidae) are the natural vectors of Leishmania spp. (Kinetoplastida: Trypanosomatidae) and phleboviruses (Bunyavirales: Phenuiviridae). In Spain, these vectors appear to be increasing their geographical distribution and have serious repercussions on public and veterinary health, encouraging studies of sand flies and their associated pathogens. An up-to-date and easily accessible compendium of current and historical data on their presence and detailed distribution is a crucial step towards the development and implementation of appropriate preventive strategies. A checklist on the presence and distribution of sand flies in Spain is compiled from data extracted from a comprehensive review of scientific literature published between 1909 and 2021 and our new records on the presence of sand flies specimens collected under the entomological surveillance of bluetongue vectors from the Spanish Ministry of Agriculture, Fishery and Food (MAPA) during the period 2004-2021. In total, 13 Spanish species of sand flies (two of them with controversial status) belonging to two genera and six subgenera are presented in this updated checklist, including new distribution data for seven species, among which several stand out as confirmed or suspected vectors of Leishmaniainfantum: Phlebotomusariasi, Ph.langeroni, Ph.mascittii, and Ph.perniciosus.

Phlebotomine sand flies (Diptera: Psychodidae) of the Maghreb region: A systematic review of distribution, morphology, and role in the transmission of the pathogens

BACKGROUND

Phlebotomine sand flies (Diptera: Psychodidae) are important vectors of various human and animal pathogens such as Bartonella bacilliformis, Phlebovirus, and parasitic protozoa of the genus Leishmania, causative agent of leishmaniases that account among most significant vector-borne diseases. The Maghreb countries Mauritania, Morocco, Algeria, Tunisia, and Libya occupy a vast area of North Africa and belong to most affected regions by these diseases. Locally varying climatic and ecological conditions support diverse sand fly fauna that includes many proven or suspected vectors. The aim of this review is to summarize often fragmented information and to provide an updated list of sand fly species of the Maghreb region with illustration of species-specific morphological features and maps of their reported distribution.

MATERIALS AND METHODS

The literature search focused on scholar databases to review information on the sand fly species distribution and their role in the disease transmissions in Mauritania, Morocco, Algeria, Tunisia, and Libya, surveying sources from the period between 1900 and 2020. Reported distribution of each species was collated using Google Earth, and distribution maps were drawn using ArcGIS software. Morphological illustrations were compiled from various published sources.

RESULTS AND CONCLUSIONS

In total, 32 species of the genera Phlebotomus (Ph.) and Sergentomyia (Se.) were reported in the Maghreb region (15 from Libya, 18 from Tunisia, 23 from Morocco, 24 from Algeria, and 9 from Mauritania). Phlebotomus mariae and Se. africana subsp. asiatica were recorded only in Morocco, Ph. mascitti, Se. hirtus, and Se. tiberiadis only in Algeria, whereas Ph. duboscqi, Se. dubia, Se. africana africana, Se. lesleyae, Se. magna, and Se. freetownensis were reported only from Mauritania. Our review has updated and summarized the geographic distribution of 26 species reported so far in Morocco, Algeria, Tunisia, and Libya, excluding Mauritania from a detailed analysis due to the unavailability of accurate distribution data. In addition, morphological differences important for species identification are summarized with particular attention to closely related species such as Ph. papatasi and Ph. bergeroti, Ph. chabaudi, and Ph. riouxi, and Se. christophersi and Se. clydei.

Identification of Nilaparvata lugens and Its Two Sibling Species (N. bakeri and N. muiri) by Direct Multiplex PCR

The brown planthopper (BPH), Nilaparvata lugens (Stål) (Hemiptera: Delphacidae), is a destructive rice pest of Asia. Currently, one important monitoring method of BPH is through black light trapping. However, two sibling species of BPH, Nilaparvata bakeri (Muri) and Nilaparvata muiri China, can also be trapped by black light, and these species feed only on gramineous weeds rather than on rice. Therefore, the accurate identification of Nilaparvata species is crucial for N. lugens forecasting and management. The traditional morphological identification method is not feasible for subadults and damaged specimens. Furthermore, this error-prone morphological identification method is time and labor intensive, with the need for expertise and experience. Here, we established a direct multiplex polymerase chain reaction (dmPCR) assay using crude tissue fluid as a template, omitting purified DNA extraction. The crude tissue fluid can be obtained by grinding specimens without any biological reagent but only using distilled water. This dmPCR assay, using three pairs of diagnostic primers, is based on internal transcribed spacers (ITS). Each primer pair amplifies a species-specific fragment of a different size, which were easily and reliably separated in a 2% agarose gel. Furthermore, the dmPCR was verified to be applicable to damaged tissue specimens, such as head, thorax, or abdomen. In conclusion, this dmPCR assay is a novel, time-saving, cost-effective, and easy-to-apply molecular diagnostic method for the identification of the above three sibling species, N. lugens, N. bakeri, and N. muiri.

Identification of wild-caught phlebotomine sand flies from Crete and Cyprus using DNA barcoding

Background

Phlebotomine sand flies (Diptera: Psychodidae) are vectors of Leishmania spp., protozoan parasites responsible for a group of neglected diseases called leishmaniases. Two sand fly genera, Phlebotomus and Sergentomyia, contain species that are present in the Mediterranean islands of Crete and Cyprus where the visceral (VL), cutaneous (CL) and canine (CanLei) leishmaniases are a public health concern. The risk of transmission of different Leishmania species can be studied in an area by monitoring their vectors. Sand fly species are traditionally identified using morphological characteristics but minute differences between individuals or populations could be overlooked leading to wrong epidemiological predictions. Molecular identification of these important vectors has become, therefore, an essential tool for research tasks concerning their geographical distribution which directly relates to leishmaniasis control efforts. DNA barcoding is a widely used molecular identification method for cataloguing animal species by sequencing a fragment of the mitochondrial gene encoding cytochrome oxidase I.

Results

DNA barcoding was used to identify individuals of five sand fly species (Phlebotomus papatasi, P. similis, P. killicki, Sergentomyia minuta, S. dentata) circulating in the islands of Crete and Cyprus during the years 2011–2014. Phlebotomus papatasi is a known vector of zoonotic CL in the Middle East and it is found in both islands. Phlebotomus similis is the suspected vector of Leishmania tropica in Greece causing anthroponotic CL. Phlebotomus killicki was collected in Cyprus for the first time. Sergentomyia minuta, found to present intraspecific diversity, is discussed for its potential as a Leishmania vector. Molecular identification was consistent with the morphological identification. It successfully identified males and females, which is difficult when using only morphological characters. A phylogenetic tree was constructed based on the barcodes acquired, representing their genetic relationships along with other species from the area studied. All individuals identified were clustered according to their species and subgenus.

Conclusions

Molecular identification of sand flies via DNA barcoding can accurately identify these medically important insects assisting traditional morphological tools, thus helping to assess their implication in Leishmania transmission.

Molecular Identification of Phlebotomus caucasicus and Phlebotomus mongolensis (Diptera: Psychodidae) in a Hyperendemic Area of Zoonotic Cutaneous Leishmaniasis in Iran

Phlebotomus caucasicus Marzinovsky and Phlebotomus mongolensis Sinton are morphologically similar sand fly species. Finding a reliable, fast, and simple method to differentiate these two sand flies is important in understanding their role in the transmission of Leishmania parasite. In our study, 20 specimens of male P. caucasicus, 4 specimens of male P. mongolensis, and 16 specimens of female of both species (Caucasicus group) were examined by polymerase chain reaction (PCR). The result shows identical patterns with a visible fragment of about 500 bp in size. In restriction fragment length polymorphism (RFLP), we observed identical patterns with TasI used as the restriction enzyme. After alignment with sequences of the ITS2 partial gene in GenBank, a perfect match was obtained for the P. mongolensis, but not for P. caucasicus whose sequence was not present in the GenBank. Based on the results of our study, the RFLP-PCR method with nucleotide gene fragment ITS2 was a rapid and reliable method for differentiating these sand fly species.

DNA sequencing confirms PCR-RFLP identification of wild caught Larroussius sand flies from Crete and Cyprus

Many Phlebotomine sand fly species (Diptera, Psychodidae) are vectors of the protozoan parasite Leishmania causing a group of diseases called the leishmaniases. The subgenus Larroussius includes sand fly vectors found in South East Mediterranean Basin responsible for Visceral (VL) and Cutaneous human leishmaniasis (CL). It is important to monitor these medically important insects in order to safely predict possible Leishmania transmission cycles. Leishmania infantum is endemic in the islands of Crete and Cyprus with increasing VL cases in humans and dogs and in Cyprus the newly introduced Leishmania donovani causes both VL and CL in humans. The morphological identification of the females of the subgenus Larroussius often presents difficulties. Morphology and COI PCR - RFLP were used to identify wild caught Larroussius sand flies belonging to Phlebotomus tobbi, P. perfiliewi, and P. neglectus species from Crete and Cyprus. The identification results were further confirmed by sequencing (DNA barcoding) and Bayesian phylogenetic analysis. COI PCR - RFLP, when correctly optimized and with respect to geographical origin, can serve as an initial patterning identification tool when large sand fly numbers need to be identified. It could accurately assign Larroussius females and males to their taxa overcoming the difficulties of morphological identification. Finally, DNA barcoding will contribute to a molecular identification database to be used for in-depth species studies.

Intragenomic variations of multicopy ITS2 marker in Agrodiaetus blue butterflies (Lepidoptera, Lycaenidae)

The eukaryotic ribosomal DNA cluster consists of multiple copies of three genes, 18S, 5. 8S and 28S rRNAs, separated by multiple copies of two internal transcribed spacers, ITS1 and ITS2. It is an important, frequently used marker in both molecular cytogenetic and molecular phylogenetic studies. Despite this, little is known about intragenomic variations within the copies of eukaryotic ribosomal DNA genes and spacers. Here we present data on intraindividual variations of ITS2 spacer in three species of Agrodiaetus Hübner, 1822 blue butterflies revealed by cloning technique. We demonstrate that a distinctly different intragenomic ITS2 pattern exists for every individual analysed. ITS2 sequences of these species show significant intragenomic variation (up to 3.68% divergence), setting them apart from each other on inferred phylogenetic tree. This variation is enough to obscure phylogenetic relationships at the species level.

Limits of a rapid identification of common Mediterranean sandflies using polymerase chain reaction-restriction fragment length polymorphism

A total of 131 phlebotomine Algerian sandflies have been processed in the present study. They belong to the species Phlebotomus bergeroti, Phlebotomus alexandri, Phlebotomus sergenti, Phlebotomus chabaudi, Phlebotomus riouxi, Phlebotomus perniciosus, Phlebotomus longicuspis, Phlebotomus perfiliewi, Phlebotomus ariasi, Phlebotomus chadlii, Sergentomyia fallax, Sergentomyia minuta, Sergentomyia antennata, Sergentomyia schwetzi, Sergentomyia clydei, Sergentomyia christophersi and Grassomyia dreyfussi. They have been characterised by sequencing of a part of the cytochrome b (cyt b), t RNA serine and NADH1 on the one hand and of the cytochrome C oxidase I of the mitochondrial DNA (mtDNA) on the other hand. Our study highlights two sympatric populations within P. sergenti in the area of its type-locality and new haplotypes of P. perniciosus and P. longicuspis without recording the specimens called lcx previously found in North Africa. We tried to use a polymerase chain reaction-restriction fragment length polymorphism method based on a combined double digestion of each marker. These method is not interesting to identify sandflies all over the Mediterranean Basin.

Comparative analysis of salivary gland transcriptomes of Phlebotomus orientalis sand flies from endemic and non-endemic foci of visceral leishmaniasis

Background

In East Africa, Phlebotomus orientalis serves as the main vector of Leishmania donovani, the causative agent of visceral leishmaniasis (VL). Phlebotomus orientalis is present at two distant localities in Ethiopia; Addis Zemen where VL is endemic and Melka Werer where transmission of VL does not occur. To find out whether the difference in epidemiology of VL is due to distant compositions of P. orientalis saliva we established colonies from Addis Zemen and Melka Werer, analyzed and compared the transcriptomes, proteomes and enzymatic activity of the salivary glands.

Methodology/Principal Findings

Two cDNA libraries were constructed from the female salivary glands of P. orientalis from Addis Zemen and Melka Werer. Clones of each P. orientalis library were randomly selected, sequenced and analyzed. In P. orientalis transcriptomes, we identified members of 13 main protein families. Phylogenetic analysis and multiple sequence alignments were performed to evaluate differences between the P. orientalis colonies and to show the relationship with other sand fly species from the subgenus Larroussius. To further compare both colonies, we investigated the humoral antigenicity and cross-reactivity of the salivary proteins and the activity of salivary apyrase and hyaluronidase.

Conclusions

This is the first report of the salivary components of P. orientalis, an important vector sand fly. Our study expanded the knowledge of salivary gland compounds of sand fly species in the subgenus Larroussius. Based on the phylogenetic analysis, we showed that P. orientalis is closely related to Phlebotomus tobbi and Phlebotomus perniciosus, whereas Phlebotomus ariasi is evolutionarily more distinct species. We also demonstrated that there is no significant difference between the transcriptomes, proteomes or enzymatic properties of the salivary components of Addis Zemen (endemic area) and Melka Werer (non-endemic area) P. orientalis colonies. Thus, the different epidemiology of VL in these Ethiopian foci cannot be attributed to the salivary gland composition.

Phlebotomus orientalis is the vector of visceral leishmaniasis (VL) caused by Leishmania donovani in Northeast Africa. Immunization with sand fly saliva or with individual salivary proteins has been shown to protect against leishmaniasis in different hosts, warranting the intensive study of salivary proteins of sand fly vectors. In our study, we characterize the salivary compounds of P. orientalis, thereby broadening the repertoire of salivary proteins of sand fly species belonging to the subgenus Larroussius. In order to find out whether there is any connection between the composition of P. orientalis saliva and the epidemiology of VL in two distinct Ethiopian foci, Addis Zemen and Melka Werer, we studied the transcriptomes, proteomes, enzymatic activities, and the main salivary antigens in two P. orientalis colonies originating from these areas. We did not detect any significant difference between the saliva of female sand flies originating in Addis Zemen (endemic area) and Melka Werer (non-endemic area). Therefore, the different epidemiology of VL in these Ethiopian foci cannot be related to the distant salivary gland protein composition. Identifying the sand fly salivary gland compounds will be useful for future research focused on characterizing suitable salivary proteins as potential anti-Leishmania vaccine candidates.

Genetic differentiation between sandfly populations of Phlebotomus chinensis and Phlebotomus sichuanensis (Diptera: Psychodidae) in China inferred by microsatellites

Background

Phlebotomus chinensis is a primary vector of visceral leishmaniasis; it occurs in various biotopes with a large geographical distribution, ranging from Yangtze River to northeast China. Phlebotomus sichuanensis, a species closely related to P. chinensis in high altitude regions, has a long term disputation on its taxonomic status. Both species occur in the current epidemic regions and are responsible for the transmission of leishmaniasis. Population genetic analysis will help to understand the population structure and infer the relationship for morphologically indistinguishable cryptic species. In this study, microsatellite markers were used for studying the genetic differentiation between P. chinensis and P. sichuanensis.

Methods

Sandflies were collected in 6 representative localities in China in 2005-2009. Ten microsatellite loci were used to estimate population genetic diversity. The intra-population genetic diversity, genetic differentiation and effective population size were estimated.

Results

All 10 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy-Weinberg disequilibrium was found in 23 out of 60 (38.33%) comparisons associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed three clusters. The cluster I included almost all specimens in the sample SCD collected at high altitude habitats in Sichuan. The other two clusters were shared by the remaining 5 populations, SCJ in Sichuan, GSZ in Gansu, SXL and SXX in Shaanxi and HNS in Henan. The diversity among these 5 populations was low (F_ST = -0.003-0.090) and no isolation by distance was detected. AMOVA analysis suggested that the variations were largely derived from individuals within populations and among individuals. Consistently, the analysis of ribosomal DNA second internal transcribed spacer (ITS2) sequence uncovered three types of variants, which corresponded with the three gene pools revealed by microsatellites.

Conclusions

The data suggested that the SCD population carried a distinct gene pool, which was differentiated from the other populations. The high altitude ecological habitats, distinctive ITS2 and herein divergence inferred by microsatellite loci support the species status of P. sichuanensis. The P. chinensis populations did not have a significant divergence from each another.

Incrimination of Phlebotomus kandelakii and Phlebotomus balcanicus as vectors of Leishmania infantum in Tbilisi, Georgia

A survey of potential vector sand flies was conducted in the neighboring suburban communities of Vake and Mtatsminda districts in an active focus of visceral Leishmaniasis (VL) in Tbilisi, Georgia. Using light and sticky-paper traps, 1,266 male and 1,179 female sand flies were collected during 2006-2008. Five Phlebotomus species of three subgenera were collected: Phlebotomus balcanicus Theodor and Phlebotomus halepensis Theodor of the subgenus Adlerius; Phlebotomus kandelakii Shchurenkova and Phlebotomus wenyoni Adler and Theodor of the subgenus Larroussius; Phlebotomus sergenti Perfil'ev of the subgenus Paraphlebotomus. Phlebotomus sergenti (35.1%) predominated in Vake, followed by P. kandelakii (33.5%), P. balcanicus (18.9%), P. halepensis (12.2%), and P. wenyoni (0.3%). In Mtatsminda, P. kandelakii (76.8%) comprised over three fourths of collected sand flies, followed by P. sergenti (12.6%), P. balcanicus (5.8%), P. halepensis (3.7%), and P. wenyoni (1.1%). The sand fly season in Georgia is exceptionally short beginning in early June, peaking in July and August, then declining to zero in early September. Of 659 female sand flies examined for Leishmania, 12 (1.8%) specimens without traces of blood were infected including 10 of 535 P. kandelakii (1.9%) and two of 40 P. balcanicus (5.0%). Six isolates were successfully cultured and characterized as Leishmania by PCR. Three isolates from P. kandelakii (2) and P. balcanicus (1) were further identified as L. infantum using sequence alignment of the 70 kDa heat-shock protein gene. Importantly, the sand fly isolates showed a high percent identity (99.8%-99.9%) to human and dog isolates from the same focus, incriminating the two sand fly species as vectors. Blood meal analysis showed that P. kandelakii preferentially feeds on dogs (76%) but also feeds on humans. The abundance, infection rate and feeding behavior of P. kandelakii and the infection rate in P. balcanicus establish these species as vectors in the Tbilisi VL focus.

Application of molecular techniques in the study of natural infection of Leishmania infantum vectors and utility of sandfly blood meal digestion for epidemiological surveys of leishmaniasis

Epidemiological studies on the distribution of leishmaniasis caused by Leishmania infantum Nicolle, 1908 (Kinetoplastida: Trypanosomatidae) have been based principally on serological surveys of the canine reservoir. This methodology is useful due to the facility of sampling, the rapidity in obtaining results, its consistency and because it allows the detection of heterogeneous foci of canine leishmaniasis (CanL) even in small areas. Other investigations have analysed Leishmania parasitism in sandflies (Diptera: Psychodidae: Phlebotominae) by using classical dissection techniques. These techniques allow the vector species to be incriminated in different foci, although they suffer from being very time consuming. Lately, studies in this field are increasingly using molecular techniques, which are faster and easier to perform. In the present work, we applied a nested-PCR in a study of natural infection of sandflies by Leishmania in three isolated farms where serological data on canine leishmaniasis of local dogs were also obtained. The analysis allowed the detection of 38.7% of females with positive nested-PCR (78%, 18% and 0%, respectively, in the different isolated farms). The positive Leishmania DNA samples were genotyped and identified as L. infantum. The results of this work provide new data for the vectorial capacity of Phlebotomus ariasi in a Pyrenean area, which can be considered at risk of becoming a new focus of CanL. The females with positive nested-PCR displayed blood in the midgut at different degrees of digestion, and/or were gravid. According to the multivariate logistic regression analysis, the risk of nested-PCR-positivity increased significantly with the degree of blood digestion (OR = 1.3; P value = 0.025). The Phlebotomus species and the presence of eggs were not statistically associated with nested-PCR positivity (P value of >0.05). The correlation of positive nested-PCR results with the presence of seropositive dogs in the farm confirms the utility of this technique in the study of the distribution and intensity of leishmaniasis foci. Also, the importance of sandfly blood-meal digestion for epidemiological surveys of leishmaniasis foci has been demonstrated.

Multilocus molecular and phylogenetic analysis of phlebotomine sand flies (Diptera: Psychodidae) from southern Italy

This study reports a combined analysis of mitochondrial and ribosomal DNA target regions of phlebotomine sand flies (Diptera: Psychodidae) from the Mediterranean region. A ∼900 bp long fragment of the mitochondrial DNA encompassing regions within cytb and nd1 gene and the complete ITS2 ribosomal region (∼500 bp) were sequenced and characterized for Phlebotomus perniciosus, Phlebotomus perfiliewi, Phlebotomus neglectus, Phlebotomus papatasi, and Sergentomyia minuta, captured in two sites of southern Italy. From one to eight mitochondrial haplotypes and from one to three ITS2 sequence types were found for the examined specimens according to the different sand fly species. The mean interspecific difference in the mitochondrial sequences was of 16.1%, with an overall intraspecific nucleotide variation from 0.1 to 2.8%. A higher interspecific difference (mean 25.1%) was recorded for the ITS2 sequence, with an overall intraspecific nucleotide variation up to 4.9%. The sequence types alignment of ITS2 region showed that all phlebotomine specimens possessed a split 5.8S rRNA, consisting of a mature 5.8S rRNA and a 2S rRNA separated by a short transcribed spacer. Phylogenetic analysis of the Phlebotomus spp. sequences, herein determined and of those available in GenBank™ were concordant in clustering P. neglectus, P. perfiliewi and P. papatasi with the same species collected from different geographic areas of the Mediterranean basin in four main clades for mtDNA and ITS2, respectively. This study demonstrates the utility of multilocus sequencing, provides a dataset for the molecular identification of the most prevalent phlebotomine sand flies in southern Europe and defines the phylogenetic relationships among species examined.

Intraspecific variability of natural populations of Phlebotomus sergenti, the main vector of Leishmania tropica

An intraspecific study of Phlebotomus sergenti was performed on populations from Turkey, Syria, Israel, and Uzbekistan by four different approaches: geometric morphometrics, RAPD analysis, internal transcribed spacer 2 (ITS2) sequencing (nuclear marker), and cytochrome B sequencing (mitochondrial marker). In RAPD analysis, distinct clades were formed in accordance with the geographical origin of the specimens. There was no distinct grouping according to place of origin within the Turkish samples from various localities in south-eastern Anatolia, which suggests a gene flow between populations separated spatially by the Amanos mountains, a mountain range of a considerable altitude. The results of ITS2 rDNA sequencing complied with the previously published intraspecific division of P. sergenti into two branches, northeastern and southwestern. However, mtDNA haplotypes formed three lineages with specimens from Turkey and Israel, sharing a common clade. A previously postulated hypothesis about a complex of sibling species within P. sergenti is therefore questionable. Cytochrome B seems to be a more discriminative marker for intraspecific variability assessment.

Morphometric and molecular analyses of the sand fly species Lutzomyia shannoni (Diptera: Psychodidae: Phlebotominae) collected from seven different geographical areas in the southeastern United States

A morphometric and molecular study of adult male and female Lutzomyia shannoni (Dyar 1929) collected at seven different locations within the southeastern United States was conducted to assess the degree of divergence between the grouped specimens from each location. The collection locations were as follows: Fort Bragg, NC; Fort Campbell, KY; Fort Rucker, AL; Ossabaw Island, GA; Patuxent National Wildlife Research Refuge, MD; Suwannee National Wildlife Refuge, FL; and Baton Rouge, LA. Forty males and forty females from each location were analyzed morphometrically from 54 and 49 character measurements, respectively. In addition, the molecular markers consisting of the partial cytochrome c oxidase subunit I (from 105 sand flies: 15 specimens/collection site) and the partial internal transcribed spacer 2 (from 42 sand flies: six specimens/collection site) were compared. Multivariate analyses indicate that the low degree of variation between the grouped specimens from each collection site prevents the separation of any collection site into an entity that could be interpreted as a distinct population. The molecular analyses were in concordance with the morphometric study as no collection location grouped into a separate population based on the two partial markers. The grouped specimens from each collection site appear to be within the normal variance of the species, indicating a single population in the southeast United States. It is recommended that additional character analyses of L. shannoni based on more molecular markers, behavioral, ecological, and physiological characteristics, be conducted before ruling out the possibility of populations or a cryptic species complex within the southeastern United States.

The sandflies (Diptera: Psychodidae, Phlebotominae) in military camps in northern Afghanistan (2007-2009), as identified by morphology and DNA 'barcoding'

As part of a continuous, standardized programme of monitoring the Leishmania vectors in German military camps in northern Afghanistan between 2007 and 2009, a detailed taxonomic analysis of the endemic sandfly fauna, as sampled using light and odour-baited traps, was conducted. Of the 10 sandfly species that were recorded, six may serve as enzootic and/or zooanthroponotic vectors of parasites causing human leishmaniasis. The use of a simple DNA-'barcoding' technique based on the mitochondrial cyt b gene, to identify the collected sandflies to species level, revealed (1) a clear discrimination between the potential vector species, (2) clustering of species within most subgenera, and (3) particularly high heterogeneity within the subgenus Paraphlebotomus (Phlebotomus alexandri being grouped with Ph. papatasi rather than with other Paraphlebotomus species). The data also indicate a high level of genetic heterogeneity within the subgenus Sergentomyia but close similarity between Sergentomyia sintoni and Sergentomyia murgabiensis. The morphological similarity of many medically important sandflies can make species identification difficult, if not impossible. The new DNA-barcoding techniques may provide powerful discriminatory tools in the future.

Genetic structures of sand fly (Diptera: Psychodidae) populations in a leishmaniasis endemic region of Turkey

The object of this study was to determine the genetic structures of three vector species, Phlebotomus tobbi, Phlebotomus papatasi, and Phlebotomus sergenti, in the Cukurova Region of Turkey, an endemic focus of cutaneous leishmaniasis. The genetic diversity indices, neutrality tests and hierarchical analysis of molecular variance (AMOVA) were performed using partial sequences of ITS2 and cytochrome b gene regions. In all species, within population genetic variation was higher than between population variation for ITS2 gene region. Fst values were low and non-significant for P. sergenti, and were higher for P. papatasi and P. tobbi indicating a weak structuring between populations. AMOVA tests suggest any substantial isolation between populations within species. AMOVA analysis of cyt b gene region revealed significant genetic structuring between populations for P. papatasi and P. sergenti. Fst values were relatively high and significant for these species indicating a certain degree of isolation between populations. However, in P. tobbi, any significant population genetic structuring was detected. Tajima's D and Fu's Fs values were negative and significant in all three species might be indicating a demographic expansion.

Morphometric and molecular differentiation of Phlebotomus (Phlebotomus) sandflies

The closely related sandfly species of the subgenus Phlebotomus namely, Phlebotomus papatasi (Scopoli, 1786), Phlebotomus duboscqi Neveu-Lemair, 1906 and Phlebotomus bergeroti Parrot, 1934 (Diptera: Psychodidae), are major vectors of Leishmania major (Kinetoplastida: Trypanosomatidae), the causative agent of cutaneous leishmaniasis in the Old World. Although allopatric in most of their distribution, the three species exist sympatrically in many places in central and eastern Sudan. Males of the three species can be distinguished using morphological characters; however, females are much harder to identify, thus complicating epidemiological studies. We carried out a morphometric and a molecular study to determine reliable morphological features and develop a polymerase chain reaction (PCR) assay for distinguishing females of these species. Males and females from each species were collected from sites in Sudan, East Africa and from one site in Mali, West Africa. Males were analysed morphologically and 20 characters and 10 character ratios were used in a stepwise discriminant analysis. This led to the identification of four characters with high discriminant loading scores sufficient for accurate male species identification. Male DNA was then used for the development of a PCR-based species diagnostic based on the second internal transcribed spacer (ITS2) of the ribosomal DNA. A set of four primers was developed to generate fragment sizes that are specific to each species and can reliably identify females as well as hybrid DNA. Both the morphometric and the molecular findings of this study have important applications for studies of the epidemiology of cutaneous leishmaniasis.

Longitudinal study on the detection of canine Leishmania infections by conjunctival swab analysis and correlation with entomological parameters

A longitudinal study was carried out on kennelled stray dogs in a canine leishmaniasis (CanL) endemic area, to evaluate early and late diagnostic performance of a non-invasive conjunctival swab (CS) nested (n)-PCR analysis for Leishmania detection in 2 cohorts of dogs, respectively. (A) Sixty-five IFAT- and CS n-PCR-negative dogs exposed to, and followed up once or twice a month during a full sand fly season (July-November 2008). In parallel, a sand fly survey was performed on site using standard sticky traps set twice a month, for a cumulative surface of 63 m(2). (B) Seventeen IFAT- and CS n-PCR-negative dogs found positive in July 2008 at the peripheral blood buffy-coat (BC) n-PCR. These dogs were examined again by BC n-PCR in September and November 2008, and before the subsequent transmission season (May 2009) along with CS n-PCR and IFAT. None of the cohort (A) dogs converted to positive CS n-PCR during the transmission season. Although approximately 2500 phlebotomine specimens were collected with peaks of 100-147 specimens/m(2) sticky trap, the cumulative density of the only proven CanL vector in the area (Phlebotomus perniciosus) was found to be very low (0.5/m(2)). All cohort (B) dogs remained substantially seronegative; BC n-PCR showed an intermittent positive trend during the period surveyed, resulting in 82% conversions to negative by the end of the study, in contrast with 71% conversions to positive at the CS n-PCR analysis. In conclusion, while CS n-PCR was not found effective for the early detection of Leishmania contacts in dogs exposed to a low pressure of vectorial transmission, this assay showed to slowly convert to positive in a high rate of dogs, in the absence of seroconversion. CS n-PCR technique can be a suitable marker for assessing Leishmania exposure in dogs as a non-invasive alternative to current serological and molecular tools.

Molecular epidemiology for vector research on leishmaniasis

Leishmaniasis is a protozoan disease caused by the genus Leishmania transmitted by female phlebotomine sand flies. Surveillance of the prevalence of Leishmania and responsive vector species in endemic and surrounding areas is important for predicting the risk and expansion of the disease. Molecular biological methods are now widely applied to epidemiological studies of infectious diseases including leishmaniasis. These techniques are used to detect natural infections of sand fly vectors with Leishmania protozoa and are becoming powerful tools due to their sensitivity and specificity. Recently, genetic analyses have been performed on sand fly species and genotyping using PCR-RFLP has been applied to the sand fly taxonomy. In addition, a molecular mass screening method has been established that enables both sand fly species and natural leishmanial infections to be identified simultaneously in hundreds of sand flies with limited effort. This paper reviews recent advances in the study of sand flies, vectors of leishmaniasis, using molecular biological approaches.

Distribution of Phlebotomus perniciosus in North-Italy: a study on 18S rDNA of phlebotomine sand flies

Leishmaniasis has recently spread, and is now endemic, in many parts of North Italy, even if it is not clear how sand flies vectors have reached this area. In order to clarify the origin of the Phlebotomus perniciosus, the main sand flies specie found in these areas, we analyzed and compared the 18S rDNA sequence from 33 out of 122 P. perniciosus collected in new endemic areas, from neighbor (<or=120 km) traditionally endemic area (Liguria) and from an area about 400 km far from the North West Italy. Based on the 18S rDNA sequence analysis, three different groups with different degrees of heterogeneity were identified. Two of them are suggestive of migration at a local scale and one, found in all sampled foci, is probably due to passive carriage. Results suggest that both migration on small scale from traditionally endemic area and passive carriage on long distance can contribute to sand flies colonization of new areas.

Genetic diversity of ribosomal RNA internal transcribed spacer sequences in Lutzomyia species from areas endemic for New World cutaneous leishmaniasis

In this study, each of 60 rRNA internal transcribed spacer (ITS) 1 and ITS2 sequences was determined from 44 individuals of 14 morphologically identified New World sand fly Lutzomyia species in Ecuador, and their interspecies and intraspecies genetic diversity was compared. Distinguishing between related species based on the ITS1 sequence was difficult because of variability, while the genetic diversity of ITS2 was distinct even among closely related species. Further, an assessment of intraspecies ITS sequence diversity in the subgenus Helcocyrtomyia revealed no correlation between sequence variation and geographic distribution. The results strongly suggested ITS2 to be a more suitable marker than ITS1 for the taxonomic analysis of Lutzomyia species including closely related species. Moreover, neither ITS sequence may be useful for the analysis of population structures in Lutzomyia species.

Characterization of internal transcribed spacer (ITS1)-ITS2 region of ribosomal RNA gene from 25 species of Culicoides biting midges (Diptera: Ceratopogonidae) in Japan

We determined nucleotide sequences of the nuclear rDNA internal transcribed spacer (ITS)1-5.8S-ITS2a-2S-ITS2 region in 103 individuals of 25 Culicoides species (Diptera: Ceratopogonidae) from 11 locations in Japan. Ribosomal RNA genes, 5.8S and 2S rDNA, were highly conserved among the species with few variations. The ITS2a region showed length variation among species. Both ITS1 and ITS2 showed highly varied sequences among species. The noticeable indel regions among ITS1 sequences are present in some Culicoides species, separating species into two types having long or short ITS1 region. However, Culicoides cylindratus Kitaoka possesses both types of ITS1 in each individual; these results seem to indicate that the ITS1-long type was the prototype and the short type was produced through deletion in many Culicoides species. One species, belonging to subgenus Avaritia, possessed an Avaritia-specific sequence in ITS1 and phylogenetically formed a monophyletic group. Geographical genotypes in a species were not clear. Species-specific sequence features were observed, enabling molecular identification of Culicoides species.

Molecular typing of sand fly species (Diptera, Psychodidae, Phlebotominae) from areas endemic for Leishmaniasis in Ecuador by PCR-RFLP of 18S ribosomal RNA gene

Surveillance of the distribution of sand fly species is important for prediction of the risk and expansion of Leishmania infection in endemic and surrounding areas. In the present study, a simple and reliable method of typing New World Lutzomyia species circulating in endemic areas in Ecuador was established by using polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) technique. PCR-RFLP of 18S ribosomal RNA (rRNA) genes with the restriction enzyme AfaI and subsequently HinfI successfully identified seven sand fly species in nine endemic areas in Ecuador. Although intraspecific genetic-diversity affecting the RFLP-patterns was detected in a species, the patterns were species specific. The method promises to be a powerful tool for the classification of New World Lutzomyia species.

Molecular evolution and phylogenetic utility of the internal transcribed spacer 2 (ITS2) in Calyptratae (Diptera: Brachycera)

The resolution potential of internal transcribed spacer 2 (ITS2) at deeper levels remains controversial. In this study, 105 ITS2 sequences of 55 species in Calyptratae were analyzed to examine the phylogenetic utility of the spacer above the subfamily level and to further understand its evolutionary characteristics. We predicted the secondary structure of each sequence using the minimum-energy algorithm and constructed two data matrixes for phylogenetic analysis. The ITS2 regions of Calyptratae display strong A-T bias and slight variation in length. The tandem and dispersed repeats embedded in the spacers possibly resulted from replication slippage or transposition. Most foldings conformed to the four-domain model. Sequence comparison in combination with the secondary structures revealed six conserved motifs. Covariation analysis from the conserved motifs indicated that the secondary structure restrains the sequence evolution of the spacer. The deep-level phylogeny derived from the ITS2 data largely agreed with the phylogenetic hypotheses from morphologic and other molecular evidence. Our analyses suggest that the accordant resolutions generated from different analyses can be used to infer deep-level phylogenetic relations.

Intraspecific variability (rDNA ITS and mtDNA Cyt b) of Phlebotomus sergenti in Spain and Morocco

Phlebotomus sergenti, the main vector of Leishmania tropica, is widely represented in Spain, whilst L. tropica is not an endemic species in this country. Nevertheless, the important human migrant flow from regions where L. tropica is endemic, the existence of its vector and the anthroponotic cycle of the parasite could lead to its establishment in our country. The vectorial capacity of P. sergenti could depend on the existence of cryptic species or races, which can only be identified by molecular methods. Our aim has been to study the molecular variation on two genes of Spanish populations of P. sergenti, comparing them with Moroccan ones. We have identified 5 ribosomal and 16 mitochondrial haplotypes. The results obtained indicate a high diversity of P. sergenti in Spain and the existence in the country of two P. sergenti lineages, a typically Spanish mitochondrial lineage and another one common in Morocco, where L. tropica is endemic in the south and emergent in the north and centre of the country.

Molecular phylogeny of black flies (Diptera: Simuliidae) from Thailand, using ITS2 rDNA

The sequences of the second internal transcribed spacer (ITS2) of ribosomal DNA (rDNA) were determined for 40 black fly species from Thailand, belonging to 4 subgenera of the genus Simulium, namely Gomphostilbia (12 species), Nevermannia (5 species), Montisimulium (1 species), Simulium sensu stricto (21 species), and an unknown subgenus with one species (Simulium baimaii). The length of the ITS2 ranged from 247 to 308 bp. All black fly species had high AT content, ranging from 71 to 83.8%. Intraindividual variation (clonal variation) occurred in 13 species, ranging from 0.3 to 1.1%. Large intrapopulation and interpopulation heterogeneities exist in S. feuerboni from the same and different locations in Doi Inthanon National Park, northern Thailand. Phylogenetic relationships among 40 black fly species were examined using PAUP (version 4.0b10) and MrBAYS (version 3.0B4). The topology of the trees revealed two major monophyletic clades. The subgenus Simulium and Simulium baimaii were placed in the first monophyletic clade, whereas the subgenera Nevermannia + Montisimulium were placed as the sister group to the subgenus Gomphostilbia in the second monophyletic clade. Our results suggest that S. baimaii belongs to the malyschevi-group or variegatum-group in the subgenus Simulium. The molecular phylogeny generally agrees with existing morphology-based phylogenies.

Evolution and phylogenetic information content of the ribosomal DNA repeat unit in the Blattodea (Insecta)

The organization, structure, and nucleotide variability of the ribosomal repeat unit was compared among families, genera, and species of cockroaches (Insecta:Blattodea). Sequence comparisons and molecular phylogenetic analyses were used to describe rDNA repeat unit variation at differing taxonomic levels. A reverse similar 1200 bp fragment of the 28S rDNA sequence was assessed for its potential utility in reconstructing higher-level phylogenetic relationships in cockroaches. Parsimony and maximum likelihood analyses of these data strongly support the expected pattern of relationships among cockroach groups. The examined 5' end of the 28S rDNA is shown to be an informative marker for larger studies of cockroach phylogeny. Comparative analysis of the nucleotide sequences of the rDNA internal transcribed spacers (ITS1 and ITS2) among closely related species of Blattella and Periplaneta reveals that ITS sequences can vary widely in primary sequence, length, and folding pattern. Secondary structure estimates for the ITS region of Blattella species indicate that variation in this spacer region can also influence the folding pattern of the 5.8S subunit. These results support the idea that ITS sequences play an important role in the stability and function of the rRNA cluster.

ITS 2 sequences heterogeneity in Phlebotomus sergenti and Phlebotomus similis (Diptera, Psychodidae): possible consequences in their ability to transmit Leishmania tropica

An intraspecific study on Phlebotomus sergenti, the main and only proven vector of Leishmania tropica among the members of the subgenus Paraphlebotomus was performed. The internal transcribed spacer 2 (ITS2) sequences of 12 populations from 10 countries (Cyprus, Egypt, Italy, Lebanon, Morocco, Pakistan, Portugal, Spain, Syria, and Turkey) were compared. Samples also included three species closely related to P. sergenti: Phlebotomus similis (three populations from Greece and Malta), Phlebotomus jacusieli and Phlebotomus kazeruni. Our results confirm the validity of the taxa morphologically characterised, and imply the revision of their distribution areas, which are explained through biogeographical events. At the Miocene time, a migration route, north of the Paratethys sea would have been followed by P. similis to colonise the north of the Caucasus, Crimea, Balkans including Greece and its islands, and western Turkey. Phlebotomus sergenti would have followed an Asiatic dispersion as well as a western migration route south of the Tethys sea to colonise North Africa and western Europe. This hypothesis seems to be well supported by high degree of variation observed in the present study, which is not related to colonisation or to intra-populational variation. Two groups can be individualised, one oriental and one western in connection with ecology, host preferences and distribution of L. tropica. We hypothesise that they could be correlated with differences in vectorial capacities.