Direct long-read RNA sequencing identifies a subset of questionable exitrons likely arising from reverse transcription artifacts

Sequencing accuracy and systematic errors of nanopore direct RNA sequencing

BACKGROUND

Direct RNA sequencing (dRNA-seq) on the Oxford Nanopore Technologies (ONT) platforms can produce reads covering up to full-length gene transcripts, while containing decipherable information about RNA base modifications and poly-A tail lengths. Although many published studies have been expanding the potential of dRNA-seq, its sequencing accuracy and error patterns remain understudied.

RESULTS

We present the first comprehensive evaluation of sequencing accuracy and characterisation of systematic errors in dRNA-seq data from diverse organisms and synthetic in vitro transcribed RNAs. We found that for sequencing kits SQK-RNA001 and SQK-RNA002, the median read accuracy ranged from 87% to 92% across species, and deletions significantly outnumbered mismatches and insertions. Due to their high abundance in the transcriptome, heteropolymers and short homopolymers were the major contributors to the overall sequencing errors. We also observed systematic biases across all species at the levels of single nucleotides and motifs. In general, cytosine/uracil-rich regions were more likely to be erroneous than guanines and adenines. By examining raw signal data, we identified the underlying signal-level features potentially associated with the error patterns and their dependency on sequence contexts. While read quality scores can be used to approximate error rates at base and read levels, failure to detect DNA adapters may be a source of errors and data loss. By comparing distinct basecallers, we reason that some sequencing errors are attributable to signal insufficiency rather than algorithmic (basecalling) artefacts. Lastly, we generated dRNA-seq data using the latest SQK-RNA004 sequencing kit released at the end of 2023 and found that although the overall read accuracy increased, the systematic errors remain largely identical compared to the previous kits.

CONCLUSIONS

As the first systematic investigation of dRNA-seq errors, this study offers a comprehensive overview of reproducible error patterns across diverse datasets, identifies potential signal-level insufficiency, and lays the foundation for error correction methods.

Exitron splicing of odor receptor genes in Drosophila

Proper expression of odor receptor genes is critical for the function of olfactory systems. In this study, we identified exitrons (exonic introns) in four of the 39 Odorant receptor (Or) genes expressed in the Drosophila antenna. Exitrons are sequences that can be spliced out from within a protein-coding exon, thereby altering the encoded protein. We focused on Or88a, which encodes a pheromone receptor, and found that exitron splicing of Or88a is conserved across five Drosophila species over 20 My of evolution. The exitron was spliced out in 15% of Or88a transcripts. Removal of this exitron creates a non-coding RNA rather than an RNA that encodes a stable protein. Our results suggest the hypothesis that in the case of Or88a, exitron splicing could act in neuronal modulation by decreasing the level of functional Or transcripts. Activation of Or88a-expressing olfactory receptor neurons via either optogenetics or pheromone stimulation increased the level of exitron-spliced transcripts, with optogenetic activation leading to a 14-fold increase. A fifth Or can also undergo an alternative splicing event that eliminates most of the canonical open reading frame. Besides these cases of alternative splicing, we found alternative polyadenylation of four Ors, and exposure of Or67c to its ligand ethyl lactate in the antenna downregulated all of its 3' isoforms. Our study reveals mechanisms by which neuronal activity could be modulated via regulation of the levels of Or isoforms.

Contrasting and Combining Transcriptome Complexity Captured by Short and Long RNA Sequencing Reads

Mapping transcriptomic variations using either short or long reads RNA sequencing is a staple of genomic research. Long reads are able to capture entire isoforms and overcome repetitive regions, while short reads still provides improved coverage and error rates. Yet how to quantitatively compare the technologies, can we combine those, and what may be the benefit of such a combined view remain open questions. We tackle these questions by first creating a pipeline to assess matched long and short reads data using a variety of transcriptome statistics. We find that across datasets, algorithms and technologies, matched short reads data detects roughly 50% more splice junctions, with 10-30% of the splice junctions included at 20% or more are missed by long reads. In contrast, long reads detect many more intron retention events, pointing to the benefit of combining the technologies. We introduce MAJIQ-L, an extension of the MAJIQ software to enable a unified view of transcriptome variations from both technologies and demonstrate its benefits. Our software can be used to assess any future long reads technology or algorithm, and combine it with short reads data for improved transcriptome analysis.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies

Aberrant skipping of coding exons in CD19 and CD22 compromises the response to immunotherapy in B-cell malignancies. Here, we showed that the MS4A1 gene encoding human CD20 also produces several messenger RNA (mRNA) isoforms with distinct 5' untranslated regions. Four variants (V1-4) were detected using RNA sequencing (RNA-seq) at distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma, only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform contained upstream open reading frames and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, whereas V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed chimeric antigen receptor T cells were able to kill both V3- and V1-expressing cells, but the bispecific T-cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on 4 postmosunetuzumab follicular lymphoma relapses and discovered that in 2 of them, the downregulation of CD20 was accompanied by a V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies.

Identification and quantification of small exon-containing isoforms in long-read RNA sequencing data

