Mosquito Gene Expression Research Articles

Climate change and vector-borne diseases: a multi-omics approach of temperature-induced changes in the mosquito.

Climate change and globalization contribute to the expansion of mosquito vectors and their associated pathogens. Long spared, temperate regions have had to deal with the emergence of arboviruses traditionally confined to tropical regions. Chikungunya virus (CHIKV) was reported for the first time in Europe in 2007, causing a localized outbreak in Italy, which then recurred repeatedly over the years in other European localities. This raises the question of climate effects, particularly temperature, on the dynamics of vector-borne viruses. The objective of this study is to improve the understanding of the molecular mechanisms set up in the vector in response to temperature. We combine three complementary approaches by examining Aedes albopictus mosquito gene expression (transcriptomics), bacterial flora (metagenomics) and CHIKV evolutionary dynamics (genomics) induced by viral infection and temperature changes. We show that temperature alters profoundly mosquito gene expression, bacterial microbiome and viral population diversity. We observe that (i) CHIKV infection upregulated most genes (mainly in immune and stress-related pathways) at 20°C but not at 28°C, (ii) CHIKV infection significantly increased the abundance of Enterobacteriaceae Serratia marcescens at 28°C and (iii) CHIKV evolutionary dynamics were different according to temperature. The substantial changes detected in the vectorial system (the vector and its bacterial microbiota, and the arbovirus) lead to temperature-specific adjustments to reach the ultimate goal of arbovirus transmission; at 20°C and 28°C, the Asian tiger mosquito Ae. albopictus was able to transmit CHIKV at the same efficiency. Therefore, CHIKV is likely to continue its expansion in the northern regions and could become a public health problem in more countries than those already affected in Europe.

Differences in gene expression in field populations of Wolbachia-infected Aedes aegypti mosquitoes with varying release histories in northern Australia.

Aedes aegypti is the principal mosquito vector of dengue, yellow fever, Zika and chikungunya viruses. The wMel strain of the endosymbiotic bacteria Wolbachia pipientis was introduced into the vector as a novel biocontrol strategy to stop transmission of these viruses. Mosquitoes with Wolbachia have been released in the field in Northern Queensland, Australia since 2011, at various locations and over several years, with populations remaining stably infected. Wolbachia infection is known to alter gene expression in its mosquito host, but whether (and how) this changes over the long-term in the context of field releases remains unknown. We sampled mosquitoes from Wolbachia-infected populations with three different release histories along a time gradient and performed RNA-seq to investigate gene expression changes in the insect host. We observed a significant impact on gene expression in Wolbachia-infected mosquitoes versus uninfected controls. Fewer genes had significantly upregulated expression in mosquitoes from the older releases (512 and 486 from the 2011 and 2013/14 release years, respectively) versus the more recent releases (1154 from the 2017 release year). Nonetheless, a fundamental signature of Wolbachia infection on host gene expression was observed across all releases, comprising upregulation of immunity (e.g. leucine-rich repeats, CLIPs) and metabolism (e.g. lipid metabolism, iron transport) genes. There was limited downregulation of gene expression in mosquitoes from the older releases (84 and 71 genes from the 2011 and 2013/14 release years, respectively), but significantly more in the most recent release (509 from the 2017 release year). Our findings indicate that at > 8 years post-introgression into field populations, Wolbachia continues to profoundly impact expression of host genes, such as those involved in insect immune response and metabolism. If Wolbachia-mediated virus blocking is underpinned by these differential gene expression changes, our results suggest it may remain stable long-term.

Profiling of Transcriptome-Wide N6-Methyladenosine (m6A) Modifications and Identifying m6A Associated Regulation in Sperm Tail Formation in Anopheles sinensis.

