Molecular Vision 2021; 27:270-282 <http://www.molvis.org/molvis/v27/270>
Received 27 August 2020 | Accepted 06 May 2021 | Published 08 May 2021

Whole mitochondrial genome analysis in Chinese patients with keratoconus

Liyan Xu, Kaili Yang, Qi Fan, Dongqing Zhao, Chenjiu Pang, Shengwei Ren

The first two authors contributed equally to this work

Henan Provincial People’s Hospital, Henan Eye Hospital, Henan Eye Institute, People’s Hospital of Zhengzhou University, Henan University People’s Hospital, Zhengzhou, China

Correspondence to: Shengwei Ren, Henan Provincial People’s Hospital, Henan Eye Hospital, Henan Eye Institute, People’s Hospital of Zhengzhou University, Henan University People’s Hospital, Zhengzhou, 450003, China; Phone: +86037165580908; email: shengweiren1984@163.com; ysgzz2018@163.com

Abstract

Purpose: Keratoconus (KC) is a corneal disorder characterized by corneal ectasia, progressive corneal thinning, and conical protrusion. This study aimed to elucidate the mitochondrial gene profile in Chinese patients with KC, analyze the mitochondrial haplogroup and heteroplasmy, and further explore the association between mitochondrial genes and KC.

Methods: Mitochondrial sequencing was conducted on 100 patients with KC and 100 matched controls. Haplogroup analysis was conducted with logistic regression analysis. The heteroplasmy was analyzed with ANOVA (ANOVA) and Student t test. Sequence kernel association tests (SKATs) were performed to analyze the association between mitochondrial genes and KC. Mtoolbox, Mitoclass.1, and APOGEE were used to estimate the impact of the identified variants in protein-coding genes. PON-mt-tRNA was used to annotate the impact of the variants in tRNA. RNAstructure was used to predict the secondary structures of native and mutated tRNAs.

Results: We identified 689 variants in patients with KC and 725 variants in controls (with 308 variants shared by both). The mitochondrial haplogroups exhibited no statistically significant differences between the two groups. Based on the heteroplasmy analysis, the number of heteroplasmic variants in the complete mitochondrial genome, RNA coding regions, and noncoding regions were statistically significantly different in the KC cases and controls (p<0.05). The heteroplasmic levels of the m.16180_16182delAA, m.16182insC, and m.14569 G>C variants in the KC cases were statistically significantly higher than those in the controls (p<0.05). The SKAT analysis showed that the COX3 and TRNH genes were statistically significantly associated with KC (p<0.05). Among the nine variants of COX3 included in the SKAT analysis (m.9300G>A, m.9316T>C, m.9327A>G, m.9355A>G, m.9468A>G, m.9612G>A, m.9804G>A, m.9957G>A, and m.9966 G>A), m.9612G>A was predicted to be deleterious by Mtoolbox. The m.9316T>C, m.9327A>G, m.9355A>G, m.9612G>A, m.9804G>A, and m.9957G>A variants were predicted to be damaging by Mitoclass.1. The m.9355A>G and m.9804G>A variants were predicted to be pathogenic by APOGEE. All identified variants located in TRNH (m.12153C>T, m.12178C>T, and m.12192G>A) were predicted to be neutral by the PON-mt-tRNA website.

Conclusions: This study presents the mitochondrial gene profile of Chinese patients with KC and demonstrated that the COX3 and TRNH genes were associated with KC.

Introduction

Keratoconus (KC) is a corneal disorder characterized by corneal ectasia, progressive corneal thinning, and conical protrusion [1], with a reported prevalence of 1.38 in 1,000 individuals [2]. The disease is progressive and requires keratoplasty treatment at an advanced stage [3]. Several studies have indicated that heredity plays an important role in the process of KC [4,5]. Genetic studies on the nuclear genome have identified numerous loci associated with KC [4,6]. However, these loci explain only a small proportion of the heritability, and additional heritable factors have yet to be discovered. As circular DNA independent from the nucleus, the mitochondrial genome also plays a vital role in the pathogenesis of KC [7,8].

Currently, the molecular pathogenesis of KC is poorly understood. Recent studies suggested that oxidative stress could be involved in the pathogenesis of the disease [9,10]. Variations in mtDNA have the potential to alter mitochondrial function, affect the generation of reactive oxygen species (ROS), and lead to increased oxidative stress [11]. Abu-Amero et al. [12] found that mitochondrial sequences were altered in Saudi patients with KC, which indicated that mtDNA mutation potentially contributed to KC progression through the oxidative stress mechanism. Moreover, mitochondrial complex I gene analysis in Indian patients with KC demonstrated that nonsynonymous mtDNA sequence variations might account for elevated ROS production, leading to oxidative stress and finally, resulting in the pathogenesis of KC [13]. Hao et al. [14] found that KC corneas had increased mtDNA damage and higher ROS levels compared to normal corneas in Chinese populations. In addition, Hao et al. [15] revealed that the mtDNA copy number might contribute to the pathogenesis of KC. Nevertheless, the mitochondrial genome profile has not been reported in Chinese patients with KC. In addition, no study has comprehensively elucidated the relationship between mtDNA variants and risks of KC in the Chinese population.

