Abstract:Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene ( ) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Comparative RNA-Seq analysis reveals pervasive tissue-specific alternative polyadenylation in Caenorhabditis elegans intestine and muscles

Transcriptome analysis of adult Caenorhabditis elegans cells reveals tissue-specific gene and isoform expression

Genome-wide Identification and Predictive Modeling of Tissue-Specific Alternative Polyadenylation.

Distinct Polyadenylation Landscapes of Diverse Human Tissues Revealed by A Modified Pa-Seq Strategy

A global view of the RNA-binding and regulatory protein landscape in Caenorhabditis elegans

A Differential Sequencing-Based Analysis of the C. Elegans Noncoding Transcriptome.

New Strains for Tissue-Specific RNAi Studies in Caenorhabditis elegans

Computational Identification of Tissue-Specific Splicing Regulatory Elements in Human Genes from RNA-Seq Data

A collection of toolkit strains reveals distinct localization and dynamics of membrane-associated transcripts in epithelia

by PAS-Seq Complex and dynamic landscape of RNA polyadenylation revealed Material Supplemental

Genome-Wide Tissue-Specific Gene Expression, Co-expression and Regulation of Co-expressed Genes in Adult Nematode Ascaris suum

The conserved landscape of RNA modifications and transcript diversity across mammalian evolution

Genome-wide Prediction of Cis-Acting RNA Elements Regulating Tissue-Specific Pre-Mrna Alternative Splicing

Tissue-specific silencing of integrated transgenes achieved through endogenous RNA interference in Caenorhabditis elegans

Profiling Caenorhabditis Elegans Non-Coding RNA Expression with a Combined Microarray

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data

Comparative RNA-Seq Analysis Revealed Tissue-Specific Splicing Variations during the Generation of the PDX Model

Complex and dynamic landscape of RNA polyadenylation revealed by PAS-Seq.

Interactions between RNA m6A modification, alternative splicing, and poly(A) tail revealed by MePAIso-seq2

Comparative genomic analysis of upstream miRNA regulatory motifs in Caenorhabditis

Identification of functional long non-coding RNAs in C. elegans