Genetic markers for study of the anopheline vectors of human malaria

Genetic diversity and phylogeographic patterns of the peacock jewel-damselfly, Rhinocypha fenestrella (Rambur, 1842)

Despite is known to have widespread distribution and the most active species of the family Chlorocyphidae, the molecular data of Rhinocypha fenestrella (Rambur, 1842) are relatively scarce. The present study is the first that examined the genetic diversity and phylogeographic pattern of the peacock jewel-damselfly R. fenestrella by sequencing the cytochrome C oxidase I (cox1) and 16S rRNA gene regions from 147 individuals representing eight populations in Malaysia. A total of 26 and 10 unique haplotypes were revealed by the cox1 and 16S rRNA genes, respectively, and 32 haplotypes were recovered by the concatenated sequences of cox1+16S. Analyses indicated that haplotype AB2 was the most frequent and the most widespread haplotype in Malaysia while haplotype AB1 was suggested as the common ancestor haplotype of the R. fenestrella that may arose from the Negeri Sembilan as discovered from cox1+16S haplotype network analysis. Overall haplotype and nucleotide diversities of the concatenated sequences were Hd = 0.8937 and Pi = 0.0028, respectively, with great genetic differentiation (FST = 0.6387) and low gene flow (Nm = 0.14). Population from Pahang presented the highest genetic diversity (Hd = 0.8889, Pi = 0.0022, Nh = 9), whereas Kedah population demonstrated the lowest diversity (Hd = 0.2842, Pi = 0.0003, Nh = 4). The concatenated sequences of cox1+16S showed genetic divergence ranging from 0.09% to 0.97%, whereas the genetic divergence for cox1 and 16S rRNA genes were 0.16% to 1.63% and 0.01% to 0.75% respectively. This study provides for the first-time insights on the intraspecific genetic diversity, phylogeographic pattern and ancestral haplotype of Rhinocypha fenestrella. The understanding of molecular data especially phylogeographic pattern can enhance the knowledge about insect origin, their diversity, and capability to disperse in particular environments.

Revolutionizing Malaria Vector Control: The Importance of Accurate Species Identification through Enhanced Molecular Capacity

Many factors, such as the resistance to pesticides and a lack of knowledge of the morphology and molecular structure of malaria vectors, have made it more challenging to eradicate malaria in numerous malaria-endemic areas of the globe. The primary goal of this review is to discuss malaria vector control methods and the significance of identifying species in vector control initiatives. This was accomplished by reviewing methods of molecular identification of malaria vectors and genetic marker classification in relation to their use for species identification. Due to its specificity and consistency, molecular identification is preferred over morphological identification of malaria vectors. Enhanced molecular capacity for species identification will improve mosquito characterization, leading to accurate control strategies/treatment targeting specific mosquito species, and thus will contribute to malaria eradication. It is crucial for disease epidemiology and surveillance to accurately identify the Plasmodium spp. that are causing malaria in patients. The capacity for disease surveillance will be significantly increased by the development of more accurate, precise, automated, and high-throughput diagnostic techniques. In conclusion, although morphological identification is quick and achievable at a reduced cost, molecular identification is preferred for specificity and sensitivity. To achieve the targeted malaria elimination goal, proper identification of vectors using accurate techniques for effective control measures should be prioritized.

Genetic Diversity and Population Structure of Anopheles funestus in Western Kenya Based on Mitochondrial DNA Marker COII

The mitochondrial marker, COII, was employed to assess the genetic structure and diversity of Anopheles funestus, a very important malaria vector in Africa that adapt and colonize different ecological niches in western Kenya. Mosquitoes were collected using mechanical aspirators in four areas (Bungoma, Port Victoria, Kombewa, and Migori) in western Kenya. Following morphological identification, PCR was used to confirm the species. The COII gene was amplified, sequenced, and analyzed to determine genetic diversity and population structure. A total of 126 (Port Victoria-38, Migori-38, Bungoma-22, and Kombewa-28) sequences of COII were used for population genetic analysis. Anopheles funestus had a high haplotype diversity (Hd = 0.97 to 0.98) but low nucleotide diversity (Π = 0.004 to 0.005). The neutrality test revealed negative Tajima's D and Fs values indicating an excess of low-frequency variation. This could be attributed to either population expansion or negative selection pressure across all the populations. No genetic or structural differentiation (Fst = -0.01) and a high level of gene flow (Gamma St, Nm = 17.99 to 35.22) were observed among the populations. Population expansion suggests the high adaptability of this species to various ecological requirements, hence sustaining its vectorial capacity and malaria transmission.

Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border

BACKGROUND

Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area.

METHODS

1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F_ST = 0.97824, 0.97343, P < 0.05), as well as lower gene flow (Nm = 0.01112, 0.01365) in the An. nivipes populations. Phylogenetic analysis and SAMOVA revealed a gene barrier restricting gene flow among three An. nivipes populations. Mantel test suggested a significant correlation between geography and gene distance in all An. nivipes populations (Z = 44,983.1865, r = 0.5575, P = 0.0070). Neutrality test and Mismatch distribution revealed a recent population expansion of An. nivipes in the Cambodia-Laos population.

