piRNAs: from biogenesis to function

Epigenetic regulation of epithelial to mesenchymal transition: a trophoblast perspective

Epigenetic changes alter expression of genes at both pre- and post-transcriptional levels without changing their DNA sequence. Accumulating evidence suggests that such changes can modify cellular behaviour and characteristics required during development and in response to various extracellular stimuli. Trophoblast cells develop from the outermost trophectoderm layer of the blastocyst and undergo many phenotypic changes as the placenta develops. One such phenotypic change is differentiation of the epithelial natured cytotrophoblasts into the mesenchymal natured extravillous trophoblasts. The extravillous trophoblasts are primarily responsible for invading into the maternal decidua and thus establishing connection with the maternal spiral arteries. Any dysregulation of this process can have adverse effects on the pregnancy outcome. Hence, tight regulation of this epithelial-mesenchymal transition is critical for successful pregnancy. This review summarizes the recent research on the epigenetic regulation of the epithelial-mesenchymal transition occurring in the trophoblast cells during placental development. The functional significance of chemical modifications of DNA and histone, which regulate transcription, as well as non-coding RNAs, which control gene expression post-transcriptionally, is discussed in relation to trophoblast biology.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera

Background

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development.

Results

Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination.

Conclusions

We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08478-9.

The function of novel small non-coding RNAs (piRNAs, tRFs) and PIWI protein in colorectal cancer

Although great research has been done to clarify the pathogenesis of colorectal cancer (CRC), it is still the third common cancer worldwide. Pathogenesis of CRC as a heterogeneous disease is correlated with mutations and epigenetic alterations that result in the inactivation of tumor-suppressive and activation of an oncogene. Small non-coding RNAs (sncRNAs), emerging as a key player in regulating the genes and protein expression, with a length less than 200 nucleotide (nt). In this review, we aimed to focus on the role and the biogenesis of PIWI-interacting RNA (piRNAs), and tRNA-derived small RNA (tRFs) and PIWI proteins in the initiation, progression, and metastasis of CRC and their molecular mechanisms to understand their function in cancers and to provide better therapeutic strategies for CRC.

The Small RNA Universe of Capitella teleta

RNAi is an evolutionarily fluid mechanism with dramatically different activities across animal phyla. One major group where there has been little investigation is annelid worms. Here, the small RNAs of the polychaete developmental model Capitella teleta are profiled across development. As is seen with nearly all animals, nearly 200 microRNAs were found with 58 high-confidence novel species. Greater miRNA diversity was associated with later stages consistent with differentiation of tissues. Outside miRNA, a distinct composition of other small RNA pathways was found. Unlike many invertebrates, an endogenous siRNA pathway was not observed, indicating pathway loss relative to basal planarians. No processively generated siRNA-class RNAs could be found arising from dsRNA precursors. This has a significant impact on RNAi technology development for this group of animals. Unlike the apparent absence of siRNAs, a significant population of piRNAs was observed. For many piRNAs, phasing and ping-pong biogenesis pathways were identified. Interestingly, piRNAs were found to be highly expressed during early development, suggesting a potential role in regulation in metamorphosis. Critically, the configuration of RNAi factors in C. teleta is found in other annelids and mollusks, suggesting that similar biology is likely to be present in the wider clade. This study is the first in providing comprehensive analysis of small RNAs in annelids.

Dysregulation of Human Somatic piRNA Expression in Parkinson's Disease Subtypes and Stages

Piwi interacting RNAs (piRNAs) are small non-coding single-stranded RNA species 20–31 nucleotides in size generated from distinct loci. In germline tissues, piRNAs are amplified via a “ping-pong cycle” to produce secondary piRNAs, which act in transposon silencing. In contrast, the role of somatic-derived piRNAs remains obscure. Here, we investigated the identity and distribution of piRNAs in human somatic tissues to determine their function and potential role in Parkinson’s disease (PD). Human datasets were curated from the Gene Expression Omnibus (GEO) database and a workflow was developed to identify piRNAs, which revealed 902 somatic piRNAs of which 527 were expressed in the brain. These were mainly derived from chromosomes 1, 11, and 19 compared to the germline tissues, which were from 15 and 19. Approximately 20% of somatic piRNAs mapped to transposon 3′ untranslated regions (UTRs), but a large proportion were sensed to the transcript in contrast to germline piRNAs. Gene set enrichment analysis suggested that somatic piRNAs function in neurodegenerative disease. piRNAs undergo dysregulation in different PD subtypes (PD and Parkinson’s disease dementia (PDD)) and stages (premotor and motor). piR-has-92056, piR-hsa-150797, piR-hsa-347751, piR-hsa-1909905, piR-hsa-2476630, and piR-hsa-2834636 from blood small extracellular vesicles were identified as novel biomarkers for PD diagnosis using a sparse partial least square discriminant analysis (sPLS-DA) (accuracy: 92%, AUC = 0.89). This study highlights a role for piRNAs in PD and provides tools for novel biomarker development.

Identification and Characterization of Polymorphisms in piRNA Regions

piRNAs are a class of noncoding RNAs that perform functions in epigenetic regulation and silencing of transposable elements, a mechanism conserved among most mammals. At present, there are more than 30,000 known piRNAs in humans, of which more than 80% are derived from intergenic regions, and approximately 20% are derived from the introns and exons of pre-mRNAs. It was observed that the expression of the piRNA profile is specific in several organs, suggesting that they play functional roles in different tissues. In addition, some studies suggest that changes in regions that encode piRNAs may have an impact on their function. To evaluate the conservation of these regions and explore the existence of a seed region, SNP and INDEL variant rates were investigated in several genomic regions and compared to piRNA region variant rates. Thus, data analysis, data collection, cleaning, treatment, and exploration were implemented using the R programming language with the help of the RStudio platform. We found that piRNA regions are highly conserved after considering INDELs and do not seem to present an identifiable seed region after considering SNPs and INDEL variants. These findings may contribute to future studies attempting to determine how polymorphisms in piRNA regions can impact diseases.

Structural insights into piRNA biogenesis

Discovered two decades ago, Piwi-interacting RNAs (piRNAs) play critical roles in gene regulation, transposon element repression, and antiviral defense. Dysregulation of piRNAs has been noted in diverse human diseases including cancers. Recently, extensive studies have revealed that many more proteins are involved in piRNA biogenesis. This review will summarize the recent progress in piRNA biogenesis and functions, especially the molecular mechanisms by which piRNA biogenesis-related proteins contribute to piRNA processing.

