# Exome sequencing in undiagnosed congenital myopathy reveals new genes and refines genes-phenotypes correlations

Yvan de Feraudy, Marie Vandroux, Norma Beatriz Romero, Raphaël Schneider, Safaa Saker, Anne Boland, Jean-François Deleuze, Jean-françois Deleuze, Valérie Biancalana, Johann Böhm, et al.

► To cite this version:

Yvan de Feraudy, Marie Vandroux, Norma Beatriz Romero, Raphaël Schneider, Safaa Saker, et al. Exome sequencing in undiagnosed congenital myopathy reveals new genes and refines genes-phenotypes correlations. Genome Medicine, 2024, 16 (1), pp.87. 10.1186/s13073-024-01353-0. hal-04908518

HAL Id: hal-04908518

https://hal.science/hal-04908518v1

Submitted on 23 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Exome sequencing in undiagnosed congenital myopathy reveals new genes and refines genes-phenotypes correlations

Yvan de Feraudy¹,²,³, Marie Vandroux¹, Norma Beatriz Romero⁴, Raphaël Schneider¹, Safaa Saker⁵, Anne Boland⁶, Jean-François Deleuze⁶, Valérie Biancalana¹,⁷, Johann Böhm¹† and Jocelyn Laporte¹*‡

**Abstract**

**Background** Congenital myopathies are severe genetic diseases with a strong impact on patient autonomy and often on survival. A large number of patients do not have a genetic diagnosis, precluding genetic counseling and appropriate clinical management. Our objective was to find novel pathogenic variants and genes associated with congenital myopathies and to decrease diagnostic odysseys and dead-end.

**Methods** To identify pathogenic variants and genes implicated in congenital myopathies, we established and conducted the MYOCAPTURE project from 2009 to 2018 to perform exome sequencing in a large cohort of 310 families partially excluded for the main known genes.

**Results** Pathogenic variants were identified in 156 families (50%), among which 123 families (40%) had a conclusive diagnosis. Only 44 (36%) of the resolved cases were linked to a known myopathy gene with the corresponding phenotype, while 55 (44%) were linked to pathogenic variants in a known myopathy gene with atypical signs, highlighting that most genetic diagnosis could not be anticipated based on clinical-histological assessments in this cohort. An important phenotypic and genetic heterogeneity was observed for the different genes and for the different congenital myopathy subtypes, respectively. In addition, we identified 14 new myopathy genes not previously associated with muscle diseases (20% of all diagnosed cases) that we previously reported in the literature, revealing novel pathomechanisms and potential therapeutic targets.

**Conclusions** Overall, this approach illustrates the importance of massive parallel gene sequencing as a comprehensive tool for establishing a molecular diagnosis for families with congenital myopathies. It also emphasizes the contribution of clinical data, histological findings on muscle biopsies, and the availability of DNA samples from additional family members to the diagnostic success rate. This study facilitated and accelerated the genetic diagnosis of congenital myopathies, improved health care for several patients, and opened novel perspectives for either repurposing of existing molecules or the development of novel treatments.

**Keywords** Congenital myopathy, Myopathy, Exome sequencing, Genetic diagnosis, Genetic heterogeneity, Phenotypic heterogeneity, Core myopathy, Nemaline myopathy, Centronuclear myopathy, Tubular aggregate myopathy

<footnote>†Johann Böhm and Jocelyn Laporte contributed equally to this work.</footnote>

<footnote>*Correspondence:</footnote>

<footnote>Jocelyn Laporte
jocelyn@igbmc.fr</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

Congenital myopathies (CM) are rare and severe genetic diseases strongly impacting patient autonomy and often survival. To date, many patients do not have a genetic diagnosis, precluding a better health care, the application of potential treatments, and genetic counseling including prenatal diagnosis [1–3]. Here we aimed to decrease diagnostic odysseys or dead-ends and to identify novel causative genes which will represent novel therapeutic targets.

CM are linked to muscle weakness and/or hypotonia, usually at birth or starting in early infancy, albeit some adult-onset cases were described [4–6]. Additional signs as cardiomyopathy, orthopedic deformities, or respiratory failure are also often noted. CM are not associated with muscle necrosis and regeneration and thus differ from muscular dystrophies. Instead, muscle biopsies from CM patients typically show structural anomalies characteristic for the different CM subtypes. The main subtypes are core myopathy (central core disease and multi-minicore disease), nemaline myopathy (NM) with protein aggregates forming rod-like structure, and centronuclear and myotubular myopathies with mislocalized nuclei and organelles. Less frequently, tubular aggregate, cylindrical spirals, hexagonally cross-linked tubular arrays, or other structural anomalies can be observed. Over the past 10 years, the number of causative CM genes increased from 19 to 47 genes and most code for myofilaments or proteins regulating calcium-coupled muscle contraction [7, 8].

Before 2012, the genetic diagnosis was performed on a gene-by-gene basis suggested by clinical and histopathological findings and limited by the knowledge of the implicated genes and sometimes hampered by the extreme size of some genes (TTN, RYR1, NEB). In 2012, gene panels targeting parts or all of the known neuromuscular genes were validated [9, 10]. Now, next-generation sequencing (NGS) as exome and genome sequencing is usually conducted in routine diagnosis laboratories for patients excluded for frequent genetic alterations through targeted PCR (i.e., DMPK) or enzymatic activity.

Here, to identify pathogenic variants and genes implicated in CM, we assembled a large cohort of 310 families classified into homogeneous CM subtypes based on clinical and histopathological findings. DNA samples from probands and relatives were analyzed by exome sequencing and in-house bioinformatics tools and candidate variants/genes were confirmed through genetic and functional investigations. This study is part of the French MYOCAPTURE NGS project aiming to identify novel genes mutated in myopathies by sequencing the exome of > 1000 individuals from families with different forms of myopathies.

## Methods

### Patient recruitment and evaluation

The present project focused on congenital myopathies. Patients were examined by clinicians from the French neuromuscular network FILNEMUS (www.filnemus.fr) or from their country of origin. All patients enrolling the study manifested either neonatal hypotonia or early- or late-onset muscle weakness, and most underwent muscle biopsies with histological analysis. Most biopsies were ascertained through the Morphology Unit of the Myology Institute (Paris, France), which receives requests from all over France, as well as from South America, while the remaining biopsies were analyzed by histologists from the patient's home countries. Patients were classified into homogeneous CM cohorts based on the clinical severity and histological findings on the muscle biopsies (Table 1). To enhance cohort homogeneity, the majority of histological classification was conducted by a single histopathologist (NBR, Myology Institute, Paris). Clinical, histological, and genetic data of patients are included in Table S1, where associated publications are mentioned when previously published.

For most patients, the implication of the main known genes at the time of inclusion (2009–2018) was excluded mainly through targeted gene sequencing or rarely through myopathy-related gene panel. For RYR1, the 5' sequence was barely tested, while RT-PCR was occasionally used. Based on these criteria, we collected DNA samples from 310 families, including 429 patients and 459 unaffected relatives, representing a varied ethnic population.

### Exome sequencing and data analysis

