Delineating the dynamic evolution from preneoplasia to invasive lung adenocarcinoma by integrating single-cell RNA sequencing and spatial transcriptomics

다국어 초록 (Multilingual Abstract)

The cell ecology and spatial niche implicated in the dynamic and sequential process of lung adenocarcinoma (LUAD) from adenocarcinoma in situ (AIS) to minimally invasive adenocarcinoma (MIA) and subsequent invasive adenocarcinoma (IAC) have not yet been elucidated. Here, we performed an integrative analysis of single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) to characterize the cell atlas of the invasion trajectory of LUAD. We found that the UBE2C + cancer cell subpopulation constantly increased during the invasive process of LUAD with remarkable elevation in IAC, and its spatial distribution was in the peripheral cancer region of the IAC, representing a more malignant phenotype. Furthermore, analysis of the TME cell type subpopulation showed a constant decrease in mast cells, monocytes, and lymphatic endothelial cells, which were implicated in the whole process of invasive LUAD, accompanied by an increase in NK cells and MALT B cells from AIS to MIA and an increase in Tregs and secretory B cells from MIA to IAC. Notably, for AIS, cancer cells, NK cells, and mast cells were colocalized in the cancer region; however, for IAC, Tregs colocalized with cancer cells. Finally, communication and interaction between cancer cells and TME cell-induced constitutive activation of TGF-β signaling were involved in the invasion of IAC. Therefore, our results reveal the specific cellular information and spatial architecture of cancer cells and TME subpopulations, as well as the cellular interaction between them, which will facilitate the identification and development of precision medicine in the invasive process of LUAD from AIS to IAC.

참고문헌 (Reference)

1 Ritchie, M. E., "limma powers differential expression analyses for RNAsequencing and microarray studies" 43 : e47-, 2015

2 Stahl, P. L., "Visualization and analysis of gene expression in tissue sections by spatial transcriptomics" 353 : 78-82, 2016

3 McInnes, L., "Umap: uniform manifold approximation and projection for dimension reduction"

4 Liu, Y., "UBE2C is upregulated by estrogen and promotes epithelialmesenchymal transition via p53 in endometrial cancer" 18 : 204-215, 2020

5 Morikawa, T., "UBE2C is a marker of unfavorable prognosis in bladder cancer after radical cystectomy" 6 : 1367-1374, 2013

6 Chiang, A. J., "UBE2C drives human cervical cancer progression and is positively modulated by mTOR" 11 : 37-, 2021

7 Hegde, P. S., "Top 10 challenges in cancer immunotherapy" 52 : 17-35, 2020

8 Maynard, A., "Therapy-Induced Evolution of Human Lung Cancer Revealed by Single-Cell RNA Sequencing" 182 : 1232-1251, 2020

9 Cao, J., "The single-cell transcriptional landscape of mammalian organogenesis" 566 : 496-502, 2019

10 Pickup, M., "The roles of TGFbeta in the tumour microenvironment" 13 : 788-799, 2013

11 Trapnell, C., "The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells" 32 : 381-386, 2014

12 Mo, C. H., "The clinicopathological significance of UBE2C in breast cancer : a study based on immunohistochemistry, microarray and RNA-sequencing data" 17 : 83-, 2017

13 Nicholson, A. G., "The 2021 WHO classification of lung tumors: impact of advances since 2015." 17 : 362-387, 2022

14 Derynck, R., "TGFbeta biology in cancer progression and immunotherapy" 18 : 9-34, 2021

15 Wu, Y., "Spatiotemporal immune landscape of colorectal cancer liver metastasis at single-cell level" 12 : 134-153, 2022

16 Thrane, K., "Spatially resolved transcriptomics enables dissection of genetic heterogeneity in stage III cutaneous malignant melanoma" 78 : 5970-5979, 2018

17 Berglund, E., "Spatial maps of prostate cancer transcriptomes reveal an unexplored landscape of heterogeneity" 9 : 2419-, 2018

18 Young, M. D., "Single-cell transcriptomes from human kidneys reveal the cellular identity of renal tumors" 361 : 594-599, 2018

19 Lu, T., "Single-cell transcriptome atlas of lung adenocarcinoma featured with ground glass nodules" 6 : 69-, 2020

20 Wu, F., "Single-cell profiling of tumor heterogeneity and the microenvironment in advanced non-small cell lung cancer" 12 : 2540-, 2021

21 Qiu, X., "Single-cell mRNA quantification and differential analysis with Census" 14 : 309-315, 2017

22 Dong, R., "Single-cell characterization of malignant phenotypes and developmental trajectories of adrenal neuroblastoma" 38 : 716-733, 2020

23 Wang, Y., "Single-cell analysis reveals spatial heterogeneity of immune cells in lung adenocarcinoma" 9 : 638374-, 2021

24 Wu, Y. Y., "SCUBE3 is an endogenous TGF-beta receptor ligand and regulates the epithelial-mesenchymal transition in lung cancer" 30 : 3682-3693, 2011

25 Aibar, S., "SCENIC: single-cell regulatory network inference and clustering" 14 : 1083-1086, 2017

26 Cable, D. M., "Robust decomposition of cell type mixtures in spatial transcriptomics" 40 : 517-526, 2022

27 Qiu, X., "Reversed graph embedding resolves complex single-cell trajectories" 14 : 979-982, 2017

28 Sinjab, A., "Resolving the spatial and cellular architecture of lung adenocarcinoma by multiregion single-cell sequencing" 11 : 2506-2523, 2021

29 Lambrechts, D., "Phenotype molding of stromal cells in the lung tumor microenvironment" 24 : 1277-1289, 2018

