Compound heterozygous mutations in electron transfer flavoprotein dehydrogenase identified in a young Chinese woman with late-onset glutaric aciduria type II

Background Glutaric aciduria type II (GA II) is an autosomal recessive disorder affecting fatty acid and amino acid metabolism. The late-onset form of GA II disorder is almost exclusively associated with mutations in the electron transfer flavoprotein dehydrogenase (ETFDH) gene. Till now, the clinical features of late-onset GA II vary widely and pose a great challenge for diagnosis. The aim of the current study is to characterize the clinical phenotypes and genetic basis of a late-onset GAII patient. Methods In this study, we described the clinical and biochemical manifestations of a 23-year-old female Chinese patient with late-onset GA II, and performed genomic DNA-based PCR amplifications and sequence analysis of ETFDH gene of the whole pedigree. We also used in-silicon tools to analyze the mutation and evaluated the pathogenicity of the mutation according to the criteria proposed by American College of Medical Genetics and Genomics (ACMG). Results The muscle biopsy of this patient revealed lipid storage myopathy. Blood biochemical test and urine organic acid analyses were consistent with GA II. Direct sequence analysis of the ETFDH gene (NM_004453) revealed compound heterozygous mutations: c.250G > A (p.A84T) on exon 3 and c.920C > G (p.S307C) on exon 8. Both mutations were classified as “pathogenic” according to ACMG criteria. Conclusions In conclusion, our study described the phenotype and genotype of a late-onset GA II patient, reiterating the importance of ETFDH gene screening in these patients. Electronic supplementary material The online version of this article (10.1186/s12944-017-0576-5) contains supplementary material, which is available to authorized users.


Background
Glutaric aciduria type II (GA II), also known as multiple Acyl-CoA dehydrogenase deficiency (MADD), is a rare autosomal recessive disorder of fatty acid, amino acid and choline metabolism caused by a defect in the alpha or beta subunit of the mitochondrial electron transfer flavoprotein (ETFA, ETFB) protein or the electron transfer flavoprotein dehydrogenase (ETFDH) protein [1]. The clinical manifestation of GA II is heterogeneous, from severe neonatal form to mild late-onset form. Children with the severe form usually show a fatal course in the neonatal period. In the mild form, the age of onset and the symptoms are variable with intermittent episodes, challenging the clinical diagnosis [2]. These patients often present with fluctuating proximal muscle weakness or episodes of rhabdomyolysis. Other manifestations include sensory neuropathy, loss of tendon reflex, fatty liver, etc. [2]. Current clinical diagnosis is mainly based on tandem mass spectrometry findings of elevated plasma levels of short, medium, and long chain acylcarnitines. However, in late-onset patients, the elevation of acylcarnitine levels may be mild and atypical, or detectable only during an acute episode. Under these circumstances, genetic screening is the key to establish a definitive diagnosis.
It has long been well-known that the mild form of GAII has exceptionally good response to riboflavin supplementation, although its true mechanism remains to be clarified. Recent work also suggests that most lateonset GA II patients with good response to riboflavin have mutations in the ETFDH gene. The ETFDH gene is located at 4q32.1, with a total of 13 exons. It encodes the electron transfer flavoprotein: ubiqionone oxidoreductase (ETF: QO) [3]. As a component in mitochondrial respiratory chain, ETF: QO forms a short pathway with ETF to transfer electrons from mitochondrial flavoprotein dehydrogenases to the ubiquinone pool [4].
So far, more than 190 mutations in the ETFDH gene, including point mutations, nonsense mutations, insertions, deletions, splicing mutations, have been reported according to the Human Gene Mutation Database (HGMD) [5]. The genotype and phenotype in ETFDH-mutant GAII patients were previously reported to be correlated. Nonsense mutations leading to truncation or missense mutations affecting protein stability in both alleles were generally associated with a severe phenotype [6,7]. However, according to a recent cohort study, in the late-onset mild forms, the correlation between genotype and phenotype (including age of onset, disease severity and response to treatment) is not well-established [8].
Herein, we describe the clinical presentation and course of a 23-year-old female who was diagnosed as late-onset GA II. We identified heterozygous mutations in ETFDH gene, including a p.A84T, located at the flavin adenine dinucleotide (FAD)-binding domain and a p.S307C mutation, located within the ubiquinone (UQ)binding domain of ETF dehydrogenase.

