Dataset Information

Single-cell transcriptomic analysis of the tumor ecosystems underlying initiation and progression of papillary thyroid carcinoma.

ABSTRACT: The tumor ecosystem of papillary thyroid carcinoma (PTC) is poorly characterized. Using single-cell RNA sequencing, we profile transcriptomes of 158,577 cells from 11 patients' paratumors, localized/advanced tumors, initially-treated/recurrent lymph nodes and radioactive iodine (RAI)-refractory distant metastases, covering comprehensive clinical courses of PTC. Our data identifies a "cancer-primed" premalignant thyrocyte population with normal morphology but altered transcriptomes. Along the developmental trajectory, we also discover three phenotypes of malignant thyrocytes (follicular-like, partial-epithelial-mesenchymal-transition-like, dedifferentiation-like), whose composition shapes bulk molecular subtypes, tumor characteristics and RAI responses. Furthermore, we uncover a distinct BRAF-like-B subtype with predominant dedifferentiation-like thyrocytes, enriched cancer-associated fibroblasts, worse prognosis and promising prospect of immunotherapy. Moreover, potential vascular-immune crosstalk in PTC provides theoretical basis for combined anti-angiogenic and immunotherapy. Together, our findings provide insight into the PTC ecosystem that suggests potential prognostic and therapeutic implications.

SUBMITTER: Pu W

PROVIDER: S-EPMC8523550 | biostudies-literature | 2021 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Single-cell transcriptomic analysis of the tumor ecosystems underlying initiation and progression of papillary thyroid carcinoma.

Pu Weilin W Shi Xiao X Yu Pengcheng P Zhang Meiying M Liu Zhiyan Z Tan Licheng L Han Peizhen P Wang Yu Y Ji Dongmei D Gan Hualei H Wei Wenjun W Lu Zhongwu Z Qu Ning N Hu Jiaqian J Hu Xiaohua X Luo Zaili Z Li Huajun H Ji Qinghai Q Wang Jiucun J Zhang Xiaoming X Wang Yu-Long YL

Nature communications 20211018 1

The tumor ecosystem of papillary thyroid carcinoma (PTC) is poorly characterized. Using single-cell RNA sequencing, we profile transcriptomes of 158,577 cells from 11 patients' paratumors, localized/advanced tumors, initially-treated/recurrent lymph nodes and radioactive iodine (RAI)-refractory distant metastases, covering comprehensive clinical courses of PTC. Our data identifies a "cancer-primed" premalignant thyrocyte population with normal morphology but altered transcriptomes. Along the dev ...[more]

PMID: 34663816

Similar Datasets

Project description:BackgroundDespite extensive research, the papillary thyroid carcinoma (PTC) ecosystem is poorly characterized and, in particular, locoregional progression. Available evidence supports that single-cell transcriptome sequencing (Sc-RNA seq) can dissect tumor ecosystems.MethodsTissue samples from one PTC patient, including matched primary tumor (Ca), lymph node (LN) metastasis, and paracancerous tissue (PCa), were subjected to Sc-RNA seq with 10×Genomics. Dual-label immunofluorescence and immunohistochemistry were used to confirm the existence of cell subtypes in a separate cohort.Results11,805 cell transcriptomes were profiled, cell landscapes of PTC were composed of malignant follicular epithelial cells (MFECs), CD8+ and CD4+ T cells, B cells, vascular endothelial cells, fibroblasts and cancer-associated fibroblasts (CAFs). Between Ca and LN ecosystems, the proportions of MFEC and interstitial cells were similar, less than 1/25(229/6,694, 361/3,895), while the proportion of normal follicular epithelial cells (NFECs) and interstitial cells was > 2 in PCa (455/171). NFECs in PCa formed a separate cluster, while MFECs in Ca and LN exhibited a profound transcriptional overlap, and the interstitial cells among these samples had an overall concordance in their identity and representation, albeit with some distinctions in terms of the cell percentage per subset. A fraction of the B cell subpopulation in Ca expressed inhibitory receptors, while their respective ligand genes were clearly transcribed in T cell and malignant epithelial cell clusters, while some CD8+ T cells in both Ca and LN produced high levels of inhibitory receptors whose respective ligands were overexpressed in some CD4+ T cells. Three CAF subtypes in Ca and LN were identified, which may be due to mutual transitions.ConclusionsOur data provide new insights into the PTC ecosystem and highlight the differences in ecosystems in PTC progression, which updates our understanding of PTC biology and will improve individualized patient treatment.

