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).

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 late-onset 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 exon-intron 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.

Table 1

Primers and amplicons for detecting ETFDH gene mutation

Name

Sequence of primers

Amplicon(bp)

ETFDH-Exon1-F

CAATCGATCTCGAAGGGCACT

659 bp

ETFDH--Exon1-R

TTAATTCTACCAACTGGGGCAAC

ETFDH--Exon2-F

GTTAAGCCTGACATGAGCTAAATTG

537 bp

ETFDH--Exon2-R

TAATTGTCGTGTTGTTGACTGAAGG

ETFDH--Exon3-F

TGTGCAAAACACAGGGAGAATTTC

512 bp

ETFDH--Exon3-R

AGCCTGGGCAACAAGAGTGAAA

ETFDH--Exon4-F

AGGAGAAACACTTGAACCCAGGA

502 bp

ETFDH--Exon4-R

AGTAAGTCCTTCAAATATCTGGGTCTC

ETFDH--Exon5-F

TGAAAGTGTGACCATCAATGTAGCA

437 bp

ETFDH--Exon5-R

TTGATGGGGGTACAACATGAGGA

ETFDH--Exon6-F

CTTTCTCATCAAGGTTGTGGCAT

405 bp

ETFDH--Exon6-R

TTGGAGAGATGGGGTTTCACTTT

ETFDH--Exon7-F

CCATTGGCAGAGGAGCTGGAT

683 bp

ETFDH--Exon7-R

CAACACTTGATACGACAATTCCAAC

ETFDH--Exon8-F

CACCGTGCCCAGCCTTCTT

477 bp

ETFDH--Exon8-R

AATGACACAGCAGCCATAAGCACT

ETFDH--Exon9-F

CCAGAGCACAAGGATTTCTTAATATG

693 bp

ETFDH--Exon9-R

ATTTGTTGTAGAGATGGGGTTTCG

ETFDH—Exon10-F

ATTTTTCAGCCTTTCCCTACAGC

551 bp

ETFDH—Exon10-R

GCCTACCAAAGTGCTGGGATTAC

ETFDH—Exon11-F

CCAACCTGGGTGACAGAGCAA

580 bp

ETFDH—Exon11-R

CTTCAACAATTTTGAGGGAAATGCT

ETFDH—Exon12-F

TGAGGGCTAGTCATATTTCTTTGGT

444 bp

ETFDH—Exon12-R

TTTCTAAGGAATGGAAGGAGATACAG

ETFDH—Exon13-F

TGAGAGGATGACTGTGAATAAGGGA

689 bp

ETFDH—Exon13-R

GAACTGAAGAGGTAGGAAGATGCTG

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].

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.

Fig. 2

Identification of the ETFDH gene mutations and the ETFDH protein structural analysis. a Pedigree of the family based on direct sequencing of the ETFDH gene. b Electropherogram of the proband. c.250G > A (p.A84T) mutation at exon 3 and c.920C > G (p.S307C) mutation at exon 8 of the ETFDH gene were identified. c The 3D models of human ETFDH structures. d The wild type and S307C mutational structures showed cyan and green cartoons. Oxygen atoms are shown in red and nitrogen atoms in blue. Carbon atoms of wild type, S307C mutation and ubiquinone are shown in cyan, green and white, respectively. The phosphorus atom of mutational Cys307 is shown in yellow-orange. Crucial residues in the binding site of ubiquinone are shown as stick and labeled. Figures were generated by PyMol

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.

GA II is mainly caused by homozygous or compound heterozygous mutations in ETFA, ETFB, or ETFDH genes. Mutations in ETFA and ETFB are usually associated with the neonatal forms, whereas ETFDH mutations often present as late-onset forms. For the late-onset form, the clinical manifestation varies considerably. So far, most GA II patients were reported from East Asia, especially mainland China. A recent study reported the clinical features and ETFDH mutation spectrum of 90 unrelated patients pooled from two major centers in China. Several hot spots were described: c.250G > A (p.A84T), c.770A > G (p.Y257C) and c.1227A > C (p.L409F) with a frequency of 12.2%, 15.0% and 12.2% in Chinese, respectively [8].

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 well-established 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.

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.

Acknowledgments

We thank Dr. Pu Xia (Fudan University) for critical reading of the manuscript. We acknowledge Dr. Wenhua Zhu (Department of Neurology, Huashan Hospital, Fudan University) for providing pictures of muscle pathology. We are deeply grateful to the family of the patient for their participation.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 81400834, 81500589, 81570799).

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

Authors’ contributions

YX analyzed the patient data and drafted the manuscript. YZ contributed to the case report section and revised the manuscript. KQZ contributed to data interpretation and critical revision of the manuscript. LL coordinated the 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.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.