Abstract

Context
Abnormal growth and short stature are observed in patients with mitochondrial disease, but it is unclear whether there is a relationship between final adult height and disease severity.
Objective
To determine whether patients with genetically confirmed mitochondrial disease are shorter than their peers and whether stature is related to disease severity.
Design
Analysis of final adult height in relation to disease severity as determined by the Newcastle Mitochondrial Disease Adult Scale (NMDAS).
Setting
UK Mitochondrial Disease Patient Cohort (Mito Cohort).
Patients
575 patients were identified with recorded height, weight, and molecular genetic diagnosis of mitochondrial disease within the Mito Cohort.
Main Outcome Measures
Adult height, body mass index (BMI), and their association with genetic subgroup and disease severity.
Results
Adults with mitochondrial disease were short, with a mean height of −0.49 SD (95% CI, −0.58 to −0.39; n = 575) compared with UK reference data. Patients were overweight, with a BMI SD of 0.52 (95% CI, 0.37 to 0.67; n = 472). The most common genetic subgroup (m.3243A>G mutation) had a height SD of −0.70 (95% CI, −0.85 to −0.54; n = 234) and a BMI SD of 0.12 (95% CI, −0.10 to 0.34; n = 212). NMDAS scores were negatively correlated with height SD (r = −0.25; 95% CI, −0.33 to −0.17; P < 0.001, n = 533). Rate of disease progression also correlated negatively with adult height (P < 0.001).
Conclusion
Final height in mitochondrial disease reflects disease severity and rate of disease progression. Mitochondrial dysfunction and associated subclinical comorbidities affect growth plate physiology.
Mitochondrial disease has an estimated prevalence of 1 in 4300 adults (1). Disease can present with multisystem involvement, but neurologic features often predominate (2). Poor growth and short stature are well recognized in affected patients, yet the precise mechanisms resulting in growth failure are poorly understood.
Mitochondria are essential organelles found in all nucleated cells. Mitochondria have multiple roles in cellular metabolism, including a fundamental role in ATP synthesis through the process of oxidative phosphorylation. Mitochondria contain their own DNA [mitochondrial DNA (mtDNA)] encoding 37 genes. Thousands of copies of mtDNA may exist in an individual cell. One of several unique features of the mitochondrial genome is heteroplasmy, the coexistence of both mutated and wild type mtDNA within the same cell or tissue. Mutated mtDNA is often well tolerated at low levels of heteroplasmy, but once a specific threshold is exceeded, cellular and organ functions are disrupted (2). There are >1100 mitochondrial proteins encoded by nuclear genes, which are involved in the replication, translation, and maintenance of mtDNA, as well as other aspects of mitochondrial function (3). Mitochondrial disease can therefore arise as a consequence of a mutation in either mtDNA or the nuclear genome. In either case, the adverse impact on ATP generation is detrimental to organs with high energy demands such as the brain, heart, liver, and skeletal muscle.
The Newcastle Mitochondrial Disease Adult Scale (NMDAS) is a validated clinical scale used to assess multisystem involvement and disease progression. NMDAS explores all aspects of mitochondrial disease, focusing on current function and current clinical assessment, including system-specific involvement and quality of life (4). NMDAS is routinely performed on patients attending outpatient review by their neurologists in Newcastle Upon Tyne every 12 to 24 months. A low NMDAS score suggests the patient is asymptomatic or has low disease burden, whereas a high score reflects severe multisystem involvement.
The increase in height observed during childhood and adolescence is under the control of numerous factors, with a key physiological endpoint being growth and multiplication of growth plate chondrocytes (5). Mitochondrial dysfunction therefore has potential to undermine the growth process. Wolny et al. (6) reported that children with mitochondrial disease are about −1.97 SD shorter than their unaffected peers, with an associated BMI reduction of −1.07 SD (7). The m.3243A>G mutation in MT-TL1 is the most common pathogenic mutation of mitochondrial disease and is commonly associated with short stature; 73% of a Chinese cohort of pediatric patients with this mutation were identified as being short (8). A similar observation was reported in a Japanese cohort, although a Dutch study found that short stature was observed only in female patients with mitochondrial disease (8–10). Patients with mitochondrial disease may not become symptomatic until adulthood, but they may still be shorter than expected, which suggests a subclinical effect the on growth plate metabolism.
We hypothesized that adults with mitochondrial disease would be shorter than their peers and that their stature would be related to disease severity and genetic subgroup.

Subjects and Methods

We undertook a retrospective analysis of height, weight, and molecular genetic data from the UK Mitochondrial Disease Patient Cohort (Mito Cohort), collected from inception in March 2009 to February 2017. All patients provided informed consent when joining the Mito Cohort, a database of patients with genetically determined or biochemical evidence of mitochondrial disease. Data analysis was approved by the Mitochondrial Disease Cohort Oversight Committee. Height and body mass index (BMI) were converted to SDs based on UK reference data (British 1990 growth reference centiles) (11). The overall disease severity of individual patients was evaluated with the NMDAS, with the most recent score used in the analysis (412).
We determined the relationship between final height and NMDAS scores in all adult patients where data were available. Additional analysis was conducted in specific groups of patients. They included patients with m.3243A>G, single large-scale mtDNA deletion, and nuclear gene defects. This analysis allowed us to determine whether there was any difference in the association of height, BMI, and NMDAS between genotypes. The m.3243A>G mutation is the single most common cause of mitochondrial disease (1). In patients with m.3243A>G mutation, age-adjusted blood heteroplasmy levels were also analyzed (13). In the single, large-scale mtDNA deletion group, we performed additional analysis by looking at Kearns-Sayre syndrome (KSS), chronic progressive external ophthalmoplegia and ptosis (CPEO), and CPEO plus (CPEO+) because those with KSS are known to have early-onset disease, whereas those with CPEO and CPEO+ present later in life (2).
Normality testing was performed before statistical analysis, with nonnormally distributed data subject to square root transformation before parametric testing. The square root of NMDAS score was therefore included in all described analyses (13). The independent samples t test was used to compare means between two groups. ANOVA was used to compare means between multiple groups; subsequent comparisons between groups were evaluated with Bonferroni-corrected t tests. The association of two continuous variables was tested via Pearson correlation. All means are displayed with 95% CIs unless otherwise stated. The statistical significance threshold was set at P < 0.05 for all tests.
We determined the relationship between final adult height and BMI between disease manifestations as defined by questions in the NMDAS score with Pearson correlation. To quantify the association between height SD, BMI SD, age, and progression of disease burden, we used the Linear and Nonlinear Mixed Effects Models package to perform a linear mixed effects analysis of the relationship between heteroplasmy level, age, and NMDAS scaled score (14). As fixed effects we entered age and height SD or BMI SD into the model. As random effects we had by-subject random slopes for the effect of age. Two-sided P values were obtained from the t distribution reported in Linear and Nonlinear Mixed Effects Models.
In our cohort we are not aware of any patients with mitochondrial disease who have been treated with GH. Short stature due to mitochondrial disease is not an indication for GH treatment in the United Kingdom (15).

