Defining piRNA primary transcripts

A comparative roadmap of PIWI-interacting RNAs across seven species reveals insights into de novo piRNA-precursor formation in mammals

PIWI-interacting RNAs (piRNAs) play a crucial role in safeguarding genome integrity by silencing mobile genetic elements. From flies to humans, piRNAs originate from long single-stranded precursors encoded by genomic piRNA clusters. How piRNA clusters form to adapt to genomic invaders and evolve to maintain protection remain key outstanding questions. Here, we generate a roadmap of piRNA clusters across seven species that highlights both similarities and variations. In mammals, we identify transcriptional readthrough as a mechanism to generate piRNAs from transposon insertions (piCs) downstream of genes (DoG). Together with the well-known stress-dependent DoG transcripts, our findings suggest a molecular mechanism for the formation of piRNA clusters in response to retroviral invasion. Finally, we identify a class of dynamic piRNA clusters in humans, underscoring unique features of human germ cell biology. Our results advance the understanding of conserved principles and species-specific variations in piRNA biology and provide tools for future studies.

A cascade of bHLH-regulated pathways programs maize anther development

The spatiotemporal development of somatic tissues of the anther lobe is necessary for successful fertile pollen production. This process is mediated by many transcription factors acting through complex, multi-layered networks. Here, our analysis of functional knockout mutants of interacting basic helix-loop-helix genes Ms23, Ms32, basic helix-loop-helix 122 (bHLH122), and bHLH51 in maize (Zea mays) established that male fertility requires all four genes, expressed sequentially in the tapetum (TP). Not only do they regulate each other, but also they encode proteins that form heterodimers that act collaboratively to guide many cellular processes at specific developmental stages. MS23 is confirmed to be the master factor, as the ms23 mutant showed the earliest developmental defect, cytologically visible in the TP, with the most drastic alterations in premeiotic gene expression observed in ms23 anthers. Notably, the male-sterile ms23, ms32, and bhlh122-1 mutants lack 24-nt phased secondary small interfering RNAs (phasiRNAs) and the precursor transcripts from the corresponding 24-PHAS loci, while the bhlh51-1 mutant has wild-type levels of both precursors and small RNA products. Multiple lines of evidence suggest that 24-nt phasiRNA biogenesis primarily occurs downstream of MS23 and MS32, both of which directly activate Dcl5 and are required for most 24-PHAS transcription, with bHLH122 playing a distinct role in 24-PHAS transcription.

Promoter hypermethylation of PIWI/piRNA pathway genes associated with diminished pachytene piRNA production in bovine hybrid male sterility

Hybrid male sterility (HMS) is a postzygotic reproductive isolation mechanism that enforces speciation. A bovine example of HMS is the yattle (also called dzo), an interspecies hybrid of taurine cattle (Bos taurus) and yak (Bos grunniens). The molecular mechanisms underlying HMS of yattle are not well understood. Epigenetic modifications of DNA methylation and P-element induced wimpy testis (PIWI)-interacting RNA (piRNAs) are important regulators in spermatogenesis. In this study, we investigated DNA methylation patterns and piRNA expression in adult testes in hybrid infertile yattle bulls and fertile cattle and yak bulls using whole genome bisulphite-seq and small RNA-seq. Promoter hypermethylation in yattle were associated with DNA methylation involved in gamete generation, piRNA metabolic processes, spermatogenesis, and spermatid development (P < 2.6 × 10^-5). Male infertility in yattle was associated with the promoter hypermethylation-associated silencing of PIWI/piRNA pathway genes including PIWIL1, DDX4, PLD6, MAEL, FKBP6, TDRD1 and TDRD5. The downstream effects of silencing these genes were diminished production of 29- to 31- nucleotide pachytene piRNAs in yattle testes. Hypermethylation events at transposable element loci (LINEs, SINEs, and LTRs) were found in yattle. LINE-derived prepachytene piRNAs increased and SINE-derived prepachytene piRNAs were reduced in yattle testes. Our data suggests that DNA methylation affects the PIWI/piRNA pathway and is involved in gene expression and pachytene piRNA production during spermatogenesis in bovine HMS. DNA hypermethylation and disruption of piRNA production contributed to unsuccessful germ cell development that may drive bovine HMS.

