An integrated single-cell and spatial transcriptomic atlas of thyroid cancer progression identifies prognostic fibroblast subpopulations
Matthew A. Loberg, George J. Xu, Sheau-Chiann Chen, Hua-Chang Chen, Claudia C. Wahoski, Kailey P. Caroland, Megan L. Tigue, Heather A. Hartmann, Jean-Nicolas Gallant, Courtney J. Phifer, Andres A. Ocampo, Dayle K. Wang, Reilly G. Fankhauser, Kirti A. Karunakaran, Chia-Chin Wu, Maxime Tarabichi, Sophia M. Shaddy, James L. Netterville, Sarah L. Rohde, Carmen C. Solórzano, Lindsay A. Bischoff, Naira Baregamian, Barbara A. Murphy, Jennifer H. Choe, Jennifer R. Wang, Eric C. Huang, Quanhu Sheng, Luciane T. Kagohara, Elizabeth M. Jaffee, Ryan H. Belcher, Ken S. Lau, Fei Ye, Ethan Lee, Vivian L. Weiss
Matthew A. Loberg, George J. Xu, Sheau-Chiann Chen, Hua-Chang Chen, Claudia C. Wahoski, Kailey P. Caroland, Megan L. Tigue, Heather A. Hartmann, Jean-Nicolas Gallant, Courtney J. Phifer, Andres A. Ocampo, Dayle K. Wang, Reilly G. Fankhauser, Kirti A. Karunakaran, Chia-Chin Wu, Maxime Tarabichi, Sophia M. Shaddy, James L. Netterville, Sarah L. Rohde, Carmen C. Solórzano, Lindsay A. Bischoff, Naira Baregamian, Barbara A. Murphy, Jennifer H. Choe, Jennifer R. Wang, Eric C. Huang, Quanhu Sheng, Luciane T. Kagohara, Elizabeth M. Jaffee, Ryan H. Belcher, Ken S. Lau, Fei Ye, Ethan Lee, Vivian L. Weiss
View: Text | PDF
Research Article Genetics Oncology

An integrated single-cell and spatial transcriptomic atlas of thyroid cancer progression identifies prognostic fibroblast subpopulations

Abstract

Although well-differentiated thyroid carcinoma (WDTC) is characterized by a robust treatment response, aggressive subtypes, such as anaplastic thyroid carcinoma (ATC), remain highly lethal. To understand thyroid cancer evolution in both children and adults, we analyzed single-cell transcriptomes of 423,733 cells from 81 samples and spatially resolved key tumor and microenvironment populations across 28 tumors with spatial transcriptomics, including rare and unique composite WDTC/ATC tumors and pediatric diffuse sclerosing thyroid carcinomas. Additionally, we identified gene signatures of stromal cell populations in 5 large thyroid cancer bulk RNA-sequencing cohorts. Through this multi-institutional effort, we defined a population of POSTN+ myofibroblast cancer-associated fibroblasts (myCAFs) that are intimately associated with invasive tumor cells and correlate with poor prognosis, lymph node metastasis, and disease progression in thyroid carcinoma. We also revealed a population of inflammatory CAFs that are distant to tumor cells and are found in the inflammatory stromal microenvironment of autoimmune thyroiditis. Together, our study provides spatial profiling of thyroid cancer evolution in samples with mixed WDTC/ATC histopathology and identifies a prognostic myCAF subtype with potential clinical utility in predicting aggressive disease in both children and adults.

Authors

Matthew A. Loberg, George J. Xu, Sheau-Chiann Chen, Hua-Chang Chen, Claudia C. Wahoski, Kailey P. Caroland, Megan L. Tigue, Heather A. Hartmann, Jean-Nicolas Gallant, Courtney J. Phifer, Andres A. Ocampo, Dayle K. Wang, Reilly G. Fankhauser, Kirti A. Karunakaran, Chia-Chin Wu, Maxime Tarabichi, Sophia M. Shaddy, James L. Netterville, Sarah L. Rohde, Carmen C. Solórzano, Lindsay A. Bischoff, Naira Baregamian, Barbara A. Murphy, Jennifer H. Choe, Jennifer R. Wang, Eric C. Huang, Quanhu Sheng, Luciane T. Kagohara, Elizabeth M. Jaffee, Ryan H. Belcher, Ken S. Lau, Fei Ye, Ethan Lee, Vivian L. Weiss

×

Figure 5

myCAF gene signatures are upregulated in malignant tumors across thyroid cancer bulk RNA-sequencing cohorts.

Options: View larger image (or click on image) Download as PowerPoint
myCAF gene signatures are upregulated in malignant tumors across thyroid...
(A) Corrplot showing Spearman’s rho correlations between single-sample GSVA (ssGSVA) scores for stromal subpopulations, ssGSVA MAP score, and expression of marker genes for stromal subpopulations in the Vanderbilt University/University of Washington bulk RNA-sequencing cohort (VUMC/UW cohort). Axes are ordered by hierarchical clustering. Boxes indicate hierarchical clustering groups. Significance levels indicate *P < 0.05, **P < 0.01, or ***P < 0.001. (B) Spearman’s rho correlation of myCAF ssGSVA score and a fibrosis score derived from H&E imaging data of TCGA PTCs (39). (C–E) Box plots showing (C) myCAF score or log2-transformed expression of CAF marker genes (D) FAP and (E) POSTN by diagnosis across 4 distinct bulk RNA-sequencing cohorts. P values calculated with Wilcoxon rank-sum test. Kruskal-Wallis test with subsequent pairwise Wilcoxon rank-sum tests with Bonferroni’s correction was used when comparing more than 2 groups. PDTC, poorly differentiated thyroid cancer; CNUH/SNUH, Chungnam/Seoul National University Hospitals.