Importance of transcript variants in transcriptome analyses

Kevin Vo,Ryan Mohamadi,Yashica Sharma,Amelia Mohamadi,Patrick E. Fields,M. A. Karim Rumi
DOI: https://doi.org/10.1101/2024.07.11.603122
2024-07-16
Abstract:RNA sequencing (RNA-Seq) has become a widely adopted genome-wide technique for investigating gene expression patterns. However, conventional RNA-Seq analyses typically rely on gene expression (GE) values that aggregate all the transcripts produced by a gene under a single identifier, overlooking the complexity of transcript variants arising from different transcription start sites and alternative splicing events. In this study, we explored the implications of neglecting transcript variants in RNA-Seq analyses. Among the 1334 transcription factor (TF) genes expressed in mouse embryonic stem (ES) or trophoblast stem (TS) cells, 652 were reported to be differentially expressed in TS cells based on GE values (365 upregulated and 287 downregulated, ≥2-fold, FDR -value ≤0.05). Intriguingly, differential gene expression analysis revealed that of the 365 upregulated genes, 883 transcript variants were expressed, with only 174 (<20%) variants exhibiting upregulation based on transcript expression (TE) values. The remaining 709 (>80%) variants were either down-regulated or showed no significant change in expression analysis. Similarly, the 287 genes reported to be downregulated expressed 856 transcript variants, with only 153 (<20%) downregulated variants and 703 (>82%) variants that were upregulated or showed no significant changes. Additionally, the 682 TF genes that did not show significant changes between ES and TS cells (GE values < 2-fold changes and/or FDR p-values >0.05) expressed 2215 transcript variants, which included 477 (>21%) that were differentially expressed (276 upregulated and 201 downregulated, ≥2-fold, FDR p-value ≤0.05). Notably, a particular gene does not express just one protein; rather its transcript variants encode multiple proteins with distinct functional domains, including non-coding regulatory RNAs. Our findings underscore the critical necessity of considering transcript variants in RNA-Seq analyses. Doing so may enable a more precise understanding of the intricate functional and regulatory landscape of genes; ignoring the variants may result in an erroneous interpretation.
Bioinformatics
What problem does this paper attempt to address?