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

A comparative analysis of corpora allata-corpora cardiaca microRNA repertoires revealed significant changes during mosquito metamorphosis

The corpora allata (CA) are a pair of endocrine glands with neural connections to the brain and close association with another neuroendocrine organ, the corpora cardiaca (CC). The CA from adult female Aedes aegypti mosquitoes synthesize fluctuating levels of juvenile hormone (JH), which have been linked to the ovarian development and are influenced by nutritional signals. In this study, we investigated the potential involvement of microRNAs (miRNAs), a type of small non-coding RNAs, in the regulation of gene expression in CA-CC complexes during mosquito reproductive development, at stages with distinct JH biosynthesis patterns. We analyzed the miRNA repertoires expressed in the CA-CC of pupae, sugar-fed and blood-fed female Ae. aegypti. In total, 156 mature miRNAs were detected in the CA-CC, with 84 displaying significant differences in expression among the three CA-CC developmental stages. There were more miRNAs that were expressed in pupae, and decreased or were absent after adult emergence, when compared with changes between CA-CC of sugar and blood-fed females. Analysis of the genes identified as potential targets for the CA-CC miRNA repertoires classified them into the broad categories of metabolism, information storage and processing, and cellular processes and signaling; with genes involved in cellular processes and signaling representing the largest portion. Among them, the signal-transduction mechanisms and intracellular trafficking, secretion and vesicular transport contained almost 55% of the genes' targets. A substantial number of miRNAs were differentially abundant in the libraries of the three developmental stages, and those changes were much more notable when pupae and adult stages were compared. We detected putative binding sites for some of the most abundant miRNAs on genes encoding JH biosynthetic enzymes and CC neuropeptides. These studies should help us to gain a better understanding of the regulation of CA-CC activity mediated by miRNAs during major developmental stages in mosquitoes.

piRNA Profiling of Dengue Virus Type 2-Infected Asian Tiger Mosquito and Midgut Tissues

The Asian tiger mosquito, Aedes albopictus, is a competent vector for the majority of arboviruses. The mosquito innate immune response is a primary determinant for arthropod-borne virus transmission, and the midgut is the first barrier to pathogen transmission. Mosquito antiviral immunity is primarily mediated by the small interfering RNA pathway. However, the roles that the P-element induced wimpy testis (PIWI)-interacting RNA (piRNA) pathway play in antiviral immunity in Ae. albopictus and its midgut still need further exploration. This study aimed to explore the profiles of both viral-derived and host-originated piRNAs in the whole body and midgut infected with Dengue virus 2 (DENV-2) in Ae. albopictus, and to elucidate gene expression profile differences of the PIWI protein family between adult females and their midguts. A deep sequencing-based method was used to identify and analyze small non-coding RNAs, especially the piRNA profiles in DENV-2-infected Ae. albopictus and its midgut. The top-ranked, differentially-expressed piRNAs were further validated using Stem-loop qRT-PCR. Bioinformatics analyses and reverse-transcription PCR (RT-PCR) methods were used to detect PIWI protein family members, and their expression profiles. DENV-2 derived piRNAs (vpiRNA, 24–30 nts) were observed in both infected Ae. albopictus and its midgut; however, only vpiRNA in the whole-body library had a weak preference for adenine at position 10 (10A) in the sense molecules as a feature of secondary piRNA. These vpiRNAs were not equally distributed, instead they were derived from a few specific regions of the genome, especially several hot spots, and displayed an obvious positive strand bias. We refer to the differentially expressed host piRNAs after DENV infection as virus-induced host endogenous piRNAs (vepiRNAs). However, we found that vepiRNAs were abundant in mosquito whole-body tissue, but deficient in the midgut. A total of eleven PIWI family genes were identified in Ae. albopictus; however, only AalPiwi5–7 and AalAgo3(1–2) were readily detected in the midgut. The characteristics of piRNAs in DENV-2-infected Ae. albopictus adult females were similar to those previously described for flavivirus infections but were not observed in the midgut. The reduced levels of vepiRNAs and incomplete expression of PIWI pathway genes in midgut samples from DENV-2-infected Ae. albopictus suggests that viral regulation of host piRNAs may not be an important factor in the midgut.

Characterization and potential role of microRNA in the Chinese dominant malaria mosquito Anopheles sinensis (Diptera: Culicidae) throughout four different life stages

BACKGROUND

microRNAs (miRNAs) are one kind of small non-coding RNAs widely distributed in insects. Many studies have shown that miRNAs play critical roles in development, differentiation, apoptosis, and innate immunity. However, there are a few reports describing miRNAs in Anopheles sinensis, the most common, and one of the dominant malaria mosquito in China. Here, we investigated the global miRNA expression profile across four different developmental stages including embryo, larval, pupal, and adult stages using Illumina Hiseq 2500 sequencing.

RESULTS

In total, 164 miRNAs were obtained out of 107.46 million raw sequencing reads. 99 of them identified as known miRNAs, and the remaining 65 miRNAs were considered as novel. By analyzing the read counts of miRNAs in all developmental stages, 95 miRNAs showed stage-specific expression (q < 0.01 and |log2 (fold change)| > 1) in consecutive stages, indicating that these miRNAs may be involved in critical physiological activity during development. Sixteen miRNAs were identified to be commonly dysregulated throughout four developmental stages. Many miRNAs showed stage-specific expression, such as asi-miR-2943 was exclusively expressed in the embryo stage, and asi-miR-1891 could not be detected in larval stage. The expression of six selected differentially expressed miRNAs identified by qRT-PCR were consistent with our sequencing results. Furthermore, 5296 and 1902 target genes were identified for the dysregulated 68 known and 27 novel miRNAs respectively by combining miRanda and RNAhybrid prediction. GO annotation and KEGG pathway analysis for the predicted genes of dysregulated miRNAs revealed that they might be involved in a broad range of biological processes related with the development, such as membrane, organic substance transport and several key pathways including protein processing in endoplasmic reticulum, propanoate metabolism and folate biosynthesis. Thirty-two key miRNAs were identified by microRNA-gene network analysis.

CONCLUSION

The present study represents the first global characterization of An. sinensis miRNAs in its four developmental stages. The presence and differential expression of An. sinensis miRNAs imply that such miRNAs may play critical roles in An. sinensis life cycle. A better understanding of the functions of these miRNAs will have great implication for the effective control of vector population and therefore interrupting malaria transmission.

Modulation of mitochondrial metabolism as a biochemical trait in blood feeding organisms: the redox vampire hypothesis redux

Hematophagous organisms undergo remarkable metabolic changes during the blood digestion process, increasing fermentative glucose metabolism, and reducing respiratory rates, both consequence of functional mitochondrial remodeling. Here, we review the pathways involved in energy metabolism and mitochondrial functionality in a comparative framework across different hematophagous species, and consider how these processes regulate redox homeostasis during blood digestion. The trend across distinct species indicate that a switch in energy metabolism might represent an important defensive mechanism to avoid the potential harmful interaction of oxidants generated from aerobic energy metabolism with products derived from blood digestion. Indeed, in insect vectors, blood feeding transiently reduces respiratory rates and oxidant production, irrespective of tissue and insect model. On the other hand, a different scenario is observed in several unrelated parasite species when exposed to blood digestion products, as respiratory rates reduce and mitochondrial oxidant production increase. The emerging picture indicates that re-wiring of energy metabolism, through reduced mitochondrial function, culminates in improved tolerance to redox insults and seems to represent a key step for hematophagous organisms to cope with the overwhelming and potentially toxic blood meal.

A mosquito sperm's journey from male ejaculate to egg: Mechanisms, molecules, and methods for exploration

