Federated analysis of autosomal recessive coding variants in 29,745 developmental disorder patients from diverse populations.

Autosomal recessive coding variants are well-known causes of rare disorders. We quantified the contribution of these variants to developmental disorders in a large, ancestrally diverse cohort comprising 29,745 trios, of whom 20.4% had genetically inferred non-European ancestries. The estimated fraction of patients attributable to exome-wide autosomal recessive coding variants ranged from ~2-19% across genetically inferred ancestry groups and was significantly correlated with average autozygosity. Established autosomal recessive developmental disorder-associated (ARDD) genes explained 84.0% of the total autosomal recessive coding burden, and 34.4% of the burden in these established genes was explained by variants not already reported as pathogenic in ClinVar. Statistical analyses identified two novel ARDD genes: KBTBD2 and ZDHHC16. This study expands our understanding of the genetic architecture of developmental disorders across diverse genetically inferred ancestry groups and suggests that improving strategies for interpreting missense variants in known ARDD genes may help diagnose more patients than discovering the remaining genes.

Genotype-first analysis in an unselected health system-based population reveals variable phenotypic severity of COL4A5 variants.

Introduction: Our knowledge of X-linked Alport Syndrome [AS] comes mostly from selected cohorts with more severe disease. Methods: We examined the phenotypic spectrum of X-linked AS in males and females with a genotype-based approach using data from the Geisinger MyCode DiscovEHR study, an unselected health system-based cohort with exome sequencing and electronic health records. Patients with COL4A5 variants reported as pathogenic (P) or likely pathogenic (LP) in ClinVar, or protein-truncating variants (PTVs), were each matched with up to 5 controls without COL4A3/4/5 variants by sociodemographics, diabetes diagnosis, and year of first outpatient encounter. AS-related phenotypes included dipstick hematuria, bilateral sensorineural hearing loss (BSHL), proteinuria, decreased eGFR, and ESKD. Results: Out of 170,856 patients, there were 29 hemizygous males (mean age 52.0 y [SD 20.0]) and 55 heterozygous females (mean age 59.3 y [SD 18.8]) with a COL4A5 P/LP variant, including 48 with the hypomorphic variant p.Gly624Asp. Overall, penetrance (having any AS phenotypic feature) was highest for non-p.Gly624Asp P/LP variants (males: 94%, females: 85%), intermediate for p.Gly624Asp (males: 77%, females: 69%), compared to controls (males: 32%; females: 50%). The proportion with ESKD was highest for males with P/LP variants (44%), intermediate for males with p.Gly624Asp (15%) and females with P/LP variants (10%), compared to controls (males: 3%, females 2%). Only 47% of individuals with COL4A5 had completed albuminuria screening, and a minority were taking renin-angiotensin aldosterone system (RAAS) inhibitors. Only 38% of males and 16% of females had a known diagnosis of Alport syndrome or thin basement membrane disease. Conclusion: In an unselected cohort, we show increased risks of AS-related phenotypes in men and women compared to matched controls, while showing a wider spectrum of severity than has been described previously and variability by genotype. Future studies are needed to determine whether early genetic diagnosis can improve outcomes in Alport Syndrome.

Systematic analysis of variants escaping nonsense-mediated decay uncovers candidate Mendelian diseases.

Protein-truncating variants (PTVs) near the 3' end of genes may escape nonsense-mediated decay (NMD). PTVs in the NMD-escape region (PTVescs) can cause Mendelian disease but are difficult to interpret given their varying impact on protein function. Previously, PTVesc burden was assessed in an epilepsy cohort, but no large-scale analysis has systematically evaluated these variants in rare disease. We performed a retrospective analysis of 29,031 neurodevelopmental disorder (NDD) parent-offspring trios referred for clinical exome sequencing to identify PTVesc de novo mutations (DNMs). We identified 1,376 PTVesc DNMs and 133 genes that were significantly enriched (binomial p < 0.001). The PTVesc-enriched genes included those with PTVescs previously described to cause dominant Mendelian disease (e.g., SEMA6B, PPM1D, and DAGLA). We annotated ClinVar variants for PTVescs and identified 948 genes with at least one high-confidence pathogenic variant. Twenty-two known Mendelian PTVesc-enriched genes had no prior evidence of PTVesc-associated disease. We found 22 additional PTVesc-enriched genes that are not well established to be associated with Mendelian disease, several of which showed phenotypic similarity between individuals harboring PTVesc variants in the same gene. Four individuals with PTVesc mutations in RAB1A had similar phenotypes including NDD and spasticity. PTVesc mutations in IRF2BP1 were found in two individuals who each had severe immunodeficiency manifesting in NDD. Three individuals with PTVesc mutations in LDB1 all had NDD and multiple congenital anomalies. Using a large-scale, systematic analysis of DNMs, we extend the mutation spectrum for known Mendelian disease-associated genes and identify potentially novel disease-associated genes.

