Introduction
Distal myopathy is a degenerative disease of the distal muscle and is a clinically and genetically heterogeneous group [7]. Distal myopathy is usually classified as Nonaka, Miyoshi, Laing, Welander, Udd and Markesbery-Griggs distal myopathy [7]. Metabolic myopathy is a myopathy that causes distal weakness [14]. Glycogen storage disease type 9D (GSD9D, OMIM 300599) is one of the metabolic myopathies caused by the lack of decomposition of glycogen due to the deficiency of Phosphorylase kinase (PhK), so that sufficient energy for muscle contraction cannot be obtained [18,20]. GSD9D belongs to a mild metabolic disorder, and is characterized by progressive muscle weakness, exercise intolerance and cramps, and extensive death in muscle tissue [1, 17, 20, 21]. Most patients develop symptoms as adults and mainly show weakness in the distal muscles and muscle atrophy [1, 17, 20, 21].
Phosphorylase kinase is composed of four homotetramers: α (PHKA1, PHKA2), β (PHKB), γ (PHKG1, PHKG2), and δ (CALM1) [2, 4, 10, 11-13, 16, 19]. The γ subunit acts as a catalyst and is regulated by the α and β subunits, and δ is calmodulin, which gives the enzyme Ca2+ sensitivity [2]. There are two types of α subunits, liver-specific PHKA2 and muscle-specific PHKA1, and both genes exist on the X chromosome [11,20]. To date, over 100 PHKA2 mutations in HGMD (http://www.hgmd.cf.ac.uk/ac/index.php) have been reported to cause GSD9A. It has also been reported in Korea [5]. However, only 7 mutations in PHKA1 have been reported to cause GSD9D and have not been reported in Korea [1, 3, 8, 10, 15, 18, 20, 21].
We investigated patients with GSD9D who visited the hospital suspected of distal hereditary motor neuronopathies (dHMN) in another hospital. dHMN does not show any damage to the sensory nerves, and only the motor nerves have abnormalities, so it shows symptoms of weakness in the distal muscles [9]. Distal myopathy also indicates weakness in the distal muscles and muscle atrophy, making dHMN difficult to distinguish clinically from distal myopathy. Therefore, it is very necessary to distinguish dHMN from distal myopathy when the patient visits with distal muscular weakness and gait disorder. Because dHMN and distal myopathy are caused by mutations in different genes, genetic testing is a good way to distinguish these two diseases.
Therefore, we performed genetic analysis in Korean GSD9D patients and found PHKA1 mutation. In addition, the patient's clinical symptoms, pathological features, and MRI results were described. Accordingly, we attempted to compare the mutations and clinical patterns of patients enrolled in this study and several previously reported GSD9D patients.
Materials and Methods
Patients
This study enrolled a X-linked recessive Korean family with GSD9D (FC975, Fig. 1A). This study also included 500 healthy controls who had no clinical features or family history of distal myopathy, which was confirmed after careful clinical and electrophysiological examinations. Written informed consent was obtained from all participants according to the protocol approved by the institutional review board for Sungkyunkwan University, Samsung Medical Center.
Fig. 1. X-linked recessive metabolic myopathy family with novel hemizygous mutation in PHKA1. (A) Pedigrees of FC975 family. Arrows indicate probands whose DNA were used for the WES (□, ○ : unaffected; ■ : affected). (B) Chromatograms of the mutation sites by capillary sequencing method. The mutation of c.3314T>C (p.I1105T) in PHKA1 is clearly shown in the mutant alleles (arrows). (C) PHKA1 protein structure and causative mutations. The present five mutations as well as previously reported mutations are indicated below the diagram. (D) Conservation of amino acid sequences in the mutation site. Multiple protein sequence alignment revealed strong conservation of amino acid sequences at the p.I1105T mutation site among different vertebrate species (Human: NP_002628.2, Mouse: NP_032858.2, Cow: XP_002700055.1, Frog: NP_001121 510.1, Zebrafish: XP_005166576.1).
Exome sequencing and filtering of variants
Exome sequencing was performed with the Human SeqCap EZ Human Exome Library v3.0 (Roche/NimbleGen, Madison, WI), and the HiSeq2500 Genome Analyzer (Illumina, San Diego, CA) for 1 sample from the proband of FC975. The University of California, Santa Cruz assembly hg19 was the reference sequence. We selected functionally significant variants (missense, nonsense, exonic indel, and splicing site variants) from the whole exome sequencing data, and then variants registered as novel or uncommon variants (minor allele frequencies ≤ 0.01) in dbSNP150 (http://www.ncbi.nlm.nih.gov), the 1,000 Genomes project database (http:// www.1000genomes.org/) were further filtered.
In silico analysis
The Sanger sequencing method confirmed the variant using the genetic analyzer ABI3130XL (Life Technologies, Foster City, CA). The genomic evolutionary rate profiling (GERP) scores were determined by the GERP program (http://mendel. stanford.edu/SidowLab/downloads/gerp/index.html). We performed conservation analysis of the protein sequences using MEGA5, version 6.06 (http://www.megasoftware.net/). In silico analyses were done with the prediction algorithms SIFT (http://sift.jcvi.org), MUpro (http://www.ics.uci.edu/~ baldig/mutation), and PolyPhen-2 (http://genetics.bwh. harvard.edu/pph2/). The SMART program predicted domains for the PHKA1 protein (http://smart.embl.de/).
Clinical and electrophysiological examinations
Clinical information included assessments of age at onset, muscle impairments, sensory loss, deep tendon reflexes, and muscle atrophy. The muscle strength of the flexor and extensor muscles was assessed manually with the standard Medical Research Council scale. The age at onset was determined by asking patients for their ages when symptoms, including distal muscle weakness, first appeared. Neurophysiological studies were done on proband. Motor and sensory conduction studies of the median, ulnar, peroneal, tibial, and sural nerves were tested, and needle electromyography was performed in the bilateral upper and lower limb muscles. In a patient, serum CK levels were measured.
Muscle biopsy and histological examination
Histopathological analyses including immunohistochemistry of the left vastus lateralis muscles were done in the proband (Ⅲ-2). Frozen 10 mm sections were examined after staining with hematoxylin and eosin (H&E), modified Gomori trichrome (GT), nicotinamide adenine dinucleotidetetrazolium reductase (NADH-tr), and myosin adenosine triphosphatase (ATPase) preincubated at pH4.3, 4.6, and 9.4. Staining with Congo red, acid phosphatase, and periodic acid-Schiff (PAS) was performed on a frozen muscle specimen of Patient.
Lower limb MRI
The patient was examined by lower limb magnetic resonance imaging (MRI) of the hip, thigh, and calf muscles at 32 years. MRI was undertaken by a 1.5T system (Siemens Vision, Erlan-gen, Germany). The imaging was done in the axial (field of view [FOV]524–32cm, slice thickness510mm, slice gap50.5–1.0mm) and coronal planes (FOV538–40cm, slice thickness54–5mm, slice gap50.5–1.0mm). The following protocol was used for all patients: T1-weighted spin-echo (SE) (repetition time [TR]5570–650 milliseconds, echo time [TE]514–20 milliseconds, 512 matrixes), T2-weighted SE (TR52,800–4,000 milliseconds, TE596–99 milliseconds, 512 matrixes), and fat-suppressed T2-weighted SE (TR53,090– 4,900 milliseconds, TE585–99 milliseconds, 512 matrixes)
Results
Identification of novel mutation in PHKA1
Whole exome sequencing (WES) was performed to find the causative mutation of the patient. The patient was screened for mutations in the distal myopathy and metabolic myopathy genes and peripheral neuropathy genes. As a result of the screening, we did not find any causative candidates within the distal myopathy and peripheral neuropathy genes. However, within the metabolic myopathy gene, a c.3314T> C (p.I1105T) missense mutation was found in exon 31 of the PHKA1 gene (Table 3). Other mutations were frequently found in several public mutation databases (dbSNP 153, 1000 Genomes project and Exome Variant Server) and 500 healthy controls. The PHKA1 p.I1105T mutation was previously not reported in any mutation database and was not found in 500 healty controls. This mutation was located below the glycosidase family 15 domain (Fig. 1C) and this region was well conserved within various other species (Fig. 1D). The GERP score of the mutation was significantly high at 5.12, and three in silico analyzes predicted that the mutation is likely to be pathogenic. (SIFT: 0.02, Damaging; Poly Phen-2: 0.975, Probably damaging; MUpro: -0.997, Decrease the stability of protein structure).
Clinical manifestations and electrophysiologic feature
Table 1 shows the clinical symptoms of the patients included in this study and the 7 previously reported cases. The patient found no abnormalities in the developmental process and found it difficult to run at age 22. Along with this, he noticed that it was difficult to walk long. At age 27, he noticed that both toes weren't moving, and he noticed that his leg muscles were weakening. When he was 32, he visited an outside hospital with symptoms of both his toes not moving. At the time, the neuroscientist's opinion was suspected of dHMN, and he came to the hospital.
Table 1. Clinical features in patients with mutation in the PHKA1 gene
ND Not done
aCreatine kinase: normal range in our laboratory < 185IU/l
However, when examined in the neurology department of our hospital, there was no loss of sensory nerves and only impaired motor function was observed. As a result of the nerve conduction test, sensory nerves showed normal SMAP and NCS in both upper and lower extremities. Normal CMAP and NCS were found in the median, ulnar nerve, and fibula nerves of the motor nerve. However, CMAP was not measured in the left and right tibial nerves of the motor nerve (Table 2). Electromyography showed reasonable results for myopathy. The CK concentration in the blood was 577 IU / l, showing a nearly three-fold increase compared to normal subjects.
Table 2. Electrophysiological features of the patient with PHKA1 mutation
Abbreviations: NP: no potential, CMAP: compound muscle action potential, MNCV: motor nerve conduction Velocity, SNAP: sensory nerve action potential, SNCV: sensory nerve conduction velocity.
The patient's type did not show any neurological abnormalities and neither father nor mother (Fig. 1A). In addition, the parents and siblings of the patient's mother showed no abnormality based on the history examination.
Histopathological findings
Muscle histological examination revealed that angulated muscle fibers and subsarcolemmal vacuoles were observed in H&E staining, and atrophy and abnormal hypertrophy of muscle fibers were also observed (Fig. 2A). In addition, PAS positive muscle fibers were observed. These pathological features were similar to those of previously reported GSD9D patients [20]. Electron microscopy revealed glycogen granules of muscle fibers and thickened Z-band and mitochondrial swelling and degradation were also observed (Fig. 2B)
Fig. 2. Histopathologic observations of vastus lateralis muscle biopsies. (A) H&E showing angulated small size of myofibers with atrophy (arrow) and myofibers (arrowhead). (B) Electron micrograph showing a lot of glycogen (arrow) and and Z-band thickening (arrowhead)
Fatty infiltration of the calf muscle
A patient's lower extremity MRI study was performed at age 33 (Fig. 3A - Fig. 3E). In the patient, no damage was observed in the hip and thigh muscles, and only moderate damage was observed in the calf. Severe gastrocnemius muscle damage was found in the right calf muscle, and soleus and gastrocnemius muscles were equally damaged in the left muscle. In the calf, the rear compartments soleus and gastrocnemius showed severe muscle atrophy and fat infiltration, while the anterior and lateral compartment muscles were not damaged.
Fig. 3. T1-weighted magnetic resonance imaging (MRI) of patient. Coronal images of thigh (A) and calf (B). Axial image of hip (C) and thigh (D), calf (E).
Discussion
This study describes the genetic and MRI data and clinical characteristics of Korean GSD9D patients with PHKA1 mutations compared to previously reported cases. A total of seven PHKA1 mutations have been reported to cause GSD9D to date [1, 3, 8, 10, 15, 18, 20, 21]. Table 1 compares the clinical symptoms of these patients. Patients enrolled in this study had similar symptoms, such as late onset age, distal weakness, and intolerance to exercise, with previously reported GSD9D patients. However, no muscle cramp was found in our patients, which was also not found in some of the previously reported patients. The onset age of GSD9D was reported to be 17-64 years (average age of 34.6 years), excluding 2 children onset, and this patient was also 22 years old.
The PHKA1 mutation was inherited from the patient's mother. Mothers showed little clinical symptoms because they had mutations as heterozygotes. The mutation was not found in the patient's older brother and the father's and mother's siblings. It was also not found in 500 healthy controls and was also not reported in the dbSNP153 and 1,000 genome databases. Therefore, this PHKA1 mutation is considered a novel mutation. Previously reported mutations were mainly located in the glycoside hydrolase family 15 domain, but our patient mutations were located under the domain with p.E1112X and IVS31+1G>C. There were no symptom differences between the mutations in the domain and the mutations below.
In addition, this study performed MRI examination of the lower limbs for the first time in patients with GSD9D. As a result of MRI, fat infiltration was not observed in the hip and thigh, and fat infiltration was observed only in the calf, distal to the lower extremity. In the calf muscle, no damage was observed in the anterior and lateral compartments, whereas atrophy and severe fat infiltration were observed in the anterior compartment soleus and gastrocnemius muscles. This phenomenon is thought to be a singularity, so the patient showed disagreement with distal myopathy.
In summary, we wanted to present clinical features to patients with new PHKA1 mutations. We report clinically typical GSD9D symptoms and detect and report distal muscle damage by MRI. Since there are no reported mutations in Korea among the 7 cases, this study is the first to report a case of 9D type Glycogen patients in Korea. This study will expand the clinical spectrum of GSD9D with the PHKA1 mutation and suggest that it will be useful for molecular diagnosis of heterogeneous metabolic myopathy.
Table 3. Functionally significant variants found in 44 metabolic myopathy genes
ahet: heterozygous, hom: homozygous and hem: hemizygous.
bAD: autosomal dominant, AR: autosomal recessive and XR: X-linked recessive
Acknowledgements
This work was supported by grants from the National Research Foundation (2016R1A5A2007009, 2018R1A4A1024 506 and 2019R1A2C1087547) and the Korean Health Technology R&D Project, Ministry of Health & Welfare (HI14C 3484 and HI16C0426), Republic of Korea.
The Conflict of Interest Statement
The authors declare that they have no conflicts of interest with the contents of this article.
References
- Bruno, C., Manfredi, G., Andreu, A. L., Shanske, S., Krishna, S., Ilse, W. K. and DiMauro, S. 1998. A splice junction mutation in the alpha(M) gene of phosphorylase kinase in a patient with myopathy. Biochem. Biophys. Res. Commun. 249, 648-651. https://doi.org/10.1006/bbrc.1998.9211
- Brushia, R. J. and Walsh, D. A. 1999. Phosphorylase kinase: the complexity of its regulation is reflected in the complexity of its structure. Front. Biosci. 4, D618-641. https://doi.org/10.2741/Brushia
- Burwinkel, B, Hu, B., Schroers, A., Clemens, P. R., Moses, S. W., Shin, Y. S., Pongratz, D., Vorgerd, M. and Kilimann, M. W. 2003. Muscle glycogenosis with low phosphorylase kinase activity: mutations in PHKA1, PHKG1 or six other candidate genes explain only a minority of cases. Eur. J. Hum. Genet. 11, 516-526. https://doi.org/10.1038/sj.ejhg.5200996
- Burwinkel, B., Maichele, A. J, Aagenaes, O., Bakker, H. D., Lerner, A., Shin, Y. S., Strachan, J. A. and Kilimann, M. W. 1997. Autosomal glycogenosis of liver and muscle due to phosphorylase kinase deficiency is caused by mutations in the phosphorylase kinase beta subunit (PHKB). Hum. Mol. Genet. 6, 1109-1115. https://doi.org/10.1093/hmg/6.7.1109
- Choi, R., Park, H. D., Kang, B., Choi, S. Y., Ki, C. S., Lee, S. Y., Kim, J. W., Song, J. and Choe, Y. H. 2016. PHKA2 mutation spectrum in Korean patients with glycogen storage disease type IX: prevalence of deletion mutations. BMC. Med. Genet. 17, 33.
- Clemens, P. R., Yamamoto, M. and Engel, A. G. 1990. Adult phosphorylase b kinase deficiency. Ann. Neurol. 28, 529-538. https://doi.org/10.1002/ana.410280410
- Dimachkie, M. M. and Barohn, R. J. 2014. Distal myopathies. Neurol. Clin. 32, 817-842. https://doi.org/10.1016/j.ncl.2014.04.004
- Echaniz-Laguna, A., Akman, H. O., Mohr, M., Tranchant, C., Talmant-Verbist, V., Rolland, M. O. and Dimauro, S. 2010. Muscle phosphorylase B kinase deficiency revisited. Neuromuscul. Disord. 20, 125-127. https://doi.org/10.1016/j.nmd.2009.11.004
- Harding, A. E. and Thomas, P. K. 1980. The clinical features of hereditary motor and sensory neuropathy types I and II. Brain 103, 259-280. https://doi.org/10.1093/brain/103.2.259
- Harmann, B., Zander, N. F. and Kilimann, M. W. 1991. Isoform diversity of phosphorylase kinase alpha and beta subunits generated by alternative RNA splicing. J. Biol. Chem. 266, 15631-15637. https://doi.org/10.1016/S0021-9258(18)98453-X
- Hendrickx, J., Coucke, P., Dams, E., Lee, P., Odievre, M., Corbeel, L., Fernandes, J. F. and Willems, P. J. 1995. Mutations in the phosphorylase kinase gene PHKA2 are responsible for X-linked liver glycogen storage disease. Hum. Mol. Genet. 4, 77-83. https://doi.org/10.1093/hmg/4.1.77
- Koller, M., Schnyder, B. and Strehler, E. E. 1990. Structural organization of the human CaMIII calmodulin gene. Biochim. Biophys. Acta. 1087, 180-189. https://doi.org/10.1016/0167-4781(90)90203-E
-
Maichele, A. J., Burwinkel, B., Maire, I., Sovik, O. and Kilimann, M. W. 1996. Mutations in the testis/liver isoform of the phosphorylase kinase
$\gamma$ subunit (PHKG2) cause autosomal liver glycogenosis in the gsd rat and in humans. Nat. Genet. 14, 337-340. https://doi.org/10.1038/ng1196-337 - Mastaglia, F. L. and Laing, N. G. 1999. Distal myopathies: clinical and molecular diagnosis and classification. J. Neurol. Neurosurg. Psychiatry 67, 703-707. https://doi.org/10.1136/jnnp.67.6.703
- Orngreen, M. C., Schelhaas, H. J., Jeppesen, T. D., Akman, H. O., Wevers, R. A., Andersen, S. T., ter Laak, H. J., van Diggelen, O. P., DiMauro, S. and Vissing, J. 2008. Is muscle glycogenolysis impaired in X-linked phosphorylase B kinase deficiency? Neurology 70, 1876-1882. https://doi.org/10.1212/01.wnl.0000289190.66955.67
- Pegues, J. C. and Friedberg, F. 1990. Multiple mRNAs encoding human calmodulin. Biochem. Biophys. Res. Commun. 172, 1145-1149. https://doi.org/10.1016/0006-291X(90)91567-C
- Preisler, N., Orngreen, M. C., Echaniz-Laguna, A., Laforet, P., Lonsdorfer-Wolf, E., Doutreleau, S., Geny, B., Akman, H. O., Dimauro, S. and Vissing, J. 2012. Muscle phosphorylase kinase deficiency: a neutral metabolic variant or a disease? Neurology 78, 265-268. https://doi.org/10.1212/WNL.0b013e31824365f9
- Schneider, A., Davidson, J. J., Wullrich, A. and Kilimann, M. W. 1993. Phosphorylase kinase deficiency in I-strain mice is associated with a frameshift mutation in the alpha subunit muscle isoform. Nat. Genet. 5, 381-385. https://doi.org/10.1038/ng1293-381
-
Wehner, M. and Kilimann, M. W. 1995. Human cDNA encoding the muscle isoform of the phosphorylase kinase
$\gamma$ subunit (PHKG1). Hum. Genet. 96, 616-618. https://doi.org/10.1007/BF00197422 - Wehner, M., Clemens, P. R., Engel, A. G. and Kilimann, M. W. 1994. Human muscle glycogenosis due to phosphorylase kinase deficiency associated with a nonsense mutation in the muscle isoform of the alpha subunit. Hum. Mol. Genet. 3, 1983-1987. https://doi.org/10.1093/hmg/3.11.1983
- Wuyts, W., Reyniers, E., Ceuterick, C., Storm, K., de Barsy, T. and Martin, J. J. 2005. Myopathy and phosphorylase kinase deficiency caused by a mutation in the PHKA1 gene. Am. J. Med. Genet. A. 133A, 82-84. https://doi.org/10.1002/ajmg.a.30517