Evaluating Performance of Different RNA Secondary Structure Prediction Programs Using Self-cleaving Ribozymes.

Accurate identification of the correct, biologically relevant RNA structures is critical to understanding various aspects of RNA biology since proper folding represents the key to the functionality of all types of RNA molecules and plays pivotal roles in many essential biological processes. Thus, a plethora of approaches have been developed to predict, identify, or solve RNA structures based on various computational, molecular, genetic, chemical, or physicochemical strategies. Purely computational approaches hold distinct advantages over all other strategies in terms of the ease of implementation, time, speed, cost, and throughput, but they strongly underperform in terms of accuracy that significantly limits their broader application. Nonetheless, the advantages of these methods led to a steady development of multiple in silico RNA secondary structure prediction approaches including recent deep learning-based programs. Here, we compared the accuracy of predictions of biologically relevant secondary structures of dozens of self-cleaving ribozyme sequences using seven in silico RNA folding prediction tools with tasks of varying complexity. We found that while many programs performed well in relatively simple tasks, their performance varied significantly in more complex RNA folding problems. However, in general, a modern deep learning method outperformed the other programs in the complex tasks in predicting the RNA secondary structures, at least based on the specific class of sequences tested, suggesting that it may represent the future of RNA structure prediction algorithms.

Common occurrence of hotspots of single strand DNA breaks at transcriptional start sites.

BACKGROUND: We recently developed two high-resolution methods for genome-wide mapping of two prominent types of DNA damage, single-strand DNA breaks (SSBs) and abasic (AP) sites and found highly complex and non-random patterns of these lesions in mammalian genomes. One salient feature of SSB and AP sites was the existence of single-nucleotide hotspots for both lesions. RESULTS: In this work, we show that SSB hotspots are enriched in the immediate vicinity of transcriptional start sites (TSSs) in multiple normal mammalian tissues, however the magnitude of enrichment varies significantly with tissue type and appears to be limited to a subset of genes. SSB hotspots around TSSs are enriched on the template strand and associate with higher expression of the corresponding genes. Interestingly, SSB hotspots appear to be at least in part generated by the base-excision repair (BER) pathway from the AP sites. CONCLUSIONS: Our results highlight complex relationship between DNA damage and regulation of gene expression and suggest an exciting possibility that SSBs at TSSs might function as sensors of DNA damage to activate genes important for DNA damage response.

Genome-wide profiles of DNA damage represent highly accurate predictors of mammalian age.

The identification of novel age-related biomarkers represents an area of intense research interest. Despite multiple studies associating DNA damage with aging, there is a glaring paucity of DNA damage-based biomarkers of age, mainly due to the lack of precise methods for genome-wide surveys of different types of DNA damage. Recently, we developed two techniques for genome-wide mapping of the most prevalent types of DNA damage, single-strand breaks and abasic sites, with nucleotide-level resolution. Herein, we explored the potential of genomic patterns of DNA damage identified by these methods as a source of novel age-related biomarkers using mice as a model system. Strikingly, we found that models based on genomic patterns of either DNA lesion could accurately predict age with higher precision than the commonly used transcriptome analysis. Interestingly, the informative patterns were limited to relatively few genes and the DNA damage levels were positively or negatively correlated with age. These findings show that previously unexplored high-resolution genomic patterns of DNA damage contain useful information that can contribute significantly to both practical applications and basic science.

Evidence for widespread existence of functional novel and non-canonical human transcripts.

BACKGROUND: Fraction of functional sequence in the human genome remains a key unresolved question in Biology and the subject of vigorous debate. While a plethora of studies have connected a significant fraction of human DNA to various biochemical processes, the classical definition of function requires evidence of effects on cellular or organismal fitness that such studies do not provide. Although multiple high-throughput reverse genetics screens have been developed to address this issue, they are limited to annotated genomic elements and suffer from non-specific effects, arguing for a strong need to develop additional functional genomics approaches. RESULTS: In this work, we established a high-throughput lentivirus-based insertional mutagenesis strategy as a forward genetics screen tool in aneuploid cells. Application of this approach to human cell lines in multiple phenotypic screens suggested the presence of many yet uncharacterized functional elements in the human genome, represented at least in part by novel exons of known and novel genes. The novel transcripts containing these exons can be massively, up to thousands-fold, induced by specific stresses, and at least some can represent bi-cistronic protein-coding mRNAs. CONCLUSIONS: Altogether, these results argue that many unannotated and non-canonical human transcripts, including those that appear as aberrant splice products, have biological relevance under specific biological conditions.

