De novo profiling of RNA viruses in Anopheles malaria vector mosquitoes from forest ecological zones in Senegal and Cambodia.

BACKGROUND: Mosquitoes are colonized by a large but mostly uncharacterized natural virome of RNA viruses, and the composition and distribution of the natural RNA virome may influence the biology and immunity of Anopheles malaria vector populations. RESULTS: Anopheles mosquitoes were sampled in malaria endemic forest village sites in Senegal and Cambodia, including Anopheles funestus, Anopheles gambiae group sp., and Anopheles coustani in Senegal, and Anopheles hyrcanus group sp., Anopheles maculatus group sp., and Anopheles dirus in Cambodia. The most frequent mosquito species sampled at both study sites are human malaria vectors. Small and long RNA sequences were depleted of mosquito host sequences, de novo assembled and clustered to yield non-redundant contigs longer than 500 nucleotides. Analysis of the assemblies by sequence similarity to known virus families yielded 115 novel virus sequences, and evidence supports a functional status for at least 86 of the novel viral contigs. Important monophyletic virus clades in the Bunyavirales and Mononegavirales orders were found in these Anopheles from Africa and Asia. The remaining non-host RNA assemblies that were unclassified by sequence similarity to known viruses were clustered by small RNA profiles, and 39 high-quality independent contigs strongly matched a pattern of classic RNAi processing of viral replication intermediates, suggesting they are entirely undescribed viruses. One thousand five hundred sixty-six additional high-quality unclassified contigs matched a pattern consistent with Piwi-interacting RNAs (piRNAs), suggesting that strand-biased piRNAs are generated from the natural virome in Anopheles. To functionally query piRNA effect, we analyzed piRNA expression in Anopheles coluzzii after infection with O'nyong nyong virus (family Togaviridae), and identified two piRNAs that appear to display specifically altered abundance upon arbovirus infection. CONCLUSIONS: Anopheles vectors of human malaria in Africa and Asia are ubiquitously colonized by RNA viruses, some of which are monophyletic but clearly diverged from other arthropod viruses. The interplay between small RNA pathways, immunity, and the virome may represent part of the homeostatic mechanism maintaining virome members in a commensal or nonpathogenic state, and could potentially influence vector competence.

Correction to: De novo profiling of RNA viruses in Anopheles malaria vector mosquitoes from forest ecological zones in Senegal and Cambodia.

Following the publication of this article [1], the authors reported that the original shading in columns 3 and 4 of Table 3, which indicated the presence or absence of viruses in each library, had been removed during typesetting.

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.

Composition and genetics of malaria vector populations in the Central African Republic.

BACKGROUND: In many African countries malaria has declined sharply due to a synergy of actions marked by the introduction of vector control strategies, but the disease remains a leading cause of morbidity and mortality in Central African Republic (CAR). An entomological study was initiated with the aim to characterize the malaria vectors in Bangui, the capital of CAR, and determine their vector competence. METHODS: A cross-sectional entomological study was conducted in 15 sites of the district of Bangui, the capital of CAR, in September-October 2013 and a second collection was done in four of those sites between November and December 2013. Mosquitoes were collected by human landing catch (HLC) indoors and outdoors and by pyrethrum spray catch of indoor-resting mosquitoes. Mosquitoes were analysed for species and multiple other attributes, including the presence of Plasmodium falciparum circumsporozoite protein or DNA, blood meal source, 2La inversion karyotype, and the L1014F kdr insecticide resistance mutation. RESULTS: Overall, 1292 anophelines were analysed, revealing a predominance of Anopheles gambiae and Anopheles funestus, with a small fraction of Anopheles coluzzii. Molecular typing of the An. gambiae complex species showed that An. gambiae was predominant (95.7 %) as compared to An. coluzzii (2.1 %), and Anopheles arabiensis was not present. In some areas the involvement of secondary vectors, such as Anopheles coustani, expands the risk of infection. By HLC sampling, An. funestus displayed a stronger endophilic preference than mosquitoes from the An. gambiae sister taxa, with a mean indoor-capture rate of 54.3 % and 67.58 % for An. gambiae sister taxa and An. funestus, respectively. Human biting rates were measured overall for each of the species with 28 or 29 bites/person/night, respectively. Both vectors displayed a strong human feeding preference as determined by blood meal source, which was not different between the different sampling sites. An. coustani appears to be highly exophilic, with 92 % of HLC samples captured outdoors. The mean CSP rate in head-thorax sections of all Anopheles was 5.09 %, and was higher in An. gambiae s.l. (7.4 %) than in An. funestus (3.3 %). CSP-positive An. coustani were also detected in outdoor HLC samples. In the mosquitoes of the An. gambiae sister taxa the kdr-w mutant allele was nearly fixed, with 92.3 % resistant homozygotes, and no susceptible homozygotes detected. CONCLUSIONS: This study collected data on anopheline populations in CAR, behaviour of vectors and transmission levels. Further studies should investigate the biting behaviour and susceptibility status of the anophelines to different insecticides to allow the establishment of appropriate vector control based on practical entomological knowledge.

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.

