Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Analysis of somatic piRNAs in the malaria mosquito Anopheles coluzzii reveals atypical classes of genic small RNAs

Piwi-interacting small RNAs (piRNA) play a key role in controlling the activity of transposable elements (TEs) in the animal germline. In diverse arthropod species, including the pathogen vectors mosquitoes, the piRNA pathway is also active in nongonadal somatic tissues, where its targets and functions are less clear. Here, we studied the features of small RNA production in head and thorax tissues of an uninfected laboratory strain of Anopheles coluzzii focusing on the 24-32-nt-long RNAs. Small RNAs derived from repetitive elements constitute a minor fraction while most small RNAs process from long noncoding RNAs (lncRNAs) and protein-coding gene mRNAs. The majority of small RNAs derived from repetitive elements and lncRNAs exhibited typical piRNAs features. By contrast, majority of protein-coding gene-derived 24-32 nt small RNAs lack the hallmarks of piRNAs and have signatures of nontemplated 3' end tailing. Most of the atypical small RNAs exhibit female-biased expression and originate from mitochondrial and nuclear genes involved in energy metabolism. We also identified atypical genic small RNAs in Anopheles gambiae somatic tissues, which further validates the noncanonical mechanism of their production. We discuss a novel mechanism of small RNA production in mosquito somatic tissues and the possible functional significance of genic small RNAs.

Comprehensive Genomic Analysis Reveals Novel Transposable Element-Derived MicroRNA Regulating Caste Differentiation in Honeybees

The honeybee (Apis mellifera) is a highly social insect whose caste differentiation is regulated by epigenetic mechanisms, representing a classic example of phenotypic plasticity in social insects. Although the importance of transposable elements (TEs) in epigenetic research is well recognized, their specific role in honeybee caste differentiation has not been fully explored. This study reveals a novel regulatory mechanism where the microRNA (miRNA) ame-mir-3721-3p, derived from ApME (Apis miniature inverted-repeat TEs), suppresses DNA methyltransferase gene DNMT3, promoting queen-like development in honeybee larvae. Genome-wide analysis identified 43 ApME elements in Apis, with ApMETm15 being particularly abundant and species-specific. These elements gave rise to 6 miRNAs, including ame-mir-3721-3p which showed notable regulatory potential. Target gene prediction and luciferase reporter assays confirmed that ame-mir-3721-3p binds to and suppresses DNMT3 expression. Spatiotemporal expression analysis indicated that ame-mir-3721-3p is significantly upregulated during the critical L3 larval stage, exhibiting a similar expression pattern to DNMT3. Larval feeding experiments with agomir demonstrated that ame-mir-3721-3p suppresses DNMT3 expression and significantly impacts the expression of genes related to the juvenile hormone and ecdysone pathways. Further physiological evidence showed that when larvae were treated with agomir-3721 during the critical caste differentiation window (L3-L4 stage), the emerging adult bees exhibited increased body size, doubled ovarian area, and significantly higher frequency of ovary development, with significant upregulation of ovarian-specific marker genes. These findings provide direct evidence for ame-mir-3721-3p's role in promoting queen-like developmental trajectories during caste differentiation, uncovering a new regulatory pathway in honeybee development and offering insights into epigenetic mechanisms in social insects.

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.

New class of RNA biomarker for endometriosis diagnosis: The potential of salivary piRNA expression

OBJECTIVES

In contrast to miRNA expression, little attention has been given to piwiRNA (piRNA) expression among endometriosis patients. The aim of the present study was to explore the human piRNAome and to investigate a potential piRNA saliva-based diagnostic signature for endometriosis.

METHODS

Data from the prospective "ENDOmiRNA" study (ClinicalTrials.gov Identifier: NCT04728152) were used. Saliva samples from 200 patients were analyzed in order to evaluate human piRNA expression using the piRNA bank. Next Generation Sequencing (NGS), barcoding of unique molecular identifiers and both Artificial Intelligence (AI) and machine learning (ML) were used. For each piRNA, sensitivity, specificity, and ROC AUC values were calculated for the diagnosis of endometriosis.

RESULTS

