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.

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.

Photorhabdus luminescens, an entomopathogenic bacterium and nematode symbiont, has homologues of the Hca and Mhp enzymes. In Escherichia coli, these enzymes catalyze the degradation of the aromatic compounds 3-phenylpropionate (3PP) and cinnamic acid (CA) and allow the use of 3PP as sole carbon source. P. luminescens is not able to use 3PP and CA as sole carbon sources but can degrade them. Hca dioxygenase is involved in this degradation pathway. P. luminescens synthesizes CA from phenylalanine via a phenylalanine ammonia-lyase (PAL) and degrades it via the not-yet-characterized biosynthetic pathway of 3,5-dihydroxy-4-isopropylstilbene (ST) antibiotic. CA induces its own synthesis by enhancing the expression of the stlA gene that codes for PAL. P. luminescens bacteria release endogenous CA into the medium at the end of exponential growth and then consume it. Hca dioxygenase is involved in the consumption of endogenous CA but is not required for ST production. This suggests that CA is consumed via at least two separate pathways in P. luminescens: the biosynthesis of ST and a pathway involving the Hca and Mhp enzymes.

Bacteria inhabit the animal digestive tract and body and are generally required for health of the organism. However, some of the bacteria could be harmful if they proliferate to a high level in the host. The mechanisms that allow the host to harbor, control and benefit from bacteria are not well understood. Here we show that a small group of bacteria that are widespread in Anopheles mosquitoes appear to be maintained at an appropriate level by the presence of an immune factor, APL1, and that loss of APL1 allows growth of only these few bacterial species.

Comparison of the proteomes of wild-type Photorhabdus luminescens and its hcaR derivative, grown in insect hemolymph, showed that hcaR disruption decreased the production of toxins (tcdA1, mcf, and pirAB) and proteins involved in oxidative stress response (SodA, AhpC, Gor). The disruption of hcaR did not affect growth rate in insects, but did delay the virulence of P. luminescens in Bombyx mori and Spodoptera littoralis larvae. This delayed virulence was associated with a lower toxemia rather than delay in bacteremia. The disruption of hcaR also increased bacterial sensitivity to hydrogen peroxide. A sodA mutant and an hcaR mutant had similar phenotypes in terms of sensitivity to hydrogen peroxide, virulence, toxin gene expression, and growth rate in insects. Thus, the two processes affected by hcaR disruption - toxemia and oxidative stress response - appear to be related. Besides, expression of toxin genes tcdA1, mcf, and pirAB was decreased by paraquat challenge. We provide here the first demonstration of the importance of toxemia for P. luminescens virulence. Our results also highlight the power of proteomic analysis for detecting unexpected links between different, concomitant processes in bacteria.

BackgroundThe 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.MethodsWe 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.ResultsIn 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.ConclusionsWe 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2009-z) contains supplementary material, which is available to authorized users.

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 (APL1A2) localizes within large cytoplasmic vesicles, is not constitutively secreted, and forms only one protein complex, while a less protective isoform (APL1A1) 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.

During a blood meal,