A piRNA-lncRNA regulatory network initiates responder and trailer piRNA formation during mosquito embryonic development

Long non-coding RNAs as promising targets for controlling disease vector mosquitoes

Hematophagous female mosquitoes are important vectors of numerous devastating human diseases, posing a major public health threat. Effective prevention and control of mosquito-borne diseases rely considerably on progress in understanding the molecular mechanisms of various life activities, and accordingly, the molecules that regulate the various life activities of mosquitoes are potential targets for implementing future vector control strategies. Many long non-coding RNAs (lncRNAs) have been identified in mosquitoes and significant progress has been made in determining their functions. Here, we present a comprehensive overview of the research advances on mosquito lncRNAs, including their molecular identification, function, and interaction with other non-coding RNAs, as well as their synergistic regulatory roles in mosquito life activities. We also highlight the potential roles of competitive endogenous RNAs in mosquito growth and development, as well as in insecticide resistance and virus-host interactions. Insights into the biological functions and mechanisms of lncRNAs in mosquito life activities, viral replication, pathogenesis, and transmission will contribute to the development of novel drugs and safe vaccines.

Non-gonadal expression of piRNAs is widespread across Arthropoda

PIWI-interacting RNAs (piRNAs) were discovered in the early 2000s and became known for their role in protecting the germline genome against mobile genetic elements. Successively, piRNAs were also detected in the somatic cells of gonads in multiple animal species. In recent years, piRNAs have been reported in non-gonadal tissues in various arthropods, contrary to the initial assumptions of piRNAs being exclusive to gonads. Here, we performed an extensive literature review, which revealed that reports on non-gonadal somatic piRNA expression are not limited to a few specific species. Instead, when multiple studies are considered collectively, it appears to be a widespread phenomenon across arthropods. Furthermore, we systematically analyzed 168 publicly available small RNA-seq datasets from diverse tissues in 17 species, which further supported the bibliographic reports that piRNAs are expressed across tissues and species in Arthropoda.

Unraveling the modulatory manner and function of circRNAs in the Asian honey bee larval guts

Circular RNAs (circRNAs) are a class of non-coding RNAs (ncRNAs) that can participate in biological processes such as gene expression, growth, and development. However, little has been explored about the function of circRNAs in the development of Apis cerana larval guts. By using our previously gained deep sequencing data from the guts of A. cerana worker larvae at 4-, 5-, and 6-day-old (Ac4, Ac5, and Ac6 groups), the expression pattern and regulatory role of circular RNAs (circRNAs) during the development process was comprehensively investigated, with a focus on differentially expressed circRNAs (DEcircRNAs) relevant to immunity pathways and developmental signaling pathways, followed by validation of the binding relationships among a key competing endogenous RNA (ceRNA) axis. Here, 224 (158) DEcircRNAs were detected in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group. It's suggested that 172 (123) parental genes of DEcircRNAs were involved in 26 (20) GO terms such as developmental process and metabolic process and 138 (136) KEGG pathways like Hippo and Wnt signaling pathways. Additionally, ceRNA network analysis indicated that 21 (11) DEcircRNAs could target seven (three) DEmiRNAs, further targeting 324 (198) DEmRNAs. These DEmRNAs can be annotated to 33 (26) GO terms and 168 (200) KEGG pathways, including 12 (16) cellular and humoral immune pathways (endocytosis, lysosome, Jak-STAT, etc.) and 10 (nine) developmental signaling pathways (Hippo, mTOR, Hedgehog, etc.). Interestingly, DEcircRNAs in these two comparison groups could target the same ace-miR-6001-y, forming complex sub-networks. The results of PCR and Sanger sequencing confirmed the back-splicing sites within four randomly selected DEcircRNAs. RT-qPCR detection of these four DEcircRNAs verified the reliability of the used transcriptome data. The results of dual-luciferase reporter assay verified the binding relationships between novel_circ_001627 and ace-miR-6001-y and between ace-miR-6001-y and apterous-like. Our data demonstrated that DEcircRNAs were likely to modulate the developmental process of the A. cerana worker larval guts via regulation of parental gene transcription and ceRNA network, and novel_circ_001627/ace-miR-6001-y/apterous-like was a potential regulatory axis in the larval gut development. Findings from this work offer a basis and a candidate ceRNA axis for illustrating the circRNA-modulated mechanisms underlying the A. cerana larval guts.

