Department of Endocrinology, Ordos Center Hospital, Ordos, Inner Mongolian, 017010, People’s Republic of China

Correspondence: Xing Li, Department of Endocrinology, Ordos Center Hospital, Sudu Street No. 6, Kangbashi District, Ordos, Inner Mongolian, 017010, People’s Republic of China, Tel/Fax +86-0477-8167761, Email [email protected]

Background: Mutations in mitochondrial tRNA (mt-tRNA) could be the origin of some type 2 diabetes mellitus (T2DM) cases, but the mechanism remained largely unknown.
Aim: The aim of this study was to assess the impact of a novel mitochondrial tRNACys/tRNATyr A5826G mutation on the development and progression of T2DM.
Methods: A four-generation Han Chinese family with maternally inherited diabetes underwent clinical, genetic and biochemical analyses. The mitochondrial DNA (mtDNA) mutations of three matrilineal relatives were screened by PCR-Sanger sequencing. Furthermore, to see whether m.A5826G mutations affected mitochondrial functions, the cybrid cell lines were derived from three subjects with m.A5826G mutation and three controls without this mutation. ATP was evaluated by luminescent cell viability assay, mitochondrial membrane potential (MMP), and reactive oxygen species (ROS) were determined by flow cytometry. The student’s two-tailed, unpaired t-test was used to assess the statistical significance between the control and mutant results.
Results: The age at onset of diabetes in this pedigree varied from 40 to 63 years, with an average of 54 years. Mutational analysis of mitochondrial genomes revealed the presence of a novel m.A5826G mutation. Interestingly, the m.A5826G mutation occurred at the conjunction between tRNACys and tRNATyr, a very conserved position that was critical for tRNAs processing and functions. Using trans-mitochondrial cybrid cells, we found that mutant cells carrying the m.A5826G showed approximately 36.5% and 22.4% reductions in ATP and MMP, respectively. By contrast, mitochondrial ROS levels increased approximately 33.3%, as compared with the wild type cells.
Conclusion: A novel m.A5826G mutation was identified in a pedigree with T2DM, and this mutation would lead to mitochondrial dysfunction. Thus, the genetic spectrum of mitochondrial diabetes was expanded by including m.A5826G mutation in tRNACys/tRNATyr, our study provided novel insight into the molecular pathogenesis, early diagnosis, prevention and clinical treatment for mitochondrial diabetes.

Introduction

Mitochondria are distinct organelles within eukaryotic cells that are involved in oxidative phosphorylation (OXPHOS). They carry their own genetic material, in humans, mitochondrial DNA (mtDNA) is a circular molecule comprising 16,569-bp that is essential for respiratory chain functions.1 It encodes 13 OXPHOS-related proteins, as well as 22 tRNAs and 2 rRNAs. MtDNA in humans shows a strictly maternal inheritance as during fertilization of the oocyte the paternal mitochondria are not retained within the fertilized oocyte.2 Mitochondrial diseases are a group of genetic disorders characterized by defects in OXPHOS and caused by mutations in nuclear DNA (nDNA) and mtDNA that encode proteins involved in mitochondrial functions.3 Mitochondrial diseases are frequently related to the damage of nerves and muscle conducting to deafness and visual loss. Among them, T2DM is one of the most common types of mitochondrial diseases. It is a complex, heritable and heterogeneous condition that manifests in a range of clinical phenotypes, including hyperglycemia, diabetic ketoacidosis, retinopathy, nephropathy, foot problem which remain big challenging for clinicians.4 However, its pathophysiology is still undetermined.

In recent years, it had been found that mutations in mtDNA may lead to T2DM concomitant with symptoms caused by damage to nerves and muscles.5–8 In particular, the well-known m.A3243G mutation in tRNALeu(UUR) was implicated in the pathogenesis and progression of T2DM in many families worldwide.9–11 This mutation affected the structural stability, aminoacylation rate and post-transcriptional modification of tRNALeu(UUR).12 In addition to the m.A3243G mutation, some case–control studies suggested that tRNATrp A5514G and tRNASer(AGY) C12337T,13 tRNAGlu A14692G,14 tRNAAla T5587C and ND5 T12338C15 mutations were associated with T2DM. In fact, mt-tRNA pathogenic mutations have structural and functional consequences, such as affecting the tRNA structure, altering 5’ or 3’ processing of tRNAs, CCA addition and leading to defects in chemical modifications. Thus, these mutations would impair the normal functions of OXPHOS system,16 thereby exacerbating the mitochondrial dysfunction that is responsible for T2DM. However, the detailed mechanisms remain largely unknown.

To address this problem, most recently, we carried out a mutational screening of mt-tRNA genes from Department of Endocrinology, Ordos Center Hospital, China. Our arm was to describe a Han Chinese family with maternally transmitted T2DM, as well as carrying out a sequence analysis of the complete mitochondrial genomes from matrilineal relatives. As a result, we identified the presence of two potential pathogenic mutations: m.A5826G which resided at a spacer between tRNACys and tRNATyr gene, together with m.C14668T in ND6 gene. Therefore, our study broadened the mutational spectrum of mitochondrial diabetes. Moreover, we evaluated the impacts of m.A5826G and m.C14668T mutations on mitochondrial functions by creating cybrid cell lines derived from the three patients with these mtDNA mutations, three healthy controls without these mutations. We also examined the levels of ATP, MMP and ROS. Consequently, we found that cells with m.A5826G and m.C14668T mutations exhibited more severe mitochondrial dysfunctions than control cell lines.

