‘Dusty core disease’ (DuCD): expanding morphological spectrum of RYR1 recessive myopathies

Several morphological phenotypes have been associated to RYR1-recessive myopathies. We recharacterized the RYR1-recessive morphological spectrum by a large monocentric study performed on 54 muscle biopsies from a large cohort of 48 genetically confirmed patients, using histoenzymology, immunohistochemistry, and ultrastructural studies. We also analysed the level of RyR1 expression in patients’ muscle biopsies. We defined “dusty cores” the irregular areas of myofibrillar disorganisation characterised by a reddish-purple granular material deposition with uneven oxidative stain and devoid of ATPase activity, which represent the characteristic lesion in muscle biopsy in 54% of patients. We named Dusty Core Disease (DuCD) the corresponding entity of congenital myopathy. Dusty cores had peculiar histological and ultrastructural characteristics compared to the other core diseases. DuCD muscle biopsies also showed nuclear centralization and type1 fibre predominance. Dusty cores were not observed in other core myopathies and centronuclear myopathies. The other morphological groups in our cohort of patients were: Central Core (CCD: 21%), Core-Rod (C&R:15%) and Type1 predominance “plus” (T1P+:10%). DuCD group was associated to an earlier disease onset, a more severe clinical phenotype and a lowest level of RyR1 expression in muscle, compared to the other groups. Variants located in the bridge solenoid and the pore domains were more frequent in DuCD patients. In conclusion, DuCD is the most frequent histopathological presentation of RYR1-recessive myopathies. Dusty cores represent the unifying morphological lesion among the DuCD pathology spectrum and are the morphological hallmark for the recessive form of disease. Electronic supplementary material The online version of this article (10.1186/s40478-018-0655-5) contains supplementary material, which is available to authorized users.


Introduction
The RYR1 gene encodes the ryanodine receptor channel 1 (RyR1), a sarcoplasmic reticulum (SR) calcium channel involved in excitation-contraction coupling through interaction with the dihydropyridine receptor (DHPR) in the T-tubule. RYR1-mutations are the most common cause of congenital myopathies with cores [14]. Variable clinical and histological phenotypes have been described in both dominant and recessive forms [36]. Dominant mutations in N-terminal and central regions (MHS/ CCD1 and MHS/CCD2 domains respectively) are classically associated to Malignant Hyperthermia Susceptibility (MHS) trait, whereas dominant mutations in C-terminal region (MHS/CCD3) are more likely associated to Central Core Disease (CCD) [40]. Moreover, heterozygous variants have been previously reported in core-rod myopathy [12,27,33] and more recently in Exertional Rhabdomyolysis (ERM) [9] and late-onset axial myopathy [15,23].
By contrast, recessive forms of disease manifest a wider spectrum of clinical and histological presentations. Even if RYR1-recessive patients could present a clinical picture similar to the dominant one, including an early onset, non-progressive, proximal muscle weakness, other peculiar clinical phenotypes of recessive forms are characterized by a more diffuse muscle weakness, ocular involvement with ptosis and/or ophthalmoplegia, and a severe bulbar and respiratory muscle weakness [19,20,35]. Nevertheless, severe clinical presentations and foetal akinesia syndrome have been reported in both dominant and recessive forms [5,12,31].
Recessive cases also present an increasing complexity, given the large amount of novel variants detected by the massive sequencing technologies, which makes the genotype-phenotype correlation more challenging. In general, hypomorphic variants (non-sense, frameshift, splice site variants), reducing or abolishing RyR1 production, seem to be more frequent in recessive forms and are associated to a more severe clinical presentation, compared to non-hypomorphic (missense, in-frame ins/ del) variants [26,40]. Nevertheless, up to now, no association between clinical severity and histopathological findings has been found [2].
We report an extensive monocentric analysis of 54 muscle biopsies from a large cohort of 48 RYR1-recessive patients to have a homogenous interpretation of morphological findings with the aim to find a closer correlation between morphology, clinical phenotype and genetic background. Here, we describe and define the "dusty core fibres" as the characteristic and unifying morphological feature present in most of RYR1-recesssive biopsies and we correlate the RyR1 expression level in patients' muscle biopsies with clinical and morphological features.