Transposon and Transgene Tribulations in Mosquitoes: A Perspective of piRNA Proportions

Mosquitoes, like Drosophila, are dipterans, the order of "true flies" characterized by a single set of two wings. Drosophila are prime model organisms for biomedical research, while mosquito researchers struggle to establish robust molecular biology in these that are arguably the most dangerous vectors of human pathogens. Both insects utilize the RNA interference (RNAi) pathway to generate small RNAs to silence transposons and viruses, yet details are emerging that several RNAi features are unique to each insect family, such as how culicine mosquitoes have evolved extreme genomic feature differences connected to their unique RNAi features. A major technical difference in the molecular genetic studies of these insects is that generating stable transgenic animals are routine in Drosophila but still variable in stability in mosquitoes, despite genomic DNA-editing advances. By comparing and contrasting the differences in the RNAi pathways of Drosophila and mosquitoes, in this review we propose a hypothesis that transgene DNAs are possibly more intensely targeted by mosquito RNAi pathways and chromatin regulatory pathways than in Drosophila. We review the latest findings on mosquito RNAi pathways, which are still much less well understood than in Drosophila, and we speculate that deeper study into how mosquitoes modulate transposons and viruses with Piwi-interacting RNAs (piRNAs) will yield clues to improving transgene DNA expression stability in transgenic mosquitoes.

Diverse Regulatory Manners and Potential Roles of lncRNAs in the Developmental Process of Asian Honey Bee (Apis cerana) Larval Guts

Long non-coding RNAs (lncRNAs) are crucial modulators in a variety of biological processes, such as gene expression, development, and immune defense. However, little is known about the function of lncRNAs in the development of Asian honey bee (Apis cerana) larval guts. Here, on the basis of our previously obtained deep-sequencing data from the 4-, 5-, and 6-day-old larval guts of A. cerana workers (Ac4, Ac5, and Ac6 groups), an in-depth transcriptome-wide investigation was conducted to decipher the expression pattern, regulatory manners, and potential roles of lncRNAs during the developmental process of A. cerana worker larval guts, followed by the verification of the relative expression of differentially expressed lncRNAs (DElncRNAs) and the targeting relationships within a competing endogenous RNA (ceRNA) axis. In the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, 527 and 498 DElncRNAs were identified, respectively, which is suggestive of the dynamic expression of lncRNAs during the developmental process of larval guts. A cis-acting analysis showed that 330 and 393 neighboring genes of the aforementioned DElncRNAs were respectively involved in 29 and 32 functional terms, such as cellular processes and metabolic processes; these neighboring genes were also respectively engaged in 246 and 246 pathways such as the Hedgehog signaling pathway and the Wnt signaling pathway. Additionally, it was found that 79 and 76 DElncRNAs as potential antisense lncRNAs may, respectively, interact with 72 and 60 sense-strand mRNAs. An investigation of competing endogenous RNA (ceRNA) networks suggested that 75 (155) DElncRNAs in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group could target 7 (5) DEmiRNAs and further bind to 334 (248) DEmRNAs, which can be annotated to 33 (29) functional terms and 186 (210) pathways, including 12 (16) cellular- and humoral-immune pathways (lysosome pathway, necroptosis, MAPK signaling pathway, etc.) and 11 (10) development-associated signaling pathways (Wnt, Hippo, AMPK, etc.). The RT-qPCR detection of five randomly selected DElncRNAs confirmed the reliability of the used sequencing data. Moreover, the results of a dual-luciferase reporter assay were indicative of the binding relationship between MSTRG.11294.1 and miR-6001-y and between miR-6001-y and ncbi_107992440. These results demonstrate that DElncRNAs are likely to modulate the developmental process of larval guts via the regulation of the source genes' transcription, interaction with mRNAs, and ceRNA networks. Our findings not only yield new insights into the developmental mechanism underlying A. cerana larval guts, but also provide a candidate ceRNA axis for further functional dissection.

Chromatin profiling identifies transcriptional readthrough as a conserved mechanism for piRNA biogenesis in mosquitoes