Results

A total of 891 participants were enrolled in the Mito Cohort in Newcastle, of whom 797 patients were aged ≥18 years by January 2017. A total of 222 patients were excluded because of a lack of recorded height (n = 143), no established genetic cause (n = 78), or a height data entry error (n = 1). A total of 575 adults aged 18 to 85 years at the time of their height measurement, were identified with genetic evidence of mitochondrial disease and growth data available (Table 1). Of these, 38.3% were male (n = 220) and 61.3% female (n = 355).
Table 1.
Cohort Demographics
Genetic ID Genetic Diagnosis Number of Patients Male Female Mean Age at Height Measurement Age Range
 
Mean Height SD 95% CI of the Mean for Height SD
 
Number of Patients With BMI Mean BMI SD 95% CI of the Mean for BMI SD
 
All cohort 575 220 355 46 18–85 −0.49 −0.58 to −0.39 472 0.52 0.37 to 0.67 
Primary mtDNA mutations           
 m.14709T>C 11 10 48 26–64 0.30 −0.57 to 1.17 −0.24 −3.00 to 2.53 
 m.3243A>G 234 95 139 42 18–79 −0.70 −0.85 to −0.55 212 0.12 −0.10 to 0.34 
 m.8344A>G 40 13 27 45 20–69 −0.41 −0.72 to −0.11 28 0.66 −0.01 to 1.34 
 Single large-scale mtDNA deletion 90 34 56 47 18–85 −0.27 −0.52 to −0.03 78 0.68 0.33 to 1.03 
 Othera 54 20 34 42 20–77 −0.33 −0.59 to −0.06 40 0.72 0.23 to 1.20 
Nuclear DNA mutations           
 OPA1 17 12 46 18–72 −0.28 −0.91 to 0.35 13 2.11 1.56 to 2.67 
 POLG 27 11 16 49 19–80 −0.09 −0.50 to 0.31 19 1.00 −0.03 to 1.96 
 RRM2B 17 12 59 31–76 −0.30 −0.88 to 0.29 14 0.97 0.34 to 1.60 
 TWINK 40 15 25 58 18–82 −0.42 −0.73 to −0.11 32 1.32 0.79 to 1.85 
 Other multiple deletion 27 12 15 63 45–77 −0.69 −1.22 to −0.16 15 1.27 0.64 to 1.91 
 Otherb 18 39 18–81 −0.66 −1.24 to −0.09 14 0.00 −1.60 to 1.61 
a
Other primary mtDNA mutations include m.11778G>A, m.12147G>A, m.12148T>C, m.12258C>A, m.12271T>C, m.12283A>G, m.12315G>A, m.12320A>G, m.13051A>G, m.13094T>C, m.13513G>A, m.14674T>C, m.15699G>C, m.16002T>C, m.16023G>A, m.1624C>T, m.3243A>T, m.3365T>C, m.3460G>A, m.4175G>A, m.4267A>G, m.4298G>A, m.4300A>G, m.5543T>C, m.5650A>G, m.5690A>G, m.618T>G, m.7497G>A, m.7587T>C, m.7989T>C, m.8839G>C, m.8851T>C, m.8993T>C, m.8993T>G, m.9176T>C, and m.9185T>C.
b
Other nuclear DNA mutations include the following: AFG3L2 encodes AFG3-like protein 2, ETFDH encodes electron transfer flavoprotein dehydrogenase, GFER encodes growth factor augmenter of liver regeneration, MTFMT encodes mitochondrial methionyl-tRNA formyltransferase, OPA1 encodes dynamin-like 120-kDa protein, PDHencodes pyruvate dehydrogenase complex, POLG encodes the α-subunit of mitochondrial polymerase-γ, RRM2B encodes ribonucleotide reductase regulatory TP53 inducible subunit M2B, SDHA encodes succinate dehydrogenase complex subunit A, SPG7 encodes paraplegin, TRIT1encodes tRNA isopentenyltransferase, TWNK (previously also known as PEO1) encodes twinkle helicase, TYMP encodes thymidine phosphorylase, and YARS2 encodes tryosyl-tRNA synthetase.
The height recorded at ≥18 years of age was assumed to be the patient’s final height. In 38 patients, the date of height measurement was not documented; for these patients final height was assumed to be their consent date for entry in the cohort, provided they were >18 years old.
Subsequent analysis was performed on all patients and specific heterogeneous groups: m.3242A>G (n = 234), nuclear gene defects (n = 146), and single large-scale mtDNA deletions (n = 90). We omitted other mtDNA mutations (n = 107) because the sample size for each genotype was insufficient.

Age

The mean age at height measurement of patients recruited to the cohort was 46.1 years (95% CI, 44.8 to 47.5; n = 538). The date height was taken was not available for 38 patients (Table 1). The mean age at height measurement in the m.3243A>G group was 42.4 years (95% CI, 40.6 to 44.3; n = 226), in the single large-scale deletion group it was 47.3 years (95% CI, 44.0 to 50.6; n = 84), and in the nuclear gene defect group it was 53.5 years (95% CI, 50.5 to 56.5; n = 129).
There was a significant difference between mean age at height measurement between the genetic subgroups, as determined by one-way ANOVA [F(2, 436) = 21.38, P < 0.001]. A multiple comparison, using Bonferroni corrected t tests, revealed that the mean age at height measurement was significantly lower in the m.3243A>G group than in the nuclear gene defect group (mean difference = −11.09; 95% CI, −15.17 to −7.00; P < 0.001) and the single large scale deletion group (mean difference = −4.85; 95% CI, −9.58 to −0.12; P = 0.042). There was also a significant difference in the age at height measurement between the single large-scale deletion group and nuclear groups (mean difference = −6.23; 95% CI, −11.42 to −1.05; P = 0.012).

Height