Materials and MethodsParticipants and Clinical Examinations

In this case–control study for genetic screening of T2DM-associated novel mtDNA mutations, a four-generation Han Chinese pedigree with diabetes, as shown in Figure 1, was ascertained via the Department of Endocrinology, Ordos Center Hospital. This pedigree was designed according to the proposed system by the National Society of Genetic Counselors (NSGC), arrow indicated the proband, and affected individuals were indicated by filled symbols.17 Members of this family were evaluated at length to identify both personal and family medical histories of diabetes, deafness, vision loss and other clinical abnormalities. Moreover, 200 healthy Han Chinese subjects, including 100 males and 100 females, aged from 33 to 50 years, with an average of 41 years were recruited from the Physical Examination Center of our hospital as controls. The inclusion criteria for these control subjects were as follows: healthy subjects without any diseases or did not have any family history of mitochondrial diseases; subjects who have any history of mitochondrial diseases were excluded. The protocol used in this investigation was in accordance with the principles expressed in 1975 Declaration of Helsinki, which was revised in 2008, and the Ethic Committee of Ordos Center Hospital approved this study (Approval No. 2023–036). The written informed consent for participating in this study as well as informed consent to have their case details published were obtained from all individuals enrolled in this study, and all information was kept in confidential in order to protect their privacy.

Figure 1 T2DM pedigree information carrying m.A5826G and m.C14668T mutations. Affected individuals (I-2; II-4, II-6, III-4 and III-8) are indicated by filled symbols. Diagonal lines (I-1, I-2 and III-8) are used to show that persons are deceased, arrow shows the proband (III-4).

The diagnosis of T2DM was based on the standard proposed by the American Diabetes Association,18 there were four possible ways to diagnose diabetes and each, in the absence of hyperglycemia symptoms, must be confirmed on a different day by any one of the other three: (1) The fasting plasma glucose concentration ≥7.0 mmol/L; (2) The 2-h plasma glucose level ≥11.1 mmol/L post the oral glucose tolerance test (OGTT). It should be noted that OGTT has poor reproducibility in adolescents, with a concordance rate of <30% between tests performed a few weeks apart.19 (3) Random plasma glucose ≥11.1 mmol/L in the presence of diabetes symptoms. (4) The concentration of hemoglobin A1c (HbA1c) >6.5%. Notably, the fasting plasma glucose was performed 0.5-h prior to vein punction and was measured by chemiluminescent immunometric assay (IMMULITE®, Siemens).

All affected subjects underwent a comprehensive physical examinations and laboratory assessments: the body mass index (BMI) was calculated as the body weight (kg) divided by the square of the height (m2). BMI ≤18.5 kg/m2 suggested the person was underweight, 18.5 kg/m2<BMI<23.9 kg/m2 indicated normal weight, 24.0 kg/m2<BMI<27.9 kg/m2 indicated overweight, BMI≥ 28.0 kg/m2 indicated obesity.20

Blood samples were collected in the morning between 7:00 Am and 10:00 Am after an overnight fast. The HbA1c was measured by high-pressure liquid chromatography (Variant II; Bio-Rad, CA, USA). For biochemical assessment, the fasting plasma glucose level was determined by regular laboratory methods (Beckman Coulter, Tokyo, Japan). In addition, OGTT was carried out by measurement of plasma glucose concentrations at 2-h after 75-g glucose administration, while the fasting plasma insulin was measured by chemiluminescent immunometric assay (IMMULITE®, Siemens).21 The degree of insulin resistance (IR) was estimated using the homeostasis model assessment of IR (HOMA-IR) index, the HOMA-IR = (insulin × glucose)/22.5. A HOMA-IR≥2.69 was regarded to be IR.22 In addition, the serum creatinine (Cr) and estimated glomerular filtration rate (eGFR) were assessed to evaluate the renal functions. Notably, eGFR = 175 × standardized Scr−1.154 × age−0.203 × 1.212 (if black) × 0.742 (if female), where eGFR was expressed as mL/min/1.73m2 of body surface area and serum Cr was expressed in µmol/L.23

The levels of hearing loss were tested in each affected subject by using pure-tone audiometry (PTA), as described previously.24 The PTA was calculated as the average of the thresholds measured at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz. The severity of hearing impairment was defined as mild (21–40 dB), moderate (41–70 dB), severe (71–90 dB), and profound (>91 dB). Moreover, the visual acuity was tested separately for each eye without refraction at 4-m using Early Treatment Diabetic Retinopathy Tumbling-E study charts (Precision Vision, La Salle, IL, USA), and the degree of visual impairment was classified on the basis of visual acuity as follows: normal (>0.3), mild (0.1–0.3), moderate (0.05–0.1), severe (0.02–0.05), and profound (<0.02).25

Analysis of mtDNA Mutations