Investigating deep feedforward neural networks for classification of transposon-derived piRNAs

PIWI-interacting RNAs (piRNAS) form an important class of non-coding RNAs that play a key role in gene expression regulation and genome integrity by silencing transposable elements. However, despite the importance of piRNAs and the large application of deep learning in computational biology, there are few studies of deep learning for piRNAs prediction. Still, current methods focus on using advanced architectures like CNN and variations. This paper presents an investigation on deep feedforward network models for classification of human transposon-derived piRNAs. We developed a lightweight predictor (when compared to other deep learning methods) and we show by practical evidence that simple neural networks can perform as well as better than complex neural networks when using the appropriate hyperparameters. For that, we train, analyze and compare the results of a multilayer perceptron with different hyperparameter choices, such as numbers of hidden layers, activation functions and optimizers, clarifying the advantages and disadvantages of each choice. Our proposed predictor reached a F-score of 0.872, outperforming other state-of-the-art methods for human transposon-derived piRNAs classification. In addition, to better access the generalization of our proposal, we also showed it achieved competitive results when classifying piRNAs of other species.

piRNAs and PIWI Proteins as Diagnostic and Prognostic Markers of Genitourinary Cancers

piRNAs (PIWI-interacting RNAs) are small non-coding RNAs capable of regulation of transposon and gene expression. piRNAs utilise multiple mechanisms to affect gene expression, which makes them potentially more powerful regulators than microRNAs. The mechanisms by which piRNAs regulate transposon and gene expression include DNA methylation, histone modifications, and mRNA degradation. Genitourinary cancers (GC) are a large group of neoplasms that differ by their incidence, clinical course, biology, and prognosis for patients. Regardless of the GC type, metastatic disease remains a key therapeutic challenge, largely affecting patients’ survival rates. Recent studies indicate that piRNAs could serve as potentially useful biomarkers allowing for early cancer detection and therapeutic interventions at the stage of non-advanced tumour, improving patient’s outcomes. Furthermore, studies in prostate cancer show that piRNAs contribute to cancer progression by affecting key oncogenic pathways such as PI3K/AKT. Here, we discuss recent findings on biogenesis, mechanisms of action and the role of piRNAs and the associated PIWI proteins in GC. We also present tools that may be useful for studies on the functioning of piRNAs in cancers.

Somatic piRNAs and Transposons are Differentially Expressed Coincident with Skeletal Muscle Atrophy and Programmed Cell Death

PIWI-interacting RNAs (piRNAs) are small single-stranded RNAs that can repress transposon expression via epigenetic silencing and transcript degradation. They have been identified predominantly in the ovary and testis, where they serve essential roles in transposon silencing in order to protect the integrity of the genome in the germline. The potential expression of piRNAs in somatic cells has been controversial. In the present study we demonstrate the expression of piRNAs derived from both genic and transposon RNAs in the intersegmental muscles (ISMs) from the tobacco hawkmoth Manduca sexta. These piRNAs are abundantly expressed, ∼27 nt long, map antisense to transposons, are oxidation resistant, exhibit a 5’ uridine bias, and amplify via the canonical ping-pong pathway. An RNA-seq analysis demonstrated that 19 piRNA pathway genes are expressed in the ISMs and are developmentally regulated. The abundance of piRNAs does not change when the muscles initiate developmentally-regulated atrophy, but are repressed coincident with the commitment of the muscles undergo programmed cell death at the end of metamorphosis. This change in piRNA expression is correlated with the repression of several retrotransposons and the induction of specific DNA transposons. The developmentally-regulated changes in the expression of piRNAs, piRNA pathway genes, and transposons are all regulated by 20-hydroxyecdysone, the steroid hormone that controls the timing of ISM death. Taken together, these data provide compelling evidence for the existence of piRNA in somatic tissues and suggest that they may play roles in developmental processes such as programmed cell death.

Comparative analysis of piRNA sequences, targets and functions in nematodes

Piwi proteins and Piwi-interacting RNAs (piRNAs) are best known for their roles in suppressing transposons and promoting fertility. Yet piRNA biogenesis and its mechanisms of action differ widely between distantly related species. To better understand the evolution of piRNAs, we characterized the piRNA pathway in C. briggsae, a sibling species of the model organism C. elegans. Our analyses define 25,883 piRNA producing-loci in C. briggsae. piRNA sequences in C. briggsae are extremely divergent from their counterparts in C. elegans, yet both species adopt similar genomic organization that drive piRNA expression. By examining production of Piwi-mediated secondary small RNAs, we identified a set of protein-coding genes that are evolutionarily conserved piRNA targets. In contrast to C. elegans, small RNAs targeting ribosomal RNAs or histone transcripts are not hyper-accumulated in C. briggsae Piwi mutants. Instead, we found that transcripts with few introns are prone to small RNA overamplification. Together our work highlights evolutionary conservation and divergence of the nematode piRNA pathway and provides insights into its role in endogenous gene regulation.

Perinatal Exposure to Nicotine Alters Sperm RNA Profiles in Rats

Perinatal exposure to smoking has been associated with childhood asthma, one of the most common pediatric conditions affecting millions of children globally. Of great interest, this disease phenotype appears heritable as it can persist across multiple generations even in the absence of persistent exposure to smoking in subsequent generations. Although the molecular mechanisms underlying childhood asthma induced by perinatal exposure to smoking or nicotine remain elusive, an epigenetic mechanism has been proposed, which is supported by the data from our earlier analyses on germline DNA methylation (5mC) and histone marks (H3 and H4 acetylation). To further investigate the potential epigenetic inheritance of childhood asthma induced by perinatal nicotine exposure, we profiled both large and small RNAs in the sperm of F1 male rats. Our data revealed that perinatal exposure to nicotine leads to alterations in the profiles of sperm-borne RNAs, including mRNAs and small RNAs, and that rosiglitazone, a PPARγ agonist, can attenuate the effect of nicotine and reverse the sperm-borne RNA profiles of F1 male rats to close to placebo control levels.

piRNA: A promising biomarker in early detection of gastrointestinal cancer