201 piRNAs were identified, none had an AUC ≥ 0.70, and only three piRNAs (piR-004153, piR001918, piR-020401) had an AUC between ≥ 0.6 and < 0.70. Seven were differentially expressed: piR-004153, piR-001918, piR-020401, piR-012864, piR-017716, piR-020326 and piR-016904. The respective correlation and accuracy to diagnose endometriosis according to the F1-score, sensitivity, specificity, and AUC ranged from 0 to 0.862 %, 0-0.961 %, 0.085-1, and 0.425-0.618. A correlation was observed between the patients' age (≥35 years) and piR-004153 (p = 0.002) and piR-017716 (p = 0.030). Among the 201 piRNAs, four were differentially expressed in patients with and without hormonal treatment: piR-004153 (p = 0.015), piR-020401 (p = 0.001), piR-012864 (p = 0.036) and piR-017716 (p = 0.009).

CONCLUSION

Our results support the link between piRNAs and endometriosis physiopathology and establish its utility as a potential diagnostic biomarker using saliva samples. Per se, piRNA expression should be analyzed along with the clinical status of a patient.

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.

The role of polyplexes in developing a green sustainable approach in agriculture

Rise in global population has increased the food demands and thus the competition among farmers to produce more and more. In the race to obtain higher productivity, farmers have resorted to injudicious farming practices that include the reckless use of nitrogenous fertilizers and intensive cropping on farmlands. Such practices have paved the path for large scale infestations of crops and plants by pests thus affecting the plant productivity and crop vigour. There are several traditional techniques to control pest infestations in plants such as the use of chemical or bio-pesticides, and integrated pest management practices which face several drawbacks. Delivery of gene/nucleic acid in plants through genetic engineering approaches is a more sustainable and effective method of protection against pests. The technology of RNA interference (RNAi) provides a sustainable solution to counter pest control problems faced by other traditional techniques. The RNAi technique involves delivery of dsDNA/dsRNA or other forms of nucleic acids into target organisms thereby bringing about gene silencing. However, RNAi is also limited to its use because of their susceptibility to degradation wherein the use of cationic polymers can provide a tangible solution. Cationic polymers form stable complexes with the nucleic acids known as "polyplexes", which may be attributed to their high positive charge densities thus protecting the exogenous nucleic acids from extracellular degradation. The current paper focuses on the utility of nucleic acids as a sustainable tool for pest control in crops and the use of cationic polymers for the efficient delivery of nucleic acids in pests thus protecting the plant from infestations.

Virus and endogenous viral element-derived small non-coding RNAs and their roles in insect-virus interaction

RNA interference pathways mediated by different types of small non-coding RNAs (siRNAs, miRNAs and piRNAs) are conserved biological responses to exotic stresses, including viral infection. Aside from the well-established siRNA pathway, the miRNA pathway and the piRNA pathway process viral sequences, exogenously or endogenously, into miRNAs and piRNAs, respectively. During the host-virus interaction, viral sequences, including both coding and non-coding sequences, can be integrated as endogenous viral elements (EVEs) and thereby become present within the germline of a non-viral organism. In recent years, significant progress has been made in characterizing the biogenesis and function of viruses and EVEs associated with snRNAs. Overall, the siRNA pathway acts as the primarily antiviral defense against a wide range of exogenous viruses; the miRNA pathways associated with viruses or EVEs function in antiviral response and host gene regulation; EVE derived piRNAs with a ping-pong signature have the potential to limit cognate viral infection.

Vector Specificity of Arbovirus Transmission

More than 25% of human infectious diseases are vector-borne diseases (VBDs). These diseases, caused by pathogens shared between animals and humans, are a growing threat to global health with more than 2.5 million annual deaths. Mosquitoes and ticks are the main vectors of arboviruses including flaviviruses, which greatly affect humans. However, all tick or mosquito species are not able to transmit all viruses, suggesting important molecular mechanisms regulating viral infection, dissemination, and transmission by vectors. Despite the large distribution of arthropods (mosquitoes and ticks) and arboviruses, only a few pairings of arthropods (family, genus, and population) and viruses (family, genus, and genotype) successfully transmit. Here, we review the factors that might limit pathogen transmission: internal (vector genetics, immune responses, microbiome including insect-specific viruses, and coinfections) and external, either biotic (adult and larvae nutrition) or abiotic (temperature, chemicals, and altitude). This review will demonstrate the dynamic nature and complexity of virus–vector interactions to help in designing appropriate practices in surveillance and prevention to reduce VBD threats.

