Correction: An Evolution-Based Screen for Genetic Differentiation between Anopheles Sister Taxa Enriches for Detection of Functional Immune Factors.

[This corrects the article DOI: 10.1371/journal.ppat.1005306.].

An Evolution-Based Screen for Genetic Differentiation between Anopheles Sister Taxa Enriches for Detection of Functional Immune Factors.

Nucleotide variation patterns across species are shaped by the processes of natural selection, including exposure to environmental pathogens. We examined patterns of genetic variation in two sister species, Anopheles gambiae and Anopheles coluzzii, both efficient natural vectors of human malaria in West Africa. We used the differentiation signature displayed by a known coordinate selective sweep of immune genes APL1 and TEP1 in A. coluzzii to design a population genetic screen trained on the sweep, classified a panel of 26 potential immune genes for concordance with the signature, and functionally tested their immune phenotypes. The screen results were strongly predictive for genes with protective immune phenotypes: genes meeting the screen criteria were significantly more likely to display a functional phenotype against malaria infection than genes not meeting the criteria (p = 0.0005). Thus, an evolution-based screen can efficiently prioritize candidate genes for labor-intensive downstream functional testing, and safely allow the elimination of genes not meeting the screen criteria. The suite of immune genes with characteristics similar to the APL1-TEP1 selective sweep appears to be more widespread in the A. coluzzii genome than previously recognized. The immune gene differentiation may be a consequence of adaptation of A. coluzzii to new pathogens encountered in its niche expansion during the separation from A. gambiae, although the role, if any of natural selection by Plasmodium is unknown. Application of the screen allowed identification of new functional immune factors, and assignment of new functions to known factors. We describe biochemical binding interactions between immune proteins that underlie functional activity for malaria infection, which highlights the interplay between pathogen specificity and the structure of immune complexes. We also find that most malaria-protective immune factors display phenotypes for either human or rodent malaria, with broad specificity a rarity.

Dynamics of the Human Infectious Reservoir for Malaria Determined by Mosquito Feeding Assays and Ultrasensitive Malaria Diagnosis in Burkina Faso.

BACKGROUND: Plasmodium falciparum gametocytes are essential for malaria transmission. Malaria control measures that aim at reducing transmission require an accurate characterization of the human infectious reservoir. METHODS: We longitudinally determined human infectiousness to mosquitoes and P. falciparum carriage by an ultrasensitive RNA-based diagnostics in 130 randomly selected inhabitants of an endemic area. RESULTS: At least 1 mosquito was infected by 32.6% (100 of 307) of the blood samples; in total, 7.6% of mosquitoes (916 of 12 079) were infected. The proportion of infectious individuals and infected mosquitoes were negatively associated with age and positively with asexual parasites (P < .001). Human infectiousness was higher at the start of the wet season and subsequently declined at the peak of the wet season (adjusted odds ratio, 0.52; P = .06) and in the dry season (0.23; P < .001). Overall, microscopy-negative individuals were responsible for 28.7% of infectious individuals (25 of 87) and 17.0% of mosquito infections (145 of 855). CONCLUSIONS: Our study reveals that the infectious reservoir peaks at the start of the wet season, with prominent roles for infections in children and submicroscopic infections. These findings have important consequences for strategies and the timing of interventions, which need to include submicroscopic infections and be implemented in the dry season.

Evolution of GOUNDRY, a cryptic subgroup of Anopheles gambiae s.l., and its impact on susceptibility to Plasmodium infection.