Gastrointestinal (GI) cancer has become the primary concern of today's society due to its aggressive nature and poor prognosis. PIWI-interacting RNAs (piRNAs), a subgroup of non-coding RNAs, are mainly expressed in the germline and have emerged as a critical regulator in gene expression and the epigenetic silencing of DNA transposable elements by interacting with PIWI proteins. piRNAs' dysregulations were reported to promote or suppress the initiation and development of different malignancies, especially gastrointestinal cancers. Recently, several studies suggested the use of piRNAs as potential cancer biomarkers associated with the progression and chemoresistance of GI cancer. Hence, this review article aims to focus on the role of piRNAs in GI cancer progression, metastasis, and their molecular mechanisms as therapeutic markers for GI cancer patients.

piRNAQuest V.2: an updated resource for searching through the piRNAome of multiple species

PIWI interacting RNAs (piRNAs) have emerged as important gene regulators in recent times. Since the release of our first version of piRNAQuest in 2014, lots of novel piRNAs have been annotated in different species other than human, mouse and rat. Such new developments in piRNA research have led us to develop an updated database piRNAQuest V.2. It consists of 92,77,689 piRNA entries for 25 new species of different phylum along with human, mouse and rat. Besides providing primary piRNA features which include their genomic location, with further information on piRNAs overlapping with repeat elements, pseudogenes and syntenic regions, etc., the novel features of this version includes (i) density based cluster prediction, (ii) piRNA expression profile across various healthy and disease systems and (iii) piRNA target prediction. The concept of density-based piRNA cluster identification is robust as it does not consider parametric distribution in its model. The piRNA expression profile for 21 disease systems including cancer have been hosted in addition to 32 tissue specific piRNA expression profile for various species. Further, the piRNA target prediction section includes both predicted and curated piRNA targets within eight disease systems and developmental stages of mouse testis. Further, users can visualize the piRNA-target duplex structure and the ping-pong signature pattern for all the ping-pong piRNA partners in different species. Overall, piRNAQuest V.2 is an updated user-friendly database which will serve as a useful resource to survey, search and retrieve information on piRNAs for multiple species. This freely accessible database is available at http://dibresources.jcbose.ac.in/zhumur/pirnaquest2.

Non-coding RNAs and chromatin: key epigenetic factors from spermatogenesis to transgenerational inheritance

Cellular fate and gene expression patterns are modulated by different epigenetic factors including non-coding RNAs (ncRNAs) and chromatin organization. Both factors are dynamic throughout male germ cell differentiation on the seminiferous tubule, despite the transcriptional inactivation in the last stages of spermatogenesis. Sperm maturation during the caput-to-cauda transit on the epididymis involves changes in chromatin organization and the soma-to-germ line transference of ncRNAs that are essential to obtain a functional sperm for fertilization and embryo development. Here, the male environment (diseases, drugs, mental stress) is crucial to modulate these epigenetic factors throughout sperm maturation, affecting the corresponding offspring. Paternal transgenerational inheritance has been directly related to sperm epigenetic changes, most of them associated with variations in the ncRNA content and chromatin marks. Our aim is to give an overview about how epigenetics, focused on ncRNAs and chromatin, is pivotal to understand spermatogenesis and sperm maturation, and how the male environment impacts the sperm epigenome modulating the offspring gene expression pattern.

PIWI interacting RNAs perspectives: a new avenues in future cancer investigations

As a currently identified small non-coding RNAs (ncRNAs) category, the PIWI-interacting RNAs (piRNAs) are crucial mediators of cell biology. The human genome comprises over 30.000 piRNA genes. Although considered a new field in cancer research, the piRNA pathway is shown by the existing evidence as an active pathway in a variety of different types of cancers with critical impacts on main aspects of cancer progression. Among the regulatory molecules that contribute to maintaining the dynamics of cancer cells, the P-element Induced WImpy testis (PIWI) proteins and piRNAs, as new players, have not been broadly studied so far. Therefore, the identification of cancer-related piRNAs and the assessment of target genes of piRNAs may lead to better cancer prevention and therapy strategies. This review articleaimed to highlight the role and function of piRNAs based on existing data. Understanding the role of piRNA in cancer may provide perspectives on their applications as particular biomarker signature in diagnosis in early stage, prognosis and therapeutic strategies.

Sperm-borne sncRNAs: potential biomarkers for semen fertility?

Semen infertility or sub-fertility, whether in humans or livestock species, remains a major concern for clinicians and technicians involved in reproduction. Indeed, they can cause tragedies in human relationships or have a dramatic overall negative impact on the sustainability of livestock breeding. Understanding and predicting semen fertility issues is therefore crucial and quality control procedures as well as biomarkers have been proposed to ensure sperm fertility. However, their predictive values appeared to be too limited and additional relevant biomarkers are still required to diagnose sub-fertility efficiently. During the last decade, the study of molecular mechanisms involved in spermatogenesis and sperm maturation highlighted the regulatory role of a variety of small non-coding RNAs (sncRNAs) and led to the discovery that sperm sncRNAs comprise both remnants from spermatogenesis and post-testicular sncRNAs acquired through interactions with extracellular vesicles along epididymis. This has led to the hypothesis that sncRNAs may be a source of relevant biomarkers, associated either with sperm functionality or embryo development. This review aims at providing a synthetic overview of the current state of knowledge regarding implication of sncRNA in spermatogenesis defects and their putative roles in sperm maturation and embryo development, as well as exploring their use as fertility biomarkers.

Epigenetic Regulation during Primordial Germ Cell Development and Differentiation

Germline development varies significantly across metazoans. However, mammalian primordial germ cell (PGC) development has key conserved landmarks, including a critical period of epigenetic reprogramming that precedes sex-specific differentiation and gametogenesis. Epigenetic alterations in the germline are of unique importance due to their potential to impact the next generation. Therefore, regulation of, and by, the non-coding genome is of utmost importance during these epigenomic events. Here, we detail the key chromatin changes that occur during mammalian PGC development and how these interact with the expression of non-coding RNAs alongside broader epitranscriptomic changes. We identify gaps in our current knowledge, in particular regarding epigenetic regulation in the human germline, and we highlight important areas of future research.

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.

Cellular abundance shapes function in piRNA-guided genome defense

