Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Feb 21;11(1):1008.
doi: 10.1038/s41467-020-14646-w.

Immuno-genomic landscape of osteosarcoma

Chia-Chin Wu #  1 , Hannah C Beird #  1 , J Andrew Livingston #  2   3 , Shailesh Advani  3 , Akash Mitra  1 , Shaolong Cao  4 , Alexandre Reuben  5 , Davis Ingram  6 , Wei-Lien Wang  6 , Zhenlin Ju  4 , Cheuk Hong Leung  7 , Heather Lin  7 , Youyun Zheng  1 , Jason Roszik  1 , Wenyi Wang  4 , Shreyaskumar Patel  2 , Robert S Benjamin  2 , Neeta Somaiah  2 , Anthony P Conley  2 , Gordon B Mills  8 , Patrick Hwu  2 , Richard Gorlick  3 , Alexander Lazar  6 , Najat C Daw  3 , Valerae Lewis  9 , P Andrew Futreal  10
Affiliations

Immuno-genomic landscape of osteosarcoma

Chia-Chin Wu et al. Nat Commun. .

Abstract

Limited clinical activity has been seen in osteosarcoma (OS) patients treated with immune checkpoint inhibitors (ICI). To gain insights into the immunogenic potential of these tumors, we conducted whole genome, RNA, and T-cell receptor sequencing, immunohistochemistry and reverse phase protein array profiling (RPPA) on OS specimens from 48 pediatric and adult patients with primary, relapsed, and metastatic OS. Median immune infiltrate level was lower than in other tumor types where ICI are effective, with concomitant low T-cell receptor clonalities. Neoantigen expression in OS was lacking and significantly associated with high levels of nonsense-mediated decay (NMD). Samples with low immune infiltrate had higher number of deleted genes while those with high immune infiltrate expressed higher levels of adaptive resistance pathways. PARP2 expression levels were significantly negatively associated with the immune infiltrate. Together, these data reveal multiple immunosuppressive features of OS and suggest immunotherapeutic opportunities in OS patients.

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing interests.

Figures