In the present study, whole mitochondrial genome sequencing was used to elucidate the mitochondrial gene profile in Chinese patients with KC. Association analysis was also conducted to illustrate the relationship between mtDNA genes and the risks of KC.

Methods

Study population

A total of 100 patients with KC were recruited from Henan Eye Hospital & Henan Eye Institute. The diagnosis of KC was based on the following inclusion criteria: Corneal topography revealed an asymmetric bowtie pattern with or without skewed axes or signs of KC detected via a slit-lamp examination, such as localized stromal thinning, conical protrusion, Vogt’s striae, Fleischer’s ring, or anterior stromal scar [5]. Cases secondary to trauma, surgery, and other diseases were excluded from the study. The mean corneal curvature (steep keratometry, Ks) of all patients with KC (serious eye for each patient) was 56.3±8.30 D. One hundred matched individuals without any ocular disorder or previous ophthalmic surgeries were enrolled as controls. All the control subjects demonstrated a clear cornea on a slit-lamp examination, and a topographic map was within normal limits. The study adhered to the tenets of the Declaration of Helsinki and the ARVO statement on human subjects.

The study protocol was approved by the ethical committees of Henan Eye Hospital [ethical approval number: HNEECKY-2019 (5)]. Parents were informed of the study purpose for cases that were less than 18 years old. Written informed consent was obtained from all subjects.

Mitochondrial genome sequencing

EDTA anti-coagulated venous blood samples were collected from each subject in the current study. And the blood samples were preserved at −80 °C cryogenic refrigerator before use. Total DNA was extracted from peripheral blood samples with QIAamp DNA Blood kits (Qiagen, Hilden, Germany) according to the manufacturer’s recommendations. The DNA samples were quantified with NanoDrop 2000 and ascertained with electrophoresis.

The mitochondria genome was sequenced by Shanghai Genesky Biotechnologies. The mtDNA was amplified with six pairs of primers with long-range PCR. The six partially overlapping PCR products (about 6 kb each) were pooled and then purified with Agencourt AMPure XP-PCR Purification kits (Beckman, Germany). Then the PCR products were fragmented into 100–500 bp pieces with ultrasound (ME220, Woburn, MA). The fragmented DNA was subjected to end-repairing, dA-tailing, and adaptor ligation with the NEBNext® DNA Library Prep Reagent Set from Illumina® (New England Biolabs, Ipswich, MA). The adaptors ligated to each sample contained different indexes, which we subsequently used for pooled samples testing. Then the PCR reaction was performed in a total volume of 50 μl consisting of 15 μl genomic DNA with indexes, 5 μl pair end primers mix (Illumina, San Diego, CA), 25μl NEBNext UltrallQ5 Master mix (New England Biolabs) and 5 μl ddH2O. The cycling condition was 98 °C for 30 s, followed by 10 cycles of amplification at 98 °C for 10 s, 65 °C for 75 s, and a final extension at 65 °C for 60 s. After PCR amplification, the amplicons were size-checked and quantitated with BioAnalyzer 2100, and then subjected to 2 × 150 bp paired-end massively parallel sequencing using Illumina NovaSeq System (Illumina).

Sequencing data analysis and functional annotation

Fastq data were acquired for a preliminary quality check. Burroughs-Wheeler Aligner (BWA) [16] was used to align the sequence data to the mitochondrial reference sequence (revised Cambridge Reference Sequence, rCRS) and generate sequence alignment files. Variant calling was performed using the Genome Analysis Toolkit (GATK) Mutect2 and Haplotype Caller [17]. Then variants were filtered according to the threshold value (mapping quality ≥20 and base quality ≥20). Haplogroup analysis was conducted on https://haplogrep.i-med.ac.at/app/index.html. Variants exhibiting frequency values within the 0.1–0.9 range were considered heteroplasmic. Moreover, the mean number of heteroplasmic mtDNA variants was obtained by the average number of individual heteroplasmic variants present in the specified location from all samples within the group. The symbols and full names of the total 37 genes in human mitochondria are listed in Appendix 1. The schematic representation of the sequencing data analytical workflow is shown in Figure 1. For associated variants located in the protein-coding gene, Mtoolbox [18], Mitoclass.1 [19], and APOGEE [20] were applied to estimate the impact of the identified variants. For associated variants located in the tRNA gene, PON-mt-tRNA [21] was used to annotate the impact of the variants. Moreover, RNAstructure [22] was used to predict the secondary structures of native and mutated tRNAs.