The piRNA pathway in mosquitoes differs substantially from other model organisms, with an expanded PIWI gene family and functions in antiviral defense. Here, we define core piRNA clusters as genomic loci that show ubiquitous piRNA expression in both somatic and germline tissues. These core piRNA clusters are enriched for non-retroviral endogenous viral elements (nrEVEs) in antisense orientation and depend on key biogenesis factors, Veneno, Tejas, Yb, and Shutdown. Combined transcriptome and chromatin state analyses identify transcriptional readthrough as a conserved mechanism for cluster-derived piRNA biogenesis in the vector mosquitoes Aedes aegypti, Aedes albopictus, Culex quinquefasciatus, and Anopheles gambiae. Comparative analyses between the two Aedes species suggest that piRNA clusters function as traps for nrEVEs, allowing adaptation to environmental challenges such as virus infection. Our systematic transcriptome and chromatin state analyses lay the foundation for studies of gene regulation, genome evolution, and piRNA function in these important vector species.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

Culex Mosquito Piwi4 Is Antiviral against Two Negative-Sense RNA Viruses

Culex spp. mosquitoes transmit several pathogens concerning public health, including West Nile virus and Saint Louis encephalitis virus. Understanding the antiviral immune system of Culex spp. mosquitoes is important for reducing the transmission of these viruses. Mosquitoes rely on RNA interference (RNAi) to control viral replication. While the siRNA pathway in mosquitoes is heavily studied, less is known about the piRNA pathway. The piRNA pathway in mosquitoes has recently been connected to mosquito antiviral immunity. In Aedes aegypti, Piwi4 has been implicated in antiviral responses. The antiviral role of the piRNA pathway in Culex spp. mosquitoes is understudied compared to Ae. aegypti. Here, we aimed to identify the role of PIWI genes and piRNAs in Culex quinquefasciatus and Culex tarsalis cells during virus infection. We examined the effect of PIWI gene silencing on virus replication of two arboviruses and three insect-specific viruses in Cx. quinquefasciatus derived cells (Hsu) and Cx. tarsalis derived (CT) cells. We show that Piwi4 is antiviral against the La Crosse orthobunyavirus (LACV) in Hsu and CT cells, and the insect-specific rhabdovirus Merida virus (MERDV) in Hsu cells. None of the silenced PIWI genes impacted replication of the two flaviviruses Usutu virus (USUV) and Calbertado virus, or the phasivirus Phasi-Charoen-like virus. We further used small RNA sequencing to determine that LACV-derived piRNAs, but not USUV-derived piRNAs were generated in Hsu cells and that PIWI gene silencing resulted in a small reduction in vpiRNAs. Finally, we determined that LACV-derived DNA was produced in Hsu cells during infection, but whether this viral DNA is required for vpiRNA production remains unclear. Overall, we expanded our knowledge on the piRNA pathway and how it relates to the antiviral response in Culex spp mosquitoes.

The DEAD-box RNA helicase Dhx15 controls glycolysis and arbovirus replication in Aedes aegypti mosquito cells

Aedes aegypti mosquitoes are responsible for the transmission of arthropod-borne (arbo)viruses including dengue and chikungunya virus (CHIKV) but in contrast to human hosts, arbovirus-infected mosquitoes are able to efficiently control virus replication to sub-pathological levels. Yet, our knowledge of the molecular interactions of arboviruses with their mosquito hosts is incomplete. Here, we aimed to identify and characterize novel host genes that control arbovirus replication in Aedes mosquitoes. RNA binding proteins (RBPs) are well-known to regulate immune signaling pathways in all kingdoms of life. We therefore performed a knockdown screen targeting 461 genes encoding predicted RBPs in Aedes aegypti Aag2 cells and identified 15 genes with antiviral activity against Sindbis virus. Amongst these, the three DEAD-box RNA helicases AAEL004419/Dhx15, AAEL008728, and AAEL004859 also acted as antiviral factors in dengue and CHIKV infections. Here, we explored the mechanism of Dhx15 in regulating an antiviral transcriptional response in mosquitoes by silencing Dhx15 in Aag2 cells followed by deep-sequencing of poly-A enriched RNAs. Dhx15 knockdown in uninfected and CHIKV-infected cells resulted in differential expression of 856 and 372 genes, respectively. Interestingly, amongst the consistently downregulated genes, glycolytic process was the most enriched gene ontology (GO) term as the expression of all core enzymes of the glycolytic pathway was reduced, suggesting that Dhx15 regulates glycolytic function. A decrease in lactate production indicated that Dhx15 silencing indeed functionally impaired glycolysis. Modified rates of glycolytic metabolism have been implicated in controlling the replication of several classes of viruses and strikingly, infection of Aag2 cells with CHIKV by itself also resulted in the decrease of several glycolytic genes. Our data suggests that Dhx15 regulates replication of CHIKV, and possibly other arboviruses, by controlling glycolysis in mosquito cells.

