A novel class of somatic small RNAs similar to germ cell pachytene PIWI-interacting small RNAs

Short communication: Differential expression of piwi1 and piwi2 genes in tissues of tambacu and zebrafish: A possible relationship with the indeterminate muscle growth

Fish skeletal muscle is a component of the human diet, and understanding the mechanisms that control muscle growth can contribute to improving production in this sector and benefits the human health. In this sense, fish such as tambacu can represent a valuable source for exploring muscle growth regulators due to the indeterminate muscle growth pattern. In this context, the genes responsible for the indeterminate and determinate muscle growth pattern of fish are little explored, with piwi genes being possible candidates involved with these growth patterns. Piwi genes are associated with the proliferation and self-renewal of germ cells, and there are descriptions of these same functions in somatic cells from different tissues. However, little is known about the function of these genes in fish somatic cells. Considering this, our objective was to analyze the expression pattern of piwi 1 and 2 genes in cardiac muscle, skeletal muscle, liver, and gonad of zebrafish (species with determinate growth) and tambacu (species with indeterminate growth). We observed a distinct expression of piwi1 and piwi2 between tambacu and zebrafish, with both genes more expressed in tambacu in all tissues evaluated. Piwi genes can represent potential candidates involved with indeterminate muscle growth control.

Non-coding RNAs regulate mitochondrial dynamics in the development of gastric cancer

Gastric cancer (GC) is a malignant cancer that reduces life expectancy worldwide. Although treatment strategies have improved, patients with GC still have poor prognoses. Hence, it is necessary to understand the molecular mechanisms of GC and to find new therapeutic targets. Mitochondrial dynamics and mitochondrial dysfunction are associated with cancer cell growth and progression. Numerous studies have reported that non-coding RNAs (ncRNAs) can participate in the occurrence and development of GC by regulating mitochondrial dynamics. Elucidating the crosstalk between ncRNAs and mitochondria would be helpful in preventing and treating GC. Herein, we review and summarize the functions of oncogenes and tumor suppressors in suppressing ncRNAs and regulating mitochondrial dynamics in GC tumor growth, proliferation, invasion and metastasis. This review provides new insights into the pathogenesis of and intervention for GC.

Unraveling mitochondrial piRNAs in mouse embryonic gonadal cells

Although mitochondria are widely studied organelles, the recent interest in the role of mitochondrial small noncoding RNAs (sncRNAs), miRNAs, and more recently, piRNAs, is providing new functional perspectives in germ cell development and differentiation. piRNAs (PIWI-interacting RNAs) are single-stranded sncRNAs of mostly about 20-35 nucleotides, generated from the processing of pre-piRNAs. We leverage next-generation sequencing data obtained from mouse primordial germ cells and somatic cells purified from early-differentiating embryonic ovaries and testis from 11.5 to 13.5 days postcoitum. Using bioinformatic tools, we elucidate (i) the origins of piRNAs as transcribed from mitochondrial DNA fragments inserted in the nucleus or from the mitochondrial genome; (ii) their levels of expression; and (iii) their potential roles, as well as their association with genomic regions encoding other sncRNAs (such as tRNAs and rRNAs) and the mitochondrial regulatory region (D-loop). Finally, our results suggest how nucleo-mitochondrial communication, both anterograde and retrograde signaling, may be mediated by mitochondria-associated piRNAs.

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.

PIWI-mediated control of tissue-specific transposons is essential for somatic cell differentiation

PIWI proteins are known as mediators of transposon silencing in animal germlines but are also found in adult pluripotent stem cells of highly regenerative animals, where they are essential for regeneration. Study of the nuclear PIWI protein SMEDWI-2 in the planarian somatic stem cell system reveals an intricate interplay between transposons and cell differentiation in which a subset of transposons is inevitably activated during cell differentiation, and the PIWI protein is required to regain control. Absence of SMEDWI-2 leads to tissue-specific transposon derepression related to cell-type-specific chromatin remodeling events and in addition causes reduced accessibility of lineage-specific genes and defective cell differentiation, resulting in fatal tissue dysfunction. Finally, we show that additional PIWI proteins provide a stem-cell-specific second layer of protection in planarian neoblasts. These findings reveal a far-reaching role of PIWI proteins and PIWI-interacting RNAs (piRNAs) in stem cell biology and cell differentiation.