Defense against genome invaders universally relies on RNA-guided immunity. Prokaryotic CRISPR-Cas and eukaryotic RNA interference pathways recognize targets by complementary base-pairing, which places the sequences of their guide RNAs at the center of self/nonself discrimination. Here, we explore the sequence space of PIWI-interacting RNAs (piRNAs), the genome defense of animals, and establish functional priority among individual sequences. Our results reveal that only the topmost abundant piRNAs are commonly present in every cell, whereas rare sequences generate cell-to-cell diversity in flies and mice. We identify a skewed distribution of sequence abundance as a hallmark of piRNA populations and show that quantitative differences of more than a 1000-fold are established by conserved mechanisms of biogenesis. Finally, our genomics analyses and direct reporter assays reveal that abundance determines function in piRNA-guided genome defense. Taken together, we identify an effective sequence space and untangle two classes of piRNAs that differ in complexity and function. The first class represents the topmost abundant sequences and drives silencing of genomic parasites. The second class sparsely covers an enormous sequence space. These rare piRNAs cannot function in every cell, every individual, or every generation but create diversity with potential for adaptation in the ongoing arms race with genome invaders.

The neurobiology of non-coding RNAs and Alzheimer's disease pathogenesis: Pathways, mechanisms and translational opportunities

In the past two decades, advances in sequencing technology and analysis of the human and mouse genome have led to the discovery of many non-protein-coding RNAs (ncRNAs) including: microRNA, small-interfering RNAs, piwi-associated small RNAs, transfer RNA-derived small RNAs, long-non-coding RNAs and circular RNAs. Compared with healthy controls, levels of some ncRNAs are significantly altered in the central nervous system and blood of patients affected by neurodegenerative disorders like Alzheimer's disease (AD). Although the mechanisms are still not fully elucidated, studies have revealed that these highly conserved ncRNAs are important modulators of gene expression, amyloid-β production, tau phosphorylation, inflammation, synaptic plasticity and neuronal survival, all features considered central to AD pathogenesis. Despite considerable difficulties due to their large heterogeneity, and the complexity of their regulatory pathways, research in this rapidly growing field suggests that ncRNAs hold great potential as biomarkers and therapeutic targets against AD. Herein, we summarize the current knowledge regarding the neurobiology of ncRNA in the context of AD pathophysiology.

Two Piwis with Ago-like functions silence somatic genes at the chromatin level

Most sRNA biogenesis mechanisms involve either RNAse III cleavage or ping-pong amplification by different Piwi proteins harbouring slicer activity. Here, we follow the question why the mechanism of transgene-induced silencing in the ciliate Paramecium needs both Dicer activity and two Ptiwi proteins. This pathway involves primary siRNAs produced from non-translatable transgenes and secondary siRNAs from targeted endogenous loci. Our data does not indicate any signatures from ping-pong amplification but Dicer cleavage of long dsRNA. Ptiwi13 and 14 prefer different sub-cellular localizations and different preferences for primary and secondary siRNAs but do not load them mutually exclusive. Both Piwis enrich for antisense RNAs and show a general preference for uridine-rich sRNAs along the entire sRNA length. In addition, Ptiwi14-loaded siRNAs show a 5´-U signature. Our data indicates both Ptiwis and 2´-O-methylation contributing to strand selection of Dicer cleaved siRNAs. This unexpected function of the two distinct vegetative Piwis extends the increasing knowledge of the diversity of Piwi functions in diverse silencing pathways. We describe an unusual mode of action of Piwi proteins extending not only the great variety of Piwi-associated RNAi pathways but moreover raising the question whether this could have been the primordial one.

Identifying piRNA targets on mRNAs in C. elegans using a deep multi-head attention network

Background

Piwi-interacting RNAs (piRNAs) are the small non-coding RNAs (ncRNAs) that silence genomic transposable elements. And researchers found out that piRNA also regulates various endogenous transcripts. However, there is no systematic understanding of the piRNA binding patterns and how piRNA targets genes. While various prediction methods have been developed for other similar ncRNAs (e.g., miRNAs), piRNA holds distinctive characteristics and requires its own computational model for binding target prediction.

Results

Recently, transcriptome-wide piRNA binding events in C. elegans were probed by PRG-1 CLASH experiments. Based on the probed piRNA-messenger RNAs (mRNAs) binding pairs, in this research, we devised the first deep learning architecture based on multi-head attention to computationally identify piRNA targeting mRNA sites. In the devised deep network, the given piRNA and mRNA segment sequences are first one-hot encoded and undergo a combined operation of convolution and squeezing-extraction to unravel motif patterns. And we incorporate a novel multi-head attention sub-network to extract the hidden piRNA binding rules that can simulate the biological piRNA target recognition process. Finally, the true piRNA–mRNA binding pairs are identified by a deep fully connected sub-network. Our model obtains a supreme discriminatory power of AUC \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=$$\end{document}= 93.3% on an independent test set and successfully extracts the verified binding pattern of a synthetic piRNA. These results demonstrated that the devised model achieves high prediction performance and suggests testable potential biological piRNA binding rules.

Conclusions

In this research, we developed the first deep learning method to identify piRNA targeting sites on C. elegans mRNAs. And the developed deep learning method is demonstrated to be of high accuracy and can provide biological insights into piRNA–mRNA binding patterns. The piRNA binding target identification network can be downloaded from http://cosbi2.ee.ncku.edu.tw/data_download/piRNA_mRNA_binding.

A Review of Discovery Profiling of PIWI-Interacting RNAs and Their Diverse Functions in Metazoans

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs (sncRNAs) that perform crucial biological functions in metazoans and defend against transposable elements (TEs) in germ lines. Recently, ubiquitously expressed piRNAs were discovered in soma and germ lines using small RNA sequencing (sRNA-seq) in humans and animals, providing new insights into the diverse functions of piRNAs. However, the role of piRNAs has not yet been fully elucidated, and sRNA-seq studies continue to reveal different piRNA activities in the genome. In this review, we summarize a set of simplified processes for piRNA analysis in order to provide a useful guide for researchers to perform piRNA research suitable for their study objectives. These processes can help expand the functional research on piRNAs from previously reported sRNA-seq results in metazoans. Ubiquitously expressed piRNAs have been discovered in the soma and germ lines in Annelida, Cnidaria, Echinodermata, Crustacea, Arthropoda, and Mollusca, but they are limited to germ lines in Chordata. The roles of piRNAs in TE silencing, gene expression regulation, epigenetic regulation, embryonic development, immune response, and associated diseases will continue to be discovered via sRNA-seq.

Molecular basis for PICS-mediated piRNA biogenesis and cell division

