The landscape of regional missense mutational intolerance quantified from 125,748 exomes

Missense variants can have a range of functional impacts depending on factors such as the specific amino acid substitution and location within the gene. To interpret their deleteriousness, studies have sought to identify regions within genes that are specifically intolerant of missense variation 1-12 . Here, we leverage the patterns of rare missense variation in 125,748 individuals in the Genome Aggregation Database (gnomAD) 13 against a null mutational model to identify transcripts that display regional differences in missense constraint. Missense-depleted regions are enriched for ClinVar 14 pathogenic variants, de novo missense variants from individuals with neurodevelopmental disorders (NDDs) 15,16 , and complex trait heritability. Following ClinGen calibration recommendations for the ACMG/AMP guidelines, we establish that regions with less than 20% of their expected missense variation achieve moderate support for pathogenicity. We create a missense deleteriousness metric (MPC) that incorporates regional constraint and outperforms other deleteriousness scores at stratifying case and control de novo missense variation, with a strong enrichment in NDDs. These results provide additional tools to aid in missense variant interpretation.

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.

Differential inclusion of NEB exons 143 and 144 provides insight into NEB-related myopathy variant interpretation and disease manifestation

Biallelic pathogenic variants in the gene encoding nebulin (NEB) are a known cause of congenital myopathy. We present two individuals with congenital myopathy and compound heterozygous variants (NM_001271208.2: c.2079C>A; p.(Cys693Ter) and c.21522+3A>G ) in NEB. Transcriptomic sequencing on patient muscle revealed that the extended splice variant c.21522+3A>G causes exon 144 skipping. Nebulin isoforms containing exon 144 are known to be mutually exclusive with isoforms containing exon 143, and these isoforms are differentially expressed during development and in adult skeletal muscles. Patients MRIs were compared to the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these patients better fits the distribution of exon 144-containing isoforms in muscle than with previously published MRI findings in NEB-related disease due to other variants. To our knowledge this is the first report hypothesizing disease pathogenesis through the alteration of isoform distributions in muscle.

A harmonized public resource of deeply sequenced diverse human genomes

Underrepresented populations are often excluded from genomic studies due in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high quality set of 4,094 whole genomes from HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also demonstrate substantial added value from this dataset compared to the prior versions of the component resources, typically combined via liftover and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared to previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

Biallelic CRELD1 variants cause a multisystem syndrome, including neurodevelopmental phenotypes, cardiac dysrhythmias, and frequent infections

PURPOSE

We sought to delineate a multisystem disorder caused by recessive cysteine-rich with epidermal growth factor-like domains 1 (CRELD1) gene variants.

METHODS

The impact of CRELD1 variants was characterized through an international collaboration utilizing next-generation DNA sequencing, gene knockdown, and protein overexpression in Xenopus tropicalis, and in vitro analysis of patient immune cells.

RESULTS

Biallelic variants in CRELD1 were found in 18 participants from 14 families. Affected individuals displayed an array of phenotypes involving developmental delay, early-onset epilepsy, and hypotonia, with about half demonstrating cardiac arrhythmias and some experiencing recurrent infections. Most harbored a frameshift in trans with a missense allele, with 1 recurrent variant, p.(Cys192Tyr), identified in 10 families. X tropicalis tadpoles with creld1 knockdown displayed developmental defects along with increased susceptibility to induced seizures compared with controls. Additionally, human CRELD1 harboring missense variants from affected individuals had reduced protein function, indicated by a diminished ability to induce craniofacial defects when overexpressed in X tropicalis. Finally, baseline analyses of peripheral blood mononuclear cells showed similar proportions of immune cell subtypes in patients compared with healthy donors.

CONCLUSION

This patient cohort, combined with experimental data, provide evidence of a multisystem clinical syndrome mediated by recessive variants in CRELD1.

Recurrent de novo SPTLC2 variant causes childhood-onset amyotrophic lateral sclerosis (ALS) by excess sphingolipid synthesis

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease of the upper and lower motor neurons with varying ages of onset, progression and pathomechanisms. Monogenic childhood-onset ALS, although rare, forms an important subgroup of ALS. We recently reported specific SPTLC1 variants resulting in sphingolipid overproduction as a cause for juvenile ALS. Here, we report six patients from six independent families with a recurrent, de novo, heterozygous variant in SPTLC2 c.778G>A [p.Glu260Lys] manifesting with juvenile ALS.

METHODS