Knockout of glutathione peroxidase 5 down-regulates the piRNAs in the caput epididymidis of aged mice

The mammalian epididymis not only plays a fundamental role in the maturation of spermatozoa, but also provides protection against various stressors. The foremost among these is the threat posed by oxidative stress, which arises from an imbalance in reactive oxygen species and can elicit damage to cellular lipids, proteins, and nucleic acids. In mice, the risk of oxidative damage to spermatozoa is mitigated through the expression and secretion of glutathione peroxidase 5 (GPX5) as a major luminal scavenger in the proximal caput epididymidal segment. Accordingly, the loss of GPX5-mediated protection leads to impaired DNA integrity in the spermatozoa of aged Gpx5^-/- mice. To explore the underlying mechanism, we have conducted transcriptomic analysis of caput epididymidal epithelial cells from aged (13 months old) Gpx5^-/- mice. This analysis revealed the dysregulation of several thousand epididymal mRNA transcripts, including the downregulation of a subgroup of piRNA pathway genes, in aged Gpx5^-/- mice. In agreement with these findings, we also observed the loss of piRNAs, which potentially bind to the P-element-induced wimpy testis (PIWI)-like proteins PIWIL1 and PIWIL2. The absence of these piRNAs was correlated with the elevated mRNA levels of their putative gene targets in the caput epididymidis of Gpx5^-/- mice. Importantly, the oxidative stress response genes tend to have more targeting piRNAs, and many of them were among the top increased genes upon the loss of GPX5. Taken together, our findings suggest the existence of a previously uncharacterized somatic piRNA pathway in the mammalian epididymis and its possible involvement in the aging and oxidative stress-mediated responses.

Dynamics of cattle sperm sncRNAs during maturation, from testis to ejaculated sperm

Background

During epididymal transit, spermatozoa go through several functional maturation steps, resulting from interactions with epididymal secretomes specific to each region. In particular, the sperm membrane is under constant remodeling, with sequential attachment and shedding of various molecules provided by the epididymal lumen fluid and epididymosomes, which also deliver sncRNA cargo to sperm. As a result, the payload of sperm sncRNAs changes during the transit from the epididymis caput to the cauda. This work was designed to study the dynamics of cattle sperm sncRNAs from spermatogenesis to final maturation.

Results

Comprehensive catalogues of sperm sncRNAs were obtained from testicular parenchyma, epididymal caput, corpus and cauda, as well as ejaculated semen from three Holstein bulls. The primary cattle sncRNA sperm content is markedly remodeled as sperm mature along the epididymis. Expression of piRNAs, which are abundant in testis parenchyma, decreases dramatically at epididymis. Conversely, sperm progressively acquires miRNAs, rsRNAs, and tsRNAs along epididymis, with regional specificities. For instance, miRNAs and tsRNAs are enriched in epididymis cauda and ejaculated sperm, while rsRNA expression peaks at epididymis corpus. In addition, epididymis corpus contains mainly 20 nt long piRNAs, instead of 30 nt in all other locations. Beyond the bulk differences in abundance of sncRNAs classes, K-means clustering was performed to study their spatiotemporal expression profile, highlighting differences in specific sncRNAs and providing insights into their putative biological role at each maturation stage. For instance, Gene Ontology analyses using miRNA targets highlighted enriched processes such as cell cycle regulation, response to stress and ubiquitination processes in testicular parenchyma, protein metabolism in epididymal sperm, and embryonic morphogenesis in ejaculated sperm.

Conclusions

Our findings confirm that the sperm sncRNAome does not simply reflect a legacy of spermatogenesis. Instead, sperm sncRNA expression shows a remarkable level of plasticity resulting probably from the combination of multiple factors such as loss of the cytoplasmic droplet, interaction with epididymosomes, and more surprisingly, the putative in situ production and/or modification of sncRNAs by sperm. Given the suggested role of sncRNA in epigenetic trans-generational inheritance, our detailed spatiotemporal analysis may pave the way for a study of sperm sncRNAs role in embryo development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00397-5.