Molecular characterization and genetic authentication assay for Anopheles 'hemocyte-like' cell lines 4a-3A and 4a-3B.

BACKGROUND: Anopheles cell lines are used in a variety of ways to better understand the major vectors of malaria in sub-Saharan Africa. Despite this, commonly used cell lines are not well characterized, and no tools are available for cell line identification and authentication. METHODS: Utilizing whole genome sequencing, genomes of 4a-3A and 4a-3B 'hemocyte-like' cell lines were characterized for insertions and deletions (indels) and SNP variation. Genomic locations of distinguishing sequence variation and species origin of the cell lines were also examined. Unique indels were targeted to develop a PCR-based cell line authentication assay. Mitotic chromosomes were examined to survey the cytogenetic landscape for chromosome structure and copy number in the cell lines. RESULTS: The 4a-3A and 4a-3B cell lines are female in origin and primarily of Anopheles coluzzii ancestry. Cytogenetic analysis indicates that the two cell lines are essentially diploid, with some relatively minor chromosome structural rearrangements. Whole-genome sequence was generated, and analysis indicated that SNPs and indels which differentiate the cell lines are clustered on the 2R chromosome in the regions of the 2Rb, 2Rc and 2Ru chromosomal inversions. A PCR-based authentication assay was developed to fingerprint three indels unique to each cell line. The assay distinguishes between 4a-3A and 4a-3B cells and also uniquely identifies two additional An. coluzzii cell lines tested, Ag55 and Sua4.0. The assay has the specificity to distinguish four cell lines and also has the sensitivity to detect cellular contamination within a sample of cultured cells. CONCLUSIONS: Genomic characterization of the 4a-3A and 4a-3B Anopheles cell lines was used to develop a simple diagnostic assay that can distinguish these cell lines within and across research laboratories. A cytogenetic survey indicated that the 4a-3A and Sua4.0 cell lines carry essentially normal diploid chromosomes, which makes them amenable to CRISPR/Cas9 genome editing. The presented simple authentication assay, coupled with screening for mycoplasma, will allow validation of the integrity of experimental resources and will promote greater experimental reproducibility of results.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii.

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences. Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Leucine-Rich Immune Factor APL1 Is Associated With Specific Modulation of Enteric Microbiome Taxa in the Asian Malaria Mosquito Anopheles stephensi.

The commensal gut microbiome is contained by the enteric epithelial barrier, but little is known about the degree of specificity of host immune barrier interactions for particular bacterial taxa. Here, we show that depletion of leucine-rich repeat immune factor APL1 in the Asian malaria mosquito Anopheles stephensi is associated with higher midgut abundance of just the family Enterobacteraceae, and not generalized dysbiosis of the microbiome. The effect is explained by the response of a narrow clade containing two main taxa related to Klebsiella and Cedecea. Analysis of field samples indicate that these two taxa are recurrent members of the wild Anopheles microbiome. Triangulation using sequence and functional data incriminated relatives of C. neteri and Cedecea NFIX57 as candidates for the Cedecea component, and K. michiganensis, K. oxytoca, and K.sp. LTGPAF-6F as candidates for the Klebsiella component. APL1 presence is associated with host ability to specifically constrain the abundance of a narrow microbiome clade of the Enterobacteraceae, and the immune factor may promote homeostasis of this clade in the enteric microbiome for host benefit.

Gene copy number and function of the APL1 immune factor changed during Anopheles evolution.

