The recurrent deep intronic pseudoexon-inducing variant COL6A1 c.930+189C>T results in a consistently severe phenotype of COL6-related dystrophy: Towards clinical trial readiness for splice-modulating therapy.

Collagen VI-related dystrophies (COL6-RDs) manifest with a spectrum of clinical phenotypes, ranging from Ullrich congenital muscular dystrophy (UCMD), presenting with prominent congenital symptoms and characterised by progressive muscle weakness, joint contractures and respiratory insufficiency, to Bethlem muscular dystrophy, with milder symptoms typically recognised later and at times resembling a limb girdle muscular dystrophy, and intermediate phenotypes falling between UCMD and Bethlem muscular dystrophy. Despite clinical and immunohistochemical features highly suggestive of COL6-RD, some patients had remained without an identified causative variant in COL6A1, COL6A2 or COL6A3. With combined muscle RNA-sequencing and whole-genome sequencing we uncovered a recurrent, de novo deep intronic variant in intron 11 of COL6A1 (c.930+189C>T) that leads to a dominantly acting in-frame pseudoexon insertion. We subsequently identified and have characterised an international cohort of forty-four patients with this COL6A1 intron 11 causative variant, one of the most common recurrent causative variants in the collagen VI genes. Patients manifest a consistently severe phenotype characterised by a paucity of early symptoms followed by an accelerated progression to a severe form of UCMD, except for one patient with somatic mosaicism for this COL6A1 intron 11 variant who manifests a milder phenotype consistent with Bethlem muscular dystrophy. Characterisation of this individual provides a robust validation for the development of our pseudoexon skipping therapy. We have previously shown that splice-modulating antisense oligomers applied in vitro effectively decreased the abundance of the mutant pseudoexon-containing COL6A1 transcripts to levels comparable to the in vivo scenario of the somatic mosaicism shown here, indicating that this therapeutic approach carries significant translational promise for ameliorating the severe form of UCMD caused by this common recurrent COL6A1 causative variant to a Bethlem muscular dystrophy phenotype.

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.

Relationship of fat mass ratio - a biomarker for lipodystrophy - with cardiometabolic traits.

Familial partial lipodystrophy (FPLD) is a heterogenous group of syndromes associated with a high prevalence of cardiometabolic diseases. Prior work has proposed DEXA-derived fat mass ratio (FMR) - defined as trunk fat percentage (trunk fat %) divided by leg fat percentage (leg fat %) - as a biomarker of FPLD, but this metric has not previously been characterized in large cohort studies. We set out to (1) understand the cardiometabolic burden of individuals with high FMR in up to 40,796 participants in the UK Biobank and 9,408 participants in the Fenland study, (2) characterize the common variant genetic underpinnings of FMR, and (3) build and test a polygenic predictor for FMR. Participants with high FMR were at higher risk for type 2 diabetes (OR = 2.30, p = 3.5 × 10-41) and MASLD/MASH (OR = 2.55, p = 4.9 × 10-7) in UK Biobank, and had higher fasting insulin (difference = +19.8 pmol/L, p = 5.7 × 10-36) and fasting triglycerides (difference = +36.1 mg/dL, p = 2.5 × 10-28) in the Fenland Study. Across FMR and its component traits, 61 conditionally independent variant-trait pairs were discovered, including 13 newly-identified pairs. A polygenic score for FMR was associated with increased risk of cardiometabolic diseases. This work establishes the cardiometabolic significance of high FMR - a biomarker for FPLD - in two large cohort studies and may prove useful in increasing diagnosis rates of patients with metabolically unhealthy fat distribution to enable treatment or a preventive therapy.

Small-molecule inhibition of glycogen synthase 1 for the treatment of Pompe disease and other glycogen storage disorders.