The fate of mosquito sperm in the female reproductive tract has been addressed sporadically and incompletely, resulting in significant gaps in our understanding of sperm-female interactions that ultimately lead to fertilization. As with other Diptera, mosquito sperm have a complex journey to their ultimate destination, the egg. After copulation, sperm spend a short time at the site of insemination where they are hyperactivated and quickly congregate near the entrance of the spermathecal ducts. Within minutes, they travel up the narrow ducts to the spermathecae, likely through the combined efforts of female transport and sperm locomotion. The female nourishes sperm and maintains them in these permanent storage organs for her entire life. When she is ready, the female coordinates the release of sperm with ovulation, and the descending egg is fertilized. Although this process has been well studied via microscopy, many questions remain regarding the molecular processes that coordinate sperm motility, movement through the reproductive tract, maintenance, and usage. In this review, we describe the current understanding of a mosquito sperm's journey to the egg, highlighting gaps in our knowledge of mosquito reproductive biology. Where insufficient information is available in mosquitoes, we describe analogous processes in other organisms, such as Drosophila melanogaster, as a basis for comparison, and we suggest future areas of research that will illuminate how sperm successfully traverse the female reproductive tract. Such studies may yield molecular targets that could be manipulated to control populations of vector species. Mol. Reprod. Dev. 83: 897-911, 2016 © 2016 Wiley Periodicals, Inc.

Multifaceted functional implications of an endogenously expressed tRNA fragment in the vector mosquito Aedes aegypti

The mosquito Aedes aegypti is the primary vector of human arboviral diseases caused by dengue, chikungunya and Zika viruses. Many studies have shown the potential roles of small RNA molecules such as microRNA, small interfering RNA and PIWI-interacting RNA in vector mosquitoes. The function of tRNA fragments (tRF), the newly discovered class of small RNAs, in mosquitoes is not known. In this study, we show that specific tRFs are expressed in significantly differential manner between males and females of Ae. aegypti strains. Specific tRFs also show differential response during developmental transition from larvae to adults, as well as after blood feeding of adult females. The expression pattern of tRFs upon blood feeding varied depending upon if the blood contained dengue virus, and also if the females were treated with antibiotic prior to feeding to cleanse of the gut bacteria. Our findings show that a single tRF derived from the precursor sequences of a tRNA-Gly was differentially expressed between males and females, developmental transitions and also upon blood feeding by females of two laboratory strains that vary in midgut susceptibility to dengue virus infection. The multifaceted functional implications of this specific tRF suggest that biogenesis of small regulatory molecules from a tRNA can have wide ranging effects on key aspects of Ae. aegypti vector biology.

The mosquito Aedes aegypti is a major vector of arboviral diseases in subtropics and tropics. The confounding effects of immature development and adult microbiome on the ability of Ae. aegypti to transmit diseases (vector competence) have gained renewed attention in the recent years. However, the molecular nature of these links/ effects remains unknown. This is major gap in knowledge regarding how vector competence is regulated at molecular level, and how that regulation may be variable among different strains of this mosquito. In this study, we investigated expression of newly discovered class of small RNAs, called tRNA fragments (tRF) in Ae. aegypti strains. Based on small RNA sequencing and bioinformatics analyses, we show that tRFs are expressed in Ae. aegypti, and they are associated with significant changes in expression between males and females, during development stages, and post blood feeding responses. A single tRF showed association with sex-biased expression, developmental regulation and in response to blood meals between Moyo-S and Moyo-R strains that differ in midgut susceptibility to dengue virus. The findings of this study are expected to guide future research efforts directed toward examining detailed regulatory mechanisms of tRFs in vector competence of Ae. aegypti to disease transmission.

Germline Cas9 expression yields highly efficient genome engineering in a major worldwide disease vector, Aedes aegypti

The development of CRISPR/Cas9 technologies has dramatically increased the accessibility and efficiency of genome editing in many organisms. In general, in vivo germline expression of Cas9 results in substantially higher activity than embryonic injection. However, no transgenic lines expressing Cas9 have been developed for the major mosquito disease vector Aedes aegypti Here, we describe the generation of multiple stable, transgenic Ae. aegypti strains expressing Cas9 in the germline, resulting in dramatic improvements in both the consistency and efficiency of genome modifications using CRISPR. Using these strains, we disrupted numerous genes important for normal morphological development, and even generated triple mutants from a single injection. We have also managed to increase the rates of homology-directed repair by more than an order of magnitude. Given the exceptional mutagenic efficiency and specificity of the Cas9 strains we engineered, they can be used for high-throughput reverse genetic screens to help functionally annotate the Ae. aegypti genome. Additionally, these strains represent a step toward the development of novel population control technologies targeting Ae. aegypti that rely on Cas9-based gene drives.

Assembling large genomes: analysis of the stick insect (Clitarchus hookeri) genome reveals a high repeat content and sex-biased genes associated with reproduction

BACKGROUND

Stick insects (Phasmatodea) have a high incidence of parthenogenesis and other alternative reproductive strategies, yet the genetic basis of reproduction is poorly understood. Phasmatodea includes nearly 3000 species, yet only the genome of Timema cristinae has been published to date. Clitarchus hookeri is a geographical parthenogenetic stick insect distributed across New Zealand. Sexual reproduction dominates in northern habitats but is replaced by parthenogenesis in the south. Here, we present a de novo genome assembly of a female C. hookeri and use it to detect candidate genes associated with gamete production and development in females and males. We also explore the factors underlying large genome size in stick insects.

RESULTS

The C. hookeri genome assembly was 4.2 Gb, similar to the flow cytometry estimate, making it the second largest insect genome sequenced and assembled to date. Like the large genome of Locusta migratoria, the genome of C. hookeri is also highly repetitive and the predicted gene models are much longer than those from most other sequenced insect genomes, largely due to longer introns. Miniature inverted repeat transposable elements (MITEs), absent in the much smaller T. cristinae genome, is the most abundant repeat type in the C. hookeri genome assembly. Mapping RNA-Seq reads from female and male gonadal transcriptomes onto the genome assembly resulted in the identification of 39,940 gene loci, 15.8% and 37.6% of which showed female-biased and male-biased expression, respectively. The genes that were over-expressed in females were mostly associated with molecular transportation, developmental process, oocyte growth and reproductive process; whereas, the male-biased genes were enriched in rhythmic process, molecular transducer activity and synapse. Several genes involved in the juvenile hormone synthesis pathway were also identified.

CONCLUSIONS

The evolution of large insect genomes such as L. migratoria and C. hookeri genomes is most likely due to the accumulation of repetitive regions and intron elongation. MITEs contributed significantly to the growth of C. hookeri genome size yet are surprisingly absent from the T. cristinae genome. Sex-biased genes identified from gonadal tissues, including genes involved in juvenile hormone synthesis, provide interesting candidates for the further study of flexible reproduction in stick insects.

Histone-derived piRNA biogenesis depends on the ping-pong partners Piwi5 and Ago3 in Aedes aegypti

The piRNA pathway is of key importance in controlling transposable elements in most animal species. In the vector mosquito Aedes aegypti, the presence of eight PIWI proteins and the accumulation of viral piRNAs upon arbovirus infection suggest additional functions of the piRNA pathway beyond genome defense. To better understand the regulatory potential of this pathway, we analyzed in detail host-derived piRNAs in A. aegypti Aag2 cells. We show that a large repertoire of protein-coding genes and non-retroviral integrated RNA virus elements are processed into genic piRNAs by different combinations of PIWI proteins. Among these, we identify a class of genes that produces piRNAs from coding sequences in an Ago3- and Piwi5-dependent fashion. We demonstrate that the replication-dependent histone gene family is a genic source of ping-pong dependent piRNAs and that histone-derived piRNAs are dynamically expressed throughout the cell cycle, suggesting a role for the piRNA pathway in the regulation of histone gene expression. Moreover, our results establish the Aag2 cell line as an accessible experimental model to study gene-derived piRNAs.

Gene Drive for Mosquito Control: Where Did It Come from and Where Are We Headed?

Mosquito-borne pathogens place an enormous burden on human health. The existing toolkit is insufficient to support ongoing vector-control efforts towards meeting disease elimination and eradication goals. The perspective that genetic approaches can potentially add a significant set of tools toward mosquito control is not new, but the recent improvements in site-specific gene editing with CRISPR/Cas9 systems have enhanced our ability to both study mosquito biology using reverse genetics and produce genetics-based tools. Cas9-mediated gene-editing is an efficient and adaptable platform for gene drive strategies, which have advantages over innundative release strategies for introgressing desirable suppression and pathogen-blocking genotypes into wild mosquito populations; until recently, an effective gene drive has been largely out of reach. Many considerations will inform the effective use of new genetic tools, including gene drives. Here we review the lengthy history of genetic advances in mosquito biology and discuss both the impact of efficient site-specific gene editing on vector biology and the resulting potential to deploy new genetic tools for the abatement of mosquito-borne disease.

