Kenneth D. Vernick

De Novo Transcriptome Sequencing in Anopheles funestus Using Illumina RNA-Seq Technology

Background Anopheles funestus is one of the primary vectors of human malaria, which causes a million deaths each year in sub-Saharan Africa. Few scientific resources are available to facilitate studies of this mosquito species and relatively little is known about its basic biology and evolution, making development and implementation of novel disease control efforts more difficult. The An. funestus genome has not been sequenced, so in order to facilitate genome-scale experimental biology, we have sequenced the adult female transcriptome of An. funestus from a newly founded colony in Burkina Faso, West Africa, using the Illumina GAIIx next generation sequencing platform. Methodology/Principal Findings We assembled short Illumina reads de novo using a novel approach involving iterative de novo assemblies and “target-based” contig clustering. We then selected a conservative set of 15,527 contigs through comparisons to four Dipteran transcriptomes as well as multiple functional and conserved protein domain databases. Comparison to the Anopheles gambiae immune system identified 339 contigs as putative immune genes, thus identifying a large portion of the immune system that can form the basis for subsequent studies of this important malaria vector. We identified 5,434 1∶1 orthologues between An. funestus and An. gambiae and found that among these 1∶1 orthologues, the protein sequence of those with putative immune function were significantly more diverged than the transcriptome as a whole. Short read alignments to the contig set revealed almost 367,000 genetic polymorphisms segregating in the An. funestus colony and demonstrated the utility of the assembled transcriptome for use in RNA-seq based measurements of gene expression. Conclusions/Significance We developed a pipeline that makes de novo transcriptome sequencing possible in virtually any organism at a very reasonable cost (

A Cryptic Subgroup of Anopheles gambiae Is Highly Susceptible to Human Malaria Parasites

6,300 in sequencing costs in our case). We anticipate that our approach could be used to develop genomic resources in a diversity of systems for which full genome sequence is currently unavailable. Our An. funestus contig set and analytical results provide a valuable resource for future studies in this non-model, but epidemiologically critical, vector insect.

Fine Pathogen Discrimination within the APL1 Gene Family Protects Anopheles gambiae against Human and Rodent Malaria Species

Collecting mosquito larvae from West African ponds has revealed a previously unknown but highly abundant genotype. 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.

Fluorescent Plasmodium berghei sporozoites and pre‐erythrocytic stages: a new tool to study mosquito and mammalian host interactions with malaria parasites

Genetically controlled resistance of Anopheles gambiae mosquitoes to Plasmodium falciparum is a common trait in the natural population, and a cluster of natural resistance loci were mapped to the Plasmodium-Resistance Island (PRI) of the A. gambiae genome. The APL1 family of leucine-rich repeat (LRR) proteins was highlighted by candidate gene studies in the PRI, and is comprised of paralogs APL1A, APL1B and APL1C that share ≥50% amino acid identity. Here, we present a functional analysis of the joint response of APL1 family members during mosquito infection with human and rodent Plasmodium species. Only paralog APL1A protected A. gambiae against infection with the human malaria parasite P. falciparum from both the field population and in vitro culture. In contrast, only paralog APL1C protected against the rodent malaria parasites P. berghei and P. yoelii. We show that anti-P. falciparum protection is mediated by the Imd/Rel2 pathway, while protection against P. berghei infection was shown to require Toll/Rel1 signaling. Further, only the short Rel2-S isoform and not the long Rel2-F isoform of Rel2 confers protection against P. falciparum. Protection correlates with the transcriptional regulation of APL1A by Rel2-S but not Rel2-F, suggesting that the Rel2-S anti-parasite phenotype results at least in part from its transcriptional control over APL1A. These results indicate that distinct members of the APL1 gene family display a mutually exclusive protective effect against different classes of Plasmodium parasites. It appears that a gene-for-pathogen-class system orients the appropriate host defenses against distinct categories of similar pathogens. It is known that insect innate immune pathways can distinguish between grossly different microbes such as Gram-positive bacteria, Gram-negative bacteria, or fungi, but the function of the APL1 paralogs reveals that mosquito innate immunity possesses a more fine-grained capacity to distinguish between classes of closely related eukaryotic pathogens than has been previously recognized.

A mosquito-specific protein family includes candidate receptors for malaria sporozoite invasion of salivary glands.