The landscape of reported VUS in multi-gene panel and genomic testing: Time for a change.

PURPOSE: Variants of uncertain significance (VUS) are a common result of diagnostic genetic testing and can be difficult to manage with potential misinterpretation and downstream costs, including time investment by clinicians. We investigated the rate of VUS reported on diagnostic testing via multi-gene panels (MGPs) and exome and genome sequencing (ES/GS) to measure the magnitude of uncertain results and explore ways to reduce their potentially detrimental impact. METHODS: Rates of inconclusive results due to VUS were collected from over 1.5 million sequencing test results from 19 clinical laboratories in North America from 2020 to 2021. RESULTS: We found a lower rate of inconclusive test results due to VUSs from ES/GS (22.5%) compared with MGPs (32.6%; P < .0001). For MGPs, the rate of inconclusive results correlated with panel size. The use of trios reduced inconclusive rates (18.9% vs 27.6%; P < .0001), whereas the use of GS compared with ES had no impact (22.2% vs 22.6%; P = ns). CONCLUSION: The high rate of VUS observed in diagnostic MGP testing warrants examining current variant reporting practices. We propose several approaches to reduce reported VUS rates, while directing clinician resources toward important VUS follow-up.

Rare pathogenic variants in WNK3 cause X-linked intellectual disability.

PURPOSE: WNK3 kinase (PRKWNK3) has been implicated in the development and function of the brain via its regulation of the cation-chloride cotransporters, but the role of WNK3 in human development is unknown. METHOD: We ascertained exome or genome sequences of individuals with rare familial or sporadic forms of intellectual disability (ID). RESULTS: We identified a total of 6 different maternally-inherited, hemizygous, 3 loss-of-function or 3 pathogenic missense variants (p.Pro204Arg, p.Leu300Ser, p.Glu607Val) in WNK3 in 14 male individuals from 6 unrelated families. Affected individuals had ID with variable presence of epilepsy and structural brain defects. WNK3 variants cosegregated with the disease in 3 different families with multiple affected individuals. This included 1 large family previously diagnosed with X-linked Prieto syndrome. WNK3 pathogenic missense variants localize to the catalytic domain and impede the inhibitory phosphorylation of the neuronal-specific chloride cotransporter KCC2 at threonine 1007, a site critically regulated during the development of synaptic inhibition. CONCLUSION: Pathogenic WNK3 variants cause a rare form of human X-linked ID with variable epilepsy and structural brain abnormalities and implicate impaired phospho-regulation of KCC2 as a pathogenic mechanism.

Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms.

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

Molecular Diagnostic Yield of Exome Sequencing in Patients With Cerebral Palsy.

Importance: Cerebral palsy is a common neurodevelopmental disorder affecting movement and posture that often co-occurs with other neurodevelopmental disorders. Individual cases of cerebral palsy are often attributed to birth asphyxia; however, recent studies indicate that asphyxia accounts for less than 10% of cerebral palsy cases. Objective: To determine the molecular diagnostic yield of exome sequencing (prevalence of pathogenic and likely pathogenic variants) in individuals with cerebral palsy. Design, Setting, and Participants: A retrospective cohort study of patients with cerebral palsy that included a clinical laboratory referral cohort with data accrued between 2012 and 2018 and a health care-based cohort with data accrued between 2007 and 2017. Exposures: Exome sequencing with copy number variant detection. Main Outcomes and Measures: The primary outcome was the molecular diagnostic yield of exome sequencing. Results: Among 1345 patients from the clinical laboratory referral cohort, the median age was 8.8 years (interquartile range, 4.4-14.7 years; range, 0.1-66 years) and 601 (45%) were female. Among 181 patients in the health care-based cohort, the median age was 41.9 years (interquartile range, 28.0-59.6 years; range, 4.8-89 years) and 96 (53%) were female. The molecular diagnostic yield of exome sequencing was 32.7% (95% CI, 30.2%-35.2%) in the clinical laboratory referral cohort and 10.5% (95% CI, 6.0%-15.0%) in the health care-based cohort. The molecular diagnostic yield ranged from 11.2% (95% CI, 6.4%-16.2%) for patients without intellectual disability, epilepsy, or autism spectrum disorder to 32.9% (95% CI, 25.7%-40.1%) for patients with all 3 comorbidities. Pathogenic and likely pathogenic variants were identified in 229 genes (29.5% of 1526 patients); 86 genes were mutated in 2 or more patients (20.1% of 1526 patients) and 10 genes with mutations were independently identified in both cohorts (2.9% of 1526 patients). Conclusions and Relevance: Among 2 cohorts of patients with cerebral palsy who underwent exome sequencing, the prevalence of pathogenic and likely pathogenic variants was 32.7% in a cohort that predominantly consisted of pediatric patients and 10.5% in a cohort that predominantly consisted of adult patients. Further research is needed to understand the clinical implications of these findings.