Clinical examination of the patients along with ancillary and genetic testing, followed by biochemical investigation of patients' blood and fibroblasts, was performed.

RESULTS

All patients presented with early-childhood-onset progressive weakness, with signs and symptoms of upper and lower motor neuron degeneration in multiple myotomes, without sensory neuropathy. These findings were supported on ancillary testing including nerve conduction studies and electromyography, muscle biopsies and muscle ultrasound studies. Biochemical investigations in plasma and fibroblasts showed elevated levels of ceramides and unrestrained de novo sphingolipid synthesis. Our studies indicate that SPTLC2 variant [c.778G>A, p.Glu260Lys] acts distinctly from hereditary sensory and autonomic neuropathy (HSAN)-causing SPTLC2 variants by causing excess canonical sphingolipid biosynthesis, similar to the recently reported SPTLC1 ALS associated pathogenic variants. Our studies also indicate that serine supplementation, which is a therapeutic in SPTLC1 and SPTCL2-associated HSAN, is expected to exacerbate the excess sphingolipid synthesis in serine palmitoyltransferase (SPT)-associated ALS.

CONCLUSIONS

SPTLC2 is the second SPT-associated gene that underlies monogenic, juvenile ALS and further establishes alterations of sphingolipid metabolism in motor neuron disease pathogenesis. Our findings also have important therapeutic implications: serine supplementation must be avoided in SPT-associated ALS, as it is expected to drive pathogenesis further.

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 GREGoR consortium. Each family's CNV data was analyzed using the seqr platform and candidate CNVs classified using the 2020 ACMG/ClinGen CNV interpretation standards. We developed additional evidence criteria to address situations not covered by the current standards. The addition of CNV calling to exome analysis identified causal CNVs for 173 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb with estimates that 44% would not have been detected by standard chromosomal microarrays. The causal CNVs consisted of 141 deletions, 15 duplications, 4 suspected complex structural variants (SVs), 3 insertions and 10 complex SVs, the latter two groups being identified by orthogonal validation methods. We interpreted 153 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.

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.

Bi-allelic variants in HMGCR cause an autosomal-recessive progressive limb-girdle muscular dystrophy

Statins are a mainstay intervention for cardiovascular disease prevention, yet their use can cause rare severe myopathy. HMG-CoA reductase, an essential enzyme in the mevalonate pathway, is the target of statins. We identified nine individuals from five unrelated families with unexplained limb-girdle like muscular dystrophy and bi-allelic variants in HMGCR via clinical and research exome sequencing. The clinical features resembled other genetic causes of muscular dystrophy with incidental high CPK levels (>1,000 U/L), proximal muscle weakness, variable age of onset, and progression leading to impaired ambulation. Muscle biopsies in most affected individuals showed non-specific dystrophic changes with non-diagnostic immunohistochemistry. Molecular modeling analyses revealed variants to be destabilizing and affecting protein oligomerization. Protein activity studies using three variants (p.Asp623Asn, p.Tyr792Cys, and p.Arg443Gln) identified in affected individuals confirmed decreased enzymatic activity and reduced protein stability. In summary, we showed that individuals with bi-allelic amorphic (i.e., null and/or hypomorphic) variants in HMGCR display phenotypes that resemble non-genetic causes of myopathy involving this reductase. This study expands our knowledge regarding the mechanisms leading to muscular dystrophy through dysregulation of the mevalonate pathway, autoimmune myopathy, and statin-induced myopathy.

Recognizable Pattern of Arthrogryposis and Congenital Myopathy Caused by the Recurrent TTN Metatranscript-only c.39974-11T > G Splice Variant

INTRODUCTION

Arthrogryposis is characterized by the presence of multiple contractures at birth and can be caused by pathogenic variants in TTN (Titin). Exons and variants that are not expressed in one of the three major isoforms of titin are referred to as "metatranscript-only" and have been considered to be only expressed during fetal development. Recently, the metatranscript-only variant (c.39974-11T > G) in TTN with a second truncating TTN variant has been linked to arthrogryposis multiplex congenita and myopathy.

METHODS

Via exome sequencing we identified the TTN c.39974-11T > G splice variant in trans with one of three truncating variants (p.Arg8922*, p.Lys32998Asnfs*63, p.Tyr10345*) in five individuals from three families. Clinical presentation and muscle ultrasound as well as MRI images were analyzed.

RESULTS