There was no significant difference in the mean height SD score (SDS) between men (mean = −0.42; SD 1.13; 95% CI, −0.57 to −0.27; n = 220) and women [mean = −0.52; SD 1.18; 95% CI, −0.65 to −0.40; n = 335; t(573) = 1.02; P = 0.307; 95% CI, −0.94 to 0.29]. The 575 adults with mitochondrial disease were shorter than the general population, with a mean SD for height of −0.49 (95% CI, −0.58 to −0.39) (Table 1Fig. 1). Furthermore, 34.6% (n = 199) of adults had a height >1 SD below the mean, and 10% (n = 57) of adults had a height >2 SD below the mean.
Figure 1.
Individual boxplots of height SD and BMI SD for genotypes. Boxplots showing the median, first quartile, third quartile, and range for (a) height SD and (b) BMI SD for the cohort and each genotype.
The mean height SD of the m.3243A>G subgroup was −0.70 (95% CI, −0.85 to −0.54; n = 234), in the single large-scale deletion group it was −0.27 (95% CI, −0.52 to −0.03; n = 90), and in the nuclear group it was −0.41 (95% CI, −0.59 to −0.22; n = 146). There was a significant difference in mean final height SDs between genetic subgroups as determined by one-way ANOVA [F(2, 467) = 5.37, P = 0.005]. Multiple comparison (Bonferroni corrected) revealed that the mean final height SD was significantly lower in the m.3243A>G group than in the single large-scale deletion group (mean difference = −0.43; 95% CI, −0.77 to −0.08; P = 0.011). There was no difference between the m.3243A>G and nuclear groups (mean difference = −0.29; 95% CI, −0.59 to 0.01; P = 0.057). There was no significant difference between the single large-scale deletion group and the nuclear deletion group (mean difference = 0.13; 95% CI, −0.51 to 0.24; P = 1.000). In the m.3243A>G group, 42% (n = 102) had a height of <1 SD below the UK reference data, and 15% (n = 36) had a height <2 SD below the UK reference data.
There is limited data on when patients with mitochondrial disease attain their final height and whether this timing is delayed. In our cohort, final adult height was not significantly lower in the patients whose heights were recorded between 18 and 25 years of age (mean = −0.73; SD 1.32; 95% CI, −1.09 to −0.38; n = 55) compared with those where height was recorded at >25 years of age [mean = −0.46; SD 1.10; 95% CI, −0.56 to −0.36; n = 482; t(535) = −1.69; P = 0.092; 95% CI, −0.59 to −0.64]. This equates to a mean height of 159.2 cm for women and 172.5 cm for men <25 years and mean height of 160.9 cm for women and 174.1 cm for men >25 years.
In an analysis of the single large-scale mtDNA deletion, patients with KSS had a mean height SD of −1.41 (95% CI, −2.31 to −0.51; n = 11), patients with CPEO had a mean height SD of 0.33 (95% CI, −0.15 to −0.81; n = 20), and those with CPEO+ had a mean height SD of −0.28 (95% CI, −0.59 to 0.02; n = 50).

BMI

In total, 82.1% (n = 472) of patients had a recorded weight, which was used to calculate BMI SD. The mean weight in men (75.4 kg; 95% CI, 72.6 to 78.3; n = 182) was higher than in women (63.6 kg; 95% CI, 61.7 to 65.5; n = 290). There was no significant difference in BMI SD between men (mean = 0.48, SD 1.84, n = 182) and women [mean = 0.54, SD 1.62, n = 290; t(470) = −0.37; P = 0.713; 95% CI, −0.353 to 0.288] (Table 1).
A mean BMI SD of 0.52 (95% CI, 0.37 to 0.67; n = 472) was above average when compared with UK standards [Fig. 1(b)]. The mean BMI SD for the m.3243A>G group was 0.12 (95% CI, −0.10 to 0.34; n = 212), for the single large-scale deletion group it was 0.68 (95% CI, 0.33 to 1.03; n = 78), and for the nuclear group it was 1.13 (95% CI, 0.80 to 1.46; n = 107). There was an observed difference in mean BMI SD [F(2, 394) = 14.04, P < 0.001] between the genetic subgroups. Multiple comparisons revealed that the m.3243A>G subgroup had a lower mean BMI SD than the single large-scale deletion group (mean difference = −0.55; 95% CI, −1.08 to −0.04; P = 0.032) and the nuclear group (mean difference = −1.01; 95% CI, −1.48 to −0.54; P < 0.001).
Further analysis of single, large-scale mtDNA deletion data showed that patients with KSS had a mean BMI SD of −0.55 (95% CI, −2.02 to 0.92; n = 9), patients with CPEO had a mean BMI SD of 1.02 (95% CI, 0.31 to 1.74; n = 18), and those with CPEO+ a mean BMI SD of 0.78 (95% CI, 0.35 to 0.78; n = 46).

NMDAS

A total of 533 adults had a documented height and NMDAS score. NMDAS scores were negatively correlated with height SDS (Pearson correlation coefficient r = −0.25; 95% CI, −0.33 to −0.17; P < 0.001; n = 533) (Fig. 2). A weaker correlation was observed between the NMDAS score and BMI SD (r = −0.10; 95% CI, −0.19 to −0.01; P = 0.033; n = 439).
Figure 2.
Individual boxplots of height SD against NMDAS score. Boxplots showing the median, first quartile, third quartile, and range for height SD (a) all cohort and (b) m.3243A>G against NMDAS score. All cohort n = 533, and m.3243A>G n = 222.
In the m.3243A>G group the NMDAS score was negatively correlated with the height SD (r = −0.39; 95% CI, −0.50 to −0.27; P< 0.001; n = 222). The NMDAS score was also negatively correlated with BMI SDS in this subgroup (r = −0.32; 95% CI, −0.44 to −0.20; P < 0.001; n = 201). In the single, large-scale mtDNA deletion group the NMDAS score was negatively correlated with height SDS (r = −0.30; 95% CI, −0.47 to −0.12; P = 0.006; n = 82); however, there was no association between BMI and NMDAS score (r = −0.18; 95% CI, −0.37 to 0.043; P = 0.130; n = 72). In the nuclear deletion group the NMDAS score was negatively correlated with height SDS (r = −0.193; 95% CI, −0.34 to −0.03; P = 0.027; n = 131). There was no correlation between BMI and NMDAS score (r = 0.12; 95% CI, −0.06 to 0.30; P = 0.262; n = 96).
Mean NMDAS score of those <25 years of age when height was recorded was 23 (95% CI, 16.3 to 29.9; n = 46), and the NMDAS score of those >25 years of age when height was recorded was 26 (95% CI, 24.6 to 28.1; n = 453). NMDAS score was negatively correlated with the height SD of patients <25 years of age when height was recorded (r = −0.43; 95% CI, −0.63 to −0.18; P = 0.003; n = 46) and >25 years when height was recorded (r = −0.22; 95% CI, −0.31 to −0.14; P < 0.001; n = 453).

Heteroplasmy

