Two alternatively-spliced human nebulin isoforms with either exon 143 or exon 144 and their developmental regulation

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 brothers with congenital myopathy and compound heterozygous variants (NC_000002.12:g.151692086G>T; NM_001271208.2: c.2079C>A; p.(Cys693Ter) and NC_000002.12:g.151533439T>C; NM_001271208.2:c.21522+3A>G) in NEB. Transcriptomic sequencing on affected individual muscles 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. Affected individuals' MRI patterns of muscle involvement were compared with the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these affected individuals 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. Our report introduces disease pathogenesis and manifestation as a result of alteration of isoform distributions in muscle.

Comprehensive phenotypic characterization of an allelic series of zebrafish models of NEB-related nemaline myopathy

Nemaline myopathy (NM) is a rare congenital neuromuscular disorder characterized by muscle weakness and hypotonia, slow gross motor development, and decreased respiratory function. Mutations in at least twelve genes, all of each encode proteins that are either components of the muscle thin filament or regulate its length and stability, have been associated with NM. Mutations in Nebulin (NEB), a giant filamentous protein localized in the sarcomere, account for more than 50% of NM cases. At present, there remains a lack of understanding of whether NEB genotype influences nebulin function and NM-patient phenotypes. In addition, there is a lack of therapeutically tractable models that can enable drug discovery and address the current unmet treatment needs of patients. To begin to address these gaps, here we have characterized five new zebrafish models of NEB-related NM. These mutants recapitulate most aspects of NEB-based NM, showing drastically reduced survival, defective muscle structure, reduced contraction force, shorter thin filaments, presence of electron-dense structures in myofibers, and thickening of the Z-disks. This study represents the first extensive investigation of an allelic series of nebulin mutants, and thus provides an initial examination in pre-clinical models of potential genotype-phenotype correlations in human NEB patients. It also represents the first utilization of a set of comprehensive outcome measures in zebrafish, including correlation between molecular analyses, structural and biophysical investigations, and phenotypic outcomes. Therefore, it provides a rich source of data for future studies exploring the NM pathomechanisms, and an ideal springboard for therapy identification and development for NEB-related NM.

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.

Transcriptomic analysis of paired healthy human skeletal muscles to identify modulators of disease severity in DMD

Muscle damage and fibro-fatty replacement of skeletal muscles is a main pathologic feature of Duchenne muscular dystrophy (DMD) with more proximal muscles affected earlier and more distal affected later in the disease course, suggesting that different skeletal muscle groups possess distinctive characteristics that influence their susceptibility to disease. To explore transcriptomic factors driving differential gene expression and modulating DMD skeletal muscle severity, we characterized the transcriptome of vastus lateralis (VL), a more proximal and susceptible muscle, relative to tibialis anterior (TA), a more distal and protected muscle, in 15 healthy individuals using bulk RNA sequencing to identify gene expression differences that may mediate their relative susceptibility to damage with loss of dystrophin. Matching single nuclei RNA sequencing data was generated for 3 of the healthy individuals, to infer cell composition in the bulk RNA sequencing dataset and to improve mapping of differentially expressed genes to their cell source of expression. A total of 3,410 differentially expressed genes were identified and mapped to cell type using single nuclei RNA sequencing of muscle, including long non-coding RNAs and protein coding genes. There was an enrichment of genes involved in calcium release from the sarcoplasmic reticulum, particularly in the myofibers and these myofiber genes were higher in the VL. There was an enrichment of genes in "Collagen-Containing Extracellular Matrix" expressed by fibroblasts, endothelial, smooth muscle and pericytes, with most genes higher in the TA, as well as genes in "Regulation Of Apoptotic Process" expressed across all cell types. Previously reported genetic modifiers were also enriched within the differentially expressed genes. We also identify 6 genes with differential isoform usage between the VL and TA. Lastly, we integrate our findings with DMD RNA sequencing data from the TA, and identify "Collagen-Containing Extracellular Matrix" and "Negative Regulation Of Apoptotic Process" as differentially expressed between DMD compared to healthy. Collectively, these findings propose novel candidate mechanisms that may mediate differential muscle susceptibility in muscular dystrophies and provide new insight into potential therapeutic targets.

Early clinical and pre-clinical therapy development in Nemaline myopathy

INTRODUCTION

Nemaline myopathies (NM) represent a group of clinically and genetically heterogeneous congenital muscle disorders with the common denominator of nemaline rods on muscle biopsy. NEB and ACTA1 are the most common causative genes. Currently, available treatments are supportive.

AREAS COVERED

We explored experimental treatments for NM, identifying at least eleven mainly pre-clinical approaches utilizing murine and/or human muscle cells. These approaches target either i) the causative gene or associated genes implicated in the same pathway; ii) pathophysiologically relevant biochemical mechanisms such as calcium/myosin regulation of muscle contraction; iii) myogenesis; iv) other therapies that improve or optimize muscle function more generally; v) and/or combinations of the above. The scope and efficiency of these attempts is diverse, ranging from gene-specific effects to those widely applicable to all NM-associated genes.