Recent discoveries of reversible N6-methyladenosine (m6A) methylation on messenger RNA (mRNA) and mapping of m6A methylomes in many species have revealed potential regulatory functions of this RNA modification by m6A players—writers, readers, and erasers. Here, we first profile transcriptome-wide m6A in female and male Anopheles sinensis and reveal that m6A is also a highly conserved modification of mRNA in mosquitoes. Distinct from mammals and yeast but similar to Arabidopsis thaliana, m6A in An. sinensis is enriched not only around the stop codon and within 3′-untranslated regions but also around the start codon and 5′-UTR. Gene ontology analysis indicates the unique distribution pattern of m6A in An. sinensis is associated with mosquito sex-specific pathways such as tRNA wobble uridine modification and phospholipid-binding in females, and peptidoglycan catabolic process, exosome and signal recognition particle, endoplasmic reticulum targeting, and RNA helicase activity in males. The positive correlation between m6A deposition and mRNA abundance indicates that m6A can play a role in regulating gene expression in mosquitoes. Furthermore, many spermatogenesis-associated genes, especially those related to mature sperm flagellum formation, are positively modulated by m6A methylation. A transcriptional regulatory network of m6A in An. sinensis is first profiled in the present study, especially in spermatogenesis, which may provide a new clue for the control of this disease-transmitting vector.

Temperature modulates immune gene expression in mosquitoes during arbovirus infection.

The principal vector of dengue, Zika and chikungunya viruses is the mosquito Aedes aegypti, with its ability to transmit pathogens influenced by ambient temperature. We use chikungunya virus (CHIKV) to understand how the mosquito transcriptome responds to arbovirus infection at different ambient temperatures. We exposed CHIKV-infected mosquitoes to 18, 28 and 32°C, and found that higher temperature correlated with higher virus levels, particularly at 3 days post infection, but lower temperature resulted in reduced virus levels. RNAseq analysis indicated significantly altered gene expression levels in CHIKV infection. The highest number of significantly differentially expressed genes was observed at 28°C, with a more muted effect at the other temperatures. At the higher temperature, the expression of many classical immune genes, including Dicer-2, was not substantially altered in response to CHIKV. The upregulation of Toll, IMD and JAK-STAT pathways was only observed at 28°C. Functional annotations suggested that genes in immune response and metabolic pathways related to energy supply and DNA replication were involved in temperature-dependent changes. Time post infection also led to substantially different gene expression profiles, and this varied with temperature. In conclusion, temperature significantly modulates mosquito gene expression in response to infection, potentially leading to impairment of immune defences at higher temperatures.

Temperature Dramatically Shapes Mosquito Gene Expression With Consequences for Mosquito-Zika Virus Interactions.

Vector-borne flaviviruses are emerging threats to human health. For successful transmission, the virus needs to efficiently enter mosquito cells and replicate within and escape several tissue barriers while mosquitoes elicit major transcriptional responses to flavivirus infection. This process will be affected not only by the specific mosquito-pathogen pairing but also by variation in key environmental variables such as temperature. Thus far, few studies have examined the molecular responses triggered by temperature and how these responses modify infection outcomes, despite substantial evidence showing strong relationships between temperature and transmission in a diversity of systems. To define the host transcriptional changes associated with temperature variation during the early infection process, we compared the transcriptome of mosquito midgut samples from mosquitoes exposed to Zika virus (ZIKV) and non-exposed mosquitoes housed at three different temperatures (20, 28, and 36°C). While the high-temperature samples did not show significant changes from those with standard rearing conditions (28°C) 48 h post-exposure, the transcriptome profile of mosquitoes housed at 20°C was dramatically different. The expression of genes most altered by the cooler temperature involved aspects of blood-meal digestion, ROS metabolism, and mosquito innate immunity. Further, we did not find significant differences in the viral RNA copy number between 24 and 48 h post-exposure at 20°C, suggesting that ZIKV replication is limited by cold-induced changes to the mosquito midgut environment. In ZIKV-exposed mosquitoes, vitellogenin, a lipid carrier protein, was most up-regulated at 20°C. Our results provide a deeper understanding of the temperature-triggered transcriptional changes in Aedes aegypti and can be used to further define the molecular mechanisms driven by environmental temperature variation.

