RNA-seq analyses of blood-induced changes in gene expression in the mosquito vector species, Aedes aegypti

Characterization of the adult Aedes aegypti early midgut peritrophic matrix proteome using LC-MS

The Aedes aegypti mosquito is the principal vector of arboviruses such as dengue, chikungunya, yellow fever, and Zika virus. These arboviruses are transmitted during adult female mosquito bloodfeeding. While these viruses must transverse the midgut to replicate, the blood meal must also reach the midgut to be digested, absorbed, or excreted, as aggregation of blood meal metabolites can be toxic to the female mosquito midgut. The midgut peritrophic matrix (PM), a semipermeable extracellular layer comprised of chitin fibrils, glycoproteins, and proteoglycans, is one such mechanism of protection for the mosquito midgut. However, this structure has not been characterized for adult female Ae. aegypti. We conducted a mass spectrometry based proteomic analysis to identify proteins that comprise or are associated with the adult female Ae. aegypti early midgut PM. Altogether, 474 unique proteins were identified, with 115 predicted as secreted. GO-term enrichment analysis revealed an abundance of serine-type proteases and several known and novel intestinal mucins. In addition, approximately 10% of the peptides identified corresponded to known salivary proteins, indicating Ae. aegypti mosquitoes extensively swallow their own salivary secretions. However, the physiological relevance of this remains unclear, and further studies are needed to determine PM proteins integral for midgut protection from blood meal derived toxicity and pathogen protection. Finally, we describe substantial discordance between previously described transcriptionally changes observed in the midgut in response to a bloodmeal and the presence of the corresponding protein in the PM. Data are available via ProteomeXchange with identifier PXD007627.

Assessment of Aedes albopictus reference genes for quantitative PCR at different stages of development

Members of the Aedes genus of mosquitoes are widely recognized as vectors of viral diseases. Ae.albopictus is its most invasive species, and are known to carry viruses such as Dengue, Chikugunya and Zika. Its emerging importance puts Ae.albopictus on the forefront of genetic interaction and evolution studies. However, a panel of suitable reference genes specific for this insect is as of now undescribed. Nine reference genes, namely ACT, eEF1-γ, eIF2α, PP2A, RPL32, RPS17, PGK1, ILK and STK were evaluated. Expression patterns of the candidate reference genes were observed in a total of seventeen sample types, separated by stage of development and age. Gene stability was inferred from obtained quantification data through three widely cited evaluation algorithms i.e. BestKeeper, geNorm, and NormFinder. No single gene showed a satisfactory degree of stability throughout all developmental stages. Therefore, we propose combinations of PGK and ILK for early embryos; RPL32 and RPS17 for late embryos, all four larval instars, and pupae samples; eEF1-γ with STK for adult males; eEF1-γ with RPS17 for non-blood fed females; and eEF1-γ with eIF2α for both blood-fed females and cell culture. The results from this study should be able to provide a more informed selection of normalizing genes during qPCR in Ae.albopictus.

Meta-Analysis of Aedes aegypti Expression Datasets: Comparing Virus Infection and Blood-Fed Transcriptomes to Identify Markers of Virus Presence

The mosquito Aedes aegypti (L.) is vector of several arboviruses including dengue, yellow fever, chikungunya, and more recently zika. Previous transcriptomic studies have been performed to elucidate altered pathways in response to viral infection. However, the intrinsic coupling between alimentation and infection were unappreciated in these studies. Feeding is required for the initial mosquito contact with the virus and these events are highly dependent. Addressing this relationship, we reinterrogated datasets of virus-infected mosquitoes with two different diet schemes (fed and unfed mosquitoes), evaluating the metabolic cross-talk during both processes. We constructed coexpression networks with the differentially expressed genes of these comparison: virus-infected versus blood-fed mosquitoes and virus-infected versus unfed mosquitoes. Our analysis identified one module with 110 genes that correlated with infection status (representing ~0.7% of the A. aegypti genome). Furthermore, we performed a machine-learning approach and summarized the infection status using only four genes (AAEL012128, AAEL014210, AAEL002477, and AAEL005350). While three of the four genes were annotated as hypothetical proteins, AAEL012128 gene is a membrane amino acid transporter correlated with viral envelope binding. This gene alone is able to discriminate all infected samples and thus should have a key role to discriminate viral infection in the A. aegypti mosquito. Moreover, validation using external datasets found this gene as differentially expressed in four transcriptomic experiments. Therefore, these genes may serve as a proxy of viral infection in the mosquito and the others 106 identified genes provides a framework to future studies.

Individual co-variation between viral RNA load and gene expression reveals novel host factors during early dengue virus infection of the Aedes aegypti midgut

Dengue virus (DENV) causes more human infections than any other mosquito-borne virus. The current lack of antiviral strategies has prompted genome-wide screens for host genes that are required for DENV infectivity. Earlier transcriptomic studies that identified DENV host factors in the primary vector Aedes aegypti used inbred laboratory colonies and/or pools of mosquitoes that erase individual variation. Here, we performed transcriptome sequencing on individual midguts in a field-derived Ae. aegypti population to identify new candidate host factors modulating DENV replication. We analyzed the transcriptomic data using an approach that accounts for individual co-variation between viral RNA load and gene expression. This approach generates a prediction about the agonist or antagonist effect of candidate genes on DENV replication based on the sign of the correlation between gene expression and viral RNA load. Using this method, we identified 39 candidate genes that went undetected by conventional pairwise comparison of gene expression levels between DENV-infected midguts and uninfected controls. Only four candidate genes were detected by both methods, emphasizing their complementarity. We demonstrated the value of our approach by functional validation of a candidate agonist gene encoding a sterol regulatory element-binding protein (SREBP), which was identified by correlation analysis but not by pairwise comparison. We confirmed that SREBP promotes DENV infection in the midgut by RNAi-mediated gene knockdown in vivo. We suggest that our approach for transcriptomic analysis can empower genome-wide screens for potential agonist or antagonist factors by leveraging inter-individual variation in gene expression. More generally, this method is applicable to a wide range of phenotypic traits displaying inter-individual variation.

Response to Blood Meal in the Fat Body of Anopheles stephensi Using Quantitative Proteomics: Toward New Vector Control Strategies Against Malaria

Malaria remains a grand challenge for disruptive innovation in global health therapeutics and diagnostics. Anopheles stephensi is one of the major vectors of malaria in Asia. Vector and transmission control are key focus areas in the fight against malaria, a field of postgenomics research where proteomics can play a substantive role. Moreover, to identify novel strategies to control the vector population, it is necessary to understand the vector life processes at a global and molecular scale. In this context, fat body is a vital organ required for vitellogenesis, vector immunity, vector physiology, and vector-parasite interaction. Given its central role in energy metabolism, vitellogenesis, and immune function, the proteome profile of the fat body and the impact of blood meal (BM) ingestion on the protein abundances of this vital organ have not been investigated so far. Therefore, using a proteomics approach, we identified the proteins expressed in the fat body of An. stephensi and their differential expression in response to BM ingestion. In all, we identified 3,218 proteins in the fat body using high-resolution mass spectrometry, of which 483 were found to be differentially expressed in response to the BM ingestion. Bioinformatics analysis of these proteins underscored their role in amino acid metabolism, vitellogenesis, lipid transport, signal peptide processing, mosquito immunity, and oxidation-reduction processes. Interestingly, we identified five novel genes, which were found to be differentially expressed upon BM ingestion. Proteins that exhibited altered expression in the present study are potential targets for vector control strategies and development of transmission blocking vaccines in the fight against malaria.

Impact of transmission cycles and vector competence on global expansion and emergence of arboviruses

Arboviruses are transmitted between arthropod vectors and vertebrate host. Arboviral infection in mosquitoes is initiated when a mosquito feeds on a viremic host. Following ingestion of a viremic blood meal by mosquitoes, virus enters midgut along with the blood, infects and replicates in midgut epithelial cells, and then escapes to the hemocoel, from where it disseminates to various secondary organs including salivary glands. Subsequently, when mosquito bites another host, a new transmission cycle is initiated. The midgut and salivary glands act as anatomical barriers to virus infection and escape. These complex interactions between the virus and vector dictate the vector competence. Thus, vector competence reflects the success in overcoming different barriers within the vector. Along with these, other intrinsic factors like midgut microbiota and immune responses, extrinsic factors like temperature and humidity, and genetic factors like vector genotype and viral genotype have been discussed in this review. Recent advancement on novel molecular tools to study vector competence is also included. Different modes of arboviral transmission like horizontal, vertical, and venereal and how these play role in sustenance and emergence of arboviruses in nature are also discussed. These factors can be exploited to reduce the susceptibility of vectors for the viruses, so as to control arboviral diseases to certain extent.

Implications of autophagy on arbovirus infection of mosquitoes

