In silico karyotyping of chromosomally polymorphic malaria mosquitoes in the Anopheles gambiae complex

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained

Highly specific PCR-RFLP assays for karyotyping the widespread 2Rb inversion in malaria vectors of the Anopheles gambiae complex

Background: Chromosomal inversion polymorphisms play a role in adaptation to heterogeneous environments. Inversion polymorphisms are implicated in the very high ecological flexibility of the three main malaria vector species of the Afrotropical Anopheles gambiae complex, facilitating the exploitation of anthropogenic environmental modifications and promoting a strong association with humans. In addition to extending the species' spatial and temporal distribution, inversions are associated with epidemiologically relevant mosquito behavior and physiology, underscoring their medical importance. We here present novel PCR-RFLP based assays strongly predictive of genotype for the cosmopolitan 2Rb inversion in An. coluzzii and An. gambiae, a development which overcomes the numerous constraints inherent to traditional cytological karyotyping. Methods: We designed PCR-RFLP genotyping assays based on tag SNPs previously computationally identified as strongly predictive (> 95%) of 2Rb genotype. We targeted those tags whose alternative allelic states destroyed or created the recognition site of a commercially available restriction enzyme, and designed assays with distinctive cleavage profiles for each inversion genotype. The assays were validated on 251 An. coluzzii and 451 An. gambiae cytologically karyotyped specimens from nine countries across Africa and one An. coluzzii laboratory colony. Results: For three tag SNPs, PCR-RFLP assays (denoted DraIII, MspAI, and TatI) reliably produced robust amplicons and clearly distinguishable electrophoretic profiles for all three inversion genotypes. Results obtained with the DraIII assay are ≥ 95% concordant with cytogenetic assignments in both species, while MspAI and TatI assays produce patterns highly concordant with cytogenetic assignments only in An. coluzzii or An. gambiae, respectively. Joint application of species-appropriate pairs of assays increased the concordance levels to > 99% in An. coluzzii and 98% in An. gambiae. Potential sources of discordance (e.g. imperfect association between tag and inversion, allelic dropout, additional polymorphisms in the restriction target site, incomplete or failed restriction digestion) are discussed. Conclusions: The availability of highly specific, cost effective and accessible molecular assays for genotyping 2Rb in An. gambiae and An. coluzzii allows karyotyping of both sexes and all developmental stages. These novel tools will accelerate deeper investigations into the role of this ecologically and epidemiologically important chromosomal inversion in vector biology.[Figure not available: see fulltext.

Inversion genotyping in the Anopheles gambiae complex using high-throughput array and sequencing platforms

Chromosomal inversion polymorphisms have special importance in the Anopheles gambiae complex of malaria vector mosquitoes, due to their role in local adaptation and range expansion. The study of inversions in natural populations is reliant on polytene chromosome analysis by expert cytogeneticists, a process that is limited by the rarity of trained specialists, low throughput, and restrictive sampling requirements. To overcome this barrier, we ascertained tag single nucleotide polymorphisms (SNPs) that are highly correlated with inversion status (inverted or standard orientation). We compared the performance of the tag SNPs using two alternative high throughput molecular genotyping approaches vs. traditional cytogenetic karyotyping of the same 960 individual An. gambiae and An. coluzzii mosquitoes sampled from Burkina Faso, West Africa. We show that both molecular approaches yield comparable results, and that either one performs as well or better than cytogenetics in terms of genotyping accuracy. Given the ability of molecular genotyping approaches to be conducted at scale and at relatively low cost without restriction on mosquito sex or developmental stage, molecular genotyping via tag SNPs has the potential to revitalize research into the role of chromosomal inversions in the behavior and ongoing adaptation of An. gambiae and An. coluzzii to environmental heterogeneities

SNP Genotyping Defines Complex Gene-Flow Boundaries Among African Malaria Vector Mosquitoes

Mosquitoes in the Anopheles gambiae complex show rapid ecological and behavioral diversification, traits that promote malaria transmission and complicate vector control efforts. A high-density, genome-wide mosquito SNP-genotyping array allowed mapping of genomic differentiation between populations and species that exhibit varying levels of reproductive isolation. Regions near centromeres or within polymorphic inversions exhibited the greatest genetic divergence, but divergence was also observed elsewhere in the genomes. Signals of natural selection within populations were overrepresented among genomic regions that are differentiated between populations, implying that differentiation is often driven by population-specific selective events. Complex genomic differentiation among speciating vector mosquito populations implies that tools for genome-wide monitoring of population structure will prove useful for the advancement of malaria eradication

A PCR-RFLP method for genotyping of inversion 2Rc in Anopheles coluzzii