Age-corrected blood mutant heteroplasmy was explored in our largest subgroup, m.3243A>G (13). Data were available for 92.3% (n = 216) of patients with m.3243A>G. NMDAS score was positively correlated with age-corrected blood heteroplasmy level (r = 0.36; 95% CI, 0.28 to 0.47, P < 0.001). The age-corrected blood heteroplasmy level was negatively correlated with height (r = −0.31; 95% CI, −0.42 to −0.18; P < 0.001) and BMI (r = −0.39; 95% CI, −0.50 to −0.26, P < 0.001) SDs.

Patient-reported disease onset

Further information with regard to patient-reported symptom onset was explored in the two largest subgroups with the m.3243A>G mutation and single large-scale mtDNA deletions. In the m.3243A>G group, onset data were available for 86.8% (n = 166) of patients, of whom 18.2% (n = 37) were asymptomatic carriers. The mean age of onset of patients in the m.3243A>G group was 28.5 years (95% CI, 26.5 to 30.6; n = 166). Age of onset was available for 96.7% (n = 87) of patients with single large-scale deletions, and the mean age of onset was 23.7 years (95% CI, 21.07 to 26.47; n = 87).
In the m.3243A>G mutation group a positive correlation was identified between age of onset and height SD (r = 0.28; 95% CI, 0.13 to 0.41; P ≤ 0.001; n = 166), indicating that later onset was associated with taller final adult height. A positive correlation was also identified between age of onset and BMI SD (r = 0.27; 95% CI, 0.12 to 0.41; P ≤ 0.001; n = 155). No correlation was identified between NMDAS score and age of onset (r = −0.10; 95% CI, −0.25 to 0.06; P = 0.220; n = 159).
In the single, large-scale mtDNA deletion group no relationship was observed between age of onset and final adult height SD (r = 0.010; 95% CI, −0.20 to 0.22; P = 0.912; n = 87), BMI SD (r = 0.05; 95% CI, −0.18. to 0.27; P = 0.678; n = 77), or NMDAS score (r = −0.17; 95% CI, −0.37 to 0.07; P = 0.136; n = 80). In patients with KSS, the mean age of onset was earlier, at 11.1 (±6.2) years of age, compared with other genetic groups.

Disease manifestations and progression

There were 400 patients with multiple NMDAS measurements: 174 in the m.3243A>G subgroup, 90 in the nuclear group, 71 with single deletions, and 65 in the other subgroup. The median number of NMDAS assessments was 4 (interquartile range = 4.00; range = 2 to 22), with a median interval between NMDAS assessments of 5.20 years (interquartile range = 5.85; range = 0, 13.3). Longitudinal linear mixed modeling, used to examine the role of height and BMI in the rate of NMDAS development, revealed a significant negative association between height SD and increased disease progression in the whole cohort (P ≤ 0.001). This association was driven mainly by the m.3243A>G subgroup (P ≤ 0.001). The nuclear gene defects and single large-scale mtDNA deletion groups showed no significant association (P = 0.197 and P= 0.144, respectively).
A similar picture was seen when BMI SD was considered. The significant association with increased disease progression in the whole cohort (P ≤ 0.001) was driven mainly by the m.3243A>G subgroup (P ≤ 0.001), although the single large-scale mtDNA deletion subgroup also showed a significant effect (P = 0.021). The nuclear subgroups showed no association (P = 0.188).
Height SD was negatively correlated with several specific NMDAS fields, namely cerebellar ataxia, cognition, pyramidal involvement, hearing impairment, gastrointestinal disturbance, diabetes, and cardiovascular involvement (range of r = −0.17 to −0.28, all P < 0.001; Table 2). Interestingly, although hearing impairment and gastrointestinal disturbance are also associated with BMI SD (r for both = −0.17, P < 0.001), BMI SD is independently associated with four other phenotypes not associated with height SD (dysphonia or dysarthria, seizures, encephalopathy, and strokelike episodes; range of r is −0.16 to 0.21, P < 0.001).
Table 2.
Correlation of Height and BMI SDs With Disease Manifestations
NMDAS Question Height SD
 
BMI SD
 
r 95% CI
 
P r 95% CI
 
P 
Psychiatric involvement −0.10 −0.19 to −0.02 0.016 545 −0.06 −0.15 to 0.03 0.193 447 
Cerebellar ataxia −0.20 −0.28 to −0.11 <0.001a 536 −0.06 −0.15 to 0.04 0.238 442 
Migraine −0.10 −0.18 to −0.01 0.026 545 −0.10 −0.19 to 0.00 0.039 447 
Cognition −0.28 −0.36 to −0.20 <0.001a 515 −0.13 −0.22 to −0.04 0.007 426 
Neuropathy −0.01 −0.10 to 0.07 0.785 538 0.07 −0.02 to 0.16 0.128 443 
Dysphonia or dysarthria −0.13 −0.22 to −0.05 0.002 539 −0.16 −0.25 to −0.07 0.001a 442 
Seizures −0.08 −0.16 to 0.00 0.065 545 −0.21 −0.30 to −0.12 <0.001a 447 
Encephalopathy −0.09 −0.18 to −0.01 0.031 545 −0.22 −0.31 to −0.13 <0.001a 447 
Strokelike episodes −0.10 −0.18 to −0.01 0.025 545 −0.22 −0.31 to −0.13 <0.001a 447 
Pyramidal −0.19 −0.27 to −0.11 <0.001a 537 −0.11 −0.20 to −0.02 0.017 443 
Extrapyramidal −0.13 −0.21 to −0.05 0.002 537 −0.02 −0.12 to 0.07 0.604 443 
Visual acuity −0.12 −0.21 to −0.04 0.005 529 0.09 0.00 to 0.19 0.049 435 
Ptosis −0.01 −0.10 to 0.07 0.747 539 0.08 −0.01 to 0.18 0.075 443 
CPEO 0.02 −0.07 to 0.10 0.701 540 0.05 −0.04 to 0.14 0.288 443 
Hearing impairment −0.24 −0.32 to −0.16 <0.001a 545 −0.17 −0.26 to −0.08 <0.001a 447 
Gastrointestinal disturbance −0.17 −0.25 to −0.08 <0.001a 545 −0.17 −0.26 to −0.08 <0.001a 447 
Myopathy −0.07 −0.15 to 0.02 0.119 540 −0.10 −0.19 to −0.01 0.031 443 
Diabetes −0.19 −0.27 to −0.10 <0.001a 545 0.02 −0.07 to 0.11 0.649 447 
Cardiovascular involvement −0.21 −0.29 to −0.12 <0.001a 495 −0.07 −0.16 to 0.03 0.170 414 
Pearson correlations of height SD and BMI SD with comorbidities as defined by individual NMDAS questions. Estimated Pearson correlation coefficients and uncorrected P values are shown. The Bonferroni corrected threshold of significance for correlations is 0.001 because 38 independent tests were performed. Disease manifestations represent answers to questions from sections II and III of the NMDAS assessment, because these are physician assessed and therefore more likely to be objective. Hearing impairment from section I was included, because this is the only assessment of this common mitochondrial disease manifestation within the NMDAS.
a
Significant results.

