Correction to: Targeted Long-Read RNA Sequencing Demonstrates Transcriptional Diversity Driven by Splice-Site Variation in
            MYBPC3

ABSTRACTBackgroundClinical sequencing has traditionally focused on genomic DNA through the use of targeted panels and exome sequencing, rather than investigating the potential transcriptomic consequences of disease-associated variants. RNA sequencing has recently been shown to be an effective additional tool for identifying disease-causing variants. We here use targeted long-read genome and transcriptome sequencing to efficiently and economically identify molecular consequences of a rare, disease-associated variant in hypertrophic cardiomyopathy (HCM).Methods and ResultsOur study, which employed both Pacific Biosciences SMRT sequencing and Oxford Nanopore Technologies MinION sequencing, as well as two RNA targeting strategies, identified alternatively-spliced isoforms that resulted from a splice-site variant containing allele in HCM. These included a predicted in-frame exon-skipping event, as well as an abundance of additional isoforms with unexpected intron-inclusion, exon-extension, and pseudo-exon events. The use of long-read RNA sequencing allowed us to not only investigate full length alternatively-spliced transcripts but also to phase them back to the variant-containing allele.ConclusionsWe suggest that targeted, long-read RNA sequencing in conjunction with genome sequencing may provide additional molecular evidence of disease for rare or de novo variants in cardiovascular disease, as well as providing new information about the consequence of these variants on downstream RNA and protein expression.

Download Full-text

TagSeqTools: a flexible and comprehensive analysis pipeline for NAD tagSeq data

10.1101/2020.03.09.982934 ◽

2020 ◽

Cited By ~ 1

Author(s):

Huan Zhong ◽

Zongwei Cai ◽

Zhu Yang ◽

Yiji Xia

Keyword(s):

Rna Sequencing ◽

Comprehensive Analysis ◽

Enzymatic Reactions ◽

Computational Tool ◽

Sequencing Data ◽

Analysis Pipeline ◽

Oxford Nanopore ◽

Long Read ◽

Identification And Characterization

AbstractNAD tagSeq has recently been developed for the identification and characterization of NAD+-capped RNAs (NAD-RNAs). This method adopts a strategy of chemo-enzymatic reactions to label the NAD-RNAs with a synthetic RNA tag before subjecting to the Oxford Nanopore direct RNA sequencing. A computational tool designed for analyzing the sequencing data of tagged RNA will facilitate the broader application of this method. Hence, we introduce TagSeqTools as a flexible, general pipeline for the identification and quantification of tagged RNAs (i.e., NAD+-capped RNAs) using long-read transcriptome sequencing data generated by NAD tagSeq method. TagSeqTools comprises two major modules, TagSeek for differentiating tagged and untagged reads, and TagSeqQuant for the quantitative and further characterization analysis of genes and isoforms. Besides, the pipeline also integrates some advanced functions to identify antisense or splicing, and supports the data reformation for visualization. Therefore, TagSeqTools provides a convenient and comprehensive workflow for researchers to analyze the data produced by the NAD tagSeq method or other tagging-based experiments using Oxford nanopore direct RNA sequencing. The pipeline is available at https://github.com/dorothyzh/TagSeqTools, under Apache License 2.0.

Download Full-text

Abstract 2724: A combination of short-read and long-read RNA sequencing reveals NOVA1’s role in telomere biology

10.1158/1538-7445.sabcs18-2724 ◽

2019 ◽

Author(s):

Andrew T. Ludlow ◽

Mohammed E. Sayed ◽

Aaron L. Slusher ◽

Mark Ribick ◽

Anisha Pancholi ◽

...

Keyword(s):

Rna Sequencing ◽

Short Read ◽

Long Read ◽

Telomere Biology

Download Full-text

Swan: a library for the analysis and visualization of long-read transcriptomes

Bioinformatics ◽

10.1093/bioinformatics/btaa836 ◽

2020 ◽

Author(s):

Fairlie Reese ◽

Ali Mortazavi

Keyword(s):

Rna Sequencing ◽

Cell Lines ◽

Intron Retention ◽

Exon Skipping ◽

Differentially Expressed ◽

Transcript Isoforms ◽

Sequencing Technologies ◽

Oxford Nanopore ◽

Long Read

