Systematic characterization of full-length RNA isoforms in human colorectal cancer at single-cell resolution

Abstract

Dysregulated RNA splicing is a well-recognized characteristic of colorectal cancer (CRC); however, its intricacies remain obscure, partly due to challenges in profiling full-length transcript variants at the single-cell level. Here, we employ high-depth long-read scRNA-seq to define the full-length transcriptome of colorectal epithelial cells in 12 CRC patients, revealing extensive isoform diversities and splicing alterations. Cancer cells exhibited increased transcript complexity, with widespread 3'-UTR shortening and reduced intron retention. Distinct splicing regulation patterns were observed between intrinsic-consensus molecular subtypes (iCMS), with iCMS3 displaying even higher splicing factor activities and more pronounced 3'-UTR shortening. Furthermore, we revealed substantial shifts in isoform usage that result in alterations of protein sequences from the same gene with distinct carcinogenic effects during tumorigenesis of CRC. Allele-specific expression analysis revealed dominant mutant allele expression in key oncogenes and tumor suppressors. Moreover, mutated PPIG was linked to widespread splicing dysregulation, and functional validation experiments confirmed its critical role in modulating RNA splicing and tumor-associated processes. Our findings highlight the transcriptomic plasticity in CRC and suggest novel candidate targets for splicing-based therapeutic strategies.

Keywords

colorectal cancer / long-read scRNA-seq / full-length RNA isoforms / RNA splicing / differential transcript usage

References

[1] Agarwal V , Lopez-Darwin S , Kelley DR et al. The landscape of alternative polyadenylation in single cells of the developing mouse embryo. Nat Commun 2021; 12: 5101.

[2] Alors-Perez E , Blázquez-Encinas R , Alcalá S et al. Dysregulated splicing factor SF3B1 unveils a dual therapeutic vulnerability to target pancreatic cancer cells and cancer stem cells with an anti-splicing drug. J Exp Clin Cancer Res 2021; 40: 382.

[3] Amarasinghe SL , Su S , Dong X et al. Opportunities and challenges in long-read sequencing data analysis. Genome Biol 2020; 21: 30.

[4] Anczukow O , Akerman M , Clery A et al. SRSF1-Regulated alternative splicing in breast cancer. Mol Cell 2015; 60: 105- 117.

[5] Bessiere D , Lacroix C , Campagne S et al. Structure-function analysis of the THAP zinc finger of THAP1, a large C2CH DNA-binding module linked to Rb/E2F pathways. J Biol Chem 2008; 283: 4352- 4363.

[6] Blencowe BJ . Alternative splicing: new insights from global analyses. Cell 2006; 126: 37- 47.

[7] Boise LH , Gonzalez-Garcia M , Postema CE et al. bcl-x, a bcl-2-related gene that functions as a dominant regulator of apoptotic cell death. Cell 1993; 74: 597- 608.

[8] Bonnal SC , Lopez-Oreja I , Valcarcel J . Roles and mechanisms of alternative splicing in cancer - implications for care. Nat Rev Clin Oncol 2020; 17: 457- 474.

[9] Bonneau D , Longy M . Mutations of the human PTEN gene. Hum Mutat 2000; 16: 109- 122.

[10] Bradley RK , Anczukow O . RNA splicing dysregulation and the hallmarks of cancer. Nat Rev Cancer 2023; 23: 135- 155.

[11] Byrne A , Beaudin AE , Olsen HE et al. Nanopore long-read RNAseq reveals widespread transcriptional variation among the surface receptors of individual B cells. Nat Commun 2017; 8: 16027.

[12] Byrne A , Cole C , Volden R et al. Realizing the potential of fulllength transcriptome sequencing. Philos Trans R Soc Lond B Biol Sci 2019; 374: 20190097.

[13] Cancer Genome Atlas Network . Comprehensive molecular characterization of human colon and rectal cancer. Nature 2012; 487: 330- 337.

[14] Castilla C , Congregado B , Chinchon D et al. Bcl-xL is overexpressed in hormone-resistant prostate cancer and promotes survival of LNCaP cells via interaction with proapoptotic Bak. Endocrinology 2006; 147: 4960- 4967.

[15] Chang MT , Bhattarai TS , Schram AM et al. Accelerating discovery of functional mutant alleles in cancer. Cancer Discov 2018; 8: 174- 183.

[16] Chen C , Zhao S , Karnad A et al. The biology and role of CD44 in cancer progression: therapeutic implications. J Hematol Oncol 2018; 11: 64.

[17] Chen B , Scurrah CR , McKinley ET et al. Differential premalignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps. Cell 2021a; 184: 6262- 6280.e26.

