Joint, multifaceted genomic analysis enables diagnosis of diverse, ultra-rare monogenic presentations.

Genomics for rare disease diagnosis has advanced at a rapid pace due to our ability to perform "N-of-1" analyses on individual patients. The increasing sizes of ultra-rare, "N-of-1" disease cohorts internationally newly enables cohort-wide analyses for new discoveries, but well-calibrated statistical genetics approaches for jointly analyzing these patients are still under development.1,2 The Undiagnosed Diseases Network (UDN) brings multiple clinical, research and experimental centers under the same umbrella across the United States to facilitate and scale N-of-1 analyses. Here, we present the first joint analysis of whole genome sequencing data of UDN patients across the network. We apply existing and introduce new, well-calibrated statistical methods for prioritizing disease genes with de novo recurrence and compound heterozygosity. We also detect pathways enriched with candidate and known diagnostic genes. Our computational analysis, coupled with a systematic clinical review, recapitulated known diagnoses and revealed new disease associations. We make our gene-level findings and variant-level information across the cohort available in a public-facing browser (https://dbmi-bgm.github.io/udn-browser/). These results show that N-of-1 efforts should be supplemented by a joint genomic analysis across cohorts.

Low-frequency inherited complement receptor variants are associated with purpura fulminans.

ABSTRACT: Extreme disease phenotypes can provide key insights into the pathophysiology of common conditions, but studying such cases is challenging due to their rarity and the limited statistical power of existing methods. Herein, we used a novel approach to pathway-based mutational burden testing, the rare variant trend test (RVTT), to investigate genetic risk factors for an extreme form of sepsis-induced coagulopathy, infectious purpura fulminans (PF). In addition to prospective patient sample collection, we electronically screened over 10.4 million medical records from 4 large hospital systems and identified historical cases of PF for which archived specimens were available to perform germline whole-exome sequencing. We found a significantly increased burden of low-frequency, putatively function-altering variants in the complement system in patients with PF compared with unselected patients with sepsis (P = .01). A multivariable logistic regression analysis found that the number of complement system variants per patient was independently associated with PF after controlling for age, sex, and disease acuity (P = .01). Functional characterization of PF-associated variants in the immunomodulatory complement receptors CR3 and CR4 revealed that they result in partial or complete loss of anti-inflammatory CR3 function and/or gain of proinflammatory CR4 function. Taken together, these findings suggest that inherited defects in CR3 and CR4 predispose to the maladaptive hyperinflammation that characterizes severe sepsis with coagulopathy.

Macrocephaly and developmental delay caused by missense variants in RAB5C.

Rab GTPases are important regulators of intracellular vesicular trafficking. RAB5C is a member of the Rab GTPase family that plays an important role in the endocytic pathway, membrane protein recycling and signaling. Here we report on 12 individuals with nine different heterozygous de novo variants in RAB5C. All but one patient with missense variants (n = 9) exhibited macrocephaly, combined with mild-to-moderate developmental delay. Patients with loss of function variants (n = 2) had an apparently more severe clinical phenotype with refractory epilepsy and intellectual disability but a normal head circumference. Four missense variants were investigated experimentally. In vitro biochemical studies revealed that all four variants were damaging, resulting in increased nucleotide exchange rate, attenuated responsivity to guanine exchange factors and heterogeneous effects on interactions with effector proteins. Studies in C. elegans confirmed that all four variants were damaging in vivo and showed defects in endocytic pathway function. The variant heterozygotes displayed phenotypes that were not observed in null heterozygotes, with two shown to be through a dominant negative mechanism. Expression of the human RAB5C variants in zebrafish embryos resulted in defective development, further underscoring the damaging effects of the RAB5C variants. Our combined bioinformatic, in vitro and in vivo experimental studies and clinical data support the association of RAB5C missense variants with a neurodevelopmental disorder characterized by macrocephaly and mild-to-moderate developmental delay through disruption of the endocytic pathway.

A novel DPH5-related diphthamide-deficiency syndrome causing embryonic lethality or profound neurodevelopmental disorder.

PURPOSE: Diphthamide is a post-translationally modified histidine essential for messenger RNA translation and ribosomal protein synthesis. We present evidence for DPH5 as a novel cause of embryonic lethality and profound neurodevelopmental delays (NDDs). METHODS: Molecular testing was performed using exome or genome sequencing. A targeted Dph5 knockin mouse (C57BL/6Ncrl-Dph5em1Mbp/Mmucd) was created for a DPH5 p.His260Arg homozygous variant identified in 1 family. Adenosine diphosphate-ribosylation assays in DPH5-knockout human and yeast cells and in silico modeling were performed for the identified DPH5 potential pathogenic variants. RESULTS: DPH5 variants p.His260Arg (homozygous), p.Asn110Ser and p.Arg207Ter (heterozygous), and p.Asn174LysfsTer10 (homozygous) were identified in 3 unrelated families with distinct overlapping craniofacial features, profound NDDs, multisystem abnormalities, and miscarriages. Dph5 p.His260Arg homozygous knockin was embryonically lethal with only 1 subviable mouse exhibiting impaired growth, craniofacial dysmorphology, and multisystem dysfunction recapitulating the human phenotype. Adenosine diphosphate-ribosylation assays showed absent to decreased function in DPH5-knockout human and yeast cells. In silico modeling of the variants showed altered DPH5 structure and disruption of its interaction with eEF2. CONCLUSION: We provide strong clinical, biochemical, and functional evidence for DPH5 as a novel cause of embryonic lethality or profound NDDs with multisystem involvement and expand diphthamide-deficiency syndromes and ribosomopathies.

Returning actionable genomic results in a research biobank: Analytic validity, clinical implementation, and resource utilization.

Over 100 million research participants around the world have had research array-based genotyping (GT) or genome sequencing (GS), but only a small fraction of these have been offered return of actionable genomic findings (gRoR). Between 2017 and 2021, we analyzed genomic results from 36,417 participants in the Mass General Brigham Biobank and offered to confirm and return pathogenic and likely pathogenic variants (PLPVs) in 59 genes. Variant verification prior to participant recontact revealed that GT falsely identified PLPVs in 44.9% of samples, and GT failed to identify 72.0% of PLPVs detected in a subset of samples that were also sequenced. GT and GS detected verified PLPVs in 1% and 2.5% of the cohort, respectively. Of 256 participants who were alerted that they carried actionable PLPVs, 37.5% actively or passively declined further disclosure. 76.3% of those carrying PLPVs were unaware that they were carrying the variant, and over half of those met published professional criteria for genetic testing but had never been tested. This gRoR protocol cost approximately $129,000 USD per year in laboratory testing and research staff support, representing $14 per participant whose DNA was analyzed or $3,224 per participant in whom a PLPV was confirmed and disclosed. These data provide logistical details around gRoR that could help other investigators planning to return genomic results.

Heterozygous loss-of-function variants significantly expand the phenotypes associated with loss of GDF11.

PURPOSE: Growth differentiation factor 11 (GDF11) is a key signaling protein required for proper development of many organ systems. Only one prior study has associated an inherited GDF11 variant with a dominant human disease in a family with variable craniofacial and vertebral abnormalities. Here, we expand the phenotypic spectrum associated with GDF11 variants and document the nature of the variants. METHODS: We present a cohort of six probands with de novo and inherited nonsense/frameshift (4/6 patients) and missense (2/6) variants in GDF11. We generated gdf11 mutant zebrafish to model loss of gdf11 phenotypes and used an overexpression screen in Drosophila to test variant functionality. RESULTS: Patients with variants in GDF11 presented with craniofacial (5/6), vertebral (5/6), neurological (6/6), visual (4/6), cardiac (3/6), auditory (3/6), and connective tissue abnormalities (3/6). gdf11 mutant zebrafish show craniofacial abnormalities and body segmentation defects that match some patient phenotypes. Expression of the patients' variants in the fly showed that one nonsense variant in GDF11 is a severe loss-of-function (LOF) allele whereas the missense variants in our cohort are partial LOF variants. CONCLUSION: GDF11 is needed for human development, particularly neuronal development, and LOF GDF11 alleles can affect the development of numerous organs and tissues.

Clinical utility of brain MRS imaging of patients with adult-onset non-cirrhotic hyperammonemia.

Adult-onset non-cirrhotic hyperammonemia (NCH) is a rare, but often fatal condition that can result in both reversible and irreversible neurological defects. Here we present five cases of adult-onset non-cirrhotic hyperammonemia wherein brain magnetic resonance spectroscopy (MRS) scans for cerebral glutamine (Gln) and myo-inositol (mI) levels helped guide clinical management. Specifically, we demonstrate that when combined with traditional brain magnetic resonance imaging (MRI) scans, cerebral Gln and mI MRS can help disentangle the reversible from irreversible neurological defects associated with hyperammonemic crisis. Specifically, we demonstrate that whereas an elevated brain MRS Gln level is associated with reversible neurological defects, markedly low mI levels are associated with a risk for irreversible neurological defects such as central pontine myelinolysis. Overall, our findings indicate the utility of brain MRS in guiding clinical care and prognosis in patients with adult-onset non-cirrhotic hyperammonemia.

Clinical Utility of Pharmacogenomic Data Collected by a Health-System Biobank to Predict and Prevent Adverse Drug Events.

INTRODUCTION: Medication-related harm represents a significant issue for patient safety and quality of care. One strategy to avoid preventable adverse drug events is to utilize patient-specific factors such as pharmacogenomics (PGx) to individualize therapy. OBJECTIVE: We measured the number of patients enrolled in a health-system biobank with actionable PGx results who received relevant medications and assessed the incidence of adverse drug events (ADEs) that might have been prevented had the PGx results been used to inform prescribing. METHODS: Patients with actionable PGx results in the following four genes with Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines were identified: HLA-A*31:01, HLA-B*15:02, TPMT, and VKORC1. The patients who received interacting medications (carbamazepine, oxcarbazepine, thiopurines, or warfarin) were identified, and electronic health records were reviewed to determine the incidence of potentially preventable ADEs. RESULTS: Of 36,424 patients with PGx results, 2327 (6.4%) were HLA-A*31:01 positive; 3543 (9.7%) were HLA-B*15:02 positive; 2893 (7.9%) were TPMT intermediate metabolizers; and 4249 (11.7%) were homozygous for the VKORC1 c.1639 G>A variant. Among patients positive for one of the HLA variants who received carbamazepine or oxcarbazepine (n = 92), four (4.3%) experienced a rash that warranted drug discontinuation. Among the TPMT intermediate metabolizers who received a thiopurine (n = 56), 11 (19.6%) experienced severe myelosuppression that warranted drug discontinuation. Among patients homozygous for the VKORC1 c.1639 G>A variant who received warfarin (n = 379), 85 (22.4%) experienced active bleeding and/or international normalized ratio (INR) > 5 that warranted drug discontinuation or dose reduction. CONCLUSION: Patients with actionable PGx results from a health-system biobank who received relevant medications experienced predictable ADEs. These ADEs may have been prevented if the patients' PGx results were available in the electronic health record with clinical decision support prior to prescribing.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases.

PURPOSE: Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful. METHODS: We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols. RESULTS: We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases. CONCLUSION: The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

A retrospective study of adult patients with noncirrhotic hyperammonemia.

Adult-onset noncirrhotic hyperammonemia (NCH) is poorly understood and has a high morbidity and mortality. To elucidate the etiology and management of NCH, we performed a retrospective analysis of 23 adults (median age 51) with NCH treated between 2014 and 2020 at two academic medical centers. Hyperammonemia was diagnosed in all cases during the evaluation of altered mental status, with 22% presenting with seizures. Peak ammonia levels were >200 µmol/L in 70% of cases. Defects in ammonia metabolism were assessed using urea cycle biochemical testing, germline genetic testing, and testing for urease-producing infectious agents. Ammonia metabolism defects in these cases appear attributable to four major sources: (a) infection with urease-producing organism (n = 5); (b) previously undiagnosed inborn errors of metabolism (IEMs) (n = 4); (c) clinical exposures causing acquired urea cycle dysfunction (n = 6); and (d) unexplained acquired urea cycle dysfunction (uaUCD) (n = 8), as evidenced by biochemical signatures of urea cycle dysfunction without a genetic or clinical exposure. Severe protein malnutrition appeared to be a reversible risk factor for uaUCD. Overall, 13% of our cohort died prior to resolution of hyperammonemia, 26% died after hyperammonemia resolution, 57% survived after having reversible neurological changes, and 4% survived with irreversible neurological changes. Renal replacement therapy for ammonia clearance was often utilized for patients with an ammonia level above 250 µmol/L and patients were frequently empirically treated with antibiotics targeting urea-splitting organisms. Our study demonstrates that acquired urea cycle dysfunction, IEMs and urease-producing infections are major sources of adult-onset NCH and highlights successful management strategies for adult-onset NCH.

De Novo Variants in the ATPase Module of MORC2 Cause a Neurodevelopmental Disorder with Growth Retardation and Variable Craniofacial Dysmorphism.

MORC2 encodes an ATPase that plays a role in chromatin remodeling, DNA repair, and transcriptional regulation. Heterozygous variants in MORC2 have been reported in individuals with autosomal-dominant Charcot-Marie-Tooth disease type 2Z and spinal muscular atrophy, and the onset of symptoms ranges from infancy to the second decade of life. Here, we present a cohort of 20 individuals referred for exome sequencing who harbor pathogenic variants in the ATPase module of MORC2. Individuals presented with a similar phenotype consisting of developmental delay, intellectual disability, growth retardation, microcephaly, and variable craniofacial dysmorphism. Weakness, hyporeflexia, and electrophysiologic abnormalities suggestive of neuropathy were frequently observed but were not the predominant feature. Five of 18 individuals for whom brain imaging was available had lesions reminiscent of those observed in Leigh syndrome, and five of six individuals who had dilated eye exams had retinal pigmentary abnormalities. Functional assays revealed that these MORC2 variants result in hyperactivation of epigenetic silencing by the HUSH complex, supporting their pathogenicity. The described set of morphological, growth, developmental, and neurological findings and medical concerns expands the spectrum of genetic disorders resulting from pathogenic variants in MORC2.

De novo mutations in TOMM70, a receptor of the mitochondrial import translocase, cause neurological impairment.

The translocase of outer mitochondrial membrane (TOMM) complex is the entry gate for virtually all mitochondrial proteins and is essential to build the mitochondrial proteome. TOMM70 is a receptor that assists mainly in mitochondrial protein import. Here, we report two individuals with de novo variants in the C-terminal region of TOMM70. While both individuals exhibited shared symptoms including hypotonia, hyper-reflexia, ataxia, dystonia and significant white matter abnormalities, there were differences between the two individuals, most prominently the age of symptom onset. Both individuals were undiagnosed despite extensive genetics workups. Individual 1 was found to have a p.Thr607Ile variant while Individual 2 was found to have a p.Ile554Phe variant in TOMM70. To functionally assess both TOMM70 variants, we replaced the Drosophila Tom70 coding region with a Kozak-mini-GAL4 transgene using CRISPR-Cas9. Homozygous mutant animals die as pupae, but lethality is rescued by the mini-GAL4-driven expression of human UAS-TOMM70 cDNA. Both modeled variants lead to significantly less rescue indicating that they are loss-of-function alleles. Similarly, RNAi-mediated knockdown of Tom70 in the developing eye causes roughening and synaptic transmission defect, common findings in neurodegenerative and mitochondrial disorders. These phenotypes were rescued by the reference, but not the variants, of TOMM70. Altogether, our data indicate that de novo loss-of-function variants in TOMM70 result in variable white matter disease and neurological phenotypes in affected individuals.

Identification of an Identical de Novo SCAMP5 Missense Variant in Four Unrelated Patients With Seizures and Severe Neurodevelopmental Delay.

Objective: To establish and broaden the phenotypic spectrum of secretory carrier membrane protein (SCAMP5) associated with epilepsy and neurodevelopmental delay. Methods: A Chinese patient was identified at the First Hospital of Peking University, and the three unrelated patients were recruited from two different countries (Italy and United States) through GeneMatcher. SCAMP5 pathogenic variants were identified by whole exome sequencing; clinical data of the patients were retrospectively collected and analyzed. Result: The onset age of seizures was ranged from 6 to 15 months. Patients had different types of seizures, including focal seizures, generalized tonic-clonic seizures and tonic seizure. One patient showed typical autism spectrum disorder (ASD) symptoms. Electroencephalogram (EEG) findings presented as focal or multifocal discharges, sometimes spreading to generalization. Brain magnetic resonance imaging (MRI) abnormalities were present in each patient. Severe intellectual disability and language and motor developmental disorders were found in our patients, with all patients having poor language development and were nonverbal at last follow-up. All but one of the patients could walk independently in childhood, but the ability to walk independently in one patient had deteriorated with age. All patients had abnormal neurological exam findings, mostly signs of extrapyramidal system involvement. Dysmorphic features were found in 2/4 patients, mainly in the face and trunk. All four unrelated patients were found to have the same heterozygous pathogenic SCAMP5 de novo variant (p. Gly180Trp). Conclusion: Epilepsy, severe developmental delay, abnormal neurological exam findings, with or without ASD or variably dysmorphic features and were common in patients with SCAMP5 variant. The onset time and type of seizure varied greatly. The EEG and brain MRI findings were not consistent, but diverse and nonspecific. The motor ability of patients with heterozygous SCAMP5 variant might have a regressive course; language development was more severely affected.

De Novo Pathogenic Variants in N-cadherin Cause a Syndromic Neurodevelopmental Disorder with Corpus Collosum, Axon, Cardiac, Ocular, and Genital Defects.

Cadherins constitute a family of transmembrane proteins that mediate calcium-dependent cell-cell adhesion. The extracellular domain of cadherins consists of extracellular cadherin (EC) domains, separated by calcium binding sites. The EC interacts with other cadherin molecules in cis and in trans to mechanically hold apposing cell surfaces together. CDH2 encodes N-cadherin, whose essential roles in neural development include neuronal migration and axon pathfinding. However, CDH2 has not yet been linked to a Mendelian neurodevelopmental disorder. Here, we report de novo heterozygous pathogenic variants (seven missense, two frameshift) in CDH2 in nine individuals with a syndromic neurodevelopmental disorder characterized by global developmental delay and/or intellectual disability, variable axon pathfinding defects (corpus callosum agenesis or hypoplasia, mirror movements, Duane anomaly), and ocular, cardiac, and genital anomalies. All seven missense variants (c.1057G>A [p.Asp353Asn]; c.1789G>A [p.Asp597Asn]; c.1789G>T [p.Asp597Tyr]; c.1802A>C [p.Asn601Thr]; c.1839C>G [p.Cys613Trp]; c.1880A>G [p.Asp627Gly]; c.2027A>G [p.Tyr676Cys]) result in substitution of highly conserved residues, and six of seven cluster within EC domains 4 and 5. Four of the substitutions affect the calcium-binding site in the EC4-EC5 interdomain. We show that cells expressing these variants in the EC4-EC5 domains have a defect in cell-cell adhesion; this defect includes impaired binding in trans with N-cadherin-WT expressed on apposing cells. The two frameshift variants (c.2563_2564delCT [p.Leu855Valfs∗4]; c.2564_2567dupTGTT [p.Leu856Phefs∗5]) are predicted to lead to a truncated cytoplasmic domain. Our study demonstrates that de novo heterozygous variants in CDH2 impair the adhesive activity of N-cadherin, resulting in a multisystemic developmental disorder, that could be named ACOG syndrome (agenesis of corpus callosum, axon pathfinding, cardiac, ocular, and genital defects).

Yield of whole exome sequencing in undiagnosed patients facing insurance coverage barriers to genetic testing.

BACKGROUND: Despite growing evidence of diagnostic yield and clinical utility of whole exome sequencing (WES) in patients with undiagnosed diseases, there remain significant cost and reimbursement barriers limiting access to such testing. The diagnostic yield and resulting clinical actions of WES for patients who previously faced insurance coverage barriers have not yet been explored. METHODS: We performed a retrospective descriptive analysis of clinical WES outcomes for patients facing insurance coverage barriers prior to clinical WES and who subsequently enrolled in the Undiagnosed Diseases Network (UDN). Clinical WES was completed as a result of participation in the UDN. Payer type, molecular diagnostic yield, and resulting clinical actions were evaluated. RESULTS: Sixty-six patients in the UDN faced insurance coverage barriers to WES at the time of enrollment (67% public payer, 26% private payer). Forty-two of 66 (64%) received insurance denial for clinician-ordered WES, 19/66 (29%) had health insurance through a payer known not to cover WES, and 5/66 (8%) had previous payer denial of other genetic tests. Clinical WES results yielded a molecular diagnosis in 23 of 66 patients (35% [78% pediatric, 65% neurologic indication]). Molecular diagnosis resulted in clinical actions in 14 of 23 patients (61%). CONCLUSIONS: These data demonstrate that a substantial proportion of patients who encountered insurance coverage barriers to WES had a clinically actionable molecular diagnosis, supporting the notion that WES has value as a covered benefit for patients who remain undiagnosed despite objective clinical findings.