Exome variant discrepancies due to reference-genome differences.

Despite release of the GRCh38 human reference genome more than seven years ago, GRCh37 remains more widely used by most research and clinical laboratories. To date, no study has quantified the impact of utilizing different reference assemblies for the identification of variants associated with rare and common diseases from large-scale exome-sequencing data. By calling variants on both the GRCh37 and GRCh38 references, we identified single-nucleotide variants (SNVs) and insertion-deletions (indels) in 1,572 exomes from participants with Mendelian diseases and their family members. We found that a total of 1.5% of SNVs and 2.0% of indels were discordant when different references were used. Notably, 76.6% of the discordant variants were clustered within discrete discordant reference patches (DISCREPs) comprising only 0.9% of loci targeted by exome sequencing. These DISCREPs were enriched for genomic elements including segmental duplications, fix patch sequences, and loci known to contain alternate haplotypes. We identified 206 genes significantly enriched for discordant variants, most of which were in DISCREPs and caused by multi-mapped reads on the reference assembly that lacked the variant call. Among these 206 genes, eight are implicated in known Mendelian diseases and 53 are associated with common phenotypes from genome-wide association studies. In addition, variant interpretations could also be influenced by the reference after lifting-over variant loci to another assembly. Overall, we identified genes and genomic loci affected by reference assembly choice, including genes associated with Mendelian disorders and complex human diseases that require careful evaluation in both research and clinical applications.

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.

TP53 germline pathogenic variant frequency in anaplastic rhabdomyosarcoma: A Children's Oncology Group report.

Rhabdomyosarcoma (RMS) is a well-described cancer in Li-Fraumeni syndrome, resulting from germline TP53 pathogenic variants (PVs). RMS exhibiting anaplasia (anRMS) are associated with a high rate of germline TP53 PVs. This study provides updated estimates of the prevalence of TP53 germline PVs in RMS (3%) and anRMS (11%) from a large cohort (n = 239) enrolled in five Children's Oncology Group (COG) clinical trials. Although the prevalence of germline TP53 PVs in patients with anRMS in this series is much lower than previously reported, this prevalence remains elevated. Germline evaluation for TP53 PVs should be strongly considered in patients with anRMS.

Long read sequencing and expression studies of AHDC1 deletions in Xia-Gibbs syndrome reveal a novel genetic regulatory mechanism.