CONCLUSIONS

Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. Based on sequences analysis of COII, a genetic barrier between Cambodia-Laos and two Indian populations was found, and a recent population expanding or selecting of An. nivipes occurred in the Cambodia-Laos population, suggesting that COII might be an effective marker for describing the molecular phylogeny of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.

New assessment of Anopheles vector species identification using MALDI-TOF MS

Background

Anopheles species identification is essential for an effective malaria vector control programme. Matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry (MS) has been developed to identify adult Anopheles species, using the legs or the cephalothorax. The protein repertoire from arthropods can vary according to compartment, but there is no general consensus regarding the anatomic part to be used.

Methods

To determine the body part of the Anopheles mosquitoes best suited for the identification of field specimens, a mass spectral library was generated with head, thorax with wings and legs of Anopheles gambiae, Anopheles arabiensis and Anopheles funestus obtained from reference centres. The MSL was evaluated using two independent panels of 52 and 40 An. gambiae field-collected in Mali and Guinea, respectively. Geographic variability was also tested using the panel from Mali and several databases containing added specimens from Mali and Senegal.

Results

Using the head and a database without specimens from the same field collection, the proportion of interpretable and correct identifications was significantly higher than using the other body parts at a threshold value of 1.7 (p < 0.0001). The thorax of engorged specimens was negatively impacted by the blood meal after frozen storage. The addition of specimens from Mali into the database significantly improved the results of Mali panel (p < 0.0001), which became comparable between head and legs. With higher identification scores, the using of the head will allow to decrease the number of technical replicates of protein extract per specimen, which represents a significant improvement for routine use of MALDI-TOF MS.

Conclusions

The using of the head of Anopheles may improve the performance of MALDI-TOF MS. Region-specific mass spectrum databases will have to be produced. Further research is needed to improve the standardization in order to share online spectral databases.

Species diversity and insecticide resistance within the Anopheles hyrcanus group in Ubon Ratchathani Province, Thailand

Background

Members of the Anopheles hyrcanus group have been incriminated as important malaria vectors. This study aims to identify the species and explore the insecticide susceptibility profile within the Anopheles hyrcanus group in Ubon Ratchathani Province, northeastern Thailand where increasing numbers of malaria cases were reported in 2014.

Methods

Between 2013 and 2015, five rounds of mosquito collections were conducted using human landing and cattle bait techniques during both the rainy and dry seasons. Anopheles mosquitoes were morphologically identified and their insecticide susceptibility status was investigated. Synergist bioassays were carried out with An. hyrcanus (s.l.) due to their resistance to all insecticides. An ITS2-PCR assay was conducted to identify to species the Hyrcanus group specimens.

Results

Out of 10,361 Anopheles females collected, representing 18 taxa in 2 subgenera, 71.8% were morphologically identified as belonging to the Hyrcanus Group (subgenus Anopheles), followed by An. barbirostris group (7.9%), An. nivipes (6.5%), An. philippinensis (5.9%) and the other 14 Anopheles species. Specimens of the Hyrcanus Group were more prevalent during the rainy season and were found to be highly zoophilic. Anopheles hyrcanus (s.l.) was active throughout the night, with an early peak of activity between 18:00 h and 21:00 h. ITS2-PCR assay conducted on 603 DNA samples from specimens within the Hyrcanus Group showed the presence of five sisters species. Anopheles peditaeniatus was the most abundant species (90.5%, n = 546), followed by An. nitidus (4.5%, n = 27), An. nigerrimus (4.3%, n = 26), An. argyropus (0.5%, n = 3), and An. sinensis (0.2%, n = 1). All An. hyrcanus (s.l.) specimens that were found resistant to insecticides (deltamethrin 0.05%, permethrin 0.75% and DDT 4% and synergist tests) belonged to An. peditaeniatus. The degree of resistance in An. peditaeniatus to each of these three insecticides was approximately 50%. Addition of PBO (Piperonyl butoxide), but not DEF (S.S.S-tributyl phosphotritioate), seemed to restore susceptibility, indicating a potential role of oxidases as a detoxifying enzyme resistance mechanism.

Conclusions

A better understanding of mosquito diversity related to host preference, biting activity and insecticide resistance status will facilitate the implementation of locally adapted vector control strategies.

Screening of mosquitoes for filarioid helminths in urban areas in south western Poland-common patterns in European Setaria tundra xenomonitoring studies