Arthropod-borne (arbo-) viruses, like all viruses, are obligate intracellular parasites that have evolved mechanisms to subvert cellular processes and evade anti-viral defenses to replicate and persist. An increasing body of research is beginning to recognize the intimate relationship between arboviruses and the cellular autophagy pathway. As a result, new therapeutic approaches that modify the autophagic response to viral infection have shown great promise. The preponderance of work thus far, however, has originated from vertebrate systems. Efforts to elucidate the role of autophagy during arbovirus infection of invertebrates have emerged, providing new insights into arbovirus-vector interactions; interactions that could be exploited for novel control strategies.

Amino acid metabolic signaling influences Aedes aegypti midgut microbiome variability

The mosquito midgut microbiota has been shown to influence vector competence for multiple human pathogens. The microbiota is highly variable in the field, and the sources of this variability are not well understood, which limits our ability to understand or predict its effects on pathogen transmission. In this work, we report significant variation in female adult midgut bacterial load between strains of A. aegypti which vary in their susceptibility to dengue virus. Composition of the midgut microbiome was similar overall between the strains, with 81-92% of reads coming from the same five bacterial families, though we did detect differences in the presence of some bacterial families including Flavobacteriaceae and Entobacteriaceae. We conducted transcriptomic analysis on the two mosquito strains that showed the greatest difference in bacterial load, and found that they differ in transcript abundance of many genes implicated in amino acid metabolism, in particular the branched chain amino acid degradation pathway. We then silenced this pathway by targeting multiple genes using RNA interference, which resulted in strain-specific bacterial proliferation, thereby eliminating the difference in midgut bacterial load between the strains. This suggests that the branched chain amino acid (BCAA) degradation pathway controls midgut bacterial load, though the mechanism underlying this remains unclear. Overall, our results indicate that amino acid metabolism can act to influence the midgut microbiota. Moreover, they suggest that genetic or physiological variation in BCAA degradation pathway activity may in part explain midgut microbiota variation in the field.

The Dose Makes the Poison: Nutritional Overload Determines the Life Traits of Blood-Feeding Arthropods

Vertebrate blood composition is heavily biased towards proteins, and hemoglobin, which is a hemeprotein, is by far the most abundant protein. Typically, hematophagous insects ingest blood volumes several times their weight before the blood meal. This barbarian feast offers an abundance of nutrients, but the degradation of blood proteins generates toxic concentrations of amino acids and heme, along with unparalleled microbiota growth. Despite this challenge, hematophagous arthropods have successfully developed mechanisms that bypass the toxicity of these molecules. While these adaptations allow hematophagous arthropods to tolerate their diet, they also constitute a unique mode of operation for cell signaling, immunity, and metabolism, the study of which may offer insights into the biology of disease vectors and may lead to novel vector-specific control methods.

Blood meal induced regulation of the chemosensory gene repertoire in the southern house mosquito

Background

The southern house mosquito, Culex quinquefasciatus, is one of the most prevalent vectors of lymphatic filariasis and flavivirus-induced encephalitis. Its vectorial capacity is directly affected by its reproductive feeding behaviors, such as host seeking, blood feeding, resting, and egg laying. In mosquitoes, these gonotrophic behaviors are odor-mediated and regulated following blood feeding. Immediately after a blood meal, female mosquitoes show reduced olfactory responsiveness and flight activity, as they enter a resting state. Insights into antennal chemosensory gene regulation at this time period can provide a foundation to identify targets involved in the state switch between host seeking and resting.

Results

This study used quantitative gene expression analyses to explore blood meal induced regulation of chemosensory gene families in the antennae of 6 days post-emergence C. quinquefasciatus females. Improved annotations for multiple chemosensory gene families, and a quantitative differential gene expression analysis between host seeking and 24 h post- blood fed females of the same age, allowed for the detection of transcripts that potentially play a role in the switch from host seeking to resting, in C. quinquefasciatus. The expression profiles of chemosensory genes varied significantly between the two treatments.

Conclusions

Annotations for chemosensory gene repertoires in C. quinquefasciatus have been manually curated and corrected for 3’ exon choice and transcript length, through sequence and transcriptome analyses. The gene expression analyses identified various molecular components of the peripheral olfactory system in C. quinquefasciatus, including odorant receptors, ionotropic receptors, odorant binding proteins and chemosensory proteins, that are regulated in response to blood feeding, and could be critical for the behavioral switch from host seeking to resting. Functional characterization of these proteins in the future can identify targets essential for the females’ gonotrophic behaviors, and can be used to design novel vector control strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3779-2) contains supplementary material, which is available to authorized users.

The midgut transcriptome of Aedes aegypti fed with saline or protein meals containing chikungunya virus reveals genes potentially involved in viral midgut escape

BACKGROUND

The mosquito Aedes aegypti is the primary vector for medically important arthropod-borne viruses, including chikungunya virus (CHIKV). Following oral acquisition, an arbovirus has to persistently infect several organs in the mosquito before becoming transmissible to another vertebrate host. A major obstacle an arbovirus has to overcome during its infection cycle inside the mosquito is the midgut escape barrier, representing the exit mechanism arboviruses utilize when disseminating from the midgut. To understand the transcriptomic basis of midgut escape and to reveal genes involved in the process, we conducted a comparative transcriptomic analysis of midgut samples from mosquitoes which had received a saline meal (SM) or a protein meal (PM) (not) containing CHIKV.

RESULTS

CHIKV which was orally acquired by a mosquito along with a SM or PM productively infected the midgut epithelium and disseminated to secondary tissues. A total of 27 RNA-Seq libraries from midguts of mosquitoes that had received PM or SM (not) containing CHIKV at 1 and 2 days post-feeding were generated and sequenced. Fewer than 80 genes responded differentially to the presence of CHIKV in midguts of mosquitoes that had acquired the virus along with SM or PM. SM feeding induced differential expression (DE) of 479 genes at day 1 and 314 genes at day 2 when compared to midguts of sugarfed mosquitoes. By comparison, PM feeding induced 6029 DE genes at day 1 and 7368 genes at day 2. Twenty-three DE genes encoding trypsins, metalloproteinases, and serine-type endopeptidases were significantly upregulated in midguts of mosquitoes at day 1 following SM or PM ingestion. Two of these genes were Ae. aegypti late trypsin (AeLT) and serine collagenase 1 precursor (AeSP1). In vitro, recombinant AeLT showed strong matrix metalloproteinase activity whereas recombinant AeSP1 did not.

CONCLUSIONS

By substituting a bloodmeal for SM, we identified midgut-expressed genes not involved in blood or protein digestion. These included genes coding for trypsins, metalloproteinases, and serine-type endopeptidases, which could be involved in facilitating midgut escape for arboviruses in Ae. aegypti. The presence of CHIKV in any of the ingested meals had relatively minor effects on the overall gene expression profiles in midguts.

RNA-Seq Comparison of Larval and Adult Malpighian Tubules of the Yellow Fever Mosquito Aedes aegypti Reveals Life Stage-Specific Changes in Renal Function

Introduction: The life history of Aedes aegypti presents diverse challenges to its diuretic system. During the larval and pupal life stages mosquitoes are aquatic. With the emergence of the adult they become terrestrial. This shifts the organism within minutes from an aquatic environment to a terrestrial environment where dehydration has to be avoided. In addition, female mosquitoes take large blood meals, which present an entirely new set of challenges to salt and water homeostasis.

Methods: To determine differences in gene expression associated with these different life stages, we performed an RNA-seq analysis of the main diuretic tissue in A. aegypti, the Malpighian tubules. We compared transcript abundance in 4th instar larvae to that of adult females and analyzed the data with a focus on transcripts that encode proteins potentially involved in diuresis, like water and solute channels as well as ion transporters. We compared our results against the model of potassium- and sodium chloride excretion in the Malpighian tubules proposed by Hine et al. (2014), which involves at least eight ion transporters and a proton-pump.

Results: We found 3,421 of a total number of 17,478 (19.6%) unique transcripts with a P < 0.05 and at least a 2.5 fold change in expression levels between the two groups. We identified two novel transporter genes that are highly expressed in the adult Malpighian tubules, which have not previously been part of the transport model in this species and may play important roles in diuresis. We also identified candidates for hypothesized sodium and chloride channels. Detoxification genes were generally higher expressed in larvae.

Significance: This study represents the first comparison of Malpighian tubule transcriptomes between larval and adult A. aegypti mosquitoes, highlighting key differences in their renal systems that arise as they transform from an aquatic filter-feeding larval stage to a terrestrial, blood-feeding adult stage.

The impact of metagenomic interplay on the mosquito redox homeostasis

