A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons.

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

Male vitellogenin regulates gametogenesis through a testis-enriched big protein in Chrysopa pallens.

In insects, vitellogenin (Vg) is generally viewed as a female-specific protein. Its primary function is to supply nutrition to developing embryos. Here, we reported Vg from the male adults of a natural predator, Chrysopa pallens. The male Vg was depleted by RNAi. Mating with Vg-deficient male downregulated female Vg expression, suppressed ovarian development and decreased reproductive output. Whole-organism transcriptome analysis after male Vg knockdown showed no differential expression of the known spermatogenesis-related regulators and seminal fluid protein genes, but a sharp downregulation of an unknown gene, which encodes a testis-enriched big protein (Vcsoo). Separate knockdown of male Vg and Vcsoo disturbed the assembly of spermatid cytoplasmic organelles in males and suppressed the expansion of ovary germarium in mated females. These results demonstrated that C. pallens male Vg signals through the downstream Vcsoo and regulates male and female reproduction.

CRISPR/Cas9 modified An. gambiae carrying kdr mutation L1014F functionally validate its contribution in insecticide resistance and combined effect with metabolic enzymes.

Insecticide resistance in Anopheles mosquitoes is a major obstacle in maintaining the momentum in reducing the malaria burden; mitigating strategies require improved understanding of the underlying mechanisms. Mutations in the target site of insecticides (the voltage gated sodium channel for the most widely used pyrethroid class) and over-expression of detoxification enzymes are commonly reported, but their relative contribution to phenotypic resistance remain poorly understood. Here we present a genome editing pipeline to introduce single nucleotide polymorphisms in An. gambiae which we have used to study the effect of the classical kdr mutation L1014F (L995F based on An. gambiae numbering), one of the most widely distributed resistance alleles. Introduction of 1014F in an otherwise fully susceptible genetic background increased levels of resistance to all tested pyrethroids and DDT ranging from 9.9-fold for permethrin to >24-fold for DDT. The introduction of the 1014F allele was sufficient to reduce mortality of mosquitoes after exposure to deltamethrin treated bednets, even as the only resistance mechanism present. When 1014F was combined with over-expression of glutathione transferase Gste2, resistance to permethrin increased further demonstrating the critical combined effect between target site resistance and detoxification enzymes in vivo. We also show that mosquitoes carrying the 1014F allele in homozygosity showed fitness disadvantages including increased mortality at the larval stage and a reduction in fecundity and adult longevity, which can have consequences for the strength of selection that will apply to this allele in the field.

Regulating the expression of gene drives is key to increasing their invasive potential and the mitigation of resistance.

Homing-based gene drives use a germline source of nuclease to copy themselves at specific target sites in a genome and bias their inheritance. Such gene drives can be designed to spread and deliberately suppress populations of malaria mosquitoes by impairing female fertility. However, strong unintended fitness costs of the drive and a propensity to generate resistant mutations can limit a gene drive's potential to spread. Alternative germline regulatory sequences in the drive element confer improved fecundity of carrier individuals and reduced propensity for target site resistance. This is explained by reduced rates of end-joining repair of DNA breaks from parentally deposited nuclease in the embryo, which can produce heritable mutations that reduce gene drive penetrance. We tracked the generation and selection of resistant mutations over the course of a gene drive invasion of a population. Improved gene drives show faster invasion dynamics, increased suppressive effect and later onset of target site resistance. Our results show that regulation of nuclease expression is as important as the choice of target site when developing a robust homing-based gene drive for population suppression.

Resistance to a CRISPR-based gene drive at an evolutionarily conserved site is revealed by mimicking genotype fixation.

