Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in Patients With Osteosarcoma

Importance

Osteosarcoma, the most common malignant bone tumor in children and adolescents, occurs in a high number of cancer predisposition syndromes that are defined by highly penetrant germline mutations. The germline genetic susceptibility to osteosarcoma outside of familial cancer syndromes remains unclear.

Objective

To investigate the germline genetic architecture of 1244 patients with osteosarcoma.

Design, Setting, and Participants

Whole-exome sequencing (n = 1104) or targeted sequencing (n = 140) of the DNA of 1244 patients with osteosarcoma from 10 participating international centers or studies was conducted from April 21, 2014, to September 1, 2017. The results were compared with the DNA of 1062 individuals without cancer assembled internally from 4 participating studies who underwent comparable whole-exome sequencing and 27 173 individuals of non-Finnish European ancestry who were identified through the Exome Aggregation Consortium (ExAC) database. In the analysis, 238 high-interest cancer-susceptibility genes were assessed followed by testing of the mutational burden across 736 additional candidate genes. Principal component analyses were used to identify 732 European patients with osteosarcoma and 994 European individuals without cancer, with outliers removed for patient-control group comparisons. Patients were subsequently compared with individuals in the ExAC group. All data were analyzed from June 1, 2017, to July 1, 2019.

Main Outcomes and Measures

The frequency of rare pathogenic or likely pathogenic genetic variants.

Results

Among 1244 patients with osteosarcoma (mean [SD] age at diagnosis, 16 [8.9] years [range, 2-80 years]; 684 patients [55.0%] were male), an analysis restricted to individuals with European ancestry indicated a significantly higher pathogenic or likely pathogenic variant burden in 238 high-interest cancer-susceptibility genes among patients with osteosarcoma compared with the control group (732 vs 994, respectively; P = 1.3 × 10-18). A pathogenic or likely pathogenic cancer-susceptibility gene variant was identified in 281 of 1004 patients with osteosarcoma (28.0%), of which nearly three-quarters had a variant that mapped to an autosomal-dominant gene or a known osteosarcoma-associated cancer predisposition syndrome gene. The frequency of a pathogenic or likely pathogenic cancer-susceptibility gene variant was 128 of 1062 individuals (12.1%) in the control group and 2527 of 27 173 individuals (9.3%) in the ExAC group. A higher than expected frequency of pathogenic or likely pathogenic variants was observed in genes not previously linked to osteosarcoma (eg, CDKN2A, MEN1, VHL, POT1, APC, MSH2, and ATRX) and in the Li-Fraumeni syndrome-associated gene, TP53.

Conclusions and Relevance

In this study, approximately one-fourth of patients with osteosarcoma unselected for family history had a highly penetrant germline mutation requiring additional follow-up analysis and possible genetic counseling with cascade testing.

Low-frequency variation near common germline susceptibility loci are associated with risk of Ewing sarcoma

BACKGROUND

Ewing sarcoma (EwS) is a rare, aggressive solid tumor of childhood, adolescence and young adulthood associated with pathognomonic EWSR1-ETS fusion oncoproteins altering transcriptional regulation. Genome-wide association studies (GWAS) have identified 6 common germline susceptibility loci but have not investigated low-frequency inherited variants with minor allele frequencies below 5% due to limited genotyped cases of this rare tumor.

METHODS

We investigated the contribution of rare and low-frequency variation to EwS susceptibility in the largest EwS genome-wide association study to date (733 EwS cases and 1,346 unaffected controls of European ancestry).

RESULTS

We identified two low-frequency variants, rs112837127 and rs2296730, on chromosome 20 that were associated with EwS risk (OR = 0.186 and 2.038, respectively; P-value < 5×10-8) and located near previously reported common susceptibility loci. After adjusting for the most associated common variant at the locus, only rs112837127 remained a statistically significant independent signal (OR = 0.200, P-value = 5.84×10-8).

CONCLUSIONS

These findings suggest rare variation residing on common haplotypes are important contributors to EwS risk.

IMPACT

Motivate future targeted sequencing studies for a comprehensive evaluation of low-frequency and rare variation around common EwS susceptibility loci.

Cytotoxic T lymphocyte-associated antigen 4 rs231775 polymorphism and osteosarcoma