Mosquitoes are exposed to oxidative challenges throughout their life cycle. The primary challenge comes from a blood meal. The blood digestion turns the midgut into an oxidative environment, which imposes pressure not only on mosquito fecundity and other physiological traits but also on the microbiota in the midgut. During evolution, mosquitoes have developed numerous oxidative defense mechanisms to maintain redox homeostasis in the midgut. In addition to antioxidants, SOD, catalase, and glutathione system, sufficient supply of the reducing agent, NADPH, is vital for a successful defense against oxidative stress. Increasing evidence indicates that in response to oxidative stress, cells reconfigure metabolic pathways to increase the generation of NADPH through NADP-reducing networks including the pentose phosphate pathway and others. The microbial homeostasis is critical for the functional contributions to various host phenotypes. The symbiotic microbiota is regulated largely by the Duox-ROS pathway in Drosophila. In mosquitoes, Duox-ROS pathway, heme-mediated signaling, antimicrobial peptide production and C-type lectins work in concert to maintain the dynamic microbial community in the midgut. Microbial mechanisms against oxidative stress in this context are not well understood. Emerging evidence that microbial metabolites trigger host oxidative response warrants further study on the metagenomic interplay in an oxidative environment like mosquito gut ecosystem. Besides the classical Drosophila model, hematophagous insects like mosquitoes provide an alternative model system to study redox homeostasis in a symbiotic metagenomic context.

Retention of duplicated long-wavelength opsins in mosquito lineages by positive selection and differential expression

Background

Opsins are light sensitive receptors associated with visual processes. Insects typically possess opsins that are stimulated by ultraviolet, short and long wavelength (LW) radiation. Six putative LW-sensitive opsins predicted in the yellow fever mosquito, Aedes aegypti and malaria mosquito, Anopheles gambiae, and eight in the southern house mosquito, Culex quinquefasciatus, suggest gene expansion in the Family Culicidae (mosquitoes) relative to other insects. Here we report the first detailed molecular and evolutionary analyses of LW opsins in three mosquito vectors, with a goal to understanding the molecular basis of opsin-mediated visual processes that could be exploited for mosquito control.

Results

Time of divergence estimates suggest that the mosquito LW opsins originated from 18 or 19 duplication events between 166.9/197.5 to 1.07/0.94 million years ago (MY) and that these likely occurred following the predicted divergence of the lineages Anophelinae and Culicinae 145–226 MY. Fitmodel analyses identified nine amino acid residues in the LW opsins that may be under positive selection. Of these, eight amino acids occur in the N and C termini and are shared among all three species, and one residue in TMIII was unique to culicine species. Alignment of 5′ non-coding regions revealed potential Conserved Non-coding Sequences (CNS) and transcription factor binding sites (TFBS) in seven pairs of LW opsin paralogs.

Conclusions

Our analyses suggest opsin gene duplication and residues possibly associated with spectral tuning of LW-sensitive photoreceptors. We explore two mechanisms - positive selection and differential expression mediated by regulatory units in CNS – that may have contributed to the retention of LW opsin genes in Culicinae and Anophelinae. We discuss the evolution of mosquito LW opsins in the context of major Earth events and possible adaptation of mosquitoes to LW-dominated photo environments, and implications for mosquito control strategies based on disrupting vision-mediated behaviors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0910-6) contains supplementary material, which is available to authorized users.

Apoptosis-related genes control autophagy and influence DENV-2 infection in the mosquito vector, Aedes aegypti

The mosquito Aedes aegypti is the primary urban vector for dengue virus (DENV) worldwide. Insight into interactions occurring between host and pathogen is important in understanding what factors contribute to vector competence. However, many of the molecular mechanisms for vector competence remain unknown. Our previous global transcriptional analysis suggested that differential expression of apoptotic proteins is involved in determining refractoriness vs susceptibility to DENV-2 infection in Ae. aegypti females following a DENV-infected blood meal. To determine whether DENV-refractory Ae. aegypti showed more robust apoptosis upon infection, we compared numbers of apoptotic cells from midguts of refractory and susceptible strains and observed increased numbers of apoptotic cells in only the refractory strain upon DENV-2 infection. Thereafter, we manipulated apoptosis through dsRNA interference of the initiator caspase, Aedronc. Unexpectedly, dsAedronc-treated females showed both decreased frequency of disseminated infection and decreased virus titer in infected individuals. Insect caspases have also previously been identified as regulators of the cellular recycling process known as autophagy. We observed activation of autophagy in midgut and fat body tissues following a blood meal, as well as programmed activation of several apoptosis-related genes, including the effector caspase, Casps7. To determine whether autophagy was affected by caspase knockdown, we silenced Aedronc and Casps7, and observed reduced activation of autophagy upon silencing. Our results provide evidence that apoptosis-related genes are also involved in regulating autophagy, and that Aedronc may play an important role in DENV-2 infection success in Ae. aegypti, possibly through its regulation of autophagy.

Comparative proteome analysis reveals that blood and sugar meals induce differential protein expression in Aedes aegypti female heads

Aedes aegypti females ingest sugar or blood to obtain the nutrients needed to maintain cellular homeostasis. During human blood ingestion, female mosquitoes may transmit different viruses such as dengue, yellow fever and, more recently, zika and chikungunya. Here, we report changes in protein expression in the heads of adult female Ae. aegypti mosquitoes in response to the ingestion of blood or sugar. Proteins extracted from the heads of Ae. aegypti fed exclusively on blood (BF) or sugar (SF) were trypsin hydrolyzed (off-gel) and analyzed by the reverse-phase nano-liquid chromatography coupled with hybrid mass spectrometry. A total of 1139 proteins were identified in female heads, representing 7.4% of the predicted proteins in Ae. aegypti genome (total = 15 419 active genes). Gene ontology annotation and categories showed that, in this insect, the head was rich in proteins involved in the metabolic process, proton transport, organelle, macromolecular complex, structural molecule activity, antioxidant activity, and catalytic activity. Our report is the first indicating that many of the annotated genes are translated into functional proteins in heads of adult female Ae. aegypti. Interestingly, we identified 8.7 times more exclusively expressed proteins involved in signal transduction, replication-transcription-translation (5.5 x), and transport (2.9 x) activity in BF than in SF groups. This paper discusses the protein profile of Ae. aegypti female heads and its implications for blood ingestion and carbohydrate intake.

Functional characterization of three trehalase genes regulating the chitin metabolism pathway in rice brown planthopper using RNA interference

RNA interference (RNAi) is an effective gene-silencing tool, and double stranded RNA (dsRNA) is considered a powerful strategy for gene function studies in insects. In the present study, we aimed to investigate the function of trehalase (TRE) genes (TRE 1-1, TRE 1-2, and TRE-2) isolated from the brown planthopper Nilaparvata lugens, a typical piercing-sucking insect in rice, and investigate their regulating roles in chitin synthesis by injecting larvae with dsRNA. The results showed that TRE1 and TRE2 had compensatory function, and the expression of each increased when the other was silenced. The total rate of insects with phenotypic deformities ranged from 19.83 to 24.36% after dsTRE injection, whereas the mortality rate ranged from 14.16 to 31.78%. The mRNA levels of genes involved in the chitin metabolism pathway in RNA-Seq and DGEP, namely hexokinase (HK), glucose-6-phosphate isomerase (G6PI) and chitinase (Cht), decreased significantly at 72 h after single dsTREs injection, whereas two transcripts of chitin synthase (CHS) genes decreased at 72 h after dsTRE1-1 and dsTREs injection. These results demonstrated that TRE silencing could affect the regulation of chitin biosynthesis and degradation, causing moulting deformities. Therefore, expression inhibitors of TREs might be effective tools for the control of planthoppers in rice.

Expression Profiles and RNAi Silencing of Inhibitor of Apoptosis Transcripts in Aedes, Anopheles, and Culex Mosquitoes (Diptera: Culicidae)

Effective mosquito control is vital to curtail the devastating health effects of many vectored diseases. RNA interference (RNAi)-mediated control of mosquitoes is an attractive alternative to conventional chemical pesticides. Previous studies have suggested that transcripts for inhibitors of apoptosis (IAPs) may be good RNAi targets. To revisit and extend previous reports, we examined the expression of Aedes aegypti (L.) IAPs (AaeIAPs) 1, 2, 5, 6, 9, and a viral IAP-associated factor (vIAF) as well as Anopheles quadrimaculatus Say and Culex quinquefasciatus Say IAP1 homologs (AquIAP1 and CquIAP1) in adult females. Expression profiles of IAPs suggested that some older female mosquitoes had significantly higher IAP mRNA levels when compared to the youngest ones. Minor differences in expression of AaeIAPs were observed in mosquitoes that imbibed a bloodmeal, but the majority of the time points (up to 48 h) were not significantly different. Although in vitro experiments with the Ae. aegypti Aag-2 cell line demonstrated that the various AaeIAPs could be effectively knocked down within one day after dsRNA treatment, only Aag-2 cells treated with dsIAP1 displayed apoptotic morphology. Gene silencing and mortality were also evaluated after topical application and microinjection of the same dsRNAs into female Ae. aegypti. In contrast to previous reports, topical administration of dsRNA against AaeIAP1 did not yield a significant reduction in gene expression or increased mortality. Knockdown of IAP1 and other IAPs by microinjection did not result in significant mortality. In toto, our findings suggest that IAPs may not be suitable RNAi targets for controlling adult mosquito populations.

Transcriptome of larvae representing the Rhipicephalus sanguineus complex