Genome-Wide Profiling of Endogenous Single-Stranded DNA Using the SSiNGLe-P1 Method.

Endogenous single-stranded DNA (essDNA) can form in a mammalian genome as the result of a variety of molecular processes and can both play important roles inside the cell as well as have detrimental consequences to genome integrity, much of which remains to be fully understood. Here, we established the SSiNGLe-P1 approach based on limited digestion by P1 endonuclease for high-throughput genome-wide identification of essDNA regions. We applied this method to profile essDNA in both human mitochondrial and nuclear genomes. In the mitochondrial genome, the profiles of essDNA provide new evidence to support the strand-displacement model of mitochondrial DNA replication. In the nuclear genome, essDNA regions were found to be enriched in certain types of functional genomic elements, particularly, the origins of DNA replication, R-loops, and to a lesser degree, in promoters. Furthermore, interestingly, many of the essDNA regions identified by SSiNGLe-P1 have not been annotated and thus could represent yet unknown functional elements.

Lessons from discovery of true ADAR RNA editing sites in a human cell line.

BACKGROUND: Conversion or editing of adenosine (A) into inosine (I) catalyzed by specialized cellular enzymes represents one of the most common post-transcriptional RNA modifications with emerging connection to disease. A-to-I conversions can happen at specific sites and lead to increase in proteome diversity and changes in RNA stability, splicing, and regulation. Such sites can be detected as adenine-to-guanine sequence changes by next-generation RNA sequencing which resulted in millions reported sites from multiple genome-wide surveys. Nonetheless, the lack of extensive independent validation in such endeavors, which is critical considering the relatively high error rate of next-generation sequencing, leads to lingering questions about the validity of the current compendiums of the editing sites and conclusions based on them. RESULTS: Strikingly, we found that the current analytical methods suffer from very high false positive rates and that a significant fraction of sites in the public databases cannot be validated. In this work, we present potential solutions to these problems and provide a comprehensive and extensively validated list of A-to-I editing sites in a human cancer cell line. Our findings demonstrate that most of true A-to-I editing sites in a human cancer cell line are located in the non-coding transcripts, the so-called RNA 'dark matter'. On the other hand, many ADAR editing events occurring in exons of human protein-coding mRNAs, including those that can recode the transcriptome, represent false positives and need to be interpreted with caution. Nonetheless, yet undiscovered authentic ADAR sites that increase the diversity of human proteome exist and warrant further identification. CONCLUSIONS: Accurate identification of human ADAR sites remains a challenging problem, particularly for the sites in exons of protein-coding mRNAs. As a result, genome-wide surveys of ADAR editome must still be accompanied by extensive Sanger validation efforts. However, given the vast number of unknown human ADAR sites, there is a need for further developments of the analytical techniques, potentially those that are based on deep learning solutions, in order to provide a quick and reliable identification of the editome in any sample.

Transcriptomics for Clinical and Experimental Biology Research: Hang on a Seq.

Sequencing the human genome empowers translational medicine, facilitating transcriptome-wide molecular diagnosis, pathway biology, and drug repositioning. Initially, microarrays are used to study the bulk transcriptome; but now short-read RNA sequencing (RNA-seq) predominates. Positioned as a superior technology, that makes the discovery of novel transcripts routine, most RNA-seq analyses are in fact modeled on the known transcriptome. Limitations of the RNA-seq methodology have emerged, while the design of, and the analysis strategies applied to, arrays have matured. An equitable comparison between these technologies is provided, highlighting advantages that modern arrays hold over RNA-seq. Array protocols more accurately quantify constitutively expressed protein coding genes across tissue replicates, and are more reliable for studying lower expressed genes. Arrays reveal long noncoding RNAs (lncRNA) are neither sparsely nor lower expressed than protein coding genes. Heterogeneous coverage of constitutively expressed genes observed with RNA-seq, undermines the validity and reproducibility of pathway analyses. The factors driving these observations, many of which are relevant to long-read or single-cell sequencing are discussed. As proposed herein, a reappreciation of bulk transcriptomic methods is required, including wider use of the modern high-density array data-to urgently revise existing anatomical RNA reference atlases and assist with more accurate study of lncRNAs.

Dynamic Patterns of Mammalian Mitochondrial DNA Replication Uncovered Using SSiNGLe-5'ES.

The proper replication of mitochondrial DNA is key to the maintenance of this crucial organelle. Multiple studies aimed at understanding the mechanisms of replication of the mitochondrial genome have been conducted in the past several decades; however, while highly informative, they were conducted using relatively low-sensitivity techniques. Here, we established a high-throughput approach based on next-generation sequencing to identify replication start sites with nucleotide-level resolution and applied it to the genome of mitochondria from different human and mouse cell types. We found complex and highly reproducible patterns of mitochondrial initiation sites, both previously annotated and newly discovered in this work, that showed differences among different cell types and species. These results suggest that the patterns of the replication initiation sites are dynamic and might reflect, in some yet unknown ways, the complexities of mitochondrial and cellular physiology. Overall, this work suggests that much remains unknown about the details of mitochondrial DNA replication in different biological states, and the method established here opens up a new avenue in the study of the replication of mitochondrial and potentially other genomes.

Traditional Chinese Medicine Treatment, Gua Sha, can Induce Subtle Molecular Changes in Gene Expression.

Objective: Here, we explored molecular changes that could potentially mediate healing effects of Gua Sha - a method employed by the Chinese traditional medicine with proven track records of safe and efficient applications dating back to ancient times as well as support from randomized controlled trials performed by modern medical studies - yet remaining almost entirely unexplored by the modern-day high-throughput methods of the -omics sciences. Methods: We investigated transcriptome changes occurring shortly after Gua Sha treatment in the whole blood of healthy volunteers using bulk RNA-seq analysis. We applied various analytical tools to identify genes with consistent expression changes in multiple individuals in response to Gua Sha and their networks. Results: We found that while the changes were very subtle and individual-specific, we could identify consistent upregulation of three histone genes. Further analysis of the potential regulatory networks of these histone genes revealed the enrichment of functions involved in the immune response and inflammation. Conclusion: The significance of these results in the context of potential effects of Gua Sha and the next steps in exploring the molecular mechanisms of action of this technique are discussed.

Evidence for Existence of Multiple Functional Human Small RNAs Derived from Transcripts of Protein-Coding Genes.

The human genome encodes a multitude of different noncoding transcripts that have been traditionally separated on the basis of their lengths into long (>200 nt) or small (<200 nt) noncoding RNAs. The functions, mechanisms of action, and biological relevance of the vast majority of both long and short noncoding transcripts remain unknown. However, according to the functional understanding of the known classes of long and small noncoding RNAs (sncRNAs) that have been shown to play crucial roles in multiple biological processes, it is generally assumed that many unannotated long and small transcripts participate in important cellular functions as well. Nevertheless, direct evidence of functionality is lacking for most noncoding transcripts, especially for sncRNAs that are often dismissed as stable degradation products of longer RNAs. Here, we developed a high-throughput assay to test the functionality of sncRNAs by overexpressing them in human cells. Surprisingly, we found that a significant fraction (>40%) of unannotated sncRNAs appear to have biological relevance. Furthermore, contrary to the expectation, the potentially functional transcripts are not highly abundant and can be derived from protein-coding mRNAs. These results strongly suggest that the small noncoding transcriptome can harbor multiple functional transcripts that warrant future studies.

Complexities of mammalian transcriptome revealed by targeted RNA enrichment techniques.

Studies using highly sensitive targeted RNA enrichment methods have shown that a large portion of the human transcriptome remains to be discovered and that most of the genome is transcribed in a complex, interleaved fashion characterized by a complex web of transcripts emanating from protein coding and noncoding loci. These results resonate with those from single-cell transcriptome profiling endeavors that reveal the existence of multiple novel, cell type-specific transcripts and clearly demonstrate that our understanding of the complexities of the human transcriptome is far from being complete. Here, we review the current status of the targeted RNA enrichment techniques, their application to the discovery of novel cell type-specific transcripts, and their impact on our understanding of the human genome and transcriptome.

Long non-coding RNAs: definitions, functions, challenges and recommendations.

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Complex genomic patterns of abasic sites in mammalian DNA revealed by a high-resolution SSiNGLe-AP method.