To track malaria parasites for biological studies within the mosquito and mammalian hosts, we constructed a stably transformed clonal line of Plasmodium berghei, PbFluspo, in which sporogonic and pre‐erythrocytic liver‐stage parasites are autonomously fluorescent. A cassette containing the structural gene for the FACS‐adapted green fluorescent protein mutant 2 (GFPmut2), expressed from the 5′ and 3′ flanking sequences of the circumsporozoite (CS) protein gene, was integrated and expressed at the endogenous CS locus. Recombinant parasites, which bear a wild‐type copy of CS, generated highly fluorescent oocysts and sporozoites that invaded mosquito salivary glands and were transmitted normally to rodent hosts. The parasites infected cultured hepatocytes in vitro, where they developed into fluorescent pre‐erythrocytic forms. Mammalian cells infected by these parasites can be separated from non‐infected cells by fluorescence activated cell sorter (FACS) analysis. These fluorescent insect and mammalian stages of P. berghei should be useful for phenotypic studies in their respective hosts, as well as for identification of new genes expressed in these parasite stages.

A switch in infected erythrocyte deformability at the maturation and blood circulation of Plasmodium falciparum transmission stages

We describe a previously unrecognized protein family from Aedes and Anopheles mosquitoes, here named SGS proteins. There are no SGS homologues in Drosophila or other eukaryotes, but SGS presence in two mosquito genera suggests that the protein family is widespread among mosquitoes. Ae. aegypti aaSGS1 mRNA and protein are salivary gland specific, and protein is localized in the basal lamina covering the anatomical regions that are preferentially invaded by malaria sporozoites. Anti‐aaSGS1 antibodies inhibited sporozoite invasion into the salivary glands in vivo, confirming aaSGS1 as a candidate sporozoite receptor. By homology to aaSGS1 we identified the complete complement of four SGS genes in An. gambiae, which were not recognized in the genome annotation. Two An. gambiae SGS genes display salivary gland specific expression like aaSGS1. Bioinformatic analysis predicts that SGS proteins possess heparin‐binding domains, and have among the highest density of tyrosine sulphation sites of all An. gambiae proteins. The major sporozoite surface proteins (CS and TRAP) also bind heparin, and interact with sulphoconjugates during liver cell invasion. Thus, we speculate that sporozoite invasion of mosquito salivary glands and subsequently the vertebrate liver may share similar mechanisms based on sulphation. Phylogenomic analysis suggests that an SGS ancestor was involved in a lateral gene transfer.

Anopheles gambiae APL1 Is a Family of Variable LRR Proteins Required for Rel1-Mediated Protection from the Malaria Parasite, Plasmodium berghei

Achievement of malaria elimination requires development of novel strategies interfering with parasite transmission, including targeting the parasite sexual stages (gametocytes). The formation of Plasmodium falciparum gametocytes in the human host takes several days during which immature gametocyte-infected erythrocytes (GIEs) sequester in host tissues. Only mature stage GIEs circulate in the peripheral blood, available to uptake by the Anopheles vector. Mechanisms underlying GIE sequestration and release in circulation are virtually unknown. We show here that mature GIEs are more deformable than immature stages using ektacytometry and microsphiltration methods, and that a switch in cellular deformability in the transition from immature to mature gametocytes is accompanied by the deassociation of parasite-derived STEVOR proteins from the infected erythrocyte membrane. We hypothesize that mechanical retention contributes to sequestration of immature GIEs and that regained deformability of mature gametocytes is associated with their release in the bloodstream and ability to circulate. These processes are proposed to play a key role in P falciparum gametocyte development in the host and to represent novel and unconventional targets for interfering with parasite transmission.

Antiviral immunity of Anopheles gambiae is highly compartmentalized, with distinct roles for RNA interference and gut microbiota