PIWI-interacting RNAs in human cancer

P-element-induced wimpy testis (PIWI) interacting RNAs (piRNAs) are a class of small regulatory RNAs mechanistically similar to but much less studied than microRNAs and small interfering RNAs. Today the best understood function of piRNAs is transposon control in animal germ cells, which has earned them the name 'guardians of the germline'. Several molecular/cellular characteristics of piRNAs, including high sequence diversity, lack of secondary structures, and target-oriented generation seem to serve this purpose. Recently, aberrant expressions of piRNAs and PIWI proteins have been implicated in a variety of malignant tumors and associated with cancer hallmarks such as cell proliferation, inhibited apoptosis, invasion, metastasis and increased stemness. Researchers have also demonstrated multiple mechanisms of piRNA-mediated target deregulation associated with cancer initiation, progression or dissemination. We review current research findings on the biogenesis, normal functions and cancer associations of piRNAs, highlighting their potentials as cancer diagnostic/prognostic biomarkers and therapeutic tools. Whenever applicable, we draw connections with other research fields to encourage intercommunity conversations. We also offer recommendations and cautions regarding the general process of cancer-related piRNA studies and the methods/tools used at each step. Finally, we call attention to some issues that, if left unsolved, might impede the future development of this field.

piRDisease v1.0: a manually curated database for piRNA associated diseases

In recent years, researches focusing on PIWI-interacting RNAs (piRNAs) have increased rapidly. It has been revealed that piRNAs have strong association with a wide range of diseases; thus, it becomes very important to understand piRNAs’ role(s) in disease diagnosis, prognosis and assessment of treatment response. We searched more than 2500 articles using keywords, such as `PIWI-interacting RNAs’ and `piRNAs’, and further scrutinized the articles to collect piRNAs-disease association data. These data are highly complex and heterogeneous due to various types of piRNA idnetifiers (IDs) and different reference genome versions. We put considerable efforts into removing redundancy and anomalies and thus homogenized the data. Finally, we developed the piRDisease database, which incorporates experimentally supported data for piRNAs’ relationship with wide range of diseases. The piRDisease (piRDisease v1.0) is a novel, comprehensive and exclusive database resource, which provides 7939 manually curated associations of experimentally supported 4796 piRNAs involved in 28 diseases. piRDisease facilitates users by providing detailed information of the piRNA in respective disease, explored by experimental support, brief description, sequence and location information. Considering piRNAs’ role(s) in wide range of diseases, it is anticipated that huge amount of data would be produced in the near future. We thus offer a submitting page, on which users or researches can contribute in to update our piRDisease database.

Functions and mechanism of noncoding RNA in the somatic cells of the testis

ncRNAs are involved in numerous biological processes by regulating gene expression and cell stability. Studies have shown that ncRNAs also contribute to spermatogenesis. Leydig cells (LCs) and Sertoli cells (SCs) are somatic cells of the testis that support spermatogenesis and are vital to male fertility. In this review, we summarized the findings from studies on ncRNAs in SCs and LCs. In SCs, ncRNAs play key roles in phagocytosis, immunoprotection and development of SCs. In LCs, ncRNAs are involved in steroidogenesis, in particular production of testosterone as well as development of LCs. Here, we discuss the possible target genes and functions of ncRNAs in both types of cells. These ncRNAs regulate the expression of target genes or mRNA coding sequence regions, resulting in a chain reaction that influences cell function. In addition, microRNAs, lncRNAs, piRNA-like RNAs (pilRNAs) and natural antisense transcripts (NATs) are discussed in this review. In summary, we suggest that these ncRNAs might act in coordination to control spermatogenesis and maintain the environmental homeostasis of the testis.

Aptamer-Functionalized Activatable DNA Tetrahedron Nanoprobe for PIWI-Interacting RNA Imaging and Regulating in Cancer Cells