CRISPR-based homing gene drives can be designed to disrupt essential genes whilst biasing their own inheritance, leading to suppression of mosquito populations in the laboratory. This class of gene drives relies on CRISPR-Cas9 cleavage of a target sequence and copying ('homing') therein of the gene drive element from the homologous chromosome. However, target site mutations that are resistant to cleavage yet maintain the function of the essential gene are expected to be strongly selected for. Targeting functionally constrained regions where mutations are not easily tolerated should lower the probability of resistance. Evolutionary conservation at the sequence level is often a reliable indicator of functional constraint, though the actual level of underlying constraint between one conserved sequence and another can vary widely. Here we generated a novel adult lethal gene drive (ALGD) in the malaria vector Anopheles gambiae, targeting an ultra-conserved target site in a haplosufficient essential gene (AGAP029113) required during mosquito development, which fulfils many of the criteria for the target of a population suppression gene drive. We then designed a selection regime to experimentally assess the likelihood of generation and subsequent selection of gene drive resistant mutations at its target site. We simulated, in a caged population, a scenario where the gene drive was approaching fixation, where selection for resistance is expected to be strongest. Continuous sampling of the target locus revealed that a single, restorative, in-frame nucleotide substitution was selected. Our findings show that ultra-conservation alone need not be predictive of a site that is refractory to target site resistance. Our strategy to evaluate resistance in vivo could help to validate candidate gene drive targets for their resilience to resistance and help to improve predictions of the invasion dynamics of gene drives in field populations.

Transcriptional variation of sensory-related genes in natural populations of Aedes albopictus.

BACKGROUND: The Asian tiger mosquito, Aedes albopictus, is a highly dangerous invasive vector of numerous medically important arboviruses including dengue, chikungunya and Zika. In four decades it has spread from tropical Southeast Asia to many parts of the world in both tropical and temperate climes. The rapid invasion process of this mosquito is supported by its high ecological and genetic plasticity across different life history traits. Our aim was to investigate whether wild populations, both native and adventive, also display transcriptional genetic variability for functions that may impact their biology, behaviour and ability to transmit arboviruses, such as sensory perception. RESULTS: Antennal transcriptome data were derived from mosquitoes from a native population from Ban Rai, Thailand and from three adventive Mediterranean populations: Athens, Greece and Arco and Trento from Italy. Clear inter-population differential transcriptional activity was observed in different gene categories related to sound perception, olfaction and viral infection. The greatest differences were detected between the native Thai and the Mediterranean populations. The two Italian populations were the most similar. Nearly one million quality filtered SNP loci were identified. CONCLUSION: The ability to express this great inter-population transcriptional variability highlights, at the functional level, the remarkable genetic flexibility of this mosquito species. We can hypothesize that the differential expression of genes, including those involved in sensory perception, in different populations may enable Ae. albopictus to exploit different environments and hosts, thus contributing to its status as a global vector of arboviruses of public health importance. The large number of SNP loci present in these transcripts represents a useful addition to the arsenal of high-resolution molecular markers and a resource that can be used to detect selective pressure and adaptive changes that may have occurred during the colonization process.

Rapid evolution of female-biased genes among four species of Anopheles malaria mosquitoes.

Understanding how phenotypic differences between males and females arise from the sex-biased expression of nearly identical genomes can reveal important insights into the biology and evolution of a species. Among Anopheles mosquito species, these phenotypic differences include vectorial capacity, as it is only females that blood feed and thus transmit human malaria. Here, we use RNA-seq data from multiple tissues of four vector species spanning the Anopheles phylogeny to explore the genomic and evolutionary properties of sex-biased genes. We find that, in these mosquitoes, in contrast to what has been found in many other organisms, female-biased genes are more rapidly evolving in sequence, expression, and genic turnover than male-biased genes. Our results suggest that this atypical pattern may be due to the combination of sex-specific life history challenges encountered by females, such as blood feeding. Furthermore, female propensity to mate only once in nature in male swarms likely diminishes sexual selection of post-reproductive traits related to sperm competition among males. We also develop a comparative framework to systematically explore tissue- and sex-specific splicing to document its conservation throughout the genus and identify a set of candidate genes for future functional analyses of sex-specific isoform usage. Finally, our data reveal that the deficit of male-biased genes on the X Chromosomes in Anopheles is a conserved feature in this genus and can be directly attributed to chromosome-wide transcriptional regulation that de-masculinizes the X in male reproductive tissues.