Discussion

Our findings in a large cohort of well-characterized patients supports earlier work indicating that patients with mitochondrial disease are shorter than their peers. In addition, we describe associations between disease severity as measured by the NMDAS and rate of disease progression with reduced final adult height. We have also noted that disease severity is associated with nutritional status, as determined by BMI, in patients with the m.3243A>G mutation, although the overall BMI SD of this patient group is well within normal limits.
It was interesting to note that patients with short stature were more likely to have diabetes and cardiovascular involvement, whereas patients with a low BMI had more severe symptoms such as seizures, encephalopathy, and strokelike episodes, as noted previously (16). Both short stature and low BMI were associated with hearing impairment and gastrointestinal disturbances.
Because adults often attain their final height before the onset of overt disease, our findings suggest that otherwise subclinical impairments in mitochondrial function can affect growth and specifically the behavior of the chondrocyte.
There are a number of possible explanations for the association between phenotype and stature. Muscle function is known to affect skeletal integrity and bone mass, but our data may reflect an impact of mitochondrial and associated muscle dysfunction on bone growth. With this in mind, it is of interest that patients with the milder Becker muscular dystrophy phenotype are taller than those with Duchenne muscular dystrophy, and in children with cerebral palsy, more profound motor dysfunction is associated with slower growth (1718).
Increasing evidence in animal models shows the importance of muscle contractions in longitudinal bone growth; in one study mice with one limb paralyzed had a significantly shorter limb length in the affected limb, demonstrating that muscle loading has a role in stimulating chondrocytes in the growth plate (19).
Although it is possible that reduced muscle mass and exercise intolerance during childhood in patients with early onset mitochondrial disease could affect the growth plate, a more direct effect of abnormal mitochondrial function on chondrocyte function can also be postulated. Hence, dysfunctional ATP synthesis within the chondrocyte itself could be a mechanism leading to poor growth and short stature. In the growth plate, hypertrophy and proliferation of chondrocytes are closely regulated, and it is well recognized that in these cells decreased ATP results in apoptosis (19). In patients with mitochondrial disease this effect may result in increased apoptosis and cell death in the growth plate, with associated suboptimal growth. An increase in apoptosis has been observed in the muscle samples of patients with primary mitochondrial disorders, and a similar process may affect the chondrocyte as a potential mechanism (20).
Another potential explanation for poor growth in these patients is that abnormal mitochondrial function affects GH production and secretion. There have been numerous case reports of GH deficiency in patients with mitochondrial disease in the nuclear, m.3243A>G, and single large-scale mtDNA deletions (21–24). Postulated mechanisms include the presence of mutated mtDNA in the hypophysis, leading to pituitary dysfunction (2125).
Primary lesions affecting the hypothalamic GH releasing factor cells have been reported in two patients with mitochondrial encephalomyopathy lactic acidosis and strokelike episode, with increased GH production noted after the administration of GH releasing factor (23). Treatment with recombinant GH improved stature in a small number of patients with mitochondrial disease, but there are no larger-scale studies (232627). Concerns have been raised about possible long-term detrimental effects of increased GH concentrations in mitochondrial disease, with acute deteriorations reported in patients with KSS and other mitochondrial rearrangements due to increased tissue energy demands after administration (222627).
Suboptimal nutrition and an inability to meet metabolic energy demands could also influence stature and body habitus, although overall BMI SD was above average. Gastrointestinal symptoms commonly reported by patients with mitochondrial disease include reduced oral intake due to poor appetite and early satiety. Dysphagia of varying severity is also well recognized, affecting 48% of patients and is very common in the m.3243A>G subgroup (928). The mechanism is believed to be dysfunctional involvement of smooth muscle, peripheral nerves, and central nervous system in swallowing. Gastroparesis, delayed gastric emptying, bacterial overgrowth, and intestinal pseudo-obstruction have also been noted in this population (16). Malabsorption is commonly experienced by patients, and mtDNA defects have been shown to impair human colonic epithelial function (29). Together, these mechanisms result in malabsorption and reduced body mass. Our data also demonstrate that there is an association between disease severity and BMI. In the m.3243A>G subgroup, disease severity as measured by NMDAS is associated with a fall in BMI. We also noted that patients with KSS, known for their early onset of disease manifestations, are shorter, with a lower BMI than those with CPEO, whose onset is often later in life (30). Body mass and muscle strength affect mechanical loading, which is essential for bone remodeling and maintaining mineral density. Those with a higher BMI have a greater bone mineral density because of mechanical loading (931).
As well as gastrointestinal symptoms, patients with mitochondrial disease may have other comorbidities that could theoretically affect growth and bone health (93233). These include endocrinopathies such as hypoparathyroidism, hypothyroidism, and diabetes. However, these endocrine disorders are variably expressed within patients who have the same genotype. Indeed, phenotypic heterogeneity and tissue-specific vulnerability are often observed across both mtDNA and nuclear gene defects. In this study, our finding that patients with mitochondrial disease caused by the m.3243A>G mutation are generally shorter than patients with other genetic mutations is entirely consistent with this phenotypic heterogeneity, but the exact mechanism of such tissue-specific changes remains elusive.
There is a growing understanding that the systemic effects of stress and inflammation in chronic disease states may affect growth (34). In adults with mitochondrial disease, hypothyroidism is uncommon, and mitochondrial diabetes usually develops after the age of 35 years and hence beyond the phase of childhood and adolescence when growth occurs (3536). Hypoparathyroidism is seen in severe childhood-onset KSS and mitochondrial encephalomyopathy lactic acidosis and strokelike episode, but in these patients short stature is caused by multisystem involvement. Patients with severe mitochondrial disease are more likely to be receiving multiple medications, which may also affect growth; however, because the majority of our patients attained final height before overt disease onset, this potential mechanism is unlikely. The prevalence of one risk factor affecting bone health in mitochondrial disease was estimated at 73% in one cohort, with 30% having four risk factors (33). Patients with mitochondrial disease are at elevated risk of osteoporosis, with more than half of the adult population known to have at least one risk factor (31).
Multiple comorbidities and associated stress could result in elevated circulating levels of growth differentiation factor 15 (37). growth differentiation factor 15 has been cited as a possible biomarker of mitochondrial disease (3839). It inhibits the production of IGF-1 in the liver, limiting body size in response to decreased energy availability. This pathophysiological process could affect growth during childhood and adolescence in the context of patients subsequently diagnosed with a mitochondrial disorder.

