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

The Anopheles leucine-rich repeat protein APL1C is a pathogen binding factor recognizing Plasmodium ookinetes and sporozoites

Leucine-rich repeat (LRR) proteins are commonly involved in innate immunity of animals and plants, including for pattern recognition of pathogen-derived elicitors. The Anopheles secreted LRR proteins APL1C and LRIM1 are required for malaria ookinete killing in conjunction with the complement-like TEP1 protein. However, the mechanism of parasite immune recognition by the mosquito remains unclear, although it is known that TEP1 lacks inherent binding specificity. Here, we find that APL1C and LRIM1 bind specifically to Plasmodium berghei ookinetes, even after depletion of TEP1 transcript and protein, consistent with a role for the LRR proteins in pathogen recognition. Moreover, APL1C does not bind to ookinetes of the human malaria parasite Plasmodium falciparum, and is not required for killing of this parasite, which correlates LRR binding specificity and immune protection. Most of the live P. berghei ookinetes that migrated into the extracellular space exposed to mosquito hemolymph, and almost all dead ookinetes, are bound by APL1C, thus associating LRR protein binding with parasite killing. We also find that APL1C binds to the surface of P. berghei sporozoites released from oocysts into the mosquito hemocoel and forms a potent barrier limiting salivary gland invasion and mosquito infectivity. Pathogen binding by APL1C provides the first functional explanation for the long-known requirement of APL1C for P. berghei ookinete killing in the mosquito midgut. We propose that secreted mosquito LRR proteins are required for pathogen discrimination and orientation of immune effector activity, potentially as functional counterparts of the immunoglobulin-based receptors used by vertebrates for antigen recognition.

Gene Conversion Explains Elevated Diversity in the Immunity Modulating APL1 Gene of the Malaria Vector Anopheles funestus

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to Plasmodium and other microbial pathogens in Anopheles mosquitoes. The APL1 gene of the malaria vector Anopheles funestus has exceptional levels of non-synonymous polymorphism across the range of An. funestus, with an average πn of 0.027 versus a genome-wide average of 0.002, and πn is consistently high in populations across Africa. Elevated APL1 diversity was consistent between the independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide resistance was observed. Although lacking the diversity of APL1, two further mosquito innate-immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. The cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, a result discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogues, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. In conclusion, we hypothesize that higher gene expression of APL1 than its paralogues is correlated with a more open chromatin formation, which enhances gene conversion and elevated diversity at this locus.

The ecdysone receptor regulates several key physiological factors in Anopheles funestus

Background

Malaria is a devastating disease, transmitted by female Anopheles mosquitoes infected with Plasmodium parasites. Current insecticide-based strategies exist to control the spread of malaria by targeting vectors. However, the increase in insecticide resistance in vector populations hinder the efficacy of these methods. It is, therefore, essential to develop novel vector control methods that efficiently target transmission reducing factors such as vector density and competence. A possible vector control candidate gene, the ecdysone receptor, regulates longevity, reproduction, immunity and other physiological processes in several insects, including malaria vectors. Anopheles funestus is a prominent vector in sub-Saharan Africa, however, the function of the ecdysone receptor in this mosquito has not previously been studied. This study aimed to determine if the ecdysone receptor depletion impacts An. funestus longevity, reproduction and susceptibility to Plasmodium falciparum infection.

Methods

RNA interference was used to reduce ecdysone receptor expression levels in An. funestus females and investigate how the above-mentioned phenotypes are influenced. Additionally, the expression levels of the ecdysone receptor, and reproduction genes lipophorin and vitellogenin receptor as well as the immune gene, leucine rich immune molecule 9 were determined in ecdysone receptor-depleted mosquitoes using quantitative polymerase chain reaction.

Results

Ecdysone receptor-depleted mosquitoes had a shorter lifespan, impaired oogenesis, were less fertile, and had reduced P. falciparum infection intensity.

Conclusions