The recent discovery of a previously unknown genetic subgroup of Anopheles gambiae sensu lato underscores our incomplete understanding of complexities of vector population demographics in Anopheles. This subgroup, named GOUNDRY, does not rest indoors as adults and is highly susceptible to Plasmodium infection in the laboratory. Initial description of GOUNDRY suggested it differed from other known Anopheles taxa in surprising and sometimes contradictory ways, raising a number of questions about its age, population size and relationship to known subgroups. To address these questions, we sequenced the complete genomes of 12 wild-caught GOUNDRY specimens and compared these genomes to a panel of Anopheles genomes. We show that GOUNDRY is most closely related to Anopheles coluzzii, and the timing of cladogenesis is not recent, substantially predating the advent of agriculture. We find a large region of the X chromosome that has swept to fixation in GOUNDRY within the last 100 years, which may be an inversion that serves as a partial barrier to contemporary gene flow. Interestingly, we show that GOUNDRY has a history of inbreeding that is significantly associated with susceptibility to Plasmodium infection in the laboratory. Our results illuminate the genomic evolution of one of probably several cryptic, ecologically specialized subgroups of Anopheles and provide a potent example of how vector population dynamics may complicate efforts to control or eradicate malaria.

Association mapping by pooled sequencing identifies TOLL 11 as a protective factor against Plasmodium falciparum in Anopheles gambiae.

BACKGROUND: The genome-wide association study (GWAS) techniques that have been used for genetic mapping in other organisms have not been successfully applied to mosquitoes, which have genetic characteristics of high nucleotide diversity, low linkage disequilibrium, and complex population stratification that render population-based GWAS essentially unfeasible at realistic sample size and marker density. METHODS: We designed a novel mapping strategy for the mosquito system that combines the power of linkage mapping with the resolution afforded by genetic association. We established founder colonies from West Africa, controlled for diversity, linkage disequilibrium and population stratification. Colonies were challenged by feeding on the infectious stage of the human malaria parasite, Plasmodium falciparum, mosquitoes were phenotyped for parasite load, and DNA pools for phenotypically similar mosquitoes were Illumina sequenced. Phenotype-genotype mapping was carried out in two stages, coarse and fine. RESULTS: In the first mapping stage, pooled sequences were analysed genome-wide for intervals displaying relativereduction in diversity between phenotype pools, and candidate genomic loci were identified for influence upon parasite infection levels. In the second mapping stage, focused genotyping of SNPs from the first mapping stage was carried out in unpooled individual mosquitoes and replicates. The second stage confirmed significant SNPs in a locus encoding two Toll-family proteins. RNAi-mediated gene silencing and infection challenge revealed that TOLL 11 protects mosquitoes against P. falciparum infection. CONCLUSIONS: We present an efficient and cost-effective method for genetic mapping using natural variation segregating in defined recent Anopheles founder colonies, and demonstrate its applicability for mapping in a complex non-model genome. This approach is a practical and preferred alternative to population-based GWAS for first-pass mapping of phenotypes in Anopheles. This design should facilitate mapping of other traits involved in physiology, epidemiology, and behaviour.

The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing.

BACKGROUND: Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anopheles coluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known. METHODS: Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility. RESULTS: A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers. CONCLUSIONS: Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Plasmodium species occurrence, temporal distribution and interaction in a child-aged population in rural Burkina Faso.

BACKGROUND: Malaria can be caused by five Plasmodium species. Due to their higher prevalence, much of the research concentrates on Plasmodium falciparum and Plasmodium vivax. In Burkina Faso, where P. falciparum co-exists with Plasmodium malariae and Plasmodium ovale, there is not much data about the prevalence of the latter two species across human population. Moreover, interactions between co-infecting Plasmodium species are not documented. The aim of the current research is to determine species-specific prevalence and temporal distribution. The potential interactions between co-infecting Plasmodium species amongst the child-aged population in Burkina Faso are also discussed. METHODS: The study took place in the Sudanese savannah zone in Burkina Faso in a rural village, Laye. Surveys were conducted during the wet season across four years, 2007 to 2010. Volunteers aged three to 15 years with parental signed consent were enrolled. Ten children per week were screened for any history of pain, fever. Parasitological data were obtained by blood slide processing. RESULTS: Three sympatric Plasmodium species were recorded during this study with an average prevalence of 70.7%. Species temporal distribution showed an increase of P. malariae parasite prevalence from 0.9% in 2007 to 13.2% in 2010. Within a season, P. falciparum occurred in the overall study period while P. malariae and P. ovale were highly prevalent after the rainy part of this period. Species-specific infection analysis showed that in a comparison of mono-infections, P. malariae gametocyte prevalence and median density were higher than those of P. falciparum (88.9% vs 34.5% and 124.0 vs 40.0 gametocytes/µl, respectively). Likewise, in P. falciparum co-infections with P. malariae or P. ovale, gametocyte prevalence was also higher than in P. falciparum mono-infection. However, in P. falciparum mixed infection with P. malariae, P. falciparum gametocyte prevalence and median density as well as asexual form density decreased compared to P. falciparum mono-infection while for P. malariae mono-infection, only asexual form density significantly vary. CONCLUSION: These data revealed high gametocyte prevalence in other Plasmodium species than P. falciparum with a significant variation of P. malariae gametocyte carriers and gametocyte density across years. Molecular tools and entomological studies are needed to highly assess species-specific contribution to malaria transmission.