Statistical analyses

Statistical analyses were performed using the SPSS 21.0 package (SPSS Inc., Chicago, IL). Age was expressed as mean ± standard deviation and analyzed with Student t test, whereas gender was analyzed with the chi-square test. Haplogroup analysis was conducted through logistic regression analysis. The number of heteroplasmic variants was compared through ANOVA (ANOVA) and Student t test. The Mann–Whitney U test was used to compare the heteroplasmic level between KC cases and controls in variants with heteroplasmic samples ≥3 in each group. Sequence kernel association test (SKAT) was conducted with the software Efficient and Parallelizable Association Container Toolbox v.3.2.6 (EPACTS v.3.2.6) [23]. SKAT, which applies a multivariate regression model to uncover genes associated with a specific disease, could evaluate the cumulative effect of rare and common variants at the gene level [23]. Variants with a genotype missing rate greater than 0.1 were included in the SKAT. P values less than 0.05 were considered statistically significant.

Results

Spectrum of mitochondrial genome variants

Whole mitochondrial sequencing was performed in 100 patients with KC (73 men and 27 women, mean age of 19.7±3.50 years) and 100 controls (72 men and 28 women, mean age of 19.1±2.40 years) to explore the profile of the mitochondria genome and the role of mtDNA genes contributing to the susceptibility to KC. No statistically significant difference was found in the mean age and gender distribution between patients with KC and controls (all p>0.05). The sequencing data were generated from all subjects with a mean read depth of 3,139.50±524.160 (range from 2,122.15 to 6,140.10). As shown in Figure 2A, 368 variants were shared between patients with KC and controls, 321 were found only in patients with KC, and 357 were found only in controls. The distribution of the variants in the protein-coding genes identified in the study population is shown in Figure 2B. The results showed that the ND5 (OMIM:516005; Gene ID:4540) gene contained 107 variants, followed by the CYTB (OMIM:516020; Gene ID:4519), COX1 (OMIM:516030; Gene ID:4512), and ND4 genes. Furthermore, we categorized the variants in the protein-coding genes in patients with KC and controls separately. A total of 68 variants were located in the ND5 gene, followed by 54 variants in the CYTB gene in patients with KC. The ND5 gene contained 68 variants in the control group, followed by 57 variants in the CYTB gene (Figure 2D, Appendix 2).

Distribution of synonymous and nonsynonymous variants

When we categorized all the variants according to their effects, 462 synonymous variants were identified, while there were 201 nonsynonymous variants (Figure 2C). Among all variants in protein-coding genes, there were more synonymous variants than nonsynonymous variants except the ATP6 (OMIM:516003; Gene ID: 4538) gene (Figure 2D). Among all the nonsynonymous variants, 135 variants were identified in patients with KC, whereas 130 were identified in controls (Figure 2E). No statistically significant difference was found by comparing the number of synonymous and nonsynonymous variants between the two groups (p=0.316). Afterward, we compared the number of nonsynonymous variants in protein-coding genes identified in the two groups. Of note, the ATP6, ATP8 (OMIM:516060; Gene ID: 4508), COX3 (OMIM:516050; Gene ID: 4514), ND2 (OMIM:516001; Gene ID: 4536), ND3 (OMIM: 516002; Gene ID: 4537), ND4L (OMIM:516004; Gene ID: 4539), ND5, ND6 (OMIM:516006; Gene ID: 4541) genes had higher number of nonsynonymous variants in KC patients, whereas the COX1, COX2 (OMIM:516040; Gene ID: 4513), CYTB, ND1 (OMIM:516000; Gene ID: 4535), and ND4 genes had lower number of nonsynonymous variants (Figure 2F, Appendix 3).

Mitochondrial haplogroup analysis

Mitochondrial haplogroups were analyzed in 100 patients with KC and 100 controls. The distribution of the mitochondrial haplogroups in the KC cases and controls is shown in Table 1. No statistically significant difference was found between KC patients and controls for all mitochondrial haplogroups tested (p>0.05).

Mitochondrial heteroplasmy analysis

For this analysis, heteroplasmic variants were compared according to different mitochondrial regions in the KC cases and controls (Figure 3). Among the protein-coding regions, RNA coding regions, and noncoding regions, the noncoding regions had the highest number of heteroplasmic variants in the two groups. By comparing the number of heteroplasmic variants in different mitochondrial regions, the average number of total heteroplasmic variants in the KC cases was higher than that in controls (p<0.05) in the complete mitochondrial genome. Regarding heteroplasmic variants in the protein-coding regions, no statistically significant differences was found between the KC cases and the controls (p>0.05). However, the number of heteroplasmic variants in the RNA coding regions was lower than that in the controls (p<0.05). The number of heteroplasmic variants in the noncoding regions was higher than that in the controls (p<0.05). After comparing the heteroplasmic levels between the KC cases and the controls, we found that the heteroplasmic levels of the m.16180_16182delAA, m.16182insC and m.14569 G>C variants in the KC cases were significantly higher than those in the controls (p<0.05, Appendix 4).