[18] Chen W-C , To MD , Westcott PMK et al. Targeting KRAS4A splicing through the RBM39/DCAF15 pathway inhibits cancer stem cells. Nat Commun 2021b; 12: 4288.

[19] Dai W , Zhou F , Tang D et al. Single-cell transcriptional profiling reveals the heterogenicity in colorectal cancer. Medicine (Baltimore) 2019; 98: e16916.

[20] De Coster W , D’Hert S , Schultz DT et al. NanoPack: visualizing and processing long-read sequencing data. Bioinformatics 2018; 34: 2666- 2669.

[21] Ding D , Huang H , Jiang W et al. Reticulocalbin-2 enhances hepatocellular carcinoma proliferation via modulating the EGFR-ERK pathway. Oncogene 2017; 36: 6691- 6700.

[22] Dong Y , Xue L , Zhang Y et al. Identification of RNA-splicing factor Lsm12 as a novel tumor-associated gene and a potent biomarker in Oral Squamous Cell Carcinoma(OSCC). J Exp Clin Cancer Res 2022; 41: 150.

[23] Espana L , Fernandez Y , Rubio N et al. Overexpression of Bcl-xL in human breast cancer cells enhances organ-selective lymph node metastasis. Breast Cancer Res Treat 2004; 87: 33- 44.

[24] Fan X , Tang D , Liao Y et al. Single-cell RNA-seq analysis of mouse preimplantation embryos by third-generation sequencing. PLoS Biol 2020; 18: e3001017.

[25] Fu X-D . Both sides of the same coin: Rac1 splicing regulation by EGF signaling. Cell Res 2017; 27: 455- 456.

[26] Gao S , Yan L , Wang R et al. Publisher Correction: Tracing the temporal-spatial transcriptome landscapes of the human fetal digestive tract using single-cell RNA-sequencing. Nat Cell Biol 2018; 20: 1227.

[27] Geng S , Peng W , Wang X et al. ARIH2 regulates the proliferation, DNA damage and chemosensitivity of gastric cancer cells by reducing the stability of p21 via ubiquitination. Cell Death Dis 2022; 13: 564.

[28] Geuens T , Bouhy D , Timmerman V . The hnRNP family:insights into their role in health and disease. Hum Genet 2016; 135: 851- 867.

[29] Gruber AJ , Zavolan M . Alternative cleavage and polyadenylation in health and disease. Nat Rev Genet 2019; 20: 599- 614.

[30] Gupta I , Collier PG , Haase B et al. Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells. Nat Biotechnol 2018; 36: 1197- 1202.

[31] Hann SR , Dixit M , Sears RC et al. The alternatively initiated c-Myc proteins differentially regulate transcription through a noncanonical DNA-binding site. Genes Dev 1994; 8: 2441- 2452.

[32] Hao Y , Hao S , Andersen-Nissen E et al. Integrated analysis of multimodal single-cell data. Cell 2021; 184: 3573- 3587.e29.

[33] Huang KK , Huang J , Wu JKL et al. Long-read transcriptome sequencing reveals abundant promoter diversity in distinct molecular subtypes of gastric cancer. Genome Biol 2021; 22: 44.

[34] Ikegaki N , Katsumata M , Minna J et al. Expression of bcl-2 in small cell lung carcinoma cells. Cancer Res 1994; 54: 6- 8.

[35] Jia R , Ajiro M , Yu L et al. Oncogenic splicing factor SRSF3 regulates ILF3 alternative splicing to promote cancer cell proliferation and transformation. RNA 2019; 25: 630- 644.

[36] Joanito I , Wirapati P , Zhao N et al. Single-cell and bulk transcriptome sequencing identifies two epithelial tumor cell states and refines the consensus molecular classification of colorectal cancer. Nat Genet 2022; 54: 963- 975.

[37] Joglekar A , Prjibelski A , Mahfouz A et al. A spatially resolved brain region- and cell type-specific isoform atlas of the postnatal mouse brain. Nat Commun 2021; 12: 463.

[38] Kahles A , Lehmann KV , Toussaint NC et al. Comprehensive analysis of alternative splicing across tumors from 8,705 patients. Cancer Cell 2018; 34: 211- 224.e6.

[39] Kim DW , Kim IJ , Kang HC et al. Mutation spectrum of the APC gene in 83 Korean FAP families. Hum Mutat 2005; 26: 281.

[40] Kovaka S , Zimin AV , Pertea GM et al. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol 2019; 20: 278.

[41] Kowalski MH , Wessels HH , Linder J et al. CPA-Perturbseq: Multiplexed single-cell characterization of alternative polyadenylation regulators. bioRxiv 2023: 2023.02.09.527751.