By incorporating two mutually exclusive factors, PID-1 and TOST-1, C. elegans PICS complex plays important roles in piRNA biogenesis, chromosome segregation and cell division. We firstly map the interaction network between PICS subunits, then uncover the mechanisms underlying the interactions between PICS subunits by solving several complex structures, including those of TOFU-6/PICS-1, ERH-2/PICS-1, and ERH-2/TOST-1. Our biochemical experiment also demonstrates that PICS exists as an octamer consisting of two copies of each subunit. Combining structural analyses with mutagenesis experiments, we identify interfacial residues of PICS subunits that are critical for maintaining intact PICS complex in vitro. Furthermore, using genetics, cell biology and imaging experiments, we find that those mutants impairing the in vitro interaction network within PICS, also lead to dysfunction of PICS in vivo, including mislocalization of PICS, and reduced levels of piRNAs or aberrant chromosome segregation and cell division. Therefore, our work provides structural insights into understanding the PICS-mediated piRNA biogenesis and cell division.

PiRNA pathway in the cardiovascular system: a novel regulator of cardiac differentiation, repair and regeneration

Piwi-interacting RNAs (piRNAs) are a novel group of small non-coding RNA molecules with lengths of 21-35 nucleotides, first identified from the germline. PiRNAs and their associated PIWI clade Argonaute proteins constitute a key part of the piRNA pathway, with the best-known biological function to silence transposable elements in germ cells. The piRNA pathway, in fact, is not exclusive to the germline. Somatic functions of piRNAs have been recorded since their first discovery. To date, involvement of the piRNA pathway has been identified within the biological functions of genome rearrangement, epigenetic regulation, protein regulation in the germline and/or the soma transcriptionally or post-transcriptionally. Emerging evidence has shown that the piRNA pathway is essential for the normal function of the cardiovascular system and that its abnormal expression is correlated with cardiovascular dysfunction, although comprehensive roles of the piRNA pathway in the cardiovascular system and underlying mechanisms remain unclear. In this review, we discuss current findings of piRNA pathway expression in cardiac cell types and their potential functions in cardiac differentiation, repair and regeneration, thus providing new insights into cardiovascular disease development associated with the piRNA pathway.

Epigenetic transgenerational inheritance, gametogenesis and germline development†

One of the most important developing cell types in any biological system is the gamete (sperm and egg). The transmission of phenotypes and optimally adapted physiology to subsequent generations is in large part controlled by gametogenesis. In contrast to genetics, the environment actively regulates epigenetics to impact the physiology and phenotype of cellular and biological systems. The integration of epigenetics and genetics is critical for all developmental biology systems at the cellular and organism level. The current review is focused on the role of epigenetics during gametogenesis for both the spermatogenesis system in the male and oogenesis system in the female. The developmental stages from the initial primordial germ cell through gametogenesis to the mature sperm and egg are presented. How environmental factors can influence the epigenetics of gametogenesis to impact the epigenetic transgenerational inheritance of phenotypic and physiological change in subsequent generations is reviewed.

Nuclear and Mitochondrial Genome, Epigenome and Gut Microbiome: Emerging Molecular Biomarkers for Parkinson's Disease

BACKGROUND

Parkinson's disease (PD) is currently the second most common neurodegenerative disorder, burdening about 10 million elderly individuals worldwide. The multifactorial nature of PD poses a difficult obstacle for understanding the mechanisms involved in its onset and progression. Currently, diagnosis depends on the appearance of clinical signs, some of which are shared among various neurologic disorders, hindering early diagnosis. There are no effective tools to prevent PD onset, detect the disease in early stages or accurately report the risk of disease progression. Hence, there is an increasing demand for biomarkers that may identify disease onset and progression, as treatment-based medicine may not be the best approach for PD. Over the last few decades, the search for molecular markers to predict susceptibility, aid in accurate diagnosis and evaluate the progress of PD have intensified, but strategies aimed to improve individualized patient care have not yet been established.

CONCLUSIONS

Genomic variation, regulation by epigenomic mechanisms, as well as the influence of the host gut microbiome seem to have a crucial role in the onset and progress of PD, thus are considered potential biomarkers. As such, the human nuclear and mitochondrial genome, epigenome, and the host gut microbiome might be the key elements to the rise of personalized medicine for PD patients.

Structural basis of PETISCO complex assembly during piRNA biogenesis in C. elegans

In this study, Perez-Borrajero et al. set out to characterize PETISCO architecture and its interaction with RNA, together with its effector proteins TOST-1 and PID-1. Using biochemical and structural biology approaches, the authors found that PETISCO forms a dimer of tetramers, in which dimerization is mediated by both PID-3 and ERH-2. Crystal structures of the PID- 3/TOFU-6 and ERH-2/PID-3 subcomplexes reveal insights into PETISCO assembly, function, and subcellular localization. Using NMR spectroscopy, the authors also characterize the mutually exclusive interplay of ERH-2 with the two effector proteins TOST-1 and PID-1.

Piwi-interacting RNAs (piRNAs) constitute a class of small RNAs that bind PIWI proteins and are essential to repress transposable elements in the animal germline, thereby promoting genome stability and maintaining fertility. C. elegans piRNAs (21U RNAs) are transcribed individually from minigenes as precursors that require 5′ and 3′ processing. This process depends on the PETISCO complex, consisting of four proteins: IFE-3, TOFU-6, PID-3, and ERH-2. We used biochemical and structural biology approaches to characterize the PETISCO architecture and its interaction with RNA, together with its effector proteins TOST-1 and PID-1. These two proteins define different PETISCO functions: PID-1 governs 21U processing, whereas TOST-1 links PETISCO to an unknown process essential for early embryogenesis. Here, we show that PETISCO forms an octameric assembly with each subunit present in two copies. Determination of structures of the TOFU-6/PID-3 and PID-3/ERH-2 subcomplexes, supported by in vivo studies of subunit interaction mutants, allows us to propose a model for the formation of the TOFU-6/PID-3/ERH-2 core complex and its functionality in germ cells and early embryos. Using NMR spectroscopy, we demonstrate that TOST-1 and PID-1 bind to a common surface on ERH-2, located opposite its PID-3 binding site, explaining how PETISCO can mediate different cellular roles.

Functional role of piRNAs in animal models and its prospects in aquaculture