Overall, this study provides the first experimental evidence that supports ecdysone receptor as a potential target in the development of vector control measures targeting An. funestus.

Functional constraints on insect immune system components govern their evolutionary trajectories

Roles of constraints in shaping evolutionary outcomes are often considered in the contexts of developmental biology and population genetics, in terms of capacities to generate new variants and how selection limits or promotes consequent phenotypic changes. Comparative genomics also recognizes the role of constraints, in terms of shaping evolution of gene and genome architectures, sequence evolutionary rates, and gene gains or losses, as well as on molecular phenotypes. Characterizing patterns of genomic change where putative functions and interactions of system components are relatively well described offers opportunities to explore whether genes with similar roles exhibit similar evolutionary trajectories. Using insect immunity as our test case system, we hypothesize that characterizing gene evolutionary histories can define distinct dynamics associated with different functional roles. We develop metrics that quantify gene evolutionary histories, employ these to characterize evolutionary features of immune gene repertoires, and explore relationships between gene family evolutionary profiles and their roles in immunity to understand how different constraints may relate to distinct dynamics. We identified three main axes of evolutionary trajectories characterized by gene duplication and synteny, maintenance/stability and sequence conservation, and loss and sequence divergence, highlighting similar and contrasting patterns across these axes amongst subsets of immune genes. Our results suggest that where and how genes participate in immune responses limit the range of possible evolutionary scenarios they exhibit. The test case study system of insect immunity highlights the potential of applying comparative genomics approaches to characterize how functional constraints on different components of biological systems govern their evolutionary trajectories.

Hemocyte RNA-Seq analysis of Indian malarial vectors Anopheles stephensi and Anopheles culicifacies: From similarities to differences

Anopheles stephensi and Anopheles culicifacies are dominant malarial vectors in urban and rural India, respectively. Both species carry significant biological differences in their behavioral adaptation and immunity, but the genetic basis of these variations are still poorly understood. Here, we uncovered the genetic differences of immune blood cells, that influence several immune-physiological responses. We generated, analyzed and compared the hemocyte RNA-Seq database of both mosquitoes. A total of 5,837,223,769 assembled bases collapsed into 7,595 and 3,791 transcripts, originating from hemocytes of laboratory-reared 3-4 days old naïve (sugar-fed) mosquitoes, Anopheles stephensi and Anopheles culicifacies respectively. Comparative GO annotation analysis revealed that both mosquito hemocytes encode similar proteins. Furthermore, while An. stephensi hemocytes showed a higher percentage of immune transcripts encoding APHAG (Autophagy), IMD (Immune deficiency pathway), PRDX (Peroxiredoxin), SCR (Scavenger receptor), IAP (Inhibitor of apoptosis), GALE (galactoside binding lectins), BGBPs (1,3 beta D glucan binding proteins), CASPs (caspases) and SRRP (Small RNA regulatory pathway), An. culicifacies hemocytes yielded a relatively higher percentage of transcripts encoding CLIP (Clip domain serine protease), FREP (Fibrinogen related proteins), PPO (Prophenol oxidase), SRPN (Serpines), ML (Myeloid differentiation 2-related lipid recognition protein), Toll path and TEP (Thioester protein), family proteins. However, a detailed comparative Interproscan analysis showed An. stephensi mosquito hemocytes encode proteins with increased repeat numbers as compared to An. culicifacies. Notably, we observed an abundance of transcripts showing significant variability of encoded proteins with repeats such as LRR (Leucine rich repeat), WD40 (W-D dipeptide), Ankyrin, Annexin, Tetratricopeptide and Mitochondrial substrate carrier repeat-containing family proteins, which may have a direct influence on species-specific immune-physiological responses. Summarily, our deep sequencing analysis unraveled that An. stephensi evolved with an expansion of repeat sequences in hemocyte proteins as compared to An. culicifacies, possibly providing an advantage for better adaptation to diverse environments.

A chromosome-level assembly of the black tiger shrimp (Penaeus monodon) genome facilitates the identification of growth-associated genes

