Lessons learned from additional research analyses of unsolved clinical exome cases

PURA and GLUT1: Sweet partners for brain health

PURA, also known as Pur-alpha, is an evolutionarily conserved DNA/RNA-binding protein crucial for various cellular processes, including DNA replication, transcriptional regulation, and translational control. Comprising three PUR domains, it engages with nucleic acids and has a role in protein-protein interactions. The manifestation of PURA syndrome, arising from mutations in the PURA gene, presents neurologically with developmental delay, hypotonia, and seizures. In our prior work from 2018, we highlighted the unique case of a PURA patient displaying hypoglycorrhachia, suggesting a potential association with GLUT1 dysfunction in this syndrome. In this current study, we expand the patient cohort with PURA mutations exhibiting hypoglycorrhachia and aim to unravel the molecular basis of this phenomenon. We established an in vitro model in HeLa cells to modulate PURA expression and investigated GLUT1 function and expression. Our findings indicate that PURA levels directly impact glucose uptake through the functioning of GLUT1, without influencing significantly GLUT1 expression. Moreover, our study reveals evidence for a possible physical interaction between PURA and GLUT1, demonstrated by colocalization and co-immunoprecipitation of both proteins. Computational analyses, employing molecular dynamics, further corroborates these findings, demonstrating that PURA:GLUT1 interactions are plausible, and that the stability of the complex is altered when PURA is truncated and/or mutated. In conclusion, our results suggest that PURA plays a pivotal role in driving the function of GLUT1 for glucose uptake, potentially forming a regulatory complex. Additional investigations are warranted to elucidate the precise mechanisms governing this complex and its significance in ensuring proper GLUT1 function.

Exome and genome sequencing in a heterogeneous population of patients with rare disease: Identifying predictors of a diagnosis

PURPOSE

Exome (ES) and genome sequencing (GS) are increasingly being utilized for individuals with rare and undiagnosed diseases; however, guidelines on their use remain limited. This study aimed to identify factors associated with diagnosis by ES and/or GS in a heterogeneous population of patients with rare and undiagnosed diseases.

METHODS

In this case control study, we reviewed data from 400 diagnosed and 400 undiagnosed randomly selected participants in the Undiagnosed Diseases Network, all of whom had undergone ES and/or GS. We analyzed factors associated with receiving a diagnosis by ES and/or GS.

RESULTS

Factors associated with a decreased odds of being diagnosed included adult symptom onset, singleton sequencing, and having undergone ES and/or GS before acceptance to the Undiagnosed Diseases Network (48%, 51%, and 32% lower odds, respectively). Factors that increased the odds of being diagnosed by ES and/or GS included having primarily neurological symptoms and having undergone prior chromosomal microarray testing (44% and 59% higher odds, respectively).

CONCLUSION

We identified several factors that were associated with receiving a diagnosis by ES and/or GS. This will ideally inform the utilization of ES and/or GS and help manage expectations of individuals and families undergoing these tests.

High number of candidate gene variants are identified as disease-causing in a period of 4 years

Advances in bioinformatic tools paired with the ongoing accumulation of genetic knowledge and periodic reanalysis of genomic sequencing data have led to an improvement in genetic diagnostic rates. Candidate gene variants (CGVs) identified during sequencing or on reanalysis but not yet implicated in human disease or associated with a phenotypically distinct condition are often not revisited, leading to missed diagnostic opportunities. Here, we revisited 33 such CGVs from our previously published study and determined that 16 of them are indeed disease-causing (novel or phenotype expansion) since their identification. These results emphasize the need to focus on previously identified CGVs during sequencing or reanalysis and the importance of sharing that information with researchers around the world, including relevant functional analysis to establish disease causality.

Non-Specific Epileptic Activity, EEG, and Brain Imaging in Loss of Function Variants in SATB1: A New Case Report and Review of the Literature