Library and exome capture of patients were performed with the SureSelect Human all exon kit (Agilent v2 or v5, Santa Clara, CA) and samples were paired-end sequenced on a HiSeq2000 (Illumina, San Diego, CA). Sequences were obtained from the CNRGH (Evry, France), the Genomeast platform (Illkirch-Graffenstaden, France), or BGI (Shenzhen, China). Sequence data were aligned to the reference genome GRCh37/hg19. Analysis of NGS data and variants calling was performed as previously described [11]. We applied segregation scenarios based on the family structure (recessive homozygous/compound heterozygous, dominant, de novo, X-linked, and other specific scenarios). Filtering of variants was also based on population frequency in both the gnomAD database and an internal NGS database. Then, we used our VaRank algorithm to weight and prioritize variant [12]. For missense variants, we used in silico predictions of pathogenicity with SIFT [13] or PolyPhen[14]. For splice variants, we used splice variant predictor tools with MaxEntScan [15], NNSplice [16], and Splice Site
<table><caption>**Table 1** Patients cohorts based on clinical and histological findings. Main histological features on muscle biopsies and diagnostic success rate per group</caption><thead><tr><th>Patient cohorts</th><th>Total number of cases: no. (%)</th><th>Number of cases with identified variants: no. (%)</th></tr></thead><tbody><tr><td>Centronuclear myopathy</td><td>82 (26)</td><td>43 (52)</td></tr><tr><td>Core myopathy</td><td>54 (17)</td><td>21 (39)</td></tr><tr><td>Nemaline myopathy</td><td>41 (13)</td><td>20 (49)</td></tr><tr><td>Tubular aggregate myopathy</td><td>34 (10)</td><td>8 (24)</td></tr><tr><td>Myopathy with protein aggregates</td><td>13 (4)</td><td>3 (23)</td></tr><tr><td>Core-rod myopathy</td><td>8 (2)</td><td>6 (75)</td></tr><tr><td>Vacuolar myopathy</td><td>7 (2)</td><td>5 (71)</td></tr><tr><td>Cylindrical spirals</td><td>4 (1)</td><td>0 (0)</td></tr><tr><td>Protein aggregates and rimmed vacuoles</td><td>3 (1)</td><td>2 (67)</td></tr><tr><td>Cap myopathy</td><td>3 (1)</td><td>1 (33)</td></tr><tr><td>Cytoplasmic bodies</td><td>3 (1)</td><td>1 (33)</td></tr><tr><td>Abnormal nuclear envelop</td><td>3 (1)</td><td>0 (0)</td></tr><tr><td>Mini-rods myopathy</td><td>2 (1)</td><td>1 (50)</td></tr><tr><td>Congenital fiber type disproportion</td><td>1 (1)</td><td>1 (100)</td></tr><tr><td>Broad A band disease</td><td>1 (1)</td><td>1 (100)</td></tr><tr><td>Myofibrillar myopathy</td><td>1 (1)</td><td>0 (0)</td></tr><tr><td>Hexagonally cross-linked crystalloid inclusions</td><td>1 (1)</td><td>0 (0)</td></tr><tr><td>Dark inclusion</td><td>1 (1)</td><td>0 (0)</td></tr><tr><td>Undefined CM</td><td>19 (6)</td><td>6 (32)</td></tr><tr><td>No biopsy</td><td>29 (9)</td><td>4 (14)</td></tr><tr><td>Total</td><td>310 (100)</td><td>123 (40)</td></tr></tbody></table>

Finder [17]. For a subset of synonymous, canonical or non-canonical splice site variants, RT-PCR from muscle RNA with subsequent cDNA sequencing was used to attest a pathogenic effect on exon splicing. For certain variants, a segregation study was required to determine their pathogenicity. In the absence of DNA samples from relatives, these variants were classified as candidates. Variants were classified according to the ACMG standards and guidelines [18].

## Other investigations

In case of homozygous variants, familial consanguinity was assessed using Somalier, measuring the relatedness between two DNA samples by comparing informative genetic variations [19]. When a pathogenic variant was confirmed in a known gene, the clinical and histological phenotype of the patient was compared to those already described in the OMIM database to determine if it corresponded to a novel phenotype. For novel genes, functional validation of the pathogenicity of the variants was performed using in vitro experiments to attest their impact on RNA stability, protein level, and protein functions, while in vivo models were used to demonstrate they could reproduce patient phenotypes. Protocols used

and results obtained have been detailed and published in the associated papers [20–29].

## Ethics

The present study was approved by the IRB Comité de Protection des Personnes Est IV (DC-2012-1693), and all participants or their legal guardians provided written informed consents to the different clinical centers involved in the project.

## Results

### Description of the cohort

A total of 310 families with genetically undiagnosed CM enrolled this study from 2009 to 2018. For 214 families, at least one relative was available, and for 92 families, a complete trio composed of the index patient and both parents was available. Overall, 429 individuals were analyzed by exome, or Sanger sequencing for segregation study. Most probands were sporadic cases (78%) with no other affected family member or known family history of a neuromuscular disorder, while familial cases with apparent X-linked, autosomal dominant, or autosomal recessive disease inheritance represented 22%. Parents were reported to be consanguineous in 25 families.
A majority of 61% of the probands was male, while 39% were female. In this international multicentric study, 175 probands were from France (57%), 102 probands (32%) from 15 countries encompassing Algeria, Argentina, Australia, Belgium, Brazil, Canada, Finland, Germany, Israel, Italy, Luxembourg, New Zealand, Spain, Turkey, Uruguay, and 32 probands from unknown origin (11%) (Fig. 1A).

Age of onset was known for 99 probands (32%), with 58 probands manifesting antenatal or neonatal disease signs, 27 with childhood/adolescence onset, and 14 with onset in adulthood. This last category had the highest proportion of families with more than one affected individual. Muscle biopsies were performed and morphologically examined for 281 probands (91%). Based on the types of histological anomalies, families were classified into homogeneous cohorts (Table 1). Centronuclear myopathy (CNM) was the largest cohort (n=82 cases, 26%), followed by core myopathy (n=54, 17%), NM (n=41, 13%), and tubular aggregate myopathy (TAM) (n=34, 10%). This distribution most probably reflects an inclusion bias as our laboratory has a known expertise in CNM and TAM.

## Overall genetic characterization of the cohort

Exome sequencing identified causal pathogenic variants for 123 families (40%) (Fig. 1B). A comparative overview