BACKGROUND: The recent reference genome assembly and annotation of the Asian malaria vector Anopheles stephensi detected only one gene encoding the leucine-rich repeat immune factor APL1, while in the Anopheles gambiae and sibling Anopheles coluzzii, APL1 factors are encoded by a family of three paralogs. The phylogeny and biological function of the unique APL1 gene in An. stephensi have not yet been specifically examined. METHODS: The APL1 locus was manually annotated to confirm the computationally predicted single APL1 gene in An. stephensi. APL1 evolution within Anopheles was explored by phylogenomic analysis. The single or paralogous APL1 genes were silenced in An. stephensi and An. coluzzii, respectively, followed by mosquito survival analysis, experimental infection with Plasmodium and expression analysis. RESULTS: APL1 is present as a single ancestral gene in most Anopheles including An. stephensi but has expanded to three paralogs in an African lineage that includes only the Anopheles gambiae species complex and Anopheles christyi. Silencing of the unique APL1 copy in An. stephensi results in significant mosquito mortality. Elevated mortality of APL1-depleted An. stephensi is rescued by antibiotic treatment, suggesting that pathology due to bacteria is the cause of mortality, and indicating that the unique APL1 gene is essential for host survival. Successful Plasmodium development in An. stephensi depends upon APL1 activity for protection from high host mortality due to bacteria. In contrast, silencing of all three APL1 paralogs in An. coluzzii does not result in elevated mortality, either with or without Plasmodium infection. Expression of the single An. stephensi APL1 gene is regulated by both the Imd and Toll immune pathways, while the two signaling pathways regulate different APL1 paralogs in the expanded APL1 locus. CONCLUSIONS: APL1 underwent loss and gain of functions concomitant with expansion from a single ancestral gene to three paralogs in one lineage of African Anopheles. We infer that activity of the unique APL1 gene promotes longevity in An. stephensi by conferring protection from or tolerance to an effect of bacterial pathology. The evolution of an expanded APL1 gene family could be a factor contributing to the exceptional levels of malaria transmission mediated by human-feeding members of the An. gambiae species complex in Africa.

Influence of genetic polymorphism on transcriptional enhancer activity in the malaria vector Anopheles coluzzii.

Enhancers are cis-regulatory elements that control most of the developmental and spatial gene expression in eukaryotes. Genetic variation of enhancer sequences is known to influence phenotypes, but the effect of enhancer variation upon enhancer functional activity and downstream phenotypes has barely been examined in any species. In the African malaria vector, Anopheles coluzzii, we identified candidate enhancers in the proximity of genes relevant for immunity, insecticide resistance, and development. The candidate enhancers were functionally validated using luciferase reporter assays, and their activity was found to be essentially independent of their physical orientation, a typical property of enhancers. All of the enhancers segregated genetically polymorphic alleles, which displayed significantly different levels of functional activity. Deletion mutagenesis and functional testing revealed a fine structure of positive and negative regulatory elements that modulate activity of the enhancer core. Enhancer polymorphisms segregate in wild A. coluzzii populations in West Africa. Thus, enhancer variants that modify target gene expression leading to likely phenotypic consequences are frequent in nature. These results demonstrate the existence of naturally polymorphic A. coluzzii enhancers, which may help explain important differences between individuals or populations for malaria transmission efficiency and vector adaptation to the environment.

Comparison of naturally acquired antibody responses against the C-terminal processing products of Plasmodium vivax Merozoite Surface Protein-1 under low transmission and unstable malaria conditions in Sri Lanka.

We report here, for the first time, a comparison of naturally acquired antibody responses to the 42 and 19 kDa C-terminal processing products of Plasmodium vivax Merozoite Surface Protein-1 assayed by ELISA using p42 and p19 baculovirus-derived recombinant proteins, respectively. Test populations comprised patients with microscopy confirmed acute P. vivax infections from two regions endemic for vivax malaria where low transmission and unstable malaria conditions prevail, and a non-endemic urban area, in Sri Lanka. The antibody prevalence to the two proteins, both at the individual and population levels, tend to respond more to p42 than to p19 in all test areas, where >14% of individuals preferentially recognized p42, compared with <2% for p19. In patients with no previous exposure to malaria, 21% preferentially recognized p42, whereas none exclusively recognized p19. A significantly lower prevalence of anti-p19 IgM, but not anti-p42 IgM, was observed among residents from endemic areas compared with their non-endemic counterparts. Individuals from both endemic areas produced significantly less anti-p19 IgM compared with anti-p42 IgM. IgG1 was the predominant IgG isotype for both antigens in all individuals. With increasing exposure to malaria in both endemic areas, anti-p19 antibody responses were dominated by the functionally important IgG1 and IgG3 isotypes, with a concurrent reduction in IgM that was lacking in the non-endemic residents. This antibody switch was also reflected for PvAMA-1 as we previously reported with the identical battery of sera. In contrast, the antibody switch for p42 was restricted to endemic residents with more extensive exposure. These results suggest that an IgM-dominated antibody response against the p42 polymorphic region in endemic residents may interfere with the development of an IgG-dominated "protective" isotype shift to p19, that may complicate vaccine development.

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).