RNAi-based immunity in insects against baculoviruses and the strategies of baculoviruses involved in siRNA and miRNA pathways to weaken the defense

Protection against viral infection in hosts concerns diverse cellular and molecular mechanisms, among which RNA interference (RNAi) response is a vital one. Small interfering RNAs (siRNAs), microRNAs (miRNAs) and PIWI interacting RNAs (piRNAs) are primary categories of small RNAs involved in RNAi response, playing significant roles in restraining viral invasion. However, during a long-term coevolution, viruses have gained the ability to evade, avoid, or suppress antiviral immunity to ensure efficient replication and transmission. Baculoviruses are enveloped, insect-pathogenic viruses with double-stranded circular DNA genomes, which encode suppressors of siRNA pathway and miRNAs targeting immune-related genes to mask the antiviral activity of their hosts. This review summarized recent findings for the RNAi-based antiviral immunity in insects as well as the strategies that baculoviruses exploit to break the shield of host siRNA pathway, and hijack cellular miRNAs or encode their own miRNAs that regulate both viral and cellular gene expression to create a favorable environment for viral infection.

Argonaute-CLIP delineates versatile, functional RNAi networks in Aedes aegypti, a major vector of human viruses

Argonaute (AGO) proteins bind small RNAs to silence complementary RNA transcripts, and they are central to RNA interference (RNAi). RNAi is critical for regulation of gene expression and antiviral defense in Aedes aegypti mosquitoes, which transmit Zika, chikungunya, dengue, and yellow fever viruses. In mosquitoes, AGO1 mediates miRNA interactions, while AGO2 mediates siRNA interactions. We applied AGO-crosslinking immunoprecipitation (AGO-CLIP) for both AGO1 and AGO2, and we developed a universal software package for CLIP analysis (CLIPflexR), identifying 230 small RNAs and 5,447 small RNA targets that comprise a comprehensive RNAi network map in mosquitoes. RNAi network maps predicted expression levels of small RNA targets in specific tissues. Additionally, this resource identified unexpected, context-dependent AGO2 target preferences, including endogenous viral elements and 3'UTRs. Finally, contrary to current thinking, mosquito AGO2 repressed imperfect targets. These findings expand our understanding of small RNA networks and have broad implications for the study of antiviral RNAi.

piRNAs as Modulators of Disease Pathogenesis

Advances in understanding disease pathogenesis correlates to modifications in gene expression within different tissues and organ systems. In depth knowledge about the dysregulation of gene expression profiles is fundamental to fully uncover mechanisms in disease development and changes in host homeostasis. The body of knowledge surrounding mammalian regulatory elements, specifically regulators of chromatin structure, transcriptional and translational activation, has considerably surged within the past decade. A set of key regulators whose function still needs to be fully elucidated are small non-coding RNAs (sncRNAs). Due to their broad range of unfolding functions in the regulation of gene expression during transcription and translation, sncRNAs are becoming vital to many cellular processes. Within the past decade, a novel class of sncRNAs called PIWI-interacting RNAs (piRNAs) have been implicated in various diseases, and understanding their complete function is of vital importance. Historically, piRNAs have been shown to be indispensable in germline integrity and stem cell development. Accumulating research evidence continue to reveal the many arms of piRNA function. Although piRNA function and biogenesis has been extensively studied in Drosophila, it is thought that they play similar roles in vertebrate species, including humans. Compounding evidence suggests that piRNAs encompass a wider functional range than small interfering RNAs (siRNAs) and microRNAs (miRNAs), which have been studied more in terms of cellular homeostasis and disease. This review aims to summarize contemporary knowledge regarding biogenesis, and homeostatic function of piRNAs and their emerging roles in the development of pathologies related to cardiomyopathies, cancer, and infectious diseases.

Diverse Defenses: A Perspective Comparing Dipteran Piwi-piRNA Pathways

