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.

The expanding clinical and genetic spectrum of DYNC1H1-related disorders.

Intracellular trafficking involves an intricate machinery of motor complexes including the dynein complex to shuttle cargo for autophagolysosomal degradation. Deficiency in dynein axonemal chains as well as cytoplasmic light and intermediate chains have been linked with ciliary dyskinesia and skeletal dysplasia. The cytoplasmic dynein 1 heavy chain protein (DYNC1H1) serves as a core complex for retrograde trafficking in neuronal axons. Dominant pathogenic variants in DYNC1H1 have been previously implicated in peripheral neuromuscular disorders (NMD) and neurodevelopmental disorders (NDD). As heavy-chain dynein is ubiquitously expressed, the apparent selectivity of heavy-chain dyneinopathy for motor neuronal phenotypes remains currently unaccounted for. Here, we aimed to evaluate the full DYNC1H1-related clinical, molecular and imaging spectrum, including multisystem features and novel phenotypes presenting throughout life. We identified 47 cases from 43 families with pathogenic heterozygous variants in DYNC1H1 (aged 0-59 years) and collected phenotypic data via a comprehensive standardized survey and clinical follow-up appointments. Most patients presented with divergent and previously unrecognized neurological and multisystem features, leading to significant delays in genetic testing and establishing the correct diagnosis. Neurological phenotypes include novel autonomic features, previously rarely described behavioral disorders, movement disorders, and periventricular lesions. Sensory neuropathy was identified in nine patients (median age of onset 10.6 years), of which five were only diagnosed after the second decade of life, and three had a progressive age-dependent sensory neuropathy. Novel multisystem features included primary immunodeficiency, bilateral sensorineural hearing loss, organ anomalies, and skeletal manifestations, resembling the phenotypic spectrum of other dyneinopathies. We also identified an age-dependent biphasic disease course with developmental regression in the first decade and, following a period of stability, neurodegenerative progression after the second decade of life. Of note, we observed several cases in whom neurodegeneration appeared to be prompted by intercurrent systemic infections with double-stranded DNA viruses (Herpesviridae) or single-stranded RNA viruses (Ross-River fever, SARS-CoV-2). Moreover, the disease course appeared to be exacerbated by viral infections regardless of age and/or severity of NDD manifestations, indicating a role of dynein in anti-viral immunity and neuronal health. In summary, our findings expand the clinical, imaging, and molecular spectrum of pathogenic DYNC1H1 variants beyond motor neuropathy disorders and suggest a life-long continuum and age-related progression due to deficient intracellular trafficking. This study will facilitate early diagnosis and improve counselling and health surveillance of affected patients.

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

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

Zebrafish and cellular models of SELENON-Related Myopathy exhibit novel embryonic and metabolic phenotypes.

SELENON-Related Myopathy (SELENON-RM) is a rare congenital myopathy caused by mutations of the SELENON gene characterized by axial muscle weakness and progressive respiratory insufficiency. Muscle histopathology commonly includes multiminicores or a dystrophic pattern but is often non-specific. The SELENON gene encodes selenoprotein N (SelN), a selenocysteine-containing redox enzyme located in the endo/sarcoplasmic reticulum membrane where it colocalizes with mitochondria-associated membranes. However, the molecular mechanism(s) by which SelN deficiency causes SELENON-RM are undetermined. A hurdle is the lack of cellular and animal models that show assayable phenotypes. Here we report deep-phenotyping of SelN-deficient zebrafish and muscle cells. SelN-deficient zebrafish exhibit changes in embryonic muscle function and swimming activity in larvae. Analysis of single cell RNAseq data in a zebrafish embryo-atlas revealed coexpression between selenon and genes involved in glutathione redox pathway. SelN-deficient zebrafish and mouse myoblasts exhibit changes in glutathione and redox homeostasis, suggesting a direct relationship with SelN function. We report changes in metabolic function abnormalities in SelN-null myotubes when compared to WT. These results suggest that SelN has functional roles during zebrafish early development and myoblast metabolism.

Integrated multi-omics approach reveals the role of striated muscle preferentially expressed protein kinase in skeletal muscle including its relationship with myospryn complex.

BACKGROUND: Autosomal-recessive mutations in SPEG (striated muscle preferentially expressed protein kinase) have been linked to centronuclear myopathy with or without dilated cardiomyopathy (CNM5). Loss of SPEG is associated with defective triad formation, abnormal excitation-contraction coupling, calcium mishandling and disruption of the focal adhesion complex in skeletal muscles. To elucidate the underlying molecular pathways, we have utilized multi-omics tools and analysis to obtain a comprehensive view of the complex biological processes and molecular functions. METHODS: Skeletal muscles from 2-month-old SPEG-deficient (Speg-CKO) and wild-type (WT) mice were used for RNA sequencing (n = 4 per genotype) to profile transcriptomics and mass spectrometry (n = 4 for WT; n = 3 for Speg-CKO mice) to profile proteomics and phosphoproteomics. In addition, interactomics was performed using the SPEG antibody on pooled muscle lysates (quadriceps, gastrocnemius and triceps) from WT and Speg-CKO mice. Based on the multi-omics results, we performed quantitative real-time PCR, co-immunoprecipitation and immunoblot to verify the findings. RESULTS: We identified that SPEG interacts with myospryn complex proteins CMYA5, FSD2 and RyR1, which are critical for triad formation, and that SPEG deficiency results in myospryn complex abnormalities (protein levels decreased to 22 ± 3% for CMYA5 [P < 0.05] and 18 ± 3% for FSD2 [P < 0.01]). Furthermore, SPEG phosphorylates RyR1 at S2902 (phosphorylation level decreased to 55 ± 15% at S2902 in Speg-CKO mice; P < 0.05), and its loss affects JPH2 phosphorylation at multiple sites (increased phosphorylation at T161 [1.90 ± 0.24-fold], S162 [1.61 ± 0.37-fold] and S165 [1.66 ± 0.13-fold]; decreased phosphorylation at S228 and S231 [39 ± 6%], S234 [50 ± 12%], S593 [48 ± 3%] and S613 [66 ± 10%]; P < 0.05 for S162 and P < 0.01 for other sites). On analysing the transcriptome, the most dysregulated pathways affected by SPEG deficiency included extracellular matrix-receptor interaction (P < 1e-15) and peroxisome proliferator-activated receptor signalling (P < 9e-14). CONCLUSIONS: We have elucidated the critical role of SPEG in the triad as it works closely with myospryn complex proteins (CMYA5, FSD2 and RyR1), it regulates phosphorylation levels of various residues in JPH2 and S2902 in RyR1, and its deficiency is associated with dysregulation of several pathways. The study identifies unique SPEG-interacting proteins and their phosphorylation functions and emphasizes the importance of using a multi-omics approach to comprehensively evaluate the molecular function of proteins involved in various genetic disorders.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.