Zika virus noncoding RNA suppresses apoptosis and is required for virus transmission by mosquitoes

Flaviviruses, including Zika virus (ZIKV), utilise host mRNA degradation machinery to produce subgenomic flaviviral RNA (sfRNA). In mammalian hosts, this noncoding RNA facilitates replication and pathogenesis of flaviviruses by inhibiting IFN-signalling, whereas the function of sfRNA in mosquitoes remains largely elusive. Herein, we conduct a series of in vitro and in vivo experiments to define the role of ZIKV sfRNA in infected Aedes aegypti employing viruses deficient in production of sfRNA. We show that sfRNA-deficient viruses have reduced ability to disseminate and reach saliva, thus implicating the role for sfRNA in productive infection and transmission. We also demonstrate that production of sfRNA alters the expression of mosquito genes related to cell death pathways, and prevents apoptosis in mosquito tissues. Inhibition of apoptosis restored replication and transmission of sfRNA-deficient mutants. Hence, we propose anti-apoptotic activity of sfRNA as the mechanism defining its role in ZIKV transmission.

One-step RT-qPCR assay for ZIKV RNA detection in Aedes aegypti samples: a protocol to study infection and gene expression during ZIKV infection

BackgroundZika virus (ZIKV) is transmitted to humans during the bite of an infected mosquito. In a scenario of globalization and climate change, the frequency of outbreaks has and will increase in areas with competent vectors, revealing a need for continuous improvement of ZIKV detection tools in vector populations. A simple, rapid and sensitive assay for viral detection is quantitative reverse transcription polymerase chain reaction (qRT-PCR), yet oligos optimized for ZIKV detection in mammalian cells and samples have repeatedly shown high background when used on mosquito ribonucleic acid (RNA). In this paper, we present a one-step qRT-PCR protocol that allows for the detection of ZIKV in mosquitoes and for the evaluation of gene expression from the same mosquito sample and RNA. This assay is a less expensive qRT-PCR approach than that most frequently used in the literature and has a much lower background, allowing confident detection.MethodsOur new oligo design to detect ZIKV RNA included in silico analysis of both viral and mosquito (Ae. aegypti and Ae. albopictus) genomes, targeting sequences conserved between Asian and African ZIKV lineages, but not matching Aedes genomes. This assay will allow researchers to avoid nonspecific amplification in insect samples due to viral integration into the mosquito genome, a phenomenon known to happen in wild and colonized populations of mosquitoes. Standard curves constructed with in vitro transcribed ZIKV RNA were used to optimize the sensitivity, efficiency and reproducibility of the assay.ResultsFinally, the assay was used with success to detect both ZIKV RNA in infected mosquitoes and to detect expression of the Defensin A gene, an antimicrobial peptide (AMP) involved in Aedes aegypti immune response to virus infection.ConclusionsThe experimental approach to detect ZIKV RNA in Aedes aegypti presented here has demonstrated to be specific, sensitive and reliable, and additionally it allows for the analysis of mosquito gene expression during ZIKV infection.

The coming of age of insulin-signaling in insects.

The coming of age of insulin-signaling in insects.

Wolbachia-a foe for mosquitoes

Mosquitoes act as vectors for a wide range of viral and parasitic infectious diseases such as malaria, dengue, Chickungunya, lymphatic filariasis, Japanese encephalitis and West Nile virus in humans as well as in animals. Although a wide range of insecticides are used to control mosquitoes, it has only resulted in development of resistance to such insecticides. The evolution of insecticide resistance and lack of vaccines for many mosquito-borne diseases have made these arthropods highly harmful vectors. Recently, a novel approach to control mosquitoes by transinfection of life shortening maternally transmitted endo-symbiont Wolbachia wMelPop strain from fruitfly Drosophila into mosquito population has been developed by researchers. The wMelPop strain up-regulated the immune gene expression in mosquitoes thereby reducing the dengue and Chickungunya viral replication in Aedes aegypti, and also it significantly reduced the Plasmodium level in Anopheles gambiae. Here, we discuss the strategy of using Wolbachia in control of vector-borne diseases of mosquitoes.