The creation and selection of mutations resistant to a gene drive over multiple generations in the malaria mosquito.

Gene drives have enormous potential for the control of insect populations of medical and agricultural relevance. By preferentially biasing their own inheritance, gene drives can rapidly introduce genetic traits even if these confer a negative fitness effect on the population. We have recently developed gene drives based on CRISPR nuclease constructs that are designed to disrupt key genes essential for female fertility in the malaria mosquito. The construct copies itself and the associated genetic disruption from one homologous chromosome to another during gamete formation, a process called homing that ensures the majority of offspring inherit the drive. Such drives have the potential to cause long-lasting, sustainable population suppression, though they are also expected to impose a large selection pressure for resistance in the mosquito. One of these population suppression gene drives showed rapid invasion of a caged population over 4 generations, establishing proof of principle for this technology. In order to assess the potential for the emergence of resistance to the gene drive in this population we allowed it to run for 25 generations and monitored the frequency of the gene drive over time. Following the initial increase of the gene drive we observed a gradual decrease in its frequency that was accompanied by the spread of small, nuclease-induced mutations at the target gene that are resistant to further cleavage and restore its functionality. Such mutations showed rates of increase consistent with positive selection in the face of the gene drive. Our findings represent the first documented example of selection for resistance to a synthetic gene drive and lead to important design recommendations and considerations in order to mitigate for resistance in future gene drive applications.

Gene drive for population genetic control: non-functional resistance and parental effects.

Gene drive is a natural process of biased inheritance that, in principle, could be used to control pest and vector populations. As with any form of pest control, attention should be paid to the possibility of resistance evolving. For nuclease-based gene drive aimed at suppressing a population, resistance could arise by changes in the target sequence that maintain function, and various strategies have been proposed to reduce the likelihood that such alleles arise. Even if these strategies are successful, it is almost inevitable that alleles will arise at the target site that are resistant to the drive but do not restore function, and the impact of such sequences on the dynamics of control has been little studied. We use population genetic modelling of a strategy targeting a female fertility gene to demonstrate that such alleles may be expected to accumulate, and thereby reduce the reproductive load on the population, if nuclease expression per se causes substantial heterozygote fitness effects or if parental (especially paternal) deposition of nuclease either reduces offspring fitness or affects the genotype of their germline. All these phenomena have been observed in synthetic drive constructs. It will, therefore, be important to allow for non-functional resistance alleles in predicting the dynamics of constructs in cage populations and the impacts of any field release.

Crystallographic analyses illustrate significant plasticity and efficient recoding of meganuclease target specificity.

The retargeting of protein-DNA specificity, outside of extremely modular DNA binding proteins such as TAL effectors, has generally proved to be quite challenging. Here, we describe structural analyses of five different extensively retargeted variants of a single homing endonuclease, that have been shown to function efficiently in ex vivo and in vivo applications. The redesigned proteins harbor mutations at up to 53 residues (18%) of their amino acid sequence, primarily distributed across the DNA binding surface, making them among the most significantly reengineered ligand-binding proteins to date. Specificity is derived from the combined contributions of DNA-contacting residues and of neighboring residues that influence local structural organization. Changes in specificity are facilitated by the ability of all those residues to readily exchange both form and function. The fidelity of recognition is not precisely correlated with the fraction or total number of residues in the protein-DNA interface that are actually involved in DNA contacts, including directional hydrogen bonds. The plasticity of the DNA-recognition surface of this protein, which allows substantial retargeting of recognition specificity without requiring significant alteration of the surrounding protein architecture, reflects the ability of the corresponding genetic elements to maintain mobility and persistence in the face of genetic drift within potential host target sites.

