Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

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.

Intellectual Functioning of Children With Isolated PRS, PRS-Plus, and Syndromic PRS.

OBJECTIVE: Describe the intelligence quotient (IQ) of children with Pierre Robin sequence (PRS). DESIGN: Prospective cohort study. SETTING: Neurodevelopmental follow-up clinic within a hospital. PATIENTS: Children with PRS (n = 45) who had been in the Neonatal Intensive Care Unit (NICU) were classified by a geneticist into 3 subgroups of isolated PRS (n = 20), PRS-plus additional medical features (n = 8), and syndromic PRS (n = 17) based on medical record review and genetic testing. MAIN OUTCOME MEASURE: Children with PRS completed IQ testing at 5 or 8 years of age with the Wechsler Preschool and Primary Scale of Intelligence, Third Edition (WPPSI-III) or Fourth Edition (WPPSI-IV) or the Wechsler Intelligence Scale for Children, Fourth Edition (WISC-IV) or Fifth Edition (WISC-V). RESULTS: IQ scores were more than 1 to 2 standard deviations below the mean for 36% of the overall sample, which was significantly greater compared to test norms (binomial test P = .001). There was a significant association between PRS subtype and IQ (Fisher's exact P = .026). While only 20% of children with isolated PRS were within 1 standard deviation below average and 35% of children with syndromic PRS were below 1 to 2 standard deviations, 75% of PRS-plus children scored lower than 1 to 2 standard deviations below the mean. CONCLUSION: PRS subgroups can help identify children at risk for cognitive delay. The majority of children with PRS-plus had low intellectual functioning, in contrast to the third of children with syndromic PRS who had low IQ and the majority of children with isolated PRS who had average or higher IQ.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance.

PURPOSE: Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis. METHODS: Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins. RESULTS: Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis. CONCLUSION: This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

DNA methylation episignatures are sensitive and specific biomarkers for detection of patients with KAT6A/KAT6B variants.

Accurate diagnosis for patients living with neurodevelopmental disorders is often met with numerous challenges, related to the ambiguity of findings and lack of specificity in genetic variants leading to pathology. Genome-wide DNA methylation analysis has been used to develop highly sensitive and specific 'episignatures' as biomarkers capable of differentiating and classifying complex neurodevelopmental disorders. In this study we describe distinct episignatures for KAT6A syndrome, caused by pathogenic variants in the lysine acetyltransferase A gene (KAT6A), and for the two neurodevelopmental disorders associated with lysine acetyl transferase B (KAT6B). We demonstrate the ability of our models to differentiate between highly overlapping episignatures, increasing the ability to effectively identify and diagnose these conditions.

MOGS-CDG: Quantitative analysis of the diagnostic Glc3 Man tetrasaccharide and clinical spectrum of six new cases.

Congenital disorders of glycosylation (CDG) are a clinically and biochemically heterogeneous subgroup of inherited metabolic disorders. Most CDG with abnormal N-glycosylation can be detected by transferrin screening, however, MOGS-CDG escapes this routine screening. Combined with the clinical heterogeneity of reported cases, diagnosing MOGS-CDG can be challenging. Here, we clinically characterize ten MOGS-CDG cases including six previously unreported individuals, showing a phenotype characterized by dysmorphic features, global developmental delay, muscular hypotonia, and seizures in all patients and in a minority vision problems and hypogammaglobulinemia. Glycomics confirmed accumulation of a Glc3 Man7 GlcNAc2 glycan in plasma. For quantification of the diagnostic Glcα1-3Glcα1-3Glcα1-2Man tetrasaccharide in urine, we developed and validated a liquid chromatography-mass spectrometry method of 2-aminobenzoic acid (2AA) labeled urinary glycans. As an internal standard, isotopically labeled 13 C6 -2AA Glc3 Man was used, while labeling efficiency was controlled by use of 12 C6 -2AA and 13 C6 -2AA labeled laminaritetraose. Recovery, linearity, intra- and interassay coefficients of variability of these labeled compounds were determined. Furthermore, Glc3 Man was specifically identified by retention time matching against authentic MOGS-CDG urine and compared with Pompe urine. Glc3 Man was increased in all six analyzed cases, ranging from 34.1 to 618.0 µmol/mmol creatinine (reference <5 µmol). In short, MOGS-CDG has a broad manifestation of symptoms but can be diagnosed with the use of a quantitative method for analysis of urinary Glc3 Man excretion.

Comparing saliva and blood for the detection of mosaic genomic abnormalities that cause syndromic intellectual disability.

We aimed to determine whether SNP-microarray genomic testing of saliva had a greater diagnostic yield than blood for pathogenic copy number variants (CNVs). We selected patients who underwent CMA testing of both blood and saliva from 23,289 blood and 21,857 saliva samples. Our cohort comprised 370 individuals who had testing of both, 224 with syndromic intellectual disability (ID) and 146 with isolated ID. Mosaic pathogenic CNVs or aneuploidy were detected in saliva but not in blood in 20/370 (4.4%). All 20 individuals had syndromic ID, accounting for 9.1% of the syndromic ID sub-cohort. Pathogenic CNVs were large in size (median of 46 Mb), and terminal in nature, with median mosaicism of 27.5% (not exceeding 40%). By contrast, non-mosaic pathogenic CNVs were 100% concordant between blood and saliva, considerably smaller in size (median of 0.65 Mb), and predominantly interstitial in location. Given that salivary microarray testing has increased diagnostic utility over blood in individuals with syndromic ID, we recommend it as a first-tier testing in this group.

Unique Capabilities of Genome Sequencing for Rare Disease Diagnosis.

Background: Causal variants underlying rare disorders may remain elusive even after expansive gene panels or exome sequencing (ES). Clinicians and researchers may then turn to genome sequencing (GS), though the added value of this technique and its optimal use remain poorly defined. We therefore investigated the advantages of GS within a phenotypically diverse cohort. Methods: GS was performed for 744 individuals with rare disease who were genetically undiagnosed. Analysis included review of single nucleotide, indel, structural, and mitochondrial variants. Results: We successfully solved 218/744 (29.3%) cases using GS, with most solves involving established disease genes (157/218, 72.0%). Of all solved cases, 148 (67.9%) had previously had non-diagnostic ES. We systematically evaluated the 218 causal variants for features requiring GS to identify and 61/218 (28.0%) met these criteria, representing 8.2% of the entire cohort. These included small structural variants (13), copy neutral inversions and complex rearrangements (8), tandem repeat expansions (6), deep intronic variants (15), and coding variants that may be more easily found using GS related to uniformity of coverage (19). Conclusion: We describe the diagnostic yield of GS in a large and diverse cohort, illustrating several types of pathogenic variation eluding ES or other techniques. Our results reveal a higher diagnostic yield of GS, supporting the utility of a genome-first approach, with consideration of GS as a secondary or tertiary test when higher-resolution structural variant analysis is needed or there is a strong clinical suspicion for a condition and prior targeted genetic testing has been negative.

The diverse pleiotropic effects of spliceosomal protein PUF60: A case series of Verheij syndrome.

Verheij syndrome (VRJS) is a rare craniofacial spliceosomopathy presenting with craniofacial dysmorphism, multiple congenital anomalies and variable neurodevelopmental delay. It is caused by single nucleotide variants (SNVs) in PUF60 or interstitial deletions of the 8q24.3 region. PUF60 encodes a splicing factor which forms part of the spliceosome. To date, 36 patients with a sole diagnosis of VRJS due to disease-causing PUF60 SNVs have been reported in peer-reviewed publications. Although the depth of their phenotyping has varied greatly, they exhibit marked phenotypic heterogeneity. We report 10 additional unrelated patients, including the first described patients of Khmer, Indian, and Vietnamese ethnicities, and the eldest patient to date, with 10 heterozygous PUF60 variants identified through exome sequencing, 8 previously unreported. All patients underwent deep phenotyping identifying variable dysmorphism, growth delay, neurodevelopmental delay, and multiple congenital anomalies, including several unique features. The eldest patient is the only reported individual with a germline variant and neither neurodevelopmental delay nor intellectual disability. In combining these detailed phenotypic data with that of previously reported patients (n = 46), we further refine the known frequencies of features associated with VRJS. These include neurodevelopmental delay/intellectual disability (98%), axial skeletal anomalies (74%), appendicular skeletal anomalies (73%), oral anomalies (68%), short stature (66%), cardiac anomalies (63%), brain malformations (48%), hearing loss (46%), microcephaly (41%), colobomata (38%), and other ocular anomalies (65%). This case series, incorporating three patients from previously unreported ethnic backgrounds, further delineates the broad pleiotropy and mutational spectrum of PUF60 pathogenic variants.

Variant-specific changes in RAC3 function disrupt corticogenesis in neurodevelopmental phenotypes.

Variants in RAC3, encoding a small GTPase RAC3 which is critical for the regulation of actin cytoskeleton and intracellular signal transduction, are associated with a rare neurodevelopmental disorder with structural brain anomalies and facial dysmorphism. We investigated a cohort of 10 unrelated participants presenting with global psychomotor delay, hypotonia, behavioural disturbances, stereotyped movements, dysmorphic features, seizures and musculoskeletal abnormalities. MRI of brain revealed a complex pattern of variable brain malformations, including callosal abnormalities, white matter thinning, grey matter heterotopia, polymicrogyria/dysgyria, brainstem anomalies and cerebellar dysplasia. These patients harboured eight distinct de novo RAC3 variants, including six novel variants (NM_005052.3): c.34G > C p.G12R, c.179G > A p.G60D, c.186_188delGGA p.E62del, c.187G > A p.D63N, c.191A > G p.Y64C and c.348G > C p.K116N. We then examined the pathophysiological significance of these novel and previously reported pathogenic variants p.P29L, p.P34R, p.A59G, p.Q61L and p.E62K. In vitro analyses revealed that all tested RAC3 variants were biochemically and biologically active to variable extent, and exhibited a spectrum of different affinities to downstream effectors including p21-activated kinase 1. We then focused on the four variants p.Q61L, p.E62del, p.D63N and p.Y64C in the Switch II region, which is essential for the biochemical activity of small GTPases and also a variation hot spot common to other Rho family genes, RAC1 and CDC42. Acute expression of the four variants in embryonic mouse brain using in utero electroporation caused defects in cortical neuron morphology and migration ending up with cluster formation during corticogenesis. Notably, defective migration by p.E62del, p.D63N and p.Y64C were rescued by a dominant negative version of p21-activated kinase 1. Our results indicate that RAC3 variants result in morphological and functional defects in cortical neurons during brain development through variant-specific mechanisms, eventually leading to heterogeneous neurodevelopmental phenotypes.

Inherited variants in CHD3 show variable expressivity in Snijders Blok-Campeau syndrome.

PURPOSE: Common diagnostic next-generation sequencing strategies are not optimized to identify inherited variants in genes associated with dominant neurodevelopmental disorders as causal when the transmitting parent is clinically unaffected, leaving a significant number of cases with neurodevelopmental disorders undiagnosed. METHODS: We characterized 21 families with inherited heterozygous missense or protein-truncating variants in CHD3, a gene in which de novo variants cause Snijders Blok-Campeau syndrome. RESULTS: Computational facial and Human Phenotype Ontology-based comparisons showed that the phenotype of probands with inherited CHD3 variants overlaps with the phenotype previously associated with de novo CHD3 variants, whereas heterozygote parents are mildly or not affected, suggesting variable expressivity. In addition, similarly reduced expression levels of CHD3 protein in cells of an affected proband and of healthy family members with a CHD3 protein-truncating variant suggested that compensation of expression from the wild-type allele is unlikely to be an underlying mechanism. Notably, most inherited CHD3 variants were maternally transmitted. CONCLUSION: Our results point to a significant role of inherited variation in Snijders Blok-Campeau syndrome, a finding that is critical for correct variant interpretation and genetic counseling and warrants further investigation toward understanding the broader contributions of such variation to the landscape of human disease.

Lessons learnt from multifaceted diagnostic approaches to the first 150 families in Victoria's Undiagnosed Diseases Program.

BACKGROUND: Clinical exome sequencing typically achieves diagnostic yields of 30%-57.5% in individuals with monogenic rare diseases. Undiagnosed diseases programmes implement strategies to improve diagnostic outcomes for these individuals. AIM: We share the lessons learnt from the first 3 years of the Undiagnosed Diseases Program-Victoria, an Australian programme embedded within a clinical genetics service in the state of Victoria with a focus on paediatric rare diseases. METHODS: We enrolled families who remained without a diagnosis after clinical genomic (panel, exome or genome) sequencing between 2016 and 2018. We used family-based exome sequencing (family ES), family-based genome sequencing (family GS), RNA sequencing (RNA-seq) and high-resolution chromosomal microarray (CMA) with research-based analysis. RESULTS: In 150 families, we achieved a diagnosis or strong candidate in 64 (42.7%) (37 in known genes with a consistent phenotype, 3 in known genes with a novel phenotype and 24 in novel disease genes). Fifty-four diagnoses or strong candidates were made by family ES, six by family GS with RNA-seq, two by high-resolution CMA and two by data reanalysis. CONCLUSION: We share our lessons learnt from the programme. Flexible implementation of multiple strategies allowed for scalability and response to the availability of new technologies. Broad implementation of family ES with research-based analysis showed promising yields post a negative clinical singleton ES. RNA-seq offered multiple benefits in family ES-negative populations. International data sharing strategies were critical in facilitating collaborations to establish novel disease-gene associations. Finally, the integrated approach of a multiskilled, multidisciplinary team was fundamental to having diverse perspectives and strategic decision-making.

Pathogenic variants in SMARCA5, a chromatin remodeler, cause a range of syndromic neurodevelopmental features.

Intellectual disability encompasses a wide spectrum of neurodevelopmental disorders, with many linked genetic loci. However, the underlying molecular mechanism for more than 50% of the patients remains elusive. We describe pathogenic variants in SMARCA5, encoding the ATPase motor of the ISWI chromatin remodeler, as a cause of a previously unidentified neurodevelopmental disorder, identifying 12 individuals with de novo or dominantly segregating rare heterozygous variants. Accompanying phenotypes include mild developmental delay, frequent postnatal short stature and microcephaly, and recurrent dysmorphic features. Loss of function of the SMARCA5 Drosophila ortholog Iswi led to smaller body size, reduced sensory dendrite complexity, and tiling defects in larvae. In adult flies, Iswi neural knockdown caused decreased brain size, aberrant mushroom body morphology, and abnormal locomotor function. Iswi loss of function was rescued by wild-type but not mutant SMARCA5. Our results demonstrate that SMARCA5 pathogenic variants cause a neurodevelopmental syndrome with mild facial dysmorphia.

A DNA repair disorder caused by de novo monoallelic DDB1 variants is associated with a neurodevelopmental syndrome.

The DNA damage-binding protein 1 (DDB1) is part of the CUL4-DDB1 ubiquitin E3 ligase complex (CRL4), which is essential for DNA repair, chromatin remodeling, DNA replication, and signal transduction. Loss-of-function variants in genes encoding the complex components CUL4 and PHIP have been reported to cause syndromic intellectual disability with hypotonia and obesity, but no phenotype has been reported in association with DDB1 variants. Here, we report eight unrelated individuals, identified through Matchmaker Exchange, with de novo monoallelic variants in DDB1, including one recurrent variant in four individuals. The affected individuals have a consistent phenotype of hypotonia, mild to moderate intellectual disability, and similar facies, including horizontal or slightly bowed eyebrows, deep-set eyes, full cheeks, a short nose, and large, fleshy and forward-facing earlobes, demonstrated in the composite face generated from the cohort. Digital anomalies, including brachydactyly and syndactyly, were common. Three older individuals have obesity. We show that cells derived from affected individuals have altered DDB1 function resulting in abnormal DNA damage signatures and histone methylation following UV-induced DNA damage. Overall, our study adds to the growing family of neurodevelopmental phenotypes mediated by disruption of the CRL4 ubiquitin ligase pathway and begins to delineate the phenotypic and molecular effects of DDB1 misregulation.

The chromatin remodeler ISWI acts during Drosophila development to regulate adult sleep.

Sleep disruptions are among the most commonly reported symptoms across neurodevelopmental disorders (NDDs), but mechanisms linking brain development to normal sleep are largely unknown. From a Drosophila screen of human NDD-associated risk genes, we identified the chromatin remodeler Imitation SWItch/SNF (ISWI) to be required for adult fly sleep. Loss of ISWI also results in disrupted circadian rhythms, memory, and social behavior, but ISWI acts in different cells and during distinct developmental times to affect each of these adult behaviors. Specifically, ISWI expression in type I neuroblasts is required for both adult sleep and formation of a learning-associated brain region. Expression in flies of the human ISWI homologs SMARCA1 and SMARCA5 differentially rescues adult phenotypes, while de novo SMARCA5 patient variants fail to rescue sleep. We propose that sleep deficits are a primary phenotype of early developmental origin in NDDs and point toward chromatin remodeling machinery as critical for sleep circuit formation.

Cognitive, academic, and behavioral functioning in school-aged children born with esophageal atresia.

PURPOSE: To characterize cognitive, academic, and behavioral functioning in children who underwent neonatal surgical repair of esophageal atresia (OA) and compare outcomes according to clinical characteristics (presence of additional congenital anomalies, longer hospitalization, and prematurity). METHODS: Intellectual, language, attention, and executive functioning were assessed in 71 5-year-olds and 72 8-year-olds born with OA. At 8 years, memory and academic skills were also assessed. Parents rated children's executive functioning and behavior via questionnaires. Outcomes were compared to normative data and within subgroups of the sample. RESULTS: Intellectual functioning varied depending on the assessment tool, with some evidence of lower than expected intellectual development in children with OA. At 5 years, children with OA showed age-appropriate language and self-regulation, but reduced verbal attention. At 8 years, the OA group had lower than expected sustained attention, divided attention, and mathematics but typical memory and literacy. Parents consistently reported increased working memory difficulties. Other executive functioning and behavioral symptoms were transiently observed. Findings did not consistently differ according to clinical characteristics. CONCLUSIONS: Children with OA may be at risk of transient and persisting cognitive difficulties, particularly in attention and working memory. Difficulties were not strongly associated with additional congenital anomalies, longer hospitalization, or prematurity. LEVEL OF EVIDENCE: Level IV.

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

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

Galactose treatment of a PGM1 patient presenting with restrictive cardiomyopathy.

We report a patient diagnosed with PGM1-CDG at 11 years of age after two biallelic likely pathogenic variants in PGM1 were found on research genomic sequencing. To our knowledge, he is the first patient with PGM1-CDG to be reported with a restrictive cardiomyopathy. Other clinical manifestations included cleft palate, asymptomatic elevated transaminases, intellectual disability and myopathy resulting in exercise intolerance. He was trialed on oral galactose therapy in increasing doses for 18 weeks to assess if there was any biochemical and clinical benefit. His galactose was continued for a further 9 months beyond the initial galactose treatment period due to improvements in exercise tolerance and myopathy. Treatment with galactose demonstrated an improvement in liver function and myopathy with improved exercise tolerance. Treatment with galactose for 15 months did not change heart function and exercise stress test results were stable.