Juvenile hormone and its receptor, methoprene-tolerant, control the dynamics of mosquito gene expression

Juvenile hormone III (JH) plays a key role in regulating the reproduction of female mosquitoes. Microarray time-course analysis revealed dynamic changes in gene expression during posteclosion (PE) development in the fat body of female Aedes aegypti. Hierarchical clustering identified three major gene clusters: 1,843 early-PE (EPE) genes maximally expressed at 6 h PE, 457 mid-PE (MPE) genes at 24 h PE, and 1,815 late-PE (LPE) genes at 66 h PE. The RNAi microarray screen for the JH receptor Methoprene-tolerant (Met) showed that 27% of EPE and 40% of MPE genes were up-regulated whereas 36% of LPE genes were down-regulated in the absence of this receptor. Met repression of EPE and MPE and activation of LPE genes were validated by an in vitro fat-body culture experiment using Met RNAi. Sequence motif analysis revealed the consensus for a 9-mer Met-binding motif, CACG(C)/TG(A)/G(T)/AG. Met-binding motif variants were overrepresented within the first 300 bases of the promoters of Met RNAi-down-regulated (LPE) genes but not in Met RNAi-up-regulated (EPE) genes. EMSAs using a combination of mutational and anti-Met antibody supershift analyses confirmed the binding properties of the Met consensus motif variants. There was a striking temporal separation of expression profiles among major functional gene groups, with carbohydrate, lipid, and xenobiotics metabolism belonging to the EPE and MPE clusters and transcription and translation to the LPE cluster. This study represents a significant advancement in the understanding of the regulation of gene expression by JH and its receptor Met during female mosquito reproduction.

Model-Based Clustering of Longitudinal Data: Application to Modeling Disease Course and Gene Expression Trajectories

We consider the problem of clustering time-dependent data. The model is a mixture of regressions, with variance–covariance matrices that are allowed to vary within the extended linear mixed model family. We discuss applications to biomedical data and analyze two longitudinal data sets: one on patients with delirium, and the other on mosquito gene expression following infection.

Transgenic expression of the Aedes aegypti CYP9J28 confers pyrethroid resistance in Drosophila melanogaster

The emergence and spread of insecticide resistance in mosquitoes, such as the major vector of dengue and yellow fever Aedes aegypti, is a major public health problem.A number of studies have been conducted to-date aiming to identify specific molecular changes that are associated with the phenotype, in an attempt to use this knowledge to control the development and spread of resistant populations. Several approaches, ranging from classical techniques to high-throughput transcriptomics, have revealed a large number of genes which are associated with the trait. However, the development of functional tools for the validation of the possible role of those genes (if any) in mosquito resistance has not kept similar pace. Here, we report the use of Drosophila melanogaster as a model for establishing the role of A. aegypti cytochrome P450s. We showed that the transgenic expression of the CYP9J28 a gene which is upregulated in pyrethroid resistant A. aegypti mosquitoes and encodes an enzyme capable of metabolizing pyrethroids, in D. melanogaster confers significant levels of pyrethroid resistance. Our study provides solid functional in vivo evidence for the role of the AaegCYP9J28 in pyrethroid resistance and shows for the first time that ectopic expression of mosquito genes in D. melanogaster can be a robust approach for validating candidate resistance genes. Further potential applications of the approach are discussed.

Alterations in the Aedes aegypti Transcriptome during Infection with West Nile, Dengue and Yellow Fever Viruses