It has been reported that PIWI-interacting RNAs (piRNAs) play critical roles in activating invasion and metastasis, evading growth suppressors, and sustaining proliferative signaling of cancer and can be regarded as a novel biomarker candidate. Thus, it is necessary to develop an effective method for imaging and regulating cancer-related piRNAs to diagnose and treat cancers. Herein, we designed aptamer-functionalized activatable DNA tetrahedron nanoprobes (apt-ADTNs) to image and regulate endogenous piRNAs in cancer cells. As proof of concept, overexpressed piRNA-36026 in MCF-7 cells was used for this study. In brief, aptamer AS1411 and piRNA-36026 antisequence with Cy5 fluorescent dye are appended from the DNA tetrahedron; then, a short oligonucleotide with black hole quencher 2 (Q-oligo) is complementary with piRNA-36026 antisequence to quench the fluorescence of Cy5. The apt-ADTNs can recognize the MCF-7 cells through aptamer AS1411, and then enter the cells. Q-oligo is detached from the apt-ADTNs because of the binding between apt-ADTNs and piRNA-36026, leading to the recovery of the Cy5 fluorescence signal. Meanwhile, the hybridization of apt-ADTNs and piRNA-36026 results in down-regulating of dissociative piRNA-36026 in cytoplasm and the subsequent apoptosis of MCF-7 cells. As the achievement of synchronously imaging and regulating piRNA-36026 in MCF-7 cells, we believe that this design holds great promise in application of diagnosis and therapy for cancer.

Theragnosis by a miR-141-3p molecular beacon: simultaneous detection and sensitization of 5-fluorouracil resistant colorectal cancer cells through the activation of the TRIM13-associated apoptotic pathway

We developed a molecular beacon targeting miR-141-3p, aberrantly increased in 5-fluorouracil-resistant colorectal cancer cells (R-CRCCs). It consists of a fluorophore-labeled oligonucleotide, antisense to miR-141-3p, and a quencher. It detected R-CRCCs and recovered the chemosensitivity of them to 5-fluorouracil by hybridization with miR-141-3p, which is applicable to cancer treatment.

Genome sequencing analysis of blood cells identifies germline haplotypes strongly associated with drug resistance in osteosarcoma patients

BACKGROUND

Osteosarcoma is the most common malignant bone tumor in children. Survival remains poor among histologically poor responders, and there is a need to identify them at diagnosis to avoid delivering ineffective therapy. Genetic variation contributes to a wide range of response and toxicity related to chemotherapy. The aim of this study is to use sequencing of blood cells to identify germline haplotypes strongly associated with drug resistance in osteosarcoma patients.

METHODS

We used sequencing data from two patient datasets, from Inova Hospital and the NCI TARGET. We explored the effect of mutation hotspots, in the form of haplotypes, associated with relapse outcome. We then mapped the single nucleotide polymorphisms (SNPs) in these haplotypes to genes and pathways. We also performed a targeted analysis of mutations in Drug Metabolizing Enzymes and Transporter (DMET) genes associated with tumor necrosis and survival.

RESULTS

We found intronic and intergenic hotspot regions from 26 genes common to both the TARGET and INOVA datasets significantly associated with relapse outcome. Among significant results were mutations in genes belonging to AKR enzyme family, cell-cell adhesion biological process and the PI3K pathways; as well as variants in SLC22 family associated with both tumor necrosis and overall survival. The SNPs from our results were confirmed using Sanger sequencing. Our results included known as well as novel SNPs and haplotypes in genes associated with drug resistance.

CONCLUSION

We show that combining next generation sequencing data from multiple datasets and defined clinical data can better identify relevant pathway associations and clinically actionable variants, as well as provide insights into drug response mechanisms.

Somatic expression of piRNA and associated machinery in the mouse identifies short, tissue-specific piRNA