The recent advances in the field of aquaculture over the last decade has helped the cultured‐fish industry production sector to identify problems and choose the best approaches to achieve high‐volume production. Understanding the emerging roles of non‐coding RNA (ncRNA) in the regulation of fish physiology and health will assist in gaining knowledge on the possible applications of ncRNAs for the advancement of aquaculture. There is information available on the practical considerations of epigenetic mechanisms like DNA methylation, histone modification and ncRNAs, such as microRNA in aquaculture, for both fish and shellfish. Among the non‐coding RNAs, PIWI‐interacting RNA (piRNA) is 24–31 bp long transcripts, which is primarily involved in silencing the germline transposons. Besides, the burgeoning reports and studies establish piRNAs' role in various aspects of biology. Till date, there are no reviews that summarize the recent findings available on piRNAs in animal models, especially on piRNAs biogenesis and biological action. To gain a better understanding and get an overview on the process of piRNA genesis among the different animals, this work reviews the literature available on the processes of piRNA biogenesis in animal models with special reference to aquatic animal model zebrafish. This review also presents a short discussion and prospects of piRNA’s application in relevance to the aquaculture industry.

piRNAs and endo-siRNAs: Small molecules with large roles in the nervous system

Since their discovery, small non-coding RNAs have emerged as powerhouses in the regulation of numerous cellular processes. In addition to guarding the integrity of the reproductive system, small non-coding RNAs play critical roles in the maintenance of the soma. Accumulating evidence indicates that small non-coding RNAs perform vital functions in the animal nervous system such as restricting the activity of deleterious transposable elements, regulating nerve regeneration, and mediating learning and memory. In this review, we provide an overview of the current understanding of the contribution of two major classes of small non-coding RNAs, piRNAs and endo-siRNAs, to the nervous system development and function, and present highlights on how the dysregulation of small non-coding RNA pathways can assist in understanding the neuropathology of human neurological disorders.

pre-piRNA trimming and 2'-O-methylation protect piRNAs from 3' tailing and degradation in C. elegans

The Piwi-interacting RNA (piRNA) pathway suppresses transposable elements and promotes fertility in diverse organisms. Maturation of piRNAs involves pre-piRNA trimming followed by 2'-O-methylation at their 3' termini. Here, we report that the 3' termini of Caenorhabditis elegans piRNAs are subject to nontemplated nucleotide addition, and piRNAs with 3' addition exhibit extensive base-pairing interaction with their target RNAs. Animals deficient for PARN-1 (pre-piRNA trimmer) and HENN-1 (2'-O-methyltransferase) accumulate piRNAs with 3' nontemplated nucleotides. In henn-1 mutants, piRNAs are shortened prior to 3' addition, whereas long isoforms of untrimmed piRNAs are preferentially modified in parn-1 mutant animals. Loss of either PARN-1 or HENN-1 results in modest reduction in steady-state levels of piRNAs. Deletion of both enzymes leads to depletion of piRNAs, desilenced piRNA targets, and impaired fecundity. Together, our findings suggest that pre-piRNA trimming and 2'-O-methylation act collaboratively to protect piRNAs from tailing and degradation.

piRNA-30473 contributes to tumorigenesis and poor prognosis by regulating m6A RNA methylation in DLBCL

The initiation and progression of diffuse large B-cell lymphoma (DLBCL) is governed by genetic and epigenetic aberrations. As the most abundant eukaryotic messenger RNA (mRNA) modification, N6-methyladenosine (m6A) is known to influence various fundamental bioprocesses by regulating the target gene; however, the function of m6A modifications in DLBCL is unclear. PIWI-interacting RNAs (piRNAs) have been indicated to be epigenetic effectors in cancer. Here, we show that high expression of piRNA-30473 supports the aggressive phenotype of DLBCL, and piRNA-30473 depletion decreases proliferation and induces cell cycle arrest in DLBCL cells. In xenograft DLBCL models, piRNA-30473 inhibition reduces tumor growth. Moreover, piRNA-30473 is significantly associated with overall survival in a univariate analysis and is statistically significant after adjusting for the National Comprehensive Cancer Network-International Prognostic Index in the multivariate analysis. Additional studies demonstrate that piRNA-30473 exerts its oncogenic role through a mechanism involving the upregulation of WTAP, an m6A mRNA methylase, and thus enhances the global m6A level. Integrating transcriptome and m6A-sequencing analyses reveals that WTAP increases the expression of its critical target gene, hexokinase 2 (HK2), by enhancing the HK2 m6A level, thereby promoting the progression of DLBCL. Together, the piRNA-30473/WTAP/HK2 axis contributes to tumorigenesis by regulating m6A RNA methylation in DLBCL. Furthermore, by comprehensively analyzing our clinical data and data sets, we discover that the m6A regulatory genes piRNA-30473 and WTAP improve survival prediction in DLBCL patients. Our study highlights the functional importance of the m6A modification in DLBCL and might assist in the development of a prognostic stratification and therapeutic approach for DLBCL.

Small RNA expression and miRNA modification dynamics in human oocytes and early embryos

Small noncoding RNAs (sRNAs) play important roles during the oocyte-to-embryo transition (OET), when the maternal phenotype is reprogrammed and the embryo genome is gradually activated. The transcriptional program driving early human development has been studied with the focus mainly on protein-coding RNAs, and expression dynamics of sRNAs remain largely unexplored. We profiled sRNAs in human oocytes and early embryos using an RNA-sequencing (RNA-seq) method suitable for low inputs of material. We show that OET in humans is temporally coupled with the transition from predominant expression of oocyte short piRNAs (os-piRNAs) in oocytes, to activation of microRNA (miRNA) expression in cleavage stage embryos. Additionally, 3′ mono- and oligoadenylation of miRNAs is markedly increased in zygotes. We hypothesize that this may modulate the function or stability of maternal miRNAs, some of which are retained throughout the first cell divisions in embryos. This study is the first of its kind elucidating the dynamics of sRNA expression and miRNA modification along a continuous trajectory of early human development and provides a valuable data set for in-depth interpretative analyses.

The Regulation and Role of piRNAs and PIWI Proteins in Cancer

P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) are regulatory small non-coding RNAs that participate in transposon inactivation, chromatin regulation, and endogenous gene regulation. Numerous genetic and epigenetic factors regulate cell proliferation and tumor metastasis. PIWI proteins and piRNAs have been revealed to function in regulating upstream or downstream of oncogenes or tumor-suppressor genes in cancer tissues. In the present review, we summarize major recent findings in uncovering the regulation and role of PIWI proteins and piRNAs in tumorigenesis and highlight some of the promising applications of specific piRNAs in cancer therapeutics and as cancer biomarkers.

