Genome Sequencing for Diagnosing Rare Diseases.

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

Implementation of Rapid Genome Sequencing for Critically Ill Infants With Complex Congenital Heart Disease.

BACKGROUND: Rapid genome sequencing (rGS) has been shown to improve care of critically ill infants. Congenital heart disease (CHD) is a leading cause of infant mortality and is often caused by genetic disorders, yet the utility of rGS has not been prospectively studied in this population. METHODS: We conducted a prospective evaluation of rGS to improve the care of infants with complex CHD in our cardiac neonatal intensive care unit. RESULTS: In a cohort of 48 infants with complex CHD, rGS diagnosed 14 genetic disorders in 13 (27%) individuals and led to changes in clinical management in 8 (62%) cases with diagnostic results. These included 2 cases in whom genetic diagnoses helped avert intensive, futile interventions before cardiac neonatal intensive care unit discharge, and 3 cases in whom eye disease was diagnosed and treated in early childhood. CONCLUSIONS: Our study provides the first prospective evaluation of rGS for infants with complex CHD to our knowledge. We found that rGS diagnosed genetic disorders in 27% of cases and led to changes in management in 62% of cases with diagnostic results. Our model of care depended on coordination between neonatologists, cardiologists, surgeons, geneticists, and genetic counselors. These findings highlight the important role of rGS in CHD and demonstrate the need for expanded study of how to implement this resource to a broader population of infants with CHD.

Identification of copy number variants with genome sequencing: Clinical experiences from the NYCKidSeq program.

Copy number variations (CNVs) play a significant role in human disease. While chromosomal microarray has traditionally been the first-tier test for CNV detection, use of genome sequencing (GS) is increasing. We report the frequency of CNVs detected with GS in a diverse pediatric cohort from the NYCKidSeq program and highlight specific examples of its clinical impact. A total of 1052 children (0-21 years) with neurodevelopmental, cardiac, and/or immunodeficiency phenotypes received GS. Phenotype-driven analysis was used, resulting in 183 (17.4%) participants with a diagnostic result. CNVs accounted for 20.2% of participants with a diagnostic result (37/183) and ranged from 0.5 kb to 16 Mb. Of participants with a diagnostic result (n = 183) and phenotypes in more than one category, 5/17 (29.4%) were solved by a CNV finding, suggesting a high prevalence of diagnostic CNVs in participants with complex phenotypes. Thirteen participants with a diagnostic CNV (35.1%) had previously uninformative genetic testing, of which nine included a chromosomal microarray. This study demonstrates the benefits of GS for reliable detection of CNVs in a pediatric cohort with variable phenotypes.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