Equivalent susceptibility of Anopheles gambiae M and S molecular forms and Anopheles arabiensis to Plasmodium falciparum infection in Burkina Faso.

BACKGROUND: The Anopheles gambiae sensu lato (s.l.) species complex in Burkina Faso consists of Anopheles arabiensis, and molecular forms M and S of Anopheles gambiae sensu stricto (s.s.). Previous studies comparing the M and S forms for level of infection with Plasmodium falciparum have yielded conflicting results. METHODS: Mosquito larvae were sampled from natural pools, reared to adulthood under controlled conditions, and challenged with natural P. falciparum by experimental feeding with blood from gametocyte carriers. Oocyst infection prevalence and intensity was determined one week after infection. DNA from carcasses was genotyped to identify species and molecular form. RESULTS: In total, 7,400 adult mosquitoes grown from wild-caught larvae were challenged with gametocytes in 29 experimental infections spanning four transmission seasons. The overall infection prevalence averaged 40.7% for A. gambiae M form, 41.4% for A. gambiae S form, and 40.1% for A. arabiensis. There was no significant difference in infection prevalence or intensity between the three population groups. Notably, infection experiments in which the population groups were challenged in parallel on the same infective blood displayed less infection difference between population groups, while infections with less balanced composition of population groups had lower statistical power and displayed apparent differences that fluctuated more often from the null average. CONCLUSION: The study clearly establishes that, at the study site in Burkina Faso, there is no difference in genetic susceptibility to P. falciparum infection between three sympatric population groups of the A. gambiae s.l. complex. Feeding the mosquito groups on the same infective blood meal greatly increases statistical power. Conversely, comparison of the different mosquito groups between, rather than within, infections yields larger apparent difference between mosquito groups, resulting from lower statistical power and greater noise, and could lead to false-positive results. In making infection comparisons between population groups, it is more accurate to compare the different groups after feeding simultaneously upon the same infective blood.

Comprehensive Ecological and Geographic Characterization of Eukaryotic and Prokaryotic Microbiomes in African Anopheles.

Exposure of mosquitoes to numerous eukaryotic and prokaryotic microbes in their associated microbiomes has probably helped drive the evolution of the innate immune system. To our knowledge, a metagenomic catalog of the eukaryotic microbiome has not been reported from any insect. Here we employ a novel approach to preferentially deplete host 18S ribosomal RNA gene amplicons to reveal the composition of the eukaryotic microbial communities of Anopheles larvae sampled in Kenya, Burkina Faso and Republic of Guinea (Conakry). We identified 453 eukaryotic operational taxonomic units (OTUs) associated with Anopheles larvae in nature, but an average of 45% of the 18S rRNA sequences clustered into OTUs that lacked a taxonomic assignment in the Silva database. Thus, the Anopheles microbiome contains a striking proportion of novel eukaryotic taxa. Using sequence similarity matching and de novo phylogenetic placement, the fraction of unassigned sequences was reduced to an average of 4%, and many unclassified OTUs were assigned as relatives of known taxa. A novel taxon of the genus Ophryocystis in the phylum Apicomplexa (which also includes Plasmodium) is widespread in Anopheles larvae from East and West Africa. Notably, Ophryocystis is present at fluctuating abundance among larval breeding sites, consistent with the expected pattern of an epidemic pathogen. Species richness of the eukaryotic microbiome was not significantly different across sites from East to West Africa, while species richness of the prokaryotic microbiome was significantly lower in West Africa. Laboratory colonies of Anopheles coluzzii harbor 26 eukaryotic OTUs, of which 38% (n = 10) are shared with wild populations, while 16 OTUs are unique to the laboratory colonies. Genetically distinct An. coluzzii colonies co-housed in the same facility maintain different prokaryotic microbiome profiles, suggesting a persistent host genetic influence on microbiome composition. These results provide a foundation to understand the role of the Anopheles eukaryotic microbiome in vector immunity and pathogen transmission. We hypothesize that prevalent apicomplexans such as Ophryocystis associated with Anopheles could induce interference or competition against Plasmodium within the vector. This and other members of the eukaryotic microbiome may offer candidates for new vector control tools.