Glycogen synthase 1 (GYS1), the rate-limiting enzyme in muscle glycogen synthesis, plays a central role in energy homeostasis and has been proposed as a therapeutic target in multiple glycogen storage diseases. Despite decades of investigation, there are no known potent, selective small-molecule inhibitors of this enzyme. Here, we report the preclinical characterization of MZ-101, a small molecule that potently inhibits GYS1 in vitro and in vivo without inhibiting GYS2, a related isoform essential for synthesizing liver glycogen. Chronic treatment with MZ-101 depleted muscle glycogen and was well tolerated in mice. Pompe disease, a glycogen storage disease caused by mutations in acid α glucosidase (GAA), results in pathological accumulation of glycogen and consequent autophagolysosomal abnormalities, metabolic dysregulation, and muscle atrophy. Enzyme replacement therapy (ERT) with recombinant GAA is the only approved treatment for Pompe disease, but it requires frequent infusions, and efficacy is limited by suboptimal skeletal muscle distribution. In a mouse model of Pompe disease, chronic oral administration of MZ-101 alone reduced glycogen buildup in skeletal muscle with comparable efficacy to ERT. In addition, treatment with MZ-101 in combination with ERT had an additive effect and could normalize muscle glycogen concentrations. Biochemical, metabolomic, and transcriptomic analyses of muscle tissue demonstrated that lowering of glycogen concentrations with MZ-101, alone or in combination with ERT, corrected the cellular pathology in this mouse model. These data suggest that substrate reduction therapy with GYS1 inhibition may be a promising therapeutic approach for Pompe disease and other glycogen storage diseases.

A cryptic pathogenic NDUFV1 variant identified by RNA-seq in a patient with normal complex I activity in muscle and transient magnetic resonance imaging changes.

Mitochondrial respiratory chain disorders (MRC) are amongst the most common group of inborn errors of metabolism. MRC, of which complex I deficiency accounts for approximately a quarter, are very diverse, causing a wide range of clinical problems and can be difficult to diagnose. We report an illustrative MRC case whose diagnosis was elusive. Clinical signs included failure to thrive caused by recurrent vomiting, hypotonia and progressive loss of motor milestones. Initial brain imaging suggested Leigh syndrome but without expected diffusion restriction. Muscle respiratory chain enzymology was unremarkable. Whole-genome sequencing identified a maternally inherited NDUFV1 missense variant [NM_007103.4 (NDUFV1):c.1157G > A; p.(Arg386His)] and a paternally inherited synonymous variant [NM_007103.4 (NDUFV1):c.1080G > A; (p.Ser360=)]. RNA sequencing demonstrated aberrant splicing. This case emphasizes the diagnostic odyssey of a patient in whom a confirmed diagnosis was elusive because of atypical features and normal muscle respiratory chain enzyme (RCE) activities, along with a synonymous variant, which are often filtered out from genomic analyses. It also illustrates the following points: (1) complete resolution of magnetic resonance imaging changes may be part of the picture in mitochondrial disease; (2) analysis for synonymous variants is important for undiagnosed patients; and (3) RNA-seq is a powerful tool to demonstrate pathogenicity of putative splicing variants.

Transcriptome variation in human tissues revealed by long-read sequencing.

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

TDP-43 represses cryptic exon inclusion in the FTD-ALS gene UNC13A.

A hallmark pathological feature of the neurodegenerative diseases amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is the depletion of RNA-binding protein TDP-43 from the nucleus of neurons in the brain and spinal cord1. A major function of TDP-43 is as a repressor of cryptic exon inclusion during RNA splicing2-4. Single nucleotide polymorphisms in UNC13A are among the strongest hits associated with FTD and ALS in human genome-wide association studies5,6, but how those variants increase risk for disease is unknown. Here we show that TDP-43 represses a cryptic exon-splicing event in UNC13A. Loss of TDP-43 from the nucleus in human brain, neuronal cell lines and motor neurons derived from induced pluripotent stem cells resulted in the inclusion of a cryptic exon in UNC13A mRNA and reduced UNC13A protein expression. The top variants associated with FTD or ALS risk in humans are located in the intron harbouring the cryptic exon, and we show that they increase UNC13A cryptic exon splicing in the face of TDP-43 dysfunction. Together, our data provide a direct functional link between one of the strongest genetic risk factors for FTD and ALS (UNC13A genetic variants), and loss of TDP-43 function.

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.

WGS and RNA Studies Diagnose Noncoding DMD Variants in Males With High Creatine Kinase.