Rhipicephalus sanguineus sensu lato (Ixodida: Ixodidae) is possibly the most widespread tick species worldwide, responsible for transmitting several vector-borne pathogens of medical and veterinary importance. Here, we explore the transcriptome of R. sanguineus s.l. larvae (Putignano strain). We sequenced total RNA from R. sanguineus s.l. larvae. A total of 15,566,986 short paired-end reads were de novo-assembled into 33,396 transcripts and then annotated and analyzed. Particular attention was paid to transcripts putatively encoding ATP-binding proteins, due to their importance as mechanisms of detoxification and acaricide resistance. Additionally, microsatellite loci were investigated, as these are useful markers for population genetic studies. The present data and analyses provide a comprehensive transcriptomic resource for R. sanguineus. The results presented here will aid further genetic and genomic studies of this important tick species.

Plasticity in Insect Olfaction: To Smell or Not to Smell?

In insects, olfaction plays a crucial role in many behavioral contexts, such as locating food, sexual partners, and oviposition sites. To successfully perform such behaviors, insects must respond to chemical stimuli at the right moment. Insects modulate their olfactory system according to their physiological state upon interaction with their environment. Here, we review the plasticity of behavioral responses to different odor types according to age, feeding state, circadian rhythm, and mating status. We also summarize what is known about the underlying neural and endocrinological mechanisms, from peripheral detection to central nervous integration, and cover neuromodulation from the molecular to the behavioral level. We describe forms of olfactory plasticity that have contributed to the evolutionary success of insects and have provided them with remarkable tools to adapt to their ever-changing environment.

Arbovirus-mosquito interactions: RNAi pathway

Arthropod-borne (arbo) viruses infect hematophagous arthropods (vectors) to maintain virus transmission between vertebrate hosts. The mosquito vector actively controls arbovirus infection to minimize its fitness costs. The RNA interference (RNAi) pathway is the major antiviral response vectors use to restrict arbovirus infections. We know this because depleting RNAi gene products profoundly impacts arbovirus replication, the antiviral RNAi pathway genes undergo positive, diversifying selection and arboviruses have evolved strategies to evade the vector's RNAi responses. The vector's RNAi defense and arbovirus countermeasures lead to an arms race that prevents potential virus-induced fitness costs yet maintains arbovirus infections needed for transmission. This review will discuss the latest findings in RNAi-arbovirus interactions in the model insect (Drosophila melanogaster) and in specific mosquito vectors.

Production of Infectious Dengue Virus in Aedes aegypti Is Dependent on the Ubiquitin Proteasome Pathway

Dengue virus (DENV) relies on host factors to complete its life cycle in its mosquito host for subsequent transmission to humans. DENV first establishes infection in the midgut of Aedes aegypti and spreads to various mosquito organs for lifelong infection. Curiously, studies have shown that infectious DENV titers peak and decrease thereafter in the midgut despite relatively stable viral genome levels. However, the mechanisms that regulate this decoupling of infectious virion production from viral RNA replication have never been determined. We show here that the ubiquitin proteasome pathway (UPP) plays an important role in regulating infectious DENV production. Using RNA interference studies, we show in vivo that knockdown of selected UPP components reduced infectious virus production without altering viral RNA replication in the midgut. Furthermore, this decoupling effect could also be observed after RNAi knockdown in the head/thorax of the mosquito, which otherwise showed direct correlation between infectious DENV titer and viral RNA levels. The dependence on the UPP for successful DENV production is further reinforced by the observed up-regulation of key UPP molecules upon DENV infection that overcome the relatively low expression of these genes after a blood meal. Collectively, our findings indicate an important role for the UPP in regulating DENV production in the mosquito vector.

DENV is a re-emerging mosquito-borne flavivirus and the lack of sustainable preventative or therapeutic measures renders it a significant health burden globally. Although suppression of mosquito populations represents the most widely used dengue control strategy, there has been limited efficacy with this method. Blocking host factors required for DENV replication in mosquitoes may thus serve as an effective anti-transmission strategy. We have recently observed that the UPP plays a critical role in regulating DENV egress from infected cells, but how the UPP contributes to DENV life cycle in mosquitoes remain ill-defined. We show here that the Aedes aegypti midgut has evolved to control persistent DENV infection by differentially regulating key genes in the UPP, without harm to itself. Using RNAi, knockdown of proteasome subunits, β1, β2 and β5, as well as other UPP-specific genes in vivo reduced the production of infectious virus without altering viral RNA replication in the mosquito. Targeting these host factors via dsRNA-mediated or chemical inactivation in the mosquito salivary glands may serve as a viable anti-dengue transmission strategy from mosquitoes to humans.

Contrasting genetic structure between mitochondrial and nuclear markers in the dengue fever mosquito from Rio de Janeiro: implications for vector control

Dengue is the most prevalent global arboviral disease that affects over 300 million people every year. Brazil has the highest number of dengue cases in the world, with the most severe epidemics in the city of Rio de Janeiro (Rio). The effective control of dengue is critically dependent on the knowledge of population genetic structuring in the primary dengue vector, the mosquito Aedes aegypti. We analyzed mitochondrial and nuclear genomewide single nucleotide polymorphism markers generated via Restriction-site Associated DNA sequencing, as well as traditional microsatellite markers in Ae. aegypti from Rio. We found four divergent mitochondrial lineages and a strong spatial structuring of mitochondrial variation, in contrast to the overall nuclear homogeneity across Rio. Despite a low overall differentiation in the nuclear genome, we detected strong spatial structure for variation in over 20 genes that have a significantly altered expression in response to insecticides, xenobiotics, and pathogens, including the novel biocontrol agent Wolbachia. Our results indicate that high genetic diversity, spatially unconstrained admixing likely mediated by male dispersal, along with locally heterogeneous genetic variation that could affect insecticide resistance and mosquito vectorial capacity, set limits to the effectiveness of measures to control dengue fever in Rio.

Regulation of Gene Expression Patterns in Mosquito Reproduction

In multicellular organisms, development, growth and reproduction require coordinated expression of numerous functional and regulatory genes. Insects, in addition to being the most speciose animal group with enormous biological and economical significance, represent outstanding model organisms for studying regulation of synchronized gene expression due to their rapid development and reproduction. Disease-transmitting female mosquitoes have adapted uniquely for ingestion and utilization of the huge blood meal required for swift reproductive events to complete egg development within a 72-h period. We investigated the network of regulatory factors mediating sequential gene expression in the fat body, a multifunctional organ analogous to the vertebrate liver and adipose tissue, of the female Aedes aegypti mosquito. Transcriptomic and bioinformatics analyses revealed that ~7500 transcripts are differentially expressed in four sequential waves during the 72-h reproductive period. A combination of RNA-interference gene-silencing and in-vitro organ culture identified the major regulators for each of these waves. Amino acids (AAs) regulate the first wave of gene activation between 3 h and 12 h post-blood meal (PBM). During the second wave, between 12 h and 36 h, most genes are highly upregulated by a synergistic action of AAs, 20-hydroxyecdysone (20E) and the Ecdysone-Receptor (EcR). Between 36 h and 48 h, the third wave of gene activation—regulated mainly by HR3—occurs. Juvenile Hormone (JH) and its receptor Methoprene-Tolerant (Met) are major regulators for the final wave between 48 h and 72 h. Each of these key regulators also has repressive effects on one or more gene sets. Our study provides a better understanding of the complexity of the regulatory mechanisms related to temporal coordination of gene expression during reproduction. We have detected the novel function of 20E/EcR responsible for transcriptional repression. This study also reveals the previously unidentified large-scale effects of HR3 and JH/Met on transcriptional regulation during the termination of vitellogenesis and remodeling of the fat body.

In addition to being vectors of devastating human diseases, mosquitoes represent outstanding model organisms for studying regulatory mechanisms of differential gene expression due to their rapid reproductive cycles. About 7500 transcripts are differentially expressed in four sequential waves during the 72-h reproductive period in the fat body, a critical reproductive organ. The major regulators for these waves of gene expression are the two very important insect hormones, 20-hydroxyecdysone (20E) and Juvenile hormone (JH), their respective receptors Ecdysone Receptor (EcR) and Methoprene-Tolerant (Met), amino acids and the orphan nuclear receptor HR3. These key regulators are responsible for activation and repression of co-regulated gene sets, at different time points, within the 72-h reproductive period. Importantly, this study, apart from providing an insight into the regulatory complexity involved in the temporal coordination of gene expression, also reveals the previously unidentified roles of 20E/EcR, JH/Met and HR3 during the 72-h period post blood meal.

Heme Signaling Impacts Global Gene Expression, Immunity and Dengue Virus Infectivity in Aedes aegypti