A dyadic approach to the delineation of diagnostic entities in clinical genomics.

The delineation of disease entities is complex, yet recent advances in the molecular characterization of diseases provide opportunities to designate diseases in a biologically valid manner. Here, we have formalized an approach to the delineation of Mendelian genetic disorders that encompasses two distinct but inter-related concepts: (1) the gene that is mutated and (2) the phenotypic descriptor, preferably a recognizably distinct phenotype. We assert that only by a combinatorial or dyadic approach taking both of these attributes into account can a unitary, distinct genetic disorder be designated. We propose that all Mendelian disorders should be designated as "GENE-related phenotype descriptor" (e.g., "CFTR-related cystic fibrosis"). This approach to delineating and naming disorders reconciles the complexity of gene-to-phenotype relationships in a simple and clear manner yet communicates the complexity and nuance of these relationships.

Uniparental disomy in a population of 32,067 clinical exome trios.

PURPOSE: Data on the clinical prevalence and spectrum of uniparental disomy (UPD) remain limited. Trio exome sequencing (ES) presents a comprehensive method for detection of UPD alongside sequence and copy-number variant analysis. METHODS: We analyzed 32,067 ES trios referred for diagnostic testing to create a profile of UPD events and their disease associations. ES single-nucleotide polymorphism (SNP) and copy-number data were used to identify both whole-chromosome and segmental UPD and to categorize whole-chromosome results as isodisomy, heterodisomy, or mixed. RESULTS: Ninety-nine whole-chromosome and 13 segmental UPD events were identified. Of these, 29 were associated with an imprinting disorder, and 16 were associated with a positive test result through homozygous sequence variants. Isodisomy was more commonly observed in large chromosomes along with a higher rate of homozygous pathogenic variants, while heterodisomy was more frequent in chromosomes associated with imprinting or trisomy mosaicism (14, 15, 16, 20, 22). CONCLUSION: Whole-chromosome UPD was observed in 0.31% of cases, resulting in a diagnostic finding in 0.14%. Only three UPD-positive cases had a diagnostic finding unrelated to the UPD. Thirteen UPD events were identified in cases with prior normal SNP chromosomal microarray results, demonstrating the additional diagnostic value of UPD detection by trio ES.

Histone H3.3 beyond cancer: Germline mutations in Histone 3 Family 3A and 3B cause a previously unidentified neurodegenerative disorder in 46 patients.

Although somatic mutations in Histone 3.3 (H3.3) are well-studied drivers of oncogenesis, the role of germline mutations remains unreported. We analyze 46 patients bearing de novo germline mutations in histone 3 family 3A (H3F3A) or H3F3B with progressive neurologic dysfunction and congenital anomalies without malignancies. Molecular modeling of all 37 variants demonstrated clear disruptions in interactions with DNA, other histones, and histone chaperone proteins. Patient histone posttranslational modifications (PTMs) analysis revealed notably aberrant local PTM patterns distinct from the somatic lysine mutations that cause global PTM dysregulation. RNA sequencing on patient cells demonstrated up-regulated gene expression related to mitosis and cell division, and cellular assays confirmed an increased proliferative capacity. A zebrafish model showed craniofacial anomalies and a defect in Foxd3-derived glia. These data suggest that the mechanism of germline mutations are distinct from cancer-associated somatic histone mutations but may converge on control of cell proliferation.

Evidence for 28 genetic disorders discovered by combining healthcare and research data.

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

Mutations disrupting neuritogenesis genes confer risk for cerebral palsy.