Fig. 1
Fig. 1. Clinical associations with somatic alterations in OS.
a, b Kaplan-Meier survival analysis (disease-free survival for a; overall survival for b) of mutation signature 8 scores which describe the contribution of mutation signature 8 to the point mutation profile. The cohort was separated into tertiles based on signature 8 scores. Group 1 (denoted as T1) is the first tertile with lowest mutation signature 8 scores. Groups 2 and 3 (denoted as T2 and T3 respctively) are the next 2 tertiles with higher mutation signature 8 scores. c Hierarchical clustering heatmap of rearrangements classification pattern identified in our cohort along with clinical features: age of diagnosis (<18 years, 18– < 50 years, and >50 years), tumor specimen type (primary, local recurrence, metastasis), and vital status (dead, alive). Structural rearrangements were classified based on their type: deletion (del), tandem duplication (tds), inversion (inv), and interchromosomal translocation (trans), and size, and clustered (denoted as c_) versus non-clustered events. d Boxplot showing the ratio of clustered rearrangements within associated with chromothriptic regions as compared across three groups of specimens based on age of diagnosis (<18 years, 18–<50 years, and >50 years). Then significance values (P-values) of the comparisons are from the Wilcoxon rank sum test are shown in asterisks. e Pearson correlation between normalized telomere length and number of copy number segments of samples. f Pearson correlation between normalized telomere length and ATRX gene expression level (in log2 scale) of samples. Samples with both ATRX and TP53 alterations, ATRX alterations alone, and TP53 alterations alone were respectively marked as red, blue, and green color.
Fig. 2
Fig. 2. Expressed neoantigens and rearrangements.
a Stacked bar chart for number of nonsynonymous mutations and neoantigens detected by WGS and capture-based RNASeq. (expressed neoantigen: nonsynonymous mutations found in both WGS and capture-based RNASeq that are predicted to be a neoantigen; expressed non-neoantigen: nonsynonymous mutations found in both WGS and capture-based RNASeq that are not predicted to be a neoantigen; unexpressed: found in WGS only). b Overlapping rearrangements between WGS and RNASeq based rearrangements where both partners are within coding regions (unexpressed: found in WGS only; expressed:found in both WGS and RNASeq).
Fig. 3
Fig. 3. Immune profiling.
a ESTIMATE immune scores in TCGA tumor types, our OS cohort (MDACC.OS), the OS cohort from ICGC (BOCA-UK.OS), the OS cohort from TARGET (TARGET.OS), and from four patients who underwent combination CTLA-4 blockade and PD-L1 blockade but exhibited no objective responses (Immunotherapy.OS). Note that TCGA-THYM is derived from the thymus made up mostly of lymphocytes and would have inflated scores. b. Unsupervised hierarchical clustering based on ssGSEA enrichment scores of each immune gene list from Charoentong et al. The three predominant clusters are referred to as “C1” and “C2”, and “C3” immune infiltrate groups. c Boxplot of the ESTIMATE immune scores of samples across the three immune clusters C1, C2, and C3. The significances of the comparisons were from the Wilcoxon rank sum test are shown in asterisks. d The geometric mean expression of GZMA and PRF1 genes (cytolytic score) shown as a boxplot across the three immune clusters C1, C2, and C3. The significances of the comparisons were from the Wilcoxon rank sum test. e Immunosuppressive pathways that are significantly deregulated between immune clusters C1 and C3. The significance (FDR) of deregulation were from the GSEA analysis. f Fold change and significance of expression difference between C1 and C3 for key genes from pathways in e: CD274(PD-L1), CTLA4, IDO1, IFNG, IFNGR1, and IL6 (these genes were not present in the gene lists used to define these immune groups). The linear mode in the limma R package was used to determine significance.
Fig. 4
Fig. 4. Copy number alterations and immune score associations.
a, b Pearson correlation of the number of deleted genes and immune scores generated from ESTIMATE for our OS cohort (a) and the TARGET cohort (b). c Boxplot of PARP2 copy number levels (in log2 level) in each of the immune infiltrate groups C1, C2, and C3 (as decribed in Fig. 3b). Wilcoxon rank sum test was used to determine significance. d Boxplot of PARP2 gene expression levels (in log2 level) in each of the immune infiltrate groups C1, C2, and C3 (as decribed in Fig. 3b). Wilcoxon rank sum test was used to determine significance.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Mirabello L, Troisi RJ, Savage SA. Osteosarcoma incidence and survival rates from 1973 to 2004: data from the Surveillance, Epidemiology, and End Results Program. Cancer. 2009;115:1531–1543. doi: 10.1002/cncr.24121. - DOI - PMC - PubMed
    1. Bielack SS, et al. Prognostic factors in high-grade osteosarcoma of the extremities or trunk: an analysis of 1,702 patients treated on neoadjuvant cooperative osteosarcoma study group protocols. J. Clin. Oncol. 2002;20:776–790. doi: 10.1200/JCO.20.3.776. - DOI - PubMed
    1. Meyers PA, et al. Osteogenic sarcoma with clinically detectable metastasis at initial presentation. J. Clin. Oncol. 1993;11:449–453. doi: 10.1200/JCO.1993.11.3.449. - DOI - PubMed
    1. Kager L, et al. Primary metastatic osteosarcoma: presentation and outcome of patients treated on neoadjuvant Cooperative Osteosarcoma Study Group protocols. J. Clin. Oncol. 2003;21:2011–2018. doi: 10.1200/JCO.2003.08.132. - DOI - PubMed
    1. Chen X, et al. Recurrent somatic structural variations contribute to tumorigenesis in pediatric osteosarcoma. Cell Rep. 2014;7:104–112. doi: 10.1016/j.celrep.2014.03.003. - DOI - PMC - PubMed

Publication types

MeSH terms

Substances