Radical remodeling of the Y chromosome in a recent radiation of malaria mosquitoes.

Y chromosomes control essential male functions in many species, including sex determination and fertility. However, because of obstacles posed by repeat-rich heterochromatin, knowledge of Y chromosome sequences is limited to a handful of model organisms, constraining our understanding of Y biology across the tree of life. Here, we leverage long single-molecule sequencing to determine the content and structure of the nonrecombining Y chromosome of the primary African malaria mosquito, Anopheles gambiae We find that the An. gambiae Y consists almost entirely of a few massively amplified, tandemly arrayed repeats, some of which can recombine with similar repeats on the X chromosome. Sex-specific genome resequencing in a recent species radiation, the An. gambiae complex, revealed rapid sequence turnover within An. gambiae and among species. Exploiting 52 sex-specific An. gambiae RNA-Seq datasets representing all developmental stages, we identified a small repertoire of Y-linked genes that lack X gametologs and are not Y-linked in any other species except An. gambiae, with the notable exception of YG2, a candidate male-determining gene. YG2 is the only gene conserved and exclusive to the Y in all species examined, yet sequence similarity to YG2 is not detectable in the genome of a more distant mosquito relative, suggesting rapid evolution of Y chromosome genes in this highly dynamic genus of malaria vectors. The extensive characterization of the An. gambiae Y provides a long-awaited foundation for studying male mosquito biology, and will inform novel mosquito control strategies based on the manipulation of Y chromosomes.

Deciphering the olfactory repertoire of the tiger mosquito Aedes albopictus.

BACKGROUND: The Asian tiger mosquito Aedes albopictus is a highly invasive species and competent vector of several arboviruses (e.g. dengue, chikungunya, Zika) and parasites (e.g. dirofilaria) of public health importance. Compared to other mosquito species, Ae. albopictus females exhibit a generalist host seeking as well as a very aggressive biting behaviour that are responsible for its high degree of nuisance. Several complex mosquito behaviours such as host seeking, feeding, mating or oviposition rely on olfactory stimuli that target a range of sensory neurons localized mainly on specialized head appendages such as antennae, maxillary palps and the mouthparts. RESULTS: With the aim to describe the Ae. albopictus olfactory repertoire we have used RNA-seq to reveal the transcriptome profiles of female antennae and maxillary palps. Male heads and whole female bodies were employed as reference for differential expression analysis. The relative transcript abundance within each tissue (TPM, transcripts per kilobase per million) and the pairwise differential abundance in the different tissues (fold change values and false discovery rates) were evaluated. Contigs upregulated in the antennae (620) and maxillary palps (268) were identified and relative GO and PFAM enrichment profiles analysed. Chemosensory genes were described: overall, 77 odorant binding proteins (OBP), 82 odorant receptors (OR), 60 ionotropic receptors (IR) and 30 gustatory receptors (GR) were identified by comparative genomics and transcriptomics. In addition, orthologs of genes expressed in the female/male maxillary palps and/or antennae and involved in thermosensation (e.g. pyrexia and arrestin1), mechanosensation (e.g. piezo and painless) and neuromodulation were classified. CONCLUSIONS: We provide here the first detailed transcriptome of the main Ae. albopictus sensory appendages, i.e. antennae and maxillary palps. A deeper knowledge of the olfactory repertoire of the tiger mosquito will help to better understand its biology and may pave the way to design new attractants/repellents.

Reprogramming homing endonuclease specificity through computational design and directed evolution.

Homing endonucleases (HEs) can be used to induce targeted genome modification to reduce the fitness of pathogen vectors such as the malaria-transmitting Anopheles gambiae and to correct deleterious mutations in genetic diseases. We describe the creation of an extensive set of HE variants with novel DNA cleavage specificities using an integrated experimental and computational approach. Using computational modeling and an improved selection strategy, which optimizes specificity in addition to activity, we engineered an endonuclease to cleave in a gene associated with Anopheles sterility and another to cleave near a mutation that causes pyruvate kinase deficiency. In the course of this work we observed unanticipated context-dependence between bases which will need to be mechanistically understood for reprogramming of specificity to succeed more generally.

