Corrigendum: Cinical, Metabolic, and Genetic Analysis and Follow-Up of Eight Patients With HIBCH Mutations Presenting With Leigh/Leigh-Like Syndrome
- 1Department of Neurology, Beijing Children's Hospital, Capital Medical University, National Center for Children's Health, Beijing, China
- 2Department of Research, Development of MILS International, Ishikawa, Japan
3-Hydroxyisobutyryl-CoA hydrolase (HIBCH, NM_014362.3) gene mutation can cause HIBCH deficiency, leading to Leigh/Leigh-like disease. To date, few case series have investigated the relationship between metabolites and clinical phenotypes or the effects of treatment, although 34 patients with HIBCH mutations from 27 families have been reported. The purpose of this study was to analyze the phenotypic spectrum, follow-up results, metabolites, and genotypes of patients with HIBCH deficiency presenting with Leigh/Leigh-like syndrome and explore specific metabolites related to disease diagnosis and prognosis through retrospective and longitudinal studies. Applying next-generation sequencing, we identified eight patients with HIBCH mutations from our cohort of 181 cases of genetically diagnosed Leigh/Leigh-like syndrome. Six novel HIBCH mutations were identified: c.977T>G [p.Leu326Arg], c.1036G>T [p.Val346Phe], c.750+1G>A, c.810-2A>C, c.469C>T [p.Arg157*], and c.236delC [p.Pro79Leufs*5]. The Newcastle Pediatric Mitochondrial Disease Scale (NPMDS) was employed to assess disease progression and clinical outcomes. The non-invasive approach of metabolite analysis showed that levels of some were associated with clinical phenotype severity. Five (5/7) patients presented with elevated C4-OH in dried blood spots, and the level was probably correlated with the NPMDS scores during the peak disease phase. 2,3-Dihydroxy-2-methylbutyrate in urine was elevated in six (6/7) patients and elevated S-(2-caboxypropyl)cysteamine in urine was found in three patients (3/3). The median age at initial presentation was 13 months (8–18 months), and the median follow-up was 2.3 years (range 1.3–7.2 years). We summarized and compared with all reported patients with HIBCH mutations. The most prominent clinical manifestations were developmental regression/delay, hypotonia, encephalopathy, and feeding difficulties. We administered drug and dietary treatment. During follow-up, five patients responded positively to treatment with a significant decrease in NPMDS scores. Our research is the largest case series of patients with HIBCH mutations.
Introduction
Leigh syndrome (LS), as the most common presentation of mitochondrial disorders in children, is a devastating neurodegenerative disease accompanied by symmetrical lesions in basal ganglia, brainstem, and with genetic and clinical heterogeneity. It sometimes presents in patients with normal laboratory findings, atypical neuroimaging or symptoms, but highly suggestive of LS, which is known as Leigh-like syndrome (Baertling et al., 2014). It is related to over 95 disease-associated genes and still in expanding scenery (Lake et al., 2016; Rahman and Rahman, 2018; Baldo and Vilarinho, 2020). LS has an estimated incidence of 1:40,000 newborns (Rahman et al., 1996). There is no research on the incidence of LS in Chinese. Currently, there is no standard treatment for LS and it is very challenging to performing large-scale clinical research, due to the variable phenotypes, genetic heterogeneity and early death.
The mitochondrial nuclear gene HIBCH (cytogenetic location: 2q32.2; GenBank accession number: NM_004092.3; OMIM* 610690), encodes the 3-hydroxyisobutyryl-CoA hydrolase (HIBCH) protein. It is mainly responsible for converting 3-hydroxyisobutryl-CoA to 3-hydroxyisobutyric acid in valine catabolism, as well as 3-hydroxypropanoyl-CoA to 3-hydroxypropionic acid, a hypothetical minor pathway for propionic acid catabolism, which feeds into the valine catabolic pathway (Figure 1A) (Shimomura et al., 1994; Peters et al., 2014). HIBCH mutations were first identified to cause HIBCH deficiency (HIBCHD, OMIM #250620) in 2007 (Loupatty et al., 2007), and it was reported that it could lead to Leigh/Leigh-like disease due to secondary oxidative phosphorylation (OXPHOS) defects (Figure 1B) (Ferdinandusse et al., 2013; Yamada et al., 2014; Soler-Alfonso et al., 2015; Stiles et al., 2015).
FIGURE 1. Schematic representation of the valine metabolic pathway involved in HIBCH (A), and mechanism of HIBCH deficiency with the resulting clinical phenotypes (B). SCPCM, S-(2-caboxypropyl)cysteamine; 23DH2MB, 2,3-dihydroxy-2-methylbutyrate.
The most comprehensive case series of HIBCH mutations was a report of five Turkish cases from two unrelated families (Schottmann et al., 2016). To date, just 31 HIBCH gene mutations have been reported, with 34 patients on a broad phenotypic spectrum ranging from early death (Brown et al., 1982; D’Gama et al., 2020) to just movement disorder with survival into adulthood, although the estimated incidence of HIBCH deficiency was 1 in ∼130,000 individuals (Stiles et al., 2015). The current understanding and detection of HIBCH mutations are insufficient. There is a lack of systematic research on the clinical, genotype, and biochemical characteristics of patients with HIBCH mutations. In the early stage, detection of HIBCH enzymatic activity and respiratory chain enzymes in patients’muscle or skin tissues can help clinicians make a diagnosis. However, enzymatic tests of biopsied tissue is an invasive method and sometimes yields inconsistent results (Loupatty et al., 2007; Reuter et al., 2014), limiting their clinical application. At present, genetic testing for disease diagnosis is more widespread. We propose that non-invasive metabolites should be investigated to assist early diagnosis and timely treatment.
Herein, we analyzed the genetic, metabolic, and clinical features and follow-up findings of 8 Chinese patients with HIBCH gene mutations identified by next-generation sequencing (NGS) combined with metabolic examination and who presented with Leigh/Leigh-like Syndrome.
Methods
Patients
Eight patients with identified HIBCH mutations who were from our large cohort of case series with 181 cases of genetically diagnosed Leigh/Leigh-like syndrome were recruited to participate in this study from October 2012 to December 2019.
Genetic Analyses
Due to different times of recruitment, we applied different NGS approaches: targeted panel sequencing (Patients 1 and 2) (Fang et al., 2017) and whole-exome sequencing (WES) (Patients 3, 4, 5, 6, 7, and 8). We performed Sanger sequencing to validate the identified HIBCH mutations and test the parental origin of available family members.
Briefly, genomic DNA was extracted from peripheral blood using TIANamp Blood DNA Kit (Tiangen, Beijing, China). Then 200 ng high-quality genomic DNA were utilized for library preparation. Genomic DNA was enzymatic sheared, end repaired, phosphorylation of the 5′ prime ends, a-tailing of the 3′ends and ligated to sequencing adapters, then PCR amplified following standard library preparation SOP, by using KAPA Hyper Plus Kits (Kapa Biosystems, Wilmington, MA, United States). The post-PCR libraries were captured by SureSelectXT Human All Exon V6 (Agilent, United States) or clinical exome analysis (6110 genes) (Agilent, United States), respectively. The final enriched libraries were sequenced with 2 × 150 bp on Illumina sequencers (Illumina, San Diego, CA, United States). The average sequencing depth for each sample was 149.3X (Patient 1), 174.9X (Patient 2), 148.7X (Patients 3 and 4), 133.9X (Patient 5), 130.5X (Patient 6), 132.02X (Patient 7), and 130.2X (Patient 8). Sequence data were aligned to the human genome reference (UCSC Genome Browser build hg19) using Burrows-Wheeler Aligner (BWA) (Li, 2012). Variant filtering and annotation was carried out with ANNOVAR software (Wang et al., 2010). Variants were screened as follows: 1) Preference to the variants related to the diseases, small INDEL, canonical splice sites, and nonsense variants. 2) Minor allele frequency (MAF) in normal populations <5% (except for known MAF ≥5% pathogenicity). 3) Preference to variants in the Human Gene Mutation Database (HGMD), ClinVar. 4) Preference to variants in the Online Mendelian Inheritance in Man database. Pathogenic variants were defined according to Standards and guidelines for the interpretation of sequence variants published by the American College of Medical Genetics (ACMG) in 2015 with Human Genome Variation Society nomenclature (Richards et al., 2015).
Additionally, The deletion at 2q32.2 including HIBCH gene was firstly identified by WES as the loss of heterogeneity in the Patient 8. Droplet-based digital PCR (ddPCR) was performed in the patient and his father to verify this deletion. Using QX200™ Droplet digital™ PCR systems, we performed the processes of sample dispersion, amplification and Quantification in turn. The number of droplets was obtained in the droplet reader based on the fluorescence amplitude (Hindson et al., 2011). And the data was analyzed by the software QuantaSoft™ analysis Pro 1.0.596.
Metabolite Measurements
Amino acids and acylcarnitines in dried blood spots (DBS) were analyzed by electrospray tandem mass spectrometry (LC-MS/MS) for C4-OH (composed of 3-hydroxyisobutyryl-carnitine, 3-hydroxy-butyryl-carnitine, and malonylcarnitine [C3-DC]). Samples were prepared as follows: a 3-mm diameter disc was punched in a microtiter plate, 100 μl internal standard solution was added, and the plate was covered. After gentle shaking at 45°C for 45 min, the solvent after elution was transferred into a new microtiter plate and analyzed immediately. An LCMS-8040 (Shimadzu, Japan) in positive mode was utilized for acylcarnitines and amino acids analysis. For MS/MS data acquisition, the peak including C4-OH (248.2>85) was observed with multiple reaction monitoring.
Seven patients’ urine metabolites, including 2,3-dihydroxy-2-methylbutyrate (23DH2MB), were analyzed by standard gas chromatography-mass spectrometry (GC-MS). Briefly, 100 μl urine was incubated with 40 units urease at 37°C for 15 min. After adding an internal standard (20 μg of n-heptadecanoic acid), the sample was centrifuged to deproteinize by adding 1,000 μl of ethanol, then the supernatant was evaporated. The residue was completely dried under a nitrogen stream for 3 min and derivatized with 100 μl of N, O-bis(trimetnylsilily) trifluorcetamide and 10 μl of trimethylchlorosilane at 90°C for 40 min (Zhang et al., 2000). A GCMS-QP2010 Ultra GC-MS system (Shimadzu, Kyoto, Japan) was used with an ultra Alloy capillary column (30 × 0.25 mm internal diameter with 0.25-μm film thickness, Frontier Lab, Fukushima, Japan). The temperature was programmed to increase from 60 to 350°C at 17°C/min. Next, 2 μl of derivated sample was injected in split mode. The mass chromatographic quantitation of urine metabolites including 23DH2MB was based on the fragment ion peaks area compared with the corresponding urine creatinine fragment ion area ratio.
Through retrospective re-analysis, quantitative urine screening for S-(2-caboxypropyl) cysteamine (SCPCM) was measured by LCMS/MS in urine of cases. This project was not routinely performed in the metabolic laboratory. The process included sample and label pretreatment. The SCPCM concentration in urine of 140 control children was determined by the detection and verification of an SCPCM standard, and the reference value was obtained using the SCPCM standard concentration curve.
Data Collection and Follow-Up
At the first visit, clinical data were collected retrospectively, consisting of demographic features, early developmental milestones before first symptoms, age at presentation, initial present, precipitating causes, metabolite findings, and radiographic data.
During the follow-up period, all eight patients were assessed every 6 months to 1 year at our neurology clinic. Each follow-up evaluated motor, cognition and linguistics functions; muscle strength and tone; feeding situation; growth index; recurrence; metabolite analyses; and response to treatment.
Follow-up MRI results were classified into four levels: progression, stable, evolution, and regression (Bonfante et al., 2016). The Newcastle Pediatric Mitochondrial Disease Scale (NPMDS) was adopted to assess disease progression (Phoenix et al., 2006), and each patient was scored with the NPMDS before medication (peak period) and after treatment (the last follow-up). NPMDS scores were correlated with the degree of disease burden. Higher and lower scores indicate a more severe phenotype and clinical improvement, respectively.
Review 34 Reported Cases With HIBCH Mutations
We reviewed the data of 34 published patients, of which 4 cases had insufficient clinical data (Charng et al., 2016; 336; Ceyhan-Birsoy, et al., 2019; Jou et al., 2019; Hu et al., 2020). The data of onset and current age, clinical symptoms, metabolites, neuroimaging, genotype, and natural history were analyzed to summarizes the similarities and differences between our patients and other reported cases, both with the HIBCH mutations.
Analysis of Other 173 Patients With Leigh/Leigh-Like Syndrome
We grouped confirmed genes of 181 patients with Leigh/Leigh-like syndrome according to pathomechanism including OXPHOS complexes (I, II, IV, V), mitochondrial DNA-related protein synthesis, mitochondrial homeostasis, substrate, mitochondrial cofactors, and inhibitors. The age of onset, predisposing factors, initial symptoms, main symptoms, brain imaging, laboratory examinations, and clinical outcomes of 173 patients were investigated to analyze similarities and differences between subjects with HIBCH gene mutations and other gene mutations associated with Leigh/Leigh-like syndrome.
Results
HIBCH Genetic Analyses
Ten mutations, including six novel mutations, were identified by NGS of the HIBCH gene with compound heterozygous states (Table 1), and according to the ACMG criteria, the ten variants were classified as likely pathogenic or pathogenic (Supplementary Table S1). No other candidate genes that could fully explain clinical symptoms were found. They were confirmed by Sanger sequencing. The paternal blood samples of Patient 1 and two siblings (Patients 3 and 4) were not available. The identified mutations were distributed on coding exons 4, 7, 9, 12, and 13 (Figure 2). Two novel missense mutations, c.977T>G (p.Leu326Arg) and c.1036G>T (p.Val346Phe), were absent from the Exome Aggregation Consortium (ExAC), 1,000 Genomes, gnomAD_genome, and ESP6500 databases, but conserved amino acid residues in multiple species and several in silico tools predicted these mutations as disease causing. Two novel splicing mutations (c.750+1G>A and c.810-2A>C) were predicted as likely pathogenic. One novel nonsense mutation c.469C>T (p.Arg157*) can disrupt protein function by changing the stop codon and was also absent from population databases (gnomAD_genome, ExAC, 1,000 Genomes, and ESP6500). Another novel mutation c.236delC, p.Pro79Leufs*5, was predicted to be loss of function due to a premature truncation. The HGMD includes c.1027C>G (p.His343Asp) and c.439-2A>G, which have been reported in the literature (Zhu et al., 2015; Yang et al., 2018; Xu et al., 2019). The former was identified in four patients of this cohort, and from our NGS data there were no region of homology of the variant site. c.743A>G (p.His248Arg) and c.452C>T (p.Ser151Leu) are included in Clinvar, and both have conserved residues predicted to be damaging and disease causing. In addition, Sanger sequencing and droplet digital PCR identified maternal c.1027C>G (p.His343Asp) and a paternal 868-kb heterozygous deletion from 2q32.2 including the HIBCH gene in Patient 5 homozygous for the c.1027C>G (p.His343Asp) HIBCH variant. Our identified 6 novel mutation and 31 other reported mutation sites are summarized in Figure 2.
FIGURE 2. Schematic representation the HIBCH gene structure showing disease-associated mutation locations and conservation of mutant amino acid residues with known protein domains. Notes: dashed line, mutations in our study; triangles, our novel mutations including missense, splicing, frameshift and nonsense mutations; solid line, mutations reported in the literature.
Metabolite Results and Laboratory Investigations
There were 5 (5/8) patients with DBS showing elevated C4-OH during either the peak or recovery phase. Six (6/7) had increased 23DH2MB levels in urine (Figure 3), including four cases in the peak stage and two in the recovery stage. Through retrospective analysis, we identified an increase in urine SCPCM levels in three cases. The data are shown in Table 2. There were no significant abnormalities of amino acid levels in DBS. In addition, three patients had mild elevation of lactic acid in blood, but no cases had elevated lactate in cerebrospinal fluid.
FIGURE 3. Urine organic acid profiles of seven patients in different phases; red arrows indicate abnormally high 23DH2MB peaks, and the black arrow represents a normal 23DH2MB peak. Patients 1, 2, 6, and 7 with elevated 23DH2MB levels in the peak stage. Patients 3 and 4 with elevated 23DH2MB levels in the recovery stage. Patient 8 with a normal 23DH2MB level after adopting a valine-restricted diet.
Additionally, we found that the NPMDS scores during the peak phase may be correlated with the value of C4-OH and SCPCM. We also observed that patients with lower C4-OH during the peak phase were more likely to have a positive response to treatment.
Demographic and Clinical Features
In this case series, eight patients (four male, four female) were from seven non-consanguineous Chinese families, and Patients three and four were siblings. The clinical features are summarized in Table 3. Two patients were ethnic minorities (Patient 6 Mongolian and Patient 8 Hmong). Seven were born at term with normal birth measurements, and the perinatal period was generally normal for all subjects. Early developmental milestones before disease onset were almost normal. All patients developed initial symptoms at ages ranging from 8 to 18 months at a median of 13 months. With regard to the initial symptoms, three cases developed encephalopathy, two cases had developmental regression, two cases showed developmental delay, and one case was dystonia. Six patients had precipitating causes, with five cases triggered by infection and one case by vaccination. Brain MRI scans of all cases at onset revealed T2 hyperintensity in the bilateral basal ganglia, and four cases also had brainstem lesions, similar to manifestation. In short, it was Leigh pattern imaging (Figure 4).
FIGURE 4. MRI manifestations of 2 patients (Patient 2, Patient 5) during disease onset. (a‐d): Patient 2 at age 10 months; hyperintensity of the globus pallidus, putamen, caudate nucleus, and cerebral peduncle on T2WI (a,b) and DWI (c,d). (e‐h): Patient 5 at age 13 months; hyperintensity on T2WI (e,f) and DWI (g,h) in the globus pallidus and cerebral peduncle.
The median age of diagnosis was 2.0 years (range 1.0–5.7 years). During the disease course, the frequently prominent clinical features were hypotonia (8/8), developmental delay (8/8), developmental regression (7/8), encephalopathy (6/8), feeding difficulties (5/8), dystonia (4/8), ataxia (4/8), metabolic acidosis (4/8), strabismus (4/8), nystagmus (3/8), and seizures (3/8). Video electroencephalogram (VEEG) in two cases (Patients 3 and 4) with tonic seizures showed bilateral central and apical areas with sharp and slow wave emission that was more remarkable on the right side. The VEEG of Patient 8 with focal seizures presented as continuous release of slow waves and multiple spikes in the right frontal and temporal regions. Two cases (Patients 3 and 8) developed status epilepticus during acute encephalopathy. Three (Patients 2, 3, and 8) developed respiratory failure during acute encephalopathy, and all required intubation and ventilatory support for ∼1 week, of which one case (Patient 3) had a failed extubation and underwent a second intubation for ventilator-assisted breathing. Two cases (Patients 2 and 8) had two episodes of acute encephalopathy both induced by infection in 1 year. During acute encephalopathy period, two cases (Patient 2 and 4) were both associated with thyroid dysfunction (T3, T4 and TSH all decreased) and one case (Patient 6) with hepatic dysfunction (ALT 375 U/L, reference range 0–40 U/L), and about half a year after the peak period, these indicators gradually returned to normal.
Treatment
Once the patients diagnosed, we fully informed the parents about treatments and recommended pharmacologic therapy and adopting a valine-restricted diet. The drugs consisted of antioxidants and OXPHOS complex cofactors including L-carnitine (50 mg/kg/day), coenzyme Q10 (10 mg/kg/day), thiamine (10 mg/kg/day), and riboflavin (10 mg/kg/day), additionally including some symptomatic drugs, such as levetiracetam, or baclofen. For the dietary treatment, we worked with nutritionists to develop a plan to limit valine intake with protein intake of 1.0–1.5 g/kg/d, recommended energy intake of 80 kcal/kg/d, and no restriction on carbohydrates. It was recommended to use maple glycosuria formula milk powder and supplement leucine and isoleucine.
Two cases (Patients 1 and 5) only used drug treatments, 3 (Patients 6, 7, and 8) initiated drug and diet therapy, and 3 cases (Patients 2, 3, and 4) abandoned therapy. The NPMDS scores in the peak phase and at the last follow-up (recovery phase) were shown in Figure 5. The treated patients presented with significant decreased scores.
FIGURE 5. NMPDS scores in the peak phase and last assessment. a, All the recruited patients (n = 8); b, Patients who received therapy (Patient 1,5,6,7,8; n = 5); c, Patients who gave up therapy (Patient 2,3,4; n = 3).
Follow-Up
The median duration between initial symptom onset and the last follow-up was 2.3 years (range 1.3–7.2 years). NPMDS scores during the peak phase and at the last evaluation, as well as MRI findings were showed in Table 3. The clinical conditions of three cases (Patients 1, 5, and 7) improved, while 2 (Patient 6 mild delay and Patient 8) were relatively stable; all five cases were in the treatment group. Conversely, all three cases (Patients 2, 3, and 4) in the non-treatment group progressed.
Patient 1 exhibited developmental regression following pneumonia at age 14 months but subsequently recovered. Unfortunately, at 4 years 5 months, she experienced acute encephalopathy due to a febrile infection, which again led to developmental regression, and she lost the skills of sitting and speaking. She received drug therapy thereafter and slowly recovered and gradually improved. The patient is currently 8 years 4 months old, with a normal growth index. She is independent and attends grade 3 at a normal elementary school but sometimes has trouble concentrating and wandering. She has some movement intolerance and walks with a slightly abnormal right lower limb posture. Despite having several fevers and other infections since initiating treatment, none have worsened her condition. Her last MRI scan showed significant improvement (Figure 6).
FIGURE 6. Patient 1: MRI performance in the peak and recovery phases. (a1,a2): Initial MRI obtained at age 2 years; hyperintensity on T2WI in the globus pallidus (a1), normal cerebral peduncle (a2). (b1,b2): MRI performed at 4 years 5 months during the acute stage; hyperintensity on T2WI in the globus pallidus with swelling on the left side (b1) and in the left cerebral peduncle (b2). (c1,c2): MRI performed at 6 years 1 month during the recovery stage; only slightly hyperintensity on T2WI in globus pallidus (c1), and Abnormal signal disappeared in the cerebral peduncle (c2).
For Patient 5, roseola infantum at age 13 months caused her condition to deteriorate, presenting with no head control, no language, and unable to sit unassisted. Following drug therapy, she gradually recovered. Now she is 2 years 4 months old, can understand basic sentences, walk a dozen steps unsupported, and speak 2 and 3 words.
Patient 7 has received drug and diet treatments for 6 months and has improved. The therapy included drugs of L-carnitine 500 mg/day, coenzyme Q10 100 mg/day, thiamine 100 mg/day, riboflavin 100 mg/kg/day, and limit protein intake 12 g/d, no restriction on carbohydrates. At present she is 2 years 7 months, can walk without assistance but occasionally falls, can understand complex sentences, and pronounce monosyllables.
Patient 6 have received the standard therapy of drug and diet plus baclofen (1 mg/kg/day) since the onset of encephalopathy for nearly 1 year. Despite four bouts of pneumonia during this time, his condition did not deteriorate. He is 3 years and 2 months old, with dystonia as the main manifestation. He can walk for a few steps with support, understand general language, and pronounce sequential words.
Patient 8 was initiated on drug therapy plus levetiracetam (30 mg/kg/day) since onset at age 11 months with partial responses. However, 10 months later he experienced a second encephalopathy episode followed by frequent seizures and status epilepticus. Since being placed on a valine-restricted diet, he has shown a dramatic response. Although he has been seizure-free for 5 months, there has been no improvement in development. He is 2 years 4 months old, with no head control, language, or eye contact.
The conditions of all three cases without long-term treatment (Patients 2, 3, and 4) deteriorated. They had poor response to treatment in the early stage, and all eventually stopped therapy. At the last follow-up, Patient 2 was 4 years old, with height 90 cm (<3rd percentile) and weight 8.5 kg (<3rd percentile). She has poor head control with no eye contact or language. Patient 3 is currently 8 years old. He has no head control and is unable to sit unsupported or actively grasp objects. He knows simple words but is nonverbal. Patient 4 is presently 6 years 4 months old; his clinical course and physical examination are similar to that of his older brother (Patient 3).
Six (6/8) patients underwent follow-up neuroimaging (Table 3). Brain MRI showed progression in three cases, two of which abandoned therapy, and 1 (Patient 8) case with therapy had lesion progression to the entire right cerebral hemisphere with evident brain volume loss. One case (Patient 3) had corpus callosum thinning and a lactic acid peak on magnetic resonance spectroscopy (MRS). Brain MRI improved in two cases and was relatively stable in one case, and they were all receiving therapy.
Analysis of 34 Reported and Our Cases With HIBCH Mutations
The patient's general condition, clinical, imaging, and metabolic characteristics were summarized in Table 4, and there were no significant statistical difference between our cases series and other reported cases. In our cohort, there were no deaths, probably due to the older age of onset, but 7 (22%) of the reported cases had been fatal. The main symptoms of total cases were developmental delay/regression, hypotonia, distonia, encephalopathy, and feeding difficulties. T2 hyperintensity in the bilateral basal ganglia occurs in an average of 95% of cases. Brain atrophy and brainstem lesions were also common neuroimaging abnormalities. All patients tested for SCPCM presented with elevated levels, suggesting a more specific metabolite of the disease, and more samples were needed for verification. In both groups, there was an almost equal proportion (63 vs. 65%) of cases with elevated C4-OH levels. However, the percentage of patients with elevated 23DH2MB appeared to be higher in our cases, but the difference was not statistically significant due to the small sample size.
TABLE 4. The clinical, neuroimaging and metabolic characteristics of our cases and the reported HIBCH mutation patients.
Analysis With 173 Patients With Leigh/Leigh-Like Syndrome
There were 43 other mutant genes identified in 173 cases with Leigh/Leigh-like syndrome (Supplementary Figure S1). We summarized the clinical, genetic, metabolic and imaging characteristics of the 173 cases. Due to baseline heterogeneity and large differences in sample size, we only performed statistical descriptions of the data (Supplementary Table S2). The symptoms of hypotonia, dystonia, encephalopathy, and feeding difficulties were more common in patients with HIBCH mutations. In metabolites, the elevated levels of blood C4-OH (62.5 vs. 2.3%) and urine 23DH2MB (85.7 vs. 5.8%) were significant in the HIBCH mutation group. Among the 173 patients group, 10 cases showed elevated 23DH2MB levels, including 8 with ECHS1 mutation, 1 with SURF1 mutation, and 1 with MT-ND6 mutation, and the latter 2 patients showed slightly elevated 23DH2MB levels and normal C4-OH levels. This implied that elevated levels of C4-OH combined with 23DH2MB in patients with clinically diagnosed Leigh/Leigh-like Syndrome may serve as specific metabolites for HIBCH mutation patients.
Discussion
We analyzed the clinical, metabolic, and genetic features of 8 patients with HIBCH mutations associated with Leigh/Leigh-like syndrome and reported their follow-up data providing additional natural history information. We also investigated the relationships between metabolite levels, disease severity, and clinical outcomes.
In our case series, earlier onset age was associated with a poorer prognosis. Infection or vaccination could also accelerate or exacerbate the onset or recurrence. There have been 34 patients from 27 unrelated families reported with HIBCH mutations, 4 of which have insufficient clinical data. Twenty-six (26/32) had initial symptoms within 2 years of age, including seven cases who died within an onset age of 1 year old (Brown et al., 1982; Ferdinandusse et al., 2013; Yamada et al., 2014; Peters et al., 2015; D'Gama et al., 2020). The significant clinical presentations were developmental delay/regression, hypotonia, encephalopathy, feeding difficulties and dystonia, and the last 4 symptoms were more specific in patients with Leigh/Leigh-like syndrome caused by HIBCH mutation in our cohort. Interestingly, seizures seemed to be related to phenotype severity, but it occurred infrequently with only 3 patients in this cohort. Among 10 (10/33) patients described elsewhere, most of had poor responses to treatment and worse prognoses (Loupatty et al., 2007; Ferdinandusse et al., 2013; Reuter et al., 2014; Candelo et al., 2019; Karimzadeh et al., 2019; D'Gama et al., 2020; Marti-Sanchez et al., 2020). The neuroimaging of all 8 patients showed symmetrical lesions in basal ganglia with/without brainstem. Lactate peaks on MRS (Soler-Alfonso et al., 2015; Candelo et al., 2019; D'Gama et al., 2020) and corpus callosum abnormalities (Brown et al., 1982; Peters et al., 2015; Tan et al., 2018; Candelo et al., 2019; D'Gama et al., 2020) are uncommon in cases with HIBCH mutation, but it may be associated with a more severe phenotype. Additionally, in our cohort there were 2 cases with thyroid dysfunction during the onset phase who had negative prognoses, and only 1 case were reported in the recent literature (Marti-Sanchez et al., 2020). Furthermore, one case had hepatic dysfunction as a rare presentation, and it was proposed that HIBCH might protect cells against methacrylyl-CoA toxicity (Ishigure et al., 2001). The follow-ups showed differences in the prognoses of these patients, which may be related to therapeutic interventions or their genotype.
It is difficult to correlate the genotypes and phenotypes of HIBCH deficiency. The variation sites reported to date are distributed in all exons except exons 1 and 2. Of these, exon 12 has the most mutation sites and cases, but phenotype severity varies (Ferdinandusse et al., 2013; Schottmann et al., 2016). Additionally, there have been no identified pathogenic mutations impacting the key residues of Glu120, Gly146, Glu169, or Asp177 (substrate binding). It was reported that truncating mutations tend to cause a more severe phenotype (Tan et al., 2018). One study observed a longer survival in cases with homozygous mutations located on the protein surface than in those with variants inside or near the catalytic site (Marti-Sanchez et al., 2020). Disease severity may be related to residual HIBCH enzyme activity (Yamada et al., 2014). The p.His343Asp mutation seemed to be associated with less severe clinical presentations in this cohort. To date, there has been no confirmed founder or hot spot mutation of the HIBCH gene. In our cases, c.1027C>G [p.His343Asp] was identified in half of patients, in addition to two of the four Chinese cases previously reported (Zhu et al., 2015; Xu et al., 2019), but it has not been described in other countries. Additionally, there is no region of homology from our NGS data, and it is at a extremely low frequency in population databases. We speculate that p.His343Asp may be a hot spot mutation in Chinese. The major limitation of our study is the absent detection of HIBCH enzyme activity in recruited patients. Instead, we assessed specific metabolites in blood and urine, so that patients could easily accept and assist in disease diagnosis.
HIBCH deficiency leads to the accumulation of 3-hydroxy-isobutyryl-CoA, but it was hypothesized that methacrylyl-CoA (Figure 1), as an upstream highly reactive intermediates, resulted in disease pathogenesis and neurotoxicity. The excessive generated methacrylyl-CoA can react with cysteamine to form SCPCM. 23DH2MB was derived from acryloyl-CoA, a compound homologous to methacrylyl-CoA, involving in valine metabolic pathway. In our cohort, elevated C4-OH and 23DH2MB levels are rare in Leigh/Leigh-like syndrome patients with other gene mutations. Our results indicate that C4-OH in the blood and 23DH2MB and SCPCM in the urine can be used as specific metabolites to assist the diagnosis and evaluate disease severity, although they are not completely specific.
The elevated C4-OH level is the direct result of the accumulation of 3-hydroxy-isobutyryl-CoA, according to acylcarnitine analyses in approximately 15 reported patients with confirmed HIBCH mutations. In this cohort, a higher C4-OH level in the peak phase was associated with a more severe phenotype, indicating poor prognosis. However, in current clinical laboratories, typical LCMS/MS cannot distinguish between 3-hydroxyisobutyryl-carnitine and 3-hydroxybutyryl carnitine isomers. This requires comprehensive clinical analysis. Sometimes C4-OH was partially elevated in the asymptomatic neonatal period (Stiles et al., 2015) and was easily overlooked. For example, Patient 8 had a mild increase in C4-OH in the asymptomatic neonatal period, but it did not attract enough attention to make a timely diagnosis and initiate treatment.
Probably due to different stages and different severity of the disease, 23DH2MB in our cases showed a higher positive rate compared with reported cases. 23DH2MB levels did appear to change in response to treatment. After Patient 8 adopted a valine-restricted diet for 1 month, 23DH2MB fell to the reference value, while C4-OH had not yet returned to a normal level. The first measurement of 23DH2MB was in normal urine with trace amounts (Thompson et al., 1975), and also detected in patients with propionate disorders. In the valine pathway, it is originated from acryloyl-CoA, but the exact mechanism is not clear. 23DH2MB abnormalities were also identified in a few Leigh syndrome patients with ECHS1 mutations (Peters et al., 2014; Yang and Yu, 2020). There was only 2 reported patient with an HIBCH mutation with elevated urinary 23DH2MB excretion (Ferdinandusse et al., 2013; Marti-Sanchez et al., 2020). We identified 6 patients with HIBCH mutations with elevated 23DH2MB levels.
SCPCM as a conjugate of the toxic metabolite methacrylyl-CoA, was firstly identified in HIBCH deficiency patients (Brown et al., 1982). It was an unusual amino acid in the urine of HIBCH mutation (Peters et al., 2015) and ECHS1 mutation patients, and cannot be detected in routine metabolic screening. In our case series, urinary SCPCM level correlated with both disease severity and prognosis. To date, there have been 5 patients identified with HIBCH mutations and elevated SCPCM levels in urine, including 3 cases in our study. SCPCM seems to be more specific for disease diagnosis, but more samples are needed for verification.
There is currently no consensus on HIBCH deficiency treatment approaches. A low-valine, carbohydrate-rich diet may be effective (Loupatty et al., 2007; Soler-Alfonso et al., 2015). This treatment can avoid excessive ATP production from aberrant valine metabolism. Other possible drug therapies involve the administration of antioxidants and supplementation with vitamins/cofactors, often termed nutraceuticals. Fewer than 10 cases have reported treatment outcomes that include improvement (Loupatty et al., 2007; Soler-Alfonso et al., 2015), progression and even death (Yamada et al., 2014). It is difficult to correlate therapeutic effects with gene mutation sites and clinical manifestations, which may be related to residual HIBCH enzyme activity. Patients in this cohort with the c.1027C>G [p.His343Asp] mutation seemed to responded well to treatment, similar to previous studies (Zhu et al., 2015; Xu et al., 2019).
In real clinical work, metabolites can be affected by the severity of disease and diet state, and false positive or negative results may occur. Genetic testing, especially targeted NGS (WES/panels), is the first-line auxiliary diagnostic approach. However, the targeted NGS testing is costly, and it can not identify possible pathogenic variants in deep intron regions and larger deletions/duplications. If the metabolites and clinical phenotypes are well defined and consistent, it may be recommended to directly perform the less costly approach of Sanger sequencing. Unfortunately, we did not perform enzymatic detection of HIBCH. We should combine clinical, metabolic and genetic analysis to diagnose the disease. Early diagnosis and treatment may delay disease progression of the disease. The main limitations of our study were the small sample size and the absence of formal clinical trials.
In conclusion, we identified 10 rare HIBCH gene variants in 8 patients with Leigh/Leigh-like syndrome, including 6 novel variants. The phenotype of Leigh/Leigh-like syndrome caused by HIBCH mutations was not significantly different from other gene mutations, except for more comment presentations with hypotonia, dystonia, and acute encephalopathy. The genotype-phenotype correlation requires further research. However, measurements of metabolites including C4-OH, 23DH2MB and SCPCM were relatively specific and also associated with disease severity, therapeutic effects, and possibly prognosis. The non-invasive tools of NGS and metabolite analyses should be considered first-line auxiliary diagnostic approaches. Patients with Leigh/Leigh-like syndrome with HIBCH mutations may have a positive prognosis, especially with early diagnosis and timely therapy.
Data Availability Statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Beijing Children’s Hospital, Capital Medical University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
FF, JW, and ZL Conceived and designed the research. XH, CZ, CR, and XY collected the data and performed the study. JW, ZL, and MX analyzed the data. JW wrote the manuscript. FF revised the manuscript. All authors reviewed and approved the final version.
Funding
This work was supported by the Capital Health Research and Development Fund (2018-2-2096), and National Natural Science Foundation of China(81541115).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank the patients and the many clinicians and laboratory scientists who contributed to this study. We also acknowledge MILS (Beijing) Medical Laboratory.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.605803/full#supplementary-material.
References
Baertling, F., Rodenburg, R. J., Schaper, J., Smeitink, J. A., Koopman, W. J., Mayatepek, E., et al. (2014). A guide to diagnosis and treatment of Leigh syndrome. J. Neurol. Neurosurg. Psychiatry 85 (3), 257–265. doi:10.1136/jnnp-2012-304426
Baldo, M. S., and Vilarinho, L. (2020). Molecular basis of Leigh syndrome: a current look. Orphanet J. Rare Dis. 15 (1), 31. doi:10.1186/s13023-020-1297-9
Bonfante, E., Koenig, M. K., Adejumo, R. B., Perinjelil, V., and Riascos, R. F. (2016). The neuroimaging of Leigh syndrome: case series and review of the literature. Pediatr. Radiol. 46, 443–451. doi:10.1007/s00247-015-3523-5
Brown, G. K., Hunt, S. M., Scholem, R., Fowler, K., Grimes, A., Mercer, J. F., et al. (1982). Beta-hydroxyisobutyryl coenzyme A deacylase deficiency: a defect in valine metabolism associated with physical malformations. Pediatrics 70 (4), 532–538. doi:10.1203/00006450-198212000-00014
Candelo, E., Cochard, L., Caicedo-Herrera, G., Granados, A. M., Gomez, J. F., Díaz-Ordoñez, L., et al. (2019). Syndromic progressive neurodegenerative disease of infancy caused by novel variants in HIBCH: report of two cases in Colombia. Intractable Rare Dis. Res. 8 (3), 187–193. doi:10.5582/irdr.2019.01014
Ceyhan-Birsoy, O., Murry, J. B., Machini, K., Lebo, M. S., Yu, T. W., Fayer, S., et al. (2019). Interpretation of genomic sequencing results in healthy and ill newborns: results from the BabySeq project. Am. J. Hum. Genet. 104 (1), 76–93. doi:10.1016/j.ajhg.2018.11.016
Charng, W., Karaca, E., Akdemir, Z. C., Gambin, T., Atik, M. M., Gu, S., et al. (2016). Exome sequencing in mostly consanguineous Arab families with neurologic disease provides a high potential molecular diagnosis rate. BMC Med. Genomics 9 (1), 42. doi:10.1186/s12920-016-0208-3
D’Gama, A. M., Brucker, W. J., Zhang, T., Gubbels, C. S., Ferdinandusse, S., Shi, J., et al. (2020). A phenotypically severe, biochemically “silent” case of HIBCH deficiency in a newborn diagnosed by rapid whole exome sequencing and enzymatic testing. Am. J. Med. Genet. 182 (4), 780–784. doi:10.1002/ajmg.a.61498
Fang, F., Liu, Z., Fang, H., Wu, J., Shen, D., Sun, S., et al. (2017). The clinical and genetic characteristics in children with mitochondrial disease in China. Sci. China Life Sci. 60, 746–757. doi:10.1007/s11427-017-9080-y
Ferdinandusse, S., Waterham, H. R., Heales, S. J., Brown, G. K., Hargreaves, I. P., Taanman, J. W., et al. (2013). HIBCH mutations can cause Leigh-like disease with combined deficiency of multiple mitochondrial respiratory chain enzymes and pyruvate dehydrogenase. Orphanet J. Rare Dis. 8, 188. doi:10.1186/1750-1172-8-188
Hindson, B. J., Ness, K. D., Masquelier, D. A., Belgrader, P., Heredia, N. J., Makarewicz, A. J., et al. (2011). High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal. Chem. 83 (22), 8604–8610. doi:10.1021/ac202028g
Hu, C. P., Li, X. H., Zhao, L., Shi, Y. Y., Wu, B. B., Zhou, S. Z., et al. (2020). Clinical and molecular characterization of pediatric mitochondrial disorders in south of China. Eur. J. Med. Genet. 63, 103898. doi:10.1016/j.ejmg.2020.103898
Ishigure, K., Shimomura, Y., Murakami, T., Kaneko, T., Takeda, S., Inoue, S., et al. (2001). Human liver disease decreases methacrylyl-CoA hydratase and beta-hydroxyisobutyryl-CoA hydrolase activities in valine catabolism. Clin. Chim. Acta 312 (1–2), 115–121. doi:10.1016/s0009-8981(01)00597-6
Jou, C., Ortigoza-Escobar, J. D., O'Callaghan, M. M., Nascimento, A., Darling, A., Pias-Peleteiro, L., et al. (2019). Muscle involvement in a large cohort of pediatric patients with genetic diagnosis of mitochondrial disease. J. Clin. Med. 8 (1), 68. doi:10.3390/jcm8010068
Karimzadeh, P., Saberi, M., Sheidaee, K., Nourbakhsh, M., and Keramatipour, M. (2019). 3-Hydroxyisobutyryl-CoA hydrolase deficiency in an Iranian child with novel HIBCH compound heterozygous mutations. Clin. Case Rep. 7 (2), 375–380. doi:10.1002/ccr3.1998
Lake, N. J., Compton, A. G., Rahman, S., and Thorburn, D. R. (2016). Leigh syndrome: one disorder, more than 75 monogenic causes. Ann. Neurol. 79 (2), 190–203. doi:10.1002/ana.24551
Li, H. (2012). Exploring single-sample SNP and INDEL calling with whole-genome de novo assembly. Bioinformatics 28, 1838–1844. doi:10.1093/bioinformatics/bts280
Loupatty, F. J., Clayton, P. T., Ruiter, J. P., Ofman, R., Ijlst, L., Brown, G. K., et al. (2007). Mutations in the gene encoding 3-hydroxyisobutyryl-CoA hydrolase results in progressive infantile neurodegeneration. Am. J. Hum. Genet. 80 (1), 195–199. doi:10.1086/510725
Marti-Sanchez, L., Baide-Mairena, H., Marcé-Grau, A., Pons, R., Skouma, A., López-Laso, E., et al. (2020). Delineating the neurological phenotype in children with defects in the ECHS1 or HIBCH gene. J. Inherit. Metab. Dis. 1–14. doi:10.1002/jimd.12288
Peters, H., Buck, N., Wanders, R., Ruiter, J., Waterham, H., Koster, J., et al. (2014). ECHS1 mutations in Leigh disease: a new inborn error of metabolism affecting valine metabolism. Brain 137, 2903–2908. doi:10.1093/brain/awu216
Peters, H., Ferdinandusse, S., Ruiter, J. P., Wanders, R. J., Boneh, A., and Pitt, J. (2015). Metabolite studies in HIBCH and ECHS1 defects: implications for screening. Mol. Genet. Metab. 115 (4), 168–173. doi:10.1016/j.ymgme.2015.06.008
Phoenix, C., Schaefer, A. M., Elson, J. L., Morava, E., Bugiani, M., Uziel, G., et al. (2006). A scale to monitor progression and treatment of mitochondrial disease in children.Neuromuscul. Disord. 16 (12), 814–820. doi:10.1016/j.nmd.2006.08.006
Rahman, J., and Rahman, S. (2018). Mitochondrial medicine in the omics era. Lancet 391 (10139), 2560–2574. doi:10.1016/S0140-6736(18)30727-X
Rahman, S., Blok, R. B., Dahl, H. H., Danks, D. M., Kirby, D. M., Chow, C. W., et al. (1996). Leigh syndrome: clinical features and biochemical and DNA abnormalities. Ann. Neurol. 39 (3), 343–351. doi:10.1002/ana.410390311
Reuter, M. S., Sass, J. O., Leis, T., Köhler, J., Mayr, J. A., Feichtinger, R. G., et al. (2014). HIBCH deficiency in a patient with phenotypic characteristics of mitochondrial disorders. Am. J. Med. Genet. A. 164A (12), 3162–3169. doi:10.1002/ajmg.a.36766
Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 17 (5), 405–424. doi:10.1038/gim.2015.30
Schottmann, G., Sarpong, A., Lorenz, C., Weinhold, N., Gill, E., Teschner, L., et al. (2016). A movement disorder with dystonia and ataxia caused by a mutation in the HIBCH gene. Mov Disord. 31 (11), 1733–1739. doi:10.1002/mds.26704
Shimomura, Y., Murakami, T., Fujitsuka, N., Nakai, N., Sato, Y., Sugiyama, S., et al. (1994). Purification and partial characterization of 3-hydroxyisobutyryl-coenzyme A hydrolase of rat liver. J. Biol. Chem. 269 (19), 14248–14253. doi:10.1016/s0021-9258(17)36781-9
Soler-Alfonso, C., Enns, G. M., Koenig, M. K., Saavedra, H., Bonfante-Mejia, E., and Northrup, H. (2015). Identification of HIBCH gene mutations causing autosomal recessive Leigh syndrome: a gene involved in valine metabolism. Pediatr. Neurol. 52 (3), 361–365. doi:10.1016/j.pediatrneurol.2014.10.023
Stiles, A. R., Ferdinandusse, S., Besse, A., Appadurai, V., Leydiker, K. B., Cambray-Forker, E. J., et al. (2015). Successful diagnosis of HIBCH deficiency from exome sequencing and positive retrospective analysis of newborn screening cards in two siblings presenting with Leigh’s disease. Mol. Genet. Metab. 115 (4), 161–167. doi:10.1016/j.ymgme.2015.05.008
Tan, H., Chen, X., Lv, W., Linpeng, S., Liang, D., and Wu, L. (2018). Truncating mutations of HIBCH tend to cause severe phenotypes in cases with HIBCH deficiency: a case report and brief literature review. J. Hum. Genet. 63 (7), 851–855. doi:10.1038/s10038-018-0461-8
Thompson, J. A., Markey, S. P., and Fennessey, P. V. (1975). Gas-chromatographic/mass-spectrometric identification and quantitation of tetronic and deoxytetronic acids in urine from normal adults and neonates. Clin. Chem. 21 (13), 1892–1898. doi:10.1093/clinchem/21.13.1892
Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 38 (16), e164. doi:10.1093/nar/gkq603
Xu, Y., Zhang, J., Yu, K., Feng, F., Sun, X., Li, C., et al. (2019). A therapeutic regimen for 3-hydroxyisobutyryl-CoA hydrolase deficiency with exercise-induced dystonia. Eur. J. Paediatr. Neurol. 23 (5), 755–759. doi:10.1016/j.ejpn.2017.11.004
Yamada, K., Naiki, M., Hoshino, S., Kitaura, Y., Kondo, Y., Nomura, N., et al. (2014). Clinical and biochemical characterization of 3-hydroxyisobutyryl-CoA hydrolase (HIBCH) deficiency that causes Leigh-like disease and ketoacidosis. Mol. Genet. Metab. Rep. 1, 455–460. doi:10.1016/j.ymgmr.2014.10.003
Yang, H. Y., Wu, L. W., Deng, X. L., Yin, F., and Yang, L. F. (2018). Diagnosis and treatment of 3-hydroxyisobutyryl-CoA hydrolase deficiency: a case report and literature review. Zhongguo Dang Dai Er Ke Za Zhi 20 (8), 647–651. doi:10.7499/j.issn.1008-8830.2018.08.009
Yang, H., and Yu, D. (2020). Clinical, biochemical and metabolic characterization of patients with short-chain enoyl-CoA hydratase(ECHS1) deficiency: two case reports and the review of the literature. BMC Pediatr. 20 (1), 50. doi:10.1186/s12887-020-1947-z
Zhang, C., Xu, K., Dave, U. P., Wang, Y., and Matsumoto, I. (2000). Inborn errors of metabolism discovered in Asian department of pediatrics and mental retardation research center. J. Chromatogr. B Biomed. Sci. Appl. 746 (1), 41–49. doi:10.1016/s0378-4347(00)00087-6
Keywords: HIBCH gene, Leigh/Leigh-like syndrome, C4-OH, 2,3-dihydroxy-2-methylbutyrate (23DH2MB), mitochondrial disorders, children
Citation: Wang J, Liu Z, Xu M, Han X, Ren C, Yang X, Zhang C and Fang F (2021) Cinical, Metabolic, and Genetic Analysis and Follow-Up of Eight Patients With HIBCH Mutations Presenting With Leigh/Leigh-Like Syndrome. Front. Pharmacol. 12:605803. doi: 10.3389/fphar.2021.605803
Received: 16 September 2020; Accepted: 22 January 2021;
Published: 08 March 2021.
Edited by:
Ruth Roberts, ApconiX, United KingdomReviewed by:
Charlotte L. Alston, Wellcome Trust Center for Mitochondrial Research (WT), United KingdomCheryl D. Cropp, Samford University, United States
Copyright © 2021 Wang, Liu, Xu, Han, Ren, Yang, Zhang and Fang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fang Fang, ZmFuZ2ZhbmdAYmNoLmNvbS5jbg==