The genomic DNA of three matrilineal relatives (II-4, II-6 and III-4), together with 200 controls were extracted by using the Puregene DNA Isolation Kits (Gentra Systems). Subsequently, the affected individuals’ DNA fragments spanning the whole mitochondrial genomes were PCR amplified in 24 overlapping fragments using sets of the light and the heavy strand primers as described previously.26 Furthermore, the ABI 3700 DNA instrument was employed to analyze the sequences by comparing with the revised Cambridge sequences (rCRS, GenBank accession number: NC_012920.1).27 The DNA STAR software package version 5.01 (Madison) was used to detect mtDNA mutations or variants.

Phylogenetic Conservation Analysis

A total of 14 species’ mtDNA sequences were used in the interspecific analysis. The conservation index (CI) was calculated by comparing the human nucleotide variants with other 13 vertebrates.28 The CI was defined as the percentage of species from the list of 13 different primate species that have the wild-type nucleotide at that position. Its value ≥75% was regarded as having functional importance.29

Generation of Cybrid Cell Lines

Cybrids were obtained by transforming platelets of patients containing the m.A5826G and m.C14668T mutations (II-4, II-6 and III-8), together with three control subjects without these mtDNA mutations (C1, C2 and C3) into mtDNA-less ρ0 human osteosarcoma 143B cells, as described previously.30 Cells were cultured in high-glucose Dulbecco’s modified Eagle medium (DMEM) (Thermo Fisher Scientific, Waltham, MA, USA) containing 12% FBS (Sigma, Aldrich, St. Louis, MO, USA) at 37°C.

To confirm the successful establishment of cybrid cells, PCR-Sanger sequencing was performed to detect the presence of m.A5826G and m.C14668T mutations in cybrids. The genomic DNA was first isolated from cell lines, and the concentration of DNA was then measured. The primer sequences for detection of m.A5826G mutation were: forward-5’-CTA ACC GGC TTT TTG CCC-3’; reversed-5’-ACC TAG AAG GTT GCC TGG CT-3’; the primer sequence for detection of m.C14468T mutation was: forward-5’-GCA TAA TTA AAC TTT ACT TC-3’; reversed-5’-AGA ATA TTG AGG CGC CAT TG-3’. The PCR primers were supplied by BGI (Shenzhen, China) and the PCR mixture included 200 µM dNTP, 10× buffer, Taq DNA polymerase and 15 mmol/L Mg2+ (Takara Biotechnology, Dalian, China). After amplification, the PCR products were purified and sequenced by ABI 3700 DNA instrument, and the data were then compared with rCRS to detect the presence of m.A5826G and m.C14668T mutations (GenBank accession number: NC_012920.1).27

ATP Measurement

The levels of ATP production in six cell lines were determined by CellTiter-Glo® Luminescent cell viability assay (Promega, G7572) in accordance with the protocol provided by the manufacturer’s instructions.31

MMP Analysis

MMP from six cell lines were examined with JC-10 Assay Kit-Flow Cytometry (Abcam) following general manufacturer’s recommendations with some modifications, as described previously.32 Cells were first plated onto 96-well cell culture plate overnight in growth medium. JC-10 dye loading solution was added for 30 min at 37°C, 5% CO2. Alternatively, plated cells were preincubated with 10 µM of the protonophore uncoupler carbonyl cyanide 3-chlorophenylhydrazone (CCCP) for 30 min at 37°C, 5% CO2 prior to staining with JC-10 dye. The fluorescent intensities for both J-aggregates and monomeric forms of JC-10 were measured at Ex/Em=490/530 and 490/590 nm with flow cytometry.

ROS Analysis

Approximate 2 × 106 cells were grown in 96-well plates, first washed twice with PBS to remove the medium. Next, adding 0.5 mL of measurement buffer containing 2 mM 2’,7’-dichlorodihydrofluorescein diacetate (DCFH-DA). Immediately after addition, flow cytometry was used to measure the fluorescence in kinetic mode at 485 nm excitation and 520 nm emission wavelengths.33

Statistical Analysis

The mitochondrial functional data were expressed as mean ± standard deviation (SD). The student’s two-tailed, unpaired t-test was used to assess the statistical significance between the control and mutant results. All analyses were performed using the SPSS software version 22.0. We regarded the p<0.05 as statistically significant.

ResultsClinical Features

We ascertained a maternally transmitted T2DM via the Department of Endocrinology, Ordos Center Hospital, China (Figure 1). Detailed information was obtained from each subject of this pedigree, as well as the family history of diabetes, deafness, myopia or renal diseases. The proband (III-4) was a 42-year-old woman who came from Ordos City of Inner Mongolia Autonomous Region. The fasting glucose level, as well as the HbA1c and OGTT results, strongly indicated that she was a diabetic carrier. She also exhibited IR and obesity (HOMA-IR = 7.019, BMI = 28 kg/m2). After the genetic counseling, we noticed that she began to develop T2DM when she was 40, she was treated with metformin (0.5 μ/kg). In addition, family history suggested that the subjects (I-2 and III-8) died several years ago due to diabetes and its complications. Meanwhile, the matrilineal relatives (II-4 and II-6) were also diabetic carriers. Notably, both of them (II-4 and II-6) had hearing and vision loss, as well as declining in kidney functions. The detailed clinical features of these affected individuals were listed in Table 1.

Table 1 Clinical and Biochemical Characterizations of Matrilineal Relatives in This Pedigree with DM

MtDNA Mutations Analysis