DNA damage plays a critical role in biology and diseases; however, how different types of DNA lesions affect cellular functions is far from clear mostly due to the paucity of high-resolution methods that can map their locations in complex genomes, such as those of mammals. Here, we present the development and validation of SSiNGLe-AP method, which can map a common type of DNA damage, abasic (AP) sites, in a genome-wide and high-resolution manner. We apply this method to six different tissues of mice with different ages and human cancer cell lines. We find a nonrandom distribution of AP sites in the mammalian genome that exhibits dynamic enrichment at specific genomic locations, including single-nucleotide hotspots, and is significantly influenced by gene expression, age and tissue type in particular. Overall, these results suggest that we are only starting to understand the true complexities in the genomic patterns of DNA damage.

Hotspots of single-strand DNA "breakome" are enriched at transcriptional start sites of genes.

Single-strand breaks (SSBs) represent one of the most common types of DNA damage, yet not much is known about the genome landscapes of this type of DNA lesions in mammalian cells. Here, we found that SSBs are more likely to occur in certain positions of the human genome-SSB hotspots-in different cells of the same cell type and in different cell types. We hypothesize that the hotspots are likely to represent biologically relevant breaks. Furthermore, we found that the hotspots had a prominent tendency to be enriched in the immediate vicinity of transcriptional start sites (TSSs). We show that these hotspots are not likely to represent technical artifacts or be caused by common mechanisms previously found to cause DNA cleavage at promoters, such as apoptotic DNA fragmentation or topoisomerase type II (TOP2) activity. Therefore, such TSS-associated hotspots could potentially be generated using a novel mechanism that could involve preferential cleavage at cytosines, and their existence is consistent with recent studies suggesting a complex relationship between DNA damage and regulation of gene expression.

Methods to Analyze the Non-Coding RNA Interactome-Recent Advances and Challenges.

Most of the human genome is transcribed to generate a multitude of non-coding RNAs. However, while these transcripts have generated an immense amount of scientific interest, their biological function remains a subject of an intense debate. Understanding mechanisms of action of non-coding RNAs is a key to addressing the issue of biological relevance of these transcripts. Based on some well-understood non-coding RNAs that function inside the cell by interacting with other molecules, it is generally believed many other non-coding transcripts could also function in a similar fashion. Therefore, development of methods that can map RNA interactome is the key to understanding functionality of the extensive cellular non-coding transcriptome. Here, we review the vast progress that has been made in the past decade in technologies that can map RNA interactions with different sites in DNA, proteins or other RNA molecules; the general approaches used to validate the existence of novel interactions; and the challenges posed by interpreting the data obtained using the interactome mapping methods.

Complex Age- and Cancer-Related Changes in Human Blood Transcriptome-Implications for Pan-Cancer Diagnostics.

Early cancer detection is the key to a positive clinical outcome. While a number of early diagnostics methods exist in clinics today, they tend to be invasive and limited to a few cancer types. Thus, a clear need exists for non-invasive diagnostics methods that can be used to detect the presence of cancer of any type. Liquid biopsy based on analysis of molecular components of peripheral blood has shown significant promise in such pan-cancer diagnostics; however, existing methods based on this approach require improvements, especially in sensitivity of early-stage cancer detection. The improvement would likely require diagnostics assays based on multiple different types of biomarkers and, thus, calls for identification of novel types of cancer-related biomarkers that can be used in liquid biopsy. Whole-blood transcriptome, especially its non-coding component, represents an obvious yet under-explored biomarker for pan-cancer detection. In this study, we show that whole transcriptome analysis using RNA-seq could indeed serve as a viable biomarker for pan-cancer detection. Furthermore, a class of long non-coding (lnc) RNAs, very long intergenic non-coding (vlinc) RNAs, demonstrated superior performance compared with protein-coding mRNAs. Finally, we show that age and presence of non-blood cancers change transcriptome in similar, yet not identical, directions and explore implications of this observation for pan-cancer diagnostics.

Chromatin lncRNA Platr10 controls stem cell pluripotency by coordinating an intrachromosomal regulatory network.