Blood-feeding mosquitoes are exposed to high levels of heme, the product of hemoglobin degradation. Heme is a pro-oxidant that influences a variety of cellular processes. We performed a global analysis of heme-regulated Aedes aegypti (yellow fever mosquito) transcriptional changes to better understand influence on mosquito physiology at the molecular level. We observed an iron- and reactive oxygen species (ROS)-independent signaling induced by heme that comprised genes related to redox metabolism. By modulating the abundance of these transcripts, heme possibly acts as a danger signaling molecule. Furthermore, heme triggered critical changes in the expression of energy metabolism and immune response genes, altering the susceptibility towards bacteria and dengue virus. These findings seem to have implications on the adaptation of mosquitoes to hematophagy and consequently on their ability to transmit diseases. Altogether, these results may also contribute to the understanding of heme cell biology in eukaryotic cells.

Patterns of differential gene expression in adult rotation-resistant and wild-type western corn rootworm digestive tracts

The western corn rootworm (WCR,Diabrotica virgifera virgifera LeConte) is an important pest of corn. Annual crop rotation between corn and soybean disrupts the corn-dependent WCR life cycle and is widely adopted to manage this pest. This strategy selected for rotation-resistant (RR) WCR with reduced ovipositional fidelity to corn. Previous studies revealed that RR-WCR adults exhibit greater tolerance of soybean diets, different gut physiology, and host-microbe interactions compared to rotation-susceptible wild types (WT). To identify the genetic mechanisms underlying these phenotypic changes, a de novo assembly of the WCR adult gut transcriptome was constructed and used for RNA-sequencing analyses of RNA libraries from different WCR phenotypes fed with corn or soybean diets. Global gene expression profiles of WT- and RR-WCR were similar when feeding on corn diets, but different when feeding on soybean. Using network-based methods, we identified gene modules transcriptionally correlated with the RR phenotype. Gene ontology enrichment analyses indicated that the functions of these modules were related to metabolic processes, immune responses, biological adhesion, and other functions/processes that appear to correlate to documented traits in RR populations. These results suggest that gut transcriptomic divergence correlated with brief soybean feeding and other physiological traits may exist between RR- and WT-WCR adults.

Evaluation of reference genes for insect olfaction studies

BACKGROUND

Quantitative reverse transcription PCR (qRT-PCR) is a robust and accessible method to assay gene expression and to infer gene regulation. Being a chain of procedures, this technique is subject to systematic error due to biological and technical limitations mainly set by the starting material and downstream procedures. Thus, rigorous data normalization is critical to grant reliability and repeatability of gene expression quantification by qRT-PCR. A number of 'housekeeping genes', involved in basic cellular functions, have been commonly used as internal controls for this normalization process. However, these genes could themselves be regulated and must therefore be tested a priori.

METHODS

We evaluated eight potential reference genes for their stability as internal controls for RT-qPCR studies of olfactory gene expression in the antennae of Rhodnius prolixus, a Chagas disease vector. The set of genes included were: α-tubulin; β-actin; Glyceraldehyde-3-phosphate dehydrogenase; Eukaryotic initiation factor 1A; Glutathione-S-transferase; Serine protease; Succinate dehydrogenase; and Glucose-6-phosphate dehydrogenase. Five experimental conditions, including changes in age,developmental stage and feeding status were tested in both sexes.

RESULTS

We show that the evaluation of candidate reference genes is necessary for each combination of sex, tissue and physiological condition analyzed in order to avoid inconsistent results and conclusions. Although, Normfinder and geNorm software yielded different results between males and females, five genes (SDH, Tub, GAPDH, Act and G6PDH) appeared in the first positions in all rankings obtained. By using gene expression data of a single olfactory coreceptor gene as an example, we demonstrated the extent of changes expected using different internal standards.

CONCLUSIONS

This work underlines the need for a rigorous selection of internal standards to grant the reliability of normalization processes in qRT-PCR studies. Furthermore, we show that particular physiological or developmental conditions require independent evaluation of a diverse set of potential reference genes.

Global transcriptional dynamics of diapause induction in non-blood-fed and blood-fed Aedes albopictus

BACKGROUND

Aedes albopictus is a vector of increasing public health concern due to its rapid global range expansion and ability to transmit Dengue virus, Chikungunya virus and a wide range of additional arboviruses. Traditional vector control strategies have been largely ineffective against Ae. albopictus and novel approaches are urgently needed. Photoperiodic diapause is a crucial ecological adaptation in a wide range of temperate insects. Therefore, targeting the molecular regulation of photoperiodic diapause or diapause-associated physiological processes could provide the basis of novel approaches to vector control.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the global transcriptional profiles of diapause induction in Ae. albopictus by performing paired-end RNA-Seq of biologically replicated libraries. We sequenced RNA from whole bodies of adult females reared under diapause-inducing and non-diapause-inducing photoperiods either with or without a blood meal. We constructed a comprehensive transcriptome assembly that incorporated previous assemblies and represents over 14,000 annotated dipteran gene models. Mapping of sequence reads to the transcriptome identified differential expression of 2,251 genes in response to diapause-inducing short-day photoperiods. In non-blood-fed females, potential regulatory elements of diapause induction were transcriptionally up-regulated, including two of the canonical circadian clock genes, timeless and cryptochrome 1. In blood-fed females, genes in metabolic pathways related to energy production and offspring provisioning were differentially expressed under diapause-inducing conditions, including the oxidative phosphorylation pathway and lipid metabolism genes.

CONCLUSIONS/SIGNIFICANCE

This study is the first to utilize powerful RNA-Seq technologies to elucidate the transcriptional basis of diapause induction in any insect. We identified candidate genes and pathways regulating diapause induction, including a conserved set of genes that are differentially expressed as part of the diapause program in a diverse group of insects. These genes provide candidates whose diapause-associated function can be further interrogated using functional genomics approaches in Ae. albopictus and other insects.

The hub protein loquacious connects the microRNA and short interfering RNA pathways in mosquitoes

Aedes aegypti mosquitoes vector several arboviruses of global health significance, including dengue viruses and chikungunya virus. RNA interference (RNAi) plays an important role in antiviral immunity, gene regulation and protection from transposable elements. Double-stranded RNA binding proteins (dsRBPs) are important for efficient RNAi; in Drosophila functional specialization of the miRNA, endo-siRNA and exo-siRNA pathway is aided by the dsRBPs Loquacious (Loqs-PB, Loqs-PD) and R2D2, respectively. However, this functional specialization has not been investigated in other dipterans. We were unable to detect Loqs-PD in Ae. aegypti; analysis of other dipteran genomes demonstrated that this isoform is not conserved outside of Drosophila. Overexpression experiments and small RNA sequencing following depletion of each dsRBP revealed that R2D2 and Loqs-PA cooperate non-redundantly in siRNA production, and that these proteins exhibit an inhibitory effect on miRNA levels. Conversely, Loqs-PB alone interacted with mosquito dicer-1 and was essential for full miRNA production. Mosquito Loqs interacted with both argonaute 1 and 2 in a manner independent of its interactions with dicer. We conclude that the functional specialization of Loqs-PD in Drosophila is a recently derived trait, and that in other dipterans, including the medically important mosquitoes, Loqs-PA participates in both the miRNA and endo-siRNA based pathways.

Exome and transcriptome sequencing of Aedes aegypti identifies a locus that confers resistance to Brugia malayi and alters the immune response

Many mosquito species are naturally polymorphic for their abilities to transmit parasites, a feature which is of great interest for controlling vector-borne disease. Aedes aegypti, the primary vector of dengue and yellow fever and a laboratory model for studying lymphatic filariasis, is genetically variable for its capacity to harbor the filarial nematode Brugia malayi. The genome of Ae. aegypti is large and repetitive, making genome resequencing difficult and expensive. We designed exome captures to target protein-coding regions of the genome, and used association mapping in a wild Kenyan population to identify a single, dominant, sex-linked locus underlying resistance. This falls in a region of the genome where a resistance locus was previously mapped in a line established in 1936, suggesting that this polymorphism has been maintained in the wild for the at least 80 years. We then crossed resistant and susceptible mosquitoes to place both alleles of the gene into a common genetic background, and used RNA-seq to measure the effect of this locus on gene expression. We found evidence for Toll, IMD, and JAK-STAT pathway activity in response to early stages of B. malayi infection when the parasites are beginning to die in the resistant genotype. We also found that resistant mosquitoes express anti-microbial peptides at the time of parasite-killing, and that this expression is suppressed in susceptible mosquitoes. Together, we have found that a single resistance locus leads to a higher immune response in resistant mosquitoes, and we identify genes in this region that may be responsible for this trait.

VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases

VectorBase is a National Institute of Allergy and Infectious Diseases supported Bioinformatics Resource Center (BRC) for invertebrate vectors of human pathogens. Now in its 11th year, VectorBase currently hosts the genomes of 35 organisms including a number of non-vectors for comparative analysis. Hosted data range from genome assemblies with annotated gene features, transcript and protein expression data to population genetics including variation and insecticide-resistance phenotypes. Here we describe improvements to our resource and the set of tools available for interrogating and accessing BRC data including the integration of Web Apollo to facilitate community annotation and providing Galaxy to support user-based workflows. VectorBase also actively supports our community through hands-on workshops and online tutorials. All information and data are freely available from our website at https://www.vectorbase.org/.

Role of microRNAs in arbovirus/vector interactions