[42] Lebrigand K , Magnone V , Barbry P et al. High throughput error corrected Nanopore single cell transcriptome sequencing. Nat Commun 2020; 11: 4025.

[43] Lee Y , Rio DC . Mechanisms and regulation of alternative premRNA splicing. Annu Rev Biochem 2015; 84: 291- 323.

[44] Lee SY , Kim MJ , Jin G et al. Somatic mutations in epidermal growth factor receptor signaling pathway genes in nonsmall cell lung cancers. J Thorac Oncol 2010; 5: 1734- 1740.

[45] Lee J , Nguyen PT , Shim HS et al. EWSR1, a multifunctional protein, regulates cellular function and aging via genetic and epigenetic pathways. Biochim Biophys Acta Mol Basis Dis 2019; 1865: 1938- 1945.

[46] Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 2018; 34: 3094- 3100.

[47] Li H , Courtois ET , Sengupta D et al. Reference component analysis of single-cell transcriptomes elucidates cellular heterogeneity in human colorectal tumors. Nat Genet 2017; 49: 708- 718.

[48] Liao DJ , Dickson RB . c-Myc in breast cancer. Endocr Relat Cancer 2000; 7: 143- 164.

[49] Liao Y , Liu Z , Zhang Y et al. Publisher Correction: Highthroughput and high-sensitivity full-length single-cell RNA-seq analysis on third-generation sequencing platform. Cell Discov 2023; 9: 15.

[50] Liu J , Li H , Shen S et al. Alternative splicing events implicated in carcinogenesis and prognosis of colorectal cancer. J Cancer 2018; 9: 1754- 1764.

[51] Liu Y , Liu X , Lin C et al. Noncoding RNAs regulate alternative splicing in Cancer. J Exp Clin Cancer Res 2021; 40: 11.

[52] Love AJ , Yu C , Petukhova NV et al. Cajal bodies and their role in plant stress and disease responses. RNA Biol 2017; 14: 779- 790.

[53] Mamlouk S , Childs LH , Aust D et al. DNA copy number changes define spatial patterns of heterogeneity in colorectal cancer. Nat Commun 2017; 8: 14093.

[54] Masciari S , Dewanwala A , Stoffel EM et al. Gastric cancer in individuals with Li-Fraumeni syndrome. Genet Med 2011; 13: 651- 657.

[55] Mayr C , Bartel DP . Widespread shortening of 3'UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell 2009; 138: 673- 684.

[56] McGeary SE , Lin KS , Shi CY et al. The biochemical basis of microRNA targeting efficacy. Science 2019; 366: eaav1741.

[57] Moseley JB , Bartolini F , Okada K et al. Regulated binding of adenomatous polyposis coli protein to actin. J Biol Chem 2007; 282: 12661- 12668.

[58] Murphy D , Parker J , Zhou M et al. Constitutively overexpressed 21 kDa protein in Hodgkin lymphoma and aggressive non-Hodgkin lymphomas identified as cytochrome B5b(CYB5B). Mol Cancer 2010; 9: 14.

[59] Ning J , Liu M , Shen J et al. Expression signature and prognostic value of CREC gene family in human colorectal cancer. BMC Cancer 2023; 23: 878.

[60] O’Brien J , Hayder H , Zayed Y et al. Overview of MicroRNA biogenesis, mechanisms of actions, and circulation. Front Endocrinol (Lausanne) 2018; 9: 402.

[61] Patel AP , Tirosh I , Trombetta JJ et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 2014; 344: 1396- 1401.

[62] Pelka K , Hofree M , Chen JH et al. Spatially organized multicellular immune hubs in human colorectal cancer. Cell 2021; 184: 4734- 4752.e20.

[63] Petitjean A , Mathe E , Kato S et al. Impact of mutant p53 functional properties on TP53 mutation patterns and tumor phenotype: lessons from recent developments in the IARC TP53 database. Hum Mutat 2007; 28: 622- 629.

[64] Philpott M , Watson J , Thakurta A et al. Nanopore sequencing of single-cell transcriptomes with scCOLOR-seq. Nat Biotechnol 2021; 39: 1517- 1520.

[65] Pros E , Saigi M , Alameda D et al. Genome-wide profiling of non-smoking-related lung cancer cells reveals common RB1 rearrangements associated with histopathologic transformation in EGFR-mutant tumors. Ann Oncol 2020; 31: 274- 282.

[66] Roussigne M , Kossida S , Lavigne AC et al. The THAP domain:a novel protein motif with similarity to the DNA-binding domain of P element transposase. Trends Biochem Sci 2003; 28: 66- 69.

