Germline Genetic Testing and Survival Outcomes Among Children With Rhabdomyosarcoma: A Report From the Children's Oncology Group

Importance

Determining the impact of germline cancer-predisposition variants (CPVs) on outcomes could inform novel approaches to testing and treating children with rhabdomyosarcoma.

Objective

To assess whether CPVs are associated with outcome among children with rhabdomyosarcoma.

Design, Setting, and Participants

In this cohort study, data were obtained for individuals, aged 0.01-23.23 years, newly diagnosed with rhabdomyosarcoma who were treated across 171 Children's Oncology Group sites from March 15, 1999, to December 8, 2017. Data analysis was performed from June 16, 2021, to May 15, 2023.

Exposure

The presence of a CPV in 24 rhabdomyosarcoma-associated cancer-predisposition genes (CPGs) or an expanded set of 63 autosomal-dominant CPGs.

Main Outcomes and Measures

Overall survival (OS) and event-free survival (EFS) were the main outcomes, using the Kaplan-Meier estimator to assess survival probabilities and the Cox proportional hazards regression model to adjust for clinical covariates. Analyses were stratified by tumor histology and the fusion status of PAX3 or PAX7 to the FOXO1 gene.

Results

In this study of 580 individuals with rhabdomyosarcoma, the median patient age was 5.9 years (range, 0.01-23.23 years), and the male-to-female ratio was 1.5 to 1 (351 [60.5%] male). For patients with CPVs in rhabdomyosarcoma-associated CPGs, EFS was 48.4% compared with 57.8% for patients without a CPV (P = .10), and OS was 53.7% compared with 65.3% for patients without a CPV (P = .06). After adjustment, patients with CPVs had significantly worse OS (adjusted hazard ratio [AHR], 2.49 [95% CI, 1.39-4.45]; P = .002), and the outcomes were not better among patients with embryonal histology (EFS: AHR, 2.25 [95% CI, 1.25-4.06]; P = .007]; OS: AHR, 2.83 [95% CI, 1.47-5.43]; P = .002]). These associations were not due to the development of a second malignant neoplasm, and importantly, patients with fusion-negative rhabdomyosarcoma who harbored a CPV had similarly inferior outcomes as patients with fusion-positive rhabdomyosarcoma without CPVs (EFS: AHR, 1.35 [95% CI, 0.71-2.59]; P = .37; OS: AHR, 1.71 [95% CI, 0.84-3.47]; P = .14). There were no significant differences in outcome by CPV status of the 63 CPG set.

Conclusions and Relevance

This cohort study identified a group of patients with embryonal rhabdomyosarcoma who had a particularly poor outcome. Other important clinical findings included that individuals with TP53 had poor outcomes independent of second malignant neoplasms and that patients with fusion-negative rhabdomyosarcoma who harbored a CPV had outcomes comparable to patients with fusion-positive rhabdomyosarcoma. These findings suggest that germline CPV testing may aid in clinical prognosis and should be considered in prospective risk-based clinical trials.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Genome Sequencing in the Parkinson Disease Clinic

Background and Objectives

Genetic variants affect both Parkinson disease (PD) risk and manifestations. Although genetic information is of potential interest to patients and clinicians, genetic testing is rarely performed during routine PD clinical care. The goal of this study was to examine interest in comprehensive genetic testing among patients with PD and document reactions to possible findings from genome sequencing in 2 academic movement disorder clinics.

Methods

In 203 subjects with PD (age = 63 years, 67% male), genome sequencing was performed and filtered using a custom panel, including 49 genes associated with PD, parkinsonism, or related disorders, as well as a 90-variant PD genetic risk score. Based on the results, 231 patients (age = 67 years, 63% male) were surveyed on interest in genetic testing and responses to vignettes covering (1) familial risk of PD (LRRK2); (2) risk of PD dementia (GBA); (3) PD genetic risk score; and (4) secondary, medically actionable variants (BRCA1).

Results

Genome sequencing revealed a LRRK2 variant in 3% and a GBA risk variant in 10% of our clinical sample. The genetic risk score was normally distributed, identifying 41 subjects with a high risk of PD. Medically actionable findings were discovered in 2 subjects (1%). In our survey, the majority (82%) responded that they would share a LRRK2 variant with relatives. Most registered unchanged or increased interest in testing when confronted with a potential risk for dementia or medically actionable findings, and most (75%) expressed interest in learning their PD genetic risk score.