Small exons are pervasive in transcriptomes across organisms, and their quantification in RNA isoforms is crucial for understanding gene functions. Although long-read RNA-seq based on Oxford Nanopore Technologies (ONT) offers the advantage of covering transcripts in full length, its lower base accuracy poses challenges for identifying individual exons, particularly microexons (≤ 30 nucleotides). Here, we systematically assess small exons quantification in synthetic and human ONT RNA-seq datasets. We demonstrate that reads containing small exons are often not properly aligned, affecting the quantification of relevant transcripts. Thus, we develop a local-realignment method for misaligned exons (MisER), which remaps reads with misaligned exons to the transcript references. Using synthetic and simulated datasets, we demonstrate the high sensitivity and specificity of MisER for the quantification of transcripts containing small exons. Moreover, MisER enabled us to identify small exons with a higher percent spliced-in index (PSI) in neural, particularly neural-regulated microexons, when comparing 14 neural to 16 non-neural tissues in humans. Our work introduces an improved quantification method for long-read RNA-seq and especially facilitates studies using ONT long-reads to elucidate the regulation of genes involving small exons.

New insights into transcriptome variation during cattle adipocyte adipogenesis by direct RNA sequencing

We performed direct RNA sequencing (DRS) together with PCR-amplified cDNA long and short read sequencing for cattle adipocyte at different stages. We proved that the DRS was with advantages to avoid artificial transcripts and questionable exitrons. Totally, we obtained 68,124 transcripts with information of alternative splicing, poly (A) length and mRNA modification. The number of transcripts for adipogenesis was expanded by alternative splicing, which lead regulation mechanisms far more complex than ever known. We detected 891 differentially expressed genes (DEGs). However, 62.78% transcripts of DEGs were not significantly differentially expressed, and 248 transcripts showed opposite changing directions with their genes. The poly (A) tail became globally shorter in differentiated adipocyte than in primary adipocyte, and had a weak negative correlation with gene/transcript expression. Moreover, the study of different mRNA modifications implied their potential roles in gene expression and alternative splicing. Overall, our study promoted better understanding of adipogenesis mechanisms in cattle adipocytes.

The transcriptional landscape of endogenous retroelements delineates esophageal adenocarcinoma subtypes

Most cancer types exhibit aberrant transcriptional activity, including derepression of retrotransposable elements (RTEs). However, the degree, specificity and potential consequences of RTE transcriptional activation may differ substantially among cancer types and subtypes. Representing one extreme of the spectrum, we characterize the transcriptional activity of RTEs in cohorts of esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) from the OCCAMS (Oesophageal Cancer Clinical and Molecular Stratification) consortium, and from TCGA (The Cancer Genome Atlas). We found exceptionally high RTE inclusion in the EAC transcriptome, driven primarily by transcription of genes incorporating intronic or adjacent RTEs, rather than by autonomous RTE transcription. Nevertheless, numerous chimeric transcripts straddling RTEs and genes, and transcripts from stand-alone RTEs, particularly KLF5- and SOX9-controlled HERVH proviruses, were overexpressed specifically in EAC. Notably, incomplete mRNA splicing and EAC-characteristic intronic RTE inclusion was mirrored by relative loss of the respective fully-spliced, functional mRNA isoforms, consistent with compromised cellular fitness. Defective RNA splicing was linked with strong transcriptional activation of a HERVH provirus on Chr Xp22.32 and defined EAC subtypes with distinct molecular features and prognosis. Our study defines distinguishable RTE transcriptional profiles of EAC, reflecting distinct underlying processes and prognosis, thus providing a framework for targeted studies.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies

Aberrant skipping of coding exons in CD19 and CD22 compromises responses to immunotherapy for B-cell malignancies. Here, we show that the MS4A1 gene encoding human CD20 also produces several mRNA isoforms with distinct 5' untranslated regions (5'-UTR). Four variants (V1-4) were detectable by RNA-seq in distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant by far. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform was found to contain upstream open reading frames (uORFs) and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching Morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, while V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed CAR T cells were able to kill both V3- and V1-expressing cells, but the bispecific T cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on four post-mosunetuzumab follicular lymphoma relapses and discovered that in two of them downregulation of CD20 was accompanied by the V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies.

Key Points

In normal & malignant human B cells, CD20 mRNA is alternatively spliced into four 5'-UTR isoforms, some of which are translation-deficient.The balance between translation-deficient and -competent isoforms modulates CD20 protein levels & responses to CD20-directed immunotherapies.

Explanation of Novelty

We discovered that in normal and malignant B-cells, CD20 mRNA is alternatively spliced to generate four distinct 5'-UTRs, including the longer translation-deficient V1 variant. Cells predominantly expressing V1 were still sensitive to CD20-targeting chimeric antigen receptor T-cells. However, they were resistant to the bispecific anti-CD3/CD20 antibody mosunetuzumab, and the shift to V1 were observed in CD20-negative post-mosunetuzumab relapses of follicular lymphoma.

Artifacts and biases of the reverse transcription reaction in RNA sequencing