The Anopheles gambiae 2La chromosome inversion is associated with susceptibility to Plasmodium falciparum in Africa.

Chromosome inversions suppress genetic recombination and establish co-adapted gene complexes, or supergenes. The 2La inversion is a widespread polymorphism in the Anopheles gambiae species complex, the major African mosquito vectors of human malaria. Here we show that alleles of the 2La inversion are associated with natural malaria infection levels in wild-captured vectors from West and East Africa. Mosquitoes carrying the more-susceptible allele (2L+a) are also behaviorally less likely to be found inside houses. Vector control tools that target indoor-resting mosquitoes, such as bednets and insecticides, are currently the cornerstone of malaria control in Africa. Populations with high levels of the 2L+a allele may form reservoirs of persistent outdoor malaria transmission requiring novel measures for surveillance and control. The 2La inversion is a major and previously unappreciated component of the natural malaria transmission system in Africa, influencing both malaria susceptibility and vector behavior.

Genetic Structure of a Local Population of the Anopheles gambiae Complex in Burkina Faso.

Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. Population heterogeneities for ecological and behavioral attributes expand and stabilize malaria transmission over space and time, and populations may change in response to vector control, urbanization and other factors. There is a need for approaches to comprehensively describe the structure and characteristics of a sympatric local mosquito population, because incomplete knowledge of vector population composition may hinder control efforts. To this end, we used a genome-wide custom SNP typing array to analyze a population collection from a single geographic region in West Africa. The combination of sample depth (n = 456) and marker density (n = 1536) unambiguously resolved population subgroups, which were also compared for their relative susceptibility to natural genotypes of Plasmodium falciparum malaria. The population subgroups display fluctuating patterns of differentiation or sharing across the genome. Analysis of linkage disequilibrium identified 19 new candidate genes for association with underlying population divergence between sister taxa, A. coluzzii (M-form) and A. gambiae (S-form).

Reticulate Speciation and Barriers to Introgression in the Anopheles gambiae Species Complex.

Speciation as a process remains a central focus of evolutionary biology, but our understanding of the genomic architecture and prevalence of speciation in the face of gene flow remains incomplete. The Anopheles gambiae species complex of malaria mosquitoes is a radiation of ecologically diverse taxa. This complex is well-suited for testing for evidence of a speciation continuum and genomic barriers to introgression because its members exhibit partially overlapping geographic distributions as well as varying levels of divergence and reproductive isolation. We sequenced 20 genomes from wild A. gambiae s.s., Anopheles coluzzii, Anopheles arabiensis, and compared these with 12 genomes from the "GOUNDRY" subgroup of A. gambiae s.l. Amidst a backdrop of strong reproductive isolation, we find strong evidence for a speciation continuum with introgression of autosomal chromosomal regions among species and subgroups. The X chromosome, however, is strongly differentiated among all taxa, pointing to a disproportionately large effect of X chromosome genes in driving speciation among anophelines. Strikingly, we find that autosomal introgression has occurred from contemporary hybridization between A. gambiae and A. arabiensis despite strong divergence (∼5× higher than autosomal divergence) and isolation on the X chromosome. In addition to the X, we find strong evidence that lowly recombining autosomal regions, especially pericentromeric regions, serve as barriers to introgression secondarily to the X. We show that speciation with gene flow results in genomic mosaicism of divergence and introgression. Such a reticulate gene pool connecting vector taxa across the speciation continuum has important implications for malaria control efforts.