Limitations

Missing data in height, weight, and NMDAS scores in some patients may result in selection bias. However, there was no significant difference between the mean NMDAS score of our cohort (n = 537) and those excluded for having no documented height (n = 36) (27 vs 33, P = 0.096).
We wondered whether the negative association between NMDAS score and stature might reflect the impact of the disease process on progress through puberty with a recorded height that was not representative of final adult height in younger patients. An analysis of the association between NMDAS score and height in older patients was similar to the data as a whole, confirming that the association was not simply caused by delayed pubertal development and growth. At present, limited information is available about growth in patients with mitochondrial disease during childhood and adolescence.
Stature is multifactorial, and in the pediatric population midparental height is often used by clinicians to determine a child’s potential final adult height and to look for a possible growth disorder. In our cohort no parental final height data were available; it would be useful to consider whether stature in those with a mitochondrial defect differed from that of unaffected parents.
Intrauterine growth is known to affect final height. In our cohort birthweight and gestation data were not available. Studies looking at birth weight in mitochondrial disease have focused on patients with pediatric mitochondrial disease and have shown that they are significantly lighter than their counterparts; in one study ≤22% were below the third percentile for weight at birth (4041). We know that infants who are small for gestational age often attain a shorter final height than their larger counterparts, particularly if no catch-up growth is demonstrated within the first years of life. Our pregnancy study has recently shown that lower birth weight was evident in the m.3243A>G group compared with other genotypes. It is therefore possible that short stature in patients with more severe mitochondrial disease is linked to smaller size at birth, and this is an area that requires further study.
Because our data are predominantly single-point data, the impact of loss of vertebral height due to osteoporosis and aging is important to consider. We addressed this limitation by correlating age at height measurement and height SDs and found no significant correlation (r = −0.01; 95% CI, −0.10 to 0.07; P = 0.80; n = 537). Nevertheless, it is possible that current height is less than the maximally attained height in a small number of older patients who have osteoporotic fracture, causing loss of vertebral height, or simply the normal age-related change in adult height.

Conclusions

This is a comprehensive description of stature in adults with mitochondrial disease. Numerous mechanisms for short stature can be postulated, including the impact of muscle function on the growth plate, abnormal chondrocyte activity and suboptimal nutrition. To date the definitive mechanism or mechanisms remain unknown.
We suspect that along with the impact of chronic disease, abnormalities of growth in mitochondrial disease also reflect abnormal chondrocyte behavior within the growth plate. More animal mechanistic studies are needed to confirm this suspicion. The role of heteroplasmy levels in chondrocytes remains to be explored.
Our findings suggest that height is a simple measure that may have a predictive value for the disease burden in adult mitochondrial disease. Additional observational analysis of serial measurements of NMDAS scores, height, and weight would help explore these associations in greater detail. Accurate growth data may help to guide clinicians and families toward the more appropriate use of targeted physiotherapy programs in the future. Short stature may be an early indicator of mitochondrial disease and should be considered within the differential diagnosis to optimize patient care.

Abbreviations:

     
  • BMI
    body mass index
  •  
  • CPEO
    chronic progressive external ophthalmoplegia and ptosis
  •  
  • CPEO+
    chronic progressive external ophthalmoplegia and ptosis plus
  •  
  • KSS
    Kearns-Sayre syndrome
  •  
  • Mito Cohort
    UK Mitochondrial Disease Patient Cohort
  •  
  • mtDNA
    mitochondrial DNA
  •  
  • NMDAS
    Newcastle Mitochondrial Disease Adult Scale
  •  
  • SDS
    standard deviation score

Acknowledgments