Vertebrate Lineages Exhibit Diverse Patterns of Transposable Element Regulation and Expression across Tissues

Transposable elements (TEs) comprise a major fraction of vertebrate genomes, yet little is known about their expression and regulation across tissues, and how this varies across major vertebrate lineages. We present the first comparative analysis integrating TE expression and TE regulatory pathway activity in somatic and gametic tissues for a diverse set of 12 vertebrates. We conduct simultaneous gene and TE expression analyses to characterize patterns of TE expression and TE regulation across vertebrates and examine relationships between these features. We find remarkable variation in the expression of genes involved in TE negative regulation across tissues and species, yet consistently high expression in germline tissues, particularly in testes. Most vertebrates show comparably high levels of TE regulatory pathway activity across gonadal tissues except for mammals, where reduced activity of TE regulatory pathways in ovarian tissues may be the result of lower relative germ cell densities. We also find that all vertebrate lineages examined exhibit remarkably high levels of TE-derived transcripts in somatic and gametic tissues, with recently active TE families showing higher expression in gametic tissues. Although most TE-derived transcripts originate from inactive ancient TE families (and are likely incapable of transposition), such high levels of TE-derived RNA in the cytoplasm may have secondary, unappreciated biological relevance.

The Role of Non-coding RNAs in the Pathogenesis of Glial Tumors

Among the many malignant neoplasms, glioblastoma (GBM) leads to one of the worst prognosis for patients and has an almost 100% recurrence rate. The only chemotherapeutic drug that is widely used for treating glioblastoma is temozolomide, a DNA alkylating agent. Its impact, however, is only minor; it increases patients' survival just by 12 to 14 months. Multiple highly selective compounds that affect specific proteins and have performed well in other types of cancer have proved ineffective against glioblastoma. Hence, there is an urgent need for novel methods that could help achieve the long-awaited progress in glioblastoma treatment. One of the potentially promising approaches is the targeting of non-coding RNAs (ncRNAs). These molecules are characterized by extremely high multifunctionality and often act as integrators by coordinating multiple key signaling pathways within the cell. Thus, the impact on ncRNAs has the potential to lead to a broader and stronger impact on cells, as opposed to the more focused action of inhibitors targeting specific proteins. In this review, we summarize the functions of long noncoding RNAs, circular RNAs, as well as microRNAs, PIWI-interacting RNAs, small nuclear and small nucleolar RNAs. We provide a classification of these transcripts and describe their role in various signaling pathways and physiological processes. We also provide examples of oncogenic and tumor suppressor ncRNAs belonging to each of these classes in the context of their involvement in the pathogenesis of gliomas and glioblastomas. In conclusion, we considered the potential use of ncRNAs as diagnostic markers and therapeutic targets for the treatment of glioblastoma.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Epigenetic roles of PIWI proteins and piRNAs in colorectal cancer

Small non‐coding RNAs (sncRNAs) are a subgroup of non‐coding RNAs, with less than 200 nucleotides length and no potential for coding proteins. PiRNAs, a member of sncRNAs, were first discovered more than a decade ago and have attracted researcher’s attention because of their gene regulatory function both in the nucleus and in the cytoplasm. Recent investigations have found that the abnormal expression of these sncRNAs is involved in many human diseases, including cancers. Colorectal cancer (CRC), as a common gastrointestinal malignancy, is one of the important causes of cancer‐related deaths through the entire world and appears to be a consequence of mutation in the genome and epigenetic alterations. The aim of this review is to realize whether there is a relationship between CRC and piRNAs or not.

The antiviral immunity of ticks against transmitted viral pathogens

Ticks, being obligate hematophagous arthropods, are exposed to various blood-borne pathogens, including arboviruses. Consequently, their feeding behavior can readily transmit economically important viral pathogens to humans and animals. With this tightly knit vector and pathogen interaction, the replication and transmission of tick-borne viruses (TBVs) must be highly regulated by their respective tick vectors to avoid any adverse effect on the ticks' biological development and viability. Knowledge about the tick-virus interface, although gaining relevant advances in recent years, is advancing at a slower pace than the scientific developments related to mosquito-virus interactions. The unique and complicated feeding behavior of ticks, compared to that of other blood-feeding arthropods, also limits the studies that would further elaborate the antiviral immunity of ticks against TBVs. Hence, knowledge of molecular and cellular immune mechanisms at the tick-virus interface, will further elucidate the successful viral replication of TBVs in ticks and their effective transmission to human and animal hosts.

SID-2 negatively regulates development likely independent of nutritional dsRNA uptake

RNA interference (RNAi) is a gene regulatory mechanism based on RNA-RNA interaction conserved through eukaryotes. Surprisingly, many animals can take-up human-made double stranded RNA (dsRNA) from the environment to initiate RNAi suggesting a mechanism for dsRNA-based information exchange between organisms and their environment. However, no naturally occurring example has been identified since the discovery of the phenomenon 22 years ago. Therefore it remains enigmatic why animals are able to take up dsRNA. Here, we explore other possible functions by performing phenotypic studies of dsRNA uptake deficient sid-2 mutants in Caenorhabditis elegans. We find that SID-2 does not have a nutritional role in feeding experiments using genetic sensitized mutants. Furthermore, we use robot assisted imaging to show that sid-2 mutants accelerate growth rate and, by maternal contribution, body length at hatching. Finally, we perform transcriptome and lipidome analysis showing that sid-2 has no effect on energy storage lipids, but affects signalling lipids and the embryo transcriptome. Overall, these results suggest that sid-2 has mild effects on development and is unlikely functioning in the nutritional uptake of dsRNA. These findings broaden our understanding of the biological role of SID-2 and motivate studies identifying the role of environmental dsRNA uptake.

Network-based visualisation reveals new insights into transposable element diversity