Systematic reviews on osteosarcoma have concluded that CTLA4 rs231775 AA genotype influences risk for the disease in the Chinese population. Remarkably, rs231775 shows different frequencies in different human populations. Therefore, it would be interesting to know whether this SNP is related to the risk of osteosarcoma in other populations. The present study aimed to evaluate the association between rs231775 and the susceptibility of osteosarcoma in the Spanish population. We performed an updated meta-analysis including a total of 538 cases and 623 controls. The genotypic association analyses showed that the CTLA4 rs231755 was associated with osteosarcoma susceptibility in the Spanish population. When meta-analysis was performed, the results displayed that CTLA4 rs231775 AA genotype was associated with the risk of developing osteosarcoma in all analyzed populations (OR=2.07; 95% CI: 1.48-2.89).The rs231775 AA genotype could be considered as a susceptibility marker in osteosarcoma Keywords: CTLA4, rs231775, +49A/G and osteogenic sarcoma.

Abstract A13: Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility

Ewing sarcoma (EWS), a pediatric tumor predominantly occurring in children of European ancestry, is characterized by the EWSR1-FLI1 fusion oncogene. We performed a genome-wide association study (GWAS) of 749 EWS cases and 1,378 unaffected individuals of European ancestry. Our study replicated previously reported susceptibility loci at 1p36.22, 10q21.3, and 15q15.1 and identified new loci at 6p25.1, 8q24.23, 20p11.22, and 20p11.23 (P-values &lt;5x10-8). Effect estimates exhibited odds ratios (ORs) in excess of 1.7, which is high for cancer GWAS, and striking in light of the rarity of EWS cases in familial cancer syndromes. In expression quantitative trait locus (eQTL) analyses, we observed altered expression patterns for plausible candidate genes at 6p25.1 with RREB1, a RAS-responsive element, and at 20p11.23 with KIZ, a centrosomal stabilization protein. The 20p11.22 locus is also near NKX2-2, a highly overexpressed gene in EWS. Interestingly, most loci reside near GGAA repeat sequences where binding of the EWSR1-FLI1 fusion protein occurs. The seven EWS susceptibility loci discovered in only 749 cases make EWS one of the most productive GWAS studied cancers when considering a locus-to-case discovery ratio. The substantive estimated ORs suggest a distinctive underlying genetic architecture for EWS in which moderate-risk SNPs likely constitute a significant fraction and may interact with EWSR1-FLI1 binding.

Citation Format: Mitchell J. Machiela, Thomas G.P. Grünewald, Didier Surdez, Stephanie Reynaud, Olivier Mirabeau, Eric Karlins, Rebeca Alba Rubio, Sakina Zaidi, Sandrine Grossetete-Lalami, Stelly Ballet, Eve Lapouble, Valérie Laurence, Jean Michon, Gaelle Pierron, Heinrich Kovar, Nathalie Gaspar, Udo Kontny, Anna González-Neira, Piero Picci, Javier Alonso, Ana Patino-Garcia, Nadège Corradini, Neal D. Freedman, Nathaniel Rothman, Casey L. Dagnall, Laurie Burdett, Kristine Jones, Michelle Manning, Kathleen Wyatt, Weiyin Zhou, Meredith Yeager, David G. Cox, Robert N. Hoover, Javed Khan, Gregory T. Armstrong, Wendy M. Leisenring, Smita Bhatia, Leslie L. Robison, Uta Dirksen, Markus Metzler, Wolfgang Hartmann, Konstantin Strauch, Thomas Kirchner, Andreas E. Kulozik, Lindsay M. Morton, Lisa Mirabello, Margaret A. Tucker, Franck Tirode, Stephen Chanock, Olivier Delattre. Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility [abstract]. In: Proceedings of the AACR Special Conference: Pediatric Cancer Research: From Basic Science to the Clinic; 2017 Dec 3-6; Atlanta, Georgia. Philadelphia (PA): AACR; Cancer Res 2018;78(19 Suppl):Abstract nr A13.

Abstract 2970: Multiple new susceptibility loci identified in genome-wide association study of Ewing sarcoma

