Targeted, High-Resolution RNA Sequencing of Non-coding Genomic Regions Associated With Neuropsychiatric Functions

Biochemical-free enrichment or depletion of RNA classes in real-time during direct RNA sequencing with RISER

The heterogeneous composition of cellular transcriptomes poses a major challenge for detecting weakly expressed RNA classes, as they can be obscured by abundant RNAs. Although biochemical protocols can enrich or deplete specified RNAs, they are time-consuming, expensive and can compromise RNA integrity. Here we introduce RISER, a biochemical-free technology for the real-time enrichment or depletion of RNA classes. RISER performs selective rejection of molecules during direct RNA sequencing by identifying RNA classes directly from nanopore signals with deep learning and communicating with the sequencing hardware in real time. By targeting the dominant messenger and mitochondrial RNA classes for depletion, RISER reduces their respective read counts by more than 85%, resulting in an increase in sequencing depth of 47% on average for long non-coding RNAs. We also apply RISER for the depletion of globin mRNA in whole blood, achieving a decrease in globin reads by more than 90% as well as an increase in non-globin reads by 16% on average. Furthermore, using a GPU or a CPU, RISER is faster than GPU-accelerated basecalling and mapping. RISER's modular and retrainable software and intuitive command-line interface allow easy adaptation to other RNA classes. RISER is available at https://github.com/comprna/riser .

Cryptic splicing of stathmin-2 and UNC13A mRNAs is a pathological hallmark of TDP-43-associated Alzheimer's disease

Nuclear clearance and cytoplasmic accumulations of the RNA-binding protein TDP-43 are pathological hallmarks in almost all patients with amyotrophic lateral sclerosis (ALS) and up to 50% of patients with frontotemporal dementia (FTD) and Alzheimer's disease. In Alzheimer's disease, TDP-43 pathology is predominantly observed in the limbic system and correlates with cognitive decline and reduced hippocampal volume. Disruption of nuclear TDP-43 function leads to abnormal RNA splicing and incorporation of erroneous cryptic exons in numerous transcripts including Stathmin-2 (STMN2, also known as SCG10) and UNC13A, recently reported in tissues from patients with ALS and FTD. Here, we identify both STMN2 and UNC13A cryptic exons in Alzheimer's disease patients, that correlate with TDP-43 pathology burden, but not with amyloid-β or tau deposits. We also demonstrate that processing of the STMN2 pre-mRNA is more sensitive to TDP-43 loss of function than UNC13A. In addition, full-length RNAs encoding STMN2 and UNC13A are suppressed in large RNA-seq datasets generated from Alzheimer's disease post-mortem brain tissue. Collectively, these results open exciting new avenues to use STMN2 and UNC13A as potential therapeutic targets in a broad range of neurodegenerative conditions with TDP-43 proteinopathy including Alzheimer's disease.

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.

A Kuhnian revolution in molecular biology: Most genes in complex organisms express regulatory RNAs

Thomas Kuhn described the progress of science as comprising occasional paradigm shifts separated by interludes of 'normal science'. The paradigm that has held sway since the inception of molecular biology is that genes (mainly) encode proteins. In parallel, theoreticians posited that mutation is random, inferred that most of the genome in complex organisms is non-functional, and asserted that somatic information is not communicated to the germline. However, many anomalies appeared, particularly in plants and animals: the strange genetic phenomena of paramutation and transvection; introns; repetitive sequences; a complex epigenome; lack of scaling of (protein-coding) genes and increase in 'noncoding' sequences with developmental complexity; genetic loci termed 'enhancers' that control spatiotemporal gene expression patterns during development; and a plethora of 'intergenic', overlapping, antisense and intronic transcripts. These observations suggest that the original conception of genetic information was deficient and that most genes in complex organisms specify regulatory RNAs, some of which convey intergenerational information. Also see the video abstract here: https://youtu.be/qxeGwahBANw.

Transcript capture and ultradeep long-read RNA sequencing (CAPLRseq) to diagnose HNPCC/Lynch syndrome

PURPOSE