The role of microRNAs (miRNAs) as small non-coding RNAs in regulation of gene expression has been recognized. They appear to be involved in regulation of a wide range of cellular pathways that affect several biological processes such as development, the immune system, survival, metabolism and host-pathogen interactions. Arthropod-borne viruses impose great economic and health risks around the world. Recent advances in miRNA biology have shed some light on the role of these small RNAs in vector-virus interactions. In this review, I will reflect on our current knowledge on the role of miRNAs in arbovirus-vector interactions and the potential avenues for their utilization in limiting virus replication and/or transmission.

Molecular analysis of Aedes aegypti classical protein tyrosine phosphatases uncovers an ortholog of mammalian PTP-1B implicated in the control of egg production in mosquitoes

Background

Protein Tyrosine Phosphatases (PTPs) are enzymes that catalyze phosphotyrosine dephosphorylation and modulate cell differentiation, growth and metabolism. In mammals, PTPs play a key role in the modulation of canonical pathways involved in metabolism and immunity. PTP1B is the prototype member of classical PTPs and a major target for treating human diseases, such as cancer, obesity and diabetes. These signaling enzymes are, hence, targets of a wide array of inhibitors. Anautogenous mosquitoes rely on blood meals to lay eggs and are vectors of the most prevalent human diseases. Identifying the mosquito ortholog of PTP1B and determining its involvement in egg production is, therefore, important in the search for a novel and crucial target for vector control.

Methodology/Principal Findings

We conducted an analysis to identify the ortholog of mammalian PTP1B in the Aedes aegypti genome. We identified eight genes coding for classical PTPs. In silico structural and functional analyses of proteins coded by such genes revealed that four of these code for catalytically active enzymes. Among the four genes coding for active PTPs, AAEL001919 exhibits the greatest degree of homology with the mammalian PTP1B. Next, we evaluated the role of this enzyme in egg formation. Blood feeding largely affects AAEL001919 expression, especially in the fat body and ovaries. These tissues are critically involved in the synthesis and storage of vitellogenin, the major yolk protein. Including the classical PTP inhibitor sodium orthovanadate or the PTP substrate DiFMUP in the blood meal decreased vitellogenin synthesis and egg production. Similarly, silencing AAEL001919 using RNA interference (RNAi) assays resulted in 30% suppression of egg production.

Conclusions/Significance

The data reported herein implicate, for the first time, a gene that codes for a classical PTP in mosquito egg formation. These findings raise the possibility that this class of enzymes may be used as novel targets to block egg formation in mosquitoes.

De novo transcriptome sequencing and sequence analysis of the malaria vector Anopheles sinensis (Diptera: Culicidae)

Background

Anopheles sinensis is the major malaria vector in China and Southeast Asia. Vector control is one of the most effective measures to prevent malaria transmission. However, there is little transcriptome information available for the malaria vector. To better understand the biological basis of malaria transmission and to develop novel and effective means of vector control, there is a need to build a transcriptome dataset for functional genomics analysis by large-scale RNA sequencing (RNA-seq).

Methods

To provide a more comprehensive and complete transcriptome of An. sinensis, eggs, larvae, pupae, male adults and female adults RNA were pooled together for cDNA preparation, sequenced using the Illumina paired-end sequencing technology and assembled into unigenes. These unigenes were then analyzed in their genome mapping, functional annotation, homology, codon usage bias and simple sequence repeats (SSRs).

Results

Approximately 51.6 million clean reads were obtained, trimmed, and assembled into 38,504 unigenes with an average length of 571 bp, an N50 of 711 bp, and an average GC content 51.26%. Among them, 98.4% of unigenes could be mapped onto the reference genome, and 69% of unigenes could be annotated with known biological functions. Homology analysis identified certain numbers of An. sinensis unigenes that showed homology or being putative 1:1 orthologues with genomes of other Dipteran species. Codon usage bias was analyzed and 1,904 SSRs were detected, which will provide effective molecular markers for the population genetics of this species.

Conclusions

Our data and analysis provide the most comprehensive transcriptomic resource and characteristics currently available for An. sinensis, and will facilitate genetic, genomic studies, and further vector control of An. sinensis.

Transcriptome sequencing and developmental regulation of gene expression in Anopheles aquasalis

BACKGROUND

Anopheles aquasalis is a major malaria vector in coastal areas of South and Central America where it breeds preferentially in brackish water. This species is very susceptible to Plasmodium vivax and it has been already incriminated as responsible vector in malaria outbreaks. There has been no high-throughput investigation into the sequencing of An. aquasalis genes, transcripts and proteins despite its epidemiological relevance. Here we describe the sequencing, assembly and annotation of the An. aquasalis transcriptome.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 419 thousand cDNA sequence reads, encompassing 164 million nucleotides, were assembled in 7544 contigs of ≥ 2 sequences, and 1999 singletons. The majority of the An. aquasalis transcripts encode proteins with their closest counterparts in another neotropical malaria vector, An. darlingi. Several analyses in different protein databases were used to annotate and predict the putative functions of the deduced An. aquasalis proteins. Larval and adult-specific transcripts were represented by 121 and 424 contig sequences, respectively. Fifty-one transcripts were only detected in blood-fed females. The data also reveal a list of transcripts up- or down-regulated in adult females after a blood meal. Transcripts associated with immunity, signaling networks and blood feeding and digestion are discussed.

CONCLUSIONS/SIGNIFICANCE

This study represents the first large-scale effort to sequence the transcriptome of An. aquasalis. It provides valuable information that will facilitate studies on the biology of this species and may lead to novel strategies to reduce malaria transmission on the South American continent. The An. aquasalis transcriptome is accessible at http://exon.niaid.nih.gov/transcriptome/An_aquasalis/Anaquexcel.xlsx.

Evidence for intercellular communication in mosquito renal tubules: a putative role of gap junctions in coordinating and regulating the rapid diuretic effects of neuropeptides

Adult female mosquitoes require a blood meal from a vertebrate host to successfully reproduce. During a single blood feeding, a female may ingest more than the equivalent of her own body mass, resulting in an acute stress to osmotic and ionic homeostasis. In response to this stress, the renal (Malpighian) tubules mediate a rapid diuresis that commences as soon as blood is ingested. The diuresis is regulated by neuropeptides (e.g., kinins, calcitonin-like peptide) that act on receptors in the Malpighian tubule epithelium. Interestingly, the expression of these receptors is discontinuous throughout the epithelium, which raises the question as to how Malpighian tubules mount such a rapid and synchronized response to neuropeptide stimulation. Here we propose a hypothesis that gap junctions functionally couple the epithelial cells of Malpighian tubules, resulting in a coordinated physiological response to the binding of neuropeptides. We review recent, relevant literature on the electrophysiology, physiology, and molecular biology of mosquito Malpighian tubules that indicate the presence of gap junctions in the epithelium. We also provide new physiological and immunochemical data that are consistent with the proposed hypothesis.

Transcriptomic evidence for a dramatic functional transition of the malpighian tubules after a blood meal in the Asian tiger mosquito Aedes albopictus

Background

The consumption of a vertebrate blood meal by adult female mosquitoes is necessary for their reproduction, but it also presents significant physiological challenges to mosquito osmoregulation and metabolism. The renal (Malpighian) tubules of mosquitoes play critical roles in the initial processing of the blood meal by excreting excess water and salts that are ingested. However, it is unclear how the tubules contribute to the metabolism and excretion of wastes (e.g., heme, ammonia) produced during the digestion of blood.

Methodology/Principal Findings

Here we used RNA-Seq to examine global changes in transcript expression in the Malpighian tubules of the highly-invasive Asian tiger mosquito Aedes albopictus during the first 24 h after consuming a blood meal. We found progressive, global changes in the transcriptome of the Malpighian tubules isolated from mosquitoes at 3 h, 12 h, and 24 h after a blood meal. Notably, a DAVID functional cluster analysis of the differentially-expressed transcripts revealed 1) a down-regulation of transcripts associated with oxidative metabolism, active transport, and mRNA translation, and 2) an up-regulation of transcripts associated with antioxidants and detoxification, proteolytic activity, amino-acid metabolism, and cytoskeletal dynamics.

Conclusions/Significance

The results suggest that blood feeding elicits a functional transition of the epithelium from one specializing in active transepithelial fluid secretion (e.g., diuresis) to one specializing in detoxification and metabolic waste excretion. Our findings provide the first insights into the putative roles of mosquito Malpighian tubules in the chronic processing of blood meals.

The Asian tiger mosquito Aedes albopictus is a vector of several medically-important arboviruses and one of the most invasive mosquito species in the world. Existing control measures for mosquitoes are presently being challenged by the emergence of resistance to insecticides that target the nervous system. Thus, it is necessary to identify novel physiological targets to guide the development of new insecticides. We recently demonstrated that the ‘kidneys’ (Malpighian tubules) of mosquitoes offer a valuable, new physiological target for insecticides. However, our understanding of how this tissue contributes to the chronic metabolic processing of blood meals by mosquitoes is limited. Here we characterize the changes in transcript expression that occur in the Malpighian tubules of adult female A. albopictus with the goal of identifying key molecular pathways that may reveal valuable targets for insecticide development. We find dramatic changes in transcript accumulation in Malpighian tubules, which 1) provide new insights into the potential functional roles of Malpighian tubules after a blood meal, and 2) reveal new potential molecular pathways and targets to guide the development of new insecticides that would disrupt the renal functions of mosquitoes.