Animals face the dual threat of virus infections hijacking cellular function and transposons proliferating in germline genomes. For insects, the deeply conserved RNA interference (RNAi) pathways and other chromatin regulators provide an important line of defense against both viruses and transposons. For example, this innate immune system displays adaptiveness to new invasions by generating cognate small RNAs for targeting gene silencing measures against the viral and genomic intruders. However, within the Dipteran clade of insects, Drosophilid fruit flies and Culicids mosquitoes have evolved several unique mechanistic aspects of their RNAi defenses to combat invading transposons and viruses, with the Piwi-piRNA arm of the RNAi pathways showing the greatest degree of novel evolution. Whereas central features of Piwi-piRNA pathways are conserved between Drosophilids and Culicids, multiple lineage-specific innovations have arisen that may reflect distinct genome composition differences and specific ecological and physiological features dividing these two branches of Dipterans. This perspective review focuses on the most recent findings illuminating the Piwi/piRNA pathway distinctions between fruit flies and mosquitoes, and raises open questions that need to be addressed in order to ameliorate human diseases caused by pathogenic viruses that mosquitoes transmit as vectors.

Immunopathology and immune homeostasis during viral infection in insects

Organisms clear infections by mounting an immune response that is normally turned off once the pathogens have been cleared. However, sometimes this immune response is not properly or timely arrested, resulting in the host damaging itself. This immune dysregulation may be referred to as immunopathology. While our knowledge of immune and metabolic pathways in insects, particularly in response to viral infections, is growing, little is known about the mechanisms that regulate this immune response and hence little is known about immunopathology in this important and diverse group of organisms. In this chapter we focus both on documenting the molecular mechanisms described involved in restoring immune homeostasis in insects after viral infections and on identifying potential mechanisms for future investigation. We argue that learning about the immunopathological consequences of an improperly regulated immune response in insects will benefit both insect and human health.

Global Analysis of Small Non-Coding RNA Populations across Tissues in the Malaria Vector, Anopheles gambiae

Malaria is a major global health problem, where the anautogenous female mosquito Anopheles gambiae serves as a major vector. In order to combat this devastating disease, understanding mosquito physiology is paramount. Numerous studies in the vector field demonstrate that small non-coding RNAs (ncRNAs) play essential roles in numerous aspects of mosquito physiology. While our previous miRNA annotation work demonstrated expression dynamics across differing tissues, miRNAs represented less than 20% of all small ncRNAs in our small RNA-Seq libraries. To this end, we systematically classified multiple small ncRNA groups across mosquito tissues. Here we (i) determined a new enriched-midgut miRNA, (ii) updated the piRNA annotation in ovaries with a genomic map of unique-mapping piRNAs, (iii) identified pan-tissue and tissue-enriched mRNA-derived small ncRNAs, and (iv) assessed AGO1- and AGO2- loading of candidate small ncRNAs. Continued research will broaden our view of small ncRNAs and greatly aide in our understanding on how these molecules contribute to mosquito physiology.

Conserved Small Nucleotidic Elements at the Origin of Concerted piRNA Biogenesis from Genes and lncRNAs

PIWI-interacting RNAs (piRNAs) target transcripts by sequence complementarity serving as guides for RNA slicing in animal germ cells. The piRNA pathway is increasingly recognized as critical for essential cellular functions such as germline development and reproduction. In the Anopheles gambiae ovary, as much as 11% of piRNAs map to protein-coding genes. Here, we show that ovarian mRNAs and long non-coding RNAs (lncRNAs) are processed into piRNAs that can direct other transcripts into the piRNA biogenesis pathway. Targeting piRNAs fuel transcripts either into the ping-pong cycle of piRNA amplification or into the machinery of phased piRNA biogenesis, thereby creating networks of inter-regulating transcripts. RNAs of the same network share related genomic repeats. These repeats give rise to piRNAs, which target other transcripts and lead to a cascade of concerted RNA slicing. While ping-pong networks are based on repeats of several hundred nucleotides, networks that rely on phased piRNA biogenesis operate through short ~40-nucleotides long repeats, which we named snetDNAs. Interestingly, snetDNAs are recurring in evolution from insects to mammals. Our study brings to light a new type of conserved regulatory pathway, the snetDNA-pathway, by which short sequences can include independent genes and lncRNAs in the same biological pathway.