To salvage marine ecosystems from fishery overexploitation, sustainable and efficient aquaculture must be emphasized. The knowledge obtained from available genome sequence of marine organisms has accelerated marine aquaculture in many cases. The black tiger shrimp (Penaeus monodon) is one of the most prominent cultured penaeid shrimps (Crustacean) with an average annual global production of half a million tons in the last decade. However, its currently available genome assemblies lack the contiguity and completeness required for accurate genome annotation due to the highly repetitive nature of the genome and technical difficulty in extracting high-quality, high-molecular weight DNA. Here, we report the first chromosome-level whole-genome assembly of P. monodon. The combination of long-read Pacific Biosciences (PacBio) and long-range Chicago and Hi-C technologies enabled a successful assembly of this first high-quality genome sequence. The final assembly covered 2.39 Gb (92.3% of the estimated genome size) and contained 44 pseudomolecules, corresponding to the haploid chromosome number. Repetitive elements occupied a substantial portion of the assembly (62.5%), the highest of the figures reported among crustacean species. The availability of this high-quality genome assembly enabled the identification of genes associated with rapid growth in the black tiger shrimp through the comparison of hepatopancreas transcriptome of slow-growing and fast-growing shrimps. The results highlighted several growth-associated genes. Our high-quality genome assembly provides an invaluable resource for genetic improvement and breeding penaeid shrimp in aquaculture. The availability of P. monodon genome enables analyses of ecological impact, environment adaptation and evolution, as well as the role of the genome to protect the ecological resources by promoting sustainable shrimp farming.

Delimiting cryptic morphological variation among human malaria vector species using convolutional neural networks

Deep learning is a powerful approach for distinguishing classes of images, and there is a growing interest in applying these methods to delimit species, particularly in the identification of mosquito vectors. Visual identification of mosquito species is the foundation of mosquito-borne disease surveillance and management, but can be hindered by cryptic morphological variation in mosquito vector species complexes such as the malaria-transmitting Anopheles gambiae complex. We sought to apply Convolutional Neural Networks (CNNs) to images of mosquitoes as a proof-of-concept to determine the feasibility of automatic classification of mosquito sex, genus, species, and strains using whole-body, 2D images of mosquitoes. We introduce a library of 1, 709 images of adult mosquitoes collected from 16 colonies of mosquito vector species and strains originating from five geographic regions, with 4 cryptic species not readily distinguishable morphologically even by trained medical entomologists. We present a methodology for image processing, data augmentation, and training and validation of a CNN. Our best CNN configuration achieved high prediction accuracies of 96.96% for species identification and 98.48% for sex. Our results demonstrate that CNNs can delimit species with cryptic morphological variation, 2 strains of a single species, and specimens from a single colony stored using two different methods. We present visualizations of the CNN feature space and predictions for interpretation of our results, and we further discuss applications of our findings for future applications in malaria mosquito surveillance.

Clusters of polymorphic transmembrane genes control resistance to schistosomes in snail vectors

Schistosomiasis is a debilitating parasitic disease infecting hundreds of millions of people. Schistosomes use aquatic snails as intermediate hosts. A promising avenue for disease control involves leveraging innate host mechanisms to reduce snail vectorial capacity. In a genome-wide association study of Biomphalaria glabrata snails, we identify genomic region PTC2 which exhibits the largest known correlation with susceptibility to parasite infection (>15 fold effect). Using new genome assemblies with substantially higher contiguity than the Biomphalaria reference genome, we show that PTC2 haplotypes are exceptionally divergent in structure and sequence. This variation includes multi-kilobase indels containing entire genes, and orthologs for which most amino acid residues are polymorphic. RNA-Seq annotation reveals that most of these genes encode single-pass transmembrane proteins, as seen in another resistance region in the same species. Such groups of hyperdiverse snail proteins may mediate host-parasite interaction at the cell surface, offering promising targets for blocking the transmission of schistosomiasis.