Abstract Motivation Long-read RNA-sequencing technologies such as PacBio and Oxford Nanopore have discovered an explosion of new transcript isoforms that are difficult to visually analyze using currently available tools. We introduce the Swan Python library, which is designed to analyze and visualize transcript models. Results Swan finds 4909 differentially expressed transcripts between cell lines HepG2 and HFFc6, including 279 that are differentially expressed even though the parent gene is not. Additionally, Swan discovers 285 reproducible exon skipping and 47 intron retention events not recorded in the GENCODE v29 annotation. Availability and implementation The Swan library for Python 3 is available on PyPi at https://pypi.org/project/swan-vis/ and on GitHub at https://github.com/mortazavilab/swan_vis.

Download Full-text

Whole Genome Sequencing and RNA Sequencing of 27 Patients with Persistent Polyclonal B-Cell Lymphocytosis Reveals a High Mutation Frequency/Overexpression of Lymphoma Associated Genes: Really a Benign Disorder?

Blood ◽

10.1182/blood-2018-99-111480 ◽

2018 ◽

Vol 132 (Supplement 1) ◽

pp. 924-924

Author(s):

Anna Stengel ◽

Alexander Höllein ◽

Wolfgang Kern ◽

Manja Meggendorfer ◽

Claudia Haferlach ◽

...

Keyword(s):

Gene Expression ◽

B Cell ◽

Whole Genome Sequencing ◽

Rna Sequencing ◽

Genome Sequencing ◽

Splice Site ◽

Mutation Frequency ◽

Whole Genome ◽

Employment Equity ◽

Equity Ownership

Abstract Background: Persistent polyclonal B-cell lymphocytosis (PPBL) is a rare disorder, occurs almost exclusively in smoking women and is characterized by a chronic polyclonal lymphocytosis with circulating binucleated lymphocytes, clonal cytogenetic abnormalities involving chromosome 3, and chromosomal instability. Outcome of PPBL patients is mostly benign, but subsequent malignancies (non-Hodgkin´s lymphomas and solid tumors) were described. Potential molecular factors leading to their development are yet unclear. Aims: Detailed molecular genetic characterization of PPBL by whole genome sequencing (WGS) and RNA sequencing (RNAseq) in comparison to the well-characterized lymphoid malignancy CLL. Patient cohorts and methods: The total cohort comprised 27 PPBL (3 male, 24 female) and 250 CLL cases (163 male, 87 female). WGS was performed for all patients: 150bp paired-end reads where generated on Illumina HiseqX and NovaSeq 6000 machines (Illumina, San Diego, CA). A mixture genomic DNA from multiple anonymous donors was used as normal controls. To remove potential germline variants, each variant was queried against the gnomAD database, variants with global population frequencies >1% where excluded. Final analysis was performed only on protein-altering and splice-site variants. For further analysis, a virtual panel of 355 lymphoid genes was selected. All reported p-values are two-sided and were considered significant at p<0.05. For gene expression analysis, estimated gene counts were normalized applying Trimmed mean of M-values (TMM) normalization method and the resulting log2 counts per million (CPMs) were used as a proxy of gene expression in each sample. Genes were kept if they were expressed (> 5 CPM) in at least 66% of the samples. Genes with FDR (false discovery rate) < 0.05 and an absolute logFC > 1.5 were considered differentially expressed (DE). Results: Median age was 46 years for PPBL patients (range: 23-67 years) and 67 years for CLL patients (range: 39-94 years). Mean number of mutations per patient was 18 for PPBL and 20 for CLL. For both entities, the majority of mutations were missense mutations (88% in PPBL vs. 81% in CLL), followed by splice-site mutations (7% vs. 10%), other mutation types were only rarely detected. In PPBL, 42 genes were found to be mutated at a frequency of >15%, including ATM (22%), CREBBP (19%), NCOR2 (19%), AHNAK2 (15%), JAK3 (15%), NOTCH2 (15%) and TRAF1 (15%), all of which have been associated with a variety of cancers. Moreover, ATM, NOTCH2 and TRAF1 mutations were described before to be associated with lymphomas. In PPBL patients, mutations in TRAF1 and ATM as well as mutations in TRAF1 and NOTCH2 were found to be mutually exclusive. For CLL patients, 29 genes showed a mutation frequency of >15%, comprising ATM (26%), KMT2D (23%), NOTCH1 (23%), LRP1B (19%), TP53 (16%) and CREBBP (15%). Comparison of the mutation frequencies between the two entities revealed several genes with significant differences: whereas mutations in CKAP5 (11% vs. 2%, p=0.022), DNMT3A (11% vs. 3%, p=0.033), MAP2 (19% vs. 4%, p=0.009), ROBO1 (15% vs. 4%, p=0.046) and TRAF1 (15% vs. 2%, p=0.006) were found to be more frequent in PPBL cases compared to CLL cases, KMT2D (4% vs. 23%, p=0.014), TDRD6 (0% vs. 14%, p=0.032) and TP53 (4% vs. 16%, p=0.048) mutations were more abundantly detected in CLL patients. Moreover, NOTCH1 was mutated more frequently in CLL cases (7% vs. 23%, p=0.082), whereas mutated NOTCH2 (known to be frequently mutated in splenic marginal zone lymphoma), was more abundant in PPBL patients (15% vs. 6%, p=0.116), although both correlations were not statistically significant. Gene expression analyses by RNAseq revealed 337 genes to be differentially expressed between the entities. 207 genes were upregulated in PPBL, including PTPRK, CXCR1, BCL11B, CEPBA, CCR4 and MYC, whereas 130 genes were found to be upregulated in CLL cases, comprising ID3, BCL2, FGF2 and FLT1. Conclusions: 1) WGS analysis identifies high frequencies of cancer/lymphoma-associated gene mutations in PPBL, including mutated ATM, NOTCH2 and TRAF1. 2) Five genes showed a higher mutation frequency compared to CLL including TRAF1,DNMT3A, CKAP5 and MAP2. 3) Lymphoma associated genes (BCL11B and MYC) were overexpressed in PPBL vs CLL. 4) Taken together our results question PPBL as a benign entity and identify molecular markers that might contribute to development of subsequent malignancies. Disclosures Stengel: MLL Munich Leukemia Laboratory: Employment. Höllein:MLL Munich Leukemia Laboratory: Employment. Kern:MLL Munich Leukemia Laboratory: Employment, Equity Ownership. Meggendorfer:MLL Munich Leukemia Laboratory: Employment. Haferlach:MLL Munich Leukemia Laboratory: Employment, Equity Ownership. Walter:MLL Munich Leukemia Laboratory: Employment. Hutter:MLL Munich Leukemia Laboratory: Employment. Haferlach:MLL Munich Leukemia Laboratory: Employment, Equity Ownership.