Improving RNAi efficiency for pest control in crop species

The application of RNAi promotes the development of novel approaches toward plant protection in a sustainable way. Genetically modified crops expressing dsRNA have been developed as commercial products with great potential in insect pest management. Alternatively, some nontransformative approaches, including foliar spray, irrigation and trunk injection, are favorable in actual utilization. In this review, we summarize the recent progress and successful cases of RNAi-based pest management strategy, explore essential implications and possibilities to improve RNAi efficiency by delivery of dsRNA through transformative and nontransformative approaches, and highlight the remaining challenges and important issues related to the application of this technology.

Expression Dynamics of Core RNAi Machinery Genes in Pea Aphids Upon Exposure to Artificially Synthesized dsRNA and miRNAs

The pea aphid is an important pest of vegetables and causes serious losses worldwide. RNA interference (RNAi) is an effective pest control tool, and three sub-pathways have been described: The miRNA pathway, siRNA pathway, and piRNA pathway. A large number of genes in miRNA pathway and piRNA pathway are found to be expanded. To study the roles of these genes, the expression of 25 core RNAi genes was screened in spatiotemporal samples, artificially synthesized dsRNA and miRNA treated samples. The 25 genes were all expressed during different development stages and in different tissues. In dsRNA-treated samples and miRNA-treated samples, the expressions of genes in these three pathways were induced, especially the expanded genes. This suggests a complex network of RNAi core genes in the three sub-pathways. Treatment of miRNA seems to induce gene expression in a dosage-dependent manner. These results increase our knowledge of the siRNA pathway and related factors from RNAi pathway in aphids and promote the use of RNAi for the control of aphid pests.

Mechanisms, Applications, and Challenges of Insect RNA Interference

The RNA interference (RNAi) triggered by short/small interfering RNA (siRNA) was discovered in nematodes and found to function in most living organisms. RNAi has been widely used as a research tool to study gene functions and has shown great potential for the development of novel pest management strategies. RNAi is highly efficient and systemic in coleopterans but highly variable or inefficient in many other insects. Differences in double-stranded RNA (dsRNA) degradation, cellular uptake, inter- and intracellular transports, processing of dsRNA to siRNA, and RNA-induced silencing complex formation influence RNAi efficiency. The basic dsRNA delivery methods include microinjection, feeding, and soaking. To improve dsRNA delivery, various new technologies, including cationic liposome-assisted, nanoparticle-enabled, symbiont-mediated, and plant-mediated deliveries, have been developed. Major challenges to widespread use of RNAi in insect pest management include variable RNAi efficiency among insects, lack of reliable dsRNA delivery methods, off-target and nontarget effects, and potential development of resistance in insect populations.

The Genome of Drosophila innubila Reveals Lineage-Specific Patterns of Selection in Immune Genes

Pathogenic microbes can exert extraordinary evolutionary pressure on their hosts. They can spread rapidly and sicken or even kill their host to promote their own proliferation. Because of this strong selective pressure, immune genes are some of the fastest evolving genes across metazoans, as highlighted in mammals and insects. Drosophila melanogaster serves as a powerful model for studying host/pathogen evolution. While Drosophila melanogaster are frequently exposed to various pathogens, little is known about D. melanogaster's ecology, or if they are representative of other Drosophila species in terms of pathogen pressure. Here, we characterize the genome of Drosophila innubila, a mushroom-feeding species highly diverged from D. melanogaster and investigate the evolution of the immune system. We find substantial differences in the rates of evolution of immune pathways between D. innubila and D. melanogaster. Contrasting what was previously found for D. melanogaster, we find little evidence of rapid evolution of the antiviral RNAi genes and high rates of evolution in the Toll pathway. This suggests that, while immune genes tend to be rapidly evolving in most species, the specific genes that are fastest evolving may depend either on the pathogens faced by the host and/or divergence in the basic architecture of the host's immune system.

Endogenous non-retroviral elements in genomes of Aedes mosquitoes and vector competence