OBJECTIVE: To describe the diagnostic utility of whole-genome sequencing and RNA studies in boys with suspected dystrophinopathy, for whom multiplex ligation-dependent probe amplification and exomic parallel sequencing failed to yield a genetic diagnosis, and to use remnant normal DMD splicing in 3 families to define critical levels of wild-type dystrophin bridging clinical spectrums of Duchenne to myalgia. METHODS: Exome, genome, and/or muscle RNA sequencing was performed for 7 males with elevated creatine kinase. PCR of muscle-derived complementary DNA (cDNA) studied consequences for DMD premessenger RNA (pre-mRNA) splicing. Quantitative Western blot was used to determine levels of dystrophin, relative to control muscle. RESULTS: Splice-altering intronic single nucleotide variants or structural rearrangements in DMD were identified in all 7 families. Four individuals, with abnormal splicing causing a premature stop codon and nonsense-mediated decay, expressed remnant levels of normally spliced DMD mRNA. Quantitative Western blot enabled correlation of wild-type dystrophin and clinical severity, with 0%-5% dystrophin conferring a Duchenne phenotype, 10% ± 2% a Becker phenotype, and 15% ± 2% dystrophin associated with myalgia without manifesting weakness. CONCLUSIONS: Whole-genome sequencing relied heavily on RNA studies to identify DMD splice-altering variants. Short-read RNA sequencing was regularly confounded by the effectiveness of nonsense-mediated mRNA decay and low read depth of the giant DMD mRNA. PCR of muscle cDNA provided a simple, yet informative approach. Highly relevant to genetic therapies for dystrophinopathies, our data align strongly with previous studies of mutant dystrophin in Becker muscular dystrophy, with the collective conclusion that a fractional increase in levels of normal dystrophin between 5% and 20% is clinically significant.

Pathogenic deep intronic MTM1 variant activates a pseudo-exon encoding a nonsense codon resulting in severe X-linked myotubular myopathy.

X-linked myotubular myopathy (XLMTM) is a severe congenital myopathy characterised by generalised weakness and respiratory insufficiency. XLMTM is associated with pathogenic variants in MTM1; a gene encoding the lipid phosphatase myotubularin. Whole genome sequencing (WGS) of an exome-negative male proband with severe hypotonia, respiratory insufficiency and centralised nuclei on muscle biopsy identified a deep intronic MTM1 variant NG_008199.1(NM_000252.2):c.1468-577A>G, which strengthened a cryptic 5' splice site (A>G substitution at the +5 position). Muscle RNA sequencing was non-diagnostic due to low read depth. Reverse transcription PCR (RT-PCR) of muscle RNA confirmed the c.1468-577A>G variant activates inclusion of a pseudo-exon encoding a premature stop codon into all detected MTM1 transcripts. Western blot analysis establishes deficiency of myotubularin protein, consistent with the severe XLMTM phenotype. We expand the genotypic spectrum of XLMTM and highlight benefits of screening non-coding regions of MTM1 in male probands with phenotypically concordant XLMTM who remain undiagnosed following exome sequencing.

Landscape of multi-nucleotide variants in 125,748 human exomes and 15,708 genomes.

Multi-nucleotide variants (MNVs), defined as two or more nearby variants existing on the same haplotype in an individual, are a clinically and biologically important class of genetic variation. However, existing tools typically do not accurately classify MNVs, and understanding of their mutational origins remains limited. Here, we systematically survey MNVs in 125,748 whole exomes and 15,708 whole genomes from the Genome Aggregation Database (gnomAD). We identify 1,792,248 MNVs across the genome with constituent variants falling within 2 bp distance of one another, including 18,756 variants with a novel combined effect on protein sequence. Finally, we estimate the relative impact of known mutational mechanisms - CpG deamination, replication error by polymerase zeta, and polymerase slippage at repeat junctions - on the generation of MNVs. Our results demonstrate the value of haplotype-aware variant annotation, and refine our understanding of genome-wide mutational mechanisms of MNVs.

Evaluating drug targets through human loss-of-function genetic variation.

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous 'knockout' humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development.

Transcript expression-aware annotation improves rare variant interpretation.

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)1, we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.