Population genetics of Anopheles coluzzii immune pathways and genes.

Natural selection is expected to drive adaptive evolution in genes involved in host-pathogen interactions. In this study, we use molecular population genetic analyses to understand how natural selection operates on the immune system of Anopheles coluzzii (formerly A. gambiae "M form"). We analyzed patterns of intraspecific and interspecific genetic variation in 20 immune-related genes and 17 nonimmune genes from a wild population of A. coluzzii and asked if patterns of genetic variation in the immune genes are consistent with pathogen-driven selection shaping the evolution of defense. We found evidence of a balanced polymorphism in CTLMA2, which encodes a C-type lectin involved in regulation of the melanization response. The two CTLMA2 haplotypes, which are distinguished by fixed amino acid differences near the predicted peptide cleavage site, are also segregating in the sister species A. gambiae ("S form") and A. arabiensis. Comparison of the two haplotypes between species indicates that they were not shared among the species through introgression, but rather that they arose before the species divergence and have been adaptively maintained as a balanced polymorphism in all three species. We additionally found that STAT-B, a retroduplicate of STAT-A, shows strong evidence of adaptive evolution that is consistent with neofunctionalization after duplication. In contrast to the striking patterns of adaptive evolution observed in these Anopheles-specific immune genes, we found no evidence of adaptive evolution in the Toll and Imd innate immune pathways that are orthologously conserved throughout insects. Genes encoding the Imd pathway exhibit high rates of amino acid divergence between Anopheles species but also display elevated amino acid diversity that is consistent with relaxed purifying selection. These results indicate that adaptive coevolution between A. coluzzii and its pathogens is more likely to involve novel or lineage-specific molecular mechanisms than the canonical humoral immune pathways.

Evidence for population-specific positive selection on immune genes of Anopheles gambiae.

Host-pathogen interactions can be powerful drivers of adaptive evolution, shaping the patterns of molecular variation at the genes involved. In this study, we sequenced alleles from 28 immune-related loci in wild samples of multiple genetic subpopulations of the African malaria mosquito Anopheles gambiae, obtaining unprecedented sample sizes and providing the first opportunity to contrast patterns of molecular evolution at immune-related loci in the recently discovered GOUNDRY population to those of the indoor-resting M and S molecular forms. In contrast to previous studies that focused on immune genes identified in laboratory studies, we centered our analysis on genes that fall within a quantitative trait locus associated with resistance to Plasmodium falciparum in natural populations of A. gambiae. Analyses of haplotypic and genetic diversity at these 28 loci revealed striking differences among populations in levels of genetic diversity and allele frequencies in coding sequence. Putative signals of positive selection were identified at 11 loci, but only one was shared by two subgroups of A. gambiae. Striking patterns of linkage disequilibrium were observed at several loci. We discuss these results with respect to ecological differences among these strata as well as potential implications for disease transmission.

A cryptic subgroup of Anopheles gambiae is highly susceptible to human malaria parasites.

Population subgroups of the African malaria vector Anopheles gambiae have not been comprehensively characterized owing to the lack of unbiased sampling methods. In the arid savanna zone of West Africa, where potential oviposition sites are scarce, widespread collection from larval pools in the peridomestic human habitat yielded a comprehensive genetic survey of local A. gambiae population subgroups, independent of adult resting behavior and ecological preference. A previously unknown subgroup of exophilic A. gambiae is sympatric with the known endophilic A. gambiae in this region. The exophilic subgroup is abundant, lacks differentiation into M and S molecular forms, and is highly susceptible to infection with wild Plasmodium falciparum. These findings might have implications for the epidemiology of malaria transmission and control.