Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Digenic Leigh syndrome on the background of the m.11778G>A Leber hereditary optic neuropathy variant.

Leigh syndrome spectrum (LSS) is a primary mitochondrial disorder defined neuropathologically by a subacute necrotizing encephalomyelopathy and characterized by bilateral basal ganglia and/or brainstem lesions. LSS is associated with variants in several mitochondrial DNA genes and more than 100 nuclear genes, most often related to mitochondrial complex I (CI) dysfunction. Rarely, LSS has been reported in association with primary Leber hereditary optic neuropathy (LHON) variants of the mitochondrial DNA, coding for CI subunits (m.3460G>A in MT-ND1, m.11778G>A in MT-ND4 and m.14484T>C in MT-ND6). The underlying mechanism by which these variants manifest as LSS, a severe neurodegenerative disease, as opposed to the LHON phenotype of isolated optic neuropathy, remains an open question. Here, we analyse the exome sequencing of six probands with LSS carrying primary LHON variants, and report digenic co-occurrence of the m.11778G > A variant with damaging heterozygous variants in nuclear disease genes encoding CI subunits as a plausible explanation. Our findings suggest a digenic mechanism of disease for m.11778G>A-associated LSS, consistent with recent reports of digenic disease in individuals manifesting with LSS due to biallelic variants in the recessive LHON-associated disease gene DNAJC30 in combination with heterozygous variants in CI subunits.

Narrowing the diagnostic gap: Genomes, episignatures, long-read sequencing, and health economic analyses in an exome-negative intellectual disability cohort.

PURPOSE: Genome sequencing (GS)-specific diagnostic rates in prospective tightly ascertained exome sequencing (ES)-negative intellectual disability (ID) cohorts have not been reported extensively. METHODS: ES, GS, epigenetic signatures, and long-read sequencing diagnoses were assessed in 74 trios with at least moderate ID. RESULTS: The ES diagnostic yield was 42 of 74 (57%). GS diagnoses were made in 9 of 32 (28%) ES-unresolved families. Repeated ES with a contemporary pipeline on the GS-diagnosed families identified 8 of 9 single-nucleotide variations/copy-number variations undetected in older ES, confirming a GS-unique diagnostic rate of 1 in 32 (3%). Episignatures contributed diagnostic information in 9% with GS corroboration in 1 of 32 (3%) and diagnostic clues in 2 of 32 (6%). A genetic etiology for ID was detected in 51 of 74 (69%) families. Twelve candidate disease genes were identified. Contemporary ES followed by GS cost US$4976 (95% CI: $3704; $6969) per diagnosis and first-line GS at a cost of $7062 (95% CI: $6210; $8475) per diagnosis. CONCLUSION: Performing GS only in ID trios would be cost equivalent to ES if GS were available at $2435, about a 60% reduction from current prices. This study demonstrates that first-line GS achieves higher diagnostic rate than contemporary ES but at a higher cost.

Leigh Syndrome: A Study of 209 Patients at the Beijing Children's Hospital.