West Nile (WNV), dengue (DENV) and yellow fever (YFV) viruses are (re)emerging, mosquito-borne flaviviruses that cause human disease and mortality worldwide. Alterations in mosquito gene expression common and unique to individual flaviviral infections are poorly understood. Here, we present a microarray analysis of the Aedes aegypti transcriptome over time during infection with DENV, WNV or YFV. We identified 203 mosquito genes that were ≥5-fold differentially up-regulated (DUR) and 202 genes that were ≥10-fold differentially down-regulated (DDR) during infection with one of the three flaviviruses. Comparative analysis revealed that the expression profile of 20 DUR genes and 15 DDR genes was quite similar between the three flaviviruses on D1 of infection, indicating a potentially conserved transcriptomic signature of flaviviral infection. Bioinformatics analysis revealed changes in expression of genes from diverse cellular processes, including ion binding, transport, metabolic processes and peptidase activity. We also demonstrate that virally-regulated gene expression is tissue-specific. The overexpression of several virally down-regulated genes decreased WNV infection in mosquito cells and Aedes aegypti mosquitoes. Among these, a pupal cuticle protein was shown to bind WNV envelope protein, leading to inhibition of infection in vitro and the prevention of lethal WNV encephalitis in mice. This work provides an extensive list of targets for controlling flaviviral infection in mosquitoes that may also be used to develop broad preventative and therapeutic measures for multiple flaviviruses.

Optimization of the Gal4-UAS system in an Anopheles gambiae cell line

The development of the bipartite Gal4-UAS system in Anopheles gambiae would improve the functional characterization of genes in this important malaria vector. Towards this aim, we used Gal4 driver plasmids to successfully activate expression of the reporter gene, luciferase, from UAS responder plasmids when cotransfected into an An. gambiae cell line. To optimize Gal4-regulated gene expression in mosquitoes, we compared the efficiency of a series of alternative Gal4 transactivators to drive reporter gene expression from responder plasmids incorporating different numbers of tandemly arrayed Gal4 binding sites or upstream activation sequences (UAS). The results indicated that the native Gal4 is only weakly active in these cells. Modified forms of Gal4, including those carrying minimal VP16 activation domains, as well as a deleted form of Gal4, give up to 20-fold greater activity than the native protein, when used in conjunction with a responder plasmid having 14 UAS repeats. The identification of Gal4-UAS vectors that are efficiently expressed in a mosquito cell line should facilitate the transfer of this versatile expression system to An. gambiae, and potentially to other insects of medical importance.

Dissecting gene expression in mosquito

Gene expression is known to vary extensively among tissues and between sexes. However, detailed descriptions of tissue- and sex-specific gene expression are available for only a few model organisms. A new study published in BMC Genomics presents such a data set for the mosquito, Anopheles gambiae, which is the vector of human malaria. In addition to providing a valuable resource for the community of mosquito researchers, the study allows comparative transcriptomic studies of dipteran insects to be extended over 250 million years of evolution, since the divergence of A. gambiae and Drosophila melanogaster.

Developmental and evolutionary basis for drought tolerance of the Anopheles gambiae embryo

During the evolution of the Diptera there is a dramatic modification of the embryonic ectoderm, whereby mosquitoes contain separate amnion and serosa lineages while higher flies such as Drosophila melanogaster contain a single amnioserosa. Whole-genome transcriptome assays were performed with isolated serosa from Anopheles gambiae embryos. These assays identified a large number of genes implicated in the production of the larval cuticle. In D. melanogaster, these genes are activated just once during embryogenesis, during late stages where they are used for the production of the larval cuticle. Evidence is presented that the serosal cells secrete a dedicated serosal cuticle, which protects A. gambiae embryos from desiccation. Detailed temporal microarray assays of mosquito gene expression profiles revealed that the cuticular genes display biphasic expression during A. gambiae embryogenesis, first in the serosa of early embryos and then again during late stages as seen in D. melanogaster. We discuss how evolutionary modifications in the well-defined dorsal–ventral patterning network led to the wholesale deployment of the cuticle biosynthesis pathway in early embryos of A. gambiae.

Comparative genomics allows the discovery of cis -regulatory elements in mosquitoes