All five patients presented with generalized muscular hypotonia, reduced muscle bulk, and congenital contractures most prominently affecting the upper limbs and distal joints. Muscular hypotonia persisted and contractures improved over time. One individual, the recipient twin in the setting of twin-to-twin transfusion syndrome, died from severe cardiac hypertrophy 1 day after birth. Ultrasound and MRI imaging studies revealed a recognizable pattern of muscle involvement with striking fibrofatty involvement of the hamstrings and calves, and relative sparing of the femoral adductors and anterior segment of the thighs.

CONCLUSION

The recurrent TTN c.39974-11T > G variant consistently causes congenital arthrogryposis and persisting myopathy providing evidence that the metatranscript-only 213 to 217 exons impact muscle elasticity during early development and beyond. There is a recognizable pattern of muscle involvement, which is distinct from other myopathies and provides valuable clues for diagnostic work-up.

Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes

Genome-wide association studies have successfully discovered thousands of common variants associated with human diseases and traits, but the landscape of rare variations in human disease has not been explored at scale. Exome-sequencing studies of population biobanks provide an opportunity to systematically evaluate the impact of rare coding variations across a wide range of phenotypes to discover genes and allelic series relevant to human health and disease. Here, we present results from systematic association analyses of 4,529 phenotypes using single-variant and gene tests of 394,841 individuals in the UK Biobank with exome-sequence data. We find that the discovery of genetic associations is tightly linked to frequency and is correlated with metrics of deleteriousness and natural selection. We highlight biological findings elucidated by these data and release the dataset as a public resource alongside the Genebass browser for rapidly exploring rare-variant association results.

seqr: A web-based analysis and collaboration tool for rare disease genomics

Exome and genome sequencing have become the tools of choice for rare disease diagnosis, leading to large amounts of data available for analyses. To identify causal variants in these datasets, powerful filtering and decision support tools that can be efficiently used by clinicians and researchers are required. To address this need, we developed seqr - an open-source, web-based tool for family-based monogenic disease analysis that allows researchers to work collaboratively to search and annotate genomic callsets. To date, seqr is being used in several research pipelines and one clinical diagnostic lab. In our own experience through the Broad Institute Center for Mendelian Genomics, seqr has enabled analyses of over 10,000 families, supporting the diagnosis of more than 3,800 individuals with rare disease and discovery of over 300 novel disease genes. Here, we describe a framework for genomic analysis in rare disease that leverages seqr's capabilities for variant filtration, annotation, and causal variant identification, as well as support for research collaboration and data sharing. The seqr platform is available as open source software, allowing low-cost participation in rare disease research, and a community effort to support diagnosis and gene discovery in rare disease.

BET1 variants establish impaired vesicular transport as a cause for muscular dystrophy with epilepsy

BET1 is required, together with its SNARE complex partners GOSR2, SEC22b, and Syntaxin-5 for fusion of endoplasmic reticulum-derived vesicles with the ER-Golgi intermediate compartment (ERGIC) and the cis-Golgi. Here, we report three individuals, from two families, with severe congenital muscular dystrophy (CMD) and biallelic variants in BET1 (P1 p.(Asp68His)/p.(Ala45Valfs*2); P2 and P3 homozygous p.(Ile51Ser)). Due to aberrant splicing and frameshifting, the variants in P1 result in low BET1 protein levels and impaired ER-to-Golgi transport. Since in silico modeling suggested that p.(Ile51Ser) interferes with binding to interaction partners other than SNARE complex subunits, we set off and identified novel BET1 interaction partners with low affinity for p.(Ile51Ser) BET1 protein compared to wild-type, among them ERGIC-53. The BET1/ERGIC-53 interaction was validated by endogenous co-immunoprecipitation with both proteins colocalizing to the ERGIC compartment. Mislocalization of ERGIC-53 was observed in P1 and P2's derived fibroblasts; while in the p.(Ile51Ser) P2 fibroblasts specifically, mutant BET1 was also mislocalized along with ERGIC-53. Thus, we establish BET1 as a novel CMD/epilepsy gene and confirm the emerging role of ER/Golgi SNAREs in CMD.

The importance of automation in genetic diagnosis: Lessons from analyzing an inherited retinal degeneration cohort with the Mendelian Analysis Toolkit (MATK)

PURPOSE

In Mendelian disease diagnosis, variant analysis is a repetitive, error-prone, and time consuming process. To address this, we have developed the Mendelian Analysis Toolkit (MATK), a configurable, automated variant ranking program.

METHODS

MATK aggregates variant information from multiple annotation sources and uses expert-designed rules with parameterized weights to produce a ranked list of potentially causal solutions. MATK performance was measured by a comparison between MATK-aided and human-domain expert analyses of 1060 families with inherited retinal degeneration (IRD), analyzed using an IRD-specific gene panel (589 individuals) and exome sequencing (471 families).

RESULTS

When comparing MATK-assisted analysis with expert curation in both the IRD-specific gene panel and exome sequencing (1060 subjects), 97.3% of potential solutions found by experts were also identified by the MATK-assisted analysis (541 solutions identified with MATK of 556 solutions found by conventional analysis). Furthermore, MATK-assisted analysis identified 114 additional potential solutions from the 504 cases unsolved by conventional analysis.

CONCLUSION

MATK expedites the process of identification of likely solving variants in Mendelian traits, and reduces variability stemming from human error and researcher bias. MATK facilitates data reanalysis to keep up with the constantly improving annotation sources and next-generation sequencing processing pipelines. The software is open source and available at https://gitlab.com/matthew_maher/mendelanalysis.

A Hidden Structural Variation in a Known IRD Gene: A Cautionary Tale of Two New Disease Candidate Genes

Rod cone dystrophy (RCD), also known as retinitis pigmentosa, is an inherited condition leading to vision loss, affecting 1/3500 people. Over 270 genes are known to be implicated in the inherited retinal degenerations (IRDs), yet genetic diagnosis for ~30% IRD of patients remains elusive despite advances in sequencing technologies. The goal of this study was to determine the genetic causality in a family with Rod-cone dystrophy (RCD). Family members were given a full ophthalmic exam at the Retinal Service at MEE and consented to genetic testing. Whole exome sequencing (WES) was performed and variants of interest were Sanger validated. Functional assays were conducted in zebrafish along with splicing assays in relevant cell lines to determine the impact on retinal function. WES identified variants in two potential candidate genes that segregated with disease: GNL3 (G Protein Nucleolar 3) c.1187+3A>C and c.1568-8C>A; and PDE4DIP (Phosphodiester 4D Interacting Protein) c.3868G>A (p.Glu1290Lys) and c.4603G>A (p.Ala1535Thr). Both genes were promising candidates based on their retinal involvement (development and interactions with IRD-associated proteins), however the functional assays did not validate either gene. Subsequent WES reanalysis with an updated bioinformatics pipeline and widened search parameters led to the detection of a 94bp duplication in PRPF31 (pre-mRNA Processing Factor 31) c.73_266dup (p.Asp56GlyfsTer33) as the causal variant. Our study demonstrates the importance of thorough functional characterization of new disease candidate genes, and the value of reanalyzing NGS sequence data, which in our case led to identification of a hidden pathogenic variant in a known IRD gene.

Childhood amyotrophic lateral sclerosis caused by excess sphingolipid synthesis

Amyotrophic lateral sclerosis (ALS) is a progressive, neurodegenerative disease of the lower and upper motor neurons with sporadic or hereditary occurrence. Age of onset, pattern of motor neuron degeneration and disease progression vary widely among individuals with ALS. Various cellular processes may drive ALS pathomechanisms, but a monogenic direct metabolic disturbance has not been causally linked to ALS. Here we show SPTLC1 variants that result in unrestrained sphingoid base synthesis cause a monogenic form of ALS. We identified four specific, dominantly acting SPTLC1 variants in seven families manifesting as childhood-onset ALS. These variants disrupt the normal homeostatic regulation of serine palmitoyltransferase (SPT) by ORMDL proteins, resulting in unregulated SPT activity and elevated levels of canonical SPT products. Notably, this is in contrast with SPTLC1 variants that shift SPT amino acid usage from serine to alanine, result in elevated levels of deoxysphingolipids and manifest with the alternate phenotype of hereditary sensory and autonomic neuropathy. We custom designed small interfering RNAs that selectively target the SPTLC1 ALS allele for degradation, leave the normal allele intact and normalize sphingolipid levels in vitro. The role of primary metabolic disturbances in ALS has been elusive; this study defines excess sphingolipid biosynthesis as a fundamental metabolic mechanism for motor neuron disease.

Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

A form of muscular dystrophy associated with pathogenic variants in JAG2

JAG2 encodes the Notch ligand Jagged2. The conserved Notch signaling pathway contributes to the development and homeostasis of multiple tissues, including skeletal muscle. We studied an international cohort of 23 individuals with genetically unsolved muscular dystrophy from 13 unrelated families. Whole-exome sequencing identified rare homozygous or compound heterozygous JAG2 variants in all 13 families. The identified bi-allelic variants include 10 missense variants that disrupt highly conserved amino acids, a nonsense variant, two frameshift variants, an in-frame deletion, and a microdeletion encompassing JAG2. Onset of muscle weakness occurred from infancy to young adulthood. Serum creatine kinase (CK) levels were normal or mildly elevated. Muscle histology was primarily dystrophic. MRI of the lower extremities revealed a distinct, slightly asymmetric pattern of muscle involvement with cores of preserved and affected muscles in quadriceps and tibialis anterior, in some cases resembling patterns seen in POGLUT1-associated muscular dystrophy. Transcriptome analysis of muscle tissue from two participants suggested misregulation of genes involved in myogenesis, including PAX7. In complementary studies, Jag2 downregulation in murine myoblasts led to downregulation of multiple components of the Notch pathway, including Megf10. Investigations in Drosophila suggested an interaction between Serrate and Drpr, the fly orthologs of JAG1/JAG2 and MEGF10, respectively. In silico analysis predicted that many Jagged2 missense variants are associated with structural changes and protein misfolding. In summary, we describe a muscular dystrophy associated with pathogenic variants in JAG2 and evidence suggests a disease mechanism related to Notch pathway dysfunction.

Familial thrombocytopenia due to a complex structural variant resulting in a WAC-ANKRD26 fusion transcript

Advances in genome sequencing have resulted in the identification of the causes for numerous rare diseases. However, many cases remain unsolved with standard molecular analyses. We describe a family presenting with a phenotype resembling inherited thrombocytopenia 2 (THC2). THC2 is generally caused by single nucleotide variants that prevent silencing of ANKRD26 expression during hematopoietic differentiation. Short-read whole-exome and genome sequencing approaches were unable to identify a causal variant in this family. Using long-read whole-genome sequencing, a large complex structural variant involving a paired-duplication inversion was identified. Through functional studies, we show that this structural variant results in a pathogenic gain-of-function WAC-ANKRD26 fusion transcript. Our findings illustrate how complex structural variants that may be missed by conventional genome sequencing approaches can cause human disease.

Pathogenic Variants in the Myosin Chaperone UNC-45B Cause Progressive Myopathy with Eccentric Cores

The myosin-directed chaperone UNC-45B is essential for sarcomeric organization and muscle function from Caenorhabditis elegans to humans. The pathological impact of UNC-45B in muscle disease remained elusive. We report ten individuals with bi-allelic variants in UNC45B who exhibit childhood-onset progressive muscle weakness. We identified a common UNC45B variant that acts as a complex hypomorph splice variant. Purified UNC-45B mutants showed changes in folding and solubility. In situ localization studies further demonstrated reduced expression of mutant UNC-45B in muscle combined with abnormal localization away from the A-band towards the Z-disk of the sarcomere. The physiological relevance of these observations was investigated in C. elegans by transgenic expression of conserved UNC-45 missense variants, which showed impaired myosin binding for one and defective muscle function for three. Together, our results demonstrate that UNC-45B impairment manifests as a chaperonopathy with progressive muscle pathology, which discovers the previously unknown conserved role of UNC-45B in myofibrillar organization.

Exome sequencing identifies novel missense and deletion variants in RTN4IP1 associated with optic atrophy, global developmental delay, epilepsy, ataxia, and choreoathetosis

Inherited optic neuropathies (IONs) are neurodegenerative disorders characterized by optic atrophy with or without extraocular manifestations. Optic atrophy-10 (OPA10) is an autosomal recessive ION recently reported to be caused by mutations in RTN4IP1, which encodes reticulon 4 interacting protein 1 (RTN4IP1), a mitochondrial ubiquinol oxydo-reductase. Here we report novel compound heterozygous mutations in RTN4IP1 in a male proband with developmental delay, epilepsy, optic atrophy, ataxia, and choreoathetosis. Workup was notable for transiently elevated lactate and lactate-to-pyruvate ratio, brain magnetic resonance imaging with optic atrophy and T2 signal abnormalities, and a nondiagnostic initial genetic workup, including chromosomal microarray and mitochondrial panel testing. Exome sequencing identified a paternally inherited missense variant (c.263T>G, p.Val88Gly) predicted to be deleterious and a maternally inherited deletion encompassing RTN4IP1. To our knowledge, this is the first report of a non-single nucleotide pathogenic variant associated with OPA10. This case highlights the expanding phenotypic spectrum of OPA10, the association between "syndromic" cases and severe RTN4IP1 mutations, and the importance of nonbiased genetic testing, such as ES, to analyze multiple genes and variants types, in patients suspected of having genetic disease.

Biallelic loss of function variants in SYT2 cause a treatable congenital onset presynaptic myasthenic syndrome

Synaptotagmins are integral synaptic vesicle membrane proteins that function as calcium sensors and regulate neurotransmitter release at the presynaptic nerve terminal. Synaptotagmin-2 (SYT2), is the major isoform expressed at the neuromuscular junction. Recently, dominant missense variants in SYT2 have been reported as a rare cause of distal motor neuropathy and myasthenic syndrome, manifesting with stable or slowly progressive distal weakness of variable severity along with presynaptic NMJ impairment. These variants are thought to have a dominant-negative effect on synaptic vesicle exocytosis, although the precise pathomechanism remains to be elucidated. Here we report seven patients of five families, with biallelic loss of function variants in SYT2, clinically manifesting with a remarkably consistent phenotype of severe congenital onset hypotonia and weakness, with variable degrees of respiratory involvement. Electrodiagnostic findings were consistent with a presynaptic congenital myasthenic syndrome (CMS) in some. Treatment with an acetylcholinesterase inhibitor pursued in three patients showed clinical improvement with increased strength and function. This series further establishes SYT2 as a CMS-disease gene and expands its clinical and genetic spectrum to include recessive loss-of-function variants, manifesting as a severe congenital onset presynaptic CMS with potential treatment implications.

Expanding the phenotypic spectrum in RDH12-associated retinal disease

Retinol dehydrogenase 12, RDH12, plays a pivotal role in the visual cycle to ensure the maintenance of normal vision. Alterations in activity of this protein result in photoreceptor death and decreased vision beginning at an early age and progressing to substantial vision loss later in life. Here we describe 11 patients with retinal degeneration that underwent next-generation sequencing (NGS) with a targeted panel of all currently known inherited retinal degeneration (IRD) genes and whole-exome sequencing to identify the genetic causality of their retinal disease. These patients display a range of phenotypic severity prompting clinical diagnoses of macular dystrophy, cone-rod dystrophy, retinitis pigmentosa, and early-onset severe retinal dystrophy all attributed to biallelic recessive mutations in RDH12. We report 15 causal alleles and expand the repertoire of known RDH12 mutations with four novel variants: c.215A > G (p.Asp72Gly); c.362T > C (p.Ile121Thr); c.440A > C (p.Asn147Thr); and c.697G > A (p.Val233Ille). The broad phenotypic spectrum observed with biallelic RDH12 mutations has been observed in other genetic forms of IRDs, but the diversity is particularly notable here given the prior association of RDH12 primarily with severe early-onset disease. This breadth emphasizes the importance of broad genetic testing for inherited retinal disorders and extends the pool of individuals who may benefit from imminent gene-targeted therapies.

Recurrent TTN metatranscript-only c.39974-11T>G splice variant associated with autosomal recessive arthrogryposis multiplex congenita and myopathy

We present eight families with arthrogryposis multiplex congenita and myopathy bearing a TTN intron 213 extended splice-site variant (NM_001267550.1:c.39974-11T>G), inherited in trans with a second pathogenic TTN variant. Muscle-derived RNA studies of three individuals confirmed mis-splicing induced by the c.39974-11T>G variant; in-frame exon 214 skipping or use of a cryptic 3' splice-site effecting a frameshift. Confounding interpretation of pathogenicity is the absence of exons 213-217 within the described skeletal muscle TTN N2A isoform. However, RNA-sequencing from 365 adult human gastrocnemius samples revealed that 56% specimens predominantly include exons 213-217 in TTN transcripts (inclusion rate ≥66%). Further, RNA-sequencing of five fetal muscle samples confirmed that 4/5 specimens predominantly include exons 213-217 (fifth sample inclusion rate 57%). Contractures improved significantly with age for four individuals, which may be linked to decreased expression of pathogenic fetal transcripts. Our study extends emerging evidence supporting a vital developmental role for TTN isoforms containing metatranscript-only exons.

ACTN2 mutations cause "Multiple structured Core Disease" (MsCD)

The identification of genes implicated in myopathies is essential for diagnosis and for revealing novel therapeutic targets. Here we characterize a novel subclass of congenital myopathy at the morphological, molecular, and functional level. Through exome sequencing, we identified de novo ACTN2 mutations, a missense and a deletion, in two unrelated patients presenting with progressive early-onset muscle weakness and respiratory involvement. Morphological and ultrastructural analyses of muscle biopsies revealed a distinctive pattern with the presence of muscle fibers containing small structured cores and jagged Z-lines. Deeper analysis of the missense mutation revealed mutant alpha-actinin-2 properly localized to the Z-line in differentiating myotubes and its level was not altered in muscle biopsy. Modelling of the disease in zebrafish and mice by exogenous expression of mutated alpha-actinin-2 recapitulated the abnormal muscle function and structure seen in the patients. Motor deficits were noted in zebrafish, and muscle force was impaired in isolated muscles from AAV-transduced mice. In both models, sarcomeric disorganization was evident, while expression of wild-type alpha-actinin-2 did not result in muscle anomalies. The murine muscles injected with mutant ACTN2 displayed cores and Z-line defects. Dominant ACTN2 mutations were previously associated with cardiomyopathies, and our data demonstrate that specific mutations in the well-known Z-line regulator alpha-actinin-2 can cause a skeletal muscle disorder.

The Genetic Landscape of Diamond-Blackfan Anemia

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.

MACF1 Mutations Encoding Highly Conserved Zinc-Binding Residues of the GAR Domain Cause Defects in Neuronal Migration and Axon Guidance

To date, mutations in 15 actin- or microtubule-associated genes have been associated with the cortical malformation lissencephaly and variable brainstem hypoplasia. During a multicenter review, we recognized a rare lissencephaly variant with a complex brainstem malformation in three unrelated children. We searched our large brain-malformation databases and found another five children with this malformation (as well as one with a less severe variant), analyzed available whole-exome or -genome sequencing data, and tested ciliogenesis in two affected individuals. The brain malformation comprised posterior predominant lissencephaly and midline crossing defects consisting of absent anterior commissure and a striking W-shaped brainstem malformation caused by small or absent pontine crossing fibers. We discovered heterozygous de novo missense variants or an in-frame deletion involving highly conserved zinc-binding residues within the GAR domain of MACF1 in the first eight subjects. We studied cilium formation and found a higher proportion of mutant cells with short cilia than of control cells with short cilia. A ninth child had similar lissencephaly but only subtle brainstem dysplasia associated with a heterozygous de novo missense variant in the spectrin repeat domain of MACF1. Thus, we report variants of the microtubule-binding GAR domain of MACF1 as the cause of a distinctive and most likely pathognomonic brain malformation. A gain-of-function or dominant-negative mechanism appears likely given that many heterozygous mutations leading to protein truncation are included in the ExAC Browser. However, three de novo variants in MACF1 have been observed in large schizophrenia cohorts.

Novel SPEG mutations in congenital myopathies: Genotype-phenotype correlations

INTRODUCTION

Centronuclear myopathies (CNMs) are a subtype of congenital myopathies (CMs) characterized by muscle weakness, predominant type 1 fibers, and increased central nuclei. SPEG (striated preferentially expressed protein kinase) mutations have recently been identified in 7 CM patients (6 with CNMs). We report 2 additional patients with SPEG mutations expanding the phenotype and evaluate genotype-phenotype correlations associated with SPEG mutations.

METHODS

Using whole exome/genome sequencing in CM families, we identified novel recessive SPEG mutations in 2 patients.

RESULTS

Patient 1, with severe muscle weakness requiring respiratory support, dilated cardiomyopathy, ophthalmoplegia, and findings of nonspecific CM on muscle biopsy carried a homozygous SPEG mutation (p.Val3062del). Patient 2, with milder muscle weakness, ophthalmoplegia, and CNM carried compound heterozygous mutations (p.Leu728Argfs*82) and (p.Val2997Glyfs*52).

CONCLUSIONS

The 2 patients add insight into genotype-phenotype correlations of SPEG-associated CMs. Clinicians should consider evaluating a CM patient for SPEG mutations even in the absence of CNM features. Muscle Nerve 59:357-362, 2019.

A Syndromic Neurodevelopmental Disorder Caused by De Novo Variants in EBF3