The Translational Status of Cancer Liquid Biopsies

Abstract

Precision oncology aims to tailor clinical decisions specifically to patients with the objective of improving treatment outcomes. This can be achieved by leveraging omics information for accurate molecular characterization of tumors. Tumor tissue biopsies are currently the main source of information for molecular profiling. However, biopsies are invasive and limited in resolving spatiotemporal heterogeneity in tumor tissues. Alternative non-invasive liquid biopsies can exploit patient’s body fluids to access multiple layers of tumor-specific biological information (genomes, epigenomes, transcriptomes, proteomes, metabolomes, circulating tumor cells, and exosomes). Analysis and integration of these large and diverse datasets using statistical and machine learning approaches can yield important insights into tumor biology and lead to discovery of new diagnostic, predictive, and prognostic biomarkers. Translation of these new diagnostic tools into standard clinical practice could transform oncology, as demonstrated by a number of liquid biopsy assays already entering clinical use. In this review, we highlight successes and challenges facing the rapidly evolving field of cancer biomarker research.

Lay Summary

Precision oncology aims to tailor clinical decisions specifically to patients with the objective of improving treatment outcomes. The discovery of biomarkers for precision oncology has been accelerated by high-throughput experimental and computational methods, which can inform fine-grained characterization of tumors for clinical decision-making. Moreover, advances in the liquid biopsy field allow non-invasive sampling of patient’s body fluids with the aim of analyzing circulating biomarkers, obviating the need for invasive tumor tissue biopsies. In this review, we highlight successes and challenges facing the rapidly evolving field of liquid biopsy cancer biomarker research.

BPLLDA: Predicting lncRNA-Disease Associations Based on Simple Paths With Limited Lengths in a Heterogeneous Network

In recent years, it has been increasingly clear that long noncoding RNAs (lncRNAs) play critical roles in many biological processes associated with human diseases. Inferring potential lncRNA-disease associations is essential to reveal the secrets behind diseases, develop novel drugs, and optimize personalized treatments. However, biological experiments to validate lncRNA-disease associations are very time-consuming and costly. Thus, it is critical to develop effective computational models. In this study, we have proposed a method called BPLLDA to predict lncRNA-disease associations based on paths of fixed lengths in a heterogeneous lncRNA-disease association network. Specifically, BPLLDA first constructs a heterogeneous lncRNA-disease network by integrating the lncRNA-disease association network, the lncRNA functional similarity network, and the disease semantic similarity network. It then infers the probability of an lncRNA-disease association based on paths connecting them and their lengths in the network. Compared to existing methods, BPLLDA has a few advantages, including not demanding negative samples and the ability to predict associations related to novel lncRNAs or novel diseases. BPLLDA was applied to a canonical lncRNA-disease association database called LncRNADisease, together with two popular methods LRLSLDA and GrwLDA. The leave-one-out cross-validation areas under the receiver operating characteristic curve of BPLLDA are 0.87117, 0.82403, and 0.78528, respectively, for predicting overall associations, associations related to novel lncRNAs, and associations related to novel diseases, higher than those of the two compared methods. In addition, cervical cancer, glioma, and non-small-cell lung cancer were selected as case studies, for which the predicted top five lncRNA-disease associations were verified by recently published literature. In summary, BPLLDA exhibits good performances in predicting novel lncRNA-disease associations and associations related to novel lncRNAs and diseases. It may contribute to the understanding of lncRNA-associated diseases like certain cancers.

Two modes of targeting transposable elements by piRNA pathway in human testis