Whereas most human genes encode multiple mRNA isoforms with distinct function, clinical workflows for assessing this heterogeneity are not readily available. This is a substantial shortcoming, considering that up to 25% of disease-causing gene variants are suspected of disrupting mRNA splicing or mRNA abundance. Long-read sequencing can readily portray mRNA isoform diversity, but its sensitivity is relatively low due to insufficient transcriptome penetration.

METHODS

We developed and applied capture-based target enrichment from patient RNA samples combined with Oxford Nanopore long-read sequencing for the analysis of 123 hereditary cancer transcripts (capture and ultradeep long-read RNA sequencing (CAPLRseq)).

RESULTS

Validating CAPLRseq, we confirmed 17 cases of hereditary non-polyposis colorectal cancer/Lynch syndrome based on the demonstration of splicing defects and loss of allele expression of mismatch repair genes MLH1, PMS2, MSH2 and MSH6. Using CAPLRseq, we reclassified two variants of uncertain significance in MSH6 and PMS2 as either likely pathogenic or benign.

CONCLUSION

Our data show that CAPLRseq is an automatable and adaptable workflow for effective transcriptome-based identification of disease variants in a clinical diagnostic setting.

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.

RNA out of the mist

RNA has long been regarded primarily as the intermediate between genes and proteins. It was a surprise then to discover that eukaryotic genes are mosaics of mRNA sequences interrupted by large tracts of transcribed but untranslated sequences, and that multicellular organisms also express many long 'intergenic' and antisense noncoding RNAs (lncRNAs). The identification of small RNAs that regulate mRNA translation and half-life did not disturb the prevailing view that animals and plant genomes are full of evolutionary debris and that their development is mainly supervised by transcription factors. Gathering evidence to the contrary involved addressing the low conservation, expression, and genetic visibility of lncRNAs, demonstrating their cell-specific roles in cell and developmental biology, and their association with chromatin-modifying complexes and phase-separated domains. The emerging picture is that most lncRNAs are the products of genetic loci termed 'enhancers', which marshal generic effector proteins to their sites of action to control cell fate decisions during development.

Long non-coding RNAs associated with infection and vaccine-induced immunity

The immune system responds to infection or vaccination through a dynamic and complex process that involves several molecular and cellular factors. Among these factors, long non-coding RNAs (lncRNAs) have emerged as significant players in all areas of biology, particularly in immunology. Most of the mammalian genome is transcribed in a highly regulated manner, generating a diversity of lncRNAs that impact the differentiation and activation of immune cells and affect innate and adaptive immunity. Here, we have reviewed the range of functions and mechanisms of lncRNAs in response to infectious disease, including pathogen recognition, interferon (IFN) response, and inflammation. We describe examples of lncRNAs exploited by pathogenic agents during infection, which indicate that lncRNAs are a fundamental part of the arms race between hosts and pathogens. We also discuss lncRNAs potentially implicated in vaccine-induced immunity and present examples of lncRNAs associated with the antibody response of subjects receiving Influenza or Yellow Fever vaccines. Elucidating the widespread involvement of lncRNAs in the immune system will improve our understanding of the factors affecting immune response to different pathogenic agents, to better prevent and treat disease.

Isoform Age - Splice Isoform Profiling Using Long-Read Technologies

Alternative splicing (AS) of RNA is a key mechanism that results in the expression of multiple transcript isoforms from single genes and leads to an increase in the complexity of both the transcriptome and proteome. Regulation of AS is critical for the correct functioning of many biological pathways, while disruption of AS can be directly pathogenic in diseases such as cancer or cause risk for complex disorders. Current short-read sequencing technologies achieve high read depth but are limited in their ability to resolve complex isoforms. In this review we examine how long-read sequencing (LRS) technologies can address this challenge by covering the entire RNA sequence in a single read and thereby distinguish isoform changes that could impact RNA regulation or protein function. Coupling LRS with technologies such as single cell sequencing, targeted sequencing and spatial transcriptomics is producing a rapidly expanding suite of technological approaches to profile alternative splicing at the isoform level with unprecedented detail. In addition, integrating LRS with genotype now allows the impact of genetic variation on isoform expression to be determined. Recent results demonstrate the potential of these techniques to elucidate the landscape of splicing, including in tissues such as the brain where AS is particularly prevalent. Finally, we also discuss how AS can impact protein function, potentially leading to novel therapeutic targets for a range of diseases.

Neuroepigenetics of psychiatric disorders: Focus on lncRNA

Understanding the pathology of psychiatric disorders is challenging due to their complexity and multifactorial origin. However, development of high-throughput technologies has allowed for better insight into their molecular signatures. Advancement of sequencing methodologies have made it possible to study not only the protein-coding but also the noncoding genome. It is now clear that besides the genetic component, different epigenetic mechanisms play major roles in the onset and development of psychiatric disorders. Among them, examining the role of long noncoding RNAs (lncRNAs) is a relatively new field. Here, we present an overview of what is currently known about the involvement of lncRNAs in schizophrenia, major depressive and bipolar disorders, as well as suicide. The diagnosis of psychiatric disorders mainly relies on clinical evaluation without using measurable biomarkers. In this regard, lncRNA may open new opportunities for development of molecular tests. However, so far only a small set of known lncRNAs have been characterized at molecular level, which means they have a long way to go before clinical implementation. Understanding how changes in lncRNAs affect the appearance and development of psychiatric disorders may lead to a more classified and objective diagnostic system, but also open up new therapeutic targets for these patients.

Single-molecule, full-length transcript isoform sequencing reveals disease-associated RNA isoforms in cardiomyocytes

Alternative splicing generates differing RNA isoforms that govern phenotypic complexity of eukaryotes. Its malfunction underlies many diseases, including cancer and cardiovascular diseases. Comparative analysis of RNA isoforms at the genome-wide scale has been difficult. Here, we establish an experimental and computational pipeline that performs de novo transcript annotation and accurately quantifies transcript isoforms from cDNA sequences with a full-length isoform detection accuracy of 97.6%. We generate a searchable, quantitative human transcriptome annotation with 31,025 known and 5,740 novel transcript isoforms ( http://steinmetzlab.embl.de/iBrowser/ ). By analyzing the isoforms in the presence of RNA Binding Motif Protein 20 (RBM20) mutations associated with aggressive dilated cardiomyopathy (DCM), we identify 121 differentially expressed transcript isoforms in 107 cardiac genes. Our approach enables quantitative dissection of complex transcript architecture instead of mere identification of inclusion or exclusion of individual exons, as exemplified by the discovery of IMMT isoforms mis-spliced by RBM20 mutations. Thereby we achieve a path to direct differential expression testing independent of an existing annotation of transcript isoforms, providing more immediate biological interpretation and higher resolution transcriptome comparisons.

Large-Scale Transcriptomics-Driven Approach Revealed Overexpression of CRNDE as a Poor Survival Prognosis Biomarker in Glioblastoma

Glioblastoma is the most common and malignant brain malignancy worldwide, with a 10-year survival of only 0.7%. Aggressive multimodal treatment is not enough to increase life expectancy and provide good quality of life for glioblastoma patients. In addition, despite decades of research, there are no established biomarkers for early disease diagnosis and monitoring of patient response to treatment. High throughput sequencing technologies allow for the identification of unique molecules from large clinically annotated datasets. Thus, the aim of our study was to identify significant molecular changes between short- and long-term glioblastoma survivors by transcriptome RNA sequencing profiling, followed by differential pathway-activation-level analysis. We used data from the publicly available repositories The Cancer Genome Atlas (TCGA; number of annotated cases = 135) and Chinese Glioma Genome Atlas (CGGA; number of annotated cases = 218), and experimental clinically annotated glioblastoma tissue samples from the Institute of Pathology, Faculty of Medicine in Ljubljana corresponding to 2-58 months overall survival (n = 16). We found one differential gene for long noncoding RNA CRNDE whose overexpression showed correlation to poor patient OS. Moreover, we identified overlapping sets of congruently regulated differential genes involved in cell growth, division, and migration, structure and dynamics of extracellular matrix, DNA methylation, and regulation through noncoding RNAs. Gene ontology analysis can provide additional information about the function of protein- and nonprotein-coding genes of interest and the processes in which they are involved. In the future, this can shape the design of more targeted therapeutic approaches.

Long-Read Sequencing Improves the Detection of Structural Variations Impacting Complex Non-Coding Elements of the Genome

The advent of long-read sequencing offers a new assessment method of detecting genomic structural variation (SV) in numerous rare genetic diseases. For autism spectrum disorders (ASD) cases where pathogenic variants fail to be found in the protein-coding genic regions along chromosomes, we proposed a scalable workflow to characterize the risk factor of SVs impacting non-coding elements of the genome. We applied whole-genome sequencing on an Emirati family having three children with ASD using long and short-read sequencing technology. A series of analytical pipelines were established to identify a set of SVs with high sensitivity and specificity. At 15-fold coverage, we observed that long-read sequencing technology (987 variants) detected a significantly higher number of SVs when compared to variants detected using short-read technology (509 variants) (p-value < 1.1020 × 10-57). Further comparison showed 97.9% of long-read sequencing variants were spanning within the 1-100 kb size range (p-value < 9.080 × 10-67) and impacting over 5000 genes. Moreover, long-read variants detected 604 non-coding RNAs (p-value < 9.02 × 10-9), comprising 58% microRNA, 31.9% lncRNA, and 9.1% snoRNA. Even at low coverage, long-read sequencing has shown to be a reliable technology in detecting SVs impacting complex elements of the genome.

On the functional relevance of spatiotemporally-specific patterns of experience-dependent long noncoding RNA expression in the brain

The majority of transcriptionally active RNA derived from the mammalian genome does not code for protein. Long noncoding RNA (lncRNA) is the most abundant form of noncoding RNA found in the brain and is involved in many aspects of cellular metabolism. Beyond their fundamental role in the nucleus as decoys for RNA-binding proteins associated with alternative splicing or as guides for the epigenetic regulation of protein-coding gene expression, recent findings indicate that activity-induced lncRNAs also regulate neural plasticity. In this review, we discuss how lncRNAs may exert molecular control over brain function beyond their known roles in the nucleus. We propose that subcellular localization is a critical feature of experience-dependent lncRNA activity in the brain, and that lncRNA-mediated control over RNA metabolism at the synapse serves to regulate local mRNA stability and translation, thereby influencing neuronal function, learning and memory.

Preparation of Mammalian Nascent RNA for Long Read Sequencing

Long read sequencing technologies now allow high-quality sequencing of RNAs (or their cDNAs) that are hundreds to thousands of nucleotides long. Long read sequences of nascent RNA provide single-nucleotide-resolution information about co-transcriptional RNA processing events-e.g., splicing, folding, and base modifications. Here, we describe how to isolate nascent RNA from mammalian cells through subcellular fractionation of chromatin-associated RNA, as well as how to deplete poly(A)+ RNA and rRNA, and, finally, how to generate a full-length cDNA library for use on long read sequencing platforms. This approach allows for an understanding of coordinated splicing status across multi-intron transcripts by revealing patterns of splicing or other RNA processing events that cannot be gained from traditional short read RNA sequencing. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Subcellular fractionation Basic Protocol 2: Nascent RNA isolation and adapter ligation Basic Protocol 3: cDNA amplicon preparation.

Improving the diversity of captured full-length isoforms using a normalized single-molecule RNA-sequencing method

Human genes form a large variety of isoforms after transcription, encoding distinct transcripts to exert different functions. Single-molecule RNA sequencing facilitates accurate identification of the isoforms by extending nucleotide read length significantly. However, the gene or isoform diversity is lowly represented by the mRNA molecules captured by single-molecule RNA sequencing. Here, we show that a cDNA normalization procedure before the library preparation for PacBio RS II sequencing captures 3.2–6.0 fold more full-length high-quality isoform species for different human samples, as compared to the non-normalized capture procedure. Many lowly expressed, functionally important isoforms can be detected. In addition, normalized PacBio RNA sequencing also resolves more allele-specific haplotype transcripts. Finally, we apply the cDNA normalization based long-read RNA sequencing method to profile the transcriptome of human gastric signet-ring cell carcinomas, identify new cancer-specific transcriptome signatures, and thus, bring out the utility of the improved protocols in gene expression studies.

Hu et al. combine cDNA normalization before library preparation with a software tool to increase the capture of RNA isoform species in single-molecule RNA sequencing. They demonstrate that this approach can detect previously unknown transcripts in gastric signet-ring cell carcinomas that are not present in non-malignant tissue.

Functional annotation of human long noncoding RNAs via molecular phenotyping

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

Nanopore Sequencing in Blood Diseases: A Wide Range of Opportunities

The molecular pathogenesis of hematological diseases is often driven by genetic and epigenetic alterations. Next-generation sequencing has considerably increased our genomic knowledge of these disorders becoming ever more widespread in clinical practice. In 2012 Oxford Nanopore Technologies (ONT) released the MinION, the first long-read nanopore-based sequencer, overcoming the main limits of short-reads sequences generation. In the last years, several nanopore sequencing approaches have been performed in various "-omic" sciences; this review focuses on the challenge to introduce ONT devices in the hematological field, showing advantages, disadvantages and future perspectives of this technology in the precision medicine era.

Non-Coding RNAs in Psychiatric Disorders and Suicidal Behavior

It is well known that only a small proportion of the human genome code for proteins; the rest belong to the family of RNAs that do not code for protein and are known as non-coding RNAs (ncRNAs). ncRNAs are further divided into two subclasses based on size: 1) long non-coding RNAs (lncRNAs; >200 nucleotides) and 2) small RNAs (<200 nucleotides). Small RNAs contain various family members that include microRNAs (miRNAs), small interfering RNAs (siRNAs), piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs), and small nuclear RNAs (snRNAs). The roles of ncRNAs, especially lncRNAs and miRNAs, are well documented in brain development, homeostasis, stress responses, and neural plasticity. It has also been reported that ncRNAs can influence the development of psychiatric disorders including schizophrenia, major depressive disorder, and bipolar disorder. More recently, their roles are being investigated in suicidal behavior. In this article, we have comprehensively reviewed the findings of lncRNA and miRNA expression changes and their functions in various psychiatric disorders including suicidal behavior. We primarily focused on studies that have been done in postmortem human brain. In addition, we have briefly reviewed the role of other small RNAs (e.g. piwiRNA, siRNA, snRNA, and snoRNAs) and their expression changes in psychiatric illnesses.