Background: Genotyping of polymorphic chromosomal inversions in malaria vectors such as An. coluzzii Coetzee & Wilkerson is important, both because they cause cryptic population structure that can mislead vector analysis and control and because they influence epidemiologically relevant eco-phenotypes. The conventional cytogenetic method of genotyping is an impediment because it is labor intensive, requires specialized training, and can be applied only to one gender and developmental stage. Here, we circumvent these limitations by developing a simple and rapid molecular method of genotyping inversion 2Rc in An. coluzzii that is both economical and field-friendly. This inversion is strongly implicated in temporal and spatial adaptations to climatic and ecological variation, particularly aridity. Methods: Using a set of tag single-nucleotide polymorphisms (SNPs) strongly correlated with inversion orientation, we identified those that overlapped restriction enzyme recognition sites and developed four polymerase chain reaction (PCR) restriction fragment length polymorphism (RFLP) assays that distinguish alternative allelic states at the tag SNPs. We assessed the performance of these assays using mosquito population samples from Burkina Faso that had been cytogenetically karyotyped as well as genotyped, using two complementary high-throughput molecular methods based on tag SNPs. Further validation was performed using mosquito population samples from additional West African (Benin, Mali, Senegal) and Central African (Cameroon) countries. Results: Of four assays tested, two were concordant with the 2Rc cytogenetic karyotype > 90% of the time in all samples. We recommend that these two assays be employed in tandem for reliable genotyping. By accepting only those genotypic assignments where both assays agree, > 99% of assignments are expected to be accurate. Conclusions: We have developed tandem PCR-RFLP assays for the accurate genotyping of inversion 2Rc in An. coluzzii. Because this approach is simple, inexpensive, and requires only basic molecular biology equipment, it is widely accessible. These provide a crucial tool for probing the molecular basis of eco-phenotypes relevant to malaria epidemiology and vector control. [Figure not available: see fulltext.

Comparative Analysis of the Global Transcriptome of Anopheles funestus from Mali, West Africa

Background: Anopheles funestus is a principal vector of malaria across much of tropical Africa and is considered one of the most efficient of its kind, yet studies of this species have lagged behind those of its broadly sympatric congener, An. gambiae. In aid of future genomic sequencing of An. funestus, we explored the whole body transcriptome, derived from mixed stage progeny of wild-caught females from Mali, West Africa. Principal Findings: Here we report the functional annotation and comparative genomics of 2,005 expressed sequence tags (ESTs) from An. funestus, which were assembled with a previous EST set from adult female salivary glands from the same mosquito. The assembled ESTs provided for a nonredundant catalog of 1,035 transcripts excluding mitochondrial sequences. Conclusions/Significance: Comparison of the An. funestus and An. gambiae transcriptomes using computational and macroarray approaches revealed a high degree of sequence identity despite an estimated 20–80 MY divergence time between lineages. A phylogenetically broader comparative genomic analysis indicated that the most rapidly evolving proteins – those involved in immunity, hematophagy, formation of extracellular structures, and hypothetical conserved proteins – are those that probably play important roles in how mosquitoes adapt to their nutritional and externa

Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii. We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures

Radical remodeling of the Y chromosome in a recent radiation of malaria mosquitoes

open28openHall A.B.; Papathanos P.-A.; Sharma A.; Cheng C.; Akbari O.S.; Assour L.; Bergman N.H.; Cagnetti A.; Crisanti A.; Dottorini T.; Fiorentini E.; Galizi R.; Hnath J.; Jiang X.; Koren S.; Nolan T.; Radune D.; Sharakhova M.V.; Steele A.; Timoshevskiy V.A.; Windbichler N.; Zhang S.; Hahn M.W.; Phillippy A.M.; Emrich S.J.; Sharakhov I.V.; Tu Z.J.; Besansky N.J.Hall, A. B.; Papathanos, P. -A.; SHARMA DHAKAL, Apsara; Cheng, C.; Akbari, O. S.; Assour, L.; Bergman, N. H.; Cagnetti, A.; Crisanti, A.; Dottorini, T.; Fiorentini, E.; Galizi, R.; Hnath, J.; Jiang, X.; Koren, S.; Nolan, T.; Radune, D.; Sharakhova, M. V.; Steele, A.; Timoshevskiy, V. A.; Windbichler, N.; Zhang, Shangu; Hahn, M. W.; Phillippy, A. M.; Emrich, S. J.; Sharakhov, I. V.; Tu, Z. J.; Besansky, N. J

Identification of Spiroplasma insolitum symbionts in Anopheles gambiae

Background: Insect symbionts have the potential to block the transmission of vector-borne diseases by their hosts. The advancement of a symbiont-based transmission blocking strategy for malaria requires the identification and study of Anopheles symbionts. Methods: High throughput 16S amplicon sequencing was used to profile the bacteria associated with Anopheles gambiae sensu lato and identify potential symbionts. The polymerase chain reaction (PCR) with specific primers were subsequently used to monitor symbiont prevalence in field populations, as well as symbiont transmission patterns. Results: We report the discovery of the bacterial symbiont, Spiroplasma, in Anopheles gambiae in Kenya. We determine that geographically dispersed Anopheles gambiae populations in Kenya are infected with Spiroplasma at low prevalence levels. Molecular phylogenetics indicates that this Anopheles gambiae associated Spiroplasma is a member of the insolitum clade. We demonstrate that this symbiont is stably maternally transmitted across at least two generations and does not significantly affect the fecundity or egg to adult survival of its host. Conclusions: In diverse insect species, Spiroplasma has been found to render their host resistant to infection by pathogens. The identification of a maternally transmitted strain of Spiroplasma in Anopheles gambiae may therefore open new lines of investigation for the development of symbiont-based strategies for blocking malaria transmission

VectorBase: a data resource for invertebrate vector genomics

VectorBase (http://www.vectorbase.org) is an NIAID-funded Bioinformatic Resource Center focused on invertebrate vectors of human pathogens. VectorBase annotates and curates vector genomes providing a web accessible integrated resource for the research community. Currently, VectorBase contains genome information for three mosquito species: Aedes aegypti, Anopheles gambiae and Culex quinquefasciatus, a body louse Pediculus humanus and a tick species Ixodes scapularis. Since our last report VectorBase has initiated a community annotation system, a microarray and gene expression repository and controlled vocabularies for anatomy and insecticide resistance. We have continued to develop both the software infrastructure and tools for interrogating the stored data