MicroRNA-275 targets sarco/endoplasmic reticulum Ca2+ adenosine triphosphatase (SERCA) to control key functions in the mosquito gut

The yellow fever mosquito Aedes aegypti is the major vector of arboviruses, causing numerous devastating human diseases, such as dengue and yellow fevers, Chikungunya and Zika. Female mosquitoes need vertebrate blood for egg development, and repeated cycles of blood feeding are tightly linked to pathogen transmission. The mosquito’s posterior midgut (gut) is involved in blood digestion and also serves as an entry point for pathogens. Thus, the mosquito gut is an important tissue to investigate. The miRNA aae-miR-275 (miR-275) has been shown to be required for normal blood digestion in the female mosquito; however, the mechanism of its action has remained unknown. Here, we demonstrate that miR-275 directly targets and positively regulates sarco/endoplasmic reticulum Ca²⁺adenosine triphosphatase, which is implicated in active transport of Ca²⁺ from the cytosol to the sarco/endoplasmic reticulum. We utilized a combination of the gut-specific yeast transcription activator protein Gal4/upstream activating sequence (Gal4/UAS) system and miRNA Tough Decoy technology to deplete the endogenous level of miR-275 in guts of transgenic mosquitoes. This gut-specific reduction of miR-275 post blood meal decreased SERCA mRNA and protein levels of the digestive enzyme late trypsin. It also resulted in a significant reduction of gut microbiota. Moreover, the decrease of miR-275 and SERCA correlated with defects in the Notch signaling pathway and assembly of the gut actin cytoskeleton. The adverse phenotypes caused by miR-275 silencing were rescued by injections of miR-275 mimic. Thus, we have discovered that miR-275 directly targets SERCA, and the maintenance of its level is critical for multiple gut functions in mosquitoes.

Female mosquitoes transmit numerous devastating human diseases. The mosquito gut, in addition to its primary function as a site of blood digestion, represents the entry point for pathogen colonization in mosquito vectors. The conserved microRNA, miR-275, was shown to be required for blood digestion and egg development. In this study, we investigated the target of miR-275 contributing to the regulation of mosquito gut functions. We achieved spatiotemporal suppression of miR-275 using a transgenic Tough Decoy RNA approach in the A. aegypti female mosquito gut. Furthermore, we have uncovered that miR-275 targets sarco/endoplasmic reticulum Ca²⁺- adenosine triphosphatase (SERCA), affecting numerous gut functions including blood digestion, production of digestive proteases, and assembly of the gut actin cytoskeleton. SERCA is essential for maintenance of Ca²⁺ homeostasis, and its disturbance, in humans, leads to cardiac hypertrophy, heart failure and cancers. Therefore, the finding that the miRNA miR-275 targets SERCA not only contributes to the knowledge of mosquito gut regulation but also significantly adds to the general understanding of mechanisms governing this critical molecule.

Rapid Evolution of Ovarian-Biased Genes in the Yellow Fever Mosquito (Aedes aegypti)

Males and females exhibit highly dimorphic phenotypes, particularly in their gonads, which is believed to be driven largely by differential gene expression. Typically, the protein sequences of genes upregulated in males, or male-biased genes, evolve rapidly as compared to female-biased and unbiased genes. To date, the specific study of gonad-biased genes remains uncommon in metazoans. Here, we identified and studied a total of 2927, 2013, and 4449 coding sequences (CDS) with ovary-biased, testis-biased, and unbiased expression, respectively, in the yellow fever mosquito Aedes aegypti The results showed that ovary-biased and unbiased CDS had higher nonsynonymous to synonymous substitution rates (dN/dS) and lower optimal codon usage (those codons that promote efficient translation) than testis-biased genes. Further, we observed higher dN/dS in ovary-biased genes than in testis-biased genes, even for genes coexpressed in nonsexual (embryo) tissues. Ovary-specific genes evolved exceptionally fast, as compared to testis- or embryo-specific genes, and exhibited higher frequency of positive selection. Genes with ovary expression were preferentially involved in olfactory binding and reception. We hypothesize that at least two potential mechanisms could explain rapid evolution of ovary-biased genes in this mosquito: (1) the evolutionary rate of ovary-biased genes may be accelerated by sexual selection (including female-female competition or male-mate choice) affecting olfactory genes during female swarming by males, and/or by adaptive evolution of olfactory signaling within the female reproductive system (e.g., sperm-ovary signaling); and/or (2) testis-biased genes may exhibit decelerated evolutionary rates due to the formation of mating plugs in the female after copulation, which limits male-male sperm competition.

Evaluation of reference genes at different developmental stages for quantitative real-time PCR in Aedes aegypti

The mosquito Aedes aegypti (Ae. aegypti) is the most notorious vector of illness-causing viruses such as Dengue, Chikugunya, and Zika. Although numerous genetic expression studies utilizing quantitative real-time PCR (qPCR) have been conducted with regards to Ae. aegypti, a panel of genes to be used suitably as references for the purpose of expression-level normalization within this epidemiologically important insect is presently lacking. Here, the usability of seven widely-utilized reference genes i.e. actin (ACT), eukaryotic elongation factor 1 alpha (eEF1α), alpha tubulin (α-tubulin), ribosomal proteins L8, L32 and S17 (RPL8, RPL32 and RPS17), and glyceraldeyde 3-phosphate dehydrogenase (GAPDH) were investigated. Expression patterns of the reference genes were observed in sixteen pre-determined developmental stages and in cell culture. Gene stability was inferred from qPCR data through three freely available algorithms i.e. BestKeeper, geNorm, and NormFinder. The consensus rankings generated from stability values provided by these programs suggest a combination of at least two genes for normalization. ACT and RPS17 are the most dependably expressed reference genes and therefore, we propose an ACT/RPS17 combination for normalization in all Ae. aegypti derived samples. GAPDH performed least desirably, and is thus not a recommended reference gene. This study emphasizes the importance of validating reference genes in Ae. aegypti for qPCR based research.

Characterization of Carbonic Anhydrase 9 in the Alimentary Canal of Aedes aegypti and Its Relationship to Homologous Mosquito Carbonic Anhydrases

In the mosquito midgut, luminal pH regulation and cellular ion transport processes are important for the digestion of food and maintenance of cellular homeostasis. pH regulation in the mosquito gut is affected by the vectorial movement of the principal ions including bicarbonate/carbonate and protons. As in all metazoans, mosquitoes employ the product of aerobic metabolism carbon dioxide in its bicarbonate/carbonate form as one of the major buffers of cellular and extracellular pH. The conversion of metabolic carbon dioxide to bicarbonate/carbonate is accomplished by a family of enzymes encoded by the carbonic anhydrase gene family. This study characterizes Aedes aegypti carbonic anhydrases using bioinformatic, molecular, and immunohistochemical methods. Our analyses show that there are fourteen Aedes aegypti carbonic anhydrase genes, two of which are expressed as splice variants. The carbonic anhydrases were classified as either integral membrane, peripheral membrane, mitochondrial, secreted, or soluble cytoplasmic proteins. Using polymerase chain reaction and Western blotting, one of the carbonic anhydrases, Aedes aegypti carbonic anhydrase 9, was analyzed and found in each life stage, male/female pupae, male/female adults, and in the female posterior midgut. Next, carbonic anhydrase 9 was analyzed in larvae and adults using confocal microscopy and was detected in the midgut regions. According to our analyses, carbonic anhydrase 9 is a soluble cytoplasmic enzyme found in the alimentary canal of larvae and adults and is expressed throughout the life cycle of the mosquito. Based on previous physiological analyses of adults and larvae, it appears AeCA9 is one of the major carbonic anhydrases involved in producing bicarbonate/carbonate which is involved in pH regulation and ion transport processes in the alimentary canal. Detailed understanding of the molecular bases of ion homeostasis in mosquitoes will provide targets for novel mosquito control strategies into the new millennium.

Transcriptome-wide microRNA and target dynamics in the fat body during the gonadotrophic cycle of Aedes aegypti