PIWI proteins and their partner small RNAs, termed piRNAs, are known to control transposable elements (TEs) in the germline. Here, we provide evidence that in humans this control is exerted in two different modes. On the one hand, production of piRNAs specifically targeting evolutionarily youngest TEs (L1HS, L1PA2-L1PA6, LTR12C, SVA) is present both at prenatal and postnatal stages of spermatogenesis and is performed without involvement of piRNA clusters. On the other hand, at postnatal stages, piRNAs deriving from pachytene clusters target "older" TEs and thus complement cluster-independent piRNA production to achieve relevant targeting of virtually all TEs expressed in postnatal testis. We also find that converging transcription of antisense-oriented genes contributes to the origin of genic postnatal prepachytene clusters. Finally, while a fraction of pachytene piRNAs was previously shown to arise from long intergenic noncoding RNAs (lincRNAs, i.e., pachytene piRNA cluster primary transcripts), we ascertain that these are a specific set of lincRNAs that both possess distinguishing epigenetic features and are expressed exclusively in testis.

Global network random walk for predicting potential human lncRNA-disease associations

There is more and more evidence that the mutation and dysregulation of long non-coding RNA (lncRNA) are associated with numerous diseases, including cancers. However, experimental methods to identify associations between lncRNAs and diseases are expensive and time-consuming. Effective computational approaches to identify disease-related lncRNAs are in high demand; and would benefit the detection of lncRNA biomarkers for disease diagnosis, treatment, and prevention. In light of some limitations of existing computational methods, we developed a global network random walk model for predicting lncRNA-disease associations (GrwLDA) to reveal the potential associations between lncRNAs and diseases. GrwLDA is a universal network-based method and does not require negative samples. This method can be applied to a disease with no known associated lncRNA (isolated disease) and to lncRNA with no known associated disease (novel lncRNA). The leave-one-out cross validation (LOOCV) method was implemented to evaluate the predicted performance of GrwLDA. As a result, GrwLDA obtained reliable AUCs of 0.9449, 0.8562, and 0.8374 for overall, novel lncRNA and isolated disease prediction, respectively, significantly outperforming previous methods. Case studies of colon, gastric, and kidney cancers were also implemented, and the top 5 disease-lncRNA associations were reported for each disease. Interestingly, 13 (out of the 15) associations were confirmed by literature mining.

Artificial "ping-pong" cascade of PIWI-interacting RNA in silkworm cells

PIWI-interacting RNAs (piRNAs) play essential roles in the defense system against selfish elements in animal germline cells by cooperating with PIWI proteins. A subset of piRNAs is predicted to be generated via the "ping-pong" cascade, which is mainly controlled by two different PIWI proteins. Here we established a cell-based artificial piRNA production system using a silkworm ovarian cultured cell line that is believed to possess a complete piRNA pathway. In addition, we took advantage of a unique silkworm sex-determining one-to-one ping-pong piRNA pair, which enabled us to precisely monitor the behavior of individual artificial piRNAs. With this novel strategy, we successfully generated artificial piRNAs against endogenous protein-coding genes via the expected back-and-forth traveling mechanism. Furthermore, we detected "primary" piRNAs from the upstream region of the artificial "ping-pong" site in the endogenous gene. This artificial piRNA production system experimentally confirms the existence of the "ping-pong" cascade of piRNAs. Also, this system will enable us to identify the factors involved in both, or each, of the "ping" and "pong" cascades and the sequence features that are required for efficient piRNA production.

Detection of piRNAs in whitespotted bamboo shark liver

Piwi-interacting RNAs (piRNAs) are 26 to 31-nt small non-coding RNAs that have been reported mostly in germ-line cells and cancer cells. However, the presence of piRNAs in the whitespotted bamboo shark liver has not yet been reported. In a previous study of microRNAs in shark liver, some piRNAs were detected from small RNAs sequenced by Solexa technology. A total of 4857 piRNAs were predicted and found in shark liver. We further selected 17 piRNAs with high and significantly differential expression between normal and regenerative liver tissues for subsequent verification by Northern blotting. Ten piRNAs were further identified, and six of these were matched to known piRNAs in piRNABank. The actual expression of six known and four novel piRNAs was validated by qRT-PCR. In addition, a total of 401 target genes of the 10 piRNAs were predicted by miRanda. Through GO and pathway function analyses, only five piRNAs could be annotated with eighteen GO annotations. The results indicated that the identified piRNAs are involved in many important biological responses, including immune inflammation, cell-specific differentiation and development, and angiogenesis. This manuscript provides the first identification of piRNAs in the liver of whitespotted bamboo shark using Solexa technology as well as further elucidation of the regulatory role of piRNAs in whitespotted bamboo shark liver. These findings may provide a useful resource and may facilitate the development of therapeutic strategies against liver damage.