Association of mitochondrial genes with KC

Given the observation that synonymous mutant alleles had a neutral effect on protein function and disease susceptibility, we selected nonsynonymous variants in the protein-coding genes for further analysis. The SKAT was applied to assess the contribution of nonsynonymous variants and variants over rRNA and tRNA genes in regard to KC susceptibility. The results revealed that the COX3 and TRNH (OMIM:590040; Gene ID: 4564) genes were statistically significantly associated with KC (p<0.05; Table 2). Nine variants in the COX3 gene and three variants in the TRNH gene were included in the SKAT analysis. Functional annotations of these variants are listed in Table 3 and Table 4. Among the nine variants located in the COX3 gene (m.9300G>A, m.9316T>C, m.9327A>G, m.9355A>G, m.9468A>G, m.9612G>A, m.9804G>A, m.9957G>A, and m.9966 G>A), m.9612G>A was predicted to be deleterious by Mtoolbox. The m.9316T>C, m.9327A>G, m.9355A>G, m.9612G>A, m.9804G>A, and m.9957G>A variants were predicted to be damaging by Mitoclass.1. The m.9355A>G and m.9804G>A variants were predicted to be pathogenic by APOGEE. All three variants located in the TRNH gene (m.12153C>T, m.12178C>T, and m.12192G>A) were predicted to be neutral via the PON-mt-tRNA website. The secondary structures of native tRNAHis and mutated tRNAHis predicted with RNAstructure are presented in Figure 4.

Discussion

In the present study, we provided a comprehensive understanding of the mitochondrial genome profile in Chinese patients with KC. Haplogroup analysis of the KC cases and controls revealed no statistically significant difference. In addition, heteroplasmy analysis demonstrated that the number of heteroplasmic variants in RNA coding regions, noncoding regions, and complete mitochondrial genome were statistically significantly different between the two groups. Heteroplasmic level analyses indicated the heteroplasmic levels of the m.16180_16182delAA, m.16182insC, and m.14569 G>C variants in the KC cases were statistically significantly higher than those in the controls. Of note, the COX3 and TRNH genes were statistically significantly associated with KC according to the gene-based SKAT analysis.

Previous studies identified several mitochondrial genome variants in KC [12,13]. However, the mitochondrial genome profile in Chinese patients with KC was unclear, and no studies have comprehensively elucidated the relationship between mtDNA genes or variants and risks of KC in the Chinese population. Therefore, we conducted whole mitochondrial sequencing to elucidate the genetic profile of mitochondria in Chinese patients with KC and analyzed the association between mitochondrial genes or variants and KC. Among the variants identified in the KC cases, the ND5 gene had the highest number of variants in mitochondrial protein-coding genes, which was consistent with the findings of Pathak et al. [13]. Generally, there were more synonymous variants than nonsynonymous variants [24]. In this study, more synonymous variants were identified in both groups. Comparison of synonymous and nonsynonymous variants between the KC cases and the controls revealed no statistically significant difference. Nevertheless, Pathak et al. [13] found a higher number of nonsynonymous variants in patients with KC than in controls in their study. The controversial results could be attributed to the differences in populations and sample sizes.

Previous studies have identified an association between mitochondrial haplogroups and various diseases, including ophthalmic diseases such as glaucoma [25,26] and Leber hereditary optic neuropathy (LHON) [27]. In the present study, we evaluated the association between mitochondrial haplogroup and KC. However, no statistically significant differences were found for all mitochondrial haplogroups tested. Abu-Amero et al. [28] found that mitochondrial haplogroups H and R were statistically significantly associated with KC in Saudi patients with KC. Thus, the differences in populations and sample sizes potentially account for the controversial results.

At present, there is a paucity of information regarding the extent of mtDNA heteroplasmy in patients with KC and controls. We found a relatively high number of heteroplasmic mtDNA variants in the noncoding regions of both groups, which might be attributed to various factors, including the variable D-loop within noncoding regions and negative selective pressures associated with mtDNA variation in protein and RNA coding regions [29,30]. Moreover, the number of heteroplasmic variants in RNA coding regions and noncoding regions were statistically significantly different between the KC cases and controls, implicating that the RNA coding and noncoding regions may influence KC. In addition, we found the heteroplasmic levels of the m.16180_16182delAA, m.16182insC, and m.14569G>C variants were higher in the KC cases, indicating that those variants might have a potential role in the disease. However, the accurate molecular mechanisms remain unclear. Further exploration of the potential role of mtDNA heteroplasmy on KC is warranted.