The mosquito Aedes aegypti is a major vector of numerous viral diseases, because it requires a blood meal to facilitate egg development. The fat body, a counterpart of mammalian liver and adipose tissues, is the metabolic center, playing a key role in reproduction. Therefore, understanding of regulatory networks controlling its functions is critical, and the role of microRNAs (miRNAs) in the process is largely unknown. We aimed to explore miRNA expression and potential targets in the female fat body of Ae. aegypti, as well as their changes postblood meal (PBM). Small RNA library analysis revealed five unique miRNA patterns sequentially expressed at five sampled time points, likely responding to, and affecting, waves of upstream hormonal signals and gene expression in the same period. To link miRNA identities with downstream targets, transcriptome-wide mRNA 3' UTR interaction sites were experimentally determined at 72 h posteclosion and 24 h PBM through Argonaute 1 cross-linking and immunoprecipitation followed by high-throughput sequencing. Several target sites were validated by means of in vitro luciferase assays with wild-type and mutated 3' UTRs for six miRNA families. With established transgenic lines, consistent results were observed with spatiotemporal knockdown of miR-8 and luciferase assays. We further investigated miRNAs potentially regulating various physiological processes based on Clusters of Orthologous Groups functional categories. Hence, the present work comprehensively elucidated miRNA expression and target dynamics in the female mosquito fat body, providing a solid foundation for future functional studies of miRNA regulation during the gonadotrophic cycle.

Transcriptome Sequencing Reveals Large-Scale Changes in Axenic Aedes aegypti Larvae

Mosquitoes host communities of microbes in their digestive tract that consist primarily of bacteria. We previously reported that Aedes aegypti larvae colonized by a native community of bacteria and gnotobiotic larvae colonized by only Escherichia coli develop very similarly into adults, whereas axenic larvae never molt and die as first instars. In this study, we extended these findings by first comparing the growth and abundance of bacteria in conventional, gnotobiotic, and axenic larvae during the first instar. Results showed that conventional and gnotobiotic larvae exhibited no differences in growth, timing of molting, or number of bacteria in their digestive tract. Axenic larvae in contrast grew minimally and never achieved the critical size associated with molting by conventional and gnotobiotic larvae. In the second part of the study we compared patterns of gene expression in conventional, gnotobiotic and axenic larvae by conducting an RNAseq analysis of gut and nongut tissues (carcass) at 22 h post-hatching. Approximately 12% of Ae. aegypti transcripts were differentially expressed in axenic versus conventional or gnotobiotic larvae. However, this profile consisted primarily of transcripts in seven categories that included the down-regulation of select peptidases in the gut and up-regulation of several genes in the gut and carcass with roles in amino acid transport, hormonal signaling, and metabolism. Overall, our results indicate that axenic larvae exhibit alterations in gene expression consistent with defects in acquisition and assimilation of nutrients required for growth.

RNA Interference for Mosquito and Mosquito-Borne Disease Control

RNA interference (RNAi) is a powerful tool to silence endogenous mosquito and mosquito-borne pathogen genes in vivo. As the number of studies utilizing RNAi in basic research grows, so too does the arsenal of physiological targets that can be developed into products that interrupt mosquito life cycles and behaviors and, thereby, relieve the burden of mosquitoes on human health and well-being. As this technology becomes more viable for use in beneficial and pest insect management in agricultural settings, it is exciting to consider its role in public health entomology. Existing and burgeoning strategies for insecticide delivery could be adapted to function as RNAi trigger delivery systems and thereby expedite transformation of RNAi from the lab to the field for mosquito control. Taken together, development of RNAi-based vector and pathogen management techniques & strategies are within reach. That said, tools for successful RNAi design, studies exploring RNAi in the context of vector control, and studies demonstrating field efficacy of RNAi trigger delivery have yet to be honed and/or developed for mosquito control.

PIWIs Go Viral: Arbovirus-Derived piRNAs in Vector Mosquitoes

Vector mosquitoes are responsible for transmission of the majority of arthropod-borne (arbo-) viruses. Virus replication in these vectors needs to be sufficiently high to permit efficient virus transfer to vertebrate hosts. The mosquito immune response therefore is a key determinant for arbovirus transmission. Mosquito antiviral immunity is primarily mediated by the small interfering RNA pathway. Besides this well-established antiviral machinery, the PIWI-interacting RNA (piRNA) pathway processes viral RNA into piRNAs. In recent years, significant progress has been made in characterizing the biogenesis and function of these viral piRNAs. In this review, we discuss these developments, identify knowledge gaps, and suggest directions for future research.

Identification of genes for engineering the male germline of Aedes aegypti and Ceratitis capitata

BACKGROUND

Synthetic biology approaches are promising new strategies for control of pest insects that transmit disease and cause agricultural damage. These strategies require characterised modular components that can direct appropriate expression of effector sequences, with components conserved across species being particularly useful. The goal of this study was to identify genes from which new potential components could be derived for manipulation of the male germline in two major pest species, the mosquito Aedes aegypti and the tephritid fruit fly Ceratitis capitata.

RESULTS

Using RNA-seq data from staged testis samples, we identified several candidate genes with testis-specific expression and suitable expression timing for use of their regulatory regions in synthetic control constructs. We also developed a novel computational pipeline to identify candidate genes with testis-specific splicing from this data; use of alternative splicing is another method for restricting expression in synthetic systems. Some of the genes identified display testis-specific expression or splicing that is conserved across species; these are particularly promising candidates for construct development.

CONCLUSIONS

In this study we have identified a set of genes with testis-specific expression or splicing. In addition to their interest from a basic biology perspective, these findings provide a basis from which to develop synthetic systems to control important pest insects via manipulation of the male germline.

Identification of Aedes aegypti Long Intergenic Non-coding RNAs and Their Association with Wolbachia and Dengue Virus Infection

Long intergenic non-coding RNAs (lincRNAs) are appearing as an important class of regulatory RNAs with a variety of biological functions. The aim of this study was to identify the lincRNA profile in the dengue vector Aedes aegypti and evaluate their potential role in host-pathogen interaction. The majority of previous RNA-Seq transcriptome studies in Ae. aegypti have focused on the expression pattern of annotated protein coding genes under different biological conditions. Here, we used 35 publically available RNA-Seq datasets with relatively high depth to screen the Ae. aegypti genome for lincRNA discovery. This led to the identification of 3,482 putative lincRNAs. These lincRNA genes displayed a slightly lower GC content and shorter transcript lengths compared to protein-encoding genes. Ae. aegypti lincRNAs also demonstrate low evolutionary sequence conservation even among closely related species such as Culex quinquefasciatus and Anopheles gambiae. We examined their expression in dengue virus serotype 2 (DENV-2) and Wolbachia infected and non-infected adult mosquitoes and Aa20 cells. The results revealed that DENV-2 infection increased the abundance of a number of host lincRNAs, from which some suppress viral replication in mosquito cells. RNAi-mediated silencing of lincRNA_1317 led to enhancement in viral replication, which possibly indicates its potential involvement in the host anti-viral defense. A number of lincRNAs were also differentially expressed in Wolbachia-infected mosquitoes. The results will facilitate future studies to unravel the function of lncRNAs in insects and may prove to be beneficial in developing new ways to control vectors or inhibit replication of viruses in them.

Developmental piRNA profiles of the invasive vector mosquito Aedes albopictus

Background

In eukaryotic organisms, Piwi-interacting RNAs (piRNAs) control the activities of mobile genetic elements and ensure genome maintenance. Recent evidence indicates that piRNAs are involved in multiple biological pathways, including transcriptional regulation of protein-coding genes, sex determination and even interactions between host and pathogens. Aedes albopictus is a major invasive species that transmits a number of viral diseases in humans. Ae. albopictus has the largest genome and the highest abundance of repetitive sequences when compared with members that belong to Culicidae with a published genome. Analysis of piRNA profiles will provide a developmental and evolutionary perspective on piRNAs in Ae. albopictus.

Methods

piRNAs were identified and characterized during the development of Ae. albopictus, and piRNA expression patterns in adult males and females as well as sugar-fed females and blood-fed females were compared.

Results