The reference transcriptome of the adult female biting midge (Culicoides sonorensis) and differential gene expression profiling during teneral, blood, and sucrose feeding conditions

Unlike other important vectors such as mosquitoes and sandflies, genetic and genomic tools for Culicoides biting midges are lacking, despite the fact that they vector a large number of arboviruses and other pathogens impacting humans and domestic animals world-wide. In North America, female Culicoides sonorensis midges are important vectors of bluetongue virus (BTV) and epizootic hemorrhagic disease virus (EHDV), orbiviruses that cause significant disease in livestock and wildlife. Libraries of tissue-specific transcripts expressed in response to feeding and oral orbivirus challenge in C. sonorensis have previously been reported, but extensive genome-wide expression profiling in the midge has not. Here, we successfully used deep sequencing technologies to construct the first adult female C. sonorensis reference transcriptome, and utilized genome-wide expression profiling to elucidate the genetic response to blood and sucrose feeding over time. The adult female midge unigene consists of 19,041 genes, of which less than 7% are differentially expressed during the course of a sucrose meal, while up to 52% of the genes respond significantly in blood-fed midges, indicating hematophagy induces complex physiological processes. Many genes that were differentially expressed during blood feeding were associated with digestion (e.g. proteases, lipases), hematophagy (e.g., salivary proteins), and vitellogenesis, revealing many major metabolic and biological factors underlying these critical processes. Additionally, key genes in the vitellogenesis pathway were identified, which provides the first glimpse into the molecular basis of anautogeny for C. sonorensis. This is the first extensive transcriptome for this genus, which will serve as a framework for future expression studies, RNAi, and provide a rich dataset contributing to the ultimate goal of informing a reference genome assembly and annotation. Moreover, this study will serve as a foundation for subsequent studies of genome-wide expression analyses during early orbivirus infection and dissecting the molecular mechanisms behind vector competence in midges.

Dual RNA-seq of parasite and host reveals gene expression dynamics during filarial worm-mosquito interactions

BACKGROUND

Parasite biology, by its very nature, cannot be understood without integrating it with that of the host, nor can the host response be adequately explained without considering the activity of the parasite. However, due to experimental limitations, molecular studies of parasite-host systems have been predominantly one-sided investigations focusing on either of the partners involved. Here, we conducted a dual RNA-seq time course analysis of filarial worm parasite and host mosquito to better understand the parasite processes underlying development in and interaction with the host tissue, from the establishment of infection to the development of infective-stage larva.

METHODOLOGY/PRINCIPAL FINDINGS

Using the Brugia malayi-Aedes aegypti system, we report parasite gene transcription dynamics, which exhibited a highly ordered developmental program consisting of a series of cyclical and state-transitioning temporal patterns. In addition, we contextualized these parasite data in relation to the concurrent dynamics of the host transcriptome. Comparative analyses using uninfected tissues and different host strains revealed the influence of parasite development on host gene transcription as well as the influence of the host environment on parasite gene transcription. We also critically evaluated the life-cycle transcriptome of B. malayi by comparing developmental stages in the mosquito relative to those in the mammalian host, providing insight into gene expression changes underpinning the mosquito-borne parasitic lifestyle of this heteroxenous parasite.

CONCLUSIONS/SIGNIFICANCE

The data presented herein provide the research community with information to design wet lab experiments and select candidates for future study to more fully dissect the whole set of molecular interactions of both organisms in this mosquito-filarial worm symbiotic relationship. Furthermore, characterization of the transcriptional program over the complete life cycle of the parasite, including stages within the mosquito, could help devise novel targets for control strategies.

Targeted genome editing in Aedes aegypti using TALENs

The Culicine mosquito, Aedes aegypti, is both a major vector of arthropod-borne viruses (arboviruses) and a genetic model organism for arbovirus transmission. TALE nucleases (TALENs), a group of artificial enzymes capable of generating site-specific DNA lesions, consist of a non-specific FokI endonuclease cleavage domain fused to an engineered DNA binding domain specific to a target site. While TALENs have become an important tool for targeted gene disruption in a variety of organisms, application to the mosquito genome is a new approach. We recently described the use of TALENs to perform heritable genetic disruptions in A. aegypti. Here, we provide detailed methods that will allow other research laboratories to capitalize on the potential of this technology for understanding mosquito gene function. We describe target site selection, transient embryo-based assays to rapidly assess TALEN activity, embryonic microinjection and downstream screening steps to identify target site mutations.

Gustatory receptor expression in the labella and tarsi of Aedes aegypti

The yellow-fever mosquito, Aedes aegypti, infects a growing number of people every year with dengue, yellow fever and chikungunya viruses. Contact chemoreception in mosquitoes influences a number of behaviors including host-selection, oviposition and feeding. While these behaviors are in many instances well documented, the molecular mechanisms mediating them are not well understood. Here we report the results of sequencing total messenger RNA in the labella and tarsi of both male and female Ae. aegypti to reveal Gustatory Receptor (GR) gene expression profiles in these major gustatory appendages. Gene expression levels in each tissue were verified by RT-qPCR. We discuss potential functions for the GRs revealed here by considering homologous GRs in other insects. Specific GRs provide molecular targets for modification of gustatory-mediated behaviors in this important disease vector.

Parasitization by Cotesia chilonis influences gene expression in fatbody and hemocytes of Chilo suppressalis

Background

During oviposition many parasitoid wasps inject various factors, such as polydnaviruses (PDVs), along with eggs that manipulate the physiology and development of their hosts. These manipulations are thought to benefit the parasites. However, the detailed mechanisms of insect host-parasitoid interactions are not fully understood at the molecular level. Based on recent findings that some parasitoids influence gene expression in their hosts, we posed the hypothesis that parasitization by a braconid wasp, Cotesia chilonis, influences the expression of genes responsible for development, metabolism and immune functions in the fatbody and hemocytes of its host, Chilo suppressalis.

Methodology/Principal Findings

We obtained 39,344,452 reads, which were assembled into 146,770 scaffolds, and 76,016 unigenes. Parasitization impacted gene expression in fatbody and hemocytes. Of these, 8096 fatbody or 5743 hemocyte unigenes were down-regulated, and 2572 fatbody or 1452 hemocyte unigenes were up-regulated. Gene ontology data showed that the majority of the differentially expressed genes are involved in enzyme-regulated activity, binding, transcription regulator activity and catalytic activity. qPCR results show that most anti-microbial peptide transcription levels were up-regulated after parasitization. Expression of bracovirus genes was detected in parasitized larvae with 19 unique sequences identified from six PDV gene families including ankyrin, CrV1 protein, cystatin, early-expressed (EP) proteins, lectin, and protein tyrosine phosphatase.

Conclusions

The current study supports our hypothesis that parasitization influences the expression of fatbody and hemocyte genes in the host, C. suppressalis. The general view is that manipulation of host metabolism and immunity benefits the development and emergence of the parasitoid offsprings. The accepted beneficial mechanisms include the direct impact of parasitoid-associated virulence factors such as venom and polydnavirus on host tissues (such as cell damage) and, more deeply, the ability of these factors to influence gene expression. We infer that insect parasitoids generally manipulate their environments, the internal milieu of their hosts.

The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector

Mosquitoes are vectors of a number of important human and animal diseases. The development of novel vector control strategies requires a thorough understanding of mosquito biology. To facilitate this, we used RNA-seq to identify novel genes and provide the first high-resolution view of the transcriptome throughout development and in response to blood feeding in a mosquito vector of human disease, Aedes aegypti, the primary vector for Dengue and yellow fever. We characterized mRNA expression at 34 distinct time points throughout Aedes development, including adult somatic and germline tissues, by using polyA+ RNA-seq. We identify a total of 14,238 novel new transcribed regions corresponding to 12,597 new loci, as well as many novel transcript isoforms of previously annotated genes. Altogether these results increase the annotated fraction of the transcribed genome into long polyA+ RNAs by more than twofold. We also identified a number of patterns of shared gene expression, as well as genes and/or exons expressed sex-specifically or sex-differentially. Expression profiles of small RNAs in ovaries, early embryos, testes, and adult male and female somatic tissues also were determined, resulting in the identification of 38 new Aedes-specific miRNAs, and ~291,000 small RNA new transcribed regions, many of which are likely to be endogenous small-interfering RNAs and Piwi-interacting RNAs. Genes of potential interest for transgene-based vector control strategies also are highlighted. Our data have been incorporated into a user-friendly genome browser located at www.Aedes.caltech.edu, with relevant links to Vectorbase (www.vectorbase.org)