Functional Analysis of Individual piRNAs in Aedes aegypti Cells and Embryos Using Antisense Oligonucleotides

In insects, PIWI-interacting (pi)RNAs fulfill versatile regulatory functions inside and outside the germline, including posttranscriptional repression of transposable elements and regulation of gene expression. Canonically, piRNAs act-and have been studied-as a conglomerate of several thousand sequences that cooperatively silence target RNAs. Interestingly, however, an increasing number of studies have demonstrated that individual piRNAs can have profound biological activity as a unique piRNA sequence. Prime examples are the tapiR1 and 2 piRNAs, which mediate target RNA degradation in the developing embryo of Aedes mosquitoes. To study such outstanding individual piRNA species, we describe here a method to interfere with RNA target silencing using antisense oligonucleotides in cell culture as well as in mosquito pre-blastoderm embryos. Although the method has been established for Aedes mosquitoes, it can likely be adapted for use in other invertebrate species as well.

Aedes aegypti Piwi4 Structural Features Are Necessary for RNA Binding and Nuclear Localization

The PIWI-interacting RNA (piRNA) pathway provides an RNA interference (RNAi) mechanism known from Drosophila studies to maintain the integrity of the germline genome by silencing transposable elements (TE). Aedes aegypti mosquitoes, which are the key vectors of several arthropod-borne viruses, exhibit an expanded repertoire of Piwi proteins involved in the piRNA pathway, suggesting functional divergence. Here, we investigate RNA-binding dynamics and subcellular localization of A. aegypti Piwi4 (AePiwi4), a Piwi protein involved in antiviral immunity and embryonic development, to better understand its function. We found that AePiwi4 PAZ (Piwi/Argonaute/Zwille), the domain that binds the 3′ ends of piRNAs, bound to mature (3′ 2′ O-methylated) and unmethylated RNAs with similar micromolar affinities (K_D = 1.7 ± 0.8 μM and K_D of 5.0 ± 2.2 μM, respectively; p = 0.05) in a sequence independent manner. Through site-directed mutagenesis studies, we identified highly conserved residues involved in RNA binding and found that subtle changes in the amino acids flanking the binding pocket across PAZ proteins have significant impacts on binding behaviors, likely by impacting the protein secondary structure. We also analyzed AePiwi4 subcellular localization in mosquito tissues. We found that the protein is both cytoplasmic and nuclear, and we identified an AePiwi4 nuclear localization signal (NLS) in the N-terminal region of the protein. Taken together, these studies provide insights on the dynamic role of AePiwi4 in RNAi and pave the way for future studies aimed at understanding Piwi interactions with diverse RNA populations.

Endogenous piRNA-guided slicing triggers responder and trailer piRNA production from viral RNA in Aedes aegypti mosquitoes

In the germline of animals, PIWI interacting (pi)RNAs protect the genome against the detrimental effects of transposon mobilization. In Drosophila, piRNA-mediated cleavage of transposon RNA triggers the production of responder piRNAs via ping-pong amplification. Responder piRNA 3' end formation by the nuclease Zucchini is coupled to the production of downstream trailer piRNAs, expanding the repertoire of transposon piRNA sequences. In Aedes aegypti mosquitoes, piRNAs are generated from viral RNA, yet, it is unknown how viral piRNA 3' ends are formed and whether viral RNA cleavage gives rise to trailer piRNA production. Here we report that in Ae. aegypti, virus- and transposon-derived piRNAs have sharp 3' ends, and are biased for downstream uridine residues, features reminiscent of Zucchini cleavage of precursor piRNAs in Drosophila. We designed a reporter system to study viral piRNA 3' end formation and found that targeting viral RNA by abundant endogenous piRNAs triggers the production of responder and trailer piRNAs. Using this reporter, we identified the Ae. aegypti orthologs of Zucchini and Nibbler, two nucleases involved in piRNA 3' end formation. Our results furthermore suggest that autonomous piRNA production from viral RNA can be triggered and expanded by an initial cleavage event guided by genome-encoded piRNAs.