Our results reveal that, despite the large genome size of Ae. albopictus, the piRNA pool of Ae. albopictus (1.2 × 10⁷) is smaller than those of Aedes aegypti (1.7 × 10⁷) and Drosophila melanogaster (1.6 × 10⁷). In Ae. albopictus, piRNAs displayed the highest abundance at the embryo stage and the lowest abundance at the pupal stage. Approximately 50 % of the piRNAs mapped to intergenic regions with no known functions. Approximately 30 % of the piRNAs mapped to repetitive elements, and 77.69 % of these repeat-derived piRNAs mapped to Class I TEs; 45.42 % of the observed piRNA reads originated from piRNA clusters, and most of the top 10 highest expressed piRNA clusters and 100 highest expressed piRNAs from each stage displayed biased expression patterns across the developmental stages. All anti-sense-derived piRNAs displayed a preference for uridine at the 5′ end; however, the sense-derived piRNAs showed adenine bias at the tenth nucleotide position and a typical ping-pong signature, suggesting that the biogenesis of piRNAs was conserved throughout development. Our results also show that 962 piRNAs displayed sex-biased expression, and 522 piRNAs showed higher expression in the blood-fed females than in the sugar-fed females.

Conclusions

Our results suggest that piRNAs, aside from silencing transposable elements in Ae. albopictus, may have a role in other biological pathways.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1815-8) contains supplementary material, which is available to authorized users.

Mosquito Peptide Hormones: Diversity, Production, and Function

Mosquitoes, like other insects, produce a diversity of peptide hormones that are processed from different precursor proteins and have a range of activities. Early studies relied on purification of bioactive peptides for hormone identification, but more recently genomic data have provided the information needed to more comprehensively identify peptide hormone genes and associated receptors. The first part of this chapter summarizes the known or predicted peptide hormones that are produced by mosquitoes. The second part of this chapter discusses the sources of these molecules and their functions.

Massive Shift in Gene Expression during Transitions between Developmental Stages of the Gall Midge, Mayetiola Destructor

Mayetiola destructor is a destructive pest of wheat and has six developmental stages. Molecular mechanisms controlling the transition between developmental stages remain unknown. Here we analyzed genes that were expressed differentially between two successive developmental stages, including larvae at 1, 3, 5, and 7 days, pupae, and adults. A total of 17,344 genes were expressed during one or more of these studied stages. Among the expressed genes, 38-68% were differently expressed between two successive stages, with roughly equal percentages of up- and down-regulated genes. Analysis of the functions of the differentially expressed genes revealed that each developmental stage had some unique types of expressed genes that are characteristic of the physiology at that stage. This is the first genome-wide analysis of genes differentially expressed in different stages in a gall midge. The large dataset of up- and down-regulated genes in each stage of the insect shall be very useful for future research to elucidate mechanisms regulating insect development and other biological processes.

Virus-derived small RNAs: molecular footprints of host-pathogen interactions

Viruses are obligatory intracellular parasites that require the host machinery to replicate. During their replication cycle, viral RNA intermediates can be recognized and degraded by different antiviral mechanisms that include RNA decay, RNA interference, and RNase L pathways. As a consequence of viral RNA degradation, infected cells can accumulate virus‐derived small RNAs at high levels compared to cellular molecules. These small RNAs are imprinted with molecular characteristics that reflect their origin. First, small RNAs can be used to reconstruct viral sequences and identify the virus from which they originated. Second, other molecular features of small RNAs such as size, polarity, and base preferences depend on the type of viral substrate and host mechanism of degradation. Thus, the pattern of small RNAs generated in infected cells can be used as a molecular footprint to identify and characterize viruses independent on sequence homology searches against known references. Hence, sequencing of small RNAs obtained from infected cells enables virus discovery and characterization using both sequence‐dependent strategies and novel pattern‐based approaches. Recent studies are helping unlock the full application of small RNA sequencing for virus discovery and characterization. WIREs RNA 2016, 7:824–837. doi: 10.1002/wrna.1361

This article is categorized under: 1

RNA Processing > Processing of Small RNAs

2

RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms

3

RNA Methods > RNA Analyses In Vitro and In Silico

High-throughput cis-regulatory element discovery in the vector mosquito Aedes aegypti

Background

Despite substantial progress in mosquito genomic and genetic research, few cis-regulatory elements (CREs), DNA sequences that control gene expression, have been identified in mosquitoes or other non-model insects. Formaldehyde-assisted isolation of regulatory elements paired with DNA sequencing, FAIRE-seq, is emerging as a powerful new high-throughput tool for global CRE discovery. FAIRE results in the preferential recovery of open chromatin DNA fragments that are not bound by nucleosomes, an evolutionarily conserved indicator of regulatory activity, which are then sequenced. Despite the power of the approach, FAIRE-seq has not yet been applied to the study of non-model insects. In this investigation, we utilized FAIRE-seq to profile open chromatin and identify likely regulatory elements throughout the genome of the human disease vector mosquito Aedes aegypti. We then assessed genetic variation in the regulatory elements of dengue virus susceptible (Moyo-S) and refractory (Moyo-R) mosquito strains.

Results

Analysis of sequence data obtained through next generation sequencing of FAIRE DNA isolated from A. aegypti embryos revealed >121,000 FAIRE peaks (FPs), many of which clustered in the 1 kb 5’ upstream flanking regions of genes known to be expressed at this stage. As expected, known transcription factor consensus binding sites were enriched in the FPs, and of these FoxA1, Hunchback, Gfi, Klf4, MYB/ph3 and Sox9 are most predominant. All of the elements tested in vivo were confirmed to drive gene expression in transgenic Drosophila reporter assays. Of the >13,000 single nucleotide polymorphisms (SNPs) recently identified in dengue virus-susceptible and refractory mosquito strains, 3365 were found to map to FPs.

Conclusion

FAIRE-seq analysis of open chromatin in A. aegypti permitted genome-wide discovery of CREs. The results of this investigation indicate that FAIRE-seq is a powerful tool for identification of regulatory DNA in the genomes of non-model organisms, including human disease vector mosquitoes.

Electronic supplementary material

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

Laterally Transferred Gene Recruited as a Venom in Parasitoid Wasps

Parasitoid wasps use venom to manipulate the immunity and metabolism of their host insects in a variety of ways to provide resources for their offspring. Yet, how genes are recruited and evolve to perform venom functions remain open questions. A recently recognized source of eukaryotic genome innovation is lateral gene transfer (LGT). Glycoside hydrolase family 19 (GH19) chitinases are widespread in bacteria, microsporidia, and plants where they are used in nutrient acquisition or defense, but have previously not been known in metazoans. In this study, a GH19 chitinase LGT is described from the unicellular microsporidia/Rozella clade into parasitoid wasps of the superfamily Chalcidoidea, where it has become recruited as a venom protein. The GH19 chitinase is present in 15 species of chalcidoid wasps representing four families, and phylogenetic analysis indicates that it was laterally transferred near or before the origin of Chalcidoidea (∼95 Ma). The GH19 chitinase gene is highly expressed in the venom gland of at least seven species, indicating a role in the complex host manipulations performed by parasitoid wasp venom. RNAi knockdown in the model parasitoid Nasonia vitripennis reveals that-following envenomation-the GH19 chitinase induces fly hosts to upregulate genes involved in an immune response to fungi. A second, independent LGT of GH19 chitinase from microsporidia into mosquitoes was also found, also supported by phylogenetic reconstructions. Besides these two LGT events, GH19 chitinase is not found in any other sequenced animal genome, or in any fungi outside the microsporidia/Rozella clade.

Calcium influx enhances neuropeptide activation of ecdysteroid hormone production by mosquito ovaries