BACKGROUND: A specific 3-dimensional intrachromosomal architecture of core stem cell factor genes is required to reprogram a somatic cell into pluripotency. As little is known about the epigenetic readers that orchestrate this architectural remodeling, we used a novel chromatin RNA in situ reverse transcription sequencing (CRIST-seq) approach to profile long noncoding RNAs (lncRNAs) in the Oct4 promoter. RESULTS: We identify Platr10 as an Oct4 - Sox2 binding lncRNA that is activated in somatic cell reprogramming. Platr10 is essential for the maintenance of pluripotency, and lack of this lncRNA causes stem cells to exit from pluripotency. In fibroblasts, ectopically expressed Platr10 functions in trans to activate core stem cell factor genes and enhance pluripotent reprogramming. Using RNA reverse transcription-associated trap sequencing (RAT-seq), we show that Platr10 interacts with multiple pluripotency-associated genes, including Oct4, Sox2, Klf4, and c-Myc, which have been extensively used to reprogram somatic cells. Mechanistically, we demonstrate that Platr10 helps orchestrate intrachromosomal promoter-enhancer looping and recruits TET1, the enzyme that actively induces DNA demethylation for the initiation of pluripotency. We further show that Platr10 contains an Oct4 binding element that interacts with the Oct4 promoter and a TET1-binding element that recruits TET1. Mutation of either of these two elements abolishes Platr10 activity. CONCLUSION: These data suggest that Platr10 functions as a novel chromatin RNA molecule to control pluripotency in trans by modulating chromatin architecture and regulating DNA methylation in the core stem cell factor network.

Very long intergenic non-coding (vlinc) RNAs directly regulate multiple genes in cis and trans.

BACKGROUND: The majority of the human genome is transcribed in the form of long non-coding (lnc) RNAs. While these transcripts have attracted considerable interest, their molecular mechanisms of function and biological significance remain controversial. One of the main reasons behind this lies in the significant challenges posed by lncRNAs requiring the development of novel methods and concepts to unravel their functionality. Existing methods often lack cross-validation and independent confirmation by different methodologies and therefore leave significant ambiguity as to the authenticity of the outcomes. Nonetheless, despite all the caveats, it appears that lncRNAs may function, at least in part, by regulating other genes via chromatin interactions. Therefore, the function of a lncRNA could be inferred from the function of genes it regulates. In this work, we present a genome-wide functional annotation strategy for lncRNAs based on identification of their regulatory networks via the integration of three distinct types of approaches: co-expression analysis, mapping of lncRNA-chromatin interactions, and assaying molecular effects of lncRNA knockdowns obtained using an inducible and highly specific CRISPR/Cas13 system. RESULTS: We applied the strategy to annotate 407 very long intergenic non-coding (vlinc) RNAs belonging to a novel widespread subclass of lncRNAs. We show that vlincRNAs indeed appear to regulate multiple genes encoding proteins predominantly involved in RNA- and development-related functions, cell cycle, and cellular adhesion via a mechanism involving proximity between vlincRNAs and their targets in the nucleus. A typical vlincRNAs can be both a positive and negative regulator and regulate multiple genes both in trans and cis. Finally, we show vlincRNAs and their regulatory networks potentially represent novel components of DNA damage response and are functionally important for the ability of cancer cells to survive genotoxic stress. CONCLUSIONS: This study provides strong evidence for the regulatory role of the vlincRNA class of lncRNAs and a potentially important role played by these transcripts in the hidden layer of RNA-based regulation in complex biological systems.

Hovlinc is a recently evolved class of ribozyme found in human lncRNA.

Although naturally occurring catalytic RNA molecules-ribozymes-have attracted a great deal of research interest, very few have been identified in humans. Here, we developed a genome-wide approach to discovering self-cleaving ribozymes and identified a naturally occurring ribozyme in humans. The secondary structure and biochemical properties of this ribozyme indicate that it belongs to an unidentified class of small, self-cleaving ribozymes. The sequence of the ribozyme exhibits a clear evolutionary path, from its appearance between ~130 and ~65 million years ago (Ma), to acquiring self-cleavage activity very recently, ~13-10 Ma, in the common ancestors of humans, chimpanzees and gorillas. The ribozyme appears to be functional in vivo and is embedded within a long noncoding RNA belonging to a class of very long intergenic noncoding RNAs. The presence of a catalytic RNA enzyme in lncRNA creates the possibility that these transcripts could function by carrying catalytic RNA domains.