Background We previously identified by genetic mapping an Anopheles gambiae chromosome region with strong influence over the outcome of malaria parasite infection in nature. Candidate gene studies in the genetic interval, including functional tests using the rodent malaria parasite Plasmodium berghei, identified a novel leucine-rich repeat gene, APL1, with functional activity against P. berghei. Principal Findings Manual reannotation now reveals APL1 to be a family of at least 3 independently transcribed genes, APL1A, APL1B, and APL1C. Functional dissection indicates that among the three known APL1 family members, APL1C alone is responsible for host defense against P. berghei. APL1C functions within the Rel1-Cactus immune signaling pathway, which regulates APL1C transcript and protein abundance. Gene silencing of APL1C completely abolishes Rel1-mediated host protection against P. berghei, and thus the presence of APL1C is required for this protection. Further highlighting the influence of this chromosome region, allelic haplotypes at the APL1 locus are genetically associated with and have high explanatory power for the success or failure of P. berghei parasite infection. Conclusions APL1C functions as a required transducer of Rel1-dependent immune signal(s) to efficiently protect mosquitoes from P. berghei infection, and allelic genetic haplotypes of the APL1 locus display distinct levels of susceptibility and resistance to P. berghei.

Exceptional Diversity, Maintenance of Polymorphism, and Recent Directional Selection on the APL1 Malaria Resistance Genes of Anopheles gambiae

Significance It is important to understand antiviral mechanisms in potential new arbovirus vectors, such as Anopheles mosquitoes, in order to assess risks associated with arbovirus spread. Using an arbovirus naturally transmitted by Anopheles, we find that important immune mechanisms involved in the first bottleneck to Anopheles infection, the midgut, have distinct effects on arbovirus or malaria. This result is, to our knowledge, the first concrete evidence of protection tradeoffs for different human pathogens in a human disease vector, and it suggests that design of genetically immune-modified mosquitoes could result in unexpected outcomes. These results also indicate that different mosquito tissues display distinct antiviral protection that probably imposes divergent selection pressures upon viral replication during different stages of the infection. Arboviruses are transmitted by mosquitoes and other arthropods to humans and animals. The risk associated with these viruses is increasing worldwide, including new emergence in Europe and the Americas. Anopheline mosquitoes are vectors of human malaria but are believed to transmit one known arbovirus, o’nyong-nyong virus, whereas Aedes mosquitoes transmit many. Anopheles interactions with viruses have been little studied, and the initial antiviral response in the midgut has not been examined. Here, we determine the antiviral immune pathways of the Anopheles gambiae midgut, the initial site of viral infection after an infective blood meal. We compare them with the responses of the post-midgut systemic compartment, which is the site of the subsequent disseminated viral infection. Normal viral infection of the midgut requires bacterial flora and is inhibited by the activities of immune deficiency (Imd), JAK/STAT, and Leu-rich repeat immune factors. We show that the exogenous siRNA pathway, thought of as the canonical mosquito antiviral pathway, plays no detectable role in antiviral defense in the midgut but only protects later in the systemic compartment. These results alter the prevailing antiviral paradigm by describing distinct protective mechanisms in different body compartments and infection stages. Importantly, the presence of the midgut bacterial flora is required for full viral infectivity to Anopheles, in contrast to malaria infection, where the presence of the midgut bacterial flora is required for protection against infection. Thus, the enteric flora controls a reciprocal protection tradeoff in the vector for resistance to different human pathogens.

Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi

The three-gene APL1 locus encodes essential components of the mosquito immune defense against malaria parasites. APL1 was originally identified because it lies within a mapped QTL conferring the vector mosquito Anopheles gambiae natural resistance to the human malaria parasite, Plasmodium falciparum, and APL1 genes have subsequently been shown to be involved in defense against several species of Plasmodium. Here, we examine molecular population genetic variation at the APL1 gene cluster in spatially and temporally diverse West African collections of A. gambiae. The locus is extremely polymorphic, showing evidence of adaptive evolutionary maintenance of genetic variation. We hypothesize that this variability aids in defense against genetically diverse pathogens, including Plasmodium. Variation at APL1 is highly structured across geographic and temporal subpopulations. In particular, diversity is exceptionally high during the rainy season, when malaria transmission rates are at their peak. Much less allelic diversity is observed during the dry season when mosquito population sizes and malaria transmission rates are low. APL1 diversity is weakly stratified by the polymorphic 2La chromosomal inversion but is very strongly subdivided between the M and S “molecular forms.” We find evidence that a recent selective sweep has occurred at the APL1 locus in M form mosquitoes only. The independently reported observation of a similar M-form restricted sweep at the Tep1 locus, whose product physically interacts with APL1C, suggests that epistatic selection may act on these two loci causing them to sweep coordinately.