Xia-Gibbs syndrome (XGS; MIM# 615829) is a rare mendelian disorder characterized by Development Delay (DD), intellectual disability (ID), and hypotonia. Individuals with XGS typically harbor de novo protein-truncating mutations in the AT-Hook DNA binding motif containing 1 (AHDC1) gene, although some missense mutations can also cause XGS. Large de novo heterozygous deletions that encompass the AHDC1 gene have also been ascribed as diagnostic for the disorder, without substantial evidence to support their pathogenicity. We analyzed 19 individuals with large contiguous deletions involving AHDC1, along with other genes. One individual bore the smallest known contiguous AHDC1 deletion (∼350 Kb), encompassing eight other genes within chr1p36.11 (Feline Gardner-Rasheed, IFI6, FAM76A, STX12, PPP1R8, THEMIS2, RPA2, SMPDL3B) and terminating within the first intron of AHDC1. The breakpoint junctions and phase of the deletion were identified using both short and long read sequencing (Oxford Nanopore). Quantification of RNA expression patterns in whole blood revealed that AHDC1 exhibited a mono-allelic expression pattern with no deficiency in overall AHDC1 expression levels, in contrast to the other deleted genes, which exhibited a 50% reduction in mRNA expression. These results suggest that AHDC1 expression in this individual is compensated by a novel regulatory mechanism and advances understanding of mutational and regulatory mechanisms in neurodevelopmental disorders.

Combination of whole exome sequencing and animal modeling identifies TMPRSS9 as a candidate gene for autism spectrum disorder.

Autism spectrum disorders are associated with some degree of developmental regression in up to 30% of all cases. Rarely, however, is the regression so extreme that a developmentally advanced young child would lose almost all ability to communicate and interact with her surroundings. We applied trio whole exome sequencing to a young woman who experienced extreme developmental regression starting at 2.5 years of age and identified compound heterozygous nonsense mutations in TMPRSS9, which encodes for polyserase-1, a transmembrane serine protease of poorly understood physiological function. Using semiquantitative polymerase chain reaction, we showed that Tmprss9 is expressed in various mouse tissues, including the brain. To study the consequences of TMPRSS9 loss of function on the mammalian brain, we generated a knockout mouse model. Through a battery of behavioral assays, we found that Tmprss9-/- mice showed decreased social interest and social recognition. We observed a borderline recognition memory deficit by novel object recognition in aged Tmprss9-/- female mice, but not in aged Tmprss9-/- male mice or younger adult Tmprss9-/- mice in both sexes. This study provides evidence to suggest that loss of function variants in TMPRSS9 are related to an autism spectrum disorder. However, the identification of more individuals with similar phenotypes and TMPRSS9 loss of function variants is required to establish a robust gene-disease relationship.

A Genocentric Approach to Discovery of Mendelian Disorders.

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

Phenotypic and protein localization heterogeneity associated with AHDC1 pathogenic protein-truncating alleles in Xia-Gibbs syndrome.

Xia-Gibbs syndrome (XGS) is a rare Mendelian disease typically caused by de novo stop-gain or frameshift mutations in the AT-hook DNA binding motif containing 1 (AHDC1) gene. Patients usually present in early infancy with hypotonia and developmental delay and later exhibit intellectual disability (ID). The overall presentation is variable, however, and the emerging clinical picture is still evolving. A detailed phenotypic analysis of 34 XGS individuals revealed five core phenotypes (delayed motor milestones, speech delay, low muscle tone, ID, and hypotonia) in more than 80% of individuals and an additional 12 features that occurred more variably. Seizures and scoliosis were more frequently associated with truncations that arise before the midpoint of the protein although the occurrence of most features could not be predicted by the mutation position. Transient expression of wild type and different patient truncated AHDC1 protein forms in human cell lines revealed abnormal patterns of nuclear localization including a diffuse distribution of a short truncated form and nucleolar aggregation in mid-protein truncated forms. Overall, both the occurrence of variable phenotypes and the different distribution of the expressed protein reflect the heterogeneity of this syndrome.

Neptune: an environment for the delivery of genomic medicine.

PURPOSE: Genomic medicine holds great promise for improving health care, but integrating searchable and actionable genetic data into electronic health records (EHRs) remains a challenge. Here we describe Neptune, a system for managing the interaction between a clinical laboratory and an EHR system during the clinical reporting process. METHODS: We developed Neptune and applied it to two clinical sequencing projects that required report customization, variant reanalysis, and EHR integration. RESULTS: Neptune has been applied for the generation and delivery of over 15,000 clinical genomic reports. This work spans two clinical tests based on targeted gene panels that contain 68 and 153 genes respectively. These projects demanded customizable clinical reports that contained a variety of genetic data types including single-nucleotide variants (SNVs), copy-number variants (CNVs), pharmacogenomics, and polygenic risk scores. Two variant reanalysis activities were also supported, highlighting this important workflow. CONCLUSION: Methods are needed for delivering structured genetic data to EHRs. This need extends beyond developing data formats to providing infrastructure that manages the reporting process itself. Neptune was successfully applied on two high-throughput clinical sequencing projects to build and deliver clinical reports to EHR systems. The software is open source and available at https://gitlab.com/bcm-hgsc/neptune .

Genetic testing in ambulatory cardiology clinics reveals high rate of findings with clinical management implications.

PURPOSE: Cardiovascular disease (CVD) is the leading cause of death in adults in the United States, yet the benefits of genetic testing are not universally accepted. METHODS: We developed the "HeartCare" panel of genes associated with CVD, evaluating high-penetrance Mendelian conditions, coronary artery disease (CAD) polygenic risk, LPA gene polymorphisms, and specific pharmacogenetic (PGx) variants. We enrolled 709 individuals from cardiology clinics at Baylor College of Medicine, and samples were analyzed in a CAP/CLIA-certified laboratory. Results were returned to the ordering physician and uploaded to the electronic medical record. RESULTS: Notably, 32% of patients had a genetic finding with clinical management implications, even after excluding PGx results, including 9% who were molecularly diagnosed with a Mendelian condition. Among surveyed physicians, 84% reported medical management changes based on these results, including specialist referrals, cardiac tests, and medication changes. LPA polymorphisms and high polygenic risk of CAD were found in 20% and 9% of patients, respectively, leading to diet, lifestyle, and other changes. Warfarin and simvastatin pharmacogenetic variants were present in roughly half of the cohort. CONCLUSION: Our results support the use of genetic information in routine cardiovascular health management and provide a roadmap for accompanying research.

Exome sequencing in children with clinically suspected maturity-onset diabetes of the young.

OBJECTIVE: Commercial gene panels identify pathogenic variants in as low as 27% of patients suspected to have MODY, suggesting the role of yet unidentified pathogenic variants. We sought to identify novel gene variants associated with MODY. RESEARCH DESIGN AND METHODS: We recruited 10 children with a clinical suspicion of MODY but non-diagnostic commercial MODY gene panels. We performed exome sequencing (ES) in them and their parents. RESULTS: Mean age at diabetes diagnosis was 10 (± 3.8) years. Six were females; 4 were non-Hispanic white, 5 Hispanic, and 1 Asian. Our variant prioritization analysis identified a pathogenic, de novo variant in INS (c.94G > A, p.Gly32Ser), confirmed by Sanger sequencing, in a proband who was previously diagnosed with "autoantibody-negative type 1 diabetes (T1D)" at 3 y/o. This rare variant, absent in the general population (gnomAD database), has been reported previously in neonatal diabetes. We also identified a frameshift deletion (c.2650delC, p.Gln884AsnfsTer57) in RFX6 in a child with a previous diagnosis of "autoantibody-negative T1D" at 12 y/o. The variant was inherited from the mother, who was diagnosed with "thin type 2 diabetes" at 25 y/o. Heterozygous protein-truncating variants in RFX6 gene have been recently reported in individuals with MODY. CONCLUSIONS: We diagnosed two patients with MODY using ES in children initially classified as "T1D". One has a likely pathogenic novel gene variant not previously associated with MODY. We demonstrate the clinical utility of ES in patients with clinical suspicion of MODY.

Identification of Polycystic Kidney Disease 1 Like 1 Gene Variants in Children With Biliary Atresia Splenic Malformation Syndrome.

Biliary atresia (BA) is the most common cause of end-stage liver disease in children and the primary indication for pediatric liver transplantation, yet underlying etiologies remain unknown. Approximately 10% of infants affected by BA exhibit various laterality defects (heterotaxy) including splenic abnormalities and complex cardiac malformations-a distinctive subgroup commonly referred to as the biliary atresia splenic malformation (BASM) syndrome. We hypothesized that genetic factors linking laterality features with the etiopathogenesis of BA in BASM patients could be identified through whole-exome sequencing (WES) of an affected cohort. DNA specimens from 67 BASM subjects, including 58 patient-parent trios, from the National Institute of Diabetes and Digestive and Kidney Diseases-supported Childhood Liver Disease Research Network (ChiLDReN) underwent WES. Candidate gene variants derived from a prespecified set of 2,016 genes associated with ciliary dysgenesis and/or dysfunction or cholestasis were prioritized according to pathogenicity, population frequency, and mode of inheritance. Five BASM subjects harbored rare and potentially deleterious biallelic variants in polycystic kidney disease 1 like 1 (PKD1L1), a gene associated with ciliary calcium signaling and embryonic laterality determination in fish, mice, and humans. Heterozygous PKD1L1 variants were found in 3 additional subjects. Immunohistochemical analysis of liver from the one BASM subject available revealed decreased PKD1L1 expression in bile duct epithelium when compared to normal livers and livers affected by other noncholestatic diseases. Conclusion: WES identified biallelic and heterozygous PKD1L1 variants of interest in 8 BASM subjects from the ChiLDReN data set; the dual roles for PKD1L1 in laterality determination and ciliary function suggest that PKD1L1 is a biologically plausible, cholangiocyte-expressed candidate gene for the BASM syndrome.

Copy-Number Variation Contributes to the Mutational Load of Bardet-Biedl Syndrome.

Bardet-Biedl syndrome (BBS) is a defining ciliopathy, notable for extensive allelic and genetic heterogeneity, almost all of which has been identified through sequencing. Recent data have suggested that copy-number variants (CNVs) also contribute to BBS. We used a custom oligonucleotide array comparative genomic hybridization (aCGH) covering 20 genes that encode intraflagellar transport (IFT) components and 74 ciliopathy loci to screen 92 unrelated individuals with BBS, irrespective of their known mutational burden. We identified 17 individuals with exon-disruptive CNVs (18.5%), including 13 different deletions in eight BBS genes (BBS1, BBS2, ARL6/BBS3, BBS4, BBS5, BBS7, BBS9, and NPHP1) and a deletion and a duplication in other ciliopathy-associated genes (ALMS1 and NPHP4, respectively). By contrast, we found a single heterozygous exon-disruptive event in a BBS-associated gene (BBS9) in 229 control subjects. Superimposing these data with resequencing revealed CNVs to (1) be sufficient to cause disease, (2) Mendelize heterozygous deleterious alleles, and (3) contribute oligogenic alleles by combining point mutations and exonic CNVs in multiple genes. Finally, we report a deletion and a splice site mutation in IFT74, inherited under a recessive paradigm, defining a candidate BBS locus. Our data suggest that CNVs contribute pathogenic alleles to a substantial fraction of BBS-affected individuals and highlight how either deletions or point mutations in discrete splice isoforms can induce hypomorphic mutations in genes otherwise intolerant to deleterious variation. Our data also suggest that CNV analyses and resequencing studies unbiased for previous mutational burden is necessary to delineate the complexity of disease architecture.

Patterns and rates of exonic de novo mutations in autism spectrum disorders.

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.

Recurrent CNVs and SNVs at the NPHP1 locus contribute pathogenic alleles to Bardet-Biedl syndrome.

Homozygosity for a recurrent 290 kb deletion of NPHP1 is the most frequent cause of isolated nephronophthisis (NPHP) in humans. A deletion of the same genomic interval has also been detected in individuals with Joubert syndrome (JBTS), and in the mouse, Nphp1 interacts genetically with Ahi1, a known JBTS locus. Given these observations, we investigated the contribution of NPHP1 in Bardet-Biedl syndrome (BBS), a ciliopathy of intermediate severity. By using a combination of array-comparative genomic hybridization, TaqMan copy number assays, and sequencing, we studied 200 families affected by BBS. We report a homozygous NPHP1 deletion CNV in a family with classical BBS that is transmitted with autosomal-recessive inheritance. Further, we identified heterozygous NPHP1 deletions in two more unrelated persons with BBS who bear primary mutations at another BBS locus. In parallel, we identified five families harboring an SNV in NPHP1 resulting in a conserved missense change, c.14G>T (p.Arg5Leu), that is enriched in our Hispanic pedigrees; in each case, affected individuals carried additional bona fide pathogenic alleles in another BBS gene. In vivo functional modeling in zebrafish embryos demonstrated that c.14G>T is a loss-of-function variant, and suppression of nphp1 in concert with each of the primary BBS loci found in our NPHP1-positive pedigrees exacerbated the severity of the phenotype. These results suggest that NPHP1 mutations are probably rare primary causes of BBS that contribute to the mutational burden of the disorder.

Meta-analysis of 49†549 individuals imputed with the 1000 Genomes Project reveals an exonic damaging variant in ANGPTL4 determining fasting TG levels.

BACKGROUND: So far, more than 170 loci have been associated with circulating lipid levels through genome-wide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels. METHODS: We used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ∼60 000 individuals in the discovery stage and ∼90 000 samples in the replication stage. RESULTS: Our study resulted in the identification of five new associations with circulating lipid levels at four loci. All four loci are within genes that can be linked biologically to lipid metabolism. One of the variants, rs116843064, is a damaging missense variant within the ANGPTL4 gene. CONCLUSIONS: This study illustrates that GWAS with high-scale imputation may still help us unravel the biological mechanism behind circulating lipid levels.

Analysis of rare, exonic variation amongst subjects with autism spectrum disorders and population controls.

We report on results from whole-exome sequencing (WES) of 1,039 subjects diagnosed with autism spectrum disorders (ASD) and 870 controls selected from the NIMH repository to be of similar ancestry to cases. The WES data came from two centers using different methods to produce sequence and to call variants from it. Therefore, an initial goal was to ensure the distribution of rare variation was similar for data from different centers. This proved straightforward by filtering called variants by fraction of missing data, read depth, and balance of alternative to reference reads. Results were evaluated using seven samples sequenced at both centers and by results from the association study. Next we addressed how the data and/or results from the centers should be combined. Gene-based analyses of association was an obvious choice, but should statistics for association be combined across centers (meta-analysis) or should data be combined and then analyzed (mega-analysis)? Because of the nature of many gene-based tests, we showed by theory and simulations that mega-analysis has better power than meta-analysis. Finally, before analyzing the data for association, we explored the impact of population structure on rare variant analysis in these data. Like other recent studies, we found evidence that population structure can confound case-control studies by the clustering of rare variants in ancestry space; yet, unlike some recent studies, for these data we found that principal component-based analyses were sufficient to control for ancestry and produce test statistics with appropriate distributions. After using a variety of gene-based tests and both meta- and mega-analysis, we found no new risk genes for ASD in this sample. Our results suggest that standard gene-based tests will require much larger samples of cases and controls before being effective for gene discovery, even for a disorder like ASD.

Assessing structural variation in a personal genome-towards a human reference diploid genome.

BACKGROUND: Characterizing large genomic variants is essential to expanding the research and clinical applications of genome sequencing. While multiple data types and methods are available to detect these structural variants (SVs), they remain less characterized than smaller variants because of SV diversity, complexity, and size. These challenges are exacerbated by the experimental and computational demands of SV analysis. Here, we characterize the SV content of a personal genome with Parliament, a publicly available consensus SV-calling infrastructure that merges multiple data types and SV detection methods. RESULTS: We demonstrate Parliament's efficacy via integrated analyses of data from whole-genome array comparative genomic hybridization, short-read next-generation sequencing, long-read (Pacific BioSciences RSII), long-insert (Illumina Nextera), and whole-genome architecture (BioNano Irys) data from the personal genome of a single subject (HS1011). From this genome, Parliament identified 31,007 genomic loci between 100 bp and 1 Mbp that are inconsistent with the hg19 reference assembly. Of these loci, 9,777 are supported as putative SVs by hybrid local assembly, long-read PacBio data, or multi-source heuristics. These SVs span 59 Mbp of the reference genome (1.8%) and include 3,801 events identified only with long-read data. The HS1011 data and complete Parliament infrastructure, including a BAM-to-SV workflow, are available on the cloud-based service DNAnexus. CONCLUSIONS: HS1011 SV analysis reveals the limits and advantages of multiple sequencing technologies, specifically the impact of long-read SV discovery. With the full Parliament infrastructure, the HS1011 data constitute a public resource for novel SV discovery, software calibration, and personal genome structural variation analysis.

Somatic mutations affect key pathways in lung adenocarcinoma.

Determining the genetic basis of cancer requires comprehensive analyses of large collections of histopathologically well-classified primary tumours. Here we report the results of a collaborative study to discover somatic mutations in 188 human lung adenocarcinomas. DNA sequencing of 623 genes with known or potential relationships to cancer revealed more than 1,000 somatic mutations across the samples. Our analysis identified 26 genes that are mutated at significantly high frequencies and thus are probably involved in carcinogenesis. The frequently mutated genes include tyrosine kinases, among them the EGFR homologue ERBB4; multiple ephrin receptor genes, notably EPHA3; vascular endothelial growth factor receptor KDR; and NTRK genes. These data provide evidence of somatic mutations in primary lung adenocarcinoma for several tumour suppressor genes involved in other cancers--including NF1, APC, RB1 and ATM--and for sequence changes in PTPRD as well as the frequently deleted gene LRP1B. The observed mutational profiles correlate with clinical features, smoking status and DNA repair defects. These results are reinforced by data integration including single nucleotide polymorphism array and gene expression array. Our findings shed further light on several important signalling pathways involved in lung adenocarcinoma, and suggest new molecular targets for treatment.

Human genetic variation determines 24-hour rhythmic gene expression and disease risk.

24-hour biological rhythms are essential to maintain physiological homeostasis. Disruption of these rhythms increases the risks of multiple diseases. The biological rhythms are known to have a genetic basis formed by core clock genes, but how individual genetic variation shapes the oscillating transcriptome and contributes to human chronophysiology and disease risk is largely unknown. Here, we mapped interactions between temporal gene expression and genotype to identify quantitative trait loci (QTLs) contributing to rhythmic gene expression. These newly identified QTLs were termed as rhythmic QTLs (rhyQTLs), which determine previously unappreciated rhythmic genes in human subpopulations with specific genotypes. Functionally, rhyQTLs and their associated rhythmic genes contribute extensively to essential chronophysiological processes, including bile acid and lipid metabolism. The identification of rhyQTLs sheds light on the genetic mechanisms of gene rhythmicity, offers mechanistic insights into variations in human disease risk, and enables precision chronotherapeutic approaches for patients.

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.