Ewing sarcoma (EWS) is a rare pediatric tumor predominantly occurring in children of European ancestry and is characterized by the pathognomonic EWSR1-FLI1 fusion oncogene. To identify germline susceptibility loci associated with EWS risk, we performed a genome-wide association study (GWAS) meta-analysis of 749 EWS cases and 1,378 unaffected individuals of European ancestry from sample collections within the Institut Curie, National Cancer Institute and the Childhood Cancer Survivor Study. Our study replicated previously reported susceptibility loci at 1p36.22, 10q21.3 and 15q15.1 as well as identified novel loci at 6p25.1, 8q24.23, 20p11.22 and 20p11.23 (P-values&lt;5×10-8). These seven EWS susceptibility loci discovered in only 749 cases make EWS one of the most productive GWAS studied cancers when considering a locus to case discovery ratio. All estimated effect estimates were high for cancer GWAS with odds ratios in excess of 1.7 observed. These high per allele effects among relatively common germline variants are striking in light of the rarity of EWS cases and lack of evidence of EWS as part of a familial cancer syndrome and therefore suggest a distinctive genetic architecture for EWS. Interestingly, in silico bioinformatics analysis identified that most EWS susceptibility loci reside near GGAA nucleotide repeat sequences where binding of the EWSR1-FLI1 fusion protein occurs. ChIP-seq analyses confirmed in vivo binding of EWSR1-FLI1, suggesting germline variation in these regions could alter EWSR1-FLI1 binding and potentially deregulate neighboring genes. To identify genes with allele specific expression differences, we carried out expression quantitative trait locus (eQTL) analyses at each identified EWS susceptibility locus. We identified eQTLs for plausible candidate genes at 6p25.1 with RREB1, a RAS-responsive element, and at 20p11.23 with KIZ, a centrosomal stabilization protein. We also noted the 20p11.22 locus is near NKX2-2, a highly overexpressed gene in EWS, although no eQTL was observed in our expression data. Furthermore, knockdown of EWSR1-FLI1 in EWS cell lines indicated a more than 2-fold difference in expression of RREB1 and NKX2-2, further supporting the role of specific regulation of these genes by EWSR1-FLI1 and suggesting RREB1 and NKX2-2 may be transcription factors involved in core regulatory circuitries of EWS. Overall, our study suggests a distinctive underlying genetic architecture for EWS in which moderate risk common germline variants interact with EWSR1-FLI1 binding to alter expression of nearby target genes.

Citation Format: Mitchell J. Machiela, Thomas G. Grünewald, Didier Surdez, Stephanie Reynaud, Olivier Mirabeau, Eric Karlins, Rebeca Alba Rubio, Sakina Zaidi, Sandrine Grossetete-Lalami, Stelly Ballet, Eve Lapouble, Valérie Laurence, Jean Michon, Gaelle Pierron, Heinrich Kovar, Nathalie Gaspar, Udo Kontny, Anna González-Neira, Piero Picci, Javier Alonso, Ana Patino-Garcia, Nadège Corradini, Neal D. Freedman, Nathaniel Rothman, Casey L. Dagnall, Laurie Burdette, Kristine Jones, Michelle Manning, Kathleen Wyatt, Weiyin Zhou, Meredith Yeager, David G. Cox, Robert N. Hoover, Javed Khan, Gregory T. Armstrong, Wendy M. Leisenring, Smita Bhatia, Leslie L. Robison, Uta Dirksen, Markus Metzler, Wolfgang Hartmann, Konstantin Strauch, Thomas Kirchner, Andreas E. Kulozik, Lindsay M. Morton, Lisa Mirabello, Margaret A. Tucker, Franck Tirode, Stephen J. Chanock, Olivier Delattre. Multiple new susceptibility loci identified in genome-wide association study of Ewing sarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2970.

Influence of genetic variants in TPMT and COMT associated with cisplatin induced hearing loss in patients with cancer: two new cohorts and a meta-analysis reveal significant heterogeneity between cohorts