Recent extensive (re)emergences of arthropod-borne viruses (arboviruses) such as chikungunya (CHIKV), zika (ZIKV) and dengue (DENV) viruses highlight the role of the epidemic vectors, Aedes aegypti and Aedes albopictus, in their spreading. Differences of vector competence to arboviruses highlight different virus/vector interactions. While both are highly competent to transmit CHIKV (Alphavirus,Togaviridae), only Ae. albopictus is considered as a secondary vector for DENV (Flavivirus, Flaviviridae). Among other factors such as environmental temperature, mosquito antiviral immunity and microbiota, the presence of non-retroviral integrated RNA virus sequences (NIRVS) in both mosquito genomes may modulate the vector competence. Here we review the current knowledge on these elements, highlighting the mechanisms by which they are produced and endogenized into Aedes genomes. Additionally, we describe their involvement in antiviral immunity as a stimulator of the RNA interference pathways and in some rare cases, as producer of viral-interfering proteins. Finally, we mention NIRVS as a tool for understanding virus/vector co-evolution. The recent discovery of endogenized elements shows that virus/vector interactions are more dynamic than previously thought, and genetic markers such as NIRVS could be one of the potential targets to reduce arbovirus transmission.

PIWI pathway against viruses in insects

Piwi-interacting RNAs (piRNAs) are an animal-specific class of small non-coding RNAs that are generated via a biogenesis pathway distinct from small interfering RNAs (siRNAs) and microRNAs (miRNAs). There are variations in piRNA biogenesis that depend on several factors, such as the cell type (germline or soma), the organism, and the purpose for which they are being produced, such as transposon-targeting, viral-targeting, or gene-derived piRNAs. Interestingly, the genes involved in the PIWI/piRNA pathway are more rapidly evolving compared with other RNA interference (RNAi) genes. In this review, the role of the piRNA pathway in the antiviral response is reviewed based on recent findings in insect models such as Drosophila, mosquitoes, midges and the silkworm, Bombyx mori. We extensively discuss the special features that characterize host-virus piRNA responses with respect to the proteins and the genes involved, the viral piRNAs' sequence characteristics, the target strand orientation biases as well as the viral piRNA target hotspots across the viral genomes. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Genome-wide profiling of piRNAs in the whitefly Bemisia tabaci reveals cluster distribution and association with begomovirus transmission

The whitefly Bemisia tabaci MEAM1 is a notorious vector capable of transmitting many plant viruses, resulting in serious crop loss and food shortage around the world. To investigate potential sRNA-mediated regulatory mechanisms in whiteflies that are affected by virus acquisition and transmission, we conducted small RNA (sRNA) deep sequencing and performed genome-wide profiling of piwi-interacting RNAs (piRNAs) in whiteflies that were fed on tomato yellow leaf curl virus (TYLCV)-infected or non-infected tomato plants for 24, 48, and 72 h. In the present study, piRNA reads ranging from 564,395 to 1,715,652 per library were identified and shown to distribute unevenly in clusters (57 to 96 per library) on the whitefly (B. tabaci MEAM1) genome. Among them, 53 piRNA clusters were common for all treatments. Comparative analysis between libraries generated from viruliferous and non-viruliferous whiteflies identified five TYLCV-induced and 24 TYLCV-suppressed piRNA clusters. Approximately 62% of piRNAs were derived from non-coding sequences including intergenic regions, introns, and untranslated regions (UTRs). The remaining 38% were derived from coding sequences (CDS) or repeat elements. Interestingly, six protein coding genes were targeted by the TYLCV-induced piRNAs. We identified a large number of piRNAs that were distributed in clusters across the whitefly genome, with 60% being derived from non-coding regions. Comparative analysis revealed that feeding on a virus-infected host caused induction and suppression of only a small number of piRNA clusters in whiteflies. Although piRNAs primarily regulate the activity of transposable elements, our results suggest that they may have additional functions in regulating protein coding genes and in insect-virus interactions.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

In silico analysis of endogenous siRNAs associated transposable elements and NATs in Schistosoma japonicum reveals their putative roles during reproductive development