PURPOSE: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. METHODS: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. RESULTS: A total of 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses (P < .001). Yield was greater for GS vs TGPs in Hispanic/Latino(a) (17.2% vs 9.5%, P < .001) and White/European American (19.8% vs 7.9%, P < .001) but not in Black/African American (11.5% vs 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. CONCLUSION: GS may yield up to twice as many diagnoses in pediatric patients compared with TGP testing but not yet across all population groups.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

Purpose: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. Methods: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. Results: 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses ( P < .001). Yield was greater for GS vs . TGPs in Hispanic/Latino(a) (17.2% vs . 9.5%, P < .001) and White/European American (19.8% vs . 7.9%, P < .001), but not in Black/African American (11.5% vs . 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs . White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. Conclusion: GS may yield up to twice as many diagnoses in pediatric patients compared to TGP testing, but not yet across all population groups.

Detection of mosaic variants using genome sequencing in a large pediatric cohort.

The increased use of next-generation sequencing has expanded our understanding of the involvement and prevalence of mosaicism in genetic disorders. We describe a total of eleven cases: nine in which mosaic variants detected by genome sequencing (GS) and/or targeted gene panels (TGPs) were considered to be causative for the proband's phenotype, and two of apparent parental mosaicism. Variants were identified in the following genes: PHACTR1, SCN8A, KCNT1, CDKL5, NEXMIF, CUX1, TSC2, GABRB2, and SMARCB1. In addition, we identified one large duplication including three genes, UBE3A, GABRB3, and MAGEL2, and one large deletion including deletion of ARFGAP1, EEF1A2, CHRNA4, and KCNQ2. All patients were enrolled in the NYCKidSeq study, a research program studying the communication of genomic information in clinical care, as well as the clinical utility and diagnostic yield of GS for children with suspected genetic disorders in diverse populations in New York City. We observed variability in the correlation between reported variant allele fraction and the severity of the patient's phenotype, although we were not able to determine the mosaicism percentage in clinically relevant tissue(s). Although our study was not sufficiently powered to assess differences in mosaicism detection between the two testing modalities, we saw a trend toward better detection by GS as compared with TGP testing. This case series supports the importance of mosaicism in childhood-onset genetic conditions and informs guidelines for laboratory and clinical interpretation of mosaic variants detected by GS.

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.

Predictive functional, statistical and structural analysis of CSNK2A1 and CSNK2B variants linked to neurodevelopmental diseases.

Okur-Chung Neurodevelopmental Syndrome (OCNDS) and Poirier-Bienvenu Neurodevelopmental Syndrome (POBINDS) were recently identified as rare neurodevelopmental disorders. OCNDS and POBINDS are associated with heterozygous mutations in the CSNK2A1 and CSNK2B genes which encode CK2α, a serine/threonine protein kinase, and CK2ß, a regulatory protein, respectively, which together can form a tetrameric enzyme called protein kinase CK2. A challenge in OCNDS and POBINDS is to understand the genetic basis of these diseases and the effect of the various CK2⍺ and CK2ß mutations. In this study we have collected all variants available to date in CSNK2A1 and CSNK2B, and identified hotspots. We have investigated CK2⍺ and CK2ß missense mutations through prediction programs which consider the evolutionary conservation, functionality and structure or these two proteins, compared these results with published experimental data on CK2α and CK2ß mutants, and suggested prediction programs that could help predict changes in functionality of CK2α mutants. We also investigated the potential effect of CK2α and CK2ß mutations on the 3D structure of the proteins and in their binding to each other. These results indicate that there are functional and structural consequences of mutation of CK2α and CK2ß, and provide a rationale for further study of OCNDS and POBINDS-associated mutations. These data contribute to understanding the genetic and functional basis of these diseases, which is needed to identify their underlying mechanisms.

A homozygous splice variant in ATP5PO, disrupts mitochondrial complex V function and causes Leigh syndrome in two unrelated families.

Mitochondrial complex V plays an important role in oxidative phosphorylation by catalyzing the generation of ATP. Most complex V subunits are nuclear encoded and not yet associated with recognized Mendelian disorders. Using exome sequencing, we identified a rare homozygous splice variant (c.87+3A>G) in ATP5PO, the complex V subunit which encodes the oligomycin sensitivity conferring protein, in three individuals from two unrelated families, with clinical suspicion of a mitochondrial disorder. These individuals had a similar, severe infantile and often lethal multi-systemic disorder that included hypotonia, developmental delay, hypertrophic cardiomyopathy, progressive epileptic encephalopathy, progressive cerebral atrophy, and white matter abnormalities on brain MRI consistent with Leigh syndrome. cDNA studies showed a predominant shortened transcript with skipping of exon 2 and low levels of the normal full-length transcript. Fibroblasts from the affected individuals demonstrated decreased ATP5PO protein, defective assembly of complex V with markedly reduced amounts of peripheral stalk proteins, and complex V hydrolytic activity. Further, expression of human ATP5PO cDNA without exon 2 (hATP5PO-∆ex2) in yeast cells deleted for yATP5 (ATP5PO homolog) was unable to rescue growth on media which requires oxidative phosphorylation when compared to the wild type construct (hATP5PO-WT), indicating that exon 2 deletion leads to a non-functional protein. Collectively, our findings support the pathogenicity of the ATP5PO c.87+3A>G variant, which significantly reduces but does not eliminate complex V activity. These data along with the recent report of an affected individual with ATP5PO variants, add to the evidence that rare biallelic variants in ATP5PO result in defective complex V assembly, function and are associated with Leigh syndrome.

De novo variants in H3-3A and H3-3B are associated with neurodevelopmental delay, dysmorphic features, and structural brain abnormalities.

The histone H3 variant H3.3, encoded by two genes H3-3A and H3-3B, can replace canonical isoforms H3.1 and H3.2. H3.3 is important in chromatin compaction, early embryonic development, and lineage commitment. The role of H3.3 in somatic cancers has been studied extensively, but its association with a congenital disorder has emerged just recently. Here we report eleven de novo missense variants and one de novo stop-loss variant in H3-3A (n = 6) and H3-3B (n = 6) from Baylor Genetics exome cohort (n = 11) and Matchmaker Exchange (n = 1), of which detailed phenotyping was conducted for 10 individuals (H3-3A = 4 and H3-3B = 6) that showed major phenotypes including global developmental delay, short stature, failure to thrive, dysmorphic facial features, structural brain abnormalities, hypotonia, and visual impairment. Three variant constructs (p.R129H, p.M121I, and p.I52N) showed significant decrease in protein expression, while one variant (p.R41C) accumulated at greater levels than wild-type control. One H3.3 variant construct (p.R129H) was found to have stronger interaction with the chaperone death domain-associated protein 6.

Neuropathological Findings in a Case of Parkinsonism and Developmental Delay Associated with a Monoallelic Variant in PLXNA1.

BACKGROUND: PLXNA1 encodes for Plexin-A, a transmembrane protein expressed in the developing nervous system. Mutations in this gene have been associated with developmental delay but have not been previously associated with the development of parkinsonism. OBJECTIVES: To describe the case of a 38-year-old patient with developmental delay who developed parkinsonism later in life. METHODS: Post-mortem exome sequencing was performed with confirmation by Sanger sequencing. Brain autopsy was also performed. RESULTS: Post-mortem exome sequencing on the proband identified a heterozygous predicted nonsense PLXNA1 variant (c.G3361T:p.Glu1121Ter). Pathology demonstrated arhinencephaly with brainstem heterotopia, diffuse Lewy body disease, and frontotemporal lobar dementia-tau. CONCLUSIONS: This case of a patient with developmental delay and parkinsonism with PLXNA1 mutation highlights a need for assessing long-term outcomes of individuals with neurodevelopmental disorders, as well as the need for genetic testing in adults. It also suggests that the link between PLXNA1 and α-synuclein should be explored in the future. © 2021 International Parkinson and Movement Disorder Society.

Clinical and genomic characterization of 8p cytogenomic disorders.

PURPOSE: To provide a detailed clinical and cytogenomic summary of individuals with chromosome 8p rearrangements of invdupdel(8p), del(8p), and dup(8p). METHODS: We enrolled 97 individuals with invdupdel(8p), del(8p), and dup(8p). Clinical and molecular data were collected to delineate and compare the clinical findings and rearrangement breakpoints. We included additional 5 individuals with dup(8p) from the literature for a total of 102 individuals. RESULTS: Eighty-one individuals had recurrent rearrangements of invdupdel(8p) (n = 49), del(8p)_distal (n = 4), del(8p)_proximal (n = 9), del(8p)_proximal&distal (n = 12), and dup(8p)_proximal (n = 7). Twenty-one individuals had nonrecurrent rearrangements. While all individuals had neurodevelopmental features, the frequency and severity of clinical findings were higher in individuals with invdupdel(8p), and with larger duplications. All individuals with GATA4 deletion had structural congenital heart defects; however, the presence of structural heart defects in some individuals with normal GATA4 copy number suggests there are other potentially contributing gene(s) on 8p. CONCLUSION: Our study may inform families and health-care providers about the associated clinical findings and severity in individuals with chromosome 8p rearrangements, and guide researchers in investigating the underlying molecular and biological mechanisms by providing detailed clinical and cytogenomic information about individuals with distinct 8p rearrangements.

Biallelic and monoallelic variants in PLXNA1 are implicated in a novel neurodevelopmental disorder with variable cerebral and eye anomalies.

PURPOSE: To investigate the effect of PLXNA1 variants on the phenotype of patients with autosomal dominant and recessive inheritance patterns and to functionally characterize the zebrafish homologs plxna1a and plxna1b during development. METHODS: We assembled ten patients from seven families with biallelic or de novo PLXNA1 variants. We describe genotype-phenotype correlations, investigated the variants by structural modeling, and used Morpholino knockdown experiments in zebrafish to characterize the embryonic role of plxna1a and plxna1b. RESULTS: Shared phenotypic features among patients include global developmental delay (9/10), brain anomalies (6/10), and eye anomalies (7/10). Notably, seizures were predominantly reported in patients with monoallelic variants. Structural modeling of missense variants in PLXNA1 suggests distortion in the native protein. Our zebrafish studies enforce an embryonic role of plxna1a and plxna1b in the development of the central nervous system and the eye. CONCLUSION: We propose that different biallelic and monoallelic variants in PLXNA1 result in a novel neurodevelopmental syndrome mainly comprising developmental delay, brain, and eye anomalies. We hypothesize that biallelic variants in the extracellular Plexin-A1 domains lead to impaired dimerization or lack of receptor molecules, whereas monoallelic variants in the intracellular Plexin-A1 domains might impair downstream signaling through a dominant-negative effect.

A novel homozygous variant in TRAPPC2L results in a neurodevelopmental disorder and disrupts TRAPP complex function.

BACKGROUND: Next-generation sequencing has facilitated the diagnosis of neurodevelopmental disorders with variable and non-specific clinical findings. Recently, a homozygous missense p.(Asp37Tyr) variant in TRAPPC2L, a core subunit of TRAPP complexes which function as tethering factors during membrane trafficking, was reported in two unrelated individuals with neurodevelopmental delay, post-infectious encephalopathy-associated developmental arrest, tetraplegia and accompanying rhabdomyolysis. METHODS: We performed whole genome sequencing on members of an Ashkenazi Jewish pedigree to identify the underlying genetic aetiology of global developmental delay/intellectual disability in three affected siblings. To assess the effect of the identified TRAPPC2L variant, we performed biochemical and cell biological functional studies on the TRAPPC2L protein. RESULTS: A rare homozygous predicted deleterious missense variant, p.(Ala2Gly), in TRAPPC2L was identified in the affected siblings and it segregated with the neurodevelopmental phenotype within the family. Using a yeast two-hybrid assay and in vitro binding, we demonstrate that the p.(Ala2Gly) variant, but not the p.(Asp37Tyr) variant, disrupted the interaction between TRAPPC2L and another core TRAPP protein, TRAPPC6a. Size exclusion chromatography suggested that this variant affects the assembly of TRAPP complexes. Employing two different membrane trafficking assays using fibroblasts from one of the affected siblings, we found a delay in traffic into and out of the Golgi. Similar to the p.(Asp37Tyr) variant, the p.(Ala2Gly) variant resulted in an increase in the levels of active RAB11. CONCLUSION: Our data fill in a gap in the knowledge of TRAPP architecture with TRAPPC2L interacting with TRAPPC6a, positioning it as a putative adaptor for other TRAPP subunits. Collectively, our findings support the pathogenicity of the TRAPPC2L p.(Ala2Gly) variant.

Early-Onset Parkinsonism Is a Manifestation of the PPP2R5D p.E200K Mutation.

PPP2R5D-related neurodevelopmental disorder is characterized by a range of neurodevelopmental and behavioral manifestations. We report the association of early-onset parkinsonism with the PPP2R5D p.E200K mutation. Clinical characterization and exome sequencing were performed on three patients, with postmortem neuropathologic examination for one patient. All patients had mild developmental delay and developed levodopa-responsive parkinsonism between the ages of 25 and 40 years. The PPP2R5D c.598G>A (p.E200K) mutation was identified in all patients. Neuropathologic examination demonstrated uneven, focally severe neuronal loss and gliosis in the substantia nigra pars compacta, without Lewy bodies. Our findings suggest the PPP2R5D p.E200K mutation to be a possible new cause of early-onset parkinsonism. ANN NEUROL 2020;88:1028-1033.

De novo heterozygous missense and loss-of-function variants in CDC42BPB are associated with a neurodevelopmental phenotype.

CDC42BPB encodes MRCKß (myotonic dystrophy-related Cdc42-binding kinase beta), a serine/threonine protein kinase, and a downstream effector of CDC42, which has recently been associated with Takenouchi-Kosaki syndrome, an autosomal dominant neurodevelopmental disorder. We identified 12 heterozygous predicted deleterious variants in CDC42BPB (9 missense, 2 frameshift, and 1 nonsense) in 14 unrelated individuals (confirmed de novo in 11/14) with neurodevelopmental disorders including developmental delay/intellectual disability, autism, hypotonia, and structural brain abnormalities including cerebellar vermis hypoplasia and agenesis/hypoplasia of the corpus callosum. The frameshift and nonsense variants in CDC42BPB are expected to be gene-disrupting and lead to haploinsufficiency via nonsense-mediated decay. All missense variants are located in highly conserved and functionally important protein domains/regions: 3 are found in the protein kinase domain, 2 are in the citron homology domain, and 4 in a 20-amino acid sequence between 2 coiled-coil regions, 2 of which are recurrent. Future studies will help to delineate the natural history and to elucidate the underlying biological mechanisms of the missense variants leading to the neurodevelopmental and behavioral phenotypes.

Biallelic variants in AGMO with diminished enzyme activity are associated with a neurodevelopmental disorder.

Alkylglycerol monooxygenase (AGMO) is the only enzyme known to cleave the O-alkyl bonds of ether lipids (alkylglycerols) which are essential components of cell membranes. A homozygous frameshift variant [p.(Glu324LysfsTer12)] in AGMO has recently been reported in two male siblings with syndromic microcephaly. In this study, we identified rare nonsense, in frame deletion, and missense biallelic variants in AGMO in two unrelated individuals with neurodevelopmental disabilities. We assessed the activity of seven disease associated AGMO variants including the four variants identified in our two affected individuals expressed in human embryonic kidney (HEK293T) cells. We demonstrated significantly diminished enzyme activity for all disease-associated variants, supporting the mechanism as decreased AGMO activity. Future mechanistic studies are necessary to understand how decreased AGMO activity leads to the neurologic manifestations.

Clinical and genetic characterization of individuals with predicted deleterious PHIP variants.

Heterozygous deleterious variants in PHIP have been associated with behavioral problems, intellectual disability/developmental delay, obesity/overweight, and dysmorphic features (BIDOD syndrome). We report an additional 10 individuals with pleckstrin homology domain-interacting protein (PHIP)-predicted deleterious variants (four frameshift, three missense, two nonsense, and one splice site; six of which are confirmed de novo). The mutation spectrum is diverse, and there is no clustering of mutations across the protein. The clinical phenotype of these individuals is consistent with previous reports and includes behavioral problems, intellectual disability, developmental delay, hypotonia, and dysmorphic features. The additional individuals we report have a lower frequency of obesity than previous reports and a higher frequency of gastrointestinal problems, social deficits, and behavioral challenges. Characterizing additional individuals with diverse mutations longitudinally will provide better natural history data to assist with medical management and educational and behavioral support.

Homozygous noncanonical splice variant in LSM1 in two siblings with multiple congenital anomalies and global developmental delay.

Two siblings, one male and one female, ages 6 and 13 yr old, have similar clinical features of global developmental delay, multiple congenital anomalies affecting the cardiac, genitourinary, and skeletal systems, and abnormal eye movements. Whole-genome sequencing revealed a homozygous splice variant (NM_014462.3:c.231+4A>C) in LSM1 that segregated with the phenotype in the family. LSM1 has a role in pre-mRNA splicing and degradation. Expression studies revealed absence of expression of the canonical isoform in the affected individuals. The Lsm1 knockout mice have a partially overlapping phenotype that affects the brain, heart, and eye. To our knowledge, LSM1 has not been associated with any human disorder; however, the tissue expression pattern, gene constraint, and the similarity of the phenotype in our patients and the knockout mice models suggest it has a role in the development of multiple organ systems in humans.