SATB1 (MIM #602075) is a relatively new gene reported only in recent years in association with neurodevelopmental disorders characterized by variable facial dysmorphisms, global developmental delay, poor or absent speech, altered electroencephalogram (EEG), and brain abnormalities on imaging. To date about thirty variants in forty-four patients/children have been described, with a heterogeneous spectrum of clinical manifestations. In the present study, we describe a new patient affected by mild intellectual disability, speech disorder, and non-specific abnormalities on EEG and neuroimaging. Family studies identified a new de novo frameshift variant c.1818delG (p.(Gln606Hisfs*101)) in SATB1. To better define genotype-phenotype associations in the different types of reported SATB1 variants, we reviewed clinical data from our patient and from the literature and compared manifestations (epileptic activity, EEG abnormalities and abnormal brain imaging) due to missense variants versus those attributable to loss-of-function/premature termination variants. Our analyses showed that the latter variants are associated with less severe, non-specific clinical features when compared with the more severe phenotypes due to missense variants. These findings provide new insights into SATB1-related disorders.

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

Genomic Balancing Act: Deciphering DNA rearrangements in the Complex Chromosomal Aberration involving 5p15.2, 2q31.1 and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long read genome sequencing unveiled the presence of four distinct junctions in the father, who is carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndrome and to better understand genetic disease mechanisms.

A comprehensive review on DDX3X liquid phase condensation in health and neurodevelopmental disorders

DEAD-box helicases are global regulators of liquid-liquid phase separation (LLPS), a process that assembles membraneless organelles inside cells. An outstanding member of the DEAD-box family is DDX3X, a multi-functional protein that plays critical roles in RNA metabolism, including RNA transcription, splicing, nucleocytoplasmic export, and translation. The diverse functions of DDX3X result from its ability to bind and remodel RNA in an ATP-dependent manner. This capacity enables the protein to act as an RNA chaperone and an RNA helicase, regulating ribonucleoprotein complex assembly. DDX3X and its orthologs from mouse, yeast (Ded1), and C. elegans (LAF-1) can undergo LLPS, driving the formation of neuronal granules, stress granules, processing bodies or P-granules. DDX3X has been related to several human conditions, including neurodevelopmental disorders, such as intellectual disability and autism spectrum disorder. Although the research into the pathogenesis of aberrant biomolecular condensation in neurodegenerative diseases is increasing rapidly, the role of LLPS in neurodevelopmental disorders is underexplored. This review summarizes current findings relevant for DDX3X phase separation in neurodevelopment and examines how disturbances in the LLPS process can be related to neurodevelopmental disorders.

Reanalysis of whole-exome sequencing (WES) data of children with neurodevelopmental disorders in a standard patient care context

This study aims to inform future genetic reanalysis management by evaluating the yield of whole-exome sequencing (WES) reanalysis in standard patient care in the Netherlands. Single-center data of 159 patients with a neurodevelopmental disorder (NDD), in which WES analysis and reanalysis were performed between January 1, 2014, and December 31, 2021, was retrospectively collected. Patients were included if they were under the age of 18 years at initial analysis and if this initial analysis did not result in a diagnosis. Demographic, phenotypic, and genotypic characteristics of patients were collected and analyzed. The primary outcomes of our study were (i) diagnostic yield at reanalysis, (ii) reasons for detecting a new possibly causal variant at reanalysis, (iii) unsolicited findings, and (iv) factors associated with positive result of reanalysis. In addition, we conducted a questionnaire study amongst the 7 genetic department in the Netherlands creating an overview of used techniques, yield, and organization of WES reanalysis. The single-center data show that in most cases, WES reanalysis was initiated by the clinical geneticist (65%) or treating physician (30%). The mean time between initial WES analysis and reanalysis was 3.7 years. A new (likely) pathogenic variant or VUS with a clear link to the phenotype was found in 20 initially negative cases, resulting in a diagnostic yield of 12.6%. In 75% of these patients, the diagnosis had clinical consequences, as for example, a screening plan for associated signs and symptoms could be devised. Most (32%) of the (likely) causal variants identified at WES reanalysis were discovered due to a newly described gene-disease association. In addition to the 12.6% diagnostic yield based on new diagnoses, reclassification of a variant of uncertain significance found at initial analysis led to a definite diagnosis in three patients. Diagnostic yield was higher in patients with dysmorphic features compared to patients without clear dysmorphic features (yield 27% vs. 6%; p = 0.001).

CONCLUSIONS

Our results show that WES reanalysis in patients with NDD in standard patient care leads to a substantial increase in genetic diagnoses. In the majority of newly diagnosed patients, the diagnosis had clinical consequences. Knowledge about the clinical impact of WES reanalysis, clinical characteristics associated with higher yield, and the yield per year after a negative WES in larger clinical cohorts is warranted to inform guidelines for genetic reanalysis. These guidelines will be of great value for pediatricians, pediatric rehabilitation specialists, and pediatric neurologists in daily care of patients with NDD.

WHAT IS KNOWN

• Whole exome sequencing can cost-effectively identify a genetic cause of intellectual disability in about 30-40% of patients. • WES reanalysis in a research setting can lead to a definitive diagnosis in 10-20% of previously exome negative cases.

WHAT IS NEW

• WES reanalysis in standard patient care resulted in a diagnostic yield of 13% in previously exome negative children with NDD. • The presence of dysmorphic features is associated with an increased diagnostic yield of WES reanalysis.

Re-evaluation and re-analysis of 152 research exomes five years after the initial report reveals clinically relevant changes in 18

Iterative re-analysis of NGS results is not well investigated for published research cohorts of rare diseases. We revisited a cohort of 152 consanguineous families with developmental disorders (NDD) reported five years ago. We re-evaluated all reported variants according to diagnostic classification guidelines or our candidate gene scoring system (AutoCaSc) and systematically scored the validity of gene-disease associations (GDA). Sequencing data was re-processed using an up-to-date pipeline for case-level re-analysis. In 28/152 (18%) families, we identified a clinically relevant change. Ten previously reported (likely) pathogenic variants were re-classified as VUS/benign. In one case, the GDA (TSEN15) validity was judged as limited, and in five cases GDAs are meanwhile established. We identified 12 new disease causing variants. Two previously reported variants were missed by our updated pipeline due to alignment or reference issues. Our results support the need to re-evaluate screening studies, not only the negative cases but including supposedly solved ones. This also applies in a diagnostic setting. We highlight that the complexity of computational re-analysis for old data should be weighed against the decreasing re-testing costs. Since extensive re-analysis per case is beyond the resources of most institutions, we recommend a screening procedure that would quickly identify the majority (83%) of new variants.

Integrative omics approaches to advance rare disease diagnostics

Over the past decade high-throughput DNA sequencing approaches, namely whole exome and whole genome sequencing became a standard procedure in Mendelian disease diagnostics. Implementation of these technologies greatly facilitated diagnostics and shifted the analysis paradigm from variant identification to prioritisation and evaluation. The diagnostic rates vary widely depending on the cohort size, heterogeneity and disease and range from around 30% to 50% leaving the majority of patients undiagnosed. Advances in omics technologies and computational analysis provide an opportunity to increase these unfavourable rates by providing evidence for disease-causing variant validation and prioritisation. This review aims to provide an overview of the current application of several omics technologies including RNA-sequencing, proteomics, metabolomics and DNA-methylation profiling for diagnostics of rare genetic diseases in general and inborn errors of metabolism in particular.

Meta-analysis of the diagnostic and clinical utility of exome and genome sequencing in pediatric and adult patients with rare diseases across diverse populations

PURPOSE

This meta-analysis aims to compare the diagnostic and clinical utility of exome sequencing (ES) vs genome sequencing (GS) in pediatric and adult patients with rare diseases across diverse populations.

METHODS

A meta-analysis was conducted to identify studies from 2011 to 2021.

RESULTS

One hundred sixty-one studies across 31 countries/regions were eligible, featuring 50,417 probands of diverse populations. Diagnostic rates of ES (0.38, 95% CI 0.36-0.40) and GS (0.34, 95% CI 0.30-0.38) were similar (P = .1). Within-cohort comparison illustrated 1.2-times odds of diagnosis by GS over ES (95% CI 0.79-1.83, P = .38). GS studies discovered a higher range of novel genes than ES studies; yet, the rate of variant of unknown significance did not differ (P = .78). Among high-quality studies, clinical utility of GS (0.77, 95% CI 0.64-0.90) was higher than that of ES (0.44, 95% CI 0.30-0.58) (P < .01).

CONCLUSION

This meta-analysis provides an important update to demonstrate the similar diagnostic rates between ES and GS and the higher clinical utility of GS over ES. With the newly published recommendations for clinical interpretation of variants found in noncoding regions of the genome and the trend of decreasing variant of unknown significance and GS cost, it is expected that GS will be more widely used in clinical settings.

The Need to Standardize the Reanalysis of Genomic Sequencing Results: Findings from Interviews with Underserved Families in Genomic Research

The reanalysis of genomic sequencing results has the potential to provide results that are of considerable medical and personal importance to recipients. Employing interviews with forty-seven predominantly medically underserved families and ethnographic observations we argue that there is pressing need to standardize the approach taken to reanalysis. Our findings highlight that study participants were unclear as to the likelihood of reanalysis happening, the process of initiating reanalysis, and whether they would receive revised results. Their reflections mirror the lack a specific focus upon reanalysis within consent and results sessions as observed in clinical settings. Mechanisms need to be put into place that standardize the approach to reanalysis in research and in clinical contexts. This would enable clinicians and genetic counsellors to communicate clearly with research participants with respect to potential for reanalysis of results and the process of reanalysis. We argue that that the role of reanalysis is too important to be referred to in an ad-hoc manner. Furthermore, the ad-hoc nature of the current process may increase health inequities given the likelihood that only those families who have the means to press for reanalysis are likely to receive it.

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

Broadening the phenotypic and molecular spectrum of FINCA syndrome: Biallelic NHLRC2 variants in 15 novel individuals

FINCA syndrome [MIM: 618278] is an autosomal recessive multisystem disorder characterized by fibrosis, neurodegeneration and cerebral angiomatosis. To date, 13 patients from nine families with biallelic NHLRC2 variants have been published. In all of them, the recurrent missense variant p.(Asp148Tyr) was detected on at least one allele. Common manifestations included lung or muscle fibrosis, respiratory distress, developmental delay, neuromuscular symptoms and seizures often followed by early death due to rapid disease progression.Here, we present 15 individuals from 12 families with an overlapping phenotype associated with nine novel NHLRC2 variants identified by exome analysis. All patients described here presented with moderate to severe global developmental delay and variable disease progression. Seizures, truncal hypotonia and movement disorders were frequently observed. Notably, we also present the first eight cases in which the recurrent p.(Asp148Tyr) variant was not detected in either homozygous or compound heterozygous state.We cloned and expressed all novel and most previously published non-truncating variants in HEK293-cells. From the results of these functional studies, we propose a potential genotype-phenotype correlation, with a greater reduction in protein expression being associated with a more severe phenotype.Taken together, our findings broaden the known phenotypic and molecular spectrum and emphasize that NHLRC2-related disease should be considered in patients presenting with intellectual disability, movement disorders, neuroregression and epilepsy with or without pulmonary involvement.

Whole genome sequencing diagnostic yield for paediatric patients with suspected genetic disorders: systematic review, meta-analysis, and GRADE assessment

BACKGROUND

About 80% of the roughly 7,000 known rare diseases are single gene disorders, about 85% of which are ultra-rare, affecting less than one in one million individuals. NGS technologies, in particular whole genome sequencing (WGS) in paediatric patients suffering from severe disorders of likely genetic origin improve the diagnostic yield allowing targeted, effective care and management. The aim of this study is to perform a systematic review and meta-analysis to assess the effectiveness of WGS, with respect to whole exome sequencing (WES) and/or usual care, for the diagnosis of suspected genetic disorders among the paediatric population.

METHODS

A systematic review of the literature was conducted querying relevant electronic databases, including MEDLINE, EMBASE, ISI Web of Science, and Scopus from January 2010 to June 2022. A random-effect meta-analysis was run to inspect the diagnostic yield of different techniques. A network meta-analysis was also performed to directly assess the comparison between WGS and WES.

RESULTS

Of the 4,927 initially retrieved articles, thirty-nine met the inclusion criteria. Overall results highlighted a significantly higher pooled diagnostic yield for WGS, 38.6% (95% CI: [32.6 - 45.0]), in respect to WES, 37.8% (95% CI: [32.9 - 42.9]) and usual care, 7.8% (95% CI: [4.4 - 13.2]). The meta-regression output suggested a higher diagnostic yield of the WGS compared to WES after controlling for the type of disease (monogenic vs non-monogenic), with a tendency to better diagnostic performances for Mendelian diseases. The network meta-analysis showed a higher diagnostic yield for WGS compared to WES (OR = 1.54, 95%CI: [1.11 - 2.12]).

CONCLUSIONS

Although whole genome sequencing for the paediatric population with suspected genetic disorders provided an accurate and early genetic diagnosis in a high proportion of cases, further research is needed for evaluating costs, effectiveness, and cost-effectiveness of WGS and achieving an informed decision-making process.

TRIAL REGISTRATION

This systematic review has not been registered.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

Neurologic, Neuropsychologic, and Neuroradiologic Features ofEBF3-Related Syndrome

Background and Objectives

Heterozygous mutations or deletions of theEBF3gene are known to cause a syndrome characterized by intellectual disability, neurodevelopmental disorders, facial dysmorphisms, hypotonia, and ataxia; the latter is quite common despite in most patients brain MRI is reported to be normal. Despite the predominant neurologic involvement ofEBF3-related syndrome, a systematic definition of neurologic, cognitive/behavioral, and neuroradiologic features is lacking.

Methods

We report on 6 patients (2 females and 4 males, age range 2–12 years), of whom 4 carrying a heterozygous point mutation of theEBF3gene and 2 with 10q26 deletion encompassing the gene, diagnosed at Carlo Besta Neurologic Institute of Milan, Italy. Clinical evaluation was performed by a pediatric neurologist and pediatric dysmorphologist; ataxia severity was rated by Scale for the Assessment and Rating of Ataxia (SARA); brain MRIs were reviewed by expert neuroradiologists; general quotient levels were obtained through standardized Griffiths Mental Development Scales. Patients carrying a 10q26.3 deletion were diagnosed by array-CGH, whereasEBF3variants were detected by whole exome sequencing.

Results

Phenotype was consistent in all patients, but with wide variability in severity. Developmental milestones were invariably delayed and resulted in an extremely variable cognitive impairment. All patients showed ataxic signs, as confirmed by SARA scores, often associated with hypotonia. Brain MRI revealed in all children a cerebellar malformation with vermis hypoplasia and a peculiar foliation anomaly characterized by a radial disposition of cerebellar folia (dandelion sign). Neurophysiologic examinations were unremarkable.

Discussion

EBF3-related syndrome has been so far described as a neurodevelopmental condition with dysmorphic traits, with limited emphasis on the neurologic features; we highlight the predominant neurologic involvement of these patients, which can be explained at least in part by the underlying cerebellar malformation. We therefore propose thatEBF3-related syndrome should be classified and treated as a congenital, nonprogressive ataxia.

A concurrent dual analysis of genomic data augments diagnoses: Experiences of 2 clinical sites in the Undiagnosed Diseases Network

PURPOSE

Next-generation sequencing (NGS) has revolutionized the diagnostic process for rare/ultrarare conditions. However, diagnosis rates differ between analytical pipelines. In the National Institutes of Health-Undiagnosed Diseases Network (UDN) study, each individual's NGS data are concurrently analyzed by the UDN sequencing core laboratory and the clinical sites. We examined the outcomes of this practice.

METHODS

A retrospective review was performed at 2 UDN clinical sites to compare the variants and diagnoses/candidate genes identified with the dual analyses of the NGS data.

RESULTS

In total, 95 individuals had 100 diagnoses/candidate genes. There was 59% concordance between the UDN sequencing core laboratories and the clinical sites in identifying diagnoses/candidate genes. The core laboratory provided more diagnoses, whereas the clinical sites prioritized more research variants/candidate genes (P < .001). The clinical sites solely identified 15% of the diagnoses/candidate genes. The differences between the 2 pipelines were more often because of variant prioritization disparities than variant detection.

CONCLUSION

The unique dual analysis of NGS data in the UDN synergistically enhances outcomes. The core laboratory provided a clinical analysis with more diagnoses and the clinical sites prioritized more research variants/candidate genes. Implementing such concurrent dual analyses in other genomic research studies and clinical settings can improve both variant detection and prioritization.

Identification of a Splicing Variant c.3813-3A&gt;G in NPHP3 by Reanalysis of Whole Exome Sequencing in a Chinese Boy with Nephronophthisis

Nephronophthisis is an autosomal recessive cystic kidney disease characterized by tubular injury and commonly results in kidney failure. We reported a case of 4-year-old Chinese boy presented with severe anemia, kidney, and liver dysfunction. Whole exome sequencing (WES) was performed to identify the candidate variant with a negative result initially. After complete collection of clinical information, reanalysis of WES identified a homozygous NPHP3 variant c.3813-3A>G (NM_153240.4). The effect on mRNA splicing of the intronic variant was predicted through software (three in silico splice tools). Furthermore, in vitro minigene assay was conducted to validate the predicted deleterious effects of the intronic variant. All of the splice prediction programs and minigene assay indicated that the variant had an impact on the normal splicing pattern of NPHP3. Our study confirmed the effect of the c.3813-3A>G variant on NPHP3 splicing in vitro, which gives additional evidence for the clinical significance of the variant and provides a basis for genetic diagnosis of nephronophthisis 3. In addition, we think that it is essential to reanalyze WES data after the complete clinical information collection to avoid missing some important candidate variants.

A biallelic frameshift indel in PPP1R35 as a cause of primary microcephaly

Protein phosphatase 1 regulatory subunit 35 (PPP1R35) encodes a centrosomal protein required for recruiting microtubule-binding elongation machinery. Several proteins in this centriole biogenesis pathway correspond to established primary microcephaly (MCPH) genes, and multiple model organism studies hypothesize PPP1R35 as a candidate MCPH gene. Here, using exome sequencing (ES) and family-based rare variant analyses, we report a homozygous, frameshifting indel deleting the canonical stop codon in the last exon of PPP1R35 [Chr7: c.753_*3delGGAAGCGTAGACCinsCG (p.Trp251Cysfs*22)]; the variant allele maps in a 3.7 Mb block of absence of heterozygosity (AOH) in a proband with severe MCPH (-4.3 SD at birth, -6.1 SD by 42 months), pachygyria, and global developmental delay from a consanguineous Turkish kindred. Droplet digital PCR (ddPCR) confirmed mutant mRNA expression in fibroblasts. In silico prediction of the translation of mutant PPP1R35 is expected to be elongated by 18 amino acids before encountering a downstream stop codon. This complex indel allele is absent in public databases (ClinVar, gnomAD, ARIC, 1000 genomes) and our in-house database of 14,000+ exomes including 1800+ Turkish exomes supporting predicted pathogenicity. Comprehensive literature searches for PPP1R35 variants yielded two probands affected with severe microcephaly (-15 SD and -12 SD) with the same homozygous indel from a single, consanguineous, Iranian family from a cohort of 404 predominantly Iranian families. The lack of heterozygous cases in two large cohorts representative of the genetic background of these two families decreased our suspicion of a founder allele and supports the contention of a recurrent mutation. We propose two potential secondary structure mutagenesis models for the origin of this variant allele mediated by hairpin formation between complementary GC rich segments flanking the stop codon via secondary structure mutagenesis.

Bridging clinical care and research in Ontario, Canada: Maximizing diagnoses from reanalysis of clinical exome sequencing data

We examined the utility of clinical and research processes in the reanalysis of publicly-funded clinical exome sequencing data in Ontario, Canada. In partnership with eight sites, we recruited 287 families with suspected rare genetic diseases tested between 2014 and 2020. Data from seven laboratories was reanalyzed with the referring clinicians. Reanalysis of clinically relevant genes identified diagnoses in 4% (13/287); four were missed by clinical testing. Translational research methods, including analysis of novel candidate genes, identified candidates in 21% (61/287). Of these, 24 families have additional evidence through data sharing to support likely diagnoses (8% of cohort). This study indicates few diagnoses are missed by clinical laboratories, the incremental gain from reanalysis of clinically-relevant genes is modest, and the highest yield comes from validation of novel disease-gene associations. Future implementation of translational research methods, including continued reporting of compelling genes of uncertain significance by clinical laboratories, should be considered to maximize diagnoses.

The different clinical facets of SYN1-related neurodevelopmental disorders

Synapsin-I (SYN1) is a presynaptic phosphoprotein crucial for synaptogenesis and synaptic plasticity. Pathogenic SYN1 variants are associated with variable X-linked neurodevelopmental disorders mainly affecting males. In this study, we expand on the clinical and molecular spectrum of the SYN1-related neurodevelopmental disorders by describing 31 novel individuals harboring 22 different SYN1 variants. We analyzed newly identified as well as previously reported individuals in order to define the frequency of key features associated with these disorders. Specifically, behavioral disturbances such as autism spectrum disorder or attention deficit hyperactivity disorder are observed in 91% of the individuals, epilepsy in 82%, intellectual disability in 77%, and developmental delay in 70%. Seizure types mainly include tonic-clonic or focal seizures with impaired awareness. The presence of reflex seizures is one of the most representative clinical manifestations related to SYN1. In more than half of the cases, seizures are triggered by contact with water, but other triggers are also frequently reported, including rubbing with a towel, fever, toothbrushing, fingernail clipping, falling asleep, and watching others showering or bathing. We additionally describe hyperpnea, emotion, lighting, using a stroboscope, digestive troubles, and defecation as possible triggers in individuals with SYN1 variants. The molecular spectrum of SYN1 variants is broad and encompasses truncating variants (frameshift, nonsense, splicing and start-loss variants) as well as non-truncating variants (missense substitutions and in-frame duplications). Genotype-phenotype correlation revealed that epileptic phenotypes are enriched in individuals with truncating variants. Furthermore, we could show for the first time that individuals with early seizures onset tend to present with severe-to-profound intellectual disability, hence highlighting the existence of an association between early seizure onset and more severe impairment of cognitive functions. Altogether, we present a detailed clinical description of the largest series of individuals with SYN1 variants reported so far and provide the first genotype-phenotype correlations for this gene. A timely molecular diagnosis and genetic counseling are cardinal for appropriate patient management and treatment.

De novo mutations within metabolism networks of amino acid/protein/energy in Chinese autistic children with intellectual disability

BACKGROUND

Autism spectrum disorder (ASD) is often accompanied by intellectual disability (ID). Despite extensive studies, however, the genetic basis for this comorbidity is still not clear. In this study, we tried to develop an analyzing pipeline for de novo mutations and possible pathways related to ID phenotype in ASD. Whole-exome sequencing (WES) was performed to screen de novo mutations and candidate genes in 79 ASD children together with their parents (trios). The de novo altering genes and relative pathways which were associated with ID phenotype were analyzed. The connection nodes (genes) of above pathways were selected, and the diagnostic value of these selected genes for ID phenotype in the study population was also evaluated.

RESULTS

We identified 89 de novo mutant genes, of which 34 genes were previously reported to be associated with ASD, including double hits in the EGF repeats of NOTCH1 gene (p.V999M and p.S1027L). Interestingly, of these 34 genes, 22 may directly affect intelligence quotient (IQ). Further analyses revealed that these IQ-related genes were enriched in protein synthesis, energy metabolism, and amino acid metabolism, and at least 9 genes (CACNA1A, ALG9, PALM2, MGAT4A, PCK2, PLEKHA1, PSME3, ADI1, and TLE3) were involved in all these three pathways. Seven patients who harbored these gene mutations showed a high prevalence of a low IQ score (< 70), a non-verbal language, and an early diagnostic age (< 4 years). Furthermore, our panel of these 9 genes reached a 10.2% diagnostic rate (5/49) in early diagnostic patients with a low IQ score and also reached a 10% diagnostic yield in those with both a low IQ score and non-verbal language (4/40).

CONCLUSION

We found some new genetic disposition for ASD accompanied with intellectual disability in this study. Our results may be helpful for etiologic research and early diagnoses of intellectual disability in ASD. Larger population studies and further mechanism studies are warranted.

Application of exome sequencing for prenatal diagnosis of fetal structural anomalies: clinical experience and lessons learned from a cohort of 1618 fetuses

Background

Exome sequencing (ES) is becoming more widely available in prenatal diagnosis. However, data on its clinical utility and integration into clinical management remain limited in practice. Herein, we report our experience implementing prenatal ES (pES) in a large cohort of fetuses with anomalies detected by ultrasonography using a hospital-based in-house multidisciplinary team (MDT) facilitated by a three-step genotype-driven followed by phenotype-driven analysis framework.

Methods

We performed pES in 1618 fetal cases with positive ultrasound findings but negative for karyotyping and chromosome microarray analysis between January 2014 and October 2021, including both retrospective (n=565) and prospective (n=1053) cohorts. The diagnostic efficiency and its correlation to organ systems involved, phenotypic spectrum, and the clinical impacts of pES results on pregnancy outcomes were analyzed.

Results

A genotype-driven followed by phenotype-driven three-step approach was carried out in all trio pES. Step 1, a genotype-driven analysis resulted in a diagnostic rate of 11.6% (187/1618). Step 2, a phenotype-driven comprehensive analysis yielded additional diagnostic findings for another 28 cases (1.7%; 28/1618). In the final step 3, data reanalyses based on new phenotypes and/or clinical requests found molecular diagnosis in 14 additional cases (0.9%; 14/1618). Altogether, 229 fetal cases (14.2%) received a molecular diagnosis, with a higher positive rate in the retrospective than the prospective cohort (17.3% vs. 12.4%, p<0.01). The diagnostic rates were highest in fetuses with skeletal anomalies (30.4%) and multiple organ involvements (25.9%), and lowest in fetuses with chest anomalies (0%). In addition, incidental and secondary findings with childhood-onset disorders were detected in 11 (0.7%) cases. Furthermore, we described the prenatal phenotypes for the first time for 27 gene-associated conditions (20.0%, 27/135) upon a systematic analysis of the diagnosed cases and expanded the phenotype spectrum for 26 (19.3%) genes where limited fetal phenotypic information was available. In the prospective cohort, the combined prenatal ultrasound and pES results had significantly impacted the clinical decisions (61.5%, 648/1053).

Conclusions

The genotype-driven approach could identify about 81.7% positive cases (11.6% of the total cohort) with the initial limited fetal phenotype information considered. The following two steps of phenotype-driven analysis and data reanalyses helped us find the causative variants in an additional 2.6% of the entire cohort (18.3% of all positive findings). Our extensive phenotype analysis on a large number of molecularly confirmed prenatal cases had greatly enriched our current knowledge on fetal phenotype-genotype correlation, which may guide more focused prenatal ultrasound in the future. This is by far the largest pES cohort study that combines a robust trio sequence data analysis, systematic phenotype-genotype correlation, and well-established MDT in a single prenatal clinical setting. This work underlines the value of pES as an essential component in prenatal diagnosis in guiding medical management and parental decision making.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01130-x.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation

BACKGROUND

The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases.

METHODS

Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype.

RESULTS

In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations.

CONCLUSIONS

Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

X-CAP improves pathogenicity prediction of stopgain variants

Stopgain substitutions are the third-largest class of monogenic human disease mutations and often examined first in patient exomes. Existing computational stopgain pathogenicity predictors, however, exhibit poor performance at the high sensitivity required for clinical use. Here, we introduce a new classifier, termed X-CAP, which uses a novel training methodology and unique feature set to improve the AUROC by 18% and decrease the false-positive rate 4-fold on large variant databases. In patient exomes, X-CAP prioritizes causal stopgains better than existing methods do, further illustrating its clinical utility. X-CAP is available at https://github.com/bejerano-lab/X-CAP .

Clinically available testing options resulting in diagnosis in post-exome clinic at one medical center

Exome sequencing (ES) became clinically available in 2011 and promised an agnostic, unbiased next-generation sequencing (NGS) platform for patients with symptoms believed to have a genetic etiology. The diagnostic yield of ES has been estimated to be between 25–40% and may be higher in specific clinical scenarios. Those who remain undiagnosed may have no molecular findings of interest on ES, variants of uncertain significance in genes that are linked to human disease, or variants of uncertain significance in candidate genes that are not definitively tied to human disease. Recent evidence suggests that a post-exome evaluation consisting of clinical re-phenotyping, functional studies of candidate variants in known genes, and variant reevaluation can lead to a diagnosis in 5–15% of additional cases. In this brief research study, we present our experience on post-exome evaluations in a cohort of patients who are believed to have a genetic etiology for their symptoms. We have reached a full or partial diagnosis in approximately 18% (6/33) of cases that have completed evaluations to date. We accomplished this by utilizing NGS-based methods that are available on a clinical basis. A sample of these cases highlights the utility of ES reanalysis with updated phenotyping allowing for the discovery of new genes, re-adjudication of known variants, incorporating updated phenotypic information, utilizing functional testing such as targeted RNA sequencing, and deploying other NGS-based testing methods such as gene panels and genome sequencing to reach a diagnosis.

An Integrated Phenotypic and Genotypic Approach Reveals a High-Risk Subtype Association for EBF3 Missense Variants Affecting the Zinc Finger Domain

OBJECTIVE

Collier/Olf/EBF (COE) transcription factors have distinct expression patterns in the developing and mature nervous system. To date, a neurological disease association has been conclusively established for only the Early B-cell Factor-3 (EBF3) COE family member through the identification of heterozygous loss-of-function variants in individuals with autism spectrum/neurodevelopmental disorders (NDD). Here, we identify a symptom severity risk association with missense variants primarily disrupting the zinc finger domain (ZNF) in EBF3-related NDD.

METHODS

A phenotypic assessment of 41 individuals was combined with a literature meta-analysis for a total of 83 individuals diagnosed with EBF3-related NDD. Quantitative diagnostic phenotypic and symptom severity scales were developed to compare EBF3 variant type and location to identify genotype-phenotype correlations. To stratify the effects of EBF3 variants disrupting either the DNA-binding domain (DBD) or the ZNF, we used in vivo fruit fly UAS-GAL4 expression and in vitro luciferase assays.

RESULTS

We show that patient symptom severity correlates with EBF3 missense variants perturbing the ZNF, which is a key protein domain required for stabilizing the interaction between EBF3 and the target DNA sequence. We found that ZNF-associated variants failed to restore viability in the fruit fly and impaired transcriptional activation. However, the recurrent variant EBF3 p.Arg209Trp in the DBD is capable of partially rescuing viability in the fly and preserved transcriptional activation.

INTERPRETATION

We describe a symptom severity risk association with ZNF perturbations and EBF3 loss-of-function in the largest reported cohort to date of EBF3-related NDD patients. This analysis should have potential predictive clinical value for newly identified patients with EBF3 gene variants. ANN NEUROL 2022;92:138-153.

De novo heterozygous variants in SLC30A7 are a candidate cause for Joubert syndrome

Joubert syndrome (JS), a well-established ciliopathy, is characterized by the distinctive molar tooth sign on brain MRI, ataxia, and neurodevelopmental features. Other manifestations can include polydactyly, accessory frenula, renal, or liver disease. Here, we report individuals meeting criteria for JS with de novo heterozygous variants in SLC30A7 (Chr1p21.2). The first individual is a female with history of unilateral postaxial polydactyly, classic molar tooth sign on MRI, macrocephaly, ataxia, ocular motor apraxia, neurodevelopmental delay, and precocious puberty. Exome sequencing detected a de novo heterozygous missense variant in SLC30A7: NM_133496.5: c.407 T > C, (p.Val136Ala). The second individual had bilateral postaxial polydactyly, molar tooth sign, macrocephaly, developmental delay, and an extra oral frenulum. A de novo deletion-insertion variant in SLC30A7, c.490_491delinsAG (p.His164Ser) was found. Both de novo variants affect highly conserved residues. Variants were not identified in known Joubert genes for either case. SLC30A7 has not yet been associated with a human phenotype. The SLC30 family of zinc transporters, like SLC30A7, permit cellular efflux of zinc, and although it is expressed in the brain its functions remain unknown. Published data from proteomic studies support SLC30A7 interaction with TCTN3, another protein associated with JS. The potential involvement of such genes in primary cilia suggest a role in Sonic Hedgehog signaling. SLC30A7 is a candidate JS-associated gene. Future work could be directed toward further characterization of SLC30A7 variants and understanding its function.

Reanalysis of exome negative patients with rare disease: a pragmatic workflow for diagnostic applications

Background

Approximately two third of patients with a rare genetic disease remain undiagnosed after exome sequencing (ES). As part of our post-test counseling procedures, patients without a conclusive diagnosis are advised to recontact their referring clinician to discuss new diagnostic opportunities in due time. We performed a systematic study of genetically undiagnosed patients 5 years after their initial negative ES report to determine the efficiency of diverse reanalysis strategies.

Methods

We revisited a cohort of 150 pediatric neurology patients originally enrolled at Radboud University Medical Center, of whom 103 initially remained genetically undiagnosed. We monitored uptake of physician-initiated routine clinical and/or genetic re-evaluation (ad hoc re-evaluation) and performed systematic reanalysis, including ES-based resequencing, of all genetically undiagnosed patients (systematic re-evaluation).

Results

Ad hoc re-evaluation was initiated for 45 of 103 patients and yielded 18 diagnoses (including 1 non-genetic). Subsequent systematic re-evaluation identified another 14 diagnoses, increasing the diagnostic yield in our cohort from 31% (47/150) to 53% (79/150). New genetic diagnoses were established by reclassification of previously identified variants (10%, 3/31), reanalysis with enhanced bioinformatic pipelines (19%, 6/31), improved coverage after resequencing (29%, 9/31), and new disease-gene associations (42%, 13/31). Crucially, our systematic study also showed that 11 of the 14 further conclusive genetic diagnoses were made in patients without a genetic diagnosis that did not recontact their referring clinician.

Conclusions

We find that upon re-evaluation of undiagnosed patients, both reanalysis of existing ES data as well as resequencing strategies are needed to identify additional genetic diagnoses. Importantly, not all patients are routinely re-evaluated in clinical care, prolonging their diagnostic trajectory, unless systematic reanalysis is facilitated. We have translated our observations into considerations for systematic and ad hoc reanalysis in routine genetic care.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01069-z.

Biallelic loss of EMC10 leads to mild to severe intellectual disability

The endoplasmic reticulum membrane protein complex subunit 10 (EMC10) is a highly conserved protein responsible for the post‐translational insertion of tail‐anchored membrane proteins into the endoplasmic reticulum in a defined topology. Two biallelic variants in EMC10 have previously been associated with a neurodevelopmental disorder. Utilizing exome sequencing and international data sharing we have identified 10 affected individuals from six independent families with five new biallelic loss‐of‐function and one previously reported recurrent EMC10 variants. This report expands the molecular and clinical spectrum of EMC10 deficiency, provides a comprehensive dysmorphological assessment and highlights an overlap between the clinical features of EMC10‐and EMC1‐related disease.

Recommendations for next generation sequencing data reanalysis of unsolved cases with suspected Mendelian disorders: A systematic review and meta-analysis

PURPOSE

The study aimed to determine the diagnostic yield, optimal timing, and methodology of next generation sequencing data reanalysis in suspected Mendelian disorders.

METHODS

We conducted a systematic review and meta-analysis of studies that conducted data reanalysis in patients with suspected Mendelian disorders. Random effects model was used to pool the estimated outcome with subgroup analysis stratified by timing, sequencing methodology, sample size, segregation, use of research validation, and artificial intelligence (AI) variant curation tools.

RESULTS

A search of PubMed, Embase, Scopus, and Web of Science between 2007 and 2021 yielded 9327 articles, of which 29 were selected. Significant heterogeneity was noted between studies. Reanalysis had an overall diagnostic yield of 0.10 (95% CI = 0.06-0.13). Literature updates accounted for most new diagnoses. Diagnostic yield was higher after 24 months, although this was not statistically significant. Increased diagnoses were obtained with research validation and data sharing. AI-based tools did not adversely affect reanalysis diagnostic rate.

CONCLUSION

Next generation sequencing data reanalysis can improve diagnostic yield. Owing to the heterogeneity of the studies, the optimal time to reanalysis and the impact of AI-based tools could not be determined with confidence. We propose standardized guidelines for future studies to reduce heterogeneity and improve the quality of the conclusions.

Systematic Collaborative Reanalysis of Genomic Data Improves Diagnostic Yield in Neurologic Rare Diseases

Many patients experiencing a rare disease remain undiagnosed even after genomic testing. Reanalysis of existing genomic data has shown to increase diagnostic yield, although there are few systematic and comprehensive reanalysis efforts that enable collaborative interpretation and future reinterpretation. The Undiagnosed Rare Disease Program of Catalonia project collated previously inconclusive good quality genomic data (panels, exomes, and genomes) and standardized phenotypic profiles from 323 families (543 individuals) with a neurologic rare disease. The data were reanalyzed systematically to identify relatedness, runs of homozygosity, consanguinity, single-nucleotide variants, insertions and deletions, and copy number variants. Data were shared and collaboratively interpreted within the consortium through a customized Genome-Phenome Analysis Platform, which also enables future data reinterpretation. Reanalysis of existing genomic data provided a diagnosis for 20.7% of the patients, including 1.8% diagnosed after the generation of additional genomic data to identify a second pathogenic heterozygous variant. Diagnostic rate was significantly higher for family-based exome/genome reanalysis compared with singleton panels. Most new diagnoses were attributable to recent gene-disease associations (50.8%), additional or improved bioinformatic analysis (19.7%), and standardized phenotyping data integrated within the Undiagnosed Rare Disease Program of Catalonia Genome-Phenome Analysis Platform functionalities (18%).

Detection and impact of genetic disease in a level IV neonatal intensive care unit

OBJECTIVE

To determine detection rates of genetic disease in a level IV neonatal intensive care unit (NICU) and cost of care.

STUDY DESIGN

We divided 2703 neonates, admitted between 2013 and 2016 to a level IV NICU, into two epochs and determined how genetic testing utilization, genetic diagnoses identified, and cost of NICU care changed over time.

RESULT

The increasing use of multi-gene panels 104 vs 184 (P = 0.02) and whole exome sequencing (WES) 9 vs 28 (P = 0.03) improved detection of genetic disease, 9% vs 12% (P < 0.01). Individuals with genetic diagnoses had higher mean NICU charges, $723,422 vs $417,013 (P < 0.01) secondary to longer lengths of stay, not genetic services.

CONCLUSION

The increased utilization of broad genetic testing improved the detection of genetic disease but contributed minimally to the cost of care while bolstering understanding of the patient's condition and prognosis.

Pathophysiological Heterogeneity of the BBSOA Neurodevelopmental Syndrome

The formation and maturation of the human brain is regulated by highly coordinated developmental events, such as neural cell proliferation, migration and differentiation. Any impairment of these interconnected multi-factorial processes can affect brain structure and function and lead to distinctive neurodevelopmental disorders. Here, we review the pathophysiology of the Bosch–Boonstra–Schaaf Optic Atrophy Syndrome (BBSOAS; OMIM 615722; ORPHA 401777), a recently described monogenic neurodevelopmental syndrome caused by the haploinsufficiency of NR2F1 gene, a key transcriptional regulator of brain development. Although intellectual disability, developmental delay and visual impairment are arguably the most common symptoms affecting BBSOAS patients, multiple additional features are often reported, including epilepsy, autistic traits and hypotonia. The presence of specific symptoms and their variable level of severity might depend on still poorly characterized genotype–phenotype correlations. We begin with an overview of the several mutations of NR2F1 identified to date, then further focuses on the main pathological features of BBSOAS patients, providing evidence—whenever possible—for the existing genotype–phenotype correlations. On the clinical side, we lay out an up-to-date list of clinical examinations and therapeutic interventions recommended for children with BBSOAS. On the experimental side, we describe state-of-the-art in vivo and in vitro studies aiming at deciphering the role of mouse Nr2f1, in physiological conditions and in pathological contexts, underlying the BBSOAS features. Furthermore, by modeling distinct NR2F1 genetic alterations in terms of dimer formation and nuclear receptor binding efficiencies, we attempt to estimate the total amounts of functional NR2F1 acting in developing brain cells in normal and pathological conditions. Finally, using the NR2F1 gene and BBSOAS as a paradigm of monogenic rare neurodevelopmental disorder, we aim to set the path for future explorations of causative links between impaired brain development and the appearance of symptoms in human neurological syndromes.

Re-analysis of genomic data: An overview of the mechanisms and complexities of clinical adoption

Re-analyzing genomic information from a patient suspected of having an underlying genetic condition can improve the diagnostic yield of sequencing tests, potentially providing significant benefits to the patient and to the health care system. Although a significant number of studies have shown the clinical potential of re-analysis, less work has been performed to characterize the mechanisms responsible for driving the increases in diagnostic yield. Complexities surrounding re-analysis have also emerged. The terminology itself represents a challenge because "re-analysis" can refer to a range of different concepts. Other challenges include the increased workload that re-analysis demands of curators, adequate reimbursement pathways for clinical and diagnostic services, and the development of systems to handle large volumes of data. Re-analysis also raises ethical implications for patients and families, most notably when re-classification of a variant alters diagnosis, treatment, and prognosis. This review highlights the possibilities and complexities associated with the re-analysis of existing clinical genomic data. We propose a terminology that builds on the foundation presented in a recent statement from the American College of Medical Genetics and Genomics and describes each re-analysis process. We identify mechanisms for increasing diagnostic yield and provide perspectives on the range of challenges that must be addressed by health care systems and individual patients.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the non-canonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable AR trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome. This article is protected by copyright. All rights reserved.

The Genetic Diagnosis of Ultrarare DEEs: An Ongoing Challenge

Epileptic encephalopathies (EEs) and developmental and epileptic encephalopathies (DEEs) are a group of severe early-onset neurodevelopmental disorders (NDDs). In recent years, next-generation equencing (NGS) technologies enabled the discovery of numerous genes involved in these conditions. However, more than 50% of patients remained undiagnosed. A major obstacle lies in the high degree of genetic heterogeneity and the wide phenotypic variability that has characterized these disorders. Interpreting a large amount of NGS data is also a crucial challenge. This study describes a dynamic diagnostic procedure used to investigate 17 patients with DEE or EE with previous negative or inconclusive genetic testing by whole-exome sequencing (WES), leading to a definite diagnosis in about 59% of participants. Biallelic mutations caused most of the diagnosed cases (50%), and a pathogenic somatic mutation resulted in 10% of the subjects. The high diagnostic yield reached highlights the relevance of the scientific approach, the importance of the reverse phenotyping strategy, and the involvement of a dedicated multidisciplinary team. The study emphasizes the role of recessive and somatic variants, new genetic mechanisms, and the complexity of genotype–phenotype associations. In older patients, WES results could end invasive diagnostic procedures and allow a more accurate transition. Finally, an early pursued diagnosis is essential for comprehensive care of patients, precision approach, knowledge of prognosis, patient and family planning, and quality of life.

Quantitative dissection of multilocus pathogenic variation in an Egyptian infant with severe neurodevelopmental disorder resulting from multiple molecular diagnoses

Genomic sequencing and clinical genomics have demonstrated that substantial subsets of atypical and/or severe disease presentations result from multilocus pathogenic variation (MPV) causing blended phenotypes. In an infant with a severe neurodevelopmental disorder, four distinct molecular diagnoses were found by exome sequencing (ES). The blended phenotype that includes brain malformation, dysmorphism, and hypotonia was dissected using the Human Phenotype Ontology (HPO). ES revealed variants in CAPN3 (c.259C > G:p.L87V), MUSK (c.1781C > T:p.A594V), NAV2 (c.1996G > A:p.G666R), and ZC4H2 (c.595A > C:p.N199H). CAPN3, MUSK, and ZC4H2 are established disease genes linked to limb-girdle muscular dystrophy (OMIM# 253600), congenital myasthenia (OMIM# 616325), and Wieacker-Wolff syndrome (WWS; OMIM# 314580), respectively. NAV2 is a retinoic-acid responsive novel disease gene candidate with biological roles in neurite outgrowth and cerebellar dysgenesis in mouse models. Using semantic similarity, we show that no gene identified by ES individually explains the proband phenotype, but rather the totality of the clinically observed disease is explained by the combination of disease-contributing effects of the identified genes. These data reveal that multilocus pathogenic variation can result in a blended phenotype with each gene affecting a different part of the nervous system and nervous system-muscle connection. We provide evidence from this n = 1 study that in patients with MPV and complex blended phenotypes resulting from multiple molecular diagnoses, quantitative HPO analysis can allow for dissection of phenotypic contribution of both established disease genes and novel disease gene candidates not yet proven to cause human disease.

A guide for the diagnosis of rare and undiagnosed disease: beyond the exome

Rare diseases affect 30 million people in the USA and more than 300–400 million worldwide, often causing chronic illness, disability, and premature death. Traditional diagnostic techniques rely heavily on heuristic approaches, coupling clinical experience from prior rare disease presentations with the medical literature. A large number of rare disease patients remain undiagnosed for years and many even die without an accurate diagnosis. In recent years, gene panels, microarrays, and exome sequencing have helped to identify the molecular cause of such rare and undiagnosed diseases. These technologies have allowed diagnoses for a sizable proportion (25–35%) of undiagnosed patients, often with actionable findings. However, a large proportion of these patients remain undiagnosed. In this review, we focus on technologies that can be adopted if exome sequencing is unrevealing. We discuss the benefits of sequencing the whole genome and the additional benefit that may be offered by long-read technology, pan-genome reference, transcriptomics, metabolomics, proteomics, and methyl profiling. We highlight computational methods to help identify regionally distant patients with similar phenotypes or similar genetic mutations. Finally, we describe approaches to automate and accelerate genomic analysis. The strategies discussed here are intended to serve as a guide for clinicians and researchers in the next steps when encountering patients with non-diagnostic exomes.

An evaluation of pipelines for DNA variant detection can guide a reanalysis protocol to increase the diagnostic ratio of genetic diseases

Clinical exome (CE) sequencing has become a first-tier diagnostic test for hereditary diseases; however, its diagnostic rate is around 30-50%. In this study, we aimed to increase the diagnostic yield of CE using a custom reanalysis algorithm. Sequencing data were available for three cohorts using two commercial protocols applied as part of the diagnostic process. Using these cohorts, we compared the performance of general and clinically relevant variant calling and the efficacy of an in-house bioinformatic protocol (FJD-pipeline) in detecting causal variants as compared to commercial protocols. On the whole, the FJD-pipeline detected 99.74% of the causal variants identified by the commercial protocol in previously solved cases. In the unsolved cases, FJD-pipeline detects more INDELs and non-exonic variants, and is able to increase the diagnostic yield in 2.5% and 3.2% in the re-analysis of 78 cancer and 62 cardiovascular cases. These results were considered to design a reanalysis, filtering and prioritization algorithm that was tested by reassessing 68 inconclusive cases of monoallelic autosomal recessive retinal dystrophies increasing the diagnosis by 4.4%. In conclusion, a guided NGS reanalysis of unsolved cases increases the diagnostic yield in genetic disorders, making it a useful diagnostic tool in medical genetics.

Novel pathogenic variants and quantitative phenotypic analyses of Robinow syndrome: WNT signaling perturbation and phenotypic variability

Robinow syndrome (RS) is a genetically heterogeneous disorder with six genes that converge on the WNT/planar cell polarity (PCP) signaling pathway implicated (DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A). RS is characterized by skeletal dysplasia and distinctive facial and physical characteristics. To further explore the genetic heterogeneity, paralog contribution, and phenotypic variability of RS, we investigated a cohort of 22 individuals clinically diagnosed with RS from 18 unrelated families. Pathogenic or likely pathogenic variants in genes associated with RS or RS phenocopies were identified in all 22 individuals, including the first variant to be reported in DVL2. We retrospectively collected medical records of 16 individuals from this cohort and extracted clinical descriptions from 52 previously published cases. We performed Human Phenotype Ontology (HPO) based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. Individuals with FZD2 variants clustered into two groups with demonstrable phenotypic differences between those with missense and truncating alleles. Probands with biallelic NXN variants clustered together with the majority of probands carrying DVL1, DVL2, and DVL3 variants, demonstrating no phenotypic distinction between the NXN-autosomal recessive and dominant forms of RS. While phenotypically similar diseases on the RS differential matched through HPO analysis, clustering using phenotype similarity score placed RS-associated phenotypes in a unique cluster containing WNT5A, FZD2, and ROR2 apart from non-RS-associated paralogs. Through human phenotype analyses of this RS cohort and OMIM clinical synopses of Mendelian disease, this study begins to tease apart specific biologic roles for non-canonical WNT-pathway proteins.

OUP accepted manuscript

The mutations in modifier genes may contribute to some inherited diseases including Wilson disease (WD). This study was designed to identify potential modifier genes that contribute to WD. A total of 10 WD patients with single or no heterozygous ATP7B mutations were recruited for whole-exome sequencing (WES). Five hundred and thirteen candidate genes, of which the genetic variants present in at least two patients, were identified. In order to clarify which proteins might be involved in copper transfer or metabolism processes, the isobaric tags for relative and absolute quantitation (iTRAQ) was performed to identify the differentially expressed proteins between normal and CuSO₄-treated cell lines. Thirteen genes/proteins were identified by both WES and iTRAQ, indicating that disease-causing variants of these genes may actually contribute to the aberrant copper ion accumulation. Additionally, the c.86C > T (p.S29L) mutation in the SLC31A2 gene (coding CTR2) has a relative higher frequency in our cohort of WD patients (6/191) than reported (0.0024 in gnomAD database) in our healthy donors (0/109), and CTR2^S29L leads to increased intracellular Cu concentration and Cu-induced apoptosis in cultured cell lines. In conclusion, the WES and iTRAQ approaches successfully identified several disease-causing variants in potential modifier genes that may be involved in the WD phenotype.

Hyper-IgE Syndrome due to an Elusive Novel Intronic Homozygous Variant in DOCK8

Rare, biallelic loss-of-function mutations in DOCK8 result in a combined immune deficiency characterized by severe and recurrent cutaneous infections, eczema, allergies, and susceptibility to malignancy, as well as impaired humoral and cellular immunity and hyper-IgE. The advent of next-generation sequencing technologies has enabled the rapid molecular diagnosis of rare monogenic diseases, including inborn errors of immunity. These advances have resulted in the implementation of gene-guided treatments, such as hematopoietic stem cell transplant for DOCK8 deficiency. However, putative disease-causing variants revealed by next-generation sequencing need rigorous validation to demonstrate pathogenicity. Here, we report the eventual diagnosis of DOCK8 deficiency in a consanguineous family due to a novel homozygous intronic deletion variant that caused aberrant exon splicing and subsequent loss of expression of DOCK8 protein. Remarkably, the causative variant was not initially detected by clinical whole-genome sequencing but was subsequently identified and validated by combining advanced genomic analysis, RNA-seq, and flow cytometry. This case highlights the need to adopt multipronged confirmatory approaches to definitively solve complex genetic cases that result from variants outside protein-coding exons and conventional splice sites.

Functional Studies of Genetic Variants Associated with Human Diseases in Notch Signaling-Related Genes Using Drosophila

Rare variants in the many genes related to Notch signaling cause diverse Mendelian diseases that affect myriad organ systems. In addition, genome- and exome-wide association studies have linked common and rare variants in Notch-related genes to common diseases and phenotypic traits. Moreover, somatic mutations in these genes have been observed in many types of cancer, some of which are classified as oncogenic and others as tumor suppressive. While functional characterization of some of these variants has been performed through experimental studies, the number of "variants of unknown significance" identified in patients with diverse conditions keeps increasing as high-throughput sequencing technologies become more commonly used in the clinic. Furthermore, as disease gene discovery efforts identify rare variants in human genes that have yet to be linked to a disease, the demand for functional characterization of variants in these "genes of unknown significance" continues to increase. In this chapter, we describe a workflow to functionally characterize a rare variant in a Notch signaling related gene that was found to be associated with late-onset Alzheimer's disease. This pipeline involves informatic analysis of the variant of interest using diverse human and model organism databases, followed by in vivo experiments in the fruit fly Drosophila melanogaster. The protocol described here can be used to study variants that affect amino acids that are not conserved between human and fly. By "humanizing" the almondex gene in Drosophila with mutant alleles and heterologous genomic rescue constructs, a missense variant in TM2D3 (TM2 Domain Containing 3) was shown to be functionally damaging. This, and similar approaches, greatly facilitate functional interpretations of genetic variants in the human genome and propel personalized medicine.

A mutation in the neonatal isoform of SCN2A causes neonatal-onset epilepsy

SCN2A (sodium channel 2A) encodes the Nav1.2 channel protein in excitatory neurons in the brain. Nav1.2 is a critical voltage-gated sodium channel of the central nervous system. Mutations in SCN2A are responsible for a broad phenotypic spectrum ranging from autism and developmental delay to severe encephalopathy with neonatal or early infantile onset. SCN2A can be spliced into two different isoforms, a neonatal (6N) and an adult (6A) form. After birth, there is an equal or higher amount of the 6N isoform, protecting the brain from the increased neuronal excitability of the infantile brain. During postnatal development, 6N is gradually replaced by 6A. In an infant carrying the novel SCN2A mutation c.643G > A (p.Ala215Thr) only in the neonatal transcript, seizures started immediately after birth. The clinical presentation evolved from a burst-suppression pattern with 30-50 tonic seizures per day to hypsarrhythmia. The first exome analysis, focusing only on common transcripts, missed the diagnosis and delayed early therapy. A reevaluation including all transcripts revealed the SCN2A variant.

Genetic Variation and the Distribution of Variant Types in the Horse

Genetic variation is a key contributor to health and disease. Understanding the link between an individual's genotype and the corresponding phenotype is a major goal of medical genetics. Whole genome sequencing (WGS) within and across populations enables highly efficient variant discovery and elucidation of the molecular nature of virtually all genetic variation. Here, we report the largest catalog of genetic variation for the horse, a species of importance as a model for human athletic and performance related traits, using WGS of 534 horses. We show the extent of agreement between two commonly used variant callers. In data from ten target breeds that represent major breed clusters in the domestic horse, we demonstrate the distribution of variants, their allele frequencies across breeds, and identify variants that are unique to a single breed. We investigate variants with no homozygotes that may be potential embryonic lethal variants, as well as variants present in all individuals that likely represent regions of the genome with errors, poor annotation or where the reference genome carries a variant. Finally, we show regions of the genome that have higher or lower levels of genetic variation compared to the genome average. This catalog can be used for variant prioritization for important equine diseases and traits, and to provide key information about regions of the genome where the assembly and/or annotation need to be improved.