In recent years, numerous studies screening mosquitoes for filarioid helminths (xenomonitoring) have been performed in Europe. The entomological monitoring of filarial nematode infections in mosquitoes by molecular xenomonitoring might serve as the measure of the rate at which humans and animals expose mosquitoes to microfilariae and the rate at which animals and humans are exposed to the bites of the infected mosquitoes. We hypothesized that combining the data obtained from molecular xenomonitoring and phenological studies of mosquitoes in the urban environment would provide insights into the transmission risk of filarial diseases. In our search for Dirofilaria spp.-infected mosquitoes, we have found Setaria tundra-infected ones instead, as in many other European studies. We have observed that cross-reactivity in PCR assays for Dirofilaria repens, Dirofilaria immitis, and S. tundra COI gene detection was the rule rather than the exception. S. tundra infections were mainly found in Aedes mosquitoes. The differences in the diurnal rhythm of Aedes and Culex mosquitoes did not seem a likely explanation for the lack of S. tundra infections in Culex mosquitoes. The similarity of S. tundra COI gene sequences found in Aedes vexans and Aedes caspius mosquitoes and in roe deer in many European studies, supported by data on Ae. vexans biology, suggested host preference as the most likely cause of the mosquito genus-biased infections. High diversity of the COI gene sequences isolated in the city of Wroclaw in south western Poland and the presence of identical or almost identical sequences in mosquitoes and roe deer across Europe suggests that S. tundra has been established in most of Europe for a very long time.

Mitochondrial variation in subpopulations of Anopheles balabacensis Baisas in Sabah, Malaysia (Diptera: Culicidae)

Anopheles balabacensis, the primary vector of Plasmodium knowlesi in Sabah, Malaysia, is both zoophilic and anthropophilic, feeding on macaques as well as humans. It is the dominant Anopheles species found in Kudat Division where it is responsible for all the cases of P. knowlesi. However there is a paucity of basic biological and ecological information on this vector. We investigated the genetic variation of this species using the sequences of cox1 (1,383 bp) and cox2 (685 bp) to gain an insight into the population genetics and inter-population gene flow in Sabah. A total of 71 An. balabacensis were collected from seven districts constituting 14 subpopulations. A total of 17, 10 and 25 haplotypes were detected in the subpopulations respectively using the cox1, cox2 and the combined sequence. Some of the haplotypes were common among the subpopulations due to gene flow occurring between them. AMOVA showed that the genetic variation was high within subpopulations as compared to between subpopulations. Mantel test results showed that the variation between subpopulations was not due to the geographical distance between them. Furthermore, Tajima's D and Fu's Fs tests showed that An. balabacensis in Sabah is experiencing population expansion and growth. High gene flow between the subpopulations was indicated by the low genetic distance and high gene diversity in the cox1, cox2 and the combined sequence. However the population at Lipasu Lama appeared to be isolated possibly due to its higher altitude at 873 m above sea level.

Medical Entomology: A Reemerging Field of Research to Better Understand Vector-Borne Infectious Diseases

In the last decade, the Chikungunya and Zika virus outbreaks have turned public attention to the possibility of the expansion of vector-borne infectious diseases worldwide. Medical entomology is focused on the study of arthropods involved in human health. We review here some of the research approaches taken by the medical entomology team of the University Hospital Institute (UHI) Méditerranée Infection of Marseille, France, with the support of recent or representative studies. We propose our approaches to technical innovations in arthropod identification and the detection of microorganisms in arthropods, the use of arthropods as epidemiological or diagnostic tools, entomological investigations around clinical cases or within specific populations, and how we have developed experimental models to decipher the interactions between arthropods, microorganisms, and humans.

Phylogenetic Relationships Among Malaria Vectors and Closely Related Species in Thailand Using Multilocus DNA Sequences

The evolutionary and taxonomic status is important for understanding speciation events and phylogenetic relationships between closely related vector and nonvector species. This information is useful for targeting important disease vector species groups for the development of novel genetic-based vector and pathogen control methods. In this study, different phylogenetic analyses were performed to reconstruct phylogenetic trees for the primary malaria vectors in Thailand based on sequence information of 4 DNA fragments from the nuclear and mitochondrial regions. The primary Anopheles species in the subgenus Cellia involved in malaria transmission in Thailand separate clearly into 3 distinct clades: the Leucosphyrus group, Minimus subgroup, and Maculatus group. The phylogenetic trees based on different reconstructed algorithms and different gene regions provided congruent phylogenetic status of the mosquito species studied. The phylogenetic relationships of malaria vector species examined followed similar patterns based on morphological characters. An estimate of the divergence time among the Anopheles species infers that they were present during the Eocene and Miocene periods (>41 million years ago). Congruent phylogenetic analysis of malaria vectors is presented with different algorithms and gene regions. The nuclear TOLL6 fragment appears useful for molecular phylogenetic, species DNA barcode, and Anopheles population genetic analyses.

Spatial distributions of Anopheles species in relation to malaria incidence at 70 localities in the highly endemic Northwest and South Pacific coast regions of Colombia

BACKGROUND

A proper identification of malaria vectors is essential for any attempt to control this disease. Between 40 and 47 Anopheles species have been recorded in Colombia, and eight species complexes have been identified in the last decade. An update of Anopheles species distribution and its relationship with malaria is required, particularly for newly identified members of species complexes.

METHODS

A cross-sectional entomological study was conducted at 70 localities in the highest malaria transmission areas in Colombia. In each locality, immature and adult mosquitoes were collected. All specimens were determined using morphological characters and confirmed used restriction profiles of Internal Transcribed Spacer 2 (PCR-RFLP-ITS2), and Cytochrome c Oxidase I (COI) sequence gene. To detect natural Plasmodium infections, enzyme-linked immunosorbent assay and nested PCR analysis were used. Distribution of Anopheles species was spatially associated with malaria incidence.