Transposable elements (TEs) are widespread across eukaryotic genomes, yet their content varies widely between different species. Factors shaping the diversity of TEs are poorly understood. Understanding the evolution of TEs is difficult because their sequences diversify rapidly and TEs are often transferred through non-conventional means such as horizontal gene transfer. We developed a method to track TE evolution using network analysis to visualise TE sequence and TE content across different genomes. We illustrate our method by first using a monopartite network to study the sequence evolution of Tc1/mariner elements across focal species. We identify a connection between two subfamilies associated with convergent acquisition of a domain from a protein-coding gene. Second, we use a bipartite network to study how TE content across species is shaped by epigenetic silencing mechanisms. We show that the presence of Piwi-interacting RNAs is associated with differences in network topology after controlling for phylogenetic effects. Together, our method demonstrates how a network-based approach can identify hitherto unknown properties of TE evolution across species.

The piRNA response to BmNPV infection in the silkworm fat body and midgut

Bombyx mori nucleopolyhedrovirus (BmNPV) is a DNA virus that causes huge losses to the silkworm industry but the piRNA responses during BmNPV infection in the silkworm remain uninvestigated. Here, silkworm piRNA profiles of uninfected and BmNPV-infected fat body and midgut were determined by high-through sequencing in the early stages of BmNPV infection. A total of 2675 and 3396 genome-derived piRNAs were identified from fat body and midgut, respectively. These genome-derived piRNAs mainly originated from unannotated instead of transposon regions in the silkworm genome. In total, 572 piRNAs were associated with 280 putative target genes in fat body and 805 piRNAs with 380 target genes in midgut. Compared to uninfected tissues, 322 and 129 piRNAs were significantly upregulated in BmNPV-infected fat body and midgut, respectively. In addition, 276 and 117 piRNAs were significantly downregulated. Moreover, differentially expressed (DE) piRNAs during BmNPV infection differed significantly between fat body and midgut. Putative DE piRNA-targeted genes were associated with "response to stimulus" and "environmental information processing" in fat body after infection with BmNPV, which may indicate an active piRNA response to BmNPV infection in fat body. This study may lay the foundation for future research of the potential roles of the piRNA pathway and specific piRNAs in BmNPV pathogenesis.

Lost and Found: Piwi and Argonaute Pathways in Flatworms

Platyhelminthes comprise one of the major phyla of invertebrate animals, inhabiting a wide range of ecosystems, and one of the most successful in adapting to parasitic life. Small non-coding RNAs have been implicated in regulating complex developmental transitions in model parasitic species. Notably, parasitic flatworms have lost Piwi RNA pathways but gained a novel Argonaute gene. Herein, we analyzed, contrasted and compared the conservation of small RNA pathways among several free-living species (a paraphyletic group traditionally known as ‘turbellarians’) and parasitic species (organized in the monophyletic clade Neodermata) to disentangle possible adaptations during the transition to parasitism. Our findings showed that complete miRNA and RNAi pathways are present in all analyzed free-living flatworms. Remarkably, whilst all ‘turbellarians’ have Piwi proteins, these were lost in parasitic Neodermantans. Moreover, two clusters of Piwi class Argonaute genes are present in all ‘turbellarians’. Interestingly, we identified a divergent Piwi class Argonaute in free living flatworms exclusively, which we named ‘Fliwi’. In addition, other key proteins of the Piwi pathways were conserved in ‘turbellarians’, while none of them were detected in Neodermatans. Besides Piwi and the canonical Argonaute proteins, a flatworm-specific class of Argonautes (FL-Ago) was identified in the analyzed species confirming its ancestrallity to all Platyhelminthes. Remarkably, this clade was expanded in parasitic Neodermatans, but not in free-living species. These phyla-specific Argonautes showed lower sequence conservation compared to other Argonaute proteins, suggesting that they might have been subjected to high evolutionary rates. However, key residues involved in the interaction with the small RNA and mRNA cleavage in the canonical Argonautes were more conserved in the FL-Agos than in the Piwi Argonautes. Whether this is related to specialized functions and adaptations to parasitism in Neodermatans remains unclear. In conclusion, differences detected in gene conservation, sequence and structure of the Argonaute family suggest tentative biological and evolutionary diversifications that are unique to Platyhelminthes. The remarkable divergencies in the small RNA pathways between free-living and parasitic flatworms indicate that they may have been involved in the adaptation to parasitism of Neodermatans.

Evolution of plasticity in production and transgenerational inheritance of small RNAs under dynamic environmental conditions

In a changing environment, small RNAs (sRNAs) play an important role in the post-transcriptional regulation of gene expression and can vary in abundance depending on the conditions experienced by an individual (phenotypic plasticity) and its parents (non-genetic inheritance). Many sRNAs are unusual in that they can be produced in two ways, either using genomic DNA as the template (primary sRNAs) or existing sRNAs as the template (secondary sRNAs). Thus, organisms can evolve rapid plastic responses to their current environment by adjusting the amplification rate of sRNA templates. sRNA levels can also be transmitted transgenerationally by the direct transfer of either sRNAs or the proteins involved in amplification. Theory is needed to describe the selective forces acting on sRNA levels, accounting for the dual nature of sRNAs as regulatory elements and templates for amplification and for the potential to transmit sRNAs and their amplification agents to offspring. Here, we develop a model to study the dynamics of sRNA production and inheritance in a fluctuating environment. We tested the selective advantage of mutants capable of sRNA-mediated phenotypic plasticity within resident populations with fixed levels of sRNA transcription. Even when the resident was allowed to evolve an optimal constant rate of sRNA production, plastic amplification rates capable of responding to environmental conditions were favored. Mechanisms allowing sRNA transcripts or amplification agents to be inherited were favored primarily when parents and offspring face similar environments and when selection acts before the optimal level of sRNA can be reached within the organism. Our study provides a clear set of testable predictions for the evolution of sRNA-related mechanisms of phenotypic plasticity and transgenerational inheritance.

The epigenome and beyond: How does non-genetic inheritance change our view of evolution?

Evidence from across the tree of life suggests that epigenetic inheritance is more common than previously thought. If epigenetic inheritance is indeed as common as the data suggest, this finding has potentially important implications for evolutionary theory and our understanding of how evolution and adaptation progress. However, we currently lack an understanding of how common various epigenetic inheritance types are, and how they impact phenotypes. In this perspective, we review the open questions that need to be addressed to fully integrate epigenetic inheritance into evolutionary theory and to develop reliable predictive models for phenotypic evolution. We posit that addressing these challenges will require the collaboration of biologists from different disciplines and a focus on the exploration of data and phenomena without preconceived limits on potential mechanisms or outcomes.