A critical step in mosquito reproduction is the ingestion of a blood meal from a vertebrate host. In mosquitoes like Aedes aegypti, blood feeding stimulates the release of ovary ecdysteroidogenic hormone (OEH) and insulin-like peptide 3 (ILP3). This induces the ovaries to produce ecdysteroid hormone (ECD), which then drives egg maturation. In many immature insects, prothoracicotropic hormone (PTTH) stimulates the prothoracic glands to produce ECD that directs molting and metamorphosis. The receptors for OEH, ILP3 and PTTH are different receptor tyrosine kinases with OEH and ILP3 signaling converging downstream in the insulin pathway and PTTH activating the mitogen-activated protein kinase pathway. Calcium (Ca(2+)) flux and cAMP have also been implicated in PTTH signaling, but the role of Ca(2+) in OEH, ILP3, and cAMP signaling in ovaries is unknown. Here, we assessed whether Ca(2+) flux affects OEH, ILP3, and cAMP activity in A. aegypti ovaries and also asked whether PTTH stimulated ovaries to produce ECD. Results indicated that Ca(2+) flux enhanced but was not essential for OEH or ILP3 activity, whereas cAMP signaling was dependent on Ca(2+) flux. Recombinant PTTH from Bombyx mori fully activated ECD production by B. mori PTGs, but exhibited no activity toward A. aegypti ovaries. Recombinant PTTH from A. aegypti also failed to stimulate either B. mori PTGs or A. aegypti ovaries to produce ECD. We discuss the implications of these results in the context of mosquito reproduction and ECD biosynthesis by insects generally.

Small RNA Profiling in Dengue Virus 2-Infected Aedes Mosquito Cells Reveals Viral piRNAs and Novel Host miRNAs

In Aedes mosquitoes, infections with arthropod-borne viruses (arboviruses) trigger or modulate the expression of various classes of viral and host-derived small RNAs, including small interfering RNAs (siRNAs), PIWI interacting RNAs (piRNAs), and microRNAs (miRNAs). Viral siRNAs are at the core of the antiviral RNA interference machinery, one of the key pathways that limit virus replication in invertebrates. Besides siRNAs, Aedes mosquitoes and cells derived from these insects produce arbovirus-derived piRNAs, the best studied examples being viruses from the Togaviridae or Bunyaviridae families. Host miRNAs modulate the expression of a large number of genes and their levels may change in response to viral infections. In addition, some viruses, mostly with a DNA genome, express their own miRNAs to regulate host and viral gene expression. Here, we perform a comprehensive analysis of both viral and host-derived small RNAs in Aedes aegypti Aag2 cells infected with dengue virus 2 (DENV), a member of the Flaviviridae family. Aag2 cells are competent in producing all three types of small RNAs and provide a powerful tool to explore the crosstalk between arboviral infection and the distinct RNA silencing pathways. Interestingly, besides the well-characterized DENV-derived siRNAs, a specific population of viral piRNAs was identified in infected Aag2 cells. Knockdown of Piwi5, Ago3 and, to a lesser extent, Piwi6 results in reduction of vpiRNA levels, providing the first genetic evidence that Aedes PIWI proteins produce DENV-derived small RNAs. In contrast, we do not find convincing evidence for the production of virus-derived miRNAs. Neither do we find that host miRNA expression is strongly changed upon DENV2 infection. Finally, our deep-sequencing analyses detect 30 novel Aedes miRNAs, complementing the repertoire of regulatory small RNAs in this important vector species.

Mosquitoes of the Aedes family transmit many important viruses, including dengue virus, between their vertebrate hosts. In the mosquito, the growth of these viruses is limited by the antiviral RNA interference pathway. Key to this pathway is a class of small non-coding RNAs known as small interfering RNAs (siRNAs). In addition, two related but distinct small RNA pathways known as the microRNA (miRNA) and the PIWI-interacting RNA (piRNA) pathway are implicated in regulating virus replication in mosquitoes. Thus, since small RNAs may critically influence the transmission of dengue virus, we set out to analyze the populations of viral and mosquito small RNAs that are produced in infected Aedes mosquito cells. We found that besides the well-known viral siRNAs, dengue virus-derived piRNAs were produced in these cells and we identified the PIWI proteins that these small RNAs rely on. In addition, we found that viral miRNAs were not expressed from the dengue virus genome and that the levels of mosquito miRNAs were barely changed upon infection. Finally, our data allowed for the identification of novel Aedes miRNAs, complementing the repertoire of these important regulatory RNAs in vector mosquitoes.

Mating-Induced Transcriptome Changes in the Reproductive Tract of Female Aedes aegypti

The Aedes aegypti mosquito is a significant public health threat, as it is the main vector of dengue and chikungunya viruses. Disease control efforts could be enhanced through reproductive manipulation of these vectors. Previous work has revealed a relationship between male seminal fluid proteins transferred to females during mating and female post-mating physiology and behavior. To better understand this interplay, we used short-read RNA sequencing to identify gene expression changes in the lower reproductive tract of females in response to mating. We characterized mRNA expression in virgin and mated females at 0, 6 and 24 hours post-mating (hpm) and identified 364 differentially abundant transcripts between mating status groups. Surprisingly, 60 transcripts were more abundant at 0hpm compared to virgin females, suggesting transfer from males. Twenty of these encode known Ae. aegypti seminal fluid proteins. Transfer and detection of male accessory gland-derived mRNA in females at 0hpm was confirmed by measurement of eGFP mRNA in females mated to eGFP-expressing males. In addition, 150 transcripts were up-regulated at 6hpm and 24hpm, while 130 transcripts were down-regulated at 6hpm and 24hpm. Gene Ontology (GO) enrichment analysis revealed that proteases, a protein class broadly known to play important roles in reproduction, were among the most enriched protein classes. RNAs associated with immune system and antimicrobial function were also up-regulated at 24hpm. Collectively, our results suggest that copulation initiates broad transcriptome changes across the mosquito female reproductive tract, “priming” her for important subsequent processes of blood feeding, egg development and immune defense. Our transcriptome analysis provides a vital foundation for future studies of the consequences of mating on female biology and will aid studies seeking to identify specific gene families, molecules and pathways that support key reproductive processes in the female mosquito.

Female post-mating behavior has important consequences for mosquito populations and their ability to transmit diseases. Male Aedes aegypti seminal fluid substances transferred during mating cause many important changes to female behavior and physiology, including blood feeding behavior, egg development, and oviposition. In an effort to understand how males induce these responses in Ae. aegypti females, we characterized the transcriptome changes that occur in the female reproductive tract at different time points after mating. We found several RNAs that are apparently transferred by the male, and 280 genes whose mRNA abundance in the female is affected by mating. The nature of the predicted products of many of these genes suggests roles in priming the reproductive tract for egg development, protecting the female against bacterial infections or processing the blood meal. This identification of mating-responsive genes provides information potentially useful for developing tools aimed at preventing disease transmission by manipulating female mosquitoes’ post-mating responses.

The evolutionary landscape of intergenic trans-splicing events in insects

To explore the landscape of intergenic trans-splicing events and characterize their functions and evolutionary dynamics, we conduct a mega-data study of a phylogeny containing eight species across five orders of class Insecta, a model system spanning 400 million years of evolution. A total of 1,627 trans-splicing events involving 2,199 genes are identified, accounting for 1.58% of the total genes. Homology analysis reveals that mod(mdg4)-like trans-splicing is the only conserved event that is consistently observed in multiple species across two orders, which represents a unique case of functional diversification involving trans-splicing. Thus, evolutionarily its potential for generating proteins with novel function is not broadly utilized by insects. Furthermore, 146 non-mod trans-spliced transcripts are found to resemble canonical genes from different species. Trans-splicing preserving the function of 'breakup' genes may serve as a general mechanism for relaxing the constraints on gene structure, with profound implications for the evolution of genes and genomes.

Gene expression changes in the salivary glands of Anopheles coluzzii elicited by Plasmodium berghei infection

Background

Malaria is a devastating infectious disease caused by Plasmodium parasites transmitted through the bites of infected Anopheles mosquitoes. Salivary glands are the only mosquito tissue invaded by Plasmodium sporozoites, being a key stage for the effective parasite transmission, making the study of Anopheles sialome highly relevant.

Methods

RNA-sequencing was used to compare differential gene expression in salivary glands of uninfected and Plasmodium berghei-infected Anopheles coluzzii mosquitoes. RNA-seq results were validated by quantitative RT-PCR. The transmembrane glucose transporter gene AGAP007752 was selected for functional analysis by RNA interference. The effect of gene silencing on infection level was evaluated. The putative function and tertiary structure of the protein was assessed.

Results