RESULTS

A total of 1736 larvae and 12,052 adult mosquitoes were determined in the 70 localities. Thirteen Anopheles species were identified. COI sequence analysis suggested 4 new lineages for Colombia: for Anopheles albimanus (An. albimanus B), Anopheles pseudopunctipennis s.l., Anopheles neivai (An. neivai nr. neivai 4), and Anopheles apicimacula. Two members of species complexes were identified, as: Anopheles nuneztovari C, and Anopheles albitarsis I. Another seven species were confirmed. Four mosquitoes were infected with Plasmodium species, An. albimanus B and An. nuneztovari C. In Northwest of Colombia, An. nuneztovari C, An. albimanus, and Anopheles darlingi were present in the municipalities with highest annual parasitic index (API) (>35 cases/1000 inhabitants). In the north of South Pacific coast, with a similar API, An. nuneztovari C were widely distributed inland, and the main species in coastal regions were An. albimanus B and An. neivai s.l. In the South Pacific coast bordering with Ecuador, 3 Anopheles species were found in municipalities with high API (15-88 cases/1000 inhabitants): An. albimanus B, Anopheles calderoni and An. neivai s.l.

CONCLUSIONS

In the highest malaria areas of Colombia, 13 Anopheles species and four new lineages were found, which highlights the need for updating the species distribution. A DNA barcode analysis allowed the taxonomic identification to be refined, particularly for species complexes, and to improve the further understanding of their relation with malaria transmission.

Emerging tools for identification of arthropod vectors

The rapid and reliable identification of arthropod vector species is an essential component of the fight against vector-borne diseases. However, owing to the lack of entomological expertise required for the morphological identification method, development of alternative and complementary tools is needed. This review describes the main methods used for arthropod identification, focusing on the emergence of protein profiling using MALDI-TOF MS technology. Sample preparation, analysis of reproducibility, database creation and blind tests for controlling accuracy of this tool for arthropod identification are described. The advantages and limitations of the MALDI-TOF MS method are illustrated by emphasizing different hematophagous arthropods, including mosquitoes and ticks, the top two main vectors of infectious diseases.

Population Genetic Structure of the Malaria Vector Anopheles sinensis (Diptera: Culicidae) Sensu Stricto and Evidence for Possible Introgression in the Republic of Korea

Anopheles sinensis Wiedemann sensu stricto (s.s.) is a dominant mosquito and considered a secondary malaria vector in the Republic of Korea (ROK). Despite the potential significance for malaria control, population genetics studies have been conducted using only mitochondrial DNA (mtDNA), and studies of the genetics of hybridization have never been attempted. In this study, 346 specimens from 23 localities were subject to experiments. Among them, 305 An. sinensis s.s. specimens from 20 localities were used for mtDNA analysis, and 346 specimens comprising 341 An. sinensis s.s. from 22 localities and five Anopheles kleini Rueda from one locality were examined in the microsatellite study. Neighbor-joining analysis of pairwise FST and RST based on microsatellite results showed that the populations are divided into two groups, as did the mtDNA results. However, the Bayesian analysis and factorial correspondence analysis plots showed three distinct clusters. Among the mtDNA and microsatellite results, only microsatellites represented small but positive and significant isolation-by-distance patterns. Both molecular markers show the Taebaek and Sobaek Mountain ranges as barriers between the northern and southern parts of the ROK. The newly recognized third group suggests possible introgressive hybridization of An. sinensis s.s. with closely related species. The slightly different composition of populations in each group based on different markers is probably because of different population dynamics in each group. These results imply that there is restricted gene flow of epidemiologically important malaria-related genes between the northern and southern parts of the ROK.

Mitochondrial DNA-based genetic diversity of Anopheles nivipes in North East India

Malaria is a major public health problem in north-east India mainly transmitted by Anopheles baimai and An. minimus while Anopheles nivipes plays an important supportive role. The genetic diversity of An. nivipes in north-east India was investigated by employing two mitochondrial DNA markers namely NADH dehydrogenase 5 (ND5) and cytochrome oxidase sub unit II (COII). High genetic diversity in An. nivipes was observed by the detection of 16 haplotypes among 30 sequences of ND5 gene and 29 haplotypes among 29 COII sequences. Anopheles nivipes of north-east India was significantly differentiated genetically with that of neighboring South-east Asian An. nivipes as revealed by pairwise F_ST values of 0.127 (p < 0.01) and 0.044 (p < 0.001) for ND5 and COII genes, respectively, suggesting geographical barriers to gene flow in this species between the two geographical areas resulting in significant population structuring.

Multilocus nuclear DNA markers reveal population structure and demography of Anopheles minimus

Utilization of multiple putatively neutral DNA markers for inferring evolutionary history of species population is considered to be the most robust approach. Molecular population genetic studies have been conducted in many species of Anopheles genus, but studies based on single nucleotide polymorphism (SNP) data are still very scarce. Anopheles minimus is one of the principal malaria vectors of Southeast (SE) Asia including the Northeastern (NE) India. Although population genetic studies with mitochondrial genetic variation data have been utilized to infer phylogeography of the SE Asian populations of this species, limited information on the population structure and demography of Indian An. minimus is available. We herewith have developed multilocus nuclear genetic approach with SNP markers located in X chromosome of An. minimus in eight Indian and two SE Asian population samples (121 individual mosquitoes in total) to infer population history and test several hypotheses on the phylogeography of this species. While the Thai population sample of An. minimus presented the highest nucleotide diversity, majority of the Indian samples were also fairly diverse. In general, An. minimus populations were moderately substructured in the distribution range covering SE Asia and NE India, largely falling under three distinct genetic clusters. Moreover, demographic expansion events could be detected in the majority of the presently studied populations of An. minimus. Additional DNA sequencing of the mitochondrial COII region in a subset of the samples (40 individual mosquitoes) corroborated the existing hypothesis of Indian An. minimus falling under the earlier reported mitochondrial lineage B.

Vectorial capacity and genetic diversity of Anopheles annularis (Diptera: Culicidae) mosquitoes in Odisha, India from 2009 to 2011

Anopheles annularis is one of the major vectors of malaria in Odisha, India. The present study was undertaken to determine the vectorial capacity and assess the genetic diversity of An. annularis collected from different endemic regions of Odisha. Mosquitoes were collected from thirteen endemic districts using standard entomological collection methods from 2009 to 2011. Sibling species of An. annularis were identified by PCR-RFLP and sequencing of D3 region of 28S ribosomal DNA (rDNA) region. Plasmodium falciparum (Pf) sporozoite rate and human blood fed percentage (HBF) were estimated by multiplex PCR using Pf and human specific primers. Genetic diversity of An. annularis was estimated by ISSR markers. Out of 1647 An. annularis collected, 1353 (82.15%) were collected by mechanical aspirators and 294 (17.85%) by light trap. 49 (2.97%) were positive for human blood and 18 (1.09%) were positive for Pf sporozoite. PCR-RFLP and sequencing analyses detected only An annularis A in the study areas. Overall genetic differentiation among An. annularis populations was moderate (FST=0.048) and showed significant correlation between genetic distance and geographic distance (r=0.882; P<0.05). Angul population proved to be genetically unique and was highly divergent FST>0.110) from other populations, suggesting low gene flow between them. The study indicated that only An. annularis A was found in Odisha with potential vectorial capacity that can play a major role in malaria transmission. ISSR markers proved to be useful molecular tools to evaluate genetic variability in An. annularis populations.

Mitochondrial DNA analyses reveal low genetic diversity in Culex quinquefasciatus from residential areas in Malaysia

The present study explored the intraspecific genetic diversity, dispersal patterns and phylogeographic relationships of Culex quinquefasciatus Say (Diptera: Culicidae) in Malaysia using reference data available in GenBank in order to reveal this species' phylogenetic relationships. A statistical parsimony network of 70 taxa aligned as 624 characters of the cytochrome c oxidase subunit I (COI) gene and 685 characters of the cytochrome c oxidase subunit II (COII) gene revealed three haplotypes (A1-A3) and four haplotypes (B1-B4), respectively. The concatenated sequences of both COI and COII genes with a total of 1309 characters revealed seven haplotypes (AB1-AB7). Analysis using tcs indicated that haplotype AB1 was the common ancestor and the most widespread haplotype in Malaysia. The genetic distance based on concatenated sequences of both COI and COII genes ranged from 0.00076 to 0.00229. Sequence alignment of Cx. quinquefasciatus from Malaysia and other countries revealed four haplotypes (AA1-AA4) by the COI gene and nine haplotypes (BB1-BB9) by the COII gene. Phylogenetic analyses demonstrated that Malaysian Cx. quinquefasciatus share the same genetic lineage as East African and Asian Cx. quinquefasciatus. This study has inferred the genetic lineages, dispersal patterns and hypothetical ancestral genotypes of Cx. quinquefasciatus.

Towards mosquito sterile insect technique programmes: exploring genetic, molecular, mechanical and behavioural methods of sex separation in mosquitoes

When considering a mosquito release programme, one of the first issues to be addressed is how to eliminate/separate the females. The greatest number of options might eventually be available for those who can use transgenic mosquitoes, but the inherent characteristics of the target species may also provide possibilities for interim measures until more efficient methods can be developed. Differences in intrinsic size, in behaviour and in development rate between females and males are often available and useful for sexing. Efficient species-specific systems for eliminating females at the embryo stage have been developed, but most have since been discarded due to lack of use. Ideal systems specifically kill female embryos using some treatment that can be manipulated during production. Such killing systems are far more efficient than using intrinsic sexual differences, but they systems require selectable genetic markers and sex-linkage created by rare random chromosomal rearrangements. While intrinsic sexual differences should not be considered as long-term candidates for the development of robust and efficient sexing approaches, in the absence of these, the accessibility and integration of less efficient systems can provide a stop-gap measure that allows rapid start up with a minimum of investment. The International Atomic Energy Agency is funding over a 5 year period (2013-2018) a new Coordinated Research Project on "Exploring Genetic, Molecular, Mechanical and Behavioural Methods of Sex Separation in Mosquitoes" to network researchers and to address the critical need of genetic sexing strains for the implementation of the sterile insect technique (using radiation-sterilised or transgenic male mosquitoes) and for insect incompatibility technique programmes against disease-transmitting mosquitoes.

Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis

BACKGROUND

Anopheles (Kerteszia) cruzii (Diptera: Culicidae) is a primary vector of human and simian malaria parasites in southern and southeastern Brazil. Earlier studies using chromosome inversions, isoenzymes and a number of molecular markers have suggested that An. cruzii is a species complex.

RESULTS

In this study, a multilocus approach using six loci, three circadian clock genes and three encoding ribosomal proteins, was carried out to investigate in more detail the genetic differentiation between the An. cruzii populations from Florianópolis-Santa Catarina (southern Brazil) and Itatiaia-Rio de Janeiro States (southeastern Brazil). The analyses were performed first comparing Florianópolis and Itatiaia, and then comparing the two putative sympatric incipient species from Itatiaia (Itatiaia A and Itatiaia B). The analysis revealed high FST values between Florianópolis and Itatiaia (considering Itatiaia A and B together) and also between the sympatric Itatiaia A and Itatiaia B, irrespective of their function. Also, using the IM program, no strong indication of migration was found between Florianópolis and Itatiaia (considering Itatiaia A and B together) using all loci together, but between Itatiaia A and Itatiaia B, the results show evidence of migration only in the direction of Itatiaia B.

CONCLUSIONS

The results of the multilocus analysis indicate that Florianópolis and Itatiaia represent different species of the An. cruzii complex that diverged around 0.6 Mya, and also that the Itatiaia sample is composed of two sympatric incipient species A and B, which diverged around 0.2 Mya. Asymmetric introgression was found between the latter two species despite strong divergence in some loci.

Molecular phylogenetic analysis of Anopheles and Cellia subgenus anophelines (Diptera: Culicidae) in temperate and tropical regions of Iran

Molecular studies on population genetics of speciation across Iran have recently started. Morphological and molecular studies have showed that 25 species of genus Anopheles are present in the country; however, relationships between vector and non-vector species as well as compatibility of morphological characters with molecular data have not been verified. Molecular phylogenetic analysis was undertaken on the Anopheles and Cellia subgenus members internal transcribed spacer 2 (ITS2) sequences submitted to GenBank among the Oriental and Palearctic members in north and southern Iran. rDNA-ITS2 sequences were extracted from the GenBank and analyzed using bioinformatics softwares: BLAST, ITS2 annotation tool (version 3.0.13), ClustalW, and MEGA5 in neighbor-joining and maximum likelihood algorithms. There are not any submitted sequences in GenBank from Iran for the following seven species: Anopheles algeriensis, Anopheles marteri, Anopheles plumbeus, Anopheles peditaeniatus, Anopheles melanoon, Anopheles subpictus, and Anopheles mongolensis; therefore, they have not been included in the study. Although these molecular-based phylogenetic trees match well enough with classical morphological taxonomy, the arrangement of species did not match with morphological classification in some cases. Correct species identification is essential for control of vector born disease such as malaria; therefore, phylogenetic methods will help to understand the relationship among the members of the target species within the genus Anopheles. It could also help us to design molecular markers for species differentiation particularly in cryptic species, which is difficult to classify them based on morphological features.

Inter and intra-specific diversity of parasites that cause lymphatic filariasis

Lymphatic filariasis is caused by three closely related nematode parasites: Wuchereria bancrofti, Brugia malayi and Brugia timori. These species have many ecological variants that differ in several aspects of their biology such as mosquito vector species, host range, periodicity, and morphology. Although the genome of B. malayi (the first genome sequenced from a parasitic nematode) has been available for more than five years, very little is known about genetic variability among the lymphatic dwelling filariae. The genetic diversity among these worms is not only interesting from a biological perspective, but it may have important practical implications for the Global Program to Eliminate Lymphatic Filariasis, as the parasites may respond differently to diagnostic tests and/or medical interventions. Therefore, better information on their genetic variability is urgently needed. With improved methods for nucleic acid extraction and recent advances in sequencing chemistry and instrumentation, this gap can be filled relatively inexpensively. Improved information on filarial genetic diversity may increase the chances of success for lymphatic filariasis elimination programs.

Identification of Anopheles daciae in Germany through ITS2 sequencing

Until the middle of the twentieth century, malaria was frequently endemic in parts of Germany; Anopheles maculipennis complex species were considered the primary vectors. Three species of this complex have been identified in Germany: A. maculipennis s.s., Anopheles messeae and Anopheles atroparvus; the last predominantly from the coastal regions of Northern Germany. Anopheles daciae is a recently described member of the A . maculipennis complex and resembles the well-characterised species A. messeae, although the two species can be distinguished through their egg morphology and sequencing of the internal transcribed spacer 2 (ITS2) region of their nuclear rDNA. In this study, we harvested larval and adult mosquito samples from five breeding sites and ten CO(2) trap collection sites in the Upper Rhine Valley of Southwestern Germany to analyse the complement of anopheline species present. Mosquito ITS2 DNA was extracted and polymerase chain reaction (PCR)-amplified using established protocols. Genomic analysis was performed by a species-diagnostic restriction fragment length polymorphism assay as well as by sequencing of PCR products; the data obtained were aligned against nucleic acid sequences from English mosquitoes retrieved from GenBank. Additionally, the larval breeding sites of A. messeae were characterised through water quality measurement. Forty-seven samples were successfully processed, of which 6 were identified as A. daciae and 41 as A. messeae. All samples of A. daciae, which has not previously been found in Central Europe, originated from one CO(2) trap collection site in Dettenheim, close to Karlsruhe, Southwestern Germany. The identification of this malarial vector in a novel area may have implications for the re-emergence of disease subsequent to climatic changes.

Resolving genetic diversity in Australasian Culex mosquitoes: incongruence between the mitochondrial cytochrome c oxidase I and nuclear acetylcholine esterase 2

Insects that vector pathogens are under constant surveillance in Australasia although the repertoire of genetic markers to distinguish what are often cryptic mosquito species remains limited. We present a comparative assessment of the second exon-intron region of the acetylcholine esterase 2 gene (ace-2) and the mitochondrial DNA cytochrome c oxidase I (COI) using two closely related Australasia mosquitoes Culex annulirostris and Culex palpalis. The COI revealed eight divergent lineages of which four were confirmed with the ace-2. We dissect out the nuclear chromosomal haplotypes of the ace-2 as well as the exon-intron regions by assessing the protein's tertiary structure to reveal a hypervariable 5'-exon that forms part of an external protein loop and displays a higher polymorphic rate than the intron. We retrace the evolutionary history of these mosquitoes by phylogenetic inference and by testing different evolutionary hypotheses. We conclude that DNA barcoding using COI may overestimate the diversity of Culex mosquitoes in Australasia and should be applied cautiously with support from the nuclear DNA such as the ace-2. Together the COI and ace-2 provide robust evidence for distinct cryptic Culex lineages--one of which correlates exactly with the southern limit of Japanese encephalitis virus activity in Australasia.

The type and mysorensis forms of the Anopheles stephensi (Diptera: Culicidae) in India exhibit identical ribosomal DNA ITS2 and domain-3 sequences

Anopheles (Cellia) stephensi Liston 1901 is one of the major malaria vectors in the Indian subcontinent, Iran, and the Middle East. Three races in this species, namely A. stephensi stephensi (type form), A. stephensi variety mysorensis, and A. stephensi intermediate form, have earlier been reported by several investigators. We describe here the sequencing of the ribosomal DNA internal transcribed spacer 2 (ITS2) and domain-3 (D3) loci of the A. stephensi type and variety mysorensis forms. We also sequenced field-collected adult specimens of this mosquito from three different regions of India. Both forms of A. stephensi showed identical ITS2 and D3 sequences. We did not find any intraspecies sequence variation among the 70 specimens sequenced in this study. In contrast to the eight ITS2 haplotypes observed among Iranian A. stephensi population, we found only one ITS2 haplotype in India. This is the first time to our knowledge that the sequence of the D3 locus of A. stephensi is being reported here. In conclusion, the type and variety mysorensis forms of A. stephensi exhibit identical nucleotide sequences at their ITS2 and D3 loci.

Cytogenetic and molecular evidence for two species in the Anopheles barbirostris complex (Diptera: Culicidae) in Thailand

Seventeen isolines of Anopheles barbirostris derived from animal-biting female mosquitoes showed three karyotypic forms: Form A (X2, Y1) in five isolines from Phetchaburi province; Form B (X1, X3, Y2) in three and eight isolines from Chiang Mai and Ubon Ratchathani provinces, respectively; Form C (X2, Y3) in one isoline from Phetchaburi province. All 17 isolines exhibited an average branch summation of seta 2-VI pupal skins ranging from 12.1-13.0 branches, which was in the limit of A. barbirostris (6-18 branches). Of the 12 human-biting isolines from Chiang Mai province, five isolines showed Form B (X2, Y2), and seven isolines exhibited a new karyotypic form designated as Form E (X2, Y5). All of 12 isolines had an average branch summation of seta 2-VI pupal skins ranging from 22.4-24.5 branches, which was in the limit of Anopheles campestris (17-58 branches). Thus, they were tentatively designated as A. campestris-like Forms B and E. Hybridization between A. campestris-like Forms B and E showed that they were genetically compatible, yielding viable progeny for several generations suggesting conspecific relationships of these two karyotypic forms. Reproductive isolation among crosses between A. campestris-like Form B and A. barbirostris Forms A, B, and C strongly suggested the existence of these two species. In addition, the very low intraspecific variation (genetic distance<0.005) of the nucleotide sequence of ITS2 of the rDNA and COI and COII of mitochondrial DNA of the seven isolines of A. campestris-like Forms B and E supported their conspecific relationship. The large sequence divergence of ITS2 (0.203-0.268), COI (0.026-0.032), and COII (0.030-0.038) from genomic DNA of A. campestris-like Forms B and E and the A. barbirostris Forms A, B, and C clearly supported cytogenetic and morphological evidence.

