- Letter to the Editor
- Open Access
Core-rod myopathy due to a novel mutation in BTB/POZ domain of KBTBD13 manifesting as late onset LGMD
Acta Neuropathologica Communications volume 6, Article number: 94 (2018)
Few genes (RYR1, NEB, ACTA1, CFL2, KBTBD13) have been associated with core-rod congenital myopathies . KBTBD13 belongs to the Kelch-repeat super-family of proteins and is implicated in the ubiquitination pathway. Dominant mutations in KBTBD13 have been associated with a peculiar form of core-rod myopathy (NEM6) so far . Childhood onset, slowly progressive proximal muscle weakness with characteristic slowness of movements and combination of nemaline rods, irregular shaped cores and unusual type2 fibres hypotrophy at muscle biopsy, were the main characteristics shared in all the affected members of the four KBTBD13 families reported in the literature .
We report on a 65 years old patient, of Sardinian origin, with atypical clinical and morphological presentation of NEM6 due to a novel mutation in KBTBD13 gene.
Patient underwent a complete clinical examination, cardiac and respiratory evaluation, biological screening for myopathies, muscle imaging studies, muscle biopsies and molecular analysis for muscle diseases.
Three open muscle biopsies were obtained from deltoid, paraspinal and quadriceps muscles. Sampling site was based on neuroimaging study. For conventional histochemical techniques 8 μ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. Frozen muscle samples for immunohistochemical analysis were analysed for membrane proteins study including antibodies against dystrophin (N-ter, Rod, C-ter), sarcoglycans, α-dystroglycan (IIH6, VIA4), laminin α2, dysferlin and caveolin-3. Antibodies were visualized using innunoperoxidase techniques. Immunofluorescence study was performed for KBTBD13 with rabbit polyclonal antibodies (ab110507 abcam 1:50). Ultrastructural study by electron microscopy (EM) was performed on small muscle specimens from the last biopsy, fixed with glutaraldehyde (2.5%, pH 7.4), post fixed with osmium tetroxide (2%), dehydrated and embedded in resin. Ultra-thin sections from 6 blocks small blocks were stained with uranyl acetate and lead citrate. 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).
Whole body CT-scan and lower limb muscle MRI were obtained following a previously described protocol according to international consensus recommendation .
Western blot (WB) study
Skeletal muscle of patient and age-matched control muscle specimens were isolated and snap frozen in liquid nitrogen. Tissue lysates were prepared on ice in specific buffer (50 mM Tris–HCl, pH 7.4, 150 mM NaCl, 1 mM EDTA, 1% Triton X-100 or 1% NP40, 1 mM PMSF), and protease inhibitor mixture. Lysates were cleared by centrifugation at 15000×g. Protein concentrations of the cleared lysates were determined by BCA protein assay kit (Thermo Fisher Scientific). Proteins (40 μg) were separated by SDS-PAGE, loaded in a 4–12% denaturing gel, transferred to Immobilon-P membrane and probed with rabbit polyclonal against KBTBD13 (ab110507 abcam 1:500) and beta-tubulin (sigma 1:1000). Reactive bands were detected using Lite Ablot Extend Long Lasting Chemiluminescent Substrate (Euroclone, Pero, Italy). Densitometry analysis was performed using Quantity One software (BioRad, Hercules, CA, USA).
Next Generation Sequencing (NGS) was performed by a customized gene panel including the following 95 genes: NEB, MYO18B, ACTA1, TPM2, TPM3, KBTBD13, KLHL40, KLHL41, LMOD3, TNNT1, CFL2, DNM2, MTM1, SPEG1, CCDC78, BIN1, RYR1, MYH2, MYH7, SEPN1, STIM1, ORAI1, PGAM2, MEGF10, MTMR14, TRIM32, FHL1, HACD1, C-term TTN, VMA21, SCN4A, CLCN1, VCP, GNE, KLHL9, HSPB8, ADSSL1, SQSTM1, MATR3, TIA1, BAG3, CRYAB, DES, FLNC, LDB3 (ZASP), HNRPDL, LMNA, FKRP, DYSF, CAPN3, FKTN, TOR1AIP1, GTDC2 (POMGNT2), GMPPB, LAMA2, COL6A1, COL6A2, COL6A3, ANO5, PYGM, CPT2, CAV3, DMD, EMD, CHKB, POMT1, POMT2, POMGNT1, FCMD, LARGE, B3GALNT2, B4GAT1 (B3GNT1), COL4A1, DAG1, ISPD, POMK (SGK196), TMEM5, SGCG, SGCA, SGCB, SGCD, DNAJB6, MYOT, TNPO3, TCAP, RAPSN, MUSK, DOK7, LAMB2, COLQ, GFPT1, CHRNE, DPAGT1, CHAT, SMN1. Target enrichment was carried out using Nextera Rapid Capture Custom Enrichment Kit (Illumina, San Diego, California, USA) with probe design by using the Design Studio software. Paired-end sequencing was performed using an Illumina MiSeq with a sequencing depth of 100×. The Illumina VariantStudio data analysis software was used to annotate the variants. Sanger sequencing was used for POGLUT1 gene screening and validation of KBTBD13 variant.
The pCMV6-AC-GFP mammalian expression vector with subcloned wild type Human KBTBD13 cDNA (kindly donated by Dr.Nyamkhishig Sambuughin of Department of Military and Emergency Medicine, Uniformed Services University, Bethesda, USA) was used as the template DNA for producing the KBTBD13G67R variant by exploiting the QuikChange XL site-directed mutagenesis kit (Agilent Technologies). Specifically, the substitution was performed by PCR amplification using the following synthetic primers carrying the desired mutation G67RF-5′-ctg cag gtg ctg cgc cgc gac cgg ccg gcg ctg-3′ and G67RR-5′-cag cgc cgg ccg gtc gcg gcg cag cac ctg cag-3′. The right insertion of the mutation was then confirmed by DNA sequencing. Plasmidic DNA was obtained by using NucleoBond Xtra (Macherey-Nagel, Duren, Germany). Then, HeLa were transfected by using Lipofectamine 3000 (ThermoFisher Scientific, Waltham, MA, USA) and plasmidic DNA (1000 ng), following manufacturer’s instructions.
Cell culture, transfection and fluorescence microscopy
HeLa cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% foetal bovine serum (FBS) and 50 μg/ml gentamycin. Cells were maintained in 5% CO2 atmosphere at 37 °C. All reagents for cell culture were from ThermoFisher Scientific (Waltham, MA, USA). HeLa cells were seeded at 21 × 105 cells/well in 6-well plates 24 h prior transfection, and transfected by using Lipofectamine 3000 (ThermoFisher Scientific), according to manufacturer’s instructions. Cells were incubated for 48 h, fixed with 4% paraformaldehyde (Sigma-Aldrich), and viewed by using a confocal microscopy (Leica, Wetzlar, Germany).
The structure of KBTBD13 protein has yet to be determined experimentally but structural data are available for other BTB-containing proteins, in some cases belonging to the BTB-Kelch family. We used Homology Modelling using SWISS-MODEL  to predict the structure of the BTB domain only of WT KBTBD13. Since the BTB domain is known to be directly involved in the formation of Kelch proteins dimers, we calculated the putative structure of the dimer formed by the BTB domains of KBTBD13, using as template the crystal structure that of the dimeric BTB domain in SPOP:Cul3 assembly (3HQI.pdb). The structure of the G67R mutant was also calculated by Homolgy Modelling. The electrostatic potential was calculated using the APBS software  and the input file required to run APBS was prepared using PDB2PQR  on the calculated WT and mutant protein structure.
The patient was born after uneventful pregnancy from healthy non-consanguineous parents. He had a 71-years-old brother with muscle hypertrophy and late onset gait disturbances who refused to be examined and tested for biological and molecular analysis. Motor milestones were normally acquired. He started to complain about slowly progressive muscular weakness of upper and lower limbs since the age of 50. First clinical examination at the age of 60 revealed generalised muscle hypertrophy, slowness of movements conferring a slow-motion bearing, proximal upper and lower limb weakness sparing distal muscles, Gowers’ sign and waddling gait (Fig. 1a). He had moderate respiratory involvement (FVC 69%) with normal cardiac evaluation. Electromyography (EMG) showed myopathic changes, and Creatine-Kinase (CK) ranged between 250 and 350 UI/L (normal values < 180 UI/L).
Muscle imaging data
Whole body muscle CT scan performed at 60 years showed a moderate to severe fibro-fatty replacement in deltoids, subscapularis and infraspinatus, serratus, latissimus dorsi, cervical paraspinalis and rectus abdominis, iliopsoas, tensor fascia lata and gluteus maximus with sparing of leg muscles. Interestingly, several muscles presented an irregular pattern of fibro-fatty replacement. A peculiar pattern involving the internal regions and sparing the external areas (“from inside-to-outside”) was observed in thigh muscles (Fig. 1b). Moreover, rectus femoris presented the “central shadow” sign and gluteus maximus showed a “zebra” pattern of fibro-fatty replacement, characterized by consecutive strips of replaced and spared muscle tissue.
The first muscle biopsy, performed in deltoid muscle at 60 years, showed fibre size variability and splitting with type1 predominance. Small clusters of subsarcolemmal rods were detected in less than 5% of muscle fibres (Fig. 1c); immunohistochemistry for membrane proteins was normal. The second muscle biopsy was performed in paraspinal muscles at the age of 61, obtained during spinal surgery because of acute lumbar disc herniation, and showed subsarcolemmal and cytoplasmic rods in several fibres (~ 30%), type1 fibre uniformity and irregular areas of myofibrillar disorganization, sometimes resembling cores (Fig. 1d,). The third muscle biopsy performed on the quadriceps muscle at 65 years, clearly showed the presence of classic eccentric cores in some fibres (< 10%) (Fig. 2a, b, c), fibre size variability and splitting, type1 uniformity and abundant lobulated fibres with irregular myofibrillar network at oxidative stains (Fig. 2c). Ultrastructural study of the last biopsy confirmed the presence of cores with abundant electrodense longitudinally smeared material (Fig. 2d, e) and detected small nemaline rods in atrophic fibres (Fig. 2f). Areas of Z-line smearing and thickening were also observed outside cores (Fig. 2g, h). Protein aggregates were not detected.
NGS allowed to identify a novel heterozygous c.199G > C (p.Gly67Arg) mutation in the BTB/POZ domain of KBTBD13 gene (Fig. 3a). Other 94 genes included in the panel were ruled. The G67R variant falls in highly conserved protein domain with MAF < 0.0001 (genomAD Database) and is predicted to be pathogenetic by bioinformatic software.
Structural and functional investigations
Transfected HeLa cells with KBTBD13G67R and KBTBD13wt vectors did not show any difference in cellular expression or localization of KBTBD13 (Fig. 3,C). Interestingly, we observed a reduction of 48.3% of the KBTBD13 expression by WB on patient’s muscle biopsy compared to healthy controls (Fig. 3b). Finally, molecular modelling of KBTBD13 protein showed that residue 67 is very close to the BTB dimerization interface and the Gly-to-Arg substitution impacts on a wide area of surface electrostatic potential around the residue 67 (Fig. 3d).
Although some clinical features were quite similar to those previously described in KBTBD13 patients, as the characteristic slowness of movements, the late age at onset in our patient represents a novel insight for KBTBD13-related myopathy, since all the previously reported patients manifested a childhood onset. Moreover, our patient showed unusual LGMD-like clinical phenotype characterized by waddling gait and muscle hypertrophy. Muscle MRI findings were different from those previously reported in KBTBD13 patients and surprisingly, the peculiar “from inside-to outside” pattern of fibro-fatty replacement was very similar to what has been recently described in POGLUT1 muscular dystrophy . For this reason, we also ruled out possible mutations in POGLUT1 gene by Sanger analysis. The “central shadow” sign of the rectus femoris has been described in COL6-related myopathies, but contrary to the “inside-to-outside” pattern observed in our patient, COL6 myopathies show a typical peripheral involvement of thigh muscles . Moreover, the full preservation of lower leg muscles observed in our patient could help in differential diagnosis with both COL6 and POGLUT1-related myopathies which typically show soleus involvement. To the best of our knowledge, “zebra” pattern of gluteus maximus has never been reported in other myopathies.
Concerning morphology, nemaline rods were constantly observed in all muscle biopsies performed in our patient, although sometimes detected only in few fibres or just by EM study. Only the third muscle biopsy revealed sharply demarked cores, never detected in previous reports, whereas the second one showed irregular-shaped areas devoid oxidative stains as described in previously cases. Interestingly, none of the biopsies showed the type2 fibre hypotrophy previously described in all reported KBTBD13 patients. Congenital myopathies are frequently characterised by prevalence/uniformity of type1 fibres, most of which are hypotrophic, and sometimes could manifest as Congenital Fibre Type Disproportion (CFTD), characterised by marked hypotrophy of all type1 fibres mixed with normal type2 fibres . Conversely, type 2 hypotrophy represent an uncommon finding among congenital myopathies, and more broadly among inherited myopathies. Type2 atrophy is frequently observed in acquired conditions as cachexia, aging or corticosteroid-induced myopathy , and among inherited myopathies it has been rarely reported . Among congenital myopathies and particularly nemaline myopathies, the type2 hypotrophy has been described only in KBTBD13 and represented a specific signature of NEM6. Conversely, our patient showed type1 prevalence or uniformity, as commonly observed in many congenital myopathies.
The clinical and histological variability observed in our patient compared to the previously reported NEM6 patients, could be related to the different localization of the KBTBD13 mutation. In fact, all the three mutations previously reported are located in Kelch domain, whereas G67R mutation detected in our patient is in the BTB/POZ domain. Possibly, the surprising reduction of KBTBD13 protein by WB analysis might be explained considering the effect of G67R mutation. KBTBD13 proteins assemble into dimers through self-association of the BTB domains. Each BTB domain interacts with the N-terminus of one Cul3 molecule forming a RING ubiquitin ligase (Cul3-RL) complex capable of ubiquitination . Our molecular modelling revealed that location of Gly67 is very close to the BTB dimerization interface; on the contrary, the location of residue 67 is rather far from Cul3-binding interface. Since the mutation described concerns the substitution of a small, hydrophobic amino acid (Gly) with another featuring a rather bulky side chain, positively charged at physiological pH (Arg), it is likely that the mutation alters the surface properties of the protein, especially the surface electrostatic potential. Electrostatic interactions across protein interfaces are known to be crucial to the formation of multimers and complexes that often represent the biologically active conformation of proteins. Our computational investigations also suggest that the G67R mutation alters the surface electrostatic potential, probably impacting the tendency of BTB dimerization in vivo. Nevertheless, this study concerns only one patient and more evidences are needed to conclusively state the effect of mutations in BTB/POZ domain of KBTBD13 gene.
In conclusion our findings broaden the clinical, histological and genetical spectrum of the ultra-rare KBTBD13-related myopathy and enlarge pathophysiological knowledges about nemaline myopathies, showing that NEM6 can have an adult onset, may show muscle hypertrophy with peculiar “inside-to-outside” MRI/CT pattern, and may lack of type2-hypotrophy at muscle biopsy. Different domains involved of KBTBD13 could lead to the phenotypic variability in NEM6.
Arnold K, Bordoli L, Kopp J, Schwede T (2006) The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 22:195–201. https://doi.org/10.1093/bioinformatics/bti770
Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA (2001) Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci U S A 98:10037–10041. https://doi.org/10.1073/pnas.181342398
Ciciliot S, Rossi AC, Dyar KA, Blaauw B, Schiaffino S (2013) Muscle type and fiber type specificity in muscle wasting. Int J Biochem Cell Biol 45:2191–2199. https://doi.org/10.1016/j.biocel.2013.05.016
Clarke NF (2011) Congenital Fiber-type disproportion. Semin Pediatr Neurol 18:264–271. https://doi.org/10.1016/j.spen.2011.10.008
Dolinsky TJ, Czodrowski P, Li H, Nielsen JE, Jensen JH, Klebe G, Baker NA (2007) PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucleic Acids Res 35:W522–W525. https://doi.org/10.1093/nar/gkm276
Garibaldi M, Tasca G, Diaz-Manera J, Ottaviani P, Laschena F, Pantoli D, Gerevini S, Fiorillo C, Maggi L, Tasca E, D’Amico A, Musumeci O, Toscano A, Bruno C, Massa R, Angelini C, Bertini E, Antonini G, Pennisi EM (2017) Muscle MRI in neutral lipid storage disease (NLSD). J Neurol 264:1334–1342. https://doi.org/10.1007/s00415-017-8498-8
Jungbluth H, Treves S, Zorzato F, Sarkozy A, Ochala J, Sewry C, Phadke R, Gautel M, Muntoni F (2018) Congenital myopathies: disorders of excitation–contraction coupling and muscle contraction. Nat Rev Neurol 14:151–167. https://doi.org/10.1038/nrneurol.2017.191
Mercuri E, Clements E, Offiah A, Pichiecchio A, Vasco G, Bianco F, Berardinelli A, Manzur A, Pane M, Messina S, Gualandi F, Ricci E, Rutherford M, Muntoni F (2010) Muscle magnetic resonance imaging involvement in muscular dystrophies with rigidity of the spine. Ann Neurol 67:201–208. https://doi.org/10.1002/ana.21846
Oldfors A (2007) Hereditary myosin myopathies. Neuromuscul Disord 17:355–367. https://doi.org/10.1016/j.nmd.2007.02.008
Olivé M, Goldfarb LG, Lee H-S, Odgerel Z, Blokhin A, Gonzalez-Mera L, Moreno D, Laing NG, Sambuughin N (2010) Nemaline myopathy type 6: clinical and myopathological features. Muscle Nerve 42:901–907. https://doi.org/10.1002/mus.21788
Sambuughin N, Swietnicki W, Techtmann S, Matrosova V, Wallace T, Goldfarb L, Maynard E (2012) KBTBD13 interacts with Cullin 3 to form a functional ubiquitin ligase. Biochem Biophys Res Commun 421:743–749. https://doi.org/10.1016/j.bbrc.2012.04.074
Sambuughin N, Yau KS, Olivé M, Duff RM, Bayarsaikhan M, Lu S, Gonzalez-Mera L, Sivadorai P, Nowak KJ, Ravenscroft G, Mastaglia FL, North KN, Ilkovski B, Kremer H, Lammens M, van Engelen BGM, Fabian V, Lamont P, Davis MR, Laing NG, Goldfarb LG (2010) Dominant mutations in KBTBD13, a member of the BTB/Kelch family, cause Nemaline myopathy with cores. Am J Hum Genet 87:842–847. https://doi.org/10.1016/j.ajhg.2010.10.020
Servián-Morilla E, Takeuchi H, Lee TV, Clarimon J, Mavillard F, Area-Gómez E, Rivas E, Nieto-González JL, Rivero MC, Cabrera-Serrano M, Gómez-Sánchez L, Martínez-López JA, Estrada B, Márquez C, Morgado Y, Suárez-Calvet X, Pita G, Bigot A, Gallardo E, Fernández-Chacón R, Hirano M, Haltiwanger RS, Jafar-Nejad H, Paradas C (2016) A POGLUT1 mutation causes a muscular dystrophy with reduced notch signaling and satellite cell loss. EMBO Molecular Medicine 8:1289–1309. https://doi.org/10.15252/emmm.201505815
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.