Identification and Characterization of Two Novel RNA Viruses from Anopheles gambiae Species Complex Mosquitoes.

Mosquitoes of the Anopheles gambiae complex display strong preference for human bloodmeals and are major malaria vectors in Africa. However, their interaction with viruses or role in arbovirus transmission during epidemics has been little examined, with the exception of O'nyong-nyong virus, closely related to Chikungunya virus. Deep-sequencing has revealed different RNA viruses in natural insect viromes, but none have been previously described in the Anopheles gambiae species complex. Here, we describe two novel insect RNA viruses, a Dicistrovirus and a Cypovirus, found in laboratory colonies of An. gambiae taxa using small-RNA deep sequencing. Sequence analysis was done with Metavisitor, an open-source bioinformatic pipeline for virus discovery and de novo genome assembly. Wild-collected Anopheles from Senegal and Cambodia were positive for the Dicistrovirus and Cypovirus, displaying high sequence identity to the laboratory-derived virus. Thus, the Dicistrovirus (Anopheles C virus, AnCV) and Cypovirus (Anopheles Cypovirus, AnCPV) are components of the natural virome of at least some anopheline species. Their possible influence on mosquito immunity or transmission of other pathogens is unknown. These natural viruses could be developed as models for the study of Anopheles-RNA virus interactions in low security laboratory settings, in an analogous manner to the use of rodent malaria parasites for studies of mosquito anti-parasite immunity.

Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation.

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.

Identification of the Anopheles gambiae ATP-binding cassette transporter superfamily genes.

The Anopheles gambiae genome sequence has been analyzed to find ATP-binding cassette protein genes based on deduced protein similarity to known family members. A nonredundant collection of 44 putative genes was identified including five genes not detected by the original Anopheles genome project machine annotation. These genes encode at least one member of all the human and Drosophila melanogaster ATP-binding protein subgroups. Like D. melanogaster, A. gambiae has subgroup ABCH genes encoding proteins different from the ABC proteins found in other complex organisms. The largest Anopheles subgroup is the ABCC genes which includes one member that can potentially encode ten different isoforms of the protein by differential splicing. As with Drosophila, the second largest Anopheles group is the ABCG subgroup with 12 genes compared to 15 genes in D. melanogaster, but only 5 genes in the human genome. In contrast, fewer ABCA and ABCB genes were identified in the mosquito genome than in the human or Drosophila genomes. Gene duplication is very evident in the Anopheles ABC genes with two groups of four genes, one group with three genes and three groups with two head to tail duplicated genes. These characteristics argue that the A. gambiae is actively using gene duplication as a mechanism to drive genetic variation in this important gene group.

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.

Diverged alleles of the Anopheles gambiae leucine-rich repeat gene APL1A display distinct protective profiles against Plasmodium falciparum.

Functional studies have demonstrated a role for the Anopheles gambiae APL1A gene in resistance against the human malaria parasite, Plasmodium falciparum. Here, we exhaustively characterize the structure of the APL1 locus and show that three structurally different APL1A alleles segregate in the Ngousso colony. Genetic association combined with RNAi-mediated gene silencing revealed that APL1A alleles display distinct protective profiles against P. falciparum. One APL1A allele is sufficient to explain the protective phenotype of APL1A observed in silencing experiments. Epitope-tagged APL1A isoforms expressed in an in vitro hemocyte-like cell system showed that under assay conditions, the most protective APL1A isoform (APL1A(2)) localizes within large cytoplasmic vesicles, is not constitutively secreted, and forms only one protein complex, while a less protective isoform (APL1A(1)) is constitutively secreted in at least two protein complexes. The tested alleles are identical to natural variants in the wild A. gambiae population, suggesting that APL1A genetic variation could be a factor underlying natural heterogeneity of vector susceptibility to P. falciparum.

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.