Schistosomiasis is a neglected tropical disease caused by trematode of the genus Schistosoma. Successful reproductive development is critical for the production of eggs, which are responsible for host pathology and disease dissemination. Endogenous small non-coding RNAs play important roles in many biological processes such as protection against foreign pathogens, cell differentiation, and chromosomal stability by regulating target gene expression at the transcriptional and post-transcriptional levels. In this study, we performed in silico analysis of endogenous small non-coding RNAs in different stages, and sex of S. japonicum focusing on endogenous small interfering RNAs (endo-siRNAs) generated from transposable elements (TEs) and natural antisense transcripts (NATs). Both total and unique siRNA populations show 18-30 nt in length, but the predominant size was 20 nt and the leading first base was adenosine. Sense TE-derived endo-siRNAs reads were higher than antisense reads at different relative positions of TEs, whereas no such difference was observed for NAT-derived endo-siRNAs. TE- and NAT-derived endo-siRNAs were more enriched in the male compared to female worms, with the higher relative expression in early phase of pairing. Putative targets of endo-siRNAs indicated more of them in males (106 and 66) than in females (6 and 23) for TE- and NAT-derived endo-siRNAs, respectively. Our preliminary study revealed vital role of endo-siRNAs during the reproductive development of S. japonicum and provide clues for putative novel targets to suppress worm reproduction and direction for effective anti-schistosomal drug development.

Contrasting Sex-and Caste-Dependent piRNA Profiles in the Transposon Depleted Haplodiploid Honeybee Apis mellifera

Protecting genome integrity against transposable elements is achieved by intricate molecular mechanisms involving PIWI proteins, their associated small RNAs (piRNAs), and epigenetic modifiers such as DNA methylation. Eusocial bees, in particular the Western honeybee, Apis mellifera, have one of the lowest contents of transposable elements in the animal kingdom, and, unlike other animals with a functional DNA methylation system, appear not to methylate their transposons. This raises the question of whether the PIWI machinery has been retained in this species. Using comparative genomics, mass spectrometry, and expressional profiling, we present seminal evidence that the piRNA system is conserved in honeybees. We show that honey bee piRNAs contain a 2'-O-methyl modification at the 3' end, and have a bias towards a 5' terminal U, which are signature features of their biogenesis. Both piRNA repertoire and expression levels are greater in reproductive individuals than in sterile workers. Haploid males, where the detrimental effects of transposons are dominant, have the greatest piRNA levels, but surprisingly, the highest expression of transposons. These results show that even in a transposon-depleted species, the piRNA system is required to guard the vulnerable haploid genome and reproductive castes against transposon-associated genomic instability. This also suggests that dosage plays an important role in the regulation of transposons and piRNAs expression in haplo-diploid systems.

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.

Association of microRNAs with Argonaute proteins in the malaria mosquito Anopheles gambiae after blood ingestion

Drastic changes in gene expression occur after adult female mosquitoes take a blood meal and use the nutrients for egg maturation. A growing body of evidence indicates that microRNAs (miRNAs) contribute to this tightly controlled tissue- and stage-specific gene expression. To investigate the role of miRNAs, we monitored miRNA expression in the mosquito Anopheles gambiae during the 72-h period immediately after blood feeding. We also measured the association of miRNAs with Argonaute 1 (Ago1) and Argonaute 2 (Ago2) to assess the functional status of individual miRNA species. Overall, 173 mature miRNAs were precipitated with Ago1 and Ago2, including 12 new miRNAs, the orthologs of which are found thus far only in other Anopheles species. Ago1 is the predominant carrier of miRNAs in Anopheles gambiae. The abundance and Ago loading of most of the mature miRNAs were relatively stable after blood ingestion. However, miRNAs of the miR-309/286/2944 cluster were considerably upregulated after blood feeding. Injection of the specific antagomir for miR-309 resulted in smaller developing oocytes and ultimately fewer eggs. In addition, the Ago association of some miRNAs was not proportional to their cellular abundance, suggesting that integration of miRNAs into the Ago complexes is regulated by additional mechanisms.

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.

Small RNA Regulators of Plant-Hemipteran Interactions: Micromanagers with Versatile Roles