Download Full-text

Accurate spliced alignment of long RNA sequencing reads

10.1101/2020.09.02.279208 ◽

2020 ◽

Author(s):

Kristoffer Sahlin ◽

Veli Mäkinen

Keyword(s):

Rna Sequencing ◽

State Of The Art ◽

Synthetic Data ◽

Biological Data ◽

Alignment Method ◽

Spliced Alignment ◽

Sequencing Technique ◽

Long Read ◽

New Challenges ◽

Novel Isoforms

AbstractLong-read RNA sequencing techniques are quickly establishing themselves as the primary sequencing technique to study the transcriptome landscape. Many such analyses are dependent upon splice alignment of reads to the genome. However, the error rate and sequencing length of long-read technologies create new challenges for accurately aligning these reads. We present an alignment method uLTRA that, on simulated and synthetic data, shows higher accuracy over state-of-the-art with substantially higher accuracy for small exons. We show several examples on biological data where uLTRA aligns to known and novel isoforms with exon structures that are not detected with other aligners. uLTRA is available at https://github.com/ksahlin/ultra.

Download Full-text

RNA sequencing data: hitchhiker's guide to expression analysis

10.7287/peerj.preprints.27283 ◽

2018 ◽

Author(s):

Koen Van Den Berge ◽

Katharina Hembach ◽

Charlotte Soneson ◽

Simone Tiberi ◽

Lieven Clement ◽

...

Keyword(s):

Gene Expression ◽

Rna Sequencing ◽

Large Scale ◽

Science Studies ◽

Rna Seq ◽

Sequencing Data ◽

Data Types ◽

The Past ◽

Long Read ◽

Statistical Approaches