Recent studies have shown that rare and common variants can contribute to the same disease [31,32]. SKAT, which applies a multivariate regression model to uncover genes associated with a specific disease, could evaluate the cumulative effect of rare and common variants at the gene level [23]. As synonymous variants have a neutral impact on protein function and disease susceptibility, we focused on nonsynonymous variants in mitochondrial protein-coding genes for further analysis. After analyzing the contribution of nonsynonymous variants in mitochondrial protein-coding genes and variants over rRNA and tRNA genes to KC susceptibility with SKAT at the gene level, we found the COX3 and TRNH genes may contribute to the occurrence of KC. The mitochondrial genome harbors 13 protein-coding genes, all of which encode the subunits of the electron transport chain, where oxidation phosphorylation (OXPHOS) occurs, and are responsible for most cell ATP production [33]. MtDNA variants often disrupt OXPHOS, compromising ATP production and releasing ROS. Cytochrome c oxidase III, encoded by the COX3 gene, is a subunit of respiratory complex IV. Notably, respiratory complex IV plays an essential role in the electron transport chain [33]. Therefore, we speculated that the abnormality of the COX3 gene might disrupt OXPHOS by altering the electron transport chain, consequently resulting in KC. The TRNH gene encodes tRNAHis which participates in protein translation. Gong et al. [34] constructed cybrids of a mutation in TRNH by transferring mitochondria from lymphoblastoid cell lines. After detecting the levels of mitochondrial ATP and membrane potential in mutant cybrids, Gong et al. found the mutant cybrids had increased ROS levels, inferring that the mutation in TRNH could lead to the increasing production of ROS. In addition, KC fibroblasts had increased basal generation levels of ROS [35]. Thus, we inferred that variants in the TRNH gene might elevate ROS levels, leading to the occurrence of KC.

In the present study, nine variants in COX3 and three variants in TRNH were included in the SKAT analysis. Among the nine variants located in the COX3 gene, the m.9300G>A variant was reported to have a potential role in some patients with sudden unexpected infant death [36]. The m.9804G>A variant has been reported in LHON, tetralogy of Fallot (TOF), mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS), and open-angle glaucoma (OAG) [37-40]. The m.9957G>A variant has been reported in cases of MELAS, non-arteritic ischemic optic neuropathy (NAION), chronic progressive external ophthalmoplegia (CPEO), hypertrophic cardiomyopathy (HCM), and gout [41-45]. Among the variants located in the TRNH gene, m.12153C>T was reported in patients with thiamine deficiency–accompanied mitochondrial myopathy [46], m.12178C>T was previously reported in Saudi patients with KC [12], and m.12192G>A was reported in cardiomyopathy (CM) and deafness [47,48]. The predicted secondary structures of m.12153C>T and m.12192G>A were similar to native tRNAHis, whereas the predicted structure of m.12178C>T was statistically significantly different from native tRNAHis. This discrepancy of predicted structures indicates that the m.12178C>T variant might critically influence the occurrence of KC. Previously, Abu-Amero et al. [12] screened the mitochondrial genome in Saudi patients with KC and found nine potentially pathogenic variants in the study population. However, only m.12178C>T was identified in the present study, which is in accordance with Abu-Amero et al. findings. Of note, this might be attributed to different populations and different sample sizes, which warrants future studies with multicenter and larger population sizes for in-depth exploration.

In conclusion, this study presented the mitochondrial gene profile of Chinese patients with KC and identified that the COX3 and TRNH genes were associated with KC.

Appendix 1. Gene symbols and names of total 37 genes in human mitochondria.

Appendix 2. Number of variants by different genome regions in KC cases and controls.

Appendix 3. Number of nonsynonymous variants by different genome regions in KC cases and controls.

Appendix 4. Comparison of heteroplasmic levels between KC cases and controls.

Acknowledgments

This study was supported by Open Program of Shandong Provincial Key Laboratory of Ophthalmology (No. 2018–04), Special program for basic research of Henan Eye Institute (No.20JCZD003) and Basic research and Cultivation Foundation for Young Teachers of Zhengzhou University (No. JC202051049). We acknowledge all the subjects for providing blood samples.