Case description
A 23-year-old female was admitted to our hospital because of exercise intolerance and general muscle weakness, especially in her lower limbs, for 3 months. The symptom of weakness in her neck and proximal limb muscles was noted to have aggravated in last 1 month. She had difficulty walking long distances, gradually she was unable to raise her head and arms away from bed. Two weeks before admission, she experienced intermittent vomiting and weight loss. Her prenatal history and early development were both normal, and her life history was unremarkable. Physical and neurological examinations showed neck and proximal muscle weakness (manual muscle testing (MMT) score 3/5 in neck, 2/5 in lower limbs, 4/5 in upper arm).
Besides her clinical history, other clinical details including blood biochemical examinations, blood acylcarnitine analysis, urine organic acids spectrum, muscle magnetic resonance imaging (MRI), CT scan of the abdomen, electrocardiogram (EEG) were also collected. Muscle biopsy was performed on the right biceps brachii and immunohistochemistry was performed using frozen sections. The patient was treated with riboflavin (120 mg/day) and carnitine supplements (1 g/day). Two days after riboflavin and carnitine treatment, she experienced a severe respiratory infection and developed rapidly progressive quadriparesis with acute respiratory failure. Despite immediate tracheal intubation, the patient succumbed to sudden heart arrest, which showed no response to chest compressions, intravenous adrenalin, or electric defibrillation.
The proband and her family were Han Chinese and lived in southern China. Her parents were not consanguineous and had no myopathic symptom. The younger brother of the proband, a 15-year old boy, had no history of neuromuscular disease. After written informed consent was obtained, blood samples of all family members were collected for genetic testing.

Mutation analysis
Genomic DNA-based PCR amplifications and sequence analysis of ETFDH gene of the whole pedigree was performed. Briefly, genomic DNA from peripheral blood leukocytes was isolated by proteinase K digestion and phenol/chloroform extraction. All of the exons and exonintron boundaries of the ETFDH (GenBank NM_004453) gene were amplified by PCR with high fidelity KOD-plus polymerase from the hyperthermophilic Archaeon Thermococcus kodakaraensis KOD1 (TOYOBO, Japan) using primers listed in Table 1. The amplified PCR products were purified from agarose gel using QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) and sequenced via the ABI3730XL sequencer (Applied Biosystems, CA, USA). Mutations were confirmed using control DNA from 200 unaffected Chinese individuals after obtaining a written informed consent.