Non-coding small RNAs (sRNAs) in plants have important roles in regulating biological processes, including development, reproduction, and stress responses. Recent research indicates significant roles for sRNA-mediated gene silencing during plant-hemipteran interactions that involve all three of these biological processes. Plant responses to hemipteran feeding are determined by changes in the host transcriptome that appear to be fine-tuned by sRNAs. The role of sRNA in plant defense responses is complex. Different forms of sRNAs, with specific modes of action, regulate changes in the host transcriptome primarily through post-transcriptional gene silencing and occasionally through translational repression. Plant genetic resistance against hemipterans provides a model to explore the regulatory roles of sRNAs in plant defense. Aphid-induced sRNA expression in resistance genotypes delivers a new paradigm in understanding the regulation of R gene-mediated resistance in host plants. Unique sRNA profiles, including changes in sRNA biogenesis and expression can also provide insights into susceptibility to insect herbivores. Activation of phytohormone-mediated defense responses against insect herbivory is another hallmark of this interaction, and recent studies have shown that regulation of phytohormone signaling is under the control of sRNAs. Hemipterans feeding on resistant plants also show changes in insect sRNA profiles, possibly influencing insect development and reproduction. Changes in insect traits such as fecundity, host range, and resistance to insecticides are impacted by sRNAs and can directly contribute to the success of certain insect biotypes. In addition to causing direct damage to the host plant, hemipteran insects are often vectors of viral pathogens. Insect anti-viral RNAi machinery is activated to limit virus accumulation, suggesting a role in insect immunity. Virus-derived long sRNAs strongly resemble insect piRNAs, leading to the speculation that the piRNA pathway is induced in response to viral infection. Evidence for robust insect RNAi machinery in several hemipteran species is of immense interest and is being actively pursued as a possible tool for insect control. RNAi-induced gene silencing following uptake of exogenous dsRNA was successfully demonstrated in several hemipterans and the presence of sid-1 like genes support the concept of a systemic response in some species.

Increased production of piRNAs from euchromatic clusters and genes in Anopheles gambiae compared with Drosophila melanogaster

BACKGROUND

Specific genomic loci, termed Piwi-interacting RNA (piRNA) clusters, manufacture piRNAs that serve as guides for the inactivation of complementary transposable elements (TEs). The piRNA pathway has been accurately detailed in Drosophila melanogaster, while it remains poorly examined in other insects. This pathway is increasingly recognized as critical for germline development and reproduction. Understanding of the piRNA functions in mosquitoes could offer an opportunity for disease vector control by the reduction of their reproductive potential.

RESULTS

To analyze the similarities and differences in this pathway between Drosophila and mosquito, we performed an in-depth analysis of the genomic loci producing piRNAs and their targets in the African malaria vector Anopheles gambiae. We identified 187 piRNA clusters in the An. gambiae genome and 155 piRNA clusters in the D. melanogaster genome. We demonstrate that many more piRNA clusters in the mosquito compared with the fruit fly are uni-directionally transcribed and are located outside pericentromeric heterochromatin. About 11 % of the An. gambiae piRNA population map to gene transcripts. This is a noticeable increase compared with the ~6 % of the piRNA population mapped to genes in D. melanogaster. A subset of the piRNA-enriched genes in An. gambiae has functions related to reproduction and development. At least 24 and 65 % of the mapped piRNAs correspond to genomic TE sequences in An. gambiae and D. melanogaster, respectively. DNA transposons and non-LTR retrotransposons are more abundant in An. gambiae, while LTR retrotransposons are more abundant in D. melanogaster. Yet, piRNAs predominantly target LTR retrotransposons in both species, which may point to a distinct feature of these elements compared to the other classes of TEs concerning their silencing by the piRNA pathway.

CONCLUSIONS

Here, we demonstrate that piRNA-producing loci have more ubiquitous distribution in the An. gambiae genome than in the genome of D. melanogaster. Also, protein-coding genes have an increased role in production of piRNAs in the germline of this mosquito. Genes involved in germline and embryonic development of An. gambiae generate a substantial portion of piRNAs, suggesting a role of the piRNA pathway in the epigenetic regulation of the reproductive processes in the African malaria vector.