In addition to commonly associated environmental factors, genomic factors may cause cerebral palsy. We performed whole-exome sequencing of 250 parent-offspring trios, and observed enrichment of damaging de novo mutations in cerebral palsy cases. Eight genes had multiple damaging de novo mutations; of these, two (TUBA1A and CTNNB1) met genome-wide significance. We identified two novel monogenic etiologies, FBXO31 and RHOB, and showed that the RHOB mutation enhances active-state Rho effector binding while the FBXO31 mutation diminishes cyclin D levels. Candidate cerebral palsy risk genes overlapped with neurodevelopmental disorder genes. Network analyses identified enrichment of Rho GTPase, extracellular matrix, focal adhesion and cytoskeleton pathways. Cerebral palsy risk genes in enriched pathways were shown to regulate neuromotor function in a Drosophila reverse genetics screen. We estimate that 14% of cases could be attributed to an excess of damaging de novo or recessive variants. These findings provide evidence for genetically mediated dysregulation of early neuronal connectivity in cerebral palsy.

Mobile element insertion detection in 89,874 clinical exomes.

PURPOSE: Exome sequencing (ES) is increasingly used for the diagnosis of rare genetic disease. However, some pathogenic sequence variants within the exome go undetected due to the technical difficulty of identifying them. Mobile element insertions (MEIs) are a known cause of genetic disease in humans but have been historically difficult to detect via ES and similar targeted sequencing methods. METHODS: We developed and applied a novel MEI detection method prospectively to samples received for clinical ES beginning in November 2017. Positive MEI findings were confirmed by an orthogonal method and reported back to the ordering provider. In this study, we examined 89,874 samples from 38,871 cases. RESULTS: Diagnostic MEIs were present in 0.03% (95% binomial test confidence interval: 0.02-0.06%) of all cases and account for 0.15% (95% binomial test confidence interval: 0.08-0.25%) of cases with a molecular diagnosis. One diagnostic MEI was a novel founder event. Most patients with pathogenic MEIs had prior genetic testing, three of whom had previous negative DNA sequencing analysis of the diagnostic gene. CONCLUSION: MEI detection from ES is a valuable diagnostic tool, reveals molecular findings that may be undetected by other sequencing assays, and increases diagnostic yield by 0.15%.

The tale of two genes: from next-generation sequencing to phenotype.

An 18-yr-old man with a history of intellectual disability, craniofacial dysmorphism, seizure disorder, and obesity was identified to carry a de novo, pathogenic variant in ASXL1 (c.4198G>T; p.E1400X) associated with the diagnosis of Bohring-Opitz syndrome based on exome sequencing. In addition, he was identified to carry a maternally inherited and likely pathogenic variant in MC4R (c.817C>T; p.Q273X) associated with monogenic obesity. Dual genetic diagnosis occurs in 4%-6% of patients and results in unique clinical phenotypes that are a function of tissue-specific gene expression, involved pathways, clinical expressivity, and penetrance. This case highlights the utility of next-generation sequencing in patients with an unusual combination of clinical presentations for several pillars of precision medicine including (1) diagnosis, (2) prognosis and outcome, (3) management and therapy, and (4) utilization of resources.

Sex-Based Analysis of De Novo Variants in Neurodevelopmental Disorders.

While genes with an excess of de novo mutations (DNMs) have been identified in children with neurodevelopmental disorders (NDDs), few studies focus on DNM patterns where the sex of affected children is examined separately. We considered ∼8,825 sequenced parent-child trios (n ∼26,475 individuals) and identify 54 genes with a DNM enrichment in males (n = 18), females (n = 17), or overlapping in both the male and female subsets (n = 19). A replication cohort of 18,778 sequenced parent-child trios (n = 56,334 individuals) confirms 25 genes (n = 3 in males, n = 7 in females, n = 15 in both male and female subsets). As expected, we observe significant enrichment on the X chromosome for females but also find autosomal genes with potential sex bias (females, CDK13, ITPR1; males, CHD8, MBD5, SYNGAP1); 6.5% of females harbor a DNM in a female-enriched gene, whereas 2.7% of males have a DNM in a male-enriched gene. Sex-biased genes are enriched in transcriptional processes and chromatin binding, primarily reside in the nucleus of cells, and have brain expression. By downsampling, we find that DNM gene discovery is greatest when studying affected females. Finally, directly comparing de novo allele counts in NDD-affected males and females identifies one replicated genome-wide significant gene (DDX3X) with locus-specific enrichment in females. Our sex-based DNM enrichment analysis identifies candidate NDD genes differentially affecting males and females and indicates that the study of females with NDDs leads to greater gene discovery consistent with the female-protective effect.

