From bulk, single-cell to spatial RNA sequencing

Abstract

RNA sequencing (RNAseq) can reveal gene fusions, splicing variants, mutations/indels in addition to differential gene expression, thus providing a more complete genetic picture than DNA sequencing. This most widely used technology in genomics tool box has evolved from classic bulk RNA sequencing (RNAseq), popular single cell RNA sequencing (scRNAseq) to newly emerged spatial RNA sequencing (spRNAseq). Bulk RNAseq studies average global gene expression, scRNAseq investigates single cell RNA biology up to 20,000 individual cells simultaneously, while spRNAseq has ability to dissect RNA activities spatially, representing next generation of RNA sequencing. This article highlights these technologies, characteristic features and suitable applications in precision oncology.

References

[1] Emrich SJ, Barbazuk WB, Li L, Schnable PS. Gene discovery and annotation using LCM-454 transcriptome sequencing. Genome Res., 2007, 17: 69-73.

[2] Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet, 2009, 10: 57-63.

[3] Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat. Rev. Genet, 2011, 12: 87-98.

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

[5] Han LO, Li XY, Cao MM, Cao Y, Zhou LH. Development and validation of an individualized diagnostic signature in thyroid cancer. Cancer Med, 2018, 7: 1135-1140.

[6] Zhou JG, . Development and validation of an RNA-seq-based prognostic signature in neuroblastoma. Front Oncol., 2019, 9: 1361.

[7] Shukla, S. et al. Development of a RNA-seq based prognostic signature in lung adenocarcinoma. J. Natl Cancer Inst. 109, https://doi.org/10.1093/jnci/djw200 (2017).

[8] Chen PL, . Analysis of immune signatures in longitudinal tumor samples yields insight into biomarkers of response and mechanisms of resistance to immune checkpoint blockade. Cancer Disco., 2016, 6: 827-837.

[9] Van Allen EM, . Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science, 2015, 350: 207-211.

[10] Gyanchandani R, . Intratumor heterogeneity affects gene expression profile test prognostic risk stratification in early breast cancer. Clin. Cancer Res, 2016, 22: 5362-5369.

[11] Biswas D, . A clonal expression biomarker associates with lung cancer mortality. Nat. Med, 2019, 25: 1540-1548.

[12] Vardiman JW, . The 2008 revision of the World Health Organization (WHO) classification of myeloid neoplasms and acute leukemia: rationale and important changes. Blood, 2009, 114: 937-951.

[13] Font-Tello A, . Association of ERG and TMPRSS2-ERG with grade, stage, and prognosis of prostate cancer is dependent on their expression levels. Prostate, 2015, 75: 1216-1226.

[14] Stransky N, Cerami E, Schalm S, Kim JL, Lengauer C. The landscape of kinase fusions in cancer. Nat. Commun., 2014, 5

[15] Roberts KG, . Genetic alterations activating kinase and cytokine receptor signaling in high-risk acute lymphoblastic leukemia. Cancer Cell, 2012, 22: 153-166.

[16] Roberts KG, . Targetable kinase-activating lesions in Ph-like acute lymphoblastic leukemia. N. Engl. J. Med, 2014, 371: 1005-1015.

[17] Shaw AT, Hsu PP, Awad MM, Engelman JA. Tyrosine kinase gene rearrangements in epithelial malignancies. Nat. Rev. Cancer, 2013, 13: 772-787.

[18] Dehghannasiri R, . Improved detection of gene fusions by applying statistical methods reveals oncogenic RNA cancer drivers. Proc. Natl Acad. Sci. USA, 2019, 116: 15524-15533.

[19] Heyer EE, . Diagnosis of fusion genes using targeted RNA sequencing. Nat. Commun., 2019, 10

[20] Gambino G, Tancredi M, Falaschi E, Aretini P, Caligo MA. Characterization of three alternative transcripts of the BRCA1 gene in patients with breast cancer and a family history of breast and/or ovarian cancer who tested negative for pathogenic mutations. Int J. Mol. Med, 2015, 35: 950-956.

[21] Cabanski CR, . cDNA hybrid capture improves transcriptome analysis on low-input and archived samples. J. Mol. Diagn., 2014, 16: 440-451.

[22] Valle L, . Germline allele-specific expression of TGFBR1 confers an increased risk of colorectal cancer. Science, 2008, 321: 1361-1365.

[23] Tang RX, . Identification of a RNA-Seq based prognostic signature with five lncRNAs for lung squamous cell carcinoma. Oncotarget, 2017, 8: 50761-50773.

[24] Siriwardhana C, Khadka VS, Chen JJ, Deng Y. Development of a miRNA-seq based prognostic signature in lung adenocarcinoma. BMC Cancer, 2019, 19

[25] Haque A, Engel J, Teichmann SA, Lonnberg T. A practical guide to single-cell RNA-sequencing for biomedical research and clinical applications. Genome Med, 2017, 9