The germline of the malaria mosquito produces abundant miRNAs, endo-siRNAs, piRNAs and 29-nt small RNAs.

BACKGROUND: Small RNAs include different classes essential for endogenous gene regulation and cellular defence against genomic parasites. However, a comprehensive analysis of the small RNA pathways in the germline of the mosquito Anopheles gambiae has never been performed despite their potential relevance to reproductive capacity in this malaria vector. RESULTS: We performed small RNA deep sequencing during larval and adult gonadogenesis and find that they predominantly express four classes of regulatory small RNAs. We identified 45 novel miRNA precursors some of which were sex-biased and gonad-enriched , nearly doubling the number of previously known miRNA loci. We also determine multiple genomic clusters of 24-30 nt Piwi-interacting RNAs (piRNAs) that map to transposable elements (TEs) and 3'UTR of protein coding genes. Unusually, many TEs and the 3'UTR of some endogenous genes produce an abundant peak of 29-nt small RNAs with piRNA-like characteristics. Moreover, both sense and antisense piRNAs from TEs in both Anopheles gambiae and Drosophila melanogaster reveal novel features of piRNA sequence bias. We also discovered endogenous small interfering RNAs (endo-siRNAs) that map to overlapping transcripts and TEs. CONCLUSIONS: This is the first description of the germline miRNome in a mosquito species and should prove a valuable resource for understanding gene regulation that underlies gametogenesis and reproductive capacity. We also provide the first evidence of a piRNA pathway that is active against transposons in the germline and our findings suggest novel piRNA sequence bias. The contribution of small RNA pathways to germline TE regulation and genome defence in general is an important finding for approaches aimed at manipulating mosquito populations through the use of selfish genetic elements.

Measuring the Impact of Genetic Heterogeneity and Chromosomal Inversions on the Efficacy of CRISPR-Cas9 Gene Drives in Different Strains of Anopheles gambiae.

The human malaria vector Anopheles gambiae is becoming increasingly resistant to insecticides, spurring the development of genetic control strategies. CRISPR-Cas9 gene drives can modify a population by creating double-stranded breaks at highly specific targets, triggering copying of the gene drive into the cut site ("homing"), ensuring its inheritance. The DNA repair mechanism responsible requires homology between the donor and recipient chromosomes, presenting challenges for the invasion of laboratory-developed gene drives into wild populations of target species An. gambiae species complex, which show high levels of genome variation. Two gene drives (vas2-5958 and zpg-7280) were introduced into three An. gambiae strains collected across Africa with 5.3-6.6% variation around the target sites, and the effect of this variation on homing was measured. Gene drive homing across different karyotypes of the 2La chromosomal inversion was also assessed. No decrease in gene drive homing was seen despite target site heterology, demonstrating the applicability of gene drives to wild populations.

Single-cell profiling of Anopheles gambiae spermatogenesis defines the onset of meiotic silencing and premeiotic overexpression of the X chromosome.

Understanding development and genetic regulation in the Anopheles gambiae germline is essential to engineer effective genetic control strategies targeting this malaria mosquito vector. These include targeting the germline to induce sterility or using regulatory sequences to drive transgene expression for applications such as gene drive. However, only very few germline-specific regulatory elements have been characterised with the majority showing leaky expression. This has been shown to considerably reduce the efficiency of current genetic control strategies, which rely on regulatory elements with more tightly restricted spatial and/or temporal expression. Meiotic silencing of the sex chromosomes limits the flexibility of transgene expression to develop effective sex-linked genetic control strategies. Here, we build on our previous study, dissecting gametogenesis into four distinct cell populations, using single-cell RNA sequencing to define eight distinct cell clusters and associated germline cell-types using available marker genes. We reveal overexpression of X-linked genes in a distinct cluster of pre-meiotic cells and document the onset of meiotic silencing of the X chromosome in a subcluster of cells in the latter stages of spermatogenesis. This study provides a comprehensive dataset, characterising the expression of distinct cell types through spermatogenesis and widening the toolkit for genetic control of malaria mosquitoes.