Piwi-interacting RNAs (piRNAs) are small non-coding RNAs that associate with PIWI proteins for transposon silencing via DNA methylation and are highly expressed and extensively studied in the germline. Mature germline piRNAs typically consist of 24-32 nucleotides, with a strong preference for a 5' uridine signature, an adenosine signature at position 10, and a 2'-O-methylation signature at the 3' end. piRNA presence in somatic tissues, however, is not well characterized and requires further systematic evaluation. In the current study, we identified piRNAs and associated machinery from mouse somatic tissues representing the three germ layers. piRNA specificity was improved by combining small RNA size selection, sodium periodate treatment enrichment for piRNA over other small RNA, and small RNA next-generation sequencing. We identify PIWIL1, PIWIL2, and PIWIL4 expression in brain, liver, kidney, and heart. Of note, somatic piRNAs are shorter in length and tissue-specific, with increased occurrence of unique piRNAs in hippocampus and liver, compared to the germline. Hippocampus contains 5,494 piRNA-like peaks, the highest expression among all tested somatic tissues, followed by cortex (1,963), kidney (580), and liver (406). The study identifies 26 piRNA sequence species and 40 piRNA locations exclusive to all examined somatic tissues. Although piRNA expression has long been considered exclusive to the germline, our results support that piRNAs are expressed in several somatic tissues that may influence piRNA functions in the soma. Once confirmed, the PIWI/piRNA system may serve as a potential tool for future research in epigenome editing to improve human health by manipulating DNA methylation.

The meiotic phosphatase GSP-2/PP1 promotes germline immortality and small RNA-mediated genome silencing

Germ cell immortality, or transgenerational maintenance of the germ line, could be promoted by mechanisms that could occur in either mitotic or meiotic germ cells. Here we report for the first time that the GSP-2 PP1/Glc7 phosphatase promotes germ cell immortality. Small RNA-induced genome silencing is known to promote germ cell immortality, and we identified a separation-of-function allele of C. elegans gsp-2 that is compromised for germ cell immortality and is also defective for small RNA-induced genome silencing and meiotic but not mitotic chromosome segregation. Previous work has shown that GSP-2 is recruited to meiotic chromosomes by LAB-1, which also promoted germ cell immortality. At the generation of sterility, gsp-2 and lab-1 mutant adults displayed germline degeneration, univalents, histone methylation and histone phosphorylation defects in oocytes, phenotypes that mirror those observed in sterile small RNA-mediated genome silencing mutants. Our data suggest that a meiosis-specific function of GSP-2 ties small RNA-mediated silencing of the epigenome to germ cell immortality. We also show that transgenerational epigenomic silencing at hemizygous genetic elements requires the GSP-2 phosphatase, suggesting a functional link to small RNAs. Given that LAB-1 localizes to the interface between homologous chromosomes during pachytene, we hypothesize that small localized discontinuities at this interface could promote genomic silencing in a manner that depends on small RNAs and the GSP-2 phosphatase.

Bioimaging of multiple piRNAs in a single breast cancer cell using molecular beacons

Simultaneous bioimaging of piR-36026 and piR-36743 using molecular beacons successfully visualized 4 different subtypes of breast cancer.

The roles of non-coding RNAs in cardiac regenerative medicine

The emergence of non-coding RNAs (ncRNAs) has challenged the central dogma of molecular biology that dictates that the decryption of genetic information starts from transcription of DNA to RNA, with subsequent translation into a protein. Large numbers of ncRNAs with biological significance have now been identified, suggesting that ncRNAs are important in their own right and their roles extend far beyond what was originally envisaged. ncRNAs do not only regulate gene expression, but are also involved in chromatin architecture and structural conformation. Several studies have pointed out that ncRNAs participate in heart disease; however, the functions of ncRNAs still remain unclear. ncRNAs are involved in cellular fate, differentiation, proliferation and tissue regeneration, hinting at their potential therapeutic applications. Here, we review the current understanding of both the biological functions and molecular mechanisms of ncRNAs in heart disease and describe some of the ncRNAs that have potential heart regeneration effects.

Paternal pachytene piRNAs are not required for fertilization, embryonic development and sperm-mediated epigenetic inheritance in mice