Owing to the maternally transmission found in this four-generation Han Chinese pedigree, which suggested that mutations in mtDNA genes were the molecular basis for T2DM and its complications, a PCR-Sanger sequencing was performed in matrilineal relatives (II-4, II-6 and III-4), and as a result, a total of 32 variants were identified (Table 2), which belonged to East Asian mitochondrial haplogroup D4.34 There were eight variants in D-loop, two variants in 12S rRNA and two variants in 16S rRNA, one mutation in tRNACys and tRNATyr (m.A5826G), and other variants occurred at OXPHOS-associated genes. Moreover, six missense variants were identified, including A6 A8701G (Ala to Thr) and A8860G (Thr to Ala), ND3 A10398G (Thr to Ala), ND6 C14668T (Tyr to Asn), CytB C14766T (Thr to Ile) and A15326G (Thr to Ala). We further used phylogenetic conservation analysis including mouse,35 bovine36 and Xenopus laevis37 to evaluate the potential pathogenicity. We found that, except for the m.A5826G and m.C14668T mutations (Figures 2 and 3), other mutations were not well conserved and may not be involved in the pathogenesis of T2DM. Furthermore, the m.A5826G and m.C14668T mutations were not detected in 200 healthy subjects, highlighting the contributions of these mutations to T2DM.

Table 2 mtDNA Variants in This Pedigree with T2DM

Figure 2 Sequence analysis of m.A5826G mutation in tRNACys/tRNATyr genes and m.C14668T mutation in ND6 gene.

Figure 3 Evolutionary conservation analysis: sequence alignment of tRNACys gene from 14 vertebrates, arrow indicates the position 1, corresponding to the m.A5826G mutation, suggesting that the m.A5826G mutation is very conserved between different species.

As shown in Figure 4, the novel m.A5826G mutation occurred at a spacer immediately to the 3’ end of tRNACys and 5’ end of tRNATyr in the light chain of mitochondrial genome.38 A position was important for tRNA steady state level and function. Furthermore, the m.C14668T mutation led to an alternation of evolutionary conserved Tyr to Asn that may cause an aberrant in ND6 metabolism. A recent case–control study revealed that m.C14668T mutation may be involved in schizophrenia, bipolar disorder, and major depressive disorder, highlighting the importance of this mutation to mitochondrial dysfunction.39