Patients' sample selection
All the muscle biopsies were analysed at the Neuromuscular Morphology Unit of Myology Institute, in Paris. More than 11000 muscle biopsies collected between 1977 and 2015 were screened. Two hundred and thirty belonged to patients with RYR1-related myopathy, of whom 154 had dominant inheritance and 76 were confirmed or suspected to be recessive. In 13 cases the second variant was not found or was not surely pathogenetic (likely pathogenetic, VUS or likely benign). Among the 63 confirmed RYR1-recessive patients, 15 cases have been excluded because of insufficient data (uncomplete clinical and/or morphological data and deteriorated muscle sample for re-analysis). Finally, our study cohort consisted on 48 confirmed RYR1-recessive patients (20 male and 28 female) from 45 unrelated families. Clinical data were obtained from a full revision of all available medical records up to the last clinical examination in all enrolled patients. Disease onset (considered as the first reported clinical sign referred to disease, including perinatal problems or delayed motor milestones), weakness distribution (proximal/ distal limb muscles, axial, facial, ocular/extraocular and bulbar muscles), contractures, spinal deformities and dysmorphisms, respiratory and cardiac involvement were considered. Clinical evaluations were performed by different clinicians over 40 years, thus we retrospectively classified patients in three main groups of clinical severity: mild (late-adult onset, walkers, mild muscle weakness, minimal or absent ocular, facial, bulbar or respiratory involvement), moderate (early onset, non-progressive, proximal or diffuse muscle weakness, associate with mild dysmorphism, spine deformities or contractures) and severe (muscle hypotonia at birth, feeding difficulties, severe respiratory involvement requiring ventilation, contractures and/or spinal deformities, diffuse muscle weakness with facial and ocular involvement).

Genetic analysis
Total RNA was extracted from each skeletal muscle sample lysed in Trizol reagent (Invitrogen, Life Technologies SAS). Complementary DNA was synthesized from 500 to 750 ng of total RNA using 0.5 μl of Transcriptor (Roche) and 0.3 lg of oligo-dT as described [26]. Seven overlapping PCR amplification spanning the entire RYR1 sequence were performed. Each fragment was sequenced as previously described [26,28]. Each variation was confirmed on DNA sample and on both paternal and maternal DNA sample to establish the transmission. Each variant was analysed by Variant effect Predictors to obtain the different prediction score such as CADD, SIFT, Polyphen and gnomAD exome and gno-mAD Genome database frequency. To better assess the functional effect of each missense variation, 3D analysis was performed on Yasara sofware [21]. Due to the large size of the RYR1 gene, we choose not to use the total RyR1 protein structure already described (5gl1/5taz) [3,8] in the first-round analysis via FoldX prediction. We choose to split the structure in 5 parts spanning the whole human RYR1 structures (amino acid 1 to 627, 628 to 1656, 1657 to 2144, 2145 to 3613 and 3614 to 5038). Then, the sequences were submitted in I-TASSER server [39] to obtain "friendly" usable RyR1 structure. Each structure prediction was matched with the RyR1 global structure (5gl1 and 5taz) [3,8]. Delta G variations were calculated to estimate protein stability. For delta G variation <0.5 kcal/mol, meaning no destabilization, study of the whole structure was realized (5gl1/5taz) [3,8]. For ACMG classification, Intervar was used with recessive transmission correction [22].