Gene expression is the fundamental level at which the result of various genetic and regulatory programs are observable. The measurement of transcriptome-wide gene expression has convincingly switched from microarrays to sequencing in a matter of years. RNA sequencing (RNA-seq) provides a quantitative and open system for profiling transcriptional outcomes on a large scale and therefore facilitates a large diversity of applications, including basic science studies, but also agricultural or clinical situations. In the past 10 years or so, much has been learned about the characteristics of the RNA-seq datasets as well as the performance of the myriad of methods developed. In this review, we give an overall view of the developments in RNA-seq data analysis, including experimental design, with an explicit focus on quantification of gene expression and statistical approaches for differential expression. We also highlight emerging data types, such as single-cell RNA-seq and gene expression profiling using long-read technologies.

Download Full-text

Nanopore-Based Direct RNA-Sequencing Reveals a High-Resolution Transcriptional Landscape of Porcine Reproductive and Respiratory Syndrome Virus

Viruses ◽

10.3390/v13122531 ◽

2021 ◽

Vol 13 (12) ◽

pp. 2531

Author(s):

Riteng Zhang ◽

Peixin Wang ◽

Xin Ma ◽

Yifan Wu ◽

Chen Luo ◽

...

Keyword(s):

High Resolution ◽

Rna Sequencing ◽

Respiratory Syndrome Virus ◽

Regulatory Sequences ◽

Early Passage ◽

Rna Molecules ◽

Transcription Process ◽

Long Read ◽

Potential Link ◽

Transcriptional Landscape

The TRS-mediated discontinuous transcription process is a hallmark of Arteriviruses. Precise assessment of the intricate subgenomic RNA (sg mRNA) populations is required to understand the kinetics of viral transcription. It is difficult to reconstruct and comprehensively quantify splicing events using short-read sequencing, making the identification of transcription-regulatory sequences (TRS) particularly problematic. Here, we applied long-read direct RNA sequencing to characterize the recombined RNA molecules produced in porcine alveolar macrophages during early passage infection of porcine reproductive and respiratory syndrome virus (PRRSV). Based on sequencing two PRRSV isolates, namely XM-2020 and GD, we revealed a high-resolution and diverse transcriptional landscape in PRRSV. The data revealed intriguing differences in subgenomic recombination types between the two PRRSVs while also demonstrating TRS-independent heterogeneous subpopulation not previously observed in Arteriviruses. We find that TRS usage is a regulated process and share the common preferred TRS in both strains. This study also identified a substantial number of TRS-mediated transcript variants, including alternative-sg mRNAs encoding the same annotated ORF, as well as putative sg mRNAs encoded nested internal ORFs, implying that the genetic information encoded in PRRSV may be more intensively expressed. Epigenetic modifications have emerged as an essential regulatory layer in gene expression. Here, we gained a deeper understanding of m5C modification in poly(A) RNA, elucidating a potential link between methylation and transcriptional regulation. Collectively, our findings provided meaningful insights for redefining the transcriptome complexity of PRRSV. This will assist in filling the research gaps and developing strategies for better control of the PRRS.

Download Full-text

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

Genome Biology ◽

10.1186/s13059-021-02411-1 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Laura Schulz ◽

Manuel Torres-Diz ◽

Mariela Cortés-López ◽

Katharina E. Hayer ◽

Mukta Asnani ◽

...

Keyword(s):

Rna Sequencing ◽

Reverse Transcription ◽

Lymphoblastic Leukemia ◽

Mrna Splicing ◽

Exon 2 ◽

Long Read ◽

Reporter Assays ◽

Cryptic Intron

AbstractResistance to CD19-directed immunotherapies in lymphoblastic leukemia has been attributed, among other factors, to several aberrant CD19 pre-mRNA splicing events, including recently reported excision of a cryptic intron embedded within CD19 exon 2. While “exitrons” are known to exist in hundreds of human transcripts, we discovered, using reporter assays and direct long-read RNA sequencing (dRNA-seq), that the CD19 exitron is an artifact of reverse transcription. Extending our analysis to publicly available datasets, we identified dozens of questionable exitrons, dubbed “falsitrons,” that appear only in cDNA-seq, but never in dRNA-seq. Our results highlight the importance of dRNA-seq for transcript isoform validation.

Download Full-text