Detection of hybrids in natural populations of theAnopheles gambiaecomplex by the rDNA-based, PCR method

Efficient Bayesian inference under the multispecies coalescent with migration

Analyses of genome sequence data have revealed pervasive interspecific gene flow and enriched our understanding of the role of gene flow in speciation and adaptation. Inference of gene flow using genomic data requires powerful statistical methods. Yet current likelihood-based methods involve heavy computation and are feasible for small datasets only. Here, we implement the multispecies-coalescent-with-migration model in the Bayesian program bpp, which can be used to test for gene flow and estimate migration rates, as well as species divergence times and population sizes. We develop Markov chain Monte Carlo algorithms for efficient sampling from the posterior, enabling the analysis of genome-scale datasets with thousands of loci. Implementation of both introgression and migration models in the same program allows us to test whether gene flow occurred continuously over time or in pulses. Analyses of genomic data from Anopheles mosquitoes demonstrate rich information in typical genomic datasets about the mode and rate of gene flow.

Gene drive in species complexes: defining target organisms

Engineered gene drives, which bias their own inheritance to increase in frequency in target populations, are being developed to control mosquito malaria vectors. Such mosquitoes can belong to complexes of both vector and nonvector species that can produce fertile interspecific hybrids, making vertical gene drive transfer (VGDT) to sibling species biologically plausible. While VGDT to other vectors could positively impact human health protection goals, VGDT to nonvectors might challenge biodiversity ones. Therefore, environmental risk assessment of gene drive use in species complexes invites more nuanced considerations of target organisms and nontarget organisms than for transgenes not intended to increase in frequency in target populations. Incorporating the concept of target species complexes offers more flexibility when assessing potential impacts from VGDT.

Evidence for natural hybridization and novel Wolbachia strain superinfections in the Anopheles gambiae complex from Guinea

Wolbachia, a widespread bacterium which can influence mosquito-borne pathogen transmission, has recently been detected within Anopheles (An.) species that are malaria vectors in Sub-Saharan Africa. Although studies have reported Wolbachia strains in the An. gambiae complex, apparent low density and prevalence rates require confirmation. In this study, wild Anopheles mosquitoes collected from two regions of Guinea were investigated. In contrast with previous studies, RNA was extracted from adult females (n = 516) to increase the chances for the detection of actively expressed Wolbachia genes, determine Wolbachia prevalence rates and estimate relative strain densities. Molecular confirmation of mosquito species and Wolbachia multilocus sequence typing (MLST) were carried out to analyse phylogenetic relationships of mosquito hosts and newly discovered Wolbachia strains. Strains were detected in An. melas (prevalence rate of 11.6%-16/138) and hybrids between An. melas and An. gambiae sensu stricto (prevalence rate of 40.0%-6/15) from Senguelen in the Maferinyah region. Furthermore, a novel high-density strain, termed wAnsX, was found in an unclassified Anopheles species. The discovery of novel Wolbachia strains (particularly in members, and hybrids, of the An. gambiae complex) provides further candidate strains that could be used for future Wolbachia-based malaria biocontrol strategies.

Is Genetic Continuity Between Anopheles sinensis (Diptera: Culicidae) and its Sibling Species Due to Gene Introgression or Incomplete Speciation?

The Anopheles mosquito Hyrcanus Group is widely distributed geographically across both Palearctic and Oriental regions and comprises 26 valid species. Although the species Anopheles sinensis Wiedemann (1828) is the most common in China and has a low potential vector rank, it has nevertheless long been thought to be an important natural malaria vector within the middle and lower reaches of the Yangtze River. A number of previous research studies have found evidence to support the occurrence of natural hybridization between An. sinensis and Anopheles kleini Rueda, 2005 (a competent malaria vector). We, therefore, collected a sample series of An. sinensis and morphologically similar species across China and undertook ribosomal and mitochondrial DNA analyses in order to assess genetic differentiation (Fst) and gene flow (Nm) amongst different groups. This enabled us to evaluate divergence times between morphologically similar species using the cytochrome oxidase I (COI) gene. The results of this study reveal significant genetic similarities between An. sinensis, An. kleini, and Anopheles belenrae Rueda, 2005 and therefore imply that correct molecular identifications will require additional molecular markers. As results also reveal the presence of gene flow between these three species, their taxonomic status will require further work. Data suggest that An. kleini is the most basal of the three species, while An. sinensis and An. belenrae share the closest genetic relationship.

Introgression of a synthetic sex ratio distortion system from Anopheles gambiae into Anopheles arabiensis