Histological study
Histoenzymological analysis was conducted on 54 muscle biopsies (4 patients had two muscle biopsies and 1 patient three muscle biopsies available in the Myology Institute Lab). Age at muscle biopsy ranged from 1 day of life (30 weeks of adjusted gestational age) to 76 years (median 16 years, . Open muscle biopsies were obtained from deltoid or quadriceps muscles in most of patients. Histological and histochemical slides were systematically re-analysed by two authors (MG and NBR) with experience in skeletal muscle morphology, blinded to clinical and molecular data. For the oldest, deteriorated or not interpretable slides, new slides were obtained from the best muscle specimen available in the lab. Conventional histological and histochemical techniques, 8-10 μm thick cryostat sections were stained with haematoxylin and eosin (HE), modified Gomori trichrome (GT), Periodic acid Schiff technique (PAS), Oil red O, reduced nicotinamide adenine dinucleotide dehydrogenase-tetrazolium reductase (NADH-TR), succinic dehydrogenase (SDH), cytochrome c oxidase (COX), and adenosine triphosphatase (ATPase) preincubated at pH 9.4, 4.63, 4.35. Digital photographs of biopsies were obtained with a Zeiss AxioCam HRc linked to a Zeiss Axioplan Bright Field Microscope and processed with the Axio Vision 4.4 software (Zeiss, Germany). Fibre type pattern was determined in ATPases reactions, and by calculating the percentage of type 1 and type 2 fibres. We considered type 1 fibres predominance when there were more than of 60% type 1 fibre in deltoid muscles, and more than 40% in quadriceps muscle. CFTD was considered when all the type1 fibres were consistently (at least 35-40%) smaller than type2 fibres in absence of other pathological findings. Centronuclear pattern was considered only when myonuclei were centrally placed in almost 50% of muscle fibres showing nuclear internalization. Rods were considered in the presence of numerous, multiple, small nemaline bodies both in cytoplasmic or subsarcolemmal areas. Central cores were considered when single or multiple sharply demarked ovoidal areas devoid of oxidative stains were observed in transversal sections of type1 muscle fibres, centrally or peripherally placed. Multiminicore were considered in the presence of boundless, small areas of decreased enzymatic activity at oxidative stains [30]. "Dusty cores" were defined as irregular areas of reddish-purple granular material deposition at GT stain corresponding to decreased or/and increased enzymatic activity at oxidative stains and devoid of ATPase activity. Patients with available ultrastructural study were finally classified considering both histological and ultrastructural features. In the five patients with two or three muscle biopsies, final morphological classification was reached considering both muscle biopsies and most relevant findings.

Immunohistochemical (IHC) study
IHC was performed in new sections from available frozen muscle samples of enrolled patients, rejecting oldest and/ or deteriorated specimens. Finally, IHC analysis was available for 23 muscle biopsies. Antibodies against Desmin (Anti-Human Desmin, Clone D33, Dako Laboratories, Denmark A/S ), Myotilin (NCL-Myotilin, Novocastra Laboratories, Newcastle Upon Tyne, United Kingdom) and αB-crystallin (CRYAB, GeneTex International Corporation, Irvine, USA) were visualized using immunoperoxidase techniques. Immunofluorescence study was performed for RyR (anti-Ryanodine Receptor, clone 34C, Sigma Laboratories, Saint Louis, Missouri, USA), DHPR (anti-CACNA1S, ab2862, abcam Laboratories, Cambridge, UK) and alpha-actinin (anti-α actinin sarcomeric, clone EA-53, Sigma Aldrich Laboratories, Saint Louis, USA) antibodies on 10-μm-thick cryosections over night at 4°C. Subsequently, sections were incubated with appropriate conjugated secondary antibodies (Alexa Fluor-488 goat anti-rabbit antibody and Jackson IR goat anti-mouse antibody) for one hour. A set of control slides was prepared with omission of the primary antibodies.

Electron microscopy (EM) study
Ultrastructural study was newly performed in all available muscle biopsies. EM images were obtained for 39 muscle biopsies. Small muscle specimens were fixed with glutaraldehyde (2.5%, pH 7.4), post fixed with osmium tetroxide (2%), dehydrated and embedded in resin. Longitudinally oriented ultra-thin sections were obtained at different level of deepness from 1 to 3 small blocks and stained with uranyl acetate and lead citrate. Ultra-thin sections of transversally oriented blocks were obtained only for the most significant findings. The grids were observed using a Philips CM120 electron microscope (80 kV; Philips Electronics NV, Eindhoven, The Netherlands) and were photo documented using a Morada camera (Soft Imaging System).

Western blot (WB) analysis
The amount of RyR1 in muscle samples was determined by quantitative Western Blot (WB) analysis using antibodies directed against RyR1 [24] and normalized to the amount of myosin heavy chain as described previously [26]. Briefly, the muscle sample (20-40mg) was homogenized in 200 mM sucrose, 20mM HEPES (pH 7.4), 0.4mMCaCl2, 200mM phenylmethylsulfonyl fluoride, 1 mM diisopropyl Fluorophosphate using a Minilys homogenizer (Bertin, France). After electrophoretic separation on a 4-20% gradient acrylamide gel (Biorad, France) and electrotransfer to Immobilon P (Biorad, France) during 4h at 0.8A to ensure a complete transfer of the loaded proteins, the membrane was incubated with anti-RyR1 antibodies and then HRP-labelled secondary antibodies (Jackson ImmunoResearch Laboratories). Signal quantification was performed using a ChemiDoc Touch apparatus (Biorad, France) and the Image Lab software (Biorad). The total amount of RyR1 in each experiment was corrected from the amount of myosin and normalized to the amount of RyR1 present in the control referred as 100% as described previously [6]. Controls (muscle biopsy from individuals non-affected by neuromuscular disease) of different age have been used: 10 days, 3 years, 23 years and 46 years, and the amount of RyR1 in the muscle biopsy of the patient was compared to the age-related control. For each muscle biopsy, at least 3 western blots have been performed, and the value for each patient is presented as mean ± SEM of the different western blots.

Statistical analysis
Data have been expressed as range, median and interquartile range (IQR) for continuous variables (age and % of RyR expression) and as absolute values and frequencies for categorical variables (gender, disease severity, age at onset, morphology and occurrence of ocular involvement). Age at disease onset has been categorized as <1 year or ≥1 year, based on the median value in the sample. Chi-Square test has been used to compare categorical variables, while Kruskal-Wallis or Mann-Whitney tests have been used to evaluate the distribution of continuous variables with respect to demographic, clinical and morphological characteristics. P values lower than 0.05 were considered to be statistically significant. Statistical

Results
Overall morphological, clinical and genetic data are summarized in Table 1.

Clinical features
Age at onset ranged from antenatal period to 66 years (median 1 year, IQR 0-7). The first symptom/sign occurred within the first year of life in 23 patients (48%), of whom five had also antenatal manifestation. Disease onset occurred after the first year of life in 25 patients (52%) of whom 4 had a relative late onset (>20 years). Motor milestones were delayed in 25 patients (52%) and 4 out of 8 severe patients never acquired ambulation. All the patients with latest onset (>45 years) had not delayed motor milestones, as well as 9 (19%) with early-infantile onset. Ocular involvement was assessed in 38 of 48 patients and showed isolated eyelid ptosis in 4 patients, isolated ophthalmoplegia (sometimes isolated upper gaze movements) in 7 and both in 9. Overall, ocular involvement occurred in 20/38 patients and prevailed in subjects with moderate/severe disease (88%) compared to those with mild disease (24%) (p<0.0001). Overall clinical severity was mild in 25 (52%), moderate in 15 (31%) and severe in 8 patients (17%). All the clinical data are summarised in Table 2.

Morphological findings
We defined "dusty core fibres" the muscle fibres with irregular areas of reddish-purple granular material deposition at GT stain corresponding to large areas of uneven myofibrillar disorganisation at oxidative stains, characterised by blended decreased or/and increased enzymatic activity and devoid ATPase activity ("dusty cores") ( Fig. 1a-f ). These dusty core fibres are different from the classic central cores observed in CCD because the last cores are well delimited and have not fuchsine deposit inside the core.
Dusty cores were the most frequent histopathological lesion, detected in 32 muscle biopsies and confirmed by EM in 30 muscle biopsies (56,6%). We named "Dusty Core Disease" (DuCD) the corresponding group of congenital core myopathy with dusty cores.
The number of muscle fibres containing dusty cores ranged greatly from few fibres up to 30-40% of fibres in muscle biopsy.
Dusty cores were consistently associated to type1 fibre predominance or uniformity and prominent nuclear internalization and centralization. Some muscle biopsies presented some supplementary peculiar features: in 3 muscle biopsies (p4, p14, p17) dusty cores were detected only in few (2-5) muscle fibres only after an extensive and careful revision; whereas in 3 others (p14, p15, p21) none or minor nuclear internalization was observed; 3 muscle biopsies (p9, p18, p24) presented a high fibre size variability with hypertrophic and atrophic fibres, fibre slitting, increased connective tissue without necrosis and regeneration (pseudo-dystrophic appearance); five muscle biopsies (p1, p11, p16, p19, p20) showed a prominent nuclear internalization and centralization (CNM-like) (Fig. 1g,h); one case (p26) showed surprisingly both dusty and central cores. These findings provide the evidence that a certain variability exists among the DuCD morphological spectrum, in which dusty cores represent the unifying lesion.Regardless of these peculiarities, in 14/31 cases (45%) dusty cores presented a "targetoid" appearance in variable amount of muscle fibres (Fig. 1i,l). The IHC study for myofibrillar proteins (desmin, myotilin, αB-crystallin) did not show any specificity in type1 predominance group. By contrast central core and dusty core groups showed an inconsistently positive immunostaining of all three antibodies inside cores (Fig. 2a-c, g-i), with the only peculiarity for desmin which showed a peripheral halo of cores in some biopsies of CCD group. Likewise, immunofluorescence for RyR1 and DHPR showed an inconstant positivity for central core and dusty core groups, even if in DuCD group positive signal by RyR1 and DHPR did not match ( Fig. 2d-f, l-n). No specificity was observed for type1 predominance group and by a-actinin immunostaining for all groups.
The EM study showed some peculiarities in the DuCD group. In longitudinal sections, the areas of sarcomeric disorganization consisted of cellular debris with sarcomeric fragments of thin and thick filaments devoid of mitochondria. In addition, two abnormal elements were consistently observed: abundant electrodense longitudinally-smeared material and thickened short Z-line darker fragments, sometime superposed (Fig. 3a,b). These features were evident also in transversal sections (Fig. 3c), inside areas of disorganization which frequently showed sub-areas of electrodense material deposition, mixed close sub-areas of less osmophilic material. In longitudinal sections, areas of sarcomeric disorganization ranged from few to more than 50 sarcomeres (Fig. 3a). In all cases of dusty core group, smaller areas of sarcomeric disorganization, resembling minicores, were noticed (Fig. 3d). The larger areas of disorganization, possibly corresponding to dusty cores at optic microscopy, had an asymmetric or «star-like» shape, with a longer longitudinal axis and a shorter transversal axis. Sometimes, these areas of disorganization occupy the entire muscle fibre in width (Fig. 3e). Interestingly in almost half of cases, several triads duplication or multiplication were detected inside areas of disorganization (Fig. 3f). Nevertheless, in all suspected minicore lesions by EM, a deeper analysis, performed with numerous serial sections, revealed the presence of bigger areas of disorganization. Finally, Z-line abnormalities were detected in the majority of cases, including duplication, streaming and smearing across 2-4 sarcomeres.
Both central core and core-rod groups showed at EM classic unstructured and, less frequently, structured cores. Nemaline rods were consistently observed in all cases of the core-rod group, in clusters of subsarcolemmal areas or scattered in sarcoplasm, frequently in perinuclear areas.
Overall, the ultrastructural study confirmed the optic microscopy classification in 37 out of 39 muscle biopsies analysed. All cases of core-rod were confirmed by EM. Surprisingly in one case of CCD and one case of DuCD no abnormalities were found at the EM study in the specimens analysed.
Concerning patients with more than one muscle biopsy available, the average time between the two biopsies was 12,6 years (range 6-22 years). The first muscle biopsy showed the presence of dusty cores in 3 out of 5 patients, confirmed also in the second muscle biopsy, two of which also showed a "targetoid" appearance and one a pseudo-dystrophic pattern. Interestingly, the first muscle biopsy in 2 cases showed only nuclear centralization in small type1 fibres without any myofibrillar disorganization detected by optic microscopy. Both cases revealed dusty cores in the subsequent muscle biopsies performed some years later.

Molecular data
Every patient has at least one variation on each allele of RYR1 gene (Additional file 1). Each variant was identified based on its low frequency (less than 0.5% in database gnomAD). Thirteen variants were not previously reported. Twenty-six hypomorphic variants were found (non-sense, out of frame deletion, splice variant). The others were in frame deletion (p.M4000del), in frame duplication (p.E1651_L1656dup), missense variants and a combination of 3 or 4 variants resulting in the mutant haplotype p.[I1571V;R3366H;Y3933C] or p.[I1571V; R3119H; R3366H; Y3933C] already described [10]. Each variant was found in a conserved region (in the vertebrate subphylum). Only L4650P, R2140W, L883P, P4501L and R4179H variations were classified as moderate, a similar type of amino acid change being observed in other species of vertebrate phylum. All the bioinformatic scores lead to consider all these variations as impacting the protein. Each variant was localized on the RyR1 3D structures and 26 variants were predicted to induce destabilisation by foldX software. Other 10 variants, localized on the RyR1 3D structures (in the so-called domains 5gl1 and 5taz), were predicted to induce the loss of nearby interaction. Five variants (N2342S, R4179H, V4842M, E4911K) were predicted to induce no particular change in RyR1 structure at the interaction or steric hindrance level [3,8]. The variant V4842M was associated with an additional splice variant. The N2342S has already been associated with a functional defect, the E4911K is located in the triadin interaction motif, the A4946V is located in a transmembrane domain and the R4179H is located in a calcium affinity domain [38]. Overall, 17 different missenses variants (40.5%) are located in the bridge solenoid (B-sol) domain (which contains the FKBP binding domain), 2 in the S2S3 domain, 4 in the Csol (5%), 3 in the J sol (7.5%), 5 in the NTD domain (12.5%), 6 in the pore domain (15%), 3 in the pVSD domain (7.5%), 1 in the RYR1&2 domain (2.5%), 2 in the SPRY1 (5%), 2 in the SPRY2 (5%), 1 in the SPRY3 (2.5%), 3 in the TaF domain (7.5%) and 1 in undefined regions (2.5%
Concerning the position of variants across the RYR1 gene, high prevalence (40,5%) of missense variants in BSol domain was found among the entire cohort of patients regardless of clinical and morphological phenotype. Nevertheless, no significant relationship was found between the position of variant and morphological or clinical phenotype. Interestingly, variants located in the BSol domain and the pore domain were more frequent in DuCD patients (10/19 and 9/12 respectively).

Discussion
After a full revision of histopathological and ultrastructural features of 54 muscle biopsies we have reclassified 48 RYR1-recessive patients in 4 main categories: DuCD (54%), CCD (21%), C&R (15%) and T1P+ (10%). Dusty cores represent the most frequent histopathological manifestation of RYR1-recessive myopathies. Dusty cores represent the characteristic and unifying morphological feature in this group and are characterised by the presence of reddish-purple granular material ("dusty") deposition at GT stain, corresponding to irregular areas of altered enzymatic activity at oxidative stains ("cores"), and devoid of ATPase activity. Unlike classic central cores, which show a sharply demarked ovoidal area devoid of oxidative stains, dusty cores have no sharply demarked borders, no round/ovoidal shape and no regular size. Moreover, reddish purple material deposition has not been reported in classical CCD due to dominant RYR1 mutations or in MYH7-related core disease. In fact, no dusty cores were observed among 154 muscle biopsies of dominant RYR1-and 10 of MYH7-related myopathies examined in our lab [30,32].
Dusty cores range from small cytoplasmic or subsarcolemmal spots to very large areas of sarcomeric disorganization, sometime occupying the entire muscle fibre in transversal sections, which is not a common finding in CCD. These findings are also detected by EM. The areas of sarcomeric disorganization are similar to the unstructured cores of CCD, but with some differences. First, dusty cores never span longitudinally the entire muscle fibres, but mostly range between 10 and 50 sarcomeres in length. Second, unlike the regular width of the core along the muscle fibre in central cores, the disorganized areas in dusty cores are irregular, enlarged in their central part, sometimes conferring an irregular longitudinal shape to the core, sometimes with "star-like" appearance. Third, composition of cores of DuCD group presents 2 main peculiarities compared to CCD: 1) in longitudinal sections, strands of longitudinally-smeared osmophilic material and darker fragments of short thickened Z-line are consistently observed and 2) in transversal sections, these elements are not homogenously distributed in the core area but are usually accumulated in certain regions, leading to sub-areas of dense material deposition close to sub-areas of less osmophilic material, probably corresponding to the peculiar aspect of mixed increased-decreased enzymatic activity at oxidative staining (Fig. 5). This phenomenon could explain the peculiar side-by-side of increased/decreased activity at oxidative stains and the unmatched immunofluorescence by RyR1 and DHPR, corresponding to the inhomogeneous composition of dusty cores. The presence of some areas of dusty cores intensely stained at NADH and devoid of mitochondria at EM, leads to speculate that material which accumulates inside dusty cores could, at least partially, belong to sarcoplasmic reticulum or other cytoplasmic structures. In all cases of DuCD, dusty cores were consistently associated to smaller areas of sarcomeric disorganization, but surprisingly none of our cases filled the morphological criteria to be classified as MmD. No muscle biopsy presented a multiminicore appearance by optic microscopy, and ultrastructural analysis revealed the presence of larger areas of disorganization. A possible explication of this finding could be that some small areas of disorganization could correspond to the peripheralside of dusty cores (Fig. 6). Previously reported cases of MmD in RYR1-recessive myopathies, had very similar oxidative stains to those observed in DuCD [11,18]. Furthermore, it has been elucidated that morphological features of RYR1-recessive MmD are quite different to those of SEPN-related MmD [13]. In particular, the previously reported cases of RYR1-recessive MmD were associated to type1 fibre predominance/uniformity, nuclear centralization and multiple large cores, as observed in our patients, whereas the SEPN-related MmD presents typically smaller and multiple boundless cores ("minicores") [13].
Most of DuCD patients of this study, had been previously interpreted and classified as CNM with some atypical features (CNM "plus") ( Table 2) because, at that time, the internalisation of nuclei was considered the most relevant abnormality in spite of the areas of myofibrillar disorganisation and purple material deposition, since these areas of disorganization seem to appear later. The most frequent associated features were: irregular areas of myofibrillar network alteration, myofibrillar-like aggregates, or fuchsine material deposition. It is noteworthy that all these morphological features have not been reported in BIN1, DNM2 or MTM1-related CNM [29]. Furthermore, dusty cores have never been detected in muscle biopsies of our series of CNM (7 BIN1, 29 DNM2, 32 MTM1). On the other hand, it's worthy of note that a prominent nuclear internalization and centralization is a constant finding in all cases of DuCD, sometimes manifesting as the main histopathological feature. Nevertheless, also in these cases, a careful revision, consistently detected the presence of dusty cores in variable amount of muscle fibres, confirmed by EM study. Atypical morphological feature in DuCD group, without nuclear internalization or without type1 predominance, may be rarely observed.
Finally, dusty cores should be considered the unifying morphological feature among DuCD spectrum disorder and represent the morphological signature of RYR1-recessive myopathy, even when detected in only few fibres.
Few patients had been previously reported in the first paper of Bevilacqua et al. [4] in which muscle biopsies shared histopathological lesions consistent with dusty cores. At that time, we were not aware that these lesions could represent a hallmark of RYR1-recessive myopathy, as 7 cases represented a small cohort and it could had been an atypical variability of findings in RYR1-recessive patients. Only the systematically revision of muscle biopsies in all RYR1-recessive cases, allows us to realise that more than 50% of cases had dusty cores (sometimes only in few muscle fibres) in muscle biopsy. For this reason, we considered appropriate to give a specific name (dusty cores) to these lesions, considering it a specific entity among core diseases (dusty core disease). In this context, we identified DuCD as a subgroup of congenital core myopathies and we re-classified many RYR1-recessive patients, from CNM (or other morphological diagnosis) to DuCD group.
Taken together, all these considerations, lead to suppose that the cases reported as MmD and CNM related to RYR1-recessive mutations, possibly represent specific variant of DuCD spectrum, in which dusty cores could be few, underestimated, or appeared later in life, as occurred in some of our patients. This speculation is also supported by the evidence that ocular involvement is almost systematically present in the DuCD group, as reported in cases with MmD and CNM-related RYR1-recessive myopathies [1,25], and in most severe patients. Most of severe cases belonged to the DuCD group, even if not statistically significant.
Interestingly the lowest level of RyR1 expression was observed in the DuCD group. These findings suggest a close relation between the RyR1 production and clinicopathological consequences. DuCD represent the last end morphological spectrum of RYR1-related myopathies. Possibly, a severe RYR1 haploinsufficiency, as observed in DuCD with respect to the non-DuCD groups, could impair the stability and integrity of excitation-contraction coupling at triads level. Triads replication observed in several DuCD, could be the expression of tentative compensation of an insufficient RyR1 production. Indeed triads replication and T-tubule dilation has been observed also in dihydropyridine receptor (DHPR)-related congenital myopathy [34], which is the voltage-gated L-type Ca2+ channel located on the T-tubule-SR interface with RyR1. As expected, the lowest level of RyR1 production is also associated to the most severe and early onset (<1 year) cases. The presence and stability of RyR1 protein is probably necessary not only for the sarcomeric structure maintenance, but also for the overall muscle function via efficient excitation-contraction coupling.
A last morphological consideration must be done with respect to the T1P+ group. Interestingly ocular involvement and moderate clinical presentation were frequent in T1P+. These findings could suggest a possible relation between T1P+ and DuCD, supporting the idea that T1P+ could represent an early-state of DuCD. This hypothesis is also supported by the observation of two patients (p12, p13) with the first muscle biopsy showing central nuclei with type1 uniformity and atrophy and the second one showing DuCD, suggesting an evolution of the myopathological lesions with the age. Moreover, all T1P+ cases had also some degree of sarcomeric disorganization at muscle biopsy. Nevertheless, the high level of RyR1 expression in T1P+ group makes this hypothesis controversial and leads to consider T1P+ group a separate entity. More data are warranted to confirm or refuse this hypothesis, as only 2 muscle biopsies of T1P+ group were available for the RyR1 expression study.
Genetic data in our cohort confirmed the high prevalence of missense variant in BSol domain in recessive forms of disease, even if lower than previously published data (40.5% vs 84%) [8]. Even if variants located in the BSol domain and the pore domain were more frequent in DuCD patients, no significant correlation has been found between the localization of variants and morphology. Finally, the RyR1 expression study leads to re-classify a number of genetic variants according to ACMG as follow: five class 3 "variant of unknown signification" (P4501L, Y4796H, R2458C, R1679H, M2120T) to class 4 ("probably pathogenic") and four class 4 variants (M4875T, V1207M, R3903Q, E4911K) to class 5 ("pathogenic variants"). Nevertheless, RYR1 haploinsufficiency might not be the only molecular mechanism to explain the pathophysiology: functional studies (e.g. calcium imaging test) could be helpful to evaluate the impact of each variant on RyR1 pathophysiology even if these approaches are too complex to be tested in routine diagnosis.
Finally, our results represent an emerging evidence in RYR1-recessive myopathies, but restricted to a specific study population, restricted in a monocentric study. More evidences are warranted to support our findings worldwide in other study populations and other laboratories.

Conclusions
In conclusion, dusty core is the most frequent histopathological presentation of RYR1-recessive myopathies. Dusty cores are the unifying morphological lesion among the DuCD spectrum pathology and represent the morphological hallmark for the recessive form of disease. DuCD is associated to earlier disease onset, severe clinical phenotype and lowest RYR1 expression in muscle.