Pachytene piRNAs are MIWI-/MILI-bound small RNAs abundantly expressed in pachytene spermatocytes and round spermatids in adult mouse testes. Miwi knockout (KO) male mice are sterile due to spermiogenic arrest. In Caenorhabditis elegans, sperm-borne piRNAs appear to have an epigenetic role during fertilization and development because progeny of individuals with piRNA-deficient gametes display a progressive loss of fertility after several generations. In mice, it remains unknown whether pachytene piRNA-deficient round spermatids can produce offspring, and whether the progeny of Miwi mutants also exhibits transgenerational, progressive fertility loss. Here, we report that Miwi KO round spermatids could fertilize both wild-type (WT) and Miwi KO oocytes through round spermatid injection, and could produce healthy and fertile offspring despite the global downregulation of both MIWI-/MILI-bound pachytene piRNAs. Progeny of ROSI-derived heterozygotes, both male and female, displayed normal fertility for at least three generations when bred with either WT or Miwi KO females. Our data indicate that aberrant MIWI-/MILI-bound pachytene piRNA profiles in spermatids do not affect fertilization, early embryonic development, or fertility of the offspring, suggesting that pachytene piRNAs might not be required for paternal transgenerational epigenetic inheritance in mice.

SpermBase: A Database for Sperm-Borne RNA Contents

Since their discovery approximately three decades ago, sperm-borne RNAs, both large/small and coding/noncoding, have been reported in multiple organisms, and some have been implicated in spermatogenesis, early development, and epigenetic inheritance. Despite these advances, isolation, quantification, and annotation of sperm-borne RNAs remain nontrivial. The yields and subspecies of sperm-borne RNAs isolated from sperm can vary drastically depending on the methods used, and no cross-species analyses of sperm RNA contents have ever been conducted using a standardized sperm RNA isolation protocol. To address these issues, we developed a simple RNA isolation method that is applicable to sperm of various species, thus allowing for reliable interspecies comparisons. Based on RNA-Seq analyses, we established SpermBase ( www.spermbase.org), a database dedicated to sperm-borne RNA profiling of multiple species. Currently, SpermBase contains large and small RNA expression data for mouse, rat, rabbit, and human total sperm and sperm heads. By analyzing large and small RNAs for conserved features, we found that many sperm-borne RNA species were conserved across all four species analyzed, and among the conserved small RNAs, sperm-borne tRNA-derived small noncoding RNAs and miRNAs can target a large number of genes known to be critical for early development.

Multiplex bioimaging of piRNA molecular pathway-regulated theragnostic effects in a single breast cancer cell using a piRNA molecular beacon

Recently, PIWI-interacting small non-coding RNAs (piRNAs) have emerged as novel cancer biomarkers candidate because of their high expression level in various cancer types and role in the control of tumor suppressor genes. In this study, a novel breast cancer theragnostics probe based on a single system targeting the piRNA-36026 (piR-36026) molecular pathway was developed using a piR-36026 molecular beacon (MB). The piR-36026 MB successfully visualized endogenous piR-36026 biogenesis, which is highly expressed in MCF7 cells (a human breast cancer cell line), and simultaneously inhibited piR-36026-mediated cancer progression in vitro and in vivo. We discovered two tumor suppressor proteins, SERPINA1 and LRAT, that were directly regulated as endogenous piR-36026 target genes in MCF7 cells. Furthermore, multiplex bioimaging of a single MCF7 cell following treatment with piR-36026 MB clearly visualized the direct molecular interaction of piRNA-36026 with SERPINA1 or LRAT and subsequent molecular therapeutic responses including caspase-3 and PI in the nucleus.

Genome-wide differential expression of genes and small RNAs in testis of two different porcine breeds and at two different ages

Some documented evidences proved small RNAs (sRNA) and targeted genes are involved in mammalian testicular development and spermatogenesis. However, the detailed molecular regulation mechanisms of them remain largely unknown so far. In this study, we obtained a total of 10,716 mRNAs, 67 miRNAs and 16,953 piRNAs which were differentially expressed between LC and LW pig breeds or between the two sexual maturity stages. Of which, we identified 16 miRNAs and 28 targeted genes possibly related to spermatogenesis; 14 miRNA and 18 targeted genes probably associated with cell adhesion related testis development. We also annotated 579 piRNAs which could potentially regulate cell death, nucleosome organization and other basic biology process, which implied that those piRNAs might be involved in sexual maturation difference. The integrated network analysis results suggested that some differentially expressed genes were involved in spermatogenesis through the ECM-receptor interaction, focal adhesion, Wnt and PI3K-Akt signaling pathways, some particular miRNAs have the negative regulation roles and some special piRNAs have the positive and negative regulation roles in testicular development. Our data provide novel insights into the molecular expression and regulation similarities and diversities of spermatogenesis and testicular development in different pig breeds at different stages of sexual maturity.