RNA-seq data showed that 2588 genes were differentially expressed in mosquitoes salivary glands in response to P. berghei infection, being 1578 upregulated and 1010 downregulated. Metabolism, Immunity, Replication/Transcription/Translation, Proteolysis and Transport were the mosquito gene functional classes more affected by parasite infection. Endopeptidase coding genes were the most abundant within the differentially expressed genes in infected salivary glands (P < 0.001). Based on its putative function and expression level, the transmembrane glucose transporter gene, AGAP007752, was selected for functional analysis by RNA interference. The results demonstrated that the number of sporozoites was 44.3 % lower in mosquitoes fed on infected mice after AGAPP007752 gene knockdown when compared to control (P < 0.01).

Conclusions

Our hypothesis is that the protein encoded by the gene AGAPP007752 may play a role on An. coluzzii salivary glands infection by Plasmodium parasite, working as a sporozoite receptor and/or promoting a favorable environment for the capacity of sporozoites.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1079-8) contains supplementary material, which is available to authorized users.

Genetic technologies for disease vectors

The first genetic technologies for insect vectors of disease were introduced 20 years ago. As of today there are 12 classes of genetic technologies used as functional genomic tools for insect vectors of important diseases. Although the applications of genetic technologies in insect disease vectors have been conducted primarily in mosquitoes, other insect systems could benefit from current technologies. While the various technological platforms are likely to function in diverse arthropods, the delivery of these technologies to cells and tissues of interest is the major technical constraint that limits their widespread adoption. Increased community resources of various types would enhance the adoption of these technologies and potentially eliminate technical limitations.

Distinct sets of PIWI proteins produce arbovirus and transposon-derived piRNAs in Aedes aegypti mosquito cells

The PIWI-interacting RNA (piRNA) pathway is essential for transposon silencing in many model organisms. Its remarkable efficiency relies on a sophisticated amplification mechanism known as the ping-pong loop. In Alphavirus-infected Aedes mosquitoes, piRNAs with sequence features that suggest ping-pong-dependent biogenesis are produced from viral RNA. The PIWI family in Aedes mosquitoes is expanded when compared to other model organisms, raising the possibility that individual PIWI proteins have functionally diversified in these insects. Here, we show that Piwi5 and Ago3, but none of the other PIWI family members, are essential for piRNA biogenesis from Sindbis virus RNA in infected Aedes aegypti cells. In contrast, the production of piRNAs from transposons relies on a more versatile set of PIWI proteins, some of which do not contribute to viral piRNA biogenesis. These results indicate that functional specialization allows distinct mosquito PIWI proteins to process RNA from different endogenous and exogenous sources.

SEX DETERMINATION. A male-determining factor in the mosquito Aedes aegypti

Sex determination in the mosquito Aedes aegypti is governed by a dominant male-determining factor (M factor) located within a Y chromosome-like region called the M locus. Here, we show that an M-locus gene, Nix, functions as an M factor in A. aegypti. Nix exhibits persistent M linkage and early embryonic expression, two characteristics required of an M factor. Nix knockout with clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 resulted in largely feminized genetic males and the production of female isoforms of two key regulators of sexual differentiation: doublesex and fruitless. Ectopic expression of Nix resulted in genetic females with nearly complete male genitalia. Thus, Nix is both required and sufficient to initiate male development. This study provides a foundation for mosquito control strategies that convert female mosquitoes into harmless males.

Ovary ecdysteroidogenic hormone requires a receptor tyrosine kinase to activate egg formation in the mosquito Aedes aegypti

Mosquitoes are major disease vectors because most species must feed on blood from a vertebrate host to produce eggs. Blood feeding by the vector mosquito Aedes aegypti triggers the release of two neurohormones, ovary ecdysteroidogenic hormone (OEH) and insulin-like peptides (ILPs), which activate multiple processes required for egg formation. ILPs function by binding to the insulin receptor, which activates downstream components in the canonical insulin signaling pathway. OEH in contrast belongs to a neuropeptide family called neuroparsins, whose receptor is unknown. Here we demonstrate that a previously orphanized receptor tyrosine kinase (RTK) from A. aegypti encoded by the gene AAEL001915 is an OEH receptor. Phylogenetic studies indicated that the protein encoded by this gene, designated AAEL001915, belongs to a clade of RTKs related to the insulin receptor, which are distinguished by an extracellular Venus flytrap module. Knockdown of AAEL001915 by RNAi disabled OEH-mediated egg formation in A. aegypti. AAEL001915 was primarily detected in the mosquito ovary in association with follicular epithelial cells. Both monomeric and dimeric AAEL001915 were detected in mosquito ovaries and transfected Drosophila S2 cells. Functional assays further indicated that OEH bound to dimeric AAEL001915, which resulted in downstream phosphorylation of Ak strain transforming factor (Akt). We hypothesize that orthologs of AAEL001915 in other insects are neuroparsin receptors.

MicroRNAs of two medically important mosquito species: Aedes aegypti and Anopheles stephensi

MicroRNAs (miRNAs) are endogenous, single-stranded small RNAs that have important regulatory functions at the post-transcriptional level. In the present study, we characterize miRNAs in two divergent mosquito species, Aedes aegypti and Anopheles stephensi, through deep sequencing of small RNAs spanning all developmental stages. We discovered eight novel miRNAs in Ae. aegypti and 20 novel miRNAs in An. stephensi, which enabled the first systematic analysis of miRNA evolution in mosquitos. We traced the phylogenetic history of all miRNAs in both species and report a rate of 0.055-0.13 miRNA net gain per million years. Most novel miRNAs originate de novo. Duplications that produced miRNA clusters and families are more common in Ae. aegypti than in An. stephensi. We also identified arm-switch as a source of new miRNAs. Expression profile analysis identified mosquito-specific miRNAs that showed strong stage-specific expression in one or both lineages. For example, the aae-miR-2941/2946 family represents the most abundant maternally deposited and zygotically transcribed miRNAs in Ae. aegypti. miR-2943 is a highly expressed zygotic miRNA in both Ae. aegypti and An. stephensi. Such information provides the basis from which to study the function of these miRNAs in biology common to all mosquitos or unique to one particular lineage.

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.

Transcriptome analysis between invasive Pomacea canaliculata and indigenous Cipangopaludina cahayensis reveals genomic divergence and diagnostic microsatellite/SSR markers

BACKGROUND

Pomacea canaliculata is an important invasive species worldwide. However, little is known about the molecular mechanisms behind species displacement, adaptational abilities, and pesticide resistance, partly because of the lack of genomic information that is available for this species. Here, the transcriptome sequences for the invasive golden apple snail P. canaliculata and the native mudsnail Cipangopaludina cahayensis were obtained by next-generation-sequencing and used to compare genomic divergence and identify molecular markers.

RESULTS

More than 46 million high quality sequencing reads were generated from P. canaliculata and C. cahayensis using Illumina paired-end sequencing technology. Our analysis indicated that 11,312 unigenes from P. canaliculata and C. cahayensis showed significant similarities to known proteins families, among which a total of 4,320 specific protein families were identified. KEGG pathway enrichment was analyzed for the unique unigenes with 17 pathways (p-value < 10(-5)) in P. canaliculata relating predominantly to lysosomes and vitamin digestion and absorption, and with 12 identified in C. cahayensis, including cancer and toxoplasmosis pathways, respectively. Our analysis also indicated that the comparatively high number of P450 genes in the P. canaliculata transcriptome may be associated with the pesticide resistance in this species. Additionally, 16,717 simple sequence repeats derived from expressed sequence tags (EST-SSRs) were identified from the 14,722 unigenes in P. canaliculata and 100 of them were examined by PCR, revealing a species-specific molecular marker that could distinguish between the morphologically similar P. canaliculata and C. cahayensis snails.

CONCLUSIONS

Here, we present the genomic resources of P. canaliculata and C. cahayensis. Differentially expressed genes in the transcriptome of P. canaliculata compared with C. cahayensis corresponded to critical metabolic pathways, and genes specifically related to environmental stress response were detected. The CYP4 family of P450 cytochromes that may be important factors in pesticide metabolism in P. canaliculata was identified. Overall, these findings will provide valuable genetic data for the further characterization of the molecular mechanisms that support the invasive and adaptive abilities of P. canaliculata.

Maternal germline-specific genes in the Asian malaria mosquito Anopheles stephensi: characterization and application for disease control