[26] Choi, J. R., Yong, K. W., Choi, J. Y. & Cowie, A. C. Single-cell RNA sequencing and its combination with protein and DNA analyses. Cells 9, https://doi.org/10.3390/cells9051130 (2020).

[27] Arozarena I, Wellbrock C. Phenotype plasticity as enabler of melanoma progression and therapy resistance. Nat. Rev. Cancer, 2019, 19: 377-391.

[28] McGranahan N, Swanton C. Clonal heterogeneity and tumor evolution: past, present, and the future. Cell, 2017, 168: 613-628.

[29] Patel AP, . Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science, 2014, 344: 1396-1401.

[30] Roerink SF, . Intra-tumour diversification in colorectal cancer at the single-cell level. Nature, 2018, 556: 457-462.

[31] Puram SV, . Single-cell transcriptomic analysis of primary and metastatic tumor ecosystems in head and neck cancer. Cell, 2017, 171: 1611-1624 e1624.

[32] Kim KT, . Application of single-cell RNA sequencing in optimizing a combinatorial therapeutic strategy in metastatic renal cell carcinoma. Genome Biol., 2016, 17

[33] Tirosh I, . Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science, 2016, 352: 189-196.

[34] Lee MC, . Single-cell analyses of transcriptional heterogeneity during drug tolerance transition in cancer cells by RNA sequencing. Proc. Natl Acad. Sci. USA, 2014, 111: E4726-E4735.

[35] Guo X, . Publisher Correction: Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing. Nat. Med, 2018, 24: 1628.

[36] Krishna S, . Stem-like CD8 T cells mediate response of adoptive cell immunotherapy against human cancer. Science, 2020, 370: 1328-1334.

[37] Zhang Q, . Landscape and dynamics of single immune cells in hepatocellular carcinoma. Cell, 2019, 179: 829-845 e820.

[38] Wang C, . CD276 expression CD276 expression enables squamous cell carcinoma stem cells to evade immune surveillance. Cell Stem Cell, 2021, 28: 1597-1613.e7.

[39] Cheng YH, . Hydro-Seq enables contamination-free high-throughput single-cell RNA-sequencing for circulating tumor cells. Nat. Commun., 2019, 10

[40] Method of the Year 2019: Single-cell multimodal omics. Nat Methods 17, 1, https://doi.org/10.1038/s41592-019-0703-5 (2020).

[41] Method of the year 2013. Nat Methods 11, 1, https://doi.org/10.1038/nmeth.2801 (2014).

[42] Stuart T, Satija R. Integrative single-cell analysis. Nat. Rev. Genet, 2019, 20: 257-272.

[43] Zhu C, Preissl S, Ren B. Single-cell multimodal omics: the power of many. Nat. Methods, 2020, 17: 11-14.

[44] Stahl PL, . Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science, 2016, 353: 78-82.

[45] Salmen F, . Barcoded solid-phase RNA capture for Spatial Transcriptomics profiling in mammalian tissue sections. Nat. Protoc., 2018, 13: 2501-2534.

[46] Kabra M, Robie AA, Rivera-Alba M, Branson S, Branson K. JAABA: interactive machine learning for automatic annotation of animal behavior. Nat. Methods, 2013, 10: 64-67.

[47] Ihle CL, . Distinct tumor microenvironments of lytic and blastic bone metastases in prostate cancer patients. J. Immunother. Cancer, 2019, 7: 293.

[48] Amaria RN, . Neoadjuvant immune checkpoint blockade in high-risk resectable melanoma. Nat. Med, 2018, 24: 1649-1654.

[49] Blank CU, . Neoadjuvant versus adjuvant ipilimumab plus nivolumab in macroscopic stage III melanoma. Nat. Med, 2018, 24: 1655-1661.

[50] Thrane K, Eriksson H, Maaskola J, Hansson J, Lundeberg J. Spatially resolved transcriptomics enables dissection of genetic heterogeneity in stage III cutaneous malignant melanoma. Cancer Res, 2018, 78: 5970-5979.

[51] Berglund E, . Spatial maps of prostate cancer transcriptomes reveal an unexplored landscape of heterogeneity. Nat. Commun., 2018, 9

[52] Moncada R, . Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat. Biotechnol., 2020, 38: 333-342.

[53] Ji AL, . Multimodal analysis of composition and spatial architecture in human squamous cell carcinoma. Cell, 2020, 182: 1661-1662.

[54] Grauel AL, . TGFbeta-blockade uncovers stromal plasticity in tumors by revealing the existence of a subset of interferon-licensed fibroblasts. Nat. Commun., 2020, 11

[55] Deng, Y. et al. Spatial epigenome sequencing at tissue scale and cellular level. Preprint at https://doi.org/10.1101/2021.03.11.434985 (2021).

[56] Su J-H, Zheng P, Kinrot SS, Bintu B, Zhuang X. Genome-scale imaging of the 3D organization and transcriptional activity of chromatin. Cell, 2020, 182: 1641-1659.e26.

Funding

U.S. Department of Health & Human Services | NIH | National Institute of Dental and Craniofacial Research (NIDCR)(R01DE029173)