Whole mitochondrial genome analysis in Chinese patients with keratoconus

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.