EXPERT OPINION

The wide range of experimental therapies currently under consideration for NM is promising. Potential translation into clinical use requires consideration of additional factors such as the potential muscle type specificity as well as the possibility of gene expression remodeling. Challenges in clinical translation include the rarity and heterogeneity of genotypes, phenotypes, and disease trajectories, as well as the lack of longitudinal natural history data and validated outcomes and biomarkers.

Exploring the Potential of Symmetric Exon Deletion to Treat Non-Ischemic Dilated Cardiomyopathy by Removing Frameshift Mutations in TTN

Non-ischemic dilated cardiomyopathy (DCM) is one of the most frequent pathologies requiring cardiac transplants. Even though the etiology of this disease is complex, frameshift mutations in the giant sarcomeric protein Titin could explain up to 25% of the familial and 18% of the sporadic cases of DCM. Many studies have shown the potential of genome editing using CRISPR/Cas9 to correct truncating mutations in sarcomeric proteins and have established the grounds for myoediting. However, these therapies are still in an immature state, with only few studies showing an efficient treatment of cardiac diseases. This publication hypothesizes that the Titin (TTN)-specific gene structure allows the application of myoediting approaches in a broad range of locations to reframe TTNtvvariants and to treat DCM patients. Additionally, to pave the way for the generation of efficient myoediting approaches for DCM, we screened and selected promising target locations in TTN. We conceptually explored the deletion of symmetric exons as a therapeutic approach to restore TTN’s reading frame in cases of frameshift mutations. We identified a set of 94 potential candidate exons of TTN that we consider particularly suitable for this therapeutic deletion. With this study, we aim to contribute to the development of new therapies to efficiently treat titinopathies and other diseases caused by mutations in genes encoding proteins with modular structures, e.g., Obscurin.

A Homozygous Deep Intronic Mutation Alters the Splicing of Nebulin Gene in a Patient With Nemaline Myopathy

Nemaline myopathy is a rare disorder affecting the muscle sarcomere. Mutations in nebulin gene (NEB) are known to be responsible for about 50% of nemaline myopathy cases. Nebulin is a giant protein which is formed integrally with the sarcomeric thin filament. This complex gene is under extensive alternative splicing giving rise to multiple isoforms. In this study, we report a 6-year-old boy presenting with general muscular weaknesses. Identification of rod-shaped structures in the patient' biopsy raised doubt about the presence of a nemaline myopathy. Next-generation sequencing was used to identify a causative mutation for the patient syndrome. A homozygous deep intronic substitution was found in the intron 144 of the NEB. The variant was predicted by in silico tools to create a new donor splice site. Molecular analysis has shown that the mutation could alter splicing events of the nebulin gene leading to a significant decrease of isoforms level. This change in the expression level of nebulin could give rise to functional consequences in the sarcomere. These results are consistent with the phenotypes observed in the patient. Such a discovery of variants in this gene will allow a better understanding of the involvement of nebulin in neuromuscular diseases and help find new treatments for the nemaline myopathy.

Is Gene-Size an Issue for the Diagnosis of Skeletal Muscle Disorders?

Human genes have a variable length. Those having a coding sequence of extraordinary length and a high number of exons were almost impossible to sequence using the traditional Sanger-based gene-by-gene approach. High-throughput sequencing has partly overcome the size-related technical issues, enabling a straightforward, rapid and relatively inexpensive analysis of large genes.

Several large genes (e.g. TTN, NEB, RYR1, DMD) are recognized as disease-causing in patients with skeletal muscle diseases. However, because of their sheer size, the clinical interpretation of variants in these genes is probably the most challenging aspect of the high-throughput genetic investigation in the field of skeletal muscle diseases.

The main aim of this review is to discuss the technical and interpretative issues related to the diagnostic investigation of large genes and to reflect upon the current state of the art and the future advancements in the field.

A long-read RNA-seq approach to identify novel transcripts of very large genes

RNA-seq is widely used for studying gene expression, but commonly used sequencing platforms produce short reads that only span up to two exon junctions per read. This makes it difficult to accurately determine the composition and phasing of exons within transcripts. Although long-read sequencing improves this issue, it is not amenable to precise quantitation, which limits its utility for differential expression studies. We used long-read isoform sequencing combined with a novel analysis approach to compare alternative splicing of large, repetitive structural genes in muscles. Analysis of muscle structural genes that produce medium (Nrap: 5 kb), large (Neb: 22 kb), and very large (Ttn: 106 kb) transcripts in cardiac muscle, and fast and slow skeletal muscles identified unannotated exons for each of these ubiquitous muscle genes. This also identified differential exon usage and phasing for these genes between the different muscle types. By mapping the in-phase transcript structures to known annotations, we also identified and quantified previously unannotated transcripts. Results were confirmed by endpoint PCR and Sanger sequencing, which revealed muscle-type-specific differential expression of these novel transcripts. The improved transcript identification and quantification shown by our approach removes previous impediments to studies aimed at quantitative differential expression of ultralong transcripts.