CRISPR-Mediated Cassette Exchange (CriMCE): A Method to Introduce and Isolate Precise Marker-Less Edits.

The introduction of small unmarked edits to the genome of insects is essential to study the molecular underpinnings of important biological traits, such as resistance to insecticides and genetic control strategies. Advances in CRISPR genome engineering have made this possible, but prohibitively laborious for most laboratories due to low rates of editing and the lack of a selectable marker. To facilitate the generation and isolation of precise marker-less edits we have developed a two-step method based on CRISPR-mediated cassette exchange (CriMCE) of a marked placeholder for a variant of interest. This strategy can be used to introduce a wider range of potential edits compared with previous approaches while consolidating the workflow. We present proof-of-principle that CriMCE is a powerful tool by engineering three single nucleotide polymorphism variants into the genome of Anopheles gambiae, with 5-41 × higher rates of editing than homology-directed repair or prime editing.

A comprehensive gene expression atlas of sex- and tissue-specificity in the malaria vector, Anopheles gambiae.

BACKGROUND: The mosquito, Anopheles gambiae, is the primary vector of human malaria, a disease responsible for millions of deaths each year. To improve strategies for controlling transmission of the causative parasite, Plasmodium falciparum, we require a thorough understanding of the developmental mechanisms, physiological processes and evolutionary pressures affecting life-history traits in the mosquito. Identifying genes expressed in particular tissues or involved in specific biological processes is an essential part of this process. RESULTS: In this study, we present transcription profiles for ~82% of annotated Anopheles genes in dissected adult male and female tissues. The sensitivity afforded by examining dissected tissues found gene activity in an additional 20% of the genome that is undetected when using whole-animal samples. The somatic and reproductive tissues we examined each displayed patterns of sexually dimorphic and tissue-specific expression. By comparing expression profiles with Drosophila melanogaster we also assessed which genes are well conserved within the Diptera versus those that are more recently evolved. CONCLUSIONS: Our expression atlas and associated publicly available database, the MozAtlas (http://www.tissue-atlas.org), provides information on the relative strength and specificity of gene expression in several somatic and reproductive tissues, isolated from a single strain grown under uniform conditions. The data will serve as a reference for other mosquito researchers by providing a simple method for identifying where genes are expressed in the adult, however, in addition our resource will also provide insights into the evolutionary diversity associated with gene expression levels among species.

Developing transgenic Anopheles mosquitoes for the sterile insect technique.

In the last 10 years the availability of the genome sequence of Anopheles gambiae and the development of a transgenic technology for several species of Anopheles mosquitoes have, in combination, helped in enabling us to gain several insights into the biology of these mosquitoes that is relevant to their capacity as vectors of the malaria parasite. While this information is anticipated to inform many novel vector control strategies, the technique most likely to benefit in the near future from the availability of a reliable transgenic technology is the sterile insect technique (SIT), which relies on releasing large numbers of sterile insects to compete for mates in the wild, leading to population suppression. Although SIT has been proven to work reliably for many insects, the construction of suitable strains, and induction of sterility, has until now been a laborious process, combining classical genetics with radiation-induced sterility. Using transgenesis to create strains of Anopheles suitable for SIT could potentially offer several advantages over current approaches, in that the basic design of transgenic constructs designed for other insects should be rapidly transferable to mosquitoes, and induction of sterility as a product of the transgenic modification could obviate the requirement for radiation and its associated deleterious effects. In this paper the progress of different transgenic approaches in constructing tools for SIT will be reviewed.