Transcriptome characterization by RNA-Seq reveals the involvement of the complement components in noise-traumatized rat cochleae

Acoustic trauma, a leading cause of sensorineural hearing loss in adults, induces a complex degenerative process in the cochlea. Although previous investigations have identified multiple stress pathways, a comprehensive analysis of cochlear responses to acoustic injury is still lacking. In the current study, we used the next-generation RNA-sequencing (RNA-Seq) technique to sequence the whole transcriptome of the normal and noise-traumatized cochlear sensory epithelia (CSE). CSE tissues were collected from rat inner ears 1d after the rats were exposed to a 120-dB (sound pressure level) noise for 2 h. The RNA-Seq generated over 176 million sequence reads for the normal CSE and over 164 million reads for the noise-traumatized CSE. Alignment of these sequences with the rat Rn4 genome revealed the expression of over 17,000 gene transcripts in the CSE, over 2000 of which were exclusively expressed in either the normal or noise-traumatized CSE. Seventy-eight gene transcripts were differentially expressed (70 upregulated and 8 downregulated) after acoustic trauma. Many of the differentially expressed genes are related to the innate immune system. Further expression analyses using quantitative real time PCR confirmed the constitutive expression of multiple complement genes in the normal organ of Corti and the changes in the expression levels of the complement factor I (Cfi) and complement component 1, s subcomponent (C1s) after acoustic trauma. Moreover, protein expression analysis revealed strong expression of Cfi and C1s proteins in the organ of Corti. Importantly, these proteins exhibited expression changes following acoustic trauma. Collectively, the results of the current investigation suggest the involvement of the complement components in cochlear responses to acoustic trauma.

Detailed transcriptome description of the neglected cestode Taenia multiceps

BACKGROUND

The larval stage of Taenia multiceps, a global cestode, encysts in the central nervous system (CNS) of sheep and other livestock. This frequently leads to their death and huge socioeconomic losses, especially in developing countries. This parasite can also cause zoonotic infections in humans, but has been largely neglected due to a lack of diagnostic techniques and studies. Recent developments in next-generation sequencing provide an opportunity to explore the transcriptome of T. multiceps.

METHODOLOGY/PRINCIPAL FINDINGS

We obtained a total of 31,282 unigenes (mean length 920 bp) using Illumina paired-end sequencing technology and a new Trinity de novo assembler without a referenced genome. Individual transcription molecules were determined by sequence-based annotations and/or domain-based annotations against public databases (Nr, UniprotKB/Swiss-Prot, COG, KEGG, UniProtKB/TrEMBL, InterPro and Pfam). We identified 26,110 (83.47%) unigenes and inferred 20,896 (66.8%) coding sequences (CDS). Further comparative transcripts analysis with other cestodes (Taenia pisiformis, Taenia solium, Echincoccus granulosus and Echincoccus multilocularis) and intestinal parasites (Trichinella spiralis, Ancylostoma caninum and Ascaris suum) showed that 5,100 common genes were shared among three Taenia tapeworms, 261 conserved genes were detected among five Taeniidae cestodes, and 109 common genes were found in four zoonotic intestinal parasites. Some of the common genes were genes required for parasite survival, involved in parasite-host interactions. In addition, we amplified two full-length CDS of unigenes from the common genes using RT-PCR.

CONCLUSIONS/SIGNIFICANCE

This study provides an extensive transcriptome of the adult stage of T. multiceps, and demonstrates that comparative transcriptomic investigations deserve to be further studied. This transcriptome dataset forms a substantial public information platform to achieve a fundamental understanding of the biology of T. multiceps, and helps in the identification of drug targets and parasite-host interaction studies.

Sindbis virus induces the production of a novel class of endogenous siRNAs in Aedes aegypti mosquitoes

Small RNA regulatory pathways are used to control the activity of transposons, regulate gene expression and resist infecting viruses. We examined the biogenesis of mRNA-derived endogenous short-interfering RNAs (endo-siRNAs) in the disease vector mosquito Aedes aegypti. Under standard conditions, mRNA-derived endo-siRNAs were produced from the bidirectional transcription of tail-tail overlapping gene pairs. Upon infection with the alphavirus, Sindbis virus (SINV), another class of mRNA-derived endo-siRNAs was observed. Genes producing SINV-induced endo-siRNAs were not enriched for overlapping partners or nearby genes, but were enriched for transcripts with long 3' untranslated regions. Endo-siRNAs from this class derived uniformly from the entire length of the target transcript, and were found to regulate the transcript levels of the genes from which they were derived. Strand-specific quantitative PCR experiments demonstrated that antisense strands of targeted mRNA genes were produced to exonic, but not intronic regions. Finally, small RNAs mapped to both sense and antisense strands of exon-exon junctions, suggesting double-stranded RNA precursors to SINV-induced endo-siRNAs may be synthesized from mature mRNA templates. These results suggest additional complexity in small RNA pathways and gene regulation in the presence of an infecting virus in disease vector mosquitoes.

RNA-Seq and molecular docking reveal multi-level pesticide resistance in the bed bug

BACKGROUND

Bed bugs (Cimex lectularius) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of C. lectularius has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance.

RESULTS

We performed a next-generation RNA sequencing (RNA-Seq) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from de novo assembled mRNA-Seq tags (30,404 ESTs) and our previous 454 pyrosequenced database (21,088 ESTs). The two-way GLMseq analysis revealed ~15,000 highly significant differentially expressed ESTs between the PR and PS strains. Among the top 5,000 differentially expressed ESTs, 109 putative defense genes (cuticular proteins, cytochrome P450s, antioxidant genes, ABC transporters, glutathione S-transferases, carboxylesterases and acetyl cholinesterase) involved in penetration resistance and metabolic resistance were identified. Tissue and development-specific expression of P450 CYP3 clan members showed high mRNA levels in the cuticle, Malpighian tubules, and midgut; and in early instar nymphs, respectively. Lastly, molecular modeling and docking of a candidate cytochrome P450 (CYP397A1V2) revealed the flexibility of the deduced protein to metabolize a broad range of insecticide substrates including DDT, deltamethrin, permethrin, and imidacloprid.

CONCLUSIONS

We developed significant molecular resources for C. lectularius putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-Seq profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in C. lectularius. Future research that is targeted towards RNA interference (RNAi) on the identified metabolic targets such as cytochrome P450s and cuticular proteins could lay the foundation for a better understanding of the genetic basis of insecticide resistance in C. lectularius.

Strain Variation in the Transcriptome of the Dengue Fever Vector, Aedes aegypti

Studies of transcriptome dynamics provide a basis for understanding functional elements of the genome and the complexity of gene regulation. The dengue vector mosquito, Aedes aegypti, exhibits great adaptability to diverse ecological conditions, is phenotypically polymorphic, and shows variation in vectorial capacity to arboviruses. Previous genome sequencing showed richness in repetitive DNA and transposable elements that can contribute to genome plasticity. Population genetic studies revealed a varying degree of worldwide genetic polymorphism. However, the extent of functional genetic polymorphism across strains is unknown. The transcriptomes of three Ae. aegypti strains, Chetumal (CTM), Rexville D-Puerto Rico (Rex-D) and Liverpool (LVP), were compared. CTM is more susceptible than Rex- D to infection by dengue virus serotype 2. A total of 4188 transcripts exhibit either no or small variation (<2-fold) among sugar-fed samples of the three strains and between sugar- and blood-fed samples within each strain, corresponding most likely to genes encoding products necessary for vital functions. Transcripts enriched in blood-fed mosquitoes encode proteins associated with catalytic activities, molecular transport, metabolism of lipids, carbohydrates and amino acids, and functions related to blood digestion and the progression of the gonotropic cycle. Significant qualitative and quantitative differences were found in individual transcripts among strains including differential representation of paralogous gene products. The majority of immunity-associated transcripts decreased in accumulation after a bloodmeal and the results are discussed in relation to the different susceptibility of CTM and Rex-D mosquitoes to DENV2 infection.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

Targeted gene expression in the transgenic Aedes aegypti using the binary Gal4-UAS system

In this study, we report the establishment of the binary Gal4/UAS system for the yellow fever mosquito Aedes aegypti. We utilized the 1.8-kb 5' upstream region of the vitellogenin gene (Vg) to genetically engineer mosquito lines with the Vg-Gal4 activator and established UAS-EGFP responder transgenic mosquito lines to evaluate the binary Gal4/UAS system. The results show that the Vg-Gal4 driver leads to a high level of tissue-, stage- and sex-specific expression of the EGFP reporter in the fat body of Vg-Gal4/UAS-EGFP hybrids after blood-meal activation. In addition, the applicability of this system to study hormonal regulation of gene expression was demonstrated in in vitro organ culture experiments in which the EGFP reporter was highly activated in isolated fat bodies of previtellogenic Vg-Gal4/UAS-EGFP females incubated in the presence of 20-hydroxyecdysone (20E). Hence, this study has opened the door for further refinement of genetic tools in mosquitoes.