Nebulin nemaline myopathy recapitulated in a compound heterozygous mouse model with both a missense and a nonsense mutation in Neb

Nemaline myopathy (NM) caused by mutations in the gene encoding nebulin (NEB) accounts for at least 50% of all NM cases worldwide, representing a significant disease burden. Most NEB-NM patients have autosomal recessive disease due to a compound heterozygous genotype. Of the few murine models developed for NEB-NM, most are Neb knockout models rather than harbouring Neb mutations. Additionally, some models have a very severe phenotype that limits their application for evaluating disease progression and potential therapies. No existing murine models possess compound heterozygous Neb mutations that reflect the genotype and resulting phenotype present in most patients. We aimed to develop a murine model that more closely matched the underlying genetics of NEB-NM, which could assist elucidation of the pathogenetic mechanisms underlying the disease. Here, we have characterised a mouse strain with compound heterozygous Neb mutations; one missense (p.Tyr2303His), affecting a conserved actin-binding site and one nonsense mutation (p.Tyr935*), introducing a premature stop codon early in the protein. Our studies reveal that this compound heterozygous model, Neb^{Y2303H, Y935X}, has striking skeletal muscle pathology including nemaline bodies. In vitro whole muscle and single myofibre physiology studies also demonstrate functional perturbations. However, no reduction in lifespan was noted. Therefore, Neb^Y2303H,Y935X mice recapitulate human NEB-NM and are a much needed addition to the NEB-NM mouse model collection. The moderate phenotype also makes this an appropriate model for studying NEB-NM pathogenesis, and could potentially be suitable for testing therapeutic applications.

Nebulin: big protein with big responsibilities

Nebulin, encoded by NEB, is a giant skeletal muscle protein of about 6669 amino acids which forms an integral part of the sarcomeric thin filament. In recent years, the nebula around this protein has been largely lifted resulting in the discovery that nebulin is critical for a number of tasks in skeletal muscle. In this review, we firstly discussed nebulin’s role as a structural component of the thin filament and the Z-disk, regulating the length and the mechanical properties of the thin filament as well as providing stability to myofibrils by interacting with structural proteins within the Z-disk. Secondly, we reviewed nebulin’s involvement in the regulation of muscle contraction, cross-bridge cycling kinetics, Ca²⁺-homeostasis and excitation contraction (EC) coupling. While its role in Ca²⁺-homeostasis and EC coupling is still poorly understood, a large number of studies have helped to improve our knowledge on how nebulin affects skeletal muscle contractile mechanics. These studies suggest that nebulin affects the number of force generating actin-myosin cross-bridges and may also affect the force that each cross-bridge produces. It may exert this effect by interacting directly with actin and myosin and/or indirectly by potentially changing the localisation and function of the regulatory complex (troponin and tropomyosin). Besides unravelling the biology of nebulin, these studies are particularly helpful in understanding the patho-mechanism of myopathies caused by NEB mutations, providing knowledge which constitutes the critical first step towards the development of therapeutic interventions. Currently, effective treatments are not available, although a number of therapeutic strategies are being investigated.

Nemaline myopathies: a current view

Nemaline myopathies are a heterogenous group of congenital myopathies caused by de novo, dominantly or recessively inherited mutations in at least twelve genes. The genes encoding skeletal α-actin (ACTA1) and nebulin (NEB) are the commonest genetic cause. Most patients have congenital onset characterized by muscle weakness and hypotonia, but the spectrum of clinical phenotypes is broad, ranging from severe neonatal presentations to onset of a milder disorder in childhood. Most patients with adult onset have an autoimmune-related myopathy with a progressive course. The wide application of massively parallel sequencing methods is increasing the number of known causative genes and broadening the range of clinical phenotypes. Nemaline myopathies are identified by the presence of structures that are rod-like or ovoid in shape with electron microscopy, and with light microscopy stain red with the modified Gömöri trichrome technique. These rods or nemaline bodies are derived from Z lines (also known as Z discs or Z disks) and have a similar lattice structure and protein content. Their shape in patients with mutations in KLHL40 and LMOD3 is distinctive and can be useful for diagnosis. The number and distribution of nemaline bodies varies between fibres and different muscles but does not correlate with severity or prognosis. Additional pathological features such as caps, cores and fibre type disproportion are associated with the same genes as those known to cause the presence of rods. Animal models are advancing the understanding of the effects of various mutations in different genes and paving the way for the development of therapies, which at present only manage symptoms and are aimed at maintaining muscle strength, joint mobility, ambulation, respiration and independence in the activities of daily living.