Anopheles stephensi is a principal vector of urban malaria on the Indian subcontinent and an emerging model for molecular and genetic studies of mosquito biology. To enhance our understanding of female mosquito reproduction, and to develop new tools for basic research and for genetic strategies to control mosquito-borne infectious diseases, we identified 79 genes that displayed previtellogenic germline-specific expression based on RNA-Seq data generated from 11 life stage-specific and sex-specific samples. Analysis of this gene set provided insights into the biology and evolution of female reproduction. Promoters from two of these candidates, vitellogenin receptor and nanos, were used in independent transgenic cassettes for the expression of artificial microRNAs against suspected mosquito maternal-effect genes, discontinuous actin hexagon and myd88. We show these promoters have early germline-specific expression and demonstrate 73% and 42% knockdown of myd88 and discontinuous actin hexagon mRNA in ovaries 48 hr after blood meal, respectively. Additionally, we demonstrate maternal-specific delivery of mRNA and protein to progeny embryos. We discuss the application of this system of maternal delivery of mRNA/miRNA/protein in research on mosquito reproduction and embryonic development, and for the development of a gene drive system based on maternal-effect dominant embryonic arrest.

piRNA pathway gene expression in the malaria vector mosquito Anopheles stephensi

The ability of transposons to mobilize to new places in a genome enables them to introgress rapidly into populations. The piRNA pathway has been characterized recently in the germ line of the fruit fly, Drosophila melanogaster, and is responsible for downregulating transposon mobility. Transposons have been used as tools in mosquitoes to genetically transform a number of species including Anopheles stephensi, a vector of human malaria. These mobile genetic elements also have been proposed as tools to drive antipathogen effector genes into wild mosquito populations to replace pathogen-susceptible insects with those engineered genetically to be resistant to or unable to transmit a pathogen. The piRNA pathway may affect the performance of such proposed genetic engineering strategies. In the present study, we identify and describe the An. stephensi orthologues of the major genes in the piRNA pathway, Ago3, Aubergine (Aub) and Piwi. Consistent with a role in protection from transposon movement, these three genes are expressed constitutively in the germ-line cells of ovaries and induced further after a blood meal.

A role for Drosophila Cyclin J in oogenesis revealed by genetic interactions with the piRNA pathway

Cyclin J (CycJ) is a poorly characterized member of the Cyclin superfamily of cyclin-dependent kinase regulators, many of which regulate the cell cycle or transcription. Although CycJ is conserved in metazoans its cellular function has not been identified and no mutant defects have been described. In Drosophila, CycJ transcript is present primarily in ovaries and very early embryos, suggesting a role in one or both of these tissues. The CycJ gene (CycJ) lies immediately downstream of armitage (armi), a gene involved in the Piwi-associated RNA (piRNA) pathways that are required for silencing transposons in the germline and adjacent somatic cells. Mutations in armi result in oogenesis defects but a role for CycJ in oogenesis has not been defined. Here we assessed oogenesis in CycJ mutants in the presence or absence of mutations in armi or other piRNA pathway genes. CycJ null ovaries appeared normal, indicating that CycJ is not essential for oogenesis under normal conditions. In contrast, armi null ovaries produced only two egg chambers per ovariole and the eggs had severe axis specification defects, as observed previously for armi and other piRNA pathway mutants. Surprisingly, the CycJ armi double mutant failed to produce any mature eggs. The double null ovaries generally had only one egg chamber per ovariole and the egg chambers frequently contained an overabundance of differentiated germline cells. Production of these compound egg chambers could be suppressed with CycJ transgenes but not with mutations in the checkpoint gene mnk, which suppress oogenesis defects in armi mutants. The CycJ null showed similar genetic interactions with the germline and somatic piRNA pathway gene piwi, and to a lesser extent with aubergine (aub), a member of the germline-specific piRNA pathway. The strong genetic interactions between CycJ and piRNA pathway genes reveal a role for CycJ in early oogenesis. Our results suggest that CycJ is required to regulate egg chamber production or maturation when piRNA pathways are compromised.

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.

Integrated proteomic and transcriptomic analysis of the Aedes aegypti eggshell

Background

Mosquito eggshells show remarkable diversity in physical properties and structure consistent with adaptations to the wide variety of environments exploited by these insects. We applied proteomic, transcriptomic, and hybridization in situ techniques to identify gene products and pathways that participate in the assembly of the Aedes aegypti eggshell. Aedes aegypti population density is low during cold and dry seasons and increases immediately after rainfall. The survival of embryos through unfavorable periods is a key factor in the persistence of their populations. The work described here supports integrated vector control approaches that target eggshell formation and result in Ae. aegypti drought-intolerant phenotypes for public health initiatives directed to reduce mosquito-borne diseases.

Results

A total of 130 proteins were identified from the combined mass spectrometric analyses of eggshell preparations.

Conclusions

Classification of proteins according to their known and putative functions revealed the complexity of the eggshell structure. Three novel Ae. aegypti vitelline membrane proteins were discovered. Odorant-binding and cysteine-rich proteins that may be structural components of the eggshell were identified. Enzymes with peroxidase, laccase and phenoloxidase activities also were identified, and their likely involvements in cross-linking reactions that stabilize the eggshell structure are discussed.

An analysis of the small RNA transcriptome of four developmental stages of the citrus red mite (Panonychus citri)

The citrus red mite (Panonychus citri) can feed on more than 112 plant species around the world. Endogenous small RNAs (sRNAs) have proved to be important components of gene regulation in many eukaryotes. Recently, many sRNAs have been shown to be involved in various biological processes, such as development in many animals, including insects; however, to date, no sRNAs have been reported in the citrus red mite. Using Illumina sequencing, several categories of sRNAs were identified, including 594 known microRNAs (miRNAs) grouped into 206 families and 31 novel miRNAs in the four developmental stages of citrus red mite. In addition, according to bioinformatics analysis and S-Poly(T) miRNA assays, the expression level of many miRNAs varied among the developmental stages. Furthermore, the prediction of miRNAs target genes and their functional annotation indicated that miRNAs are involved in the regulation of multiple pathways in the citrus red mite. As the first report of the sRNA world in citrus red mite, the present study furthers our understanding of the roles played by sRNAs in the development of citrus red mite and the data may help to develop methods of controlling the pests in the field.

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.

BLIND ordering of large-scale transcriptomic developmental timecourses

RNA-Seq enables the efficient transcriptome sequencing of many samples from small amounts of material, but the analysis of these data remains challenging. In particular, in developmental studies, RNA-Seq is challenged by the morphological staging of samples, such as embryos, since these often lack clear markers at any particular stage. In such cases, the automatic identification of the stage of a sample would enable previously infeasible experimental designs. Here we present the 'basic linear index determination of transcriptomes' (BLIND) method for ordering samples comprising different developmental stages. The method is an implementation of a traveling salesman algorithm to order the transcriptomes according to their inter-relationships as defined by principal components analysis. To establish the direction of the ordered samples, we show that an appropriate indicator is the entropy of transcriptomic gene expression levels, which increases over developmental time. Using BLIND, we correctly recover the annotated order of previously published embryonic transcriptomic timecourses for frog, mosquito, fly and zebrafish. We further demonstrate the efficacy of BLIND by collecting 59 embryos of the sponge Amphimedon queenslandica and ordering their transcriptomes according to developmental stage. BLIND is thus useful in establishing the temporal order of samples within large datasets and is of particular relevance to the study of organisms with asynchronous development and when morphological staging is difficult.

Identification of germline transcriptional regulatory elements in Aedes aegypti

The mosquito Aedes aegypti is the principal vector for the yellow fever and dengue viruses, and is also responsible for recent outbreaks of the alphavirus chikungunya. Vector control strategies utilizing engineered gene drive systems are being developed as a means of replacing wild, pathogen transmitting mosquitoes with individuals refractory to disease transmission, or bringing about population suppression. Several of these systems, including Medea, UD(MEL), and site-specific nucleases, which can be used to drive genes into populations or bring about population suppression, utilize transcriptional regulatory elements that drive germline-specific expression. Here we report the identification of multiple regulatory elements able to drive gene expression specifically in the female germline, or in the male and female germline, in the mosquito Aedes aegypti. These elements can also be used as tools with which to probe the roles of specific genes in germline function and in the early embryo, through overexpression or RNA interference.

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.