Project description:BackgroundThe tumor microenvironment (TME) plays a pivotal role in cancer progression in papillary thyroid carcinoma (PTC), yet the composition and the phenotype of cells within the TME in bilateral PTC are poorly understood.MethodsWe performed unbiased transcriptome-wide single-cell RNA sequencing (scRNA-seq) analysis on 29,561 cells from 3 pairs of bilateral PTC and 1 non-tumor thyroid sample. The results of the analysis were validated by a large-scale bulk transcriptomic dataset deposited in The Cancer Genome Atlas (TCGA) database.ResultsOur integrative analysis of thyroid follicular cells revealed 42 signaling pathways enriched in malignant follicular cells, including cytokine-cytokine receptor interaction, PI3K/Akt signaling pathway, mitogen-activated protein kinase (MAPK) signaling pathway, and tumor necrosis factor (TNF) signaling pathway. A 6-gene signature (CXCL3, CXCL1, IL1A, CCL5, TNFRSF12A, and IL18) in the cytokine-cytokine receptor interaction pathway was constructed to predict the prognosis of patients with PTC, with high risk scores being associated with decreased overall survival [hazard ratio (HR) = 3.863, 95% CI = 2.233-6.682, p < 0.001]. Gene set variation analysis (GSVA) indicated that the pathways enriched in bilateral PTC were significantly different, indicating great heterogeneity in bilateral PTC, even with the same BRAF V600E mutation. Comprehensive analysis of T cells revealed that the proportion of CD8+ tissue-resident memory T cells expressing IFNG decreased in tumor samples with advanced N stage. Within the myeloid compartment, the ratio of suppressive M2-like to pro-inflammatory M1-like macrophages increased with advanced disease stage, which was confirmed in the bulk dataset using transcriptomic profiles. In addition, we also identified numerous biologically critical interactions among myeloid cells, T cells, and follicular cells, which were related to T-cell recruitment, M2-like macrophage polarization, malignant follicular cell progression, and T-cell inhibitory signaling.ConclusionOur integrative analyses revealed great inter-tumor heterogeneity within the TME in bilateral PTC, which will offer assistance for precise diagnosis and treatment.

Project description:Papillary thyroid carcinoma (PTC) is the fastest-growing disease caused by numerous molecular alterations in addition to previously reported DNA mutations. There is a compelling need to identify novel transcriptomic alterations that are associated with the pathogenesis of PTC with potential diagnostic and prognostic implications. Methods: We gathered and compared 242 expression profiles between paired PTC and adjacent normal tissues and identified and validated the coding and long non-coding RNAs (lncRNAs) associated with the extrathyroidal extension (ETE) of 655 PTC patients in two independent cohorts, followed by predicting their interactions with drugs. Co-expression, RNA interaction, Kaplan-Meier survival and multivariate Cox proportional regression analyses were performed to identify dysregulated lncRNAs and genes that correlated with clinical outcomes of PTC. Alternative splicing (AS), RNA circularization, and editing were also compared between transcriptomes to expand the repertoire of molecular alterations in PTC. Results: Numerous genes related to cellular microenvironment and steroid hormone response were associated with the ETE of PTC. Drug susceptibility predictions of the expression signature revealed two highly ranked compounds, 6-bromoindirubin-3'-oxime and lovastatin. Co-expression and RNA interaction analysis revealed the essential role of lncRNAs in PTC pathogenesis by modulating extracellular matrix and cell adhesion. Eight genes and two novel lncRNAs were identified that correlated with the aggressive nature and disease-free survival of PTC. Furthermore, this study provided the transcriptome-wide landscape of circRNAs in PTC and uncovered dissimilar expression profiles among circRNAs originating from the same host gene, suggesting the functional complexity of circRNAs in PTC carcinogenesis. The newly identified AS events in the SERPINA1 and FN1 genes may improve the sensitivity and specificity of these diagnostic biomarkers. Conclusions: Our study uncovered a comprehensive transcriptomic signature associated with the carcinogenesis and aggressive behavior of PTC, as well as presents a catalog of 10 potential biomarkers, which would facilitate PTC prognosis and development of new therapeutic strategies for this cancer.

Dataset Information

Single-cell transcriptomic analysis of the tumor ecosystems underlying initiation and progression of papillary thyroid carcinoma.

Publications

Single-cell transcriptomic analysis of the tumor ecosystems underlying initiation and progression of papillary thyroid carcinoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets