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.

A harmonized public resource of deeply sequenced diverse human genomes.

Underrepresented populations are often excluded from genomic studies owing 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 4094 whole genomes from 80 populations in the 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 show substantial added value from this data set compared with 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 with 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.

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.

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.

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.

Germline GATA1s-generating mutations predispose to leukemia with acquired trisomy 21 and Down syndrome-like phenotype.

Individuals with Down syndrome are at increased risk of myeloid leukemia in early childhood, which is associated with acquisition of GATA1 mutations that generate a short GATA1 isoform called GATA1s. Germline GATA1s-generating mutations result in congenital anemia in males. We report on 2 unrelated families that harbor germline GATA1s-generating mutations in which several members developed acute megakaryoblastic leukemia in early childhood. All evaluable leukemias had acquired trisomy 21 or tetrasomy 21. The leukemia characteristics overlapped with those of myeloid leukemia associated with Down syndrome, including age of onset at younger than 4 years, unique immunophenotype, complex karyotype, gene expression patterns, and drug sensitivity. These findings demonstrate that the combination of trisomy 21 and GATA1s-generating mutations results in a unique myeloid leukemia independent of whether the GATA1 mutation or trisomy 21 is the primary or secondary event and suggest that there is a unique functional cooperation between GATA1s and trisomy 21 in leukemogenesis. The family histories also indicate that germline GATA1s-generating mutations should be included among those associated with familial predisposition for myelodysplastic syndrome and leukemia.

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.

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.

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.