Analysis of the Mutational Landscape of Osteosarcomas Identifies Genes Related to Metastasis and Prognosis and Disrupted Biological Pathways of Immune Response and Bone Development
Authors:
- Sara Ferreira PiresSPSara Ferreira Pires
- Juliana Sobral de BarrosJBJuliana Sobral de Barros
- Silvia Souza da CostaSCSilvia Souza da Costa
- Gabriel Bandeira Do CarmoGCGabriel Bandeira Do Carmo
- Marília de Oliveira ScliarMSMarília de Oliveira Scliar
- André van Helvoort LengertALAndré van Helvoort Lengert
- Érica BoldriniÉBÉrica Boldrini
- Sandra Regini Morini da SilvaSSSandra Regini Morini da Silva
- Daniel Onofre VidalDVDaniel Onofre Vidal
- Mariana MaschiettoMariana Maschietto
- Ana Cristina Victorino KrepischiAKAna Cristina Victorino Krepischi
Open Access
- 21 June 2023
- journal article
- research article
- Published by MDPI AG in International Journal of Molecular Sciences
- Vol. 24 (13) , 10463
Abstract
Osteosarcoma (OS) is the most prevalent type of bone tumor, but slow progress has been achieved in disentangling the full set of genomic events involved in its initiation and progression. We assessed by NGS the mutational spectrum of 28 primary OSs from Brazilian patients, and identified 445 potentially deleterious SNVs/indels and 1176 copy number alterations (CNAs). TP53 was the most recurrently mutated gene, with an overall rate of ~60%, considering SNVs/indels and CNAs. The most frequent CNAs (~60%) were gains at 1q21.2q21.3, 6p21.1, and 8q13.3q24.22, and losses at 10q26 and 13q14.3q21.1. Seven cases presented CNA patterns reminiscent of complex events (chromothripsis and chromoanasynthesis). Putative RB1 and TP53 germline variants were found in five samples associated with metastasis at diagnosis along with complex genomic patterns of CNAs. PTPRQ, KNL1, ZFHX4, and DMD alterations were prevalent in metastatic or deceased patients, being potentially indicative of poor prognosis. TNFRSF11B, involved in skeletal system development and maintenance, emerged as a candidate for osteosarcomagenesis due to its biological function and a high frequency of copy number gains. A protein–protein network enrichment highlighted biological pathways involved in immunity and bone development. Our findings reinforced the high genomic OS instability and heterogeneity, and led to the identification of novel disrupted genes deserving further evaluation as biomarkers due to their association with poor outcomes.Keywords
Funding Information
- Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (88887.816566/2023-00, 305806/2019-0, 311141/2021-8, 140258/2020-5, 2013/08028-1, 2014/10250-7, 2018/21047-9)
- Conselho Nacional de Desenvolvimento Científico e Tecnológico (88887.816566/2023-00, 305806/2019-0, 311141/2021-8, 140258/2020-5, 2013/08028-1, 2014/10250-7, 2018/21047-9)
- Fundação de Amparo à Pesquisa do Estado de São Paulo (88887.816566/2023-00, 305806/2019-0, 311141/2021-8, 140258/2020-5, 2013/08028-1, 2014/10250-7, 2018/21047-9)
This publication has 87 references indexed in Scilit:
- The Cancer Genome Atlas Pan-Cancer analysis projectNature Genetics, 2013
- The Genetics of OsteosarcomaSarcoma, 2012
- The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics DataCancer Discovery, 2012
- Detecting false-positive signals in exome sequencingHuman Mutation, 2012
- Using Epidemiology and Genomics to Understand Osteosarcoma EtiologySarcoma, 2011
- Copy Number Variation in Human Health, Disease, and EvolutionAnnual Review of Genomics and Human Genetics, 2009
- The Sequence Alignment/Map format and SAMtoolsBioinformatics, 2009
- Fast and accurate short read alignment with Burrows–Wheeler transformBioinformatics, 2009
- Osteosarcoma incidence and survival rates from 1973 to 2004Cancer, 2009
- Functions of RANKL/RANK/OPG in bone modeling and remodelingArchives of Biochemistry and Biophysics, 2008