Early B cell factor 3 (EBF3) is a member of the highly evolutionarily conserved Collier/Olf/EBF (COE) family of transcription factors. Prior studies on invertebrate and vertebrate animals have shown that EBF3 homologs are essential for survival and that loss-of-function mutations are associated with a range of nervous system developmental defects, including perturbation of neuronal development and migration. Interestingly, aristaless-related homeobox (ARX), a homeobox-containing transcription factor critical for the regulation of nervous system development, transcriptionally represses EBF3 expression. However, human neurodevelopmental disorders related to EBF3 have not been reported. Here, we describe three individuals who are affected by global developmental delay, intellectual disability, and expressive speech disorder and carry de novo variants in EBF3. Associated features seen in these individuals include congenital hypotonia, structural CNS malformations, ataxia, and genitourinary abnormalities. The de novo variants affect a single conserved residue in a zinc finger motif crucial for DNA binding and are deleterious in a fly model. Our findings indicate that mutations in EBF3 cause a genetic neurodevelopmental syndrome and suggest that loss of EBF3 function might mediate a subset of neurologic phenotypes shared by ARX-related disorders, including intellectual disability, abnormal genitalia, and structural CNS malformations.

De Novo Truncating Variants in ASXL2 Are Associated with a Unique and Recognizable Clinical Phenotype

The ASXL genes (ASXL1, ASXL2, and ASXL3) participate in body patterning during embryogenesis and encode proteins involved in epigenetic regulation and assembly of transcription factors to specific genomic loci. Germline de novo truncating variants in ASXL1 and ASXL3 have been respectively implicated in causing Bohring-Opitz and Bainbridge-Ropers syndromes, which result in overlapping features of severe intellectual disability and dysmorphic features. ASXL2 has not yet been associated with a human Mendelian disorder. In this study, we performed whole-exome sequencing in six unrelated probands with developmental delay, macrocephaly, and dysmorphic features. All six had de novo truncating variants in ASXL2. A careful review enabled the recognition of a specific phenotype consisting of macrocephaly, prominent eyes, arched eyebrows, hypertelorism, a glabellar nevus flammeus, neonatal feeding difficulties, hypotonia, and developmental disabilities. Although overlapping features with Bohring-Opitz and Bainbridge-Ropers syndromes exist, features that distinguish the ASXL2-associated condition from ASXL1- and ASXL3-related disorders are macrocephaly, absence of growth retardation, and more variability in the degree of intellectual disabilities. We were also able to demonstrate with mRNA studies that these variants are likely to exert a dominant-negative effect, given that both alleles are expressed in blood and the mutated ASXL2 transcripts escape nonsense-mediated decay. In conclusion, de novo truncating variants in ASXL2 underlie a neurodevelopmental syndrome with a clinically recognizable phenotype. This report expands the germline disorders that are linked to the ASXL genes.

Disruption of Golgi morphology and altered protein glycosylation in PLA2G6-associated neurodegeneration

BACKGROUND

Mutations in PLA2G6, which encodes the calcium-independent phospholipase A2 group VI, cause neurodegeneration and diffuse cortical Lewy body formation by a yet undefined mechanism. We assessed whether altered protein glycosylation due to abnormal Golgi morphology might be a factor in the pathology of this disease.

METHODS

Three patients presented with PLA2G6-associated neurodegeneration (PLAN); two had infantile neuroaxonal dystrophy (INAD) and one had adult-onset dystonia-parkinsonism. We analysed protein N-linked and O-linked glycosylation in cerebrospinal fluid, plasma, urine, and cultured skin fibroblasts using high performance liquid chromatography (HPLC) and matrix-assisted laser desorption ionization--time of flight/mass spectrometry (MALDI-TOF/MS). We also assessed sialylation and Golgi morphology in cultured fibroblasts by immunofluorescence and performed rescue experiments using a lentiviral vector.

RESULTS

The patients with INAD had PLA2G6 mutations NM_003560.2: c.[950G>T];[426-1077dup] and c.[1799G>A];[2221C>T] and the patient with dystonia-parkinsonism had PLA2G6 mutations NM_003560.2: c.[609G>A];[2222G>A]. All three patients had altered Golgi morphology and abnormalities of protein O-linked glycosylation and sialylation in cultured fibroblasts that were rescued by lentiviral overexpression of wild type PLA2G6.

CONCLUSIONS

Our findings add altered Golgi morphology, O-linked glycosylation and sialylation defects to the phenotypical spectrum of PLAN; these pathways are essential for correct processing and distribution of proteins. Lewy body and Tau pathology, two neuropathological features of PLAN, could emerge from these defects. Therefore, Golgi morphology, O-linked glycosylation and sialylation may play a role in the pathogenesis of PLAN and perhaps other neurodegenerative disorders.