Figure 4 Secondary structure of tRNACys and tRNATyr, arrows indicate the location of m.A5826G mutation, the secondary structure of mt-tRNA genes are derived from tRNA dB database (http://mttrna.bioinf.uni-leipzig.de/mtDataOutput/).

Reduced in ATP Production

Mitochondrial ATP production is the main energy source for intracellular metabolic pathways.40 To see whether m.A5826G and m.C14668T mutations affected mitochondrial function, we measured ATP levels in three subjects with these mutations and three controls, as shown in Figure 5A, the mutant cell lines (relative level: 63.5 ± 8.2) exhibited approximately 36.5% reductions in ATP production as compared with controls (relative level: 100 ± 11.3) (p < 0.0001).

Figure 5 Analyses of mitochondrial functions in cybrids. (A) ATP analysis; (B) Analysis of MMP; (C) ROS measurement.

Decreased in MMP

The MMP generated by proton pumps (Complexes I, III and IV) is an essential component in the process of energy storage during OXPHOS.41 As can be seen in Figure 5B, mutant cells with m.A5826G and m.C14668T (relative fluorescence level: 78.6 ± 6.2) showed a markedly decreased in MMP level as compared with controls (relative fluorescence level: 100 ± 9.3) (p < 0.0001).

ROS Increased

Mitochondria are the main source of ROS in the cell and, particularly, the enzymes in the electron transport chain are the major contributors to this phenomenon.42 As shown in Figure 5C, we noticed that mutant cells with m.A5826G and m.C14668T (relative level: 133.3 ± 10.7) exhibited much higher levels of ROS when comparing with the controls (relative level: 100 ± 4.3) (p=0.0018).

Discussion

The followings were the major findings of the current study: (1) m.A5826G mutation occurred at the conjunction between tRNACys and tRNATyr, which was critical for tRNA steady-state level and function; (2) m.C14668T mutation was located at a well-conserved position that may affect mitochondrial ND6 metabolism; (3) Cell lines carrying m.A5826G and m.C14668T mutations exhibited much more severe of mitochondrial damage, ATP and MMP dropped significantly, whereas ROS increased; (4) The novel m.A5826G and m.C14468T mutations were involved in the pathogenesis of T2DM; (5) The incomplete penetrance of T2DM suggested that m.A5826G and m.C14468T mutations were not sufficient to produce clinical phenotypes; hence, other risk factors such as environmental factors, epigenetic modification may contribute to T2DM expression.

In recent years, mitochondrial genome mutations had been the target of study in the pathogenesis and progression of different chronic metabolic diseases.43,44 Mitochondria can store between two and ten copies of mtDNA. The lack of histones in mtDNA led to a mutation rate that was six to 17 times higher than that of nDNA. Furthermore, it had been reported that mitochondrial metabolism was involved in the processes that controlled insulin release from pancreatic β-cells,45 in addition, mitochondrial dysfunction was intricately linked to T2DM through several pathways, which included diminished energy production, oxidative stress, the altered metabolism of fatty acids, inflammation, and aberrations in cell signaling.46

To see the contributions of mtDNA mutations to T2DM and its complications, a four-generation Han Chinese family with maternally transmitted T2DM was assessed by clinical and biochemical analyses. Clinically, the age at onset of T2DM varied from 40 to 63 years, with an average of 54 years. Despite the fact that several affected individuals developed a wide range of clinical phenotypes including hearing and vision impairment, as well as declined renal functions, the pattern of T2DM transmission was consistent with maternally inheritance, suggesting that mtDNA mutations were the molecular basis for T2DM in this pedigree.

Sequence analysis of the entire mitochondrial genome from matrilineal relatives revealed a set of polymorphisms, in particular, the m.A5826G and m.C14668T were the two identified as likely involved in the pathogenesis. In fact, the novel homoplasmic m.A5826G mutation was localized at the junction between the tRNACys and tRNATyr at the light strand.38 Moreover, the m.A5826G mutation disrupted the evolutionary conserved A1-U73 Watson–Crick base-pairing. Interestingly, the homoplasmic m.A4401G mutation also occurred at the junction between tRNAMet and tRNAGln decreased the 5’ end processing efficiency of tRNAGln and tRNAMet precursors; furthermore, the m.A4401G reduced the steady-state level of tRNAGln and tRNAMet and was regarded to be a pathogenic mutation for essential hypertension.47 Thus, we proposed that the m.A5826G mutation, which was similar to the m.A4401G mutation, may also affect the tRNACys and tRNATyr 3’ end precursors and lead to the failure in tRNAs metabolism.

Human mtDNA encoded 13 proteins that were essential for OXPHOS. The only protein encoded in the mtDNA light-strand, ND6 played a critical role in the proper assembly Complex I.48 Mutation in the ND6 gene caused severe mitochondrial respiratory dysfunction and several mitochondrial genetic diseases.49 Due to the central role of Complex I in the pathological ROS production, the m.C14668T mutation was believed to be involved in the pathogenesis of T2DM in this pedigree. In fact, the homoplasmic m.C14668T mutation caused the well-conserved Tyr to Asn, was first described to be associated with schizophrenia, bipolar disorder and major depressive disorder,39 and the conserved Tyr to Asn transition perturbed the structural and function of ND6, thereby affecting Complex I deficiency, as in the case of T2DM associated ND6 T14502C mutation.11

Trans-mitochondrial cybrids were created by fusing cells devoid of mtDNA (ρ0 cells) with cytoplasts from different individuals, so the resultant cybrids had uniform nuclear background but different mtDNA.50 With the development and use of the cybrid model, questions related to the importance of the mtDNA variants and mitochondrial–nuclear interactions can be addressed.51 In fact, cybrid cells were already used to model human mitochondrial defects in the pathogenesis of T2DM. T2DM cybrid cells recapitulated diabetes’ pathological features, including excessive oxidative stress, cytochrome oxidase defects, decreased MMP, impaired energy metabolism, and alterations in mitochondrial dynamics.52 However, the greatest limiting factor to generation cybrid cell lines may turn out to be mtDNA heteroplasmy. The traditional Sanger-based sequencing approach was not designed to detect low abundance heteroplasmic mutations; if the mutational load was less than 20–30%, the ability to detect the mtDNA change was unreliable.53

In this study, we created cybrid cells containing homoplasmic m.A5826G and m.C14668T mutations, as well as control subjects without these mtDNA mutations. Approximately 36.5% and 31.8% in ATP and MMP reductions were observed in mutant cells with these mutations, which were below a proposed threshold level to develop a clinical phenotype.54 Because a certain mt-tRNA mutation can be complemented at the translational level by wild-type tRNA genes to ensure a normal rate of mitochondrial protein synthesis. Thus, genetic and functional complementation of mitochondrial defects can occur at the following levels: transcription, translation, enzyme assembly and activity, mitochondrial respiration and ATP synthesis, cell activity and organ function, and all of these different compensatory mechanisms can combine to result in the observed phenotypic threshold effect.55

Indeed, ATP was required for DNA replication, biosynthesis, protein assembly, and biochemical transport. Defects in ATP synthesis revealed an impairment of OXPHOS system.56,57 Furthermore, MMP was critical for maintaining the physiological function of the respiratory chain to generate ATP.58 Thus, loss of MMP would lead to mitochondria-mediated apoptosis.59 As a result, a decrease in ATP and MMP enhanced ROS production that was responsible for the accumulation of ROS-associated damages in DNA, proteins and lipids. Furthermore, mitochondria ROS were closely involved in cellular activity, such as cell differentiation, multiplication and death. This might affect the life and function of pancreatic β-cells.60 When cells misplaced, a large amount of ROS contributed to apoptosis condition of pancreatic β-cells, resulting in limited regeneration ability. Thus, in the presence of increased β-cell loss, either by apoptosis or dedifferentiation, a scenario would arise in which β-cells were lost with insufficient replacement, resulting in inadequate insulin secretion to maintain proper glucose homeostasis,61,62 and contributing to T2DM progression.63

Conclusion

In conclusion, our study suggests that there is a type of mitochondrial T2DM disease related to tRNACys/tRNATyr A5826G and ND6 C14668T mutations. The m.A5826G and m.C14668T should be considered risk polymorphisms for T2DM mitochondrial disease. Furthermore, the incomplete penetrance of T2DM suggested that m.A5826G and m.C14468T mutations were not sufficient to produce clinical phenotypes; hence, other risk factors such as environmental factors or epigenetic modification may contribute to T2DM expression. However, our study only enrolled one pedigree, and further investigation including more patients was needed to verify the conclusions.

Data Sharing Statement

The datasets generated and/or analyzed during the current study are not publicly available due to the original design and validation of the project but are available from the corresponding author on reasonable request.

Ethics Statement

This study was approved by the Ethics Committee of Ordos Center Hospital, and all participants provided their informed consent to participate in this study.

Acknowledgments

We thanked the patients and control subjects for participating in this study.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Funding

There is no funding to report.

Disclosure

The authors declare that they have no conflicts of interest in this work.

References

1. Anderson S, Bankier AT, Barrell BG, et al. Sequence and organization of the human mitochondrial genome. Nature. 1981;290(5806):457–465. doi:10.1038/290457a0

2. Shoubridge EA, Wai T. Mitochondrial DNA and the mammalian oocyte. Curr Top Dev Biol. 2007;77:87–111. doi:10.1016/S0070-2153(06)77004-1

3. Gorman GS, Chinnery PF, DiMauro S, et al. Mitochondrial diseases. Nat Rev Dis Primers. 2016;2:16080. doi:10.1038/nrdp.2016.80

4. Naidoo K, Khathi A. The potential role of gossypetin in the treatment of diabetes mellitus and its associated complications: a review. Int J Mol Sci. 2023;24(24):17609. doi:10.3390/ijms242417609

5. Rötig A, Bonnefont JP, Munnich A. Mitochondrial diabetes mellitus. Diabetes Metab. 1996;22(5):291–298.

6. Alharbi KK, Alsaikhan AS, Alshammary AF, et al. Screening of mitochondrial mutations in Saudi women diagnosed with gestational diabetes mellitus: a non-replicative case-control study. Saudi J Biol Sci. 2022;29(1):360–365. doi:10.1016/j.sjbs.2021.08.102

7. Montgomery MK. Mitochondrial dysfunction and diabetes: is mitochondrial transfer a friend or foe? Biology. 2019;8(2):33. doi:10.3390/biology8020033

8. So S, Lee S, Lee Y, et al. Dysfunctional pancreatic cells differentiated from induced pluripotent stem cells with mitochondrial DNA mutations. BMB Rep. 2022;55(9):453–458. doi:10.5483/BMBRep.2022.55.9.023

9. Tabebi M, Safi W, Felhi R, et al. The first concurrent detection of mitochondrial DNA m.3243A>G mutation, deletion, and depletion in a family with mitochondrial diabetes. Mol Genet Genomic Med. 2020;8(7):e1292. doi:10.1002/mgg3.1292

10. Fukao T, Kondo M, Yamamoto T, et al. Comparison of mitochondrial A3243G mutation loads in easily accessible samples from a family with maternally inherited diabetes and deafness. Mol Med Rep. 2009;2(1):69–72. doi:10.3892/mmr_00000063

11. Ding Y, Zhang S, Guo Q, et al. Mitochondrial Diabetes is Associated with tRNALeu(UUR) A3243G and ND6 T14502C Mutations. Diabetes Metab Syndr Obes. 2022;15:1687–1701. doi:10.2147/DMSO.S363978

12. Finsterer J. Genetic, pathogenetic, and phenotypic implications of the mitochondrial A3243G tRNALeu(UUR) mutation. Acta Neurol Scand. 2007;116(1):1–14. doi:10.1111/j.1600-0404.2007.00836.x

13. Yang L, Guo Q, Leng J, et al. Late onset of type 2 diabetes is associated with mitochondrial tRNATrp A5514G and tRNASer(AGY) C12237T mutations. J Clin Lab Anal. 2022;36(1):e24102. doi:10.1002/jcla.24102

14. Wang M, Liu H, Zheng J, et al. A deafness- and diabetes-associated tRNA mutation causes deficient pseudouridinylation at position 55 in tRNAGlu and mitochondrial dysfunction. J Biol Chem. 2016;291(40):21029–21041. doi:10.1074/jbc.M116.739482

15. Jiang Z, Cai X, Kong J, et al. Maternally transmitted diabetes mellitus may be associated with mitochondrial ND5 T12338C and tRNAAla T5587C variants. Ir J Med Sci. 2022;191(6):2625–2633. doi:10.1007/s11845-021-02911-w

16. Yarham JW, Elson JL, Blakely EL, et al. Mitochondrial tRNA mutations and disease. Wiley Interdiscip Rev RNA. 2010;1(2):304–324. doi:10.1002/wrna.27

17. Bennett RL, French KS, Resta RG, et al. Standardized human pedigree nomenclature: update and assessment of the recommendations of the national society of genetic counselors. J Genet Couns. 2008;17(5):424–433. doi:10.1007/s10897-008-9169-9

18. American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2010;33:S62- S69.

19. Libman IM, Barinas-Mitchell E, Bartucci A, et al. Reproducibility of the oral glucose tolerance test in overweight children. J Clin Endocrinol Metab. 2008;93(11):4231–4237. doi:10.1210/jc.2008-0801

20. Hypertension Study Group of Chinese Society of Cardiology of Chinese Medical Association. Chinese expert consensus on obesity-related hypertension management, Zhonghua Xin Xue Guan Bing Za Zhi. 2016;44(3):212–219. Article in Chinese. doi:10.3760/cma.j.issn.0253-3758.2016.03.006

21. Chailurkit LO, Jongjaroenprasert W, Chanprasertyothin S, et al. Insulin and C-peptide levels, pancreatic beta cell function, and insulin resistance across glucose tolerance status in Thais. J Clin Lab Anal. 2007;21(2):85–90. doi:10.1002/jcla.20138

22. Lin KH, Liou TL, Hsiao LC, et al. Clinical and biochemical indicators of homeostasis model assessment-estimated insulin resistance in postmenopausal women. J Chin Med Assoc. 2011;74(10):442–447. doi:10.1016/j.jcma.2011.08.014

23. Levey AS, Stevens LA, Schmid CH, et al. CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration). A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–612. doi:10.7326/0003-4819-150-9-200905050-00006