Expression dynamics, relationships, and transcriptional regulations of diverse transcripts in mouse spermatogenic cells

Among all tissues of the metazoa, the transcritpome of testis displays the highest diversity and specificity. However, its composition and dynamics during spermatogenesis have not been fully understood. Here, we have identified 20,639 message RNAs (mRNAs), 7,168 long non-coding RNAs (lncRNAs) and 15,101 circular RNAs (circRNAs) in mouse spermatogenic cells, and found many of them were specifically expressed in testes. lncRNAs are significantly more testis-specific than mRNAs. At all stages, mRNAs are generally more abundant than lncRNAs, and linear transcripts are more abundant than circRNAs. We showed that the productions of circRNAs and piRNAs were highly regulated instead of random processes. Based on the results of a small-scale functional screening experiment using cultured mouse spermatogonial stem cells, many evolutionarily conserved lncRNAs are likely to play roles in spermatogenesis. Typical classes of transcription factor binding sites are enriched in the promoters of testis-specific m/lncRNA genes. Target genes of CREM and RFX2, 2 key TFs for spermatogenesis, were further validated by using ChIP-chip assays and RNA-seq on RFX2-knockout spermatogenic cells. Our results contribute to the current understanding of the transcriptomic complexity of spermatogenic cells and provide a valuable resource from which many candidate genes may be selected for further functional studies.

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.

The Landscape of long noncoding RNA classification

Advances in the depth and quality of transcriptome sequencing have revealed many new classes of long noncoding RNAs (lncRNAs). lncRNA classification has mushroomed to accommodate these new findings, even though the real dimensions and complexity of the noncoding transcriptome remain unknown. Although evidence of functionality of specific lncRNAs continues to accumulate, conflicting, confusing, and overlapping terminology has fostered ambiguity and lack of clarity in the field in general. The lack of fundamental conceptual unambiguous classification framework results in a number of challenges in the annotation and interpretation of noncoding transcriptome data. It also might undermine integration of the new genomic methods and datasets in an effort to unravel the function of lncRNA. Here, we review existing lncRNA classifications, nomenclature, and terminology. Then, we describe the conceptual guidelines that have emerged for their classification and functional annotation based on expanding and more comprehensive use of large systems biology-based datasets.

Reduced pachytene piRNAs and translation underlie spermiogenic arrest in Maelstrom mutant mice

Pachytene piRNAs are a class of Piwi-interacting small RNAs abundant in spermatids of the adult mouse testis. They are processed from piRNA primary transcripts by a poorly understood mechanism and, unlike fetal transposon-derived piRNAs, lack complementary targets in the spermatid transcriptome. We report that immunopurified complexes of a conserved piRNA pathway protein Maelstrom (MAEL) are enriched in MIWI (Piwi partner of pachytene piRNAs), Tudor-domain proteins and processing intermediates of pachytene piRNA primary transcripts. We provide evidence of functional significance of these complexes in Mael129 knockout mice that exhibit spermiogenic arrest with acrosome and flagellum malformation. Mael129-null mutant testes possess low levels of piRNAs derived from MAEL-associated piRNA precursors and exhibit reduced translation of numerous spermiogenic mRNAs including those encoding acrosome and flagellum proteins. These translation defects in haploid round spermatids are likely indirect, as neither MAEL nor piRNA precursors associate with polyribosomes, and they may arise from an imbalance between pachytene piRNAs and MIWI.