Age-adjusted association of homologous recombination genes with ovarian cancer using clinical exomes as controls.

BACKGROUND: Genes in the homologous recombination pathway have shown varying results in the literature regarding ovarian cancer (OC) association. Recent case-control studies have used allele counts alone to quantify genetic associations with cancer. METHODS: A retrospective case-control study was performed on 6,182 women with OC referred for hereditary cancer multi-gene panel testing (cases) and 4,690 mothers from trios who were referred for whole-exome sequencing (controls). We present age-adjusted odds ratios (ORAdj) to determine association of OC with pathogenic variants (PVs) in homologous recombination genes. RESULTS: Significant associations with OC were observed in BRCA1, BRCA2, RAD51C and RAD51D. Other homologous recombination genes, BARD1, NBN, and PALB2, were not significantly associated with OC. ATM and CHEK2 were only significantly associated with OC by crude odds ratio (ORCrude) or by ORAdj, respectively. However, there was no significant difference between ORCrude and ORAdj for these two genes. The significant association of PVs in BRIP1 by ORCrude (2.05, CI = 1.11 to 3.94, P = 0.03) was not observed by ORAdj (0.87, CI = 0.41 to 1.93, P = 0.73). Interestingly, the confidence intervals of the two effect sizes were significantly different (P = 0.04). CONCLUSION: The lack of association of PVs in BRIP1 with OC by ORAdj is inconsistent with some previous literature and current management recommendations, highlighted by the significantly older age of OC onset for BRIP1 PV carriers compared to non-carriers. By reporting ORAdj, this study presents associations that reflect more informed genetic contributions to OC when compared to traditional count-based methods.

Missense Variants in the Histone Acetyltransferase Complex Component Gene TRRAP Cause Autism and Syndromic Intellectual Disability.

Acetylation of the lysine residues in histones and other DNA-binding proteins plays a major role in regulation of eukaryotic gene expression. This process is controlled by histone acetyltransferases (HATs/KATs) found in multiprotein complexes that are recruited to chromatin by the scaffolding subunit transformation/transcription domain-associated protein (TRRAP). TRRAP is evolutionarily conserved and is among the top five genes intolerant to missense variation. Through an international collaboration, 17 distinct de novo or apparently de novo variants were identified in TRRAP in 24 individuals. A strong genotype-phenotype correlation was observed with two distinct clinical spectra. The first is a complex, multi-systemic syndrome associated with various malformations of the brain, heart, kidneys, and genitourinary system and characterized by a wide range of intellectual functioning; a number of affected individuals have intellectual disability (ID) and markedly impaired basic life functions. Individuals with this phenotype had missense variants clustering around the c.3127G>A p.(Ala1043Thr) variant identified in five individuals. The second spectrum manifested with autism spectrum disorder (ASD) and/or ID and epilepsy. Facial dysmorphism was seen in both groups and included upslanted palpebral fissures, epicanthus, telecanthus, a wide nasal bridge and ridge, a broad and smooth philtrum, and a thin upper lip. RNA sequencing analysis of skin fibroblasts derived from affected individuals skin fibroblasts showed significant changes in the expression of several genes implicated in neuronal function and ion transport. Thus, we describe here the clinical spectrum associated with TRRAP pathogenic missense variants, and we suggest a genotype-phenotype correlation useful for clinical evaluation of the pathogenicity of the variants.

De novo variants in HK1 associated with neurodevelopmental abnormalities and visual impairment.

Hexokinase 1 (HK1) phosphorylates glucose to glucose-6-phosphate, the first rate-limiting step in glycolysis. Homozygous and heterozygous variants in HK1 have been shown to cause autosomal recessive non-spherocytic hemolytic anemia, autosomal recessive Russe type hereditary motor and sensory neuropathy, and autosomal dominant retinitis pigmentosa (adRP). We report seven patients from six unrelated families with a neurodevelopmental disorder associated with developmental delay, intellectual disability, structural brain abnormality, and visual impairments in whom we identified four novel, de novo missense variants in the N-terminal half of HK1. Hexokinase activity in red blood cells of two patients was normal, suggesting that the disease mechanism is not due to loss of hexokinase enzymatic activity.