24. Zhang J, Lu B, Xia WW, et al. The mitochondrial transfer RNAAsp A7551G mutation may contribute to the clinical expression of deafness associated with the A1555G mutation in a pedigree with hearing impairment. Mol Med Rep. 2019;19(3):1797–1802. doi:10.3892/mmr.2018.9790

25. Wang H, Barket B, Du S, et al. The prevalence and correlates of vision impairment and glasses ownership among ethnic minority and Han schoolchildren in rural China. PLoS One. 2021;16(8):e0256565. doi:10.1371/journal.pone.0256565

26. Ding Y, Ye YF, Li MY, et al. Mitochondrial tRNAAla 5601C>T variant may affect the clinical expression of the LHON‑related ND4 11778G>A mutation in a family. Mol Med Rep. 2020;21(1):201–208. doi:10.3892/mmr.2019.10844

27. Andrews RM, Kubacka I, Chinnery PF, et al. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999;23(2):147. doi:10.1038/13779

28. Ding Y, Yu J, Guo Q, et al. Molecular characterization of two Chinese pedigrees with maternally inherited hypertension. J Gene Med. 2021;23(4):e3328. doi:10.1002/jgm.3328

29. Ruiz-Pesini E, Wallace DC. Evidence for adaptive selection acting on the tRNA and rRNA genes of human mitochondrial DNA. Hum Mutat. 2006;27(11):1072–1081. doi:10.1002/humu.20378

30. Chomyn A, Lai ST, Shakeley R, et al. Platelet-mediated transformation of mtDNA-less human cells: analysis of phenotypic variability among clones from normal individuals--and complementation behavior of the tRNALys mutation causing myoclonic epilepsy and ragged red fibers. Am J Hum Genet. 1994;54(6):966–974.

31. Gong S, Peng Y, Jiang P, et al. A deafness-associated tRNAHis mutation alters the mitochondrial function, ROS production and membrane potential. Nucleic Acids Res. 2014;42(12):8039–8048. doi:10.1093/nar/gku466

32. Jia Z, Zhang Y, Li Q, et al. A coronary artery disease-associated tRNAThr mutation altered mitochondrial function, apoptosis and angiogenesis. Nucleic Acids Res. 2019;47(4):2056–2074. doi:10.1093/nar/gky1241

33. Ding Y, Xia BH, Zhang CJ, et al. Mitochondrial tRNALeu(UUR) C3275T, tRNAGln T4363C and tRNALys A8343G mutations may be associated with PCOS and metabolic syndrome. Gene. 2018;642:299–306. doi:10.1016/j.gene.2017.11.049

34. Kong QP, Bandelt HJ, Sun C, et al. Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations. Hum Mol Genet. 2006;15(13):2076–2086. doi:10.1093/hmg/ddl130

35. Bibb MJ, Van Etten RA, Wright CT, et al. Sequence and gene organization of mouse mitochondrial DNA. Cell. 1981;26(2 Pt 2):167–180. doi:10.1016/0092-8674(81)90300-7

36. Gadaleta G, Pepe G, De Candia G, et al. The complete nucleotide sequence of the Rattus norvegicus mitochondrial genome: cryptic signals revealed by comparative analysis between vertebrates. J Mol Evol. 1989;28(6):497–516. doi:10.1007/BF02602930

37. Roe BA, Ma DP, Wilson RK, et al. The complete nucleotide sequence of the xenopus laevis mitochondrial genome. J Biol Chem. 1985;260(17):9759–9774.

38. Florentz C, Sohm B, Tryoen-Tóth P, et al. Human mitochondrial tRNAs in health and disease. Cell Mol Life Sci. 2003;60(7):1356–1375. doi:10.1007/s00018-003-2343-1

39. Rollins B, Martin MV, Sequeira PA, et al. Mitochondrial variants in schizophrenia, bipolar disorder, and major depressive disorder. PLoS One. 2009;4(3):e4913. doi:10.1371/journal.pone.0004913

40. Paumard P, Vaillier J, Coulary B, et al. The ATP synthase is involved in generating mitochondrial cristae morphology. EMBO J. 2002;21(3):221–230. doi:10.1093/emboj/21.3.221

41. Zorova LD, Popkov VA, Plotnikov EY, et al. Mitochondrial membrane potential. Anal Biochem. 2018;552:50–59. doi:10.1016/j.ab.2017.07.009