Getting the Entire Message: Progress in Isoform Sequencing

The advent of second-generation sequencing and its application to RNA sequencing have revolutionized the field of genomics by allowing quantification of gene expression, as well as the definition of transcription start/end sites, exons, splice sites and RNA editing sites. However, due to the sequencing of fragments of cDNAs, these methods have not given a reliable picture of complete RNA isoforms. Third-generation sequencing has filled this gap and allows end-to-end sequencing of entire RNA/cDNA molecules. This approach to transcriptomics has been a "niche" technology for a couple of years but now is becoming mainstream with many different applications. Here, we review the background and progress made to date in this rapidly growing field. We start by reviewing the progressive realization that alternative splicing is omnipresent. We then focus on long-noncoding RNA isoforms and the distinct combination patterns of exons in noncoding and coding genes. We consider the implications of the recent technologies of direct RNA sequencing and single-cell isoform RNA sequencing. Finally, we discuss the parameters that define the success of long-read RNA sequencing experiments and strategies commonly used to make the most of such data.

TGStools: A Bioinformatics Suit to Facilitate Transcriptome Analysis of Long Reads from Third Generation Sequencing Platform

Recent analyses show that transcriptome sequencing can be utilized as a diagnostic tool for rare Mendelian diseases. The third generation sequencing de novo detects long reads of thousands of base pairs, thus greatly expanding the isoform discovery and identification of novel long noncoding RNAs. In this study, we developed TGStools, a bioinformatics suite to facilitate routine tasks such as characterizing full-length transcripts, detecting shifted types of alternative splicing, and long noncoding RNAs (lncRNAs) identification in transcriptome analysis. It also prioritizes the transcripts with a visualization framework that automatically integrates rich annotation with known genomic features. TGStools is a Python package freely available at Github.