Dataset Information

Single-Cell Transcriptomic Analysis of Ecosystems in Papillary Thyroid Carcinoma Progression

ABSTRACT:

Background

Despite 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.

Methods

Tissue 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.

Results

11,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.

Conclusions

Our 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.

SUBMITTER: Yan T

PROVIDER: S-EPMC8591202 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

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:BackgroundPapillary thyroid carcinoma (PTC) is the main pathological type of thyroid carcinoma (TC). Gender is a prominent background parameter for patients with PTC. Here, we aimed to delineate the differences in cell clusters and immune microenvironment in relation to gender in PTC.MethodsWe generated 6720, 14,666, and 33,373 single-cell transcriptomes that were pooled from the tissues of four male patients with PTC, seven female patients with PTC, and three patients with nodular goiter, respectively. We performed single-cell RNA-sequencing (scRNA-seq) based on BD Rhapsody and characterized the first single-cell transcriptomic landscape of PTC involving gender. The differential cell clusters and their gene profiles were identified and analyzed via a multi-resolution network in male and female patients. The interactions of fibroblasts and endothelial cells with malignant epithelial cells and the difference in the immune infiltration of B and T lymphocytes according to gender were assessed.ResultsMalignant epithelial cells were divided into two distinct subsets in male and female patients with PTC. Moreover, significant differences involving inferred copy-number variations (CNVs), gene profiles, and cell differentiation were detected between male and female patients. Regarding the interactions of fibroblasts and endothelial cells with malignant epithelial cells, members of the human leukocyte antigen (HLA) family and their receptors were considered as typical in female patients with PTC, while transforming growth factor beta 1 (TGFB1) and its receptors were typical of male patients with PTC. The characteristics of B cells, including cell clusters, cell differentiation, and dominant gene sets, were significantly different between genders.ConclusionsOur data revealed the detailed differences in cell clusters and immune microenvironment in PTC according to gender at the single-cell level, which provided new insights into the understanding of the impact of gender on PTC.

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.