The second internal transcribed spacer of nuclear ribosomal DNA as a tool for Latin American anopheline taxonomy - a critical review

Among the molecular markers commonly used for mosquito taxonomy, the internal transcribed spacer 2 (ITS2) of the ribosomal DNA is useful for distinguishing among closely-related species. Here we review 178 GenBank accession numbers matching ITS2 sequences of Latin American anophelines. Among those, we found 105 unique sequences corresponding to 35 species. Overall the ITS2 sequences distinguish anopheline species, however, information on intraspecific and geographic variations is scarce. Intraspecific variations ranged from 0.2% to 19% and our analysis indicates that misidentification and/or sequencing errors could be responsible for some of the high values of divergence. Research in Latin American malaria vector taxonomy profited from molecular data provided by single or few field capture mosquitoes. However we propose that caution should be taken and minimum requirements considered in the design of additional studies. Future studies in this field should consider that: (1) voucher specimens, assigned to the DNA sequences, need to be deposited in collections, (2) intraspecific variations should be thoroughly evaluated, (3) ITS2 and other molecular markers, considered as a group, will provide more reliable information, (4) biological data about vector populations are missing and should be prioritized, (5) the molecular markers are most powerful when coupled with traditional taxonomic tools.

Identification of two cryptic species in the Anopheles (Cellia) annularis complex using ribosomal DNA PCR-RFLP

Anopheles (Cellia) annularis Van der Wulp is a complex of two sibling species provisionally designated as species A and B and can only be differentiated on the basis of the paracentric inversion in the ovarian polytene chromosomes. To analyze the distribution of these two species and to develop a molecular method for the identification of these two cryptic species, we sequenced the ribosomal DNA internal transcribed spacer 2 (ITS2) and domain 3 (D3) of A. annularis specimens collected from Sonapur (Assam), Jabalpur (Madhya Pradesh), Ranchi (Jharkhand), and Ghaziabad (Uttar Pradesh). We did not find any sequence variation among the specimens collected from Assam, Madhya Pradesh, and Jharkhand states, whereas two types of sequences were obtained from the specimens collected from the state of Uttar Pradesh, which correspond to species A and B of the A. annularis complex. Species A was more prevalent among the all four regions studied. The ITS2 sequence of species A showed unique restriction sites for MvaI and Eco24I, while species B displayed HinfI and NruI sites. Similarly, the D3 sequence of species A showed unique restriction site for Alw26I, while species B showed a unique KpnI site. In this study, we report for the first time the development of ribosomal DNA polymerase chain reaction-restriction fragment length polymorphism methods for identifying these two cryptic species of the Annularis complex.

PCR-RFLP method for the identification of four members of the Anopheles annularis group of mosquitoes (Diptera: Culicidae)

The Anopheles annularis group of mosquitoes is widely distributed in Southeast Asia and may be locally important as malaria vectors. Members of this group are morphologically very similar and often difficult to distinguish, particularly A. nivipes and A. philippinensis. We report the sequence analysis of the rDNA internal transcribed spacer 2 (ITS2) and Domain-3 (D3) regions of the four members of the A. annularis group -A. nivipes, A. philippinensis, A. annularis and A. pallidus - and a method for their molecular identification. No intraspecies sequence variation was detected among the specimens, while interspecific sequence differences were greater for ITS2 than the D3 regions. Comparison of the D3 sequences of the four species revealed two SmaI restriction sites in A. nivipes, but only one site in A. philippinensis, A. annularis and A. pallidus. The ApaI site was present in both A. philippinensis and A. pallidus, while an NcoI site was present in A. pallidus only. Restriction digestion of the PCR products of D3 fragment individually with SmaI, ApaI and NcoI produced a distinctive pattern for all the four species. We present, for the first time, a PCR-RFLP method to distinguish the four members of the A. annularis group of mosquitoes.

A Novel Real-Time PCR Assay for the Speciation of Medically Important Ticks

The identification of ticks using morphological characters is a well-established practice, however specimens that are small or damaged are often difficult to speciate. A novel, rapid real-time PCR assay, which targets the second internal transcribed spacer (ITS2) region in the nuclear ribosomal DNA gene, was developed for identification of four tick species of utmost medical importance in the United States: Ixodes scapularis, I. pacificus, Dermacentor variabilis, and Amblyomma americanum. Computational analyses of public databases and DNA sequencing studies revealed regions that could be specifically targeted with oligonucleotides optimized for TaqMan chemistry. The oligonucleotide sets designed in this study are specific at both the genus and species levels, and are sensitive at 0.1-1 pg of total tick DNA.