of the clinical, histological, and genetic findings of all families with confirmed molecular diagnosis is given in Table S1. For 33 families (10%), we detected sequence variants in convincing candidate genes, but DNA samples from relatives were not available, precluding segregation analyses. In our cohort, there was a correlation between the number of sequenced family members and the diagnosis rate. Increasing the number of sequenced family members enhanced the identification of the pathogenic variants, with a majority of families (51%) being diagnosed when at least three DNA samples were available (Fig S1). More specifically, the proportion of diagnosed families was higher for trios (48%) compared to duos (29%). Similarly, the proportion of families with identified pathogenic variants increased when the number of affected individuals per families increased, with a majority of families (55%) being diagnosed when at least two affected individuals were available (Fig S2). Muscle biopsies with histological examination were another major factor contributing to molecular diagnosis. Indeed, causative variants were identified in 42% of the families with biopsy and only in 14% of the families without biopsy (Fig S3). This may be partially due to the fact that patients with the most severe clinical phenotypes are more often biopsied; however, these results support the importance of complementary investigations in addition to clinical assessment, as histopathology (or MRI, not implemented

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Overview of patients enrolling the MYOCAPTURE cohort. A Geographic origin of genetic samples. Europe: European countries other than France. B Diagnosis ratio. Diagnosis: patients with confirmed pathogenic variant(s). Candidate: patients with suspected pathogenic variant(s) in one candidate gene</figcaption></figure>
heterozygous, hemizygous, or compound heterozygous pathogenic variants carried de novo pathogenic variants (15%). This knowledge should be useful for genetic information of relatives. Patients from the 16 diagnosed consanguineous families did not always exhibit homozygous pathogenic variants, as one proband had a dominant pathogenic variant in the *STIM1* gene (family 105). This example confirms that the occurrence of de novo mutations is equally probable in consanguineous and non-consanguineous families and that consanguineous families should not only be screened for homozygous variants. Conversely, consanguinity was revealed in one family with homozygous pathogenic variant through the use of the Somalier software (family 33). In total and across all patient cohorts, we identified 202 pathogenic variants (Table S1) with a majority of missense variants (53%) and frameshift variants (23%) (Fig. 2B). The majority of the pathogenic variants are unique. Only 12 (6%) are shared by a maximum of 2 patients each, with 9 of them involving the *RYR1* gene. Thus, in this cohort, there was no strong mutation hotspot or founder effect. Finally, two probands (family 8 and 49) displaying different histological presentations (respectively NM and cap myopathy) were, in fact, third-degree relatives. They had pathogenic variants in different genes (respectively in *ACTA1* and *MYPN*). This highlights the possibility of the occurrence of two distinct and unrelated myopathies within large families.

## Main congenital myopathy genes and impact of exome sequencing

Among the 123 families with confirmed molecular diagnosis, 47 different causative genes were identified. Strikingly, 4 genes were over-represented with diseased-associated variants in 60 of these families (49%). Specifically, the *RYR1* gene was implicated in the highest number of cases (n=29 cases, 24%), followed by *NEB* (n=15, 12%), *TTN* (n=8, 7%), and *ACTA1* (n=8, 7%) (Fig. 3).

Sixteen of the 123 diagnosed families (13%) carried pathogenic variants that may easily escape routine diagnosis. Seven pathogenic variants were identified in non-canonical splice sites in *MTM1* (n=1), *RYR1* (n=1), *UNC45B* (n=1), and *NEB* (n=4). Five synonymous *NEB*, *RYR1*, and *TTN* variants were found to impact on splicing in three families with NM, one family with core myopathy and one family with CNM, respectively. Three families with NM had mosaic pathogenic variants in the *ACTA1* gene that was detectable on muscle DNA and only barely detectable on lymphocyte DNA [30]. Finally, heterozygous pathogenic *MTM1* variants were detected in two symptomatic female carriers with CNM, and one *RYR1* missense variant previously considered as polymorphism has now been reclassified as a pathogenic variant in a CNM family [31]. Overall, these examples emphasize the importance to thoroughly examine non-coding regions, use prediction programs to evaluate

<figure><img src="image_4.png" /><figcaption>**Fig. 3** Distribution of known and novel myopathy genes. In the bar chart, the percentages indicate the ratio of families carrying pathogenic variant(s) in a given gene. New genes identified within the MYOCAPTURE project are shown in orange. For known genes, percentages of patients with classical phenotype (blue) or with new phenotype (yellow) are indicated. The global distribution of pathogenic variants is depicted in the upper diagram</figcaption></figure>
potential pathogenic effects of intronic and synonymous variants, perform high-throughput sequencing on DNA from the affected tissue or adjust the bioinformatics pipeline filters to detect variants in a low number of reads, and consider X-linked disorders also in females.

Among the 123 diagnosed families, 11 (9%) had pathogenic variants that were not detected by previous myopathy-related gene panels or targeted sequencing of known CM genes. This was essentially due to the presence of pathogenic variants in the large *RYR1* gene (*n*=8), which was previously only sequenced for the most 3'-terminal exons, or to mosaic pathogenic variants in *ACTA1* (*n*=3), which were not detected by routine genetic tests on lymphocyte DNA. As a result, the use of exome sequencing on muscle DNA became necessary for one of these patients to identify the causative gene.

## Known myopathy genes linked to classical phenotypes

For 36% of the families (*n*=44) with confirmed molecular diagnosis, we identified pathogenic variants in known CM genes concurring with the clinical and histopathological characteristics. The majority of these cases encompassed pathogenic *RYR1* variants in core myopathy, *NEB* and *ACTA1* in nemaline myopathy, and pathogenic *MTM1* variants in CNM (Fig. 3). As examples, eight patients with core myopathy had pathogenic variants in *RYR1*, eight NM cases harbored the typical compound heterozygous combination of a nonsense or a frameshift variant with a splice site variant affecting an in-frame exon of *NEB*. We also detected heterozygous missense variants in the highly conserved *ACTA1* gene in five patients with classical NM and nemaline rods on muscle sections. Also, one CNM case with severe neonatal hypotonia and the presence of abnormally centralized nuclei in the myofiber carried heterozygous missense variant affecting the PH domain of DNM2 [32].

## Genes linked to atypical phenotypes

In addition to pathogenic variants in known CM genes associated with a classical phenotype, we also found that the majority (51%) of pathogenic variants were located in previously known neuromuscular genes, including 44% in patients displaying atypical phenotypes at the time of inclusion (Fig. 3). They were divided into three categories:

Pathogenic variants in known genes associated with a different muscle disorder: importantly, 46 of the 123 diagnosed families (37%) had pathogenic variants in genes previously associated with other forms of muscular diseases, with 40 families showing pathogenic variants in genes associated with other forms of myopathy (Table 2). Among these, the most frequently identified genes were *RYR1* in patients with CNM

(n=17) [33, 34] or dusty core disease (n=2) [35], and *TTN* (n=6) in patients with CNM [36]. These novel genotype-phenotype associations were reported by us or others during the course of this study. The remaining six families had pathogenic variants in genes usually associated with muscular dystrophies (*HNRNPDL* [37], *TOR1AIP1* [38], *SGCG*, *COL6A2*, *CAV3*, *TRIP4* [39, 40]). These variants were identical or affecting the same residues as described in previously reported dystrophy cases, highlighting the genetic and clinical heterogeneity of CMs.

Pathogenic variants in known genes associated with other neuromuscular or cardiac diseases: interestingly, three families had novel pathogenic variants in genes usually associated with neuropathies: *ASCC1* (n=2) associated with a severe myopathy with arthrogryposis and bone fractures [41] and *HSPB8* (n=1) linked to a vacuolar adult-onset myopathy [42]. In the latter case, electromyography showed no evidence of an underlying neuronopathy or polyneuropathy. In addition, four families with muscle biopsies showing characteristic intermyofibrillar network and centralized nuclei or core-like structures had pathogenic variants in *CACNA1S*, a gene usually associated with hypokalemic periodic paralysis or malignant hyperthermia [11]. Moreover, four families with CM carried pathogenic variants in *GFPT1*, *CHRNA1*, or *CHRND*, all three previously associated with congenital myasthenia often involving muscle weakness similar to CM. In contrast to the reported *GFPT1*, *CHRNA1*, and *CHRND*-related myasthenia families, the muscle biopsies from our patients displayed either protein aggregates and rimmed vacuoles or CNM-typical nuclear centralization. Lastly, pathogenic variants in *MYPN* previously described in patients with cardiomyopathy was found in two patients with nemaline and cap myopathies [43].

Pathogenic variants in known genes associated with other syndromes: four families presented a pathogenic variant in a gene not previously known to cause specific neuromuscular or cardiac diseases. In one family, we detected a known homozygous pathogenic variant in *IDUA*, associated with a mild form of mucopolysaccharidosis known as Scheie syndrome [44]. In addition to Scheie-typical features as glaucoma, clouded cornea, carpal tunnel syndrome, and aortic valve disease, our patients also presented with distal retractions and cores on the muscle biopsies. In a second family with CNM, a homozygous pathogenic variant was found in the *OPA1* gene, which encodes a mitochondrial protein previously linked to optic atrophy. In the third family, two deceased fetuses diagnosed with CNM had a hemizygous
<table><caption>**Table 2** Genes linked to atypical phenotypes. Genes and associated phenotypes, according to the OMIM database, are indicated in the first two columns. Starting from the first row, the genes are grouped into myopathy/ muscle dystrophy genes (n = 18), cardiomyopathy genes (n = 1), congenital myasthenia genes (n = 3), channelopathy genes (n = 1), neuropathy genes (n = 2), and genes associated with other syndromes (n = 4). The new phenotypes found in the MYOCAPTURE cohort are indicated in the third column. EM electron microscopy</caption><thead><tr><th>Gene</th><th>Usual phenotype(s) (OMIM)</th><th>New phenotype(s)</th></tr></thead><tbody><tr><td colspan="3">Congenital myopathy or muscle dystrophy</td></tr><tr><td>ACTA1</td><td>Nemaline myopathy</td><td>Nemaline myopathy with enlarged perinuclear space (EM)</td></tr><tr><td>CAV3</td><td>Distal myopathy</td><td>Tubular aggregates myopathy</td></tr><tr><td>CFL2</td><td>Nemaline myopathy</td><td>Centronuclear myopathy</td></tr><tr><td>COL6A2</td><td>Bethlem myopathy/Ulrich congenital muscular dystrophy</td><td>Cores myopathy</td></tr><tr><td>FLNC</td><td>Distal myopathy/myofibrillar myopathy</td><td>Nemaline myopathy</td></tr><tr><td>HNRNPA2B1</td><td>Inclusion body myopathy with early-onset Paget disease with or without frontotemporal dementia 2</td><td>Early onset form of oculopharyngeal muscular dystrophy</td></tr><tr><td>HNRNPDL</td><td>Muscular dystrophy, limb-girdle</td><td>Autophagic rimmed vacuolar myopathies</td></tr><tr><td>TPM3</td><td>CAP myopathy/congenital myopathy with fiber type disproportion/nemaline myopathy</td><td>Cores myopathy / Centronuclear myopathy / Caps and rods myopathy</td></tr><tr><td>FHL1</td><td>Emery-Dreifuss muscular dystrophy/reducing body myopathy/myopathy, X-linked, with postural atrophy/scapuloponeal myopathy</td><td>Centronuclear myopathy</td></tr><tr><td>MTM1</td><td>Centronuclear myopathy</td><td>Congenital fiber-type disproportion</td></tr><tr><td>MYH7</td><td>Laing distal myopathy/myosin storage myopathy/myopathic scapuloponeal syndrome</td><td>Central core myopathy</td></tr><tr><td>RYR1</td><td>Congenital myopathy with susceptibility to malignant hyperthermia/core myopathy/King-Denborough syndrome</td><td>Centronuclear myopathy / Dusty core disease</td></tr><tr><td>SGCG</td><td>Muscular dystrophy, limb-girdle, autosomal recessive 5</td><td>Myopathy with broad A band disease</td></tr><tr><td>SPEG</td><td>Centronuclear myopathy</td><td>Myopathy with no centralized nuclei</td></tr><tr><td>TOR1AIP1</td><td>Muscular dystrophy, autosomal recessive, with rigid spine and distal joint contractures</td><td>Centronuclear myopathy</td></tr><tr><td>TRIP4</td><td>Congenital muscular dystrophy, Davignon-Chauveau type/spinal muscular atrophy with congenital bones fractures 1</td><td>Centronuclear myopathy with cores</td></tr><tr><td>TTN</td><td>Muscular dystrophy, limb-girdle/myofibrillar myopathy/Salih myopathy/tibial muscular dystrophy tardive</td><td>Centronuclear myopathy</td></tr><tr><td colspan="3">Cardiomyopathy</td></tr><tr><td>MYPN</td><td>Cardiomyopathy</td><td>Caps myopathy</td></tr><tr><td colspan="3">Congenital myasthenia</td></tr><tr><td>CHRNA1</td><td>Congenital myasthenic syndrome type 1A or 1B</td><td>Centronuclear myopathy</td></tr><tr><td>CHRNĐ</td><td>Congenital myasthenic syndrome type 3A, 3B or 3C</td><td>Centronuclear myopathy</td></tr><tr><td>GFPT1</td><td>Congenital myasthenia with tubular aggregates</td><td>Rimmed vacuole myopathy</td></tr><tr><td colspan="3">Channelopathy</td></tr><tr><td>CACNA1S</td><td>Hypokalemic periodic paralysis type 1</td><td>Myopathy with myofibrillar disorganization</td></tr><tr><td colspan="3">Neuropathy</td></tr><tr><td>ASCC1</td><td>Spinal muscular atrophy with congenital bone fractures</td><td>Congenital myopathy</td></tr><tr><td>HSPB8</td><td>CMT type 2L/neuronopathy, distal hereditary, type IIA</td><td>Vacuolar myopathy</td></tr><tr><td colspan="3">Other syndromes</td></tr><tr><td>EXOSC3</td><td>Short stature, hearing loss, retinitis pigmentosa, and distinctive facies (SHRF)</td><td>Myopathy with caps and mini-rods</td></tr><tr><td>FOXP3</td><td>Immunodysregulation, polyendocrinopathy, and enteropathy, X-linked</td><td>Centronuclear myopathy</td></tr><tr><td>IDUA</td><td>Mucopolysaccharidosis Ih, Ih/s, Is</td><td>Central core myopathy</td></tr><tr><td>OPA1</td><td>Behr syndrome/optic atrophy/optic atrophy plus syndrome</td><td>Centronuclear myopathy</td></tr></tbody></table>
<table><caption>**Table 3** New genes found implicated in myopathies. New genes and associated phenotypes are indicated in the first three columns. The references for the characterization of the pathogenic variants found in the MYOCAPTURE cohort are indicated in the fourth column</caption><thead><tr><td>New genes</td><td>Clinical description</td><td>Histological phenotype</td><td>Reference(s)</td></tr></thead><tbody><tr><td>STIM1</td><td>Generalized slowly progressive muscle weakness with variable age of onset, Stormorken syndrome</td><td>Tubular aggregate myopathy</td><td>[20, 46, 47]</td></tr><tr><td>ORAI1</td><td>Generalized slowly progressive muscle weakness with variable age of onset</td><td>Tubular aggregate myopathy</td><td>[21, 48, 49]</td></tr><tr><td>CASQ1</td><td>Generalized slowly progressive muscle weakness, post-exercise myalgia</td><td>Tubular aggregate myopathy</td><td>[22, 50]</td></tr><tr><td>PYROXD1</td><td>Generalized slowly progressive muscle weakness</td><td>Myopathy with internalized nuclei and myofibrillar disorganization</td><td>[23, 51]</td></tr><tr><td>MAP3K20</td><td>Generalized slowly progressive muscle weakness accompanied by decreased vital capacities</td><td>Centronuclear myopathy</td><td>[24]</td></tr><tr><td>UNC45B</td><td>Generalized slowly progressive muscle weakness accompanied by decreased vital capacities</td><td>Myopathy with eccentric cores, internalized nuclei and myofibrillar disorganization</td><td>[25]</td></tr><tr><td>MYO18B</td><td>Congenital myopathy with cardiomyopathy and dysmorphic features</td><td>Nemaline myopathy</td><td>[26, 52]</td></tr><tr><td>ACTN2</td><td>Generalized slowly progressive muscle weakness accompanied by decreased vital capacities</td><td>Core myopathy with jagged Z-lines</td><td>[27]</td></tr><tr><td>GGPS1</td><td>Muscular dystrophy with hearing loss and ovarian insufficiency</td><td>Muscular dystrophy</td><td>[28]</td></tr><tr><td>SRPK3/TTN<sup>a</sup></td><td>Generalized slowly progressive muscle weakness with cardiomyopathy</td><td>Core myopathy</td><td>[29]</td></tr></tbody></table>

pathogenic variant in the *FOXP3* gene linked to a syndrome with immunodysregulation, polyendocrinopathy, and enteropathy. In the last family, one patient with a pathogenic *EXOSC3* variant presented with pontocerebellar hypoplasia. The muscle biopsy showed not only marked fiber type grouping, suggesting a neurogenic process, but also the presence of cap-like structures and mini-rods [45]. Overall, except for the two *FOXP3* fetuses, the phenotypes of the listed families conformed with the identified pathogenic variants and genes, and our investigations additionally uncovered an abnormal skeletal muscle structure that may have been overlooked in the previously reported families.

## New genes implicated in myopathies

As discussed above, we reported the first pathogenic variants causing a primary myopathy in four genes previously linked to other neuromuscular or cardiac diseases: *ASCC1*, *HSPB8*, *CACNA1S*, and *MYPN*. Moreover, in the MYOCAPTURE project, pathogenic variants were found in 10 new genes not previously linked to a neuromuscular genetic disease in 16 families, including *ACTN2*, *CASQ1*, *GGPS1*, *MAP3K20/ZAK*, *ORAI1*, *MYO18B*, *PYROXD1*, *SRPK3*, *STIM1*, and *UNC45B* (Tables 2 and 3). The genetic and functional validations evidencing the pathogenic impact of the identified variants in cell and

animal models were published previously. An overview is provided in Table 3.

In three families with TAM, we found dominant missense variants in *STIM1*, encoding a reticular Ca²⁺ sensor [20, 46, 47, 53]. Together with the plasma membrane Ca²⁺ channel ORAI1, STIM1 forms the key element of the ubiquitous store-operated Ca²⁺ entry (SOCE) pathway, an essential mechanism regulating Ca²⁺ homeostasis in all cell types. The pathogenic *STIM1* variants affected conserved amino acids in the Ca²⁺-binding EF-hands, and functional experiments evidenced a gain-of-function effect leading to SOCE overactivation and excessive extracellular Ca²⁺ entry. As second and third TAM genes, we and others identified pathogenic variants in *ORAI1* [21, 47–49] and in the muscle-specific reticular Ca²⁺ buffer calsequestrin (*CASQ1*) [22, 50].

In two families, pathogenic variants were found in the *PYROXD1* gene, resulting in CM with muscle fibers exhibiting multiple internal nuclei and cores, along with myofibrillar disorganization [23, 51]. PYROXD1 acts as an oxidoreductase, and tests in the yeast model have shown that the identified variants impair the enzymatic activity and enhance the sensitivity of the cells to oxidative stress. For one family, pathogenic *MAP3K20/ZAK* variants were associated with CNM and rimmed sarcolemma, implicating a novel kinase of unknown function in muscle [24]. *UNC45B* encodes a myosin-directed chaperone essential for sarcomeric organization and

<footnote><sup>a</sup> Digenism</footnote>
muscle function. Pathogenic variants in this gene were found in one family with CM, presenting a muscle biopsy showing internalized nuclei and myofibrillar disorganization [25]. In a single family with severe NM and cardiomyopathy, we detected a homozygous truncating variant in *MYO18B*, coding for an unconventional myosin [26, 52]. In one family, a pathogenic *ACTN2* variant induced core myopathy with jagged z-lines. *ACTN2* encodes alpha-actinin 2, a known regulator of the Z-line and sarcomeric structure [27]. In addition, a pathogenic *GGPS1* variant was identified in one family with muscle dystrophy, hearing loss, and ovarian insufficiency syndrome, suggesting an impairment in geranylgeranylation of small GTPases with multisystemic effects [28]. Lastly, digenic inheritance with co-occurrence of heterozygous truncating variants in *TTN* and on the X-chromosomal *SRPK3* gene, encoding a protein kinase, was found in three families with core myopathy [29].

## Widening the genetic and clinical heterogeneity

Through genetic analysis of this cohort, we were able to widen the known genetic and clinical-histological heterogeneity of CM. Regarding the phenotypic heterogeneity associated with a specific gene, we observed a genetic overlap among the core, nemaline, and CNM cohorts (*CACNA1S*, *RYR1*, *TPM3*, and *TTN*), whereas less common structural anomalies (e.g., TAM) tend to involve other genes (Fig. 4). Such overlap appeared mainly due to the fact that single patient depicted a mix of several structural defects previously associated with distinct cohorts (e.g., cores with mislocalized nuclei linked to pathogenic *RYR1* or *CACNA1S* variants), rather than different patients mutated in the same gene displaying different and specific histopathological hallmarks (e.g., only cores or only central nuclei). Particularly noteworthy is the observation that patients with pathogenic variants in *TTN*, *ACTA1*, *NEB*, or *TPM3* exhibit the widest diversity of histopathological hallmarks in muscle biopsies.

The genetic heterogeneity, i.e., several genes linked to a similar phenotype, is mainly due to (1) the identification of novel myopathy genes and (2) the identification of novel genotype-phenotype associations. For example, at the beginning of the MYOCAPTURE project in 2009, no gene was associated with TAM. This project and other studies identified *STIM1*, *ORAI1*, and *CASQ1* as causative genes [20–22, 48–50, 53]. Another example is the CNM cohort where many genes have been linked, probably in part due to the presence of mislocalized nuclei in the biopsy of patients.

All these findings should be valuable for geneticists, histopathologists, and clinicians for the future diagnosis of different forms of myopathies.

## Discussion

The primary goal of this MYOCAPTURE project was the identification of novel genes associated with CM and the shortening of diagnostic odysseys and dead-ends. It was mainly based on exome sequencing of 310 families with genetically undiagnosed CM, classified in homo geneous cohorts based on clinical and histopathological findings. Pathogenic variants were identified in 50% of the families, of which 79% obtained a definite diagnosis. For the remaining 21%, we identified sequence variants with likely pathogenic effect in muscle genes. Strik ingly, only 36% of the patients with resolved diagnosis were linked to a known myopathy gene with a concord ant phenotype, pointing out that clinical and histological assessments were not sufficient to anticipate molecular diagnosis in the majority of cases. This included genes previously associated with other myopathies, neuromus cular diseases, or other syndromic disorders. Moreover, the MYOCAPTURE project highlighted the wide genetic and phenotypic heterogeneity characterizing congenital myopathies and successfully identified 14 novel genes not previously associated with muscle diseases.

## Impact of massive parallel sequencing on the diagnosis of congenital myopathies

Patients included in the MYOCAPTURE cohort from 2009 to 2018 were devoid of a genetic diagnosis despite more or less extensive genetic studies. Using exome sequencing, we identified causative variants in 50% of the cases. To confirm the implication of the identified vari ants and genes, we performed complementary analyses such as segregation studies on DNA samples from rela tives or sequencing of muscle RNA. However, a subset of patients was singletons without available parental DNA for segregation studies or biological material for RNA extraction, so that a definite molecular diagno sis could not be delivered. This success rate is similar to other studies using large gene panels or exome sequenc ing in muscle-related diseases as limb-girdle muscular dystrophies, myopathies, or arthrogryposis multiplex congenital [54–57]. We demonstrate the massive paral lel analysis of all coding genes through exome sequenc ing is more efficient than previous genetic investigations on a gene-by-gene basis or gene panels to pinpoint the genetic defect. Pathogenic variants in *RYR1*, *TTN*, *NEB*, and *ACTA1* were prevalent in our CM cohort, and sev eral reasons account for this. First of all, all four are major CM genes coding for essential proteins of muscle func tion and structural maintenance, and pathogenic variants are likely to interfere with normal myofiber contractil ity. Second, the first three belong to the largest genes of the human genome and have a higher probability to
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Congenital myopathy-causing genes. In total, pathogenic variants in 48 genes have been identified in the MYOCAPTURE cohort. The diagram depicts the distribution of these genes based on histological phenotypes and the intersections between phenotypes. The new genes not previously linked to a congenital myopathy are indicated in bold and underlined. CFTD congenital fiber type disorder, CNM centronuclear myopathy, TAM tubular aggregate myopathy</figcaption></figure>

To our own surprise, only 36% of the families with resolved molecular diagnosis harbored pathogenic variants in known myopathy genes corresponding to their clinical presentation, highlighting a strong genetic and clinical-histological heterogeneity. Conversely, 64% of the pathogenic variants would have been missed with diagnostic panels targeting only a distinct subtype of myopathies. Larger panels encompassing all genes implicated in congenital myopathies or neuromuscular disorders would presumably have found these variants, but they are inapt to discover novel genes. Within the MYOCAPTURE project, we found 14 novel genes (20%

of resolved cases) not previously associated with muscle disorders and thereby provided a molecular diagnosis to patients and families that remained without conclusive result after gene panel sequencing. Nonetheless, exome sequencing does not cover most intronic and intergenic regions and only partially and indirectly uncovers structural anomalies as large deletions, inversions, or chromosomal rearrangements impacting on gene expression or mRNA splicing. To overcome these limitations and to ensure a uniform coverage of all coding and non-coding areas, genome sequencing in combination with muscle RNA-seq, and/or other omics approaches can be used [58]. Also, we sequenced DNA extracted mainly from the blood, which precluded in most cases the identification of mosaic mutations present only in a few tissues including muscles. Another limitation of this study is the lack of MRI data that may have allowed a better diagnosis orientation [59, 60]. For example, characteristic involvement patterns of specific muscles and muscle groups have been documented in NM, enabling the differentiation of NEB-, ACTA1-, and TPM3-related NM [61].

Despite all limitations, the MYOCAPTURE project resolved a large number of undiagnosed cases and provided a significant help to the affected families. Indeed, numerous patients obtained the clinical diagnosis of a
myopathy several years and decades ago, but still awaited
molecular diagnosis. As an example, a patient with reces-
sive *RYR1* variants underwent extensive genetic tests
and gene-by-gene investigations over 15 years without
conclusive result, and another patient with pathogenic
*CACNA1S* variant was clinically followed over 60 years
before enrolling our exome sequencing program (family
1 in Schartner et al. [11]). The knowledge of the causa-
tive gene enables an adapted genetic counseling for the
patients and his or her relatives and also allows a prog-
nosis on disease development in certain cases. It can
also suggest complementary examinations and improve
clinical follow-up as heart monitoring in *MYH7* or *TTN*-
related congenital myopathies or as adaptation of anes-
thetics for *RYR1* patients. Notably, two unrelated families
(family 72 and 76) had a triplet of *RYR1* variants in cis,
which was previously linked to malignant hyperthermia
susceptibility [62]. In some cases, treatments such as
acetylcholinesterase inhibitors may be considered, par-
ticularly in patients with pathogenic variants in congeni-
tal myasthenia genes (*CHRNA1*, *CHRND*, *DPAGT1*, and
*GFPT1*). Moreover, patients with resolved genetic diag-
nosis can now be included into clinical trials.

Refinement of myopathy classification reveals extensive phenotypic and genetic heterogeneity

We found a genetic diagnosis linking a typical pheno-
type with a gene previously associated to a myopathy in
only 36% of resolved cases. The other two-third of cases
included genes previously linked to other myopathies
(37%; Fig. 3), to congenital myasthenia syndrome (3.5%),
to other syndromes (3.5%), or to novel genes (20%). *TTN*,
*RYR1*, and *TPM3* are the genes linked to the highest
diversity of phenotypes. Overall, we cannot exclude an
inclusion bias leading to either the enrichment of patients
with a rather variable spectrum of clinical and histologi-
cal features in some cohorts or to cohorts with less strin-
gent classification criteria. Most biopsies were examined
by the Neuromuscular Morphology Unit of the Institute
of Myology in Paris to prevent high heterogeneity in his-
topathological classification. However, other centers have
also included patients in the MYOCAPTURE project and
may not have followed the same inclusion criteria. For
example, a large number of patients enrolling the MYO-
CAPTURE project exhibited nuclear mislocalization on
muscle biopsies and were classified as CNM by the refer-
ring clinicians. Although CNM is genetically heteroge-
neous with six causative genes reported to date (*BIN1*,
*DNM2*, *MTM1*, *RYR1*, *SPEG*, *TTN*), exome sequencing
identified pathogenic variants in a markedly higher diver-
sity of genes. Some histopathologists may define CNM by
the presence of a large number of hypotrophic myofibers

with central nuclei, while others may consider a few
internal nuclei in normally sized myofibers [33, 34, 63].

In some cases, additional investigations might have
been required to fully validate the classification of
patients into a specific cohort at the timepoint of inclu-
sion into the MYOCAPTURE project. By way of example,
electron microscopy was not performed for all patients
(for example to validate the classification of TAM by the
detection of tubular membrane aggregates). Indeed, fol-
lowing the genetic diagnosis, some patients could have
been re-classified into a different phenotypic cohort.

A major finding was an overlap of histopathological
hallmarks in a subset of patients, as the presence of cores
with mis-localized nuclei and myofibrillar disorganiza-
tion for *CACNA1S* patients [11], core-like area, central
nuclei and myofibrillar disorganization for *TTN* [36], or
nemaline rods, cap structures, and internal nuclei for
*TPM3* [64].

Conversely, other patients displayed no specific struc-
tural defects on the biopsy. Noteworthy, neonates may
not fully recapitulate the structural hallmarks at birth as
some develop or change with age, challenging the histo-
pathological classification. Lastly, this study also focused
on less explored cohorts, as TAM, leading to the dis-
covery of novel genes. In conclusion, exome sequenc-
ing of our large cohort encompassing 310 families
revealed a larger than expected phenotypic and genetic
heterogeneity.

Novel myopathy genes suggesting new pathomechanisms
and therapeutic targets

Through this exome sequencing project, we identified
ten new myopathy genes not previously linked to any
neuromuscular diseases and four new myopathy genes
that were previously associated with other neuromuscu-
lar or cardiac diseases. The validation of the implication
of these genes resulted from (1) the sampling of several
families sharing similar clinical and histopathological
findings through international collaborations, (2) confir-
mation of the impact of pathogenic variants on the RNA,
protein and protein function through in vitro experi-
ments, and (3) the demonstration that the impairment
of the in vivo gene function correlated with the patient
phenotypes.

These novel genes revealed novel pathomechanism
underlying different CM. Genes like *CACNA1S* or
*ACTN2*, for which the physiological functions are well
known, fit with a pathomechanism implicating impaired
excitation-contraction coupling or sarcomere disor-
ganization, respectively [27, 65-67]. We also uncovered
novel genes with novel functions in skeletal muscle as the
oxidoreductase PYROXD1. How and to what extend the
impaired enzymatic activity in our patients contributes to the development of the muscle phenotype is currently not understood, but our functional experiments suggested that the myofibers are more sensitive to oxidative stress [23]. In addition, a novel pathway was implicated in myopathies with the discovery of pathogenic variants in *STIM1*, *ORAII*, and *CASQ1* that are key players of $Ca^{2+}$ homeostasis through the SOCE mechanism [68, 69]. The strict control of $Ca^{2+}$ entry, storage, and release takes a central role in skeletal muscle physiology and contraction [70].

Noteworthy, the first example of digenism in myopathies was provided with the identification of the co-segregation of deleterious variants in *SRPK3* and *TTN* in patients with progressive childhood-onset skeletal muscle myopathy with cores and centralized nuclei [29]. Such a scenario with pathogenic variants in two independent genes needs to be considered if NGS does not disclose a single implicated gene.

Novel myopathy genes represent novel therapeutic targets that may be more amenable to therapies compared to previously known genes. In case of pathogenic loss-of-function variants, gene replacement strategies could be envisaged even in the absence of a detailed knowledge on the pathomechanism. Inversely, the targeted down-regulation of genes, mRNA, or proteins can be applied for pathogenic variants involving a gain-of-function. The implication of novel pathological functions triggers the development of novel therapeutic approaches, while the implication of several proteins in the same pathway sustains the validation of common therapies. For example, repurposing of therapies targeting excitation-contraction coupling and $Ca^{2+}$ homeostasis can be envisaged for *CACNA1S* or the SOCE pathway [70], or anti-oxidant modulators for *PYROXD1*-related myopathy. Also, novel drugs and approaches can be developed through pre-clinical studies, as *ORAII* RNA silencing for *STIM1*-related TAM [71].

## Conclusions

Our results confirm the efficiency of exome sequencing to provide genetic diagnoses in patients with undiagnosed congenital myopathies, especially when combined with clinical and histopathological data, and access to additional family members. Most of the resolved cases were linked to a known myopathy gene with atypical signs, highlighting that most genetic diagnoses could not be found with targeted gene sequencing restricted to clinical-histological assessments. In addition, we identified 14 new myopathy genes, revealing novel pathomechanisms and potential therapeutic targets. Finally, we observed an important phenotypic and genotypic

heterogeneity. Overall, these results facilitated and accelerated the genetic diagnosis of unresolved families, improved health care for several patients, and opened novel perspectives for either repurposing of existing molecules or the development of novel treatments.

## Abbreviations

CNM Centronuclear myopathy

CM Congenital myopathy

NM Nemaline myopathy

SOCE Store-operated $Ca^{2+}$ entry

TAM Tubular aggregate myopathy

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s13073-024-01353-0.

Additional file 1: Table S1. Clinical, histological and pathogenic variants for probands from 123 families with signs of congenital myopathy. Variants were classified according to the ACMG standards and guidelines [18]. Genes are listed in alphabetical order. -: data not available, F: female, M: male, LL: lower limbs, UL: upper limbs, N: normal, CFTD: congenital fiber type disproportion, CNM: centronuclear myopathy, TAM: tubular aggregate myopathy, SD: splice donor. Homozygous pathogenic variants are shown in bold. ¹The probands of families 8 and 49 are third-degree relatives.

Additional file 2: Figs. S1–S3. Fig S1: Diagnostic yield of exome sequencing as a function of the number of included family members. Diagnosis: patients with pathogenic variant(s). Candidate: patients with suspected pathogenic variant(s) in one candidate gene. Fig S2: Diagnostic yield of exome sequencing as a function of the number of included affected family members. Diagnosis: patients with causal pathogenic variant(s). Candidate: patients with suspected pathogenic variant(s) in one candidate gene. Fig S3: Diagnostic yield of exome sequencing with or without muscle biopsy. Diagnosis: patients with causal pathogenic variant(s). Candidate: patients with suspected pathogenic variant(s) in one candidate gene.

## Acknowledgements

We thank patients and families for participation, clinicians, and the FILNEMUS network for providing DNA samples and phenotypic information, Vincent Laugel for scientific collaboration, Osorio Abath Neto and Kirsley Chennen for informatic analysis, and other partners of MYOCAPTURE consortium (Gisèle Bone, Nicolas Levy, Isabelle Richard, France Leturcq).

## Authors' contributions

JL conceived the project. NBR, VB, JB, and JL supervised the project. NBR and SS prepared and provided the materials. AB and JFD performed sequencing. YDF, MV, RS, VB, JB, and JL performed the analysis. JB and JL acquired fundings. YDF and JL wrote the manuscript with inputs from other co-authors. All authors read and approved the final manuscript.

## Funding

This work of the Interdisciplinary Thematic Institute IMCBio, as part of the ITI 2021–2028 program of the University of Strasbourg, CNRS and Inserm, was supported by IdEx Unistra (ANR-10-IDEX-0002), and by SFRI-STRAT'US project (ANR 20-SFRI-0012) and EUR IMCBio (ANR-17-EURE-0023) under the framework of the French Investments for the Future Program, the France Génomique National infrastructure funded as part of the Investissements d'Avenir program (ANR-10-INBS-09) and by GIS Maladies Rares within the frame of the "MYOCAPTURE" sequencing project, Agence Nationale pour la Recherche (CM-WES: ANR-11-BSV1-026-01), the Fondation pour la Recherche Médicale (DBI20131228569), the Horizon Europe program from European Union (COMPASS-NMD: 101080874), and the Association Française contre les Myopathies AFM-Téléthon (17088, 20323, 21267, 22734, 23933). YF is a fellow from Inserm poste d'accueil.
## Availability of data and materials

Due to privacy, ethical, and legal issues in France and the European Union and according to the content of our IRB approved project, exome data supporting the findings of this study are not publicly available. The described pathogenic variants and the associated phenotypes supporting the conclusions of this article are available in the ClinVar repository under the submissions SUB14382184 and SUB14484902 (https://www.ncbi.nlm.nih.gov/clinvar). The datasets combining pathogenic variants and phenotypes of the patients are included in the supporting files in a machine-readable format (.xlxs).

## Declarations

### Ethics approval and consent to participate

The present study was approved by the IRB Comité de Protection des Personnes Est IV (DC-2012-1693), and all participants or legal guardians provided written informed consents to the different clinical centers involved in the project. Our research conforms to the principles of the Helsinki Declaration.

### Consent for publication

Not applicable.

### Competing interests

Safaa Saker is employee of a for-profit organization (Genethon). The other authors declare that they have no competing interests.

### Author details

¹IGBMC, Inserm U1258, Cnrs UMR7104, Université de Strasbourg, 1 Rue Laurent Fries, Illkirch 67404, France. ²Department of Pediatric Neurology, CHU Strasbourg, Strasbourg, France. ³Centre de Référence Neuromusculaire Nord-Est-Ile de France, Strasbourg, France. ⁴Myology Institute, Neuromuscular Morphology Unit, Sorbonne Université, INSERM, GHU Pitié-Salpêtrière, Paris, France. ⁵Genethon, DNA and Cell Bank, Evry 91000, France. ⁶Centre National de Recherche en Génomique Humaine (CNRGH), Université Paris-Saclay, CEA, Evry 91057, France. ⁷Laboratoire de Diagnostic Génétique CHRU de Strasbourg, Strasbourg 67091, France.

Received: 30 October 2023 Accepted: 30 May 2024

Published online: 09 July 2024

## References

1. Amburgey K, McNamara N, Bennett LR, McCormick ME, Acsadi G, Dowling JJ. Prevalence of congenital myopathies in a representative pediatric united states population. Ann Neurol. 2011;70:662–5.

2. North KN, Wang CH, Clarke N, Jungbluth H, Vainzof M, Dowling JJ, et al. Approach to the diagnosis of congenital myopathies. Neuromuscul Disord NMD. 2014;24:97–116.

3. Wang CH, Dowling JJ, North K, Schroth MK, Sejersen T, Shapiro F, et al. Consensus statement on standard of care for congenital myopathies. J Child Neurol. 2012;27:363–82.

4. Nicolau S, Liewluck T, Tracy JA, Laughlin RS, Milone M. Congenital myopathies in the adult neuromuscular clinic. Neurol Genet. 2019;5:e341.

5. Papadimas GK, Kirou S, Kararizou E, Papadopoulos C. Update on congenital myopathies in adulthood. Int J Mol Sci. 2020;21:3694.

6. Pinto MJ, Passos BA, Grangeia A, Guimarães J, Braz L. Congenital myopathies in adults: a diagnosis not to overlook. Acta Neurol Scand. 2022;146:152–9.

7. Kaplan J-C, Hamroun D. The 2013 version of the gene table of monogenic neuromuscular disorders (nuclear genome). Neuromuscul Disord NMD. 2012;22:1108–35.

8. Benarroch L, Bonne G, Rivier F, Hamroun D. The 2023 version of the gene table of neuromuscular disorders (nuclear genome). Neuromuscul Disord NMD. 2023;33:76–117.

9. Vasli N, Laporte J. Impacts of massively parallel sequencing for genetic diagnosis of neuromuscular disorders. Acta Neuropathol (Berl). 2013;125:173–85.

10. Vasli N, Böhm J, Le Gras S, Muller J, Pizot C, Jost B, et al. Next generation sequencing for molecular diagnosis of neuromuscular diseases. Acta Neuropathol (Berl). 2012;124:273–83.

11. Schartner V, Romero NB, Donkervoort S, Treves S, Munot P, Pierson TM, et al. Dihydropyridine receptor (DHPR, CACNA1S) congenital myopathy. Acta Neuropathol (Berl). 2017;133:517–33.

12. Geoffroy V, Pizot C, Redin C, Piton A, Vasli N, Stoetzel C, et al. VaRank: a simple and powerful tool for ranking genetic variants. PeerJ. 2015;3:e796.

13. Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31:3812–4.

14. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.

15. Yeo G, Burge CB. Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals. J Comput Biol J Comput Mol Cell Biol. 2004;11:377–94.

16. Reese MG, Eeckman FH, Kulp D, Haussler D. Improved splice site detection in Genie. J Comput Biol J Comput Mol Cell Biol. 1997;4:311–23.

17. Shapiro MB, Senapathy P. RNA splice junctions of different classes of eukaryotes: sequence statistics and functional implications in gene expression. Nucleic Acids Res. 1987;15:7155–74.

18. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetic and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24.

19. Pedersen BS, Bhetariya PJ, Brown J, Kravitz SN, Marth G, Jensen RL, et al. Somalier: rapid relatedness estimation for cancer and germline studies using efficient genome sketches. Genome Med. 2020;12:62.

20. Böhm J, Chevessier F, Maues De Paula A, Koch C, Attarian S, Feger C, et al. Constitutive activation of the calcium sensor STIM1 causes tubular-aggregate myopathy. Am J Hum Genet. 2013;92:271–8.

21. Böhm J, Bulla M, Urquhart JE, Malfatti E, Williams SG, O'Sullivan J, et al. ORAI1 Mutations with distinct channel gating defects in tubular aggregate myopathy. Hum Mutat. 2017;38:426–38.

22. Böhm J, Lornage X, Chevessier F, Birck C, Zanotti S, Cudia P, et al. CASQ1 mutations impair calsequestrin polymerization and cause tubular aggregate myopathy. Acta Neuropathol (Berl). 2018;135:149–51.

23. O'Grady GL, Best HA, Sztal TE, Schartner V, Sanjuan-Vazquez M, Donkervoort S, et al. Variants in the oxidoreductase PYROXD1 cause early-onset myopathy with internalized nuclei and myofibrillar disorganization. Am J Hum Genet. 2016;99:1086–105.

24. Vasli N, Harris E, Karamchandani J, Bareke E, Majewski J, Romero NB, et al. Recessive mutations in the kinase ZAK cause a congenital myopathy with fibre type disproportion. Brain J Neurol. 2017;140:37–48.

25. Donkervoort S, Kutzner CE, Hu Y, Lornage X, Rendu J, Stojkovic T, et al. Pathogenic variants in the myosin chaperone UNC-45B cause progressive myopathy with eccentric cores. Am J Hum Genet. 2020;107:1078–95.

26. Malfatti E, Böhm J, Lacène E, Beuvin M, Romero NB, Laporte J. A premature stop codon in MYO18B is associated with severe nemaline myopathy with cardiomyopathy. J Neuromuscul Dis. 2015;2:219–27.

27. Lornage X, Romero NB, Grosgogeat CA, Malfatti E, Donkervoort S, Marchetti MM, et al. ACTN2 mutations cause "Multiple structured Core Disease" (MsCD). Acta Neuropathol (Berl). 2019;137:501–19.

28. Foley AR, Zou Y, Dunford JE, Rooney J, Chandra G, Xiong H, et al. GGPS1 mutations cause muscular dystrophy/hearing loss/ovarian insufficiency syndrome. Ann Neurol. 2020;88:332–47.

29. Töpf A, Cox DT, Zaharieva I, Di Leo V, Sarparanta J, Harald Jonson P, et al. Digenic inheritance involving a muscle specific protein kinase and the giant titin protein causes a skeletal muscle myopathy. Nat Genet. 2024;56:395–407.

30. Lornage X, Quijano-Roy S, Amthor H, Carlier RY, Monnier N, Deleuze J-F, et al. Asymmetric muscle weakness due to ACTA1 mosaic mutations. Neurology. 2020;95:e3406–11.

31. Biancalana V, Rendu J, Chaussenot A, Mecili H, Bieth E, Fradin M, et al. A recurrent RYR1 mutation associated with early-onset hypotonia and benign disease course. Acta Neuropathol Commun. 2021;9:155.

32. Biancalana V, Romero NB, Thuestad IJ, Ignatius J, Kataja J, Gardberg M, et al. Some DNM2 mutations cause extremely severe congenital myopathy and phenocopy myotubular myopathy. Acta Neuropathol Commun. 2018;6:93.

33. Bevilacqua JA, Monnier N, Bitoun M, Eymard B, Ferreiro A, Monges S, et al. Recessive RYR1 mutations cause unusual congenital myopathy
with prominent nuclear internalization and large areas of myofibrillar disorganization. Neuropathol Appl Neurobiol. 2011;37:271–84.

34. Wilmshurst JM, Lillis S, Zhou H, Pillay K, Henderson H, Kress W, et al. RYR1 mutations are a common cause of congenital myopathies with central nuclei. Ann Neurol. 2010;68:717–26.

35. Garibaldi M, Rendu J, Brocard J, Lacene E, Fauré J, Brochier G, et al. "Dusty core disease" (DuCD): expanding morphological spectrum of RYR1 recessive myopathies. Acta Neuropathol Commun. 2019;7:3.

36. Ceyhan-Birsoy O, Agrawal PB, Hidalgo C, Schmitz-Abe K, DeChene ET, Swanson LC, et al. Recessive truncating titin gene, TTN, mutations presenting as centronuclear myopathy. Neurology. 2013;81:1205–14.

37. Berardo A, Lornage X, Johari M, Evangelista T, Cejas C, Barroso F, et al. HNRNPDL-related muscular dystrophy: expanding the clinical, morphological and MRI phenotypes. J Neurol. 2019;266:2524–34.

38. Lornage X, Mallaret M, Silva-Rojas R, Biancalana V, Giovannini D, Dieterich K, et al. Selective loss of a LAP1 isoform causes a muscle-specific nuclear envelopathy. Neurogenetics. 2021;22:33–41.

39. Villar-Quiles RN, Catervi F, Cabet E, Juntas-Morales R, Genetti CA, Gidaro T, et al. ASC-1 is a cell cycle regulator associated with severe and mild forms of myopathy. Ann Neurol. 2020;87:217–32.

40. Davignon L, Chauveau C, Julien C, Dill C, Duband-Goulet I, Cabet E, et al. The transcription coactivator ASC-1 is a regulator of skeletal myogenesis, and its deficiency causes a novel form of congenital muscle disease. Hum Mol Genet. 2016;25:1559–73.

41. Böhm J, Malfatti E, Oates E, Jones K, Brochier G, Boland A, et al. Novel ASCC1 mutations causing prenatal-onset muscle weakness with arthro-gyposis and congenital bone fractures. J Med Genet. 2019;56:617–21.

42. Echaniz-Laguna A, Lornage X, Lannes B, Schneider R, Bierry G, Dondaine N, et al. HSPB8 haploinsufficiency causes dominant adult-onset axial and distal myopathy. Acta Neuropathol (Berl). 2017;134:163–5.

43. Lornage X, Malfatti E, Chéraud C, Schneider R, Biancalana V, Cuisset J-M, et al. Recessive MYPN mutations cause cap myopathy with occasional nemaline rods. Ann Neurol. 2017;81:467–73.

44. Scott HS, Litjens T, Nelson PV, Thompson PR, Brooks DA, Hopwood JJ, et al. Identification of mutations in the alpha-L-iduronidase gene (IDUA) that cause Hurler and Scheie syndromes. Am J Hum Genet. 1993;53:973–86.

45. Pinto MM, Monges S, Malfatti E, Lubieniecki F, Lornage X, Alias L, et al. Sarcomeric disorganization and nemaline bodies in muscle biopsies of patients with EXOSC3-related type 1 pontocerebellar hypoplasia. Muscle Nerve. 2019;59:137–41.

46. Morin G, Biancalana V, Echaniz-Laguna A, Noury J-B, Lornage X, Moggio M, et al. Tubular aggregate myopathy and Stormorken syndrome: mutation spectrum and genotype/phenotype correlation. Hum Mutat. 2020;41:17–37.

47. Nesin V, Wiley G, Kousi M, Ong E-C, Lehmann T, Nicholl DJ, et al. Activating mutations in STIM1 and ORAI1 cause overlapping syndromes of tubular myopathy and congenital miosis. Proc Natl Acad Sci U S A. 2014;111:4197–202.

48. Garibaldi M, Fattori F, Riva B, Labasse C, Brochier G, Ottaviani P, et al. A novel gain-of-function mutation in ORAI1 causes late-onset tubular aggregate myopathy and congenital miosis. Clin Genet. 2017;91:780–6.

49. Endo Y, Noguchi S, Hara Y, Hayashi YK, Motomura K, Miyatake S, et al. Dominant mutations in ORAI1 cause tubular aggregate myopathy with hypocalcemia via constitutive activation of store-operated Ca²⁺ channels. Hum Mol Genet. 2015;24:637–48.

50. Barone V, Del Re V, Gamberucci A, Polverino V, Galli L, Rossi D, et al. Identification and characterization of three novel mutations in the CASQ1 gene in four patients with tubular aggregate myopathy. Hum Mutat. 2017;38:1761–73.

51. Lornage X, Schartner V, Balbueno I, Biancalana V, Willis T, Echaniz-Laguna A, et al. Clinical, histological, and genetic characterization of PYROXD1-related myopathy. Acta Neuropathol Commun. 2019;7:138.

52. Alazami AM, Kentab AY, Faqei E, Mohamed JY, Alkhalidi H, Hijazi H, et al. A novel syndrome of Klippel-Feil anomaly, myopathy, and characteristic facies is linked to a null mutation in MYO18B. J Med Genet. 2015;52:400–4.

53. Misceo D, Holmgren A, Louch WE, Holme PA, Mizobuchi M, Morales RJ, et al. A dominant STIM1 mutation causes Stormorken syndrome. Hum Mutat. 2014;35:556–64.

54. Töpf A, Johnson K, Bates A, Phillips L, Chao KR, England EM, et al. Sequential targeted exome sequencing of 1001 patients affected by unexplained limb-girdle weakness. Genet Med Off J Am Coll Med Genet. 2020;22:1478–88.

55. Ghaoui R, Cooper ST, Lek M, Jones K, Corbett A, Reddel SW, et al. Use of whole-exome sequencing for diagnosis of limb-girdle muscular dystrophy: outcomes and lessons learned. JAMA Neurol. 2015;72:1424–32.

56. Savarese M, Di Fruscio G, Torella A, Fiorillo C, Magri F, Fanin M, et al. The genetic basis of undiagnosed muscular dystrophies and myopathies: results from 504 patients. Neurology. 2016;87:71–6.

57. Laquerriere A, Jaber D, Abiusi E, Maluenda J, Mejlachowicz D, Vivanti A, et al. Phenotypic spectrum and genomics of undiagnosed arthrogryposis multiplex congenita. J Med Genet. 2022;59:559–67.

58. Cummings BB, Marshall JL, Tukiainen T, Lek M, Donkervoort S, Foley AR, et al. Improving genetic diagnosis in Mendelian disease with transcriptome sequencing. Sci Transl Med. 2017;9:eaa15209.

59. Wattjes MP, Kley RA, Fischer D. Neuromuscular imaging in inherited muscle diseases. Eur Radiol. 2010;20:2447–60.

60. Mercuri E, Pichiecchio A, Allsop J, Messina S, Pane M, Muntoni F. Muscle MRI in inherited neuromuscular disorders: past, present, and future. J Magn Reson Imaging JMRI. 2007;25:433–40.

61. Perry L, Stimpson G, Singh L, Morrow JM, Shah S, Baranello G, et al. Muscle magnetic resonance imaging involvement patterns in nemaline myopathies. Ann Clin Transl Neurol. 2023;10:1219–29.

62. Kraeva N, Heytens L, Jungbluth H, Treves S, Voermans N, Kamsteeg E, et al. Compound RYR1 heterozygosity resulting in a complex phenotype of malignant hyperthermia susceptibility and a core myopathy. Neuromuscul Disord NMD. 2015;25:567–76.

63. Romero NB. Centronuclear myopathies: a widening concept. Neuromuscul Disord NMD. 2010;20:223–8.

64. Malfatti E, Schaeffer U, Chapon F, Yang Y, Eymard B, Xu R, et al. Combined cap disease and nemaline myopathy in the same patient caused by an autosomal dominant mutation in the TPM3 gene. Neuromuscul Disord NMD. 2013;23:992–7.

65. Rossi D, Catallo MR, Pierantozzi E, Sorrentino V. Mutations in proteins involved in E-C coupling and SOCE and congenital myopathies. J Gen Physiol. 2022;154:e202213115.

66. Schartner V, Laporte J, Böhm J. Abnormal excitation-contraction coupling and calcium homeostasis in myopathies and cardiomyopathies. J Neuromuscul Dis. 2019;6:289–305.

67. Jungbluth H, Treves S, Zorzato F, Sarkozy A, Ochala J, Sewry C, et al. Congenital myopathies: disorders of excitation-contraction coupling and muscle contraction. Nat Rev Neurol. 2018;14:151–67.

68. Protasi F, Girolami B, Roccabianca S, Rossi D. Store-operated calcium entry: from physiology to tubular aggregate myopathy. Curr Opin Pharmacol. 2023;68:102347.

69. Silva-Rojas R, Laporte J, Böhm J. STIM1/ORAI1 Loss-of-function and gain-of-function mutations inversely impact on SOCE and calcium homeostasis and cause multi-systemic mirror diseases. Front Physiol. 2020;11:604941.

70. Riva B, Pessolano E, Quaglia E, Cordero-Sanchez C, Bhela IP, Topf A, et al. STIM1 and ORAI1 mutations leading to tubular aggregate myopathies are sensitive to the Store-operated Ca²⁺-entry modulators CIC-37 and CIC-39. Cell Calcium. 2022;105:102605.

71. Silva-Rojas R, Pérez-Guàrdia L, Lafabrie E, Moulaert D, Laporte J, Böhm J. Silencing of the Ca²⁺ channel ORAI1 improves the multi-systemic phenotype of tubular aggregate myopathy (TAM) and Stormorken Syndrome (STRMK) in mice. Int J Mol Sci. 2022;23:6968.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.