The discovery and mapping of cis-regulatory elements is important for understanding regulation of gene transcription in mosquito vectors of human diseases. Genome sequence data are available for 3 species, Aedes aegypti, Anopheles gambiae, and Culex quinquefasciatus (Diptera: Culicidae), representing 2 subfamilies (Culicinae and Anophelinae) that are estimated to have diverged 145 to 200 million years ago. Comparative genomics tools were used to screen genomic DNA fragments located in the 5'-end flanking regions of orthologous genes. These analyses resulted in the identification of 137 sequences, designated "mosquito motifs," 7 to 9 nucleotides in length, representing 18 families of putative cis-regulatory elements conserved significantly among the 3 species when compared to the fruit fly, Drosophila melanogaster. Forty-one of the motifs were implicated previously in experiments as sites for binding transcription factors or functioning in the regulation of mosquito gene expression. Further analyses revealed associations between specific motifs and expression profiles, particularly in those genes that show increased or decreased mRNA abundance in females following a blood meal, and those accumulating transcription products exclusively or preferentially in the midgut, fat bodies, or ovaries. These results validate the methodology and support a relationship between the discovered motifs and the conservation of hematophagy in mosquitoes.

Life cycle transcriptome of the malaria mosquito Anopheles gambiae and comparison with the fruitfly Drosophila melanogaster

The African mosquito Anopheles gambiae is the major vector of human malaria. We report a genome-wide survey of mosquito gene expression profiles clustered temporally into developmental programs and spatially into adult tissue-specific patterns. Global expression analysis shows that genes that belong to related functional categories or that encode the same or functionally linked protein domains are associated with characteristic developmental programs or tissue patterns. Comparative analysis of our data together with data published from Drosophila melanogaster reveal an overall strong and positive correlation of developmental expression between orthologous genes. The degree of correlation varies, depending on association of orthologs with certain developmental programs or functional groups. Interestingly, the similarity of gene expression is not correlated with the coding sequence similarity of orthologs, indicating that expression profiles and coding sequences evolve independently. In addition to providing a comprehensive view of temporal and spatial gene expression during the A. gambiae life cycle, this large-scale comparative transcriptomic analysis has detected important evolutionary features of insect transcriptomes.

Sodium channel gene expression in mosquitoes, Aedes albopictus (S.)

Abstract A mosquito strain of Aedes albopictus, HAmAalG0, from Huntsville, Alabama, USA, showed a normal susceptibility and low tolerance to permethrin and resmethrin (pyrethroid insecticides) compared to a susceptible Ikaken strain, even though these pyrethroid insecticides have been used in the field for a long period of time in Alabama. Recently, we treated HAmAalG0 in the laboratory with permethrin for five generations and detected no significant change in the level of resistance to permethrin in the selected mosquitoes, HAmAalG5, compared with the parental strain HAmAalG0. We then examined the allelic expression at the L‐to‐F kdr site of the sodium channel gene in the Aedes mosquitoes to address our hypothesis that the L‐to‐F kdr mutation was not present in HAmAalG0 and HAmAalG5 mosquitoes. We found that every tested individual in Ikaken, HAmAalG0, and HAmAalG5 populations expressed a codon of CTA at the L‐to‐F kdr site encoding Leu, strongly corresponding to their susceptibility to insecticides.

Plasmodium yoelii: The effect of second blood meal and anti-sporozoite antibodies on development and gene expression in the mosquito vector, Anopheles stephensi

The sporogonic development of the malaria parasite takes place in the mosquito and a wide range of factors modulates it. Among those, the contents of the blood meal can influence the parasite development directly or indirectly through the mosquito response to the infection. We have studied the effect of a second blood meal in previously infected mosquitoes and the effect of anti-sporozoite immune serum on parasite development and mosquito response to the infection. The prevalence and intensity of infection and gene expression of both Plasmodium yoelii and Anopheles stephensi was analyzed. We verified that a second blood meal and its immune status interfere with parasite development and with Plasmodium and mosquito gene expression.