We are extremely grateful to members of the Mitochondrial Disease Patient Cohort and staff of the National Health Service (NHS) Highly Specialised Mitochondrial Disease Service laboratory in Newcastle. We thank Douglas Jerry for the data management of the Mito Cohort.
Financial Support: This work was supported by the Wellcome Centre for Mitochondrial Research (203105/Z/16/Z), Newcastle University Centre for Ageing and Vitality (supported by the Biotechnology and Biological Sciences Research Council and Medical Research Council L016354), UK National Institute for Health Research (NIHR) Biomedical Research Centre for Ageing and Age-Related Disease award to the Newcastle Upon Tyne Hospitals NHS Foundation Trust, NIHR, and the UK NHS Specialist Commissioners, which funds the Rare Mitochondrial Disorders of Adults and Children Diagnostic Service in Newcastle Upon Tyne (http://www.newcastle-mitochondria.com/). This work also received infrastructure support from the UK Mitochondrial Disease Patient Cohort: A Natural History Study and Patient Registry (REC ref 13/NE/0326), NIHR Biomedical Research Centre, Newcastle and North Tyneside Comprehensive Local Research Network. Y.S.N. holds an NIHR Clinical Lectureship in Neurology (CL-2016-01-003) and was funded by the Medical Research Council Centre of Neuromuscular Diseases for his PhD study (MR/K000608/1). S.P. is funded by Wellcome Trust Fellowship (204709/Z/16/Z). The views expressed are those of the authors and not necessarily those of the NHS, the NIHR, or the Department of Health.
Disclosure Summary: The authors have nothing to disclose.

References

1.
Gorman
 
GS
Schaefer
 
AM
Ng
 
Y
Gomez
 
N
Blakely
 
EL
Alston
 
CL
Feeney
 
C
Horvath
 
R
Yu-Wai-Man
 
P
Chinnery
 
PF
Taylor
 
RW
Turnbull
 
DM
McFarland
 
R
Prevalence of nuclear and mitochondrial DNA mutations related to adult mitochondrial disease
Ann Neurol
 . 
2015
;
77
(
5
):
753
759
.
2.
Gorman
 
GS
Chinnery
 
PF
DiMauro
 
S
Hirano
 
M
Koga
 
Y
McFarland
R
Suomalainen
 
A
Thorburn
 
DR
Zeviani
 
M
Turnbull
 
DM
Mitochondrial diseases
Nat Rev Dis Primers
 . 
2016
;
2
:
16080
.
3.
Calvo
 
SE
Clauser
 
KR
Mootha
 
VK
MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins
Nucleic Acids Res
 . 
2016
;
44
(
D1
):
D1251
D1257
.
4.
Schaefer
 
AM
Phoenix
 
C
Elson
 
JL
McFarland
 
R
Chinnery
 
PF
Turnbull
 
DM
Mitochondrial disease in adults: a scale to monitor progression and treatment
Neurology
 . 
2006
;
66
(
12
):
1932
1934
.
5.
Mackie
 
EJ
Tatarczuch
 
L
Mirams
 
M
The skeleton: a multi-functional complex organ: the growth plate chondrocyte and endochondral ossification
J Endocrinol
 . 
2011
;
211
(
2
):
109
121
.
6.
Wolny
 
S
McFarland
 
R
Chinnery
 
P
Cheetham
 
T
Abnormal growth in mitochondrial disease
Acta Paediatr
 . 
2009
;
98
(
3
):
553
554
.
7.
Castro-Gago
 
M
Gómez-Lado
 
C
Pérez-Gay
 
L
Eirís-Puñal
 
J
Abnormal growth in mitochondrial disease
Acta Paediatr
 . 
2010
;
99
(
6
):
796
.
8.
Xia
 
C-Y
Liu
 
Y
Liu
 
H
Zhang
 
Y-C
Ma
 
YN
Qi
 
Y
Clinical and molecular characteristics in 100 Chinese pediatric patients with m.3243A>G mutation in mitochondrial DNA
Chin Med J (Engl)
 . 
2016
;
129
(
16
):
1945
1949
.
9.
de Laat
 
P
Zweers
 
HE
Knuijt
 
S
Smeitink
 
JA
Wanten
 
GJ
Janssen
MC
Dysphagia, malnutrition and gastrointestinal problems in patients with mitochondrial disease caused by the m3243A&gt;G mutation
Neth J Med
 . 
2015
;
73
(
1
):
30
36
.
10.
Yatsuga
 
S
Povalko
 
N
Nishioka
 
J
Katayama
 
K
Kakimoto
 
N
Matsuishi
 
T
Kakuma
 
T
Koga
 
Y
Taro Matsuoka for MELAS Study Group in Japan
MELAS: a nationwide prospective cohort study of 96 patients in Japan
Biochim Biophys Acta
 . 
2012
;
1820
(
5
):
619
624
.
11.
Cole
 
TJ
Freeman
 
JV
Preece
 
MA
British 1990 growth reference centiles for weight, height, body mass index and head circumference fitted by maximum penalized likelihood
Stat Med
 . 
1998
;
17
(
4
):
407
429
.
12.
Schafer
 
A
Manual for Newcastle Mitochondrial Disease Adult Scale.
  
Newcastle Upon Tyne, UK
Wellcome Centre Mitochondrial Research
2002
. Accessed August 2017. http://www.newcastle-mitochondria.com/clinical-professional-home-page/clinical-publications/nmdas/.
13.
Grady
 
JP
Pickett
 
SJ
Ng
 
YS
Alston
 
CL
Blakely
 
EL
Hardy
 
SA
Feeney
 
CL
Bright
 
AA
Schaefer
 
AM
Gorman
 
GS
McNally
 
RJ
Taylor
 
RW
Turnbull
 
DM
McFarland
 
R
mtDNA heteroplasmy level and copy number indicate disease burden in m.3243A>G mitochondrial disease
EMBO Mol Med
 . 
2018
;
10
(
6
):
e8262
.
14.
Pinheiro JBD DS, Sarkar D; R Core Team. nlme: linear and nonlinear mixed effects models. Updated R package version 3.1-127. 2016. Accessed 18 August 2018. http://F-EMAIL01.xnuth.nhs.uk:32224/?dmVyPTEuMDAxJiYwYTQ4MjQ1ZWNjZTU4N2YxNz01QjkyODRDRl84MjE1NV8xNDcxXzEmJjc5MTU1MThhY2Y0NTYzZT0xMjIyJiZ1cmw9Q1JBTiUyRVItcHJvamVjdCUyRW9yZyUyRnBhY2thZ2UlM0RubG1l.
15.
National Institute for Health and Care Excellence
Human growth hormone (somatropin) for the treatment of growth failure in children
2010
https://www.nice.org.uk/guidance/ta188.
16.
Ng
 
YS
Feeney
 
C
Schaefer
 
AM
Holmes
 
CE
Hynd
 
P
Alston
 
CL
Grady
 
JP
Roberts
 
M
Maguire
 
M
Bright
 
A
Taylor
 
RW
Yiannakou
Y
McFarland
 
R
Turnbull
 
DM
Gorman
 
GS
Pseudo-obstruction, stroke, and mitochondrial dysfunction: a lethal combination
Ann Neurol
 . 
2016
;
80
(
5
):
686
692
.
17.
Oftedal
 
S
Davies
 
PSW
Boyd
 
RN
Stevenson
 
RD
Ware
 
RS
Keawutan
 
P
Benfer
 
KA
Bell
 
KL
Longitudinal growth, diet, and physical activity in young children with cerebral palsy
Pediatrics
 . 
2016
;
138
(
4
):
138
.
18.
Matsumoto
 
M
Awano
 
H
Lee
 
T
Takeshima
 
Y
Matsuo
 
M
Iijima
 
K
Patients with Duchenne muscular dystrophy are significantly shorter than those with Becker muscular dystrophy, with the higher incidence of short stature in Dp71 mutated subgroup
Neuromuscul Disord
 . 
2017
;
27
(
11
):
1023
1028
.
19.
Killion
 
CH
Mitchell
 
EH
Duke
 
CG
Serra
 
R
Mechanical loading regulates organization of the actin cytoskeleton and column formation in postnatal growth plate
Mol Biol Cell
 . 
2017
;
28
(
14
):
1862
1870
.
20.
Mirabella
 
M
Di Giovanni
 
S
Silvestri
 
G
Tonali
 
P
Servidei
 
S
Apoptosis in mitochondrial encephalomyopathies with mitochondrial DNA mutations: a potential pathogenic mechanism
Brain
 . 
2000
;
123
(
Pt 1
):
93
104
.
21.
Matsuzaki
 
M
Izumi
 
T
Ebato
 
K
Suzuki
 
H
Shishikura
 
K
Osawa
 
M
Fukuyama
 
Y
Shimizu
 
N
[Hypothalamic GH deficiency and gelastic seizures in a 10-year-old girl with MELAS]
No To Hattatsu
 . 
1991
;
23
(
4
):
411
416
.
22.
Obara-Moszynska
 
M
Maceluch
 
J
Bobkowski
 
W
Baszko
 
A
Jaremba
 
O
Krawczynski
 
MR
Niedziela
 
M
A novel mitochondrial DNA deletion in a patient with Kearns-Sayre syndrome: a late-onset of the fatal cardiac conduction deficit and cardiomyopathy accompanying long-term rGH treatment
BMC Pediatr
 . 
2013
;
13
(
1
):
27
.
23.
Matsuzaki
 
M
Izumi
 
T
Shishikura
 
K
Suzuki
 
H
Hirayama
 
Y
Hypothalamic growth hormone deficiency and supplementary GH therapy in two patients with mitochondrial myopathy, encephalopathy, lactic acidosis and stroke-like episodes
Neuropediatrics
 . 
2002
;
33
(
5
):
271
273
.
24.
Broomfield
 
A
Sweeney
 
MG
Woodward
 
CE
Fratter
 
C
Morris
 
AM
Leonard
 
JV
Abulhoul
 
L
Grunewald
 
S
Clayton
 
PT
Hanna
 
MG
Poulton
 
J
Rahman
 
S
Paediatric single mitochondrial DNA deletion disorders: an overlapping spectrum of disease
J Inherit Metab Dis
 . 
2015
;
38
(
3
):
445
457
.
25.
Shiraiwa
 
N
Ishii
 
A
Iwamoto
 
H
Mizusawa
 
H
Kagawa
 
Y
Ohta
 
S
Content of mutant mitochondrial DNA and organ dysfunction in a patient with a MELAS subgroup of mitochondrial encephalomyopathies
J Neurol Sci
 . 
1993
;
120
(
2
):
174
179
.
26.
Barberi
 
S
Bozzola
 
E
Berardinelli
 
A
Meazza
 
C
Bozzola
 
M
Long-term growth hormone therapy in mitochondrial cytopathy
Horm Res
 . 
2004
;
62
(
2
):
103
106
.
27.
Romano
 
S
Samara
 
D
Crosnier
 
H
Valayannopoulos
 
V
Polak
 
M
Chrétien
 
D
Rötig
 
A
Munnich
 
A
Brauner
 
R
de Lonlay
 
P
Variable outcome of growth hormone administration in respiratory chain deficiency
Mol Genet Metab
 . 
2008
;
93
(
2
):
195
199
.
28.
Read
 
JL
Whittaker
 
RG
Miller
 
N
Clark
 
S
Taylor
 
R
McFarland
 
R
Turnbull
 
D
Prevalence and severity of voice and swallowing difficulties in mitochondrial disease
Int J Lang Commun Disord
 . 
2012
;
47
(
1
):
106
111
.
29.
Nooteboom
 
M
Johnson
 
R
Taylor
 
RW
Wright
 
NA
Lightowlers
RN
Kirkwood
 
TB
Mathers
 
JC
Turnbull
 
DM
Greaves
 
LC
Age-associated mitochondrial DNA mutations lead to small but significant changes in cell proliferation and apoptosis in human colonic crypts
Aging Cell
 . 
2010
;
9
(
1
):
96
99
.
30.
Alston
 
CL
Schaefer
 
AM
Raman
 
P
Solaroli
 
N
Krishnan
 
KJ
Blakely
 
EL
He
 
L
Craig
 
K
Roberts
 
M
Vyas
 
A
Nixon
 
J
Horvath
 
R
Turnbull
 
DM
Karlsson
 
A
Gorman
 
GS
Taylor
 
RW
Late-onset respiratory failure due to TK2 mutations causing multiple mtDNA deletions
Neurology
 . 
2013
;
81
(
23
):
2051
2053
.
31.
Evans
 
AL
Paggiosi
 
MA
Eastell
 
R
Walsh
 
JS
Bone density, microstructure and strength in obese and normal weight men and women in younger and older adulthood
J Bone Miner Res
 . 
2015
;
30
(
5
):
920
928
.
32.
Chow
 
J
Rahman
 
J
Achermann
 
JC
Dattani
 
MT
Rahman
 
S
Mitochondrial disease and endocrine dysfunction
Nat Rev Endocrinol
 . 
2017
;
13
(
2
):
92
104
.
33.
Gandhi
 
SS
Muraresku
 
C
McCormick
 
EM
Falk
 
MJ
McCormack
SE
Risk factors for poor bone health in primary mitochondrial disease
J Inherit Metab Dis
 . 
2017
;
40
(
5
):
673
683
.
34.
Wong
 
SC
Dobie
 
R
Altowati
 
MA
Werther
 
GA
Farquharson
 
C
Ahmed
 
SF
Growth and the growth hormone-insulin like growth factor 1 axis in children with chronic inflammation: current evidence, gaps in knowledge, and future directions
Endocr Rev
 . 
2016
;
37
(
1
):
62
110
.
35.
Schaefer
 
AM
Walker
 
M
Turnbull
 
DM
Taylor
 
RW
Endocrine disorders in mitochondrial disease
Mol Cell Endocrinol
 . 
2013
;
379
(
1–2
):
2
11
.
36.
Whittaker
 
RG
Schaefer
 
AM
McFarland
 
R
Taylor
 
RW
Walker
 
M
Turnbull
 
DM
Prevalence and progression of diabetes in mitochondrial disease
Diabetologia
 . 
2007
;
50
(
10
):
2085
2089
.
37.
Wang
 
T
Liu
 
J
McDonald
 
C
Lupino
 
K
Zhai
 
X
Wilkins
 
BJ
Hakonarson
 
H
Pei
 
L
GDF15 is a heart-derived hormone that regulates body growth
EMBO Mol Med
 . 
2017
;
9
(
8
):
1150
1164
.
38.
Yatsuga
 
S
Fujita
 
Y
Ishii
 
A
Fukumoto
 
Y
Arahata
 
H
Kakuma
 
T
Kojima
 
T
Ito
 
M
Tanaka
 
M
Saiki
 
R
Koga
 
Y
Growth differentiation factor 15 as a useful biomarker for mitochondrial disorders
Ann Neurol
 . 
2015
;
78
(
5
):
814
823
.
39.
Davis
 
RL
Liang
 
C
Sue
 
CM
A comparison of current serum biomarkers as diagnostic indicators of mitochondrial diseases
Neurology
 . 
2016
;
86
(
21
):
2010
2015
.
40.
von Kleist-Retzow
 
J-C
Cormier-Daire
 
V
Viot
 
G
Goldenberg
 
A
Mardach
 
B
Amiel
 
J
Saada
 
P
Dumez
 
Y
Brunelle
 
F
Saudubray
JM
Chrétien
 
D
Rötig
 
A
Rustin
 
P
Munnich
 
A
De Lonlay
 
P
Antenatal manifestations of mitochondrial respiratory chain deficiency
J Pediatr
 . 
2003
;
143
(
2
):
208
212
.
41.
Tavares
 
MV
Santos
 
MJ
Domingues
 
AP
Pratas
 
J
Mendes
 
C
Simões
 
M
Moura
 
P
Diogo
 
L
Grazina
 
M
Antenatal manifestations of mitochondrial disorders
J Inherit Metab Dis
 . 
2013
;
36
(
5
):
805
811
.
This article has been published under the terms of the Creative Commons Attribution License (CC BY; https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Copyright for this article is retained by the author(s).