References

  1. Rabinowitz YS. Keratoconus. Surv Ophthalmol. 1998; 42:297-319. [PMID: 9493273]
  2. Hashemi H, Heydarian S, Hooshmand E, Saatchi M, Yekta A, Aghamirsalim M, Valadkhan M, Mortazavi M, Hashemi A, Khabazkhoob M. The Prevalence and Risk Factors for Keratoconus: A Systematic Review and Meta-Analysis. Cornea. 2020; 39:263-70. [PMID: 31498247]
  3. Gomes JA, Tan D, Rapuano CJ, Belin MW, Ambrosio R, , Jr Guell JL, Malecaze F, Nishida K, Sangwan VS. Group of Panelists for the Global Delphi Panel of K, Ectatic D. Global consensus on keratoconus and ectatic diseases. Cornea. 2015; 34:359-69. [PMID: 25738235]
  4. Valgaeren H, Koppen C, Van Camp G. A new perspective on the genetics of keratoconus: why have we not been more successful? Ophthalmic Genet. 2018; 39:158-74. [PMID: 29111844]
  5. Mas Tur V, MacGregor C, Jayaswal R, O’Brart D, Maycock N. A review of keratoconus: Diagnosis, pathophysiology, and genetics. Surv Ophthalmol. 2017; 62:770-83. [PMID: 28688894]
  6. Lucas SEM, Burdon KP. Genetic and Environmental Risk Factors for Keratoconus. Annu Rev Vis Sci. 2020; 6:25-46. [PMID: 32320633]
  7. Loukovitis E, Sfakianakis K, Syrmakesi P, Tsotridou E, Orfanidou M, Bakaloudi DR, Stoila M, Kozei A, Koronis S, Zachariadis Z, Tranos P, Kozeis N, Balidis M, Gatzioufas Z, Fiska A, Anogeianakis G. Genetic Aspects of Keratoconus: A Literature Review Exploring Potential Genetic Contributions and Possible Genetic Relationships with Comorbidities. Ophthalmol Ther. 2018; 7:263-92. [PMID: 30191404]
  8. Abu-Amero KK, Al-Muammar AM, Kondkar AA. Genetics of Keratoconus: Where Do We Stand? J Ophthalmol. 2014; 2014:1-11. [PMID: 25254113]
  9. Nita M, Grzybowski A. The Role of the Reactive Oxygen Species and Oxidative Stress in the Pathomechanism of the Age-Related Ocular Diseases and Other Pathologies of the Anterior and Posterior Eye Segments in Adults. Oxid Med Cell Longev. 2016; 2016:3164734 [PMID: 26881021]
  10. Wojcik KA, Kaminska A, Blasiak J, Szaflik J, Szaflik JP. Oxidative stress in the pathogenesis of keratoconus and Fuchs endothelial corneal dystrophy. Int J Mol Sci. 2013; 14:19294-308. [PMID: 24065107]
  11. Vallabh NA, Romano V, Willoughby CE. Mitochondrial dysfunction and oxidative stress in corneal disease. Mitochondrion. 2017; 36:103-13. [PMID: 28549842]
  12. Abu-Amero KK, Azad TA, Kalantan H, Sultan T, Al-Muammar AM. Mitochondrial sequence changes in keratoconus patients. Invest Ophthalmol Vis Sci. 2014; 55:1706-10. [PMID: 24569587]
  13. Pathak D, Nayak B, Singh M, Sharma N, Tandon R, Sinha R, Titiyal JS, Dada R. Mitochondrial complex 1 gene analysis in keratoconus. Mol Vis. 2011; 17:1514-25. [PMID: 21691575]
  14. Hao XD, Chen ZL, Qu ML, Zhao XW, Li SX, Chen P. Decreased Integrity, Content, and Increased Transcript Level of Mitochondrial DNA Are Associated with Keratoconus. PLoS One. 2016; 11:e0165580 [PMID: 27783701]
  15. Hao XD, Chen P, Wang Y, Li SX, Xie LX. Mitochondrial DNA copy number, but not haplogroup is associated with keratoconus in Han Chinese population. Exp Eye Res. 2015; 132:59-63. [PMID: 25613073]
  16. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25:1754-60. [PMID: 19451168]
  17. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20:1297-303. [PMID: 20644199]
  18. Calabrese C, Simone D, Diroma MA, Santorsola M, Gutta C, Gasparre G, Picardi E, Pesole G, Attimonelli M. MToolBox: a highly automated pipeline for heteroplasmy annotation and prioritization analysis of human mitochondrial variants in high-throughput sequencing. Bioinformatics. 2014; 30:3115-7. [PMID: 25028726]
  19. Martin-Navarro A, Gaudioso-Simon A, Alvarez-Jarreta J, Montoya J, Mayordomo E, Ruiz-Pesini E. Machine learning classifier for identification of damaging missense mutations exclusive to human mitochondrial DNA-encoded polypeptides. BMC Bioinformatics. 2017; 18:158 [PMID: 28270093]
  20. Castellana S, Fusilli C, Mazzoccoli G, Biagini T, Capocefalo D, Carella M, Vescovi AL, Mazza T. High-confidence assessment of functional impact of human mitochondrial non-synonymous genome variations by APOGEE. PLOS Comput Biol. 2017; 13:e1005628 [PMID: 28640805]
  21. Niroula A, Vihinen M. PON-mt-tRNA: a multifactorial probability-based method for classification of mitochondrial tRNA variations. Nucleic Acids Res. 2016; 44:2020-7. [PMID: 26843426]
  22. Reuter JS, Mathews DH. RNAstructure: software for RNA secondary structure prediction and analysis. BMC Bioinformatics. 2010; 11:129 [PMID: 20230624]
  23. Ionita-Laza I, Lee S, Makarov V, Buxbaum JD, Lin X. Sequence kernel association tests for the combined effect of rare and common variants. Am J Hum Genet. 2013; 92:841-53. [PMID: 23684009]
  24. Zeng Z, Bromberg Y. Predicting Functional Effects of Synonymous Variants: A Systematic Review and Perspectives. Front Genet. 2019; 10:914 [PMID: 31649718]
  25. Abu-Amero KK, Morales J, Bosley TM, Mohamed GH, Cabrera VM. The role of mitochondrial haplogroups in glaucoma: a study in an Arab population. Mol Vis. 2008; 14:518-22. [PMID: 18385785]
  26. Abu-Amero KK, Cabrera VM, Larruga JM, Osman EA, Gonzalez AM, Al-Obeidan SA. Eurasian and Sub-Saharan African mitochondrial DNA haplogroup influences pseudoexfoliation glaucoma development in Saudi patients. Mol Vis. 2011; 17:543-7. [PMID: 21364909]
  27. Abu-Amero KK, Bosley TM. Mitochondrial abnormalities in patients with LHON-like optic neuropathies. Invest Ophthalmol Vis Sci. 2006; 47:4211-20. [PMID: 17003408]
  28. Abu-Amero KK, Azad TA, Sultan T, Kalantan H, Kondkar AA, Al-Muammar AM. Association of mitochondrial haplogroups H and R with keratoconus in Saudi Arabian patients. Invest Ophthalmol Vis Sci. 2014; 55:2827-31. [PMID: 24722698]
  29. Ye K, Lu J, Ma F, Keinan A, Gu Z. Extensive pathogenicity of mitochondrial heteroplasmy in healthy human individuals. Proc Natl Acad Sci USA. 2014; 111:10654-9. [PMID: 25002485]
  30. Stoneking M. Hypervariable sites in the mtDNA control region are mutational hotspots. Am J Hum Genet. 2000; 67:1029-32. [PMID: 10968778]
  31. Asselbergs FW, Guo Y, van Iperen EP, Sivapalaratnam S, Tragante V, Lanktree MB, Lange LA, Almoguera B, Appelman YE, Barnard J, Baumert J, Beitelshees AL, Bhangale TR, Chen YD, Gaunt TR, Gong Y, Hopewell JC, Johnson T, Kleber ME, Langaee TY, Li M, Li YR, Liu K, McDonough CW, Meijs MF, Middelberg RP, Musunuru K, Nelson CP, O’Connell JR, Padmanabhan S, Pankow JS, Pankratz N, Rafelt S, Rajagopalan R, Romaine SP, Schork NJ, Shaffer J, Shen H, Smith EN, Tischfield SE, van der Most PJ, van Vliet-Ostaptchouk JV, Verweij N, Volcik KA, Zhang L, Bailey KR, Bailey KM, Bauer F, Boer JM, Braund PS, Burt A, Burton PR, Buxbaum SG, Chen W, Cooper-Dehoff RM, Cupples LA, deJong JS, Delles C, Duggan D, Fornage M, Furlong CE, Glazer N, Gums JG, Hastie C, Holmes MV, Illig T, Kirkland SA, Kivimaki M, Klein R, Klein BE, Kooperberg C, Kottke-Marchant K, Kumari M, LaCroix AZ, Mallela L, Murugesan G, Ordovas J, Ouwehand WH, Post WS, Saxena R, Scharnagl H, Schreiner PJ, Shah T, Shields DC, Shimbo D, Srinivasan SR, Stolk RP, Swerdlow DI, Taylor HA, , Jr Topol EJ, Toskala E, van Pelt JL, van Setten J, Yusuf S, Whittaker JC, Zwinderman AH. Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci. Am J Hum Genet. 2012; 91:823-38. [PMID: 23063622]
  32. Diogo D, Kurreeman F, Stahl EA, Liao KP, Gupta N, Greenberg JD, Rivas MA, Hickey B, Flannick J, Thomson B, Guiducci C, Ripke S, Adzhubey I, Barton A, Kremer JM, Alfredsson L. Rare, low-frequency, and common variants in the protein-coding sequence of biological candidate genes from GWASs contribute to risk of rheumatoid arthritis. Am J Hum Genet. 2013; 92:15-27. [PMID: 23261300]
  33. Sousa JS, D’Imprima E, Vonck J. Mitochondrial Respiratory Chain Complexes. Subcell Biochem. 2018; 87:167-227. [PMID: 29464561]
  34. Gong S, Peng Y, Jiang P, Wang M, Fan M, Wang X, Zhou H, Li H, Yan Q, Huang T, Guan MX. A deafness-associated tRNAHis mutation alters the mitochondrial function, ROS production and membrane potential. Nucleic Acids Res. 2014; 42:8039-48. [PMID: 24920829]
  35. Chwa M, Atilano SR, Hertzog D, Zheng H, Langberg J, Kim DW, Kenney MC. Hypersensitive response to oxidative stress in keratoconus corneal fibroblasts. Invest Ophthalmol Vis Sci. 2008; 49:4361-9. [PMID: 18515575]
  36. Opdal SH, Vege A, Egeland T, Musse MA, Rognum TO. Possible role of mtDNA mutations in sudden infant death. Pediatr Neurol. 2002; 27:23-9. [PMID: 12160969]
  37. Dogulu CF, Kansu T, Seyrantepe V, Ozguc M, Topaloglu H, Johns DR. Mitochondrial DNA analysis in the Turkish Leber’s hereditary optic neuropathy population. Eye (Lond). 2001; 15:183-8. [PMID: 11339587]
  38. Tansel T, Pacal F, Ustek D. A novel ATP8 gene mutation in an infant with tetralogy of Fallot. Cardiol Young. 2014; 24:531-3. [PMID: 23735083]
  39. Wani AA, Ahanger SH, Bapat SA, Rangrez AY, Hingankar N, Suresh CG, Barnabas S, Patole MS, Shouche YS. Analysis of mitochondrial DNA sequences in childhood encephalomyopathies reveals new disease-associated variants. PLoS One. 2007; 2:e942 [PMID: 17895983]
  40. Inagaki Y, Mashima Y, Fuse N, Ohtake Y, Fujimaki T, Fukuchi T. Glaucoma Gene Research G. Mitochondrial DNA mutations with Leber’s hereditary optic neuropathy in Japanese patients with open-angle glaucoma. Jpn J Ophthalmol. 2006; 50:128-34. [PMID: 16604388]
  41. Manfredi G, Schon EA, Moraes CT, Bonilla E, Berry GT, Sladky JT, DiMauro S. A new mutation associated with MELAS is located in a mitochondrial DNA polypeptide-coding gene. Neuromuscul Disord. 1995; 5:391-8. [PMID: 7496173]
  42. Abu-Amero KK, Bosley TM, Bohlega S, Hansen E. Mitochondrial T9957C mutation in association with NAION and seizures but not MELAS. Ophthalmic Genet. 2005; 26:31-6. [PMID: 15823923]
  43. Liu CH, Liou CW, Liu CH, Kuo HC, Chu CC, Huang CC. Chronic progressive external ophthalmoplegia with T9957C mitochondrial DNA mutation in a Taiwanese patient. Acta Neurol Taiwan. 2011; 20:53-8. [PMID: 21249588]
  44. Arbustini E, Fasani R, Morbini P, Diegoli M, Grasso M, Dal Bello B, Marangoni E, Banfi P, Banchieri N, Bellini O, Comi G, Narula J, Campana C, Gavazzi A, Danesino C, Vigano M. Coexistence of mitochondrial DNA and beta myosin heavy chain mutations in hypertrophic cardiomyopathy with late congestive heart failure. Heart. 1998; 80:548-58. [PMID: 10065021]
  45. Tseng CC, Chen CJ, Yen JH, Huang HY, Chang JG, Chang SJ, Liao WT. Next-generation sequencing profiling of mitochondrial genomes in gout. Arthritis Res Ther. 2018; 20:137 [PMID: 29976239]
  46. Morovvati S, Nakagawa M, Sato Y, Hamada K, Higuchi I, Osame M. Phenotypes and mitochondrial DNA substitutions in families with A3243G mutation. Acta Neurol Scand. 2002; 106:104-8. [PMID: 12100370]
  47. Shin WS, Tanaka M, Suzuki J, Hemmi C, Toyo-oka T. A novel homoplasmic mutation in mtDNA with a single evolutionary origin as a risk factor for cardiomyopathy. Am J Hum Genet. 2000; 67:1617-20. [PMID: 11038324]
  48. Ding Y, Teng YS, Zhuo GC, Xia BH, Leng JH. The Mitochondrial tRNAHis G12192A Mutation May Modulate the Clinical Expression of Deafness-Associated tRNAThr G15927A Mutation in a Chinese Pedigree. Curr Mol Med. 2019; 19:136-46. [PMID: 30854964]