30 Tumeh, P. C., "PD-1 blockade induces responses by inhibiting adaptive immune resistance" 515 : 568-571, 2014

31 Hafemeister, C., "Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression" 20 : 296-, 2019

32 Chen, Z., "Non-small-cell lung cancers : a heterogeneous set of diseases" 14 : 535-546, 2014

33 Browaeys, R., "NicheNet : modeling intercellular communication by linking ligands to target genes" 17 : 159-162, 2020

34 Wang, Y., "Myeloid-derived suppressor cells impair B cell responses in lung cancer through IL-7 and STAT5" 201 : 278-295, 2018

35 Sutherland, K. D., "Multiple cells-of-origin of mutant K-Ras-induced mouse lung adenocarcinoma" 111 : 4952-4957, 2014

36 Hu, X., "Multi-region exome sequencing reveals genomic evolution from preneoplasia to lung adenocarcinoma" 10 : 2978-, 2019

37 Li, Q., "Measuring reproducibility of highthroughput experiments" 5 : 1752-1779, 2011

38 Barresi, V., "Lymphatic vessel density and its prognostic value in stage I colorectal carcinoma" 64 : 6-12, 2011

39 Yotsukura, M., "Long-term prognosis of patients with resected adenocarcinoma in situ and minimally invasive adenocarcinoma of the lung" 16 : 1312-1320, 2021

40 Moncada, R., "Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas" 38 : 333-342, 2020

41 Kueckelhaus, J., "Inferring spatially transient gene expression pattern from spatial transcriptomic studies"

42 Huynh-Thu, V. A., "Inferring regulatory networks from expression data using tree-based methods" 5 : e12776-, 2010

43 Jin, S., "Inference and analysis of cell-cell communication using CellChat" 12 : 1088-, 2021

44 Mainardi, S., "Identification of cancer initiating cells in K-Ras driven lung adenocarcinoma" 111 : 255-260, 2014

45 Zhang, C., "Genomic landscape and immune microenvironment features of preinvasive and early invasive lung adenocarcinoma" 14 : 1912-1923, 2019

46 Li, Y., "Genomic characterisation of pulmonary subsolid nodules : mutational landscape and radiological features" 55 : 1901409-, 2020

47 Chen, H., "Genomic and immune profiling of pre-invasive lung adenocarcinoma" 10 : 5472-, 2019

48 Cortal, A., "Gene signature extraction and cell identity recognition at the single-cell level with Cell-ID" 39 : 1095-1102, 2021

49 Subramanian, A., "Gene set enrichment analysis : a knowledge-based approach for interpreting genome-wide expression profiles" 102 : 15545-15550, 2005

50 Hanzelmann, S., "GSVA: gene set variation analysis for microarray and RNA-seq data" 14 : 7-, 2013

51 Vincent, D. B., "Fast unfolding of communities in large networks" 2008 : 2008

52 Xu, X., "Evidence for type II cells as cells of origin of K-Ras-induced distal lung adenocarcinoma" 109 : 4910-4915, 2012

53 Sarode, P., "Epithelial cell plasticity defines heterogeneity in lung cancer" 65 : 109463-, 2020

54 Rowbotham, S. P., "Diverse cells at the origin of lung adenocarcinoma" 111 : 4745-4746, 2014

55 Zhang, M., "Dissecting transcriptional heterogeneity in primary gastric adenocarcinoma by single cell RNA sequencing" 70 : 464-475, 2021

56 Guo, J., "Deregulation of UBE2C-mediated autophagy repression aggravates NSCLC progression" 7 : 49-, 2018

57 Li, J., "Depletion of UBE2C reduces ovarian cancer malignancy and reverses cisplatin resistance via downregulating CDK1" 523 : 434-440, 2020

58 Xing, X., "Decoding the multicellular ecosystem of lung adenocarcinoma manifested as pulmonary subsolid nodules by single-cell RNA sequencing" 7 : eabd9738-, 2021

59 Armingol, E., "Deciphering cell-cell interactions and communication from gene expression" 22 : 71-88, 2021

60 Wang, Z., "Deciphering cell lineage specification of human lung adenocarcinoma with single-cell RNA sequencing" 12 : 6500-, 2021

61 Li, B., "Cumulus provides cloud-based data analysis for large-scale single-cell and single-nucleus RNA-seq" 17 : 793-798, 2020

62 Stuart, T., "Comprehensive integration of single-cell data" 177 : 1888-1902, 2019

63 Yang, Y., "Combined consideration of tumor-associated immune cell density and immune checkpoint expression in the peritumoral microenvironment for prognostic stratification of non-small-cell lung cancer patients" 13 : 811007-, 2022

64 Jia, M., "Clinicopathologic features and genetic alterations in adenocarcinoma in situ and minimally invasive adenocarcinoma of the lung : long-term follow-up study of 121 Asian Patients" 27 : 3052-3063, 2020

65 Teng, M. W., "Classifying cancers based on T-cell infiltration and PD-L1" 75 : 2139-2145, 2015

66 Irvine, A. F., "Characterising cancerassociated fibroblast heterogeneity in non-small cell lung cancer : a systematic review and meta-analysis" 11 : 3727-, 2021

67 Lou, E., "Cellular and molecular networking within the ecosystem of cancer cell communication via tunneling nanotubes" 6 : 95-, 2018

68 Pallasch, F. B., "Angiotensin inhibition, TGF-beta and EMT in cancer" 12 : 2785-, 2020

69 Martin-Gayo, E., "A reproducibility-based computational framework identifies an inducible, enhanced antiviral state in dendritic cells from HIV-1 elite controllers" 19 : 10-, 2018