I-PpoI is a homing endonuclease that has a high cleavage activity and specificity for a conserved sequence within the ribosomal rDNA repeats, located in a single cluster on the Anopheles gambiae X chromosome. This property has been exploited to develop a synthetic sex ratio distortion system in this mosquito species. When I-PpoI is expressed from a transgene during spermatogenesis in mosquitoes, the paternal X chromosome is shredded and only Y chromosome-bearing sperm are viable, resulting in a male-biased sex ratio of >95% in the progeny. These distorter male mosquitoes can efficiently suppress caged wild-type populations, providing a powerful tool for vector control strategies. Given that malaria mosquito vectors belong to a species complex comprising at least two major vectors, we investigated whether the sex distorter I-PpoI, originally integrated in the A. gambiae genome, could be transferred via introgression to the sibling vector species Anopheles arabiensis. In compliance with Haldane’s rule, F1 hybrid male sterility is known to occur in all intercrosses among members of the Anopheles gambiae complex. A scheme based on genetic crosses and transgene selection was used to bypass F1 hybrid male sterility and introgress the sex distorter I-PpoI into the A. arabiensis genetic background. Our data suggest that this sex distortion technique can be successfully applied to target A. arabiensis mosquitoes.

Cross-Species Y Chromosome Function Between Malaria Vectors of the Anopheles gambiae Species Complex

Y chromosome function, structure and evolution is poorly understood in many species, including the Anopheles genus of mosquitoes-an emerging model system for studying speciation that also represents the major vectors of malaria. While the Anopheline Y had previously been implicated in male mating behavior, recent data from the Anopheles gambiae complex suggests that, apart from the putative primary sex-determiner, no other genes are conserved on the Y. Studying the functional basis of the evolutionary divergence of the Y chromosome in the gambiae complex is complicated by complete F1 male hybrid sterility. Here, we used an F1 × F0 crossing scheme to overcome a severe bottleneck of male hybrid incompatibilities that enabled us to experimentally purify a genetically labeled A. gambiae Y chromosome in an A. arabiensis background. Whole genome sequencing (WGS) confirmed that the A. gambiae Y retained its original sequence content in the A. arabiensis genomic background. In contrast to comparable experiments in Drosophila, we find that the presence of a heterospecific Y chromosome has no significant effect on the expression of A. arabiensis genes, and transcriptional differences can be explained almost exclusively as a direct consequence of transcripts arising from sequence elements present on the A. gambiae Y chromosome itself. We find that Y hybrids show no obvious fertility defects, and no substantial reduction in male competitiveness. Our results demonstrate that, despite their radically different structure, Y chromosomes of these two species of the gambiae complex that diverged an estimated 1.85 MYA function interchangeably, thus indicating that the Y chromosome does not harbor loci contributing to hybrid incompatibility. Therefore, Y chromosome gene flow between members of the gambiae complex is possible even at their current level of divergence. Importantly, this also suggests that malaria control interventions based on sex-distorting Y drive would be transferable, whether intentionally or contingent, between the major malaria vector species.

Identification of field caught Anopheles gambiae s.s. and Anopheles arabiensis by TaqMan single nucleotide polymorphism genotyping

Background

Identification of Anopheles gambiae s.s. and Anopheles arabiensis from field-collected Anopheles gambiae s.l. is often necessary in basic and applied research, and in operational control programmes. The currently accepted method involves use of standard polymerase chain reaction amplification of ribosomal DNA (rDNA) from the 3' 28S to 5' intergenic spacer region of the genome, and visual confirmation of amplicons of predicted size on agarose gels, after electrophoresis. This report describes development and evaluation of an automated, quantitative PCR method based upon TaqMan™ single nucleotide polymorphism (SNP) genotyping.

Methods

Standard PCR, and TaqMan SNP genotyping with newly designed primers and fluorophore-labeled probes hybridizing to sequences of complementary rDNA specific for either An. gambiae s.s. or An. arabiensis, were conducted in three experiments involving field-collected An. gambiae s.l. from western Kenya, and defined laboratory strains. DNA extraction was from a single leg, sonicated for five minutes in buffer in wells of 96-well PCR plates.

Results

TaqMan SNP genotyping showed a reaction success rate, sensitivity, and species specificity comparable to that of standard PCR. In an extensive field study, only 29 of 3,041 (0.95%) were determined to be hybrids by TaqMan (i.e., having rDNA sequences from both species), however, all but one were An. arabiensis by standard PCR, suggesting an acceptably low (ca. 1%) error rate for TaqMan genotyping in mistakenly identifying species hybrids.

Conclusion

TaqMan SNP genotyping proved to be a sensitive and rapid method for identification of An. gambiae s.l. and An. arabiensis, with a high success rate, specific results, and congruence with the standard PCR method.