Abstract:Background: Solving the structure of mRNA transcripts is a major challenge for both research and molecular diagnostic purposes. Current approaches based on short-read RNA sequencing and RT-PCR techniques cannot fully explore the complexity of transcript structure. The emergence of third-generation long-read sequencing addresses this problem by solving this sequence directly. However, genes with low expression levels are difficult to study with the whole transcriptome sequencing approach. To fix this technical limitation, we propose a novel method to capture transcripts of a gene panel using a targeted enrichment approach suitable for Pacific Biosciences and Oxford Nanopore Technologies platforms. Results: We designed a set of probes to capture transcripts of a panel of genes involved in hereditary breast and ovarian cancer syndrome. We present SOSTAR (iSofOrmS annoTAtoR), a versatile pipeline to assemble, quantify and annotate isoforms from long read sequencing using a new tool specially designed for this application. The significant enrichment of transcripts by our capture protocol, together with the SOSTAR annotation, allowed the identification of 1,231 unique transcripts within the gene panel from the eight patients sequenced. The structure of these transcripts was annotated with a resolution of one base relative to a reference transcript. All major alternative splicing events of the BRCA1 and BRCA2 genes described in the literature were found. Complex splicing events such as pseudoexons were correctly annotated. SOSTAR enabled the identification of abnormal transcripts in the positive controls. In addition, a case of unexplained inheritance in a family with a history of breast and ovarian cancer was solved by identifying an SVA retrotransposon in intron 13 of the BRCA1 gene. Conclusions: We have validated a new protocol for the enrichment of transcripts of interest using probes adapted to the ONT and PacBio platforms. This protocol allows a complete description of the alternative structures of transcripts, the estimation of their expression and the identification of aberrant transcripts in a single experiment. This proof-of-concept opens new possibilities for RNA structure exploration in both research and molecular diagnostics.

Transipedia.org: k-mer-based exploration of large RNA sequencing datasets and application to cancer data

Exploring a large cancer cell line RNA-sequencing dataset with k-mers

Fleximer: Accurate Quantification of RNA-Seq via Variable-Length k-mers

Effective requesting method to detect fusion transcripts in chronic myelomonocytic leukemia RNA-seq

RNA-Skim: a Rapid Method for RNA-Seq Quantification at Transcript Level.

Accurate quantification of single-cell and single-nucleus RNA-seq transcripts using distinguishing flanking k-mers

Reference-free inferring of transcriptomic events in cancer cells on single-cell data

A Survey of K-mer Methods and Applications in Bioinformatics

A cloud-based workflow to quantify transcript-expression levels in public cancer compendia

Indexing and real-time user-friendly queries in terabyte-sized complex genomic datasets with kmindex and ORA

Molecular indexing enables quantitative targeted RNA sequencing and reveals poor efficiencies in standard library preparations

Guide to k-mer approaches for genomics across the tree of life

FusionQ: a Novel Approach for Gene Fusion Detection and Quantification from Paired-End RNA-Seq.

RNASEQR—a Streamlined and Accurate RNA-seq Sequence Analysis Program

DORQ-seq: high-throughput quantification of femtomol tRNA pools by combination of cDNA hybridization and Deep sequencing

Towards the Latent Transcriptome

Leveraging Basecaller's Move Table to Generate a Lightweight k-mer Model

Improved functions for nonlinear sequence comparison using SEEKR

Single-Cell Multimodal Analytical Approach by Integrating Raman Optical Tweezers and RNA Sequencing

Ensemble of nucleic acid absolute quantitation modules for copy number variation detection and RNA profiling

Fine mapping of RNA isoform diversity using an innovative targeted long-read RNA sequencing protocol with novel dedicated bioinformatics pipeline