Ancestral vinclozolin exposure alters the epigenetic transgenerational inheritance of sperm small noncoding RNAs

Exposure to the agricultural fungicide vinclozolin during gestation promotes a higher incidence of various diseases in the subsequent unexposed F3 and F4 generations. This phenomenon is termed epigenetic transgenerational inheritance and has been shown to in part involve alterations in DNA methylation, but the role of other epigenetic mechanisms remains unknown. The current study investigated the alterations in small noncoding RNA (sncRNA) in the sperm from F3 generation control and vinclozolin lineage rats. Over 200 differentially expressed sncRNAs were identified and the tRNA-derived sncRNAs, namely 5' halves of mature tRNAs (5' halves), displayed the most dramatic changes. Gene targets of the altered miRNAs and tRNA 5' halves revealed associations between the altered sncRNAs and differentially DNA methylated regions. Dysregulated sncRNAs appear to correlate with mRNA profiles associated with the previously observed vinclozolin-induced disease phenotypes. Data suggest potential connections between sperm-borne RNAs and the vinclozolin-induced epigenetic transgenerational inheritance phenomenon.

piRNAs derived from ancient viral processed pseudogenes as transgenerational sequence-specific immune memory in mammals

Endogenous bornavirus-like nucleoprotein elements (EBLNs) are sequences within vertebrate genomes derived from reverse transcription and integration of ancient bornaviral nucleoprotein mRNA via the host retrotransposon machinery. While species with EBLNs appear relatively resistant to bornaviral disease, the nature of this association is unclear. We hypothesized that EBLNs could give rise to antiviral interfering RNA in the form of PIWI-interacting RNAs (piRNAs), a class of small RNA known to silence transposons but not exogenous viruses. We found that in both rodents and primates, which acquired their EBLNs independently some 25-40 million years ago, EBLNs are present within piRNA-generating regions of the genome far more often than expected by chance alone (ℙ = 8 × 10(-3)-6 × 10(-8)). Three of the seven human EBLNs fall within annotated piRNA clusters and two marmoset EBLNs give rise to bona fide piRNAs. In both rats and mice, at least two of the five EBLNs give rise to abundant piRNAs in the male gonad. While no EBLNs are syntenic between rodent and primate, some of the piRNA clusters containing EBLNs are; thus we deduce that EBLNs were integrated into existing piRNA clusters. All true piRNAs derived from EBLNs are antisense relative to the proposed ancient bornaviral nucleoprotein mRNA. These observations are consistent with a role for EBLN-derived piRNA-like RNAs in interfering with ancient bornaviral infection. They raise the hypothesis that retrotransposon-dependent virus-to-host gene flow could engender RNA-mediated, sequence-specific antiviral immune memory in metazoans analogous to the CRISPR/Cas system in prokaryotes.

Unique somatic and malignant expression patterns implicate PIWI-interacting RNAs in cancer-type specific biology

Human PIWI-interacting RNAs (piRNAs) are known to be expressed in germline cells, functionally silencing LINEs and SINEs. Their expression patterns in somatic tissues are largely uncharted. We analyzed 6,260 human piRNA transcriptomes derived from non-malignant and tumour tissues from 11 organs. We discovered that only 273 of the 20,831 known piRNAs are expressed in somatic non-malignant tissues. However, expression patterns of these piRNAs were able to distinguish tissue-of-origin. A total of 522 piRNAs are expressed in corresponding tumour tissues, largely distinguishing tumour from non-malignant tissues in a cancer-type specific manner. Most expressed piRNAs mapped to known transcripts, contrary to “piRNA clusters” reported in germline cells. We showed that piRNA expression can delineate clinical features, such as histological subgroups, disease stages, and survival. PiRNAs common to many cancer types might represent a core gene-set that facilitates cancer growth, while piRNAs unique to individual cancer types likely contribute to cancer-specific biology.