Bioinformatic analysis
Two in silicon tools were used for the analysis of mutation pathogenicity. The machine learning method predicts variants according to Bayesian methods (PolyPhen2) [9], or mathematical operations (SIFT) [10]. Multiple sequence alignment was performed using ClustalX [11]. The X-ray crystal structure of pig ETF:QO (Protein Data Bank entry 2GMH, resolution 2.5 Å) [12] was used as the template to construct human wild type and S307C mutational ETFDH structure models. The sequence of human ETFDH was retrieved from UniProt (http://www.uniprot.org/uniprot/ Q16134). The 3D models of human ETFDH complexed with substrates (ubiquinone, FAD, and 4Fe4S) were generated using MODELLER (version 9.12) [13].
Muscle MRI of lower limbs showed high signal intensity areas in the bilateral lower limb muscles on T2-weighted MR imaging, indicating diffuse muscle injury (Fig. 1a). CT scans of the abdomen indicated severe lipid accumulation in the liver (fatty liver) as suggested by significantly lower density compared with that of the spleen (Fig. 1b).
Histochemical staining of muscle frozen sections revealed vacuolar myopathy, with numerous fine vacuoles found mainly in type 1 fibers (Fig. 1c). The Oil Red O stain showed diffuse lipid droplets accumulation, predominantly in type 1 fibers (Fig. 1d). Nicotinamide adenine dinucleotide (NADH) staining indicated the myofibrillar network was disorganized in type 1 fibers, and no excessive glycogen content was observed on periodic acid-schiff (PAS) staining ( Fig. 1f and g). The pathological findings were compatible with lipid storage myopathy. Histological findings of a control case were shown in Additional file 1: Figure S1.

Molecular studies
Direct sequencing of all 13 exons of ETFDH gene revealed compound heterozygous mutations in the proband: c.250G > A (p.A84T) on exon 3 and c.920C > G (p.S307C) on exon 8 (Fig. 2b). Her father carried a c.250G > A (p.A84T) on exon 3, while her mother carried a c.920C > G (p.S307C) on exon 8. Her younger brother, a 17-year old boy, inherited c.920C > G (p.S307C) from his mother. The missense mutation c.920C > G (p.S307C) was not identified in a screening of 200 control chromosomes in Chinese individuals. No mutation was identified in ETFA or ETFB gene.

Bioinformatic and structural analysis
Both silicon tools, namely PolyPhen-2 and SIFT, supported that p.S307C of ETFDH was a deleterious mutation ( Table 3). The human and pig ETFDH proteins share 95% sequence identity, which makes homology modeling very reasonable. The results showed a minimal difference between the wild type and mutant structures, even in the ubiquinone binding pocket for p.S307C mutation (Fig. 2d). It should be mentioned that, the structure of amino acid serine (S) was very close to cysteine (C). The only difference in the chemical structure was the Sulfur atom instead of the Oxygen atom.  [8].

Discussion
According to American College of Medical Genetics and Genomics (ACMG), the criterion for classification of pathogenic or likely pathogenic variants is weighted as very strong (PVS1), strong (PS1-4), moderate (PM1-6), or supporting (PP1-5) [5]. Several independent studies have consistently reported that the c.250G > A (p.A84T) mutation was the most common mutation in GAII patients from Southern China, as well as those from Taiwan, Hong Kong, Thailand and Singapore (PS1)   [14,15]. Consistently, the current proband and her father, both carrying the c.250G > A (p.A84T) mutation, lived in Zhejiang province, geographically classified as Southern China. The p.A84T mutation was previously reported to cause reduced ETFDH protein expression in vivo (PS3) [16]. The prevalence of p.A84T in GAII patients was significantly increased compared with the prevalence in East Asians (PS4) [15]. The p.A84T mutation itself was a mutational hot spot and located in a critical and well-established functional domain (FAD domain) (PM1). Thus, it was classified as "pathogenic" according to the criteria proposed by ACMG (3PS + 1 PM) [5].
In contrast, the c.920C > G (p.S307C) mutation was relatively new. The p.S307C mutation was previously reported and associated with GA II (PS1) [17,18]. The p.S307C mutation was located in a critical and wellestablished functional domain (UQ-binding domain) (PM1) [18]. Furthermore, this variant was neither found in Exome Aggregation Consortium (ExAC) nor 1000 Genomes Project (PM2). Taking co-segregation into account in the pedigree, the p.S307C mutation was in trans with a pathogenic variant p.A84T (PM3). The deleterious effect was supported by multiple lines of computational evidence such as PolyPhen2 and SIFT (PP3). The phenotype of the patient is highly consistent with the phenotype caused by the ETHDH gene (PP4). Based on these findings, it should also be classified as "pathogenic" according to the criteria proposed by ACMG(1PS + 3 PM + 2PP) [5].
The c.920C > G (p.S307C) mutation was previously described in only two patients. Both patients were Han Chinese, whose age of onset were much older than our patient [17,18]. The features of these three patients carrying the p.S307C mutation were summarized in Additional file 2: Table S1. One patient initially manifested as Bent spine syndrome. Both patients had limb and lip numbness, and EMG showed no sensory nerve action potential, suggesting severe sensory neuropathy besides muscle weakness [17,18]. However, in the current case, progressive muscle weakness was the major symptom, while sensory nerve action potential was preserved. CT scan of the abdomen in our patient indicated severe lipid accumulation in the liver. Some GA II patients were previously reported to have fatty liver [19,20]. In one case, liver lipid accumulation relieved after riboflavin supplementation [20]. Since the patient died unexpectedly, we were unable to see an effect of riboflavin supplementation on fatty liver.
The ETF:QO protein acts in conjunction with ETF to govern electron transport from several flavoprotein acyl-CoA dehydrogenases to the main respiratory chain in mitochondria [12]. Residues A84 and S307 are located in the flavin adenine dinucleotide (FAD)-binding domain and the ubiquinone (UQ)-binding domain of ETF dehydrogenase, respectively. Both of the codons are highly evolutionarily conserved. Previous studies have suggested that A84T mutation could affect the interaction between FAD and ETF:QO, thus, could reduce the efficiency of electron transfer [21]. S307, together with Y304, G305 and G306, belong to the β-sheet 11, which wrap around the C5 methyl, O4 carbonyl, and C3 methoxy groups of UQ. In particular, the O4 atom of UQ makes hydrogen bonds to the backbone nitrogen of G306 and carbonyl oxygen of G305 [12]. Based on the structural studies, it is reasonable to assume that the p.S307C mutation, which is beside G306 and G305, might disturb the formation of the hydrogen bonds between UQ and G306/G305. Further in vivo studies are still needed to confirm the pathogenicity of S307C mutation.

Conclusion
In summary, we described the phenotype and genotype of a late-onset GA II patient. Clinicians should be alert to young patients with unexplained muscle weakness and high CK levels. ETFDH gene screening is necessary for establishing the diagnosis and genetic counseling if muscle biopsy showed lipid accumulation.

Additional files
Additional file 1: Figure S1. Histological findings of a control case. research. AK contributed to bioinformatic analysis. LYC carried out the molecular genetic study. YHW and ZQL designed the study, revised and prepared the final version of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate This study was conducted in accordance with the Declaration of Helsinki and approved by the local ethics committee (Tongji Hospital of Tongji University). Written informed consent was obtained from each subject participating in the study after full explanation of the purpose and all procedures.

Consent for publication
Not applicable.