Treatment with cisplatin-containing chemotherapy regimens causes hearing loss in 40–60% of cancer patients. It has been suggested that genetic variants in the genes encoding thiopurine S-methyltransferase (TPMT) and catechol O-methyltransferase (COMT) can predict the development of cisplatin-induced ototoxicity and may explain interindividual variability in sensitivity to cisplatin-induced hearing loss. Two recently published studies however, sought to validate these findings and showed inconsistent results. The aim of this study was to evaluate the role of polymorphisms in the TPMT and COMT genes in cisplatin-induced ototoxicity. Therefore we investigated two independent cohorts of 110 Dutch and 38 Spanish patients with osteosarcoma and performed a meta-analysis including all previously published studies resulting in a total population of 664 patients with cancer. With this largest meta-analysis performed to date, we show that the influence of TPMT and COMT on the development of cisplatin-induced hearing loss may be less important than previously suggested.

Imputation and subset-based association analysis across different cancer types identifies multiple independent risk loci in the TERT-CLPTM1L region on chromosome 5p15.33

Genome-wide association studies (GWAS) have mapped risk alleles for at least 10 distinct cancers to a small region of 63 000 bp on chromosome 5p15.33. This region harbors the TERT and CLPTM1L genes; the former encodes the catalytic subunit of telomerase reverse transcriptase and the latter may play a role in apoptosis. To investigate further the genetic architecture of common susceptibility alleles in this region, we conducted an agnostic subset-based meta-analysis (association analysis based on subsets) across six distinct cancers in 34 248 cases and 45 036 controls. Based on sequential conditional analysis, we identified as many as six independent risk loci marked by common single-nucleotide polymorphisms: five in the TERT gene (Region 1: rs7726159, P = 2.10 × 10(-39); Region 3: rs2853677, P = 3.30 × 10(-36) and PConditional = 2.36 × 10(-8); Region 4: rs2736098, P = 3.87 × 10(-12) and PConditional = 5.19 × 10(-6), Region 5: rs13172201, P = 0.041 and PConditional = 2.04 × 10(-6); and Region 6: rs10069690, P = 7.49 × 10(-15) and PConditional = 5.35 × 10(-7)) and one in the neighboring CLPTM1L gene (Region 2: rs451360; P = 1.90 × 10(-18) and PConditional = 7.06 × 10(-16)). Between three and five cancers mapped to each independent locus with both risk-enhancing and protective effects. Allele-specific effects on DNA methylation were seen for a subset of risk loci, indicating that methylation and subsequent effects on gene expression may contribute to the biology of risk variants on 5p15.33. Our results provide strong support for extensive pleiotropy across this region of 5p15.33, to an extent not previously observed in other cancer susceptibility loci.

Common variants near TARDBP and EGR2 are associated with susceptibility to Ewing sarcoma

Ewing sarcoma, a pediatric tumor characterized by EWSR1-ETS fusions, is predominantly observed in populations of European ancestry. We performed a genome-wide association study (GWAS) of 401 French individuals with Ewing sarcoma, 684 unaffected French individuals and 3,668 unaffected individuals of European descent and living in the United States. We identified candidate risk loci at 1p36.22, 10q21 and 15q15. We replicated these loci in two independent sets of cases and controls. Joint analysis identified associations with rs9430161 (P = 1.4 × 10(-20); odds ratio (OR) = 2.2) located 25 kb upstream of TARDBP, rs224278 (P = 4.0 × 10(-17); OR = 1.7) located 5 kb upstream of EGR2 and, to a lesser extent, rs4924410 at 15q15 (P = 6.6 × 10(-9); OR = 1.5). The major risk haplotypes were less prevalent in Africans, suggesting that these loci could contribute to geographical differences in Ewing sarcoma incidence. TARDBP shares structural similarities with EWSR1 and FUS, which encode RNA binding proteins, and EGR2 is a target gene of EWSR1-ETS. Variants at these loci were associated with expression levels of TARDBP, ADO (encoding cysteamine dioxygenase) and EGR2.

Abstract 5748: Height, birth-weight, and osteosarcoma risk in a large meta-analysis

Osteosarcoma (OS) is the most commonly diagnosed primary malignancy of bone. Its primary incidence peak occurs in adolescents and young adults, suggesting that perinatal factors and/or pubertal growth may have a role in etiology. In young patients, OS occurs most frequently at sites of rapid bone growth, the peak incidence occurs earlier in females than males, and males have a higher incidence. Several studies have investigated the relationship between height and birth-weight and OS, but the data are conflicting, therefore, we conducted a meta-analysis of the association of height and birth-weight with OS risk.