42. Zorov DB, Juhaszova M, Sollott SJ. Mitochondrial ROS-induced ROS release: an update and review. Biochim Biophys Acta. 2006;1757(5–6):509–517. doi:10.1016/j.bbabio.2006.04.029

43. Swan EJ, Salem RM, Sandholm N, et al.; GENIE Consortium. Genetic risk factors affecting mitochondrial function are associated with kidney disease in people with Type 1 diabetes. Diabet Med. 2015;32(8):1104–1109. doi:10.1111/dme.12763

44. Langenberg C, Lotta LA. Genomic insights into the causes of type 2 diabetes. Lancet. 2018;391(10138):2463–2474. doi:10.1016/S0140-6736(18)31132-2

45. Mulder H. Transcribing β-cell mitochondria in health and disease. Mol Metab. 2017;6(9):1040–1051. doi:10.1016/j.molmet.2017.05.014

46. Hesselink MK, Schrauwen-Hinderling V, Schrauwen P. Skeletal muscle mitochondria as a target to prevent or treat type 2 diabetes mellitus. Nat Rev Endocrinol. 2016;12(11):633–645. doi:10.1038/nrendo.2016.104

47. Li R, Liu Y, Li Z, et al. Failures in mitochondrial tRNAMet and tRNAGln metabolism caused by the novel 4401A>G mutation are involved in essential hypertension in a Han Chinese family. Hypertension. 2009;54(2):329–337. doi:10.1161/HYPERTENSIONAHA.109.129270

48. Yin Z, Burger N, Kula-Alwar D, et al. Structural basis for a complex I mutation that blocks pathological ROS production. Nat Commun. 2021;12(1):707. doi:10.1038/s41467-021-20942-w

49. McManus MJ, Picard M, Chen HW, et al. Mitochondrial DNA variation dictates expressivity and progression of nuclear DNA mutations causing cardiomyopathy. Cell Metab. 2019;29(1):78–90.e5. doi:10.1016/j.cmet.2018.08.002

50. Vithayathil SA, Ma Y, Kaipparettu BA. Transmitochondrial cybrids: tools for functional studies of mutant mitochondria. Methods Mol Biol. 2012;837:219–230. doi:10.1007/978-1-61779-504-6_15

51. Paliwal S, Fiumera AC, Fiumera HL. Mitochondrial-nuclear epistasis contributes to phenotypic variation and coadaptation in natural isolates of Saccharomyces cerevisiae. Genetics. 2014;198(3):1251–1265. doi:10.1534/genetics.114.168575

52. Chang YH, Lin HY, Shen FC, et al. The causal role of mitochondrial dynamics in regulating innate immunity in diabetes. Front Endocrinol (Lausanne). 2020;11:445. doi:10.3389/fendo.2020.00445

53. Wilkins HM, Carl SM, Swerdlow RH. Cytoplasmic hybrid (cybrid) cell lines as a practical model for mitochondriopathies. Redox Biol. 2014;2:619–631. doi:10.1016/j.redox.2014.03.006

54. Houstek J, Klement P, Floryk D, et al. A novel deficiency of mitochondrial ATPase of nuclear origin. Hum Mol Genet. 1999;8(11):1967–1974. doi:10.1093/hmg/8.11.1967

55. Rossignol R, Faustin B, Rocher C, et al. Mitochondrial threshold effects. Biochem J. 2003;370(Pt 3):751–762. doi:10.1042/BJ20021594

56. Bertram R, Gram Pedersen M, Luciani DS, et al. A simplified model for mitochondrial ATP production. J Theor Biol. 2006;243(4):575–586. doi:10.1016/j.jtbi.2006.07.019

57. Giacomello M, Pyakurel A, Glytsou C, et al. The cell biology of mitochondrial membrane dynamics. Nat Rev Mol Cell Biol. 2020;21(4):204–224. doi:10.1038/s41580-020-0210-7

58. Widlansky ME, Wang J, Shenouda SM, et al. Altered mitochondrial membrane potential, mass, and morphology in the mononuclear cells of humans with type 2 diabetes. Transl Res. 2010;156(1):15–25. doi:10.1016/j.trsl.2010.04.001

59. Xiong S, Mu T, Wang G, et al. Mitochondria-mediated apoptosis in mammals. Protein Cell. 2014;5:737–749. doi:10.1007/s13238-014-0089-1

60. Circu ML, Aw TY. Reactive oxygen species, cellular redox systems, and apoptosis. Free Radic Biol Med. 2010;48(6):749–762. doi:10.1016/j.freeradbiomed.2009.12.02

61. Maassen JA, ‘T Hart LM, Van Essen E, et al. Mitochondrial diabetes: molecular mechanisms and clinical presentation. Diabetes. 2004;53(Suppl 1):S103–9. doi:10.2337/diabetes.53.2007.s103

62. Anello M, Lupi R, Spampinato D, et al. Functional and morphological alterations of mitochondria in pancreatic beta cells from type 2 diabetic patients. Diabetologia. 2005;48(2):282–289. doi:10.1007/s00125-004-1627-9

63. Nicholas LM, Valtat B, Medina A, et al. Mitochondrial transcription factor B2 is essential for mitochondrial and cellular function in pancreatic β-cells. Mol Metab. 2017;6(7):651–663. doi:10.1016/j.molmet.2017.05.005