[67] Sakamoto Y , Sereewattanawoot S , Suzuki A . A new era of long-read sequencing for cancer genomics. J Hum Genet 2020; 65: 3- 10.

[68] Seiler M , Peng S , Agrawal AA et al. Somatic mutational landscape of splicing factor genes and their functional consequences across 33 cancer types. Cell Rep 2018; 23: 282- 296.e4.

[69] Singh M , Al-Eryani G , Carswell S et al. High-throughput targeted long-read single cell sequencing reveals the clonal and transcriptional landscape of lymphocytes. Nat Commun 2019; 10: 3120.

[70] Smillie CS , Biton M , Ordovas-Montanes J et al. Intra- and inter-cellular rewiring of the human colon during ulcerative colitis. Cell 2019; 178: 714- 730.e22.

[71] Stark R , Grzelak M , Hadfield J . RNA sequencing: the teenage years. Nat Rev Genet 2019; 20: 631- 656.

[72] Subramanian A , Tamayo P , Mootha VK et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 2005; 102: 15545- 15550.

[73] Sveen A , Kilpinen S , Ruusulehto A et al. Aberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes. Oncogene 2016; 35: 2413- 2427.

[74] Talevich E , Shain AH , Botton T et al. CNVkit: Genome-Wide copy number detection and visualization from targeted DNA sequencing. PLoS Comput Biol 2016; 12: e1004873.

[75] Tardaguila M , de la Fuente L , Marti C et al. SQANTI: extensive characterization of long-read transcript sequences for quality control in full-length transcriptome identification and quantification. Genome Res 2018; 28: 396- 411.

[76] Tekath T , Dugas M . Differential transcript usage analysis of bulk and single-cell RNA-seq data with DTUrtle. Bioinformatics 2021; 37: 3781- 3787.

[77] Tian L , Jabbari JS , Thijssen R et al. Comprehensive characterization of single-cell full-length isoforms in human and mouse with long-read sequencing. Genome Biol 2021; 22: 310.

[78] Trincado JL , Entizne JC , Hysenaj G et al. SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions. Genome Biol 2018; 19: 40.

[79] Trisciuoglio D , Tupone MG , Desideri M et al. BCL-XL overexpression promotes tumor progression-associated properties. Cell Death Dis 2017; 8: 3216.

[80] Urbanski LM , Leclair N , Anczukow O . Alternative-splicing defects in cancer: Splicing regulators and their downstream targets, guiding the way to novel cancer therapeutics. Wiley Interdiscip Rev RNA 2018; 9: e1476.

[81] Veiga DFT , Nesta A , Zhao Y et al. A comprehensive long-read isoform analysis platform and sequencing resource for breast cancer. Sci Adv 2022; 8: eabg6711.

[82] Wang ET , Sandberg R , Luo S et al. Alternative isoform regulation in human tissue transcriptomes. Nature 2008; 456: 470- 476.

[83] Wang Y , Song W , Wang J et al. Single-cell transcriptome analysis reveals differential nutrient absorption functions in human intestine. J Exp Med 2020; 217: e20191130.

[84] Wright CJ , Smith CWJ , Jiggins CD . Alternative splicing as a source of phenotypic diversity. Nat Rev Genet 2022; 23: 697- 710.

[85] Wurth L , Papasaikas P , Olmeda D et al. UNR/CSDE1 drives a post-transcriptional program to promote melanoma invasion and metastasis. Cancer Cell 2016; 30: 694- 707.

[86] Xiong Y , Deng Y , Wang K et al. Profiles of alternative splicing in colorectal cancer and their clinical significance: A study based on large-scale sequencing data. EBioMedicine 2018; 36: 183- 195.

[87] Xu T , Verhagen M , Joosten R et al. Alternative splicing downstream of EMT enhances phenotypic plasticity and malignant behaviour in colon cancer. Elife 2022; 11: e82006.

[88] You Y , Prawer YDJ , De Paoli-Iseppi R et al. Identification of cell barcodes from long-read single-cell RNA-seq with BLAZE. Genome Biol 2023; 24: 66.

[89] Yu G , Wang LG , Han Y et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 2012; 16: 284- 287.

[90] Yurgelun MB , Kulke MH , Fuchs CS et al. Cancer susceptibility gene mutations in individuals with colorectal cancer. J Clin Oncol 2017; 35: 1086- 1095.

[91] Zong Z , Li H , Yi C et al. Genome-Wide profiling of prognostic alternative splicing signature in colorectal cancer. Front Oncol 2018; 8: 537.

RIGHTS & PERMISSIONS

The Author(s). Published by Oxford University Press on behalf of Higher Education Press.