We identified 18 studies of OS that included height and/or birth weight. Due to heterogeneity in study findings, the published results of all could not be directly combined. The individual patient data from 6 published studies were obtained from the primary investigators, and 2 unpublished case datasets, resulting in 1067 OS cases with height data, and 434 OS cases with birth-weight data. Since many studies only collected only case information, a control population was created based on the 2000 US National Center for Health Statistics Growth Charts; for each case, 1000 controls were simulated matched on age and gender for height analyses, and matched on gender for birth-weight analyses. Unconditional logistic regression models were used to estimate odds ratios (ORs) for the association between height at diagnosis or birth-weight and risk of OS for each individual study, adjusting for potential confounders. Birth-weight categories were defined based on the lowest 10% and highest 10% of the control distribution. Study specific ORs were combined using random-effects models.

Compared to those of average birth-weight (2665-4045g), patients with low birth-weight (≤2664g) or with high birth-weight (≥4046g) had no significant increased risk of OS [OR 1.09, 95% confidence intervals (CI) 0.78-1.53, and OR 1.41, 95% CI 0.93-2.13, respectively]. By gender, females with high birth-weight had an increased risk of OS (OR 1.49, 95% CI 1.00-2.22); there were no significant associations in males. Patients who were taller (51st to 89th percentile of height for age and gender) than average (≤50th percentile of height) and very tall (≥90th percentile) had a statistically significant increased risk of OS (OR 1.40, 95% CI 1.13-1.73, and OR 2.63, 95% CI 1.98-3.49, respectively; Ptrend &lt;0.0001). A significant trend between increasing height and OS risk was observed in both males and females.

This is the largest aggregated dataset to investigate the associations between height and birth-weight with OS. High birth-weight was associated with increased risk of OS in females, although overall the results did not reach statistical significance. A taller-than-average stature for a person's age and gender significantly increased the risk of OS. This suggests that rapid and/or sustained bone growth during puberty, and possibly in utero, contributes to OS etiology.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 5748.

Prognostic factors and outcomes for osteosarcoma: an international collaboration

We aimed to evaluate the prognostic significance of traditional clinical predictors in osteosarcoma through an international collaboration of 10 teams of investigators (2680 patients) who participated. In multivariate models the mortality risk increased with older age, presence of metastatic disease at diagnosis, development of local recurrence when the patient was first seen, use of amputation instead of limb salvage/wide resection, employment of unusual treatments, use of chemotherapeutic regimens other than anthracycline and platinum and use of methotrexate. It was also influenced by the site of the tumour. The risk of metastasis increased when metastatic disease was present at the time the patient was first seen and also increased with use of amputation or unusual treatment combinations or chemotherapy regimens not including anthracycline and platinum. Local recurrence risk was higher in older patients, in those who had local recurrence when first seen and when no anthracycline and platinum were used in chemotherapy. Results were similar when limited to patients seen after 1990 and treated with surgery plus combination chemotherapy. This large-scale international collaboration identifies strong predictors of major clinical outcomes in osteosarcoma.

CDH11 expression is associated with survival in patients with osteosarcoma

Previous studies have shown that cadherin-11 (CDH11) may be involved in the metastatic process of osteosarcoma. The correlation of the expression levels of CDH11 in osteosarcoma samples with the risk of disease progression and metastasis was examined. Real time qRT-PCR was used to quantify CDH11 expression in a set of newly established osteosarcoma cell lines, 11 primaries and five metastases, compared to the levels in 12 normal osteoblast cell lines established from healthy bone, and also in a set of 10 snap-frozen osteosarcoma samples. In all cases long term clinical follow-up data was available. The CDH11 expression level decreased gradually from the osteoblast to the primary cell lines (p=0.2184) and further to those established from the tumor metastases (p=0.0275). Importantly, the level of CDH11 expression correlated significantly (p=0.01) with patient survival (Kaplan-Meier survival analysis) in both sample sets (p=0.0128 for the cell lines, p=0.0492 for the biopsies). In conclusion, the results indicate that CDH11 may be useful as a prognostic marker of disease progression and survival in osteosarcoma.