OBJECTIVE: Leigh syndrome (LS) is a heterogeneous neurodegenerative disease and the most frequent pediatric manifestation of mitochondrial disease. In the largest patient collection to date, this study aimed to provide new insights into the clinical and genetic spectrum of LS, defect-specific associations, and predictors of disease course and survival. METHODS: Clinical, metabolic, neuroimaging, onset, and survival data were collected from the medical records of 209 patients referred to the Beijing Children's Hospital with symmetrical basal ganglia and/or brainstem neuroimaging changes indicative of LS by 30 centers from the Chinese network of mitochondrial disease (mitoC-NET) between January 2013 and July 2021 for exploratory analysis. RESULTS: Pathogenic variants were identified in 52 genes, most frequently MT-ATP6, SURF1, and PDHA1. Maternally inherited variants accounted for 42% (heteroplasmy level ≥90% in 64%). Phenotypes spanned 92 Human Phenotype Ontology terms. Elevated serum lactate (144/195), global developmental delay (142/209), and developmental regression (103/209) were most frequent. Discriminating neuroimaging and/or clinical features were identified for MT-ATP6 (m.9176T>C), MT-ND5, PDHA1, SUCLG1, and SURF1. Poorest survival was associated with MT-ND5, MT-ATP6 (m.8993T>C and m.9176T>C), SURF1, and ALDH5A1 (≤50% 3 year's survival), in contrast to milder defects with specific treatment (ECHS1 and SLC19A3, 100% 3 year's survival). INTERPRETATION: Our data define phenotype, onset, and survival of LS in a defect-specific manner, identifying features discriminating between genetic defects and predictive of disease outcome. These findings are essential to early diagnosis, in optimizing family counseling, and to the design and monitoring of future clinical trials, the next frontier of LS research. ANN NEUROL 2022;91:466-482.

DNAJC30 defect: a frequent cause of recessive Leber hereditary optic neuropathy and Leigh syndrome.

The recent description of biallelic DNAJC30 variants in Leber hereditary optic neuropathy (LHON) and Leigh syndrome challenged the longstanding assumption for LHON to be exclusively maternally inherited and broadened the genetic spectrum of Leigh syndrome, the most frequent paediatric mitochondrial disease. Herein, we characterize 28 so far unreported individuals from 26 families carrying a homozygous DNAJC30 p.Tyr51Cys founder variant, 24 manifesting with LHON, two manifesting with Leigh syndrome, and two remaining asymptomatic. This collection of unreported variant carriers confirms sex-dependent incomplete penetrance of the homozygous variant given a significant male predominance of disease and the report of asymptomatic homozygous variant carriers. The autosomal recessive LHON patients demonstrate an earlier age of disease onset and a higher rate of idebenone-treated and spontaneous recovery of vision in comparison to reported figures for maternally inherited disease. Moreover, the report of two additional patients with childhood- or adult-onset Leigh syndrome further evidences the association of DNAJC30 with Leigh syndrome, previously only reported in a single childhood-onset case.

Oral Coenzyme Q10 supplementation leads to better preservation of kidney function in steroid-resistant nephrotic syndrome due to primary Coenzyme Q10 deficiency.

Primary Coenzyme Q10 (CoQ10) deficiency is an ultra-rare disorder caused by defects in genes involved in CoQ10 biosynthesis leading to multidrug-resistant nephrotic syndrome as the hallmark kidney manifestation. Promising early results have been reported anecdotally with oral CoQ10 supplementation. However, the long-term efficacy and optimal prescription remain to be established. In a global effort, we collected and analyzed information from 116 patients who received CoQ10 supplements for primary CoQ10 deficiency due to biallelic pathogenic variants in either the COQ2, COQ6 or COQ8B genes. Median duration of follow up on treatment was two years. The effect of treatment on proteinuria was assessed, and kidney survival was analyzed in 41 patients younger than 18 years with chronic kidney disease stage 1-4 at the start of treatment compared with that of an untreated cohort matched by genotype, age, kidney function, and proteinuria. CoQ10 supplementation was associated with a substantial and significant sustained reduction of proteinuria by 88% at 12 months. Complete remission of proteinuria was more frequently observed in COQ6 disease. CoQ10 supplementation led to significantly better preservation of kidney function (5-year kidney failure-free survival 62% vs. 19%) with an improvement in general condition and neurological manifestations. Side effects of treatment were uncommon and mild. Thus, our findings indicate that all patients diagnosed with primary CoQ10 deficiency should receive early and life-long CoQ10 supplementation to decelerate the progression of kidney disease and prevent further damage to other organs.

Variation of the clinical spectrum and genotype-phenotype associations in Coenzyme Q10 deficiency associated glomerulopathy.

Primary Coenzyme Q10 deficiency is a rare mitochondriopathy with a wide spectrum of organ involvement, including steroid-resistant nephrotic syndrome mainly associated with disease-causing variants in the genes COQ2, COQ6 or COQ8B. We performed a systematic literature review, PodoNet, mitoNET, and CCGKDD registries queries and an online survey, collecting comprehensive clinical and genetic data of 251 patients spanning 173 published (47 updated) and 78 new cases. Kidney disease was first diagnosed at median age 1.0, 1.2 and 9.8 years in individuals with disease-causing variants in COQ2, COQ6 and COQ8B, respectively. Isolated kidney involvement at diagnosis occurred in 34% of COQ2, 10.8% of COQ6 and 70.7% of COQ8B variant individuals. Classic infantile multiorgan involvement comprised 22% of the COQ2 variant cohort while 47% of them developed neurological symptoms at median age 2.7 years. The association of steroid-resistant nephrotic syndrome and sensorineural hearing loss was confirmed as the distinctive phenotype of COQ6 variants, with hearing impairment manifesting at average age three years. None of the patients with COQ8B variants, but 50% of patients with COQ2 and COQ6 variants progressed to kidney failure by age five. At adult age, kidney survival was equally poor (20-25%) across all disorders. A number of sequence variants, including putative local founder mutations, had divergent clinical presentations, in terms of onset age, kidney and non-kidney manifestations and kidney survival. Milder kidney phenotype was present in those with biallelic truncating variants within the COQ8B variant cohort. Thus, significant intra- and inter-familial phenotype variability was observed, suggesting both genetic and non-genetic modifiers of disease severity.

The genetics of mitochondrial disease: dissecting mitochondrial pathology using multi-omic pipelines.

Mitochondria play essential roles in numerous metabolic pathways including the synthesis of adenosine triphosphate through oxidative phosphorylation. Clinically, mitochondrial diseases occur when there is mitochondrial dysfunction - manifesting at any age and affecting any organ system; tissues with high energy requirements, such as muscle and the brain, are often affected. The clinical heterogeneity is parallel to the degree of genetic heterogeneity associated with mitochondrial dysfunction. Around 10% of human genes are predicted to have a mitochondrial function, and defects in over 300 genes are reported to cause mitochondrial disease. Some involve the mitochondrial genome (mtDNA), but the vast majority occur within the nuclear genome. Except for a few specific genetic defects, there remains no cure for mitochondrial diseases, which means that a genetic diagnosis is imperative for genetic counselling and the provision of reproductive options for at-risk families. Next-generation sequencing strategies, particularly exome and whole-genome sequencing, have revolutionised mitochondrial diagnostics such that the traditional muscle biopsy has largely been replaced with a minimally-invasive blood sample for an unbiased approach to genetic diagnosis. Where these genomic approaches have not identified a causative defect, or where there is insufficient support for pathogenicity, additional functional investigations are required. The application of supplementary 'omics' technologies, including transcriptomics, proteomics, and metabolomics, has the potential to greatly improve diagnostic strategies. This review aims to demonstrate that whilst a molecular diagnosis can be achieved for many cases through next-generation sequencing of blood DNA, the use of patient tissues and an integrated, multidisciplinary multi-omics approach is pivotal for the diagnosis of more challenging cases. Moreover, the analysis of clinically relevant tissues from affected individuals remains crucial for understanding the molecular mechanisms underlying mitochondrial pathology. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Expanding the clinical and genetic spectrum of FDXR deficiency by functional validation of variants of uncertain significance.

Ferrodoxin reductase (FDXR) deficiency is a mitochondrial disease described in recent years primarily in association with optic atrophy, acoustic neuropathy, and developmental delays. Here, we identified seven unpublished patients with FDXR deficiency belonging to six independent families. These patients show a broad clinical spectrum ranging from Leigh syndrome with early demise and severe infantile-onset encephalopathy, to milder movement disorders. In total nine individual pathogenic variants, of which seven were novel, were identified in FDXR using whole exome sequencing in suspected mitochondrial disease patients. Over 80% of these variants are missense, a challenging variant class in which to determine pathogenic consequence, especially in the setting of nonspecific phenotypes and in the absence of a reliable biomarker, necessitating functional validation. Here we implement an Arh1-null yeast model to confirm the pathogenicity of variants of uncertain significance in FDXR, bypassing the requirement for patient-derived material.

The diagnosis of inborn errors of metabolism by an integrative "multi-omics" approach: A perspective encompassing genomics, transcriptomics, and proteomics.

Given the rapidly decreasing cost and increasing speed and accessibility of massively parallel technologies, the integration of comprehensive genomic, transcriptomic, and proteomic data into a "multi-omics" diagnostic pipeline is within reach. Even though genomic analysis has the capability to reveal all possible perturbations in our genetic code, analysis typically reaches a diagnosis in just 35% of cases, with a diagnostic gap arising due to limitations in prioritization and interpretation of detected variants. Here we review the utility of complementing genetic data with transcriptomic data and give a perspective for the introduction of proteomics into the diagnostic pipeline. Together these methodologies enable comprehensive capture of the functional consequence of variants, unobtainable by the analysis of each methodology in isolation. This facilitates functional annotation and reprioritization of candidate genes and variants-a promising approach to shed light on the underlying molecular cause of a patient's disease, increasing diagnostic rate, and allowing actionability in clinical practice.

Biallelic DMXL2 mutations impair autophagy and cause Ohtahara syndrome with progressive course.

Ohtahara syndrome, early infantile epileptic encephalopathy with a suppression burst EEG pattern, is an aetiologically heterogeneous condition starting in the first weeks or months of life with intractable seizures and profound developmental disability. Using whole exome sequencing, we identified biallelic DMXL2 mutations in three sibling pairs with Ohtahara syndrome, belonging to three unrelated families. Siblings in Family 1 were compound heterozygous for the c.5135C>T (p.Ala1712Val) missense substitution and the c.4478C>G (p.Ser1493*) nonsense substitution; in Family 2 were homozygous for the c.4478C>A (p.Ser1493*) nonsense substitution and in Family 3 were homozygous for the c.7518-1G>A (p.Trp2507Argfs*4) substitution. The severe developmental and epileptic encephalopathy manifested from the first day of life and was associated with deafness, mild peripheral polyneuropathy and dysmorphic features. Early brain MRI investigations in the first months of life revealed thin corpus callosum with brain hypomyelination in all. Follow-up MRI scans in three patients revealed progressive moderate brain shrinkage with leukoencephalopathy. Five patients died within the first 9 years of life and none achieved developmental, communicative or motor skills following birth. These clinical findings are consistent with a developmental brain disorder that begins in the prenatal brain, prevents neural connections from reaching the expected stages at birth, and follows a progressive course. DMXL2 is highly expressed in the brain and at synaptic terminals, regulates v-ATPase assembly and activity and participates in intracellular signalling pathways; however, its functional role is far from complete elucidation. Expression analysis in patient-derived skin fibroblasts demonstrated absence of the DMXL2 protein, revealing a loss of function phenotype. Patients' fibroblasts also exhibited an increased LysoTracker® signal associated with decreased endolysosomal markers and degradative processes. Defective endolysosomal homeostasis was accompanied by impaired autophagy, revealed by lower LC3II signal, accumulation of polyubiquitinated proteins, and autophagy receptor p62, with morphological alterations of the autolysosomal structures on electron microscopy. Altered lysosomal homeostasis and defective autophagy were recapitulated in Dmxl2-silenced mouse hippocampal neurons, which exhibited impaired neurite elongation and synaptic loss. Impaired lysosomal function and autophagy caused by biallelic DMXL2 mutations affect neuronal development and synapse formation and result in Ohtahara syndrome with profound developmental impairment and reduced life expectancy.

Assessment of the evidence yield for the calibrated PP3/BP4 computational recommendations.

Purpose: To investigate the number of rare missense variants observed in human genome sequences by ACMG/AMP PP3/BP4 evidence strength, following the calibrated PP3/BP4 computational recommendations. Methods: Missense variants from the genome sequences of 300 probands from the Rare Genomes Project with suspected rare disease were analyzed using computational prediction tools able to reach PP3_Strong and BP4_Moderate evidence strengths (BayesDel, MutPred2, REVEL, and VEST4). The numbers of variants at each evidence strength were analyzed across disease-associated genes and genome-wide. Results: From a median of 75.5 rare (≤1% allele frequency) missense variants in disease-associated genes per proband, a median of one reached PP3_Strong, 3-5 PP3_Moderate, and 3-5 PP3_Supporting. Most were allocated BP4 evidence (median 41-49 per proband) or were indeterminate (median 17.5-19 per proband). Extending the analysis to all protein-coding genes genome-wide, the number of PP3_Strong variants increased approximately 2.6-fold compared to disease-associated genes, with a median per proband of 1-3 PP3_Strong, 8-16 PP3_Moderate, and 10-17 PP3_Supporting. Conclusion: A small number of variants per proband reached PP3_Strong and PP3_Moderate in 3,424 disease-associated genes, and though not the intended use of the recommendations, also genome-wide. Use of PP3/BP4 evidence as recommended from calibrated computational prediction tools in the clinical diagnostic laboratory is unlikely to inappropriately contribute to the classification of an excessive number of variants as Pathogenic or Likely Pathogenic by ACMG/AMP rules.

Exploring penetrance of clinically relevant variants in over 800,000 humans from the Genome Aggregation Database.

Incomplete penetrance, or absence of disease phenotype in an individual with a disease-associated variant, is a major challenge in variant interpretation. Studying individuals with apparent incomplete penetrance can shed light on underlying drivers of altered phenotype penetrance. Here, we investigate clinically relevant variants from ClinVar in 807,162 individuals from the Genome Aggregation Database (gnomAD), demonstrating improved representation in gnomAD version 4. We then conduct a comprehensive case-by-case assessment of 734 predicted loss of function variants (pLoF) in 77 genes associated with severe, early-onset, highly penetrant haploinsufficient disease. We identified explanations for the presumed lack of disease manifestation in 701 of the variants (95%). Individuals with unexplained lack of disease manifestation in this set of disorders rarely occur, underscoring the need and power of deep case-by-case assessment presented here to minimize false assignments of disease risk, particularly in unaffected individuals with higher rates of secondary properties that result in rescue.

Inferring compound heterozygosity from large-scale exome sequencing data.

Recessive diseases arise when both copies of a gene are impacted by a damaging genetic variant. When a patient carries two potentially causal variants in a gene, accurate diagnosis requires determining that these variants occur on different copies of the chromosome (that is, are in trans) rather than on the same copy (that is, in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. Here we developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in the Genome Aggregation Database (v2, n = 125,748 exomes). Our approach estimates phase with 96% accuracy, both in trio data and in patients with Mendelian conditions and presumed causal compound heterozygous variants. We provide a public resource of phasing estimates for coding variants and counts per gene of rare variants in trans that can aid interpretation of rare co-occurring variants in the context of recessive disease.

De novo variants in the non-coding spliceosomal snRNA gene RNU4-2 are a frequent cause of syndromic neurodevelopmental disorders.

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes1. Increasingly, large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a novel syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome2. We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 119 individuals with NDD. The vast majority of individuals (77.3%) have the same highly recurrent single base-pair insertion (n.64_65insT). We estimate that variants in this region explain 0.41% of individuals with NDD. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to its contiguous counterpart RNU4-1 and other U4 homologs, supporting RNU4-2's role as the primary U4 transcript in the brain. Overall, this work underscores the importance of non-coding genes in rare disorders. It will provide a diagnosis to thousands of individuals with NDD worldwide and pave the way for the development of effective treatments for these individuals.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.