Discussion

Our results highlight broad interest in comprehensive genetic testing among patients with PD and may facilitate integration of genome sequencing in clinical practice.

Accounting for population structure in genetic studies of cystic fibrosis

CFTR F508del (c.1521_1523delCTT, p.Phe508delPhe) is the most common pathogenic allele underlying cystic fibrosis (CF), and its frequency varies in a geographic cline across Europe. We hypothesized that genetic variation associated with this cline is overrepresented in a large cohort (N > 5,000) of persons with CF who underwent whole-genome sequencing and that this pattern could result in spurious associations between variants correlated with both the F508del genotype and CF-related outcomes. Using principal-component (PC) analyses, we showed that variation in the CFTR region disproportionately contributes to a PC explaining a relatively high proportion of genetic variance. Variation near CFTR was correlated with population structure among persons with CF, and this correlation was driven by a subset of the sample inferred to have European ancestry. We performed genome-wide association studies comparing persons with CF with one versus two copies of the F508del allele; this allowed us to identify genetic variation associated with the F508del allele and to determine that standard PC-adjustment strategies eliminated the significant association signals. Our results suggest that PC adjustment can adequately prevent spurious associations between genetic variants and CF-related traits and are therefore effective tools to control for population structure even when population structure is confounded with disease severity and a common pathogenic variant.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV₁] and its ratio to forced vital capacity [FEV₁/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV₁ and FEV₁/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10^-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Germline Cancer Predisposition Variants in Pediatric Rhabdomyosarcoma: A Report From the Children's Oncology Group

BACKGROUND

Several cancer-susceptibility syndromes are reported to underlie pediatric rhabdomyosarcoma (RMS); however, to our knowledge there have been no systematic efforts to characterize the heterogeneous genetic etiologies of this often-fatal malignancy.

METHODS

We performed exome-sequencing on germline DNA from 615 patients with newly diagnosed RMS consented through the Children's Oncology Group. We compared the prevalence of cancer predisposition variants in 63 autosomal-dominant cancer predisposition genes in these patients with population controls (n = 9963). All statistical tests were 2-sided.

RESULTS

We identified germline cancer predisposition variants in 45 RMS patients (7.3%; all FOXO1 fusion negative) across 15 autosomal dominant genes, which was statistically significantly enriched compared with controls (1.4%, P = 1.3 × 10-22). Specifically, 73.3% of the predisposition variants were found in predisposition syndrome genes previously associated with pediatric RMS risk, such as Li-Fraumeni syndrome (TP53) and neurofibromatosis type I (NF1). Notably, 5 patients had well-described oncogenic missense variants in HRAS (p.G12V and p.G12S) associated with Costello syndrome. Also, genetic etiology differed with histology, as germline variants were more frequent in embryonal vs alveolar RMS patients (10.0% vs 3.0%, P = .02). Although patients with a cancer predisposition variant tended to be younger at diagnosis (P = 9.9 × 10-4), 40.0% of germline variants were identified in those older than 3 years of age, which is in contrast to current genetic testing recommendations based on early age at diagnosis.

CONCLUSIONS

These findings demonstrate that genetic risk of RMS results from germline predisposition variants associated with a wide spectrum of cancer susceptibility syndromes. Germline genetic testing for children with RMS should be informed by RMS subtypes and not be limited to only young patients.

Whole genome sequence analysis of pulmonary function and COPD in 19,996 multi-ethnic participants

Chronic obstructive pulmonary disease (COPD), diagnosed by reduced lung function, is a leading cause of morbidity and mortality. We performed whole genome sequence (WGS) analysis of lung function and COPD in a multi-ethnic sample of 11,497 participants from population- and family-based studies, and 8499 individuals from COPD-enriched studies in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We identify at genome-wide significance 10 known GWAS loci and 22 distinct, previously unreported loci, including two common variant signals from stratified analysis of African Americans. Four novel common variants within the regions of PIAS1, RGN (two variants) and FTO show evidence of replication in the UK Biobank (European ancestry n ~ 320,000), while colocalization analyses leveraging multi-omic data from GTEx and TOPMed identify potential molecular mechanisms underlying four of the 22 novel loci. Our study demonstrates the value of performing WGS analyses and multi-omic follow-up in cohorts of diverse ancestry.

Exome sequencing reveals novel genetic loci influencing obesity-related traits in Hispanic children

OBJECTIVE

To perform whole exome sequencing in 928 Hispanic children and identify variants and genes associated with childhood obesity.

METHODS

Single-nucleotide variants (SNVs) were identified from Illumina whole exome sequencing data using integrated read mapping, variant calling, and an annotation pipeline (Mercury). Association analyses of 74 obesity-related traits and exonic variants were performed using SeqMeta software. Rare autosomal variants were analyzed using gene-based association analyses, and common autosomal variants were analyzed at the SNV level.

RESULTS

(1) Rare exonic variants in 10 genes and 16 common SNVs in 11 genes that were associated with obesity traits in a cohort of Hispanic children were identified, (2) novel rare variants in peroxisome biogenesis factor 1 (PEX1) associated with several obesity traits (weight, weight z score, BMI, BMI z score, waist circumference, fat mass, trunk fat mass) were discovered, and (3) previously reported SNVs associated with childhood obesity were replicated.

CONCLUSIONS

Convergence of whole exome sequencing, a family-based design, and extensive phenotyping discovered novel rare and common variants associated with childhood obesity. Linking PEX1 to obesity phenotypes poses a novel mechanism of peroxisomal biogenesis and metabolism underlying the development of childhood obesity.

A t(5;16) translocation is the likely driver of a syndrome with ambiguous genitalia, facial dysmorphism, intellectual disability, and speech delay

Genetic studies grounded on monogenic paradigms have accelerated both gene discovery and molecular diagnosis. At the same time, complex genomic rearrangements are also appreciated as potent drivers of disease pathology. Here, we report two male siblings with a dysmorphic face, ambiguous genitalia, intellectual disability, and speech delay. Through quad-based whole-exome sequencing and concomitant molecular cytogenetic testing, we identified two copy-number variants (CNVs) in both affected individuals likely arising from a balanced translocation: a 13.5-Mb duplication on Chromosome 16 (16q23.1 → 16qter) and a 7.7-Mb deletion on Chromosome 5 (5p15.31 → 5pter), as well as a hemizygous missense variant in CXorf36 (also known as DIA1R). The 5p terminal deletion has been associated previously with speech delay, whereas craniofacial dysmorphia and genital/urinary anomalies have been reported in patients with a terminal duplication of 16q. However, dosage changes in either genomic region alone could not account for the overall clinical presentation in our family; functional testing of CXorf36 in zebrafish did not induce defects in neurogenesis or the craniofacial skeleton. Notably, literature and database analysis revealed a similar dosage disruption in two siblings with extensive phenotypic overlap with our patients. Taken together, our data suggest that dosage perturbation of genes within the two chromosomal regions likely drives the syndromic manifestations of our patients and highlight how multiple genetic lesions can contribute to complex clinical pathologies.

Whole exome sequencing identifies novel genes for fetal hemoglobin response to hydroxyurea in children with sickle cell anemia

Hydroxyurea has proven efficacy in children and adults with sickle cell anemia (SCA), but with considerable inter-individual variability in the amount of fetal hemoglobin (HbF) produced. Sibling and twin studies indicate that some of that drug response variation is heritable. To test the hypothesis that genetic modifiers influence pharmacological induction of HbF, we investigated phenotype-genotype associations using whole exome sequencing of children with SCA treated prospectively with hydroxyurea to maximum tolerated dose (MTD). We analyzed 171 unrelated patients enrolled in two prospective clinical trials, all treated with dose escalation to MTD. We examined two MTD drug response phenotypes: HbF (final %HbF minus baseline %HbF), and final %HbF. Analyzing individual genetic variants, we identified multiple low frequency and common variants associated with HbF induction by hydroxyurea. A validation cohort of 130 pediatric sickle cell patients treated to MTD with hydroxyurea was genotyped for 13 non-synonymous variants with the strongest association with HbF response to hydroxyurea in the discovery cohort. A coding variant in Spalt-like transcription factor, or SALL2, was associated with higher final HbF in this second independent replication sample and SALL2 represents an outstanding novel candidate gene for further investigation. These findings may help focus future functional studies and provide new insights into the pharmacological HbF upregulation by hydroxyurea in patients with SCA.