RNA sequencing has spurred a significant number of research areas in recent years. Most protocols rely on synthesizing a more stable complementary DNA (cDNA) copy of the RNA molecule during the reverse transcription reaction. The resulting cDNA pool is often wrongfully assumed to be quantitatively and molecularly similar to the original RNA input. Sadly, biases and artifacts confound the resulting cDNA mixture. These issues are often overlooked or ignored in the literature by those that rely on the reverse transcription process. In this review, we confront the reader with intra- and intersample biases and artifacts caused by the reverse transcription reaction during RNA sequencing experiments. To fight the reader's despair, we also provide solutions to most issues and inform on good RNA sequencing practices. We hope the reader can use this review to their advantage, thereby contributing to scientifically sound RNA studies.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

Landscape of exitrons in gastric cancer

Background

Exitron is a new type of non-canonical alternative splicing. Accumulating evidence implies exitron may have pathological function and contribute to another source of anti-tumor immunogenicity in various cancers. Its role in gastric cancer remains poorly understood. Large-scale, multi-omics analysis could comprehensively characterize the landscape of exitrons in gastric cancer, reveal undiscovered mechanism and hopefully identify molecular biomarkers for predicting immunotherapy response.

Methods

We collected datasets from five studies for analysis. RNA sequencing was used for exitron identification. Somatic mutations were detected by whole exome sequencing. Neopeptides were confirmed by proteome mass spectrometry.

Findings

42174 gastric cancer-specific exitrons (GCSEs) were identified in 632 patients. GCSEs were clinically relevant to gender, age, Lauren type, tumor stage and prognosis. Tissue specificity test and pathogenic exitron prediction revealed their unique functional impact. GCSEs were mutually exclusive with mutations and demonstrated both unique and complementary function against TP53 mutation in gastric cancer. We further established splicing regulatory network to reveal upstream regulation of exitron splicing. We also evaluated the immunogenicity and diagnostic potential of GCSEs. Evidence of GCSEs-derived neopeptide expression was validated by whole proteome mass spectrometry. PD-1 and Siglecs were significantly increased in high neoantigen load patients. But exitron-related biomarkers failed to predict immunotherapy response, possibly due to small sample size and insufficient sequencing depth.

Interpretation

The present study provided a comprehensive multidimensional landscape of gastric cancer exitrons and underscores insights into underexplored mechanism in gastric cancer pathology.

Funding

The Guangdong Provincial Key Laboratory of Precision Medicine for Gastroinstestinal Cancer (2020B121201004), the Guangdong Provincial Major Talents Project (No. 2019JC05Y361) and National Natural Science Foundation of China (grant number:82172960 and 81872013).

High-throughput mutagenesis identifies mutations and RNA-binding proteins controlling CD19 splicing and CART-19 therapy resistance

Following CART-19 immunotherapy for B-cell acute lymphoblastic leukaemia (B-ALL), many patients relapse due to loss of the cognate CD19 epitope. Since epitope loss can be caused by aberrant CD19 exon 2 processing, we herein investigate the regulatory code that controls CD19 splicing. We combine high-throughput mutagenesis with mathematical modelling to quantitatively disentangle the effects of all mutations in the region comprising CD19 exons 1-3. Thereupon, we identify ~200 single point mutations that alter CD19 splicing and thus could predispose B-ALL patients to developing CART-19 resistance. Furthermore, we report almost 100 previously unknown splice isoforms that emerge from cryptic splice sites and likely encode non-functional CD19 proteins. We further identify cis-regulatory elements and trans-acting RNA-binding proteins that control CD19 splicing (e.g., PTBP1 and SF3B4) and validate that loss of these factors leads to pervasive CD19 mis-splicing. Our dataset represents a comprehensive resource for identifying predictive biomarkers for CART-19 therapy.

Multiple alternative splicing events in CD19 mRNA have been associated with resistance/relapse to CD19 CAR-T therapy in patients with B cell malignancies. Here, by combining patient data and a high-throughput mutagenesis screen, the authors identify single point mutations and RNA-binding proteins that can control CD19 splicing and be associated with CD19 CAR-T therapy resistance.

ScisorWiz: Visualizing Differential Isoform Expression in Single-Cell Long-Read Data

Summary

RNA isoforms contribute to the diverse functionality of the proteins they encode within the cell. Visualizing how isoform expression differs across cell types and brain regions can inform our understanding of disease and gain or loss of functionality caused by alternative splicing with potential negative impacts. However, the extent to which this occurs in specific cell types and brain regions is largely unknown. This is the kind of information that ScisorWiz plots can provide in an informative and easily communicable manner. ScisorWiz affords its user the opportunity to visualize specific genes across any number of cell types, and provides various sorting options for the user to gain different ways to understand their data. ScisorWiz provides a clear picture of differential isoform expression through various clustering methods and highlights features such as alternative exons and single-nucleotide variants. Tools like ScisorWiz are key for interpreting single-cell isoform sequencing data. This tool applies to any single-cell long-read RNA sequencing data in any cell type, tissue or species.

Availability and implementation

Source code is available at http://github.com/ans4013/ScisorWiz. No new data were generated for this publication. Data used to generate figures was sourced from GEO accession token GSE158450 and available on GitHub as example data.