- Henan Provincial People’s Hospital, Henan Eye Hospital, People’s Hospital of Zhengzhou University, Zhengzhou, China
Background: Mitochondrial DNA (mtDNA) variants have been implicated in keratoconus (KC). The present study aimed to characterize the mtDNA heteroplasmy profile in KC and explore the association of mitochondrial heteroplasmic levels with KC.
Methods: Mitochondrial sequencing of peripheral blood samples and corneal tomography were conducted in 300 KC cases and 300 matched controls. The number of heteroplasmic and homoplasmic variants was calculated across the mitochondrial genome. Spearman’s correlation was used to analyze the correlation between the number of heteroplasmic variants and age. The association of mtDNA heteroplasmic level with KC was analyzed by logistic regression analysis. Moreover, the relationship between mitochondrial heteroplasmic levels and clinical parameters was determined by linear regression analysis.
Results: The distribution of mtDNA heteroplasmic variants showed the highest number of heteroplasmic variants in the non-coding region, while the COX3 gene exhibited the highest number in protein-coding genes. Comparisons of the number of heteroplasmic and homoplasmic non-synonymous variants in protein-coding genes revealed no significant differences between KC cases and controls (all p > 0.05). In addition, the number of heteroplasmic variants was positively associated with age in all subjects (r = 0.085, p = 0.037). The logistic regression analyses indicated that the heteroplasmic levels of m.16180_16181delAA was associated with KC (p < 0.005). Linear regression analyses demonstrated that the heteroplasmic levels of m.16180_16181delAA and m.302A>C were not correlated with thinnest corneal thickness (TCT), steep keratometry (Ks), and flat keratometry (Kf) (all p > 0.05) in KC cases and controls separately.
Conclusion: The current study characterized the mtDNA heteroplasmy profile in KC, and revealed that the heteroplasmic levels of m.16180_16181delAA were associated with KC.
1 Introduction
Keratoconus (KC) is a progressive corneal ectatic disorder characterized by corneal thinning, irregular astigmatism, and vision deterioration (Rabinowitz, 1998). The disease is typically bilateral and usually happens in adolescence (Lucas and Burdon, 2020). It is also one of the most frequent causes of corneal transplantation (Sarezky et al., 2017). Although the exact cause of KC is still unclear in most patients, the genetic, environmental, and behavioral factors were reported to be associated with its pathogenesis (Ferrari and Rama, 2020; Hao et al., 2022; Santodomingo-Rubido et al., 2022). Recently, several studies have reported that oxidative stress plays an important role in the development of KC (Navel et al., 2021; Vallabh et al., 2017), and mitochondria have been implicated in the pathogenesis of KC due to its crucial role in oxidative stress (Roy et al., 2019).
The human mitochondrial genome comprises circular double-stranded DNA that encodes 37 genes, of which thirteen genes are subunits of respiratory complexes, twenty-two genes encode mitochondrial tRNAs, and a further two encode rRNA. Generally, the mutation rate within the mitochondrial genome is much higher than the nuclear genome (Hahn and Zuryn, 2019). In addition, the mutations of mitochondrial DNA (mtDNA) could emerge as homoplasmy or heteroplasmy. The mtDNA homoplasmy represents a uniform type of mtDNA in an individual (Parakatselaki and Ladoukakis, 2021), while the mtDNA heteroplasmy represents two or more types of mtDNA simultaneously existing in the same individual (Stefano et al., 2017). A recent study reported that mitochondrial heteroplasmy is pathogenic and could lead to a greater disease burden in several diseases (Wallace and Chalkia, 2013).
Previous studies have identified multiple mtDNA variants in KC. Eighty-four mtDNA variants, including two novel frameshift mutations in the mitochondrial complex I gene, have been reported in Indian patients with KC (Pathak et al., 2011). Similarly, Abu-Amero et al. (2014b) reported ten non-synonymous mtDNA mutations in KC patients from Saudi Arabia, with one non-synonymous variant heteroplasmic. Our previous study observed an increased number of non-synonymous mtDNA variants in KC patients, though not statistically significant (Xu et al., 2021). Currently, there is a paucity of information regarding mtDNA heteroplasmy in KC from China. Therefore, we aimed to characterize the mtDNA heteroplasmy profile and explore the association of mitochondrial heteroplasmic levels with KC in the present study.
2 Materials and methods
2.1 Study population
A total of 300 KC patients (221 males and 79 females, with a mean age of 20.69 ± 4.68 years) were consecutively recruited between June 2018 and June 2021 in Henan Eye Hospital. The ratio of male to female was 2.80:1 in this study. The ratio was consistent with our previous study (Yang et al., 2021) indicating that the ratio was 2.65:1 in central China, and the study population was more representative. In addition, 300 controls without KC (224 males and 76 females, with a mean age of 20.44 ± 4.21 years) were also enrolled. KC patients were diagnosed based on the following criteria: corneal tomography revealing an asymmetric bowtie pattern with or without skewed axes or keratoconus sign detected by slit lamp examination, such as localized stromal thinning, conical protrusion, Vogt’s striae, Fleischer’s ring or anterior stromal scar (Mas Tur et al., 2017). While for the control group, the slit lamp examination showed normal cornea, and the elevation map in the tomographic map was within the normal limits. All the subjects were unrelated individuals and matched by age and sex with KC cases.
2.2 Clinical examination
All the subjects underwent clinical examinations. Basic characteristics were collected through medical records. The slit lamp examinations and corneal tomography measurements were performed by an experienced operator. Corneal tomography was carried out based on the Scheimplfug technique using Pentacam HR (Oculus, Wetzlar, Germany). Moreover, the thinnest corneal thickness (TCT), steep keratometry (Ks), and flat keratometry (Kf) were collected. Considering the characteristic of higher Ks in KC, the eye with a higher Ks in each subject was included in the current analysis.
2.3 Mitochondrial DNA sequencing and bioinformatics analysis
Total DNA was extracted from peripheral blood samples with QIAamp DNA Blood kits (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The mitochondrial genome was amplified by long-range PCR using human mitochondrial genome specific primers. The detailed procedure of mtDNA sequencing has been described in our previous article (Xu et al., 2021). After sequencing on Illumina NovaSeq System (Illumina, San Diego, CA, United States), Fastq data were obtained and aligned to the mitochondrial reference sequence (revised Cambridge Reference Sequence, rCRS) by Burroughs-Wheeler Aligner (BWA) (Li and Durbin, 2009). The quality control of the sequencing data was conducted for each sample. The mean Q20 of the samples in R1 and R2 were 0.947 ± 0.047 and 0.969 ± 0.038 separately. The mean Q30 of the samples in R1 and R2 were 0.921 ± 0.051 and 0.843 ± 0.077 separately. The Genome Analysis Toolkit (GATK) (McKenna et al., 2010) was used to identify variants. Variants with a frequency range of 0.1–0.9 were considered heteroplasmic, while variants with a frequency of ≤0.1 or a frequency of ≥0.9 were considered homoplasmic. The ratio of the mutant allele was used to represent the heteroplasmic level. As previously recommended (Cosemans et al., 2023), only heteroplasmic variants with a prevalence of at least 10% in the study population were included for further analysis.
2.4 Statistical analyses
Statistical analyses were performed using SPSS 21.0 (SPSS Inc., Chicago, IL, United States). Quantitative variables were expressed using mean ± standard deviation (SD) and analyzed by the Student’s t-test. The distribution of sex was analyzed by the chi-squared test. The number of heteroplasmic and homoplasmic mtDNA variants was calculated in protein-coding genes (including ATP6, ATP8, COX1, COX2, COX3, CYTB, ND1, ND2, ND3, ND4, ND4L, ND5, and ND6), rRNA genes, tRNA genes, and non-coding region. A comparison of heteroplasmic and homoplasmic non-synonymous variants numbers in protein-coding genes was carried out by Fisher’s exact test or chi-squared test. The Spearman’s correlation was used to demonstrate the correlation between the number of heteroplasmic variants and age. Logistic regression analysis was adopted to explore the association of mtDNA heteroplasmic level with KC, and the β coefficient, stand error (SE), and odds ratio (OR) values were recorded. An adjusted significance level of p = 0.05/10 = 0.005 was applied. In addition, linear regression analysis was used to explore the relationship between mtDNA heteroplasmic levels and clinical parameters, and the β coefficient, SE, and t values were recorded.
3 Results
3.1 Clinical characteristics of study subjects
The clinical characteristics were compared between 300 KC cases and 300 controls. As is shown in Table 1, there were no significant differences in age and sex between KC cases and controls (p > 0.05). The TCT in KC cases was significantly lower than that in controls (p < 0.05). In contrast, the Ks and Kf in KC cases were significantly higher than those in controls (p < 0.05).
3.2 Distribution of heteroplasmic and homoplasmic mtDNA variants
The distribution of heteroplasmic and homoplasmic mtDNA variants across the mitochondrial genome was shown in Figure 1. Firstly, the heteroplasmic and homoplasmic variants were calculated in all the subjects (Supplementary Table S1). A total of 311 heteroplasmic variants and 1936 homoplasmic variants were identified. The non-coding region harbored the highest (125 out of 311) heteroplasmic variants among different mitochondrial regions. While in the protein-coding genes, the COX3 gene exhibited the highest heteroplasmic variants (Figure 1A). Then the heteroplasmic and homoplasmic mtDNA variants were calculated in KC cases and controls separately (Supplementary Tables S2, S3). In KC cases, nearly half of the heteroplasmic variants (92 out of 200) were observed in the non-coding region (Figure 1B). Similarly, 91 out of 199 variants were located in the non-coding region in the controls (Figure 1C). Among the protein-coding genes, the COX3 gene contained 13 heteroplasmic variants in KC cases, followed by 12 heteroplasmic variants in the ATP6 region. While in the control group, 13 heteroplasmic variants were observed in the COX3 gene and CYTB gene separately, followed by 10 heteroplasmic variants in the ND5 region.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1251951/fgene-14-1251951-HTML/image_m/fgene-14-1251951-g001.jpg)
FIGURE 1. The distribution of heteroplasmic and homoplasmic mtDNA variants across the mitochondrial genome. [(A): Distribution of heteroplasmic mtDNA variants in all subjects; (B): Distribution of heteroplasmic mtDNA variants in KC cases; (C): Distribution of heteroplasmic mtDNA variants in controls; (D): Distribution of homoplasmic mtDNA variants in all subjects; (E): Distribution of homoplasmic mtDNA variants in KC cases; (F): Distribution of homoplasmic mtDNA variants in controls].
In addition, the number of heteroplasmic and homoplasmic non-synonymous variants in protein-coding genes was compared between KC cases and controls (Figure 2; Supplementary Table S4). Finally, no statistically significant differences were found between the two groups. The correlation between the number of heteroplasmic variants and age was further analyzed in all subjects, KC cases, and controls separately (Figure 3). The number of heteroplasmic variants was significantly positively associated with age in all subjects (r = 0.085, p = 0.037), while the number of heteroplasmic variants was not associated with age in KC cases and controls separately (all p > 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1251951/fgene-14-1251951-HTML/image_m/fgene-14-1251951-g002.jpg)
FIGURE 2. Comparison of the number of heteroplasmic and homoplasmic nonsynonymous variants in protein-coding genes between KC cases and controls. [(A): ATP6; (B): ATP8; (C): COX1; (D): COX2; (E): COX3; (F): CYTB; (G): ND1; (H): ND2; (I): ND3; (J): ND4; (K): ND4L; (L): ND5; (M): ND6; (N): Total protein-coding genes].
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1251951/fgene-14-1251951-HTML/image_m/fgene-14-1251951-g003.jpg)
FIGURE 3. Correlation of mtDNA heteroplasmic variants number with age. [(A): All subjects; (B): KC cases; (C): Controls].
3.3 Association of mtDNA heteroplasmic level with KC
The association of mtDNA heteroplasmic level with KC was explored by logistic regression with adjusting age and sex. Ten heteroplasmic variants were included in the analysis. The prevalence of the ten heteroplasmic variants in our study population was shown in Supplementary Table S5. The heteroplasmic levels of variants in KC cases and controls were shown in Supplementary Table S6, and the associations of heteroplasmic levels with KC were shown in Table 2. The results showed that the heteroplasmic level of m.16180_16181delAA were associated with KC after Bonferroni correction (p = 0.004, OR = 0.001, 95%CI: 0–0.097).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1251951/fgene-14-1251951-HTML/image_m/fgene-14-1251951-t002.jpg)
TABLE 2. Logistic regression analysis of mitochondrial heteroplasmic level and KC with adjusting age and sex.
3.4 Correlation of mtDNA heteroplasmic level with clinical characteristics
Moreover, linear regression analyses were performed between m.302A>C,m.16180_16181delAA and clinical parameters including TCT, Ks, and Kf in KC cases and controls separately. As is shown in Table 3, the heteroplasmic levels of m.302A>C and m.16180_16181delAA were not correlated with TCT, Ks, and Kf (all p > 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1251951/fgene-14-1251951-HTML/image_m/fgene-14-1251951-t003.jpg)
TABLE 3. Liner regression analyses between heteroplasmic variants and TCT, Ks, Kf in KC cases and controls.
4 Discussion
The pathogenesis of KC is complex, and multiple studies have indicated an association of the mitochondrial genome with KC (Pathak et al., 2011; Abu-Amero et al., 2014b; Vallabh et al., 2017). In this study, we presented an extensive profile of mitochondrial heteroplasmy in KC. The results showed that the COX3 gene harbored the highest number of heteroplasmic variants among different protein-coding genes. In addition, a positive association was found between the number of heteroplasmic variants and age in all the subjects. Logistic regression analyses indicated that the heteroplasmic level of m.16180_16181delAA was associated with KC.
KC is a corneal disorder with complex etiology. Currently, growing evidence suggests that oxidative stress may play an important role in its pathogenesis (Navel et al., 2021; Monteiro de Barros and Chakravarti, 2022). Oxidative stress refers to the imbalance of two opposite and antagonistic forces, production of reactive oxygen species (ROS) and antioxidants (Sies, 2020). Actually, mitochondria are the most important source of ROS in most mammalian cells (Hernansanz-Agustin and Enriquez, 2021). Therefore, the mitochondrial dysfunction might lead to the imbalance of oxidative stress, finally resulting in the occurrence of diseases. Recently, multiple studies indicated an association of mitochondrial genome with KC, and several mtDNA variants have been identified in KC (Atilano et al., 2005; Abu-Amero et al., 2014a). As is reported previously, most of the pathogenic mtDNA variants are heteroplasmic, and mitochondrial heteroplasmy tends to show high pathogenicity (Ye et al., 2014; Nandakumar et al., 2021). Mitochondrial heteroplasmy is the co-existence of multiple mtDNA variants in a single source, and it is now generally accepted to play important roles in many diseases (Parakatselaki and Ladoukakis, 2021). However, the profile of mitochondrial heteroplasmy in KC has not been characterized up to now. The present study elucidated the genetic profile of mitochondrial heteroplasmy in KC. A comparison of heteroplasmic variants number among different regions showed that the heteroplasmic variants were mostly located in the non-coding region, consistent with Fendt et al. (2020) who reported that most heteroplasmic variants were located in the non-coding region in oral squamous cell carcinoma. While in the protein-coding genes, the COX3 gene exhibited the highest heteroplasmic variants. Since most mtDNA variants implicated in diseases are heteroplasmic, we inferred that the COX3 gene might have important roles in KC. Although the heteroplasmic levels of variants in COX3 gene exhibited no association with KC in the current study, our previous study identified an association between the COX3 gene and KC through gene-based SKAT analysis (Xu et al., 2021). Therefore, we speculated that each variant in COX3 gene might have a minor effect on KC, the combined effects of multiple variants resulted in the occurrence of the disease. Abu-Amero et al. (2014b) detected ten potentially pathogenic non-synonymous mtDNA mutations in KC, of which one non-synonymous variant was heteroplasmic, and the other nine were homoplasmic. Similarly, we also found a lower number of heteroplasmic non-synonymous variants in KC. However, there were no significant differences between heteroplasmic and homoplasmic non-synonymous variants. It has been reported that mtDNA heteroplasmy is associated with age (Nandakumar et al., 2021; Stewart and Chinnery, 2021). Therefore, we analyzed the correlation between the number of heteroplasmic variants and age. The results showed that the number of heteroplasmic variants was positively correlated with age in all subjects. However, there was no correlation when KC cases and controls were analyzed separately, which might be attributed to different phenotypes and sample sizes.
Different mitochondrial heteroplasmy levels have been implicated in many diseases. As is reported by Geng et al. (2019), the heteroplasmic levels in peripheral blood leukocytes were closely associated with clinical manifestations and valuable for evaluating the clinical severity of the m.3243A>G mutation. Besides, Atilano et al. (2021) found the m.13095T>C and m.13105A>G heteroplasmic levels were higher in age-related macular degeneration, with the higher heteroplasmic levels possibly representing potential biomarkers. In the present study, the association of heteroplasmic levels with KC was analyzed. The results showed that the heteroplasmic level of m.16180_16181delAA was associated with KC after adjusting the age and sex. In detail, the lower heteroplasmic level of m.16180_16181delAA owned a higher risk of KC. The m.16180_16181delAA was located in the non-coding region. Currently, multiple mtDNA variants in non-coding regions were reported to be associated with certain diseases. Deng et al. (2021) analyzed the correlation of mtDNA variants in the D-loop region with polycystic ovary syndrome, and found that the variants m.G207A, m.16036insGG, and m.16049insG were associated with decreased risk of the disease. Wu et al. (2018) identified several mtDNA variants in the D-loop region associated with Parkinson’s disease. However, the linkage of m.16180_16181delAA with other diseases hasn’t been reported. The lower TCT, higher Ks, and higher Kf were typical clinical characteristics of KC in corneal tomography (Mas Tur et al., 2017). However, the genetic correlation between those clinical parameters and the mitochondrial genome was not clear. Considering the association between the heteroplasmic levels of mtDNA variants and KC, we speculated that there might exist a correlation between the mitochondrial heteroplasmic levels and TCT, Ks, and Kf. The association analysis indicated that the heteroplasmic level of m.16180_16181delAA was associated with KC after Bonferroni correction (p = 0.004), and the heteroplasmic level of m.302A>C was associated with KC before correction (p = 0.009). Therefore, we explored the correlation between m.16180_16181delAA, m.302A>C heteroplasmic levels and TCT, Ks, Kf. Nevertheless, our findings showed that the heteroplasmic levels in m.16180_16181delAA and m.302A>C were not correlated with TCT, Ks and Kf in KC cases and controls, indicating that the differences of m.16180_16181delAA and m.302A>C heteroplasmic levels might not be attributed to TCT, Ks, and Kf. It is reported that the mitochondria dysfunction has been observed in KC. Given that the mitochondrial heteroplasmy played important roles in maintaining the mitochondrial function (Stefano et al., 2017), we inferred that the heteroplasmic levels of m.16180_16181delAA and m.302A>C might affect the risk of KC through influencing mitochondrial function. The mitochondrial dysfunction results in the change of oxidative stress, thus leading to the occurrence of KC.
Several limitations should be noted in the present study. Firstly, we only included peripheral blood samples to explore the mtDNA heteroplasmy in KC due to the difficulty of acquiring the cornea. Secondly, the mtDNA content, smoking status, and BMI were not adjusted in the analysis. Those influencing factors would be collected and analyzed in our following search. Thirdly, the mitochondrial function of heteroplasmic variants and the molecular basis of the relationship between mtDNA heteroplasmy and KC were not explored in the present study. Further studies will be conducted in the future.
In conclusion, the present study characterized the profile of mtDNA heteroplasmy in KC, and revealed that the heteroplasmic level of m.16180_16181delAA was associated with KC. The data implied that mitochondrial heteroplasmy might be involved in the pathogenesis of KC. Further research is required to better understand the complex interactions between mitochondrial heteroplasmy and KC.
Data availability statement
The datasets presented in this article are not readily available because data obtained from mitochondrial sequencing are sensitive. Requests to access the datasets can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by the Ethics Committee of Henan Eye Hospital [HNEECKY-2019 (5)]. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
SR conceived and designed this experiment. LX and KY collected the DNA samples and analyzed the data. QF and YG helped collect the patients’ information and analyzed the data. LX wrote the manuscript. KY and SR revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the Henan Natural Science Foundation (No. 222300420536), Special Program for Basic Research of Henan Eye Hospital (No. 20JCZD003), Henan Young Health Science and Technology Innovation Outstanding Program (No. YXKC2020023), Henan Provincial Building Key Program (No. SBGJ202002028), Henan Provincial Medical Science and Technology Joint Program (Nos. LHGJ20210080 and LHGJ20220707), Henan Provincial Scientific and Technological Research Project (Nos. 222102310307 and 222102310599), Youth Special Program for Basic Research of Henan Eye Hospital (Nos. 21JCQN006 and 21JCQN008).
Acknowledgments
We acknowledge all the subjects for providing blood samples.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1251951/full#supplementary-material
References
Abu-Amero, K. K., Al-Muammar, A. M., and Kondkar, A. A. (2014a). Genetics of keratoconus: where do we stand? J. Ophthalmol. 2014, 641708–641711. doi:10.1155/2014/641708
Abu-Amero, K. K., Azad, T. A., Kalantan, H., Sultan, T., and Al-Muammar, A. M. (2014b). Mitochondrial sequence changes in keratoconus patients. Invest. Ophthalmol. Vis. Sci. 55, 1706–1710. doi:10.1167/iovs.14-13938
Atilano, S. R., Coskun, P., Chwa, M., Jordan, N., Reddy, V., Le, K., et al. (2005). Accumulation of mitochondrial DNA damage in keratoconus corneas. Invest. Ophthalmol. Vis. Sci. 46, 1256–1263. doi:10.1167/iovs.04-1395
Atilano, S. R., Udar, N., Satalich, T. A., Udar, V., Chwa, M., and Kenney, M. C. (2021). Low frequency mitochondrial DNA heteroplasmy SNPs in blood, retina, and [RPE+choroid] of age-related macular degeneration subjects. PloS one 16, e0246114. doi:10.1371/journal.pone.0246114
Cosemans, C., Wang, C., Martens, D. S., Janssen, B. G., Vanpoucke, C., Lefebvre, W., et al. (2023). In utero exposure to air pollutants and mitochondrial heteroplasmy in neonates. Environ. Sci. Technol. 57, 350–359. doi:10.1021/acs.est.2c02556
Deng, X., Ji, D., Li, X., Xu, Y., Cao, Y., Zou, W., et al. (2021). Polymorphisms and haplotype of mitochondrial DNA D-loop region are associated with polycystic ovary syndrome in a Chinese population. Mitochondrion 57, 173–181. doi:10.1016/j.mito.2020.12.006
Fendt, L., Fazzini, F., Weissensteiner, H., Bruckmoser, E., Schonherr, S., Schafer, G., et al. (2020). Profiling of mitochondrial DNA heteroplasmy in a prospective oral squamous cell carcinoma study. Cancers (Basel) 12 (7), 1933. doi:10.3390/cancers12071933
Ferrari, G., and Rama, P. (2020). The keratoconus enigma: A review with emphasis on pathogenesis. Ocul. Surf. 18, 363–373. doi:10.1016/j.jtos.2020.03.006
Geng, X., Zhang, Y., Yan, J., Chu, C., Gao, F., Jiang, Z., et al. (2019). Mitochondrial DNA mutation m.3243A>G is associated with altered mitochondrial function in peripheral blood mononuclear cells, with heteroplasmy levels and with clinical phenotypes. Diabet. Med. 36, 776–783. doi:10.1111/dme.13874
Hahn, A., and Zuryn, S. (2019). Mitochondrial genome (mtDNA) mutations that generate reactive oxygen species. Antioxidants (Basel) 8 (9), 392. doi:10.3390/antiox8090392
Hao, X. D., Yao, Y. Z., Xu, K. G., Dong, B., Xu, W. H., and Zhang, J. J. (2022). Insufficient dose of ERCC8 protein caused by a frameshift mutation is associated with keratoconus with congenital cataracts. Invest. Ophthalmol. Vis. Sci. 63 (13), 1. doi:10.1167/iovs.63.13.1
Hernansanz-Agustin, P., and Enriquez, J. A. (2021). Generation of reactive oxygen species by mitochondria. Antioxidants (Basel) 10 (3), 415. doi:10.3390/antiox10030415
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324
Lucas, S. E. M., and Burdon, K. P. (2020). Genetic and environmental risk factors for keratoconus. Annu. Rev. Vis. Sci. 6, 25–46. doi:10.1146/annurev-vision-121219-081723
Mas Tur, V., MacGregor, C., Jayaswal, R., O'Brart, D., and Maycock, N. (2017). A review of keratoconus: diagnosis, pathophysiology, and genetics. Surv. Ophthalmol. 62, 770–783. doi:10.1016/j.survophthal.2017.06.009
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Monteiro de Barros, M. R., and Chakravarti, S. (2022). Pathogenesis of keratoconus: NRF2-antioxidant, extracellular matrix and cellular dysfunctions. Exp. Eye Res. 219, 109062. doi:10.1016/j.exer.2022.109062
Nandakumar, P., Tian, C., O'Connell, J., andMe Research, T., Hinds, D., Paterson, A. D., et al. (2021). Nuclear genome-wide associations with mitochondrial heteroplasmy. Sci. Adv. 7 (12), eabe7520. doi:10.1126/sciadv.abe7520
Navel, V., Malecaze, J., Pereira, B., Baker, J. S., Malecaze, F., Sapin, V., et al. (2021). Oxidative and antioxidative stress markers in keratoconus: a systematic review and meta-analysis. Acta Ophthalmol. 99, e777–e794. doi:10.1111/aos.14714
Parakatselaki, M. E., and Ladoukakis, E. D. (2021). mtDNA heteroplasmy: origin, detection, significance, and evolutionary consequences. Life (Basel) 11 (7), 633. doi:10.3390/life11070633
Pathak, D., Nayak, B., Singh, M., Sharma, N., Tandon, R., Sinha, R., et al. (2011). Mitochondrial complex 1 gene analysis in keratoconus. Mol. Vis. 17, 1514–1525. doi:10.1016/j.fsigss.2011.09.012
Rabinowitz, Y. S. (1998). Keratoconus. Surv. Ophthalmol. 42, 297–319. doi:10.1016/s0039-6257(97)00119-7
Roy, S., Yadav, S., Dasgupta, T., Chawla, S., Tandon, R., and Ghosh, S. (2019). Interplay between hereditary and environmental factors to establish an in vitro disease model of keratoconus. Drug Discov. Today 24, 403–416. doi:10.1016/j.drudis.2018.10.017
Santodomingo-Rubido, J., Carracedo, G., Suzaki, A., Villa-Collar, C., Vincent, S. J., and Wolffsohn, J. S. (2022). Keratoconus: an updated review. Cont. Lens Anterior Eye 45, 101559. doi:10.1016/j.clae.2021.101559
Sarezky, D., Orlin, S. E., Pan, W., and VanderBeek, B. L. (2017). Trends in corneal transplantation in keratoconus. Cornea 36, 131–137. doi:10.1097/ICO.0000000000001083
Sies, H. (2020). Oxidative stress: concept and some practical aspects. Antioxidants (Basel) 9 (9), 852. doi:10.3390/antiox9090852
Stefano, G. B., Bjenning, C., Wang, F., Wang, N., and Kream, R. M. (2017). Mitochondrial heteroplasmy. Adv. Exp. Med. Biol. 982, 577–594. doi:10.1007/978-3-319-55330-6_30
Stewart, J. B., and Chinnery, P. F. (2021). Extreme heterogeneity of human mitochondrial DNA from organelles to populations. Nat. Rev. Genet. 22, 106–118. doi:10.1038/s41576-020-00284-x
Vallabh, N. A., Romano, V., and Willoughby, C. E. (2017). Mitochondrial dysfunction and oxidative stress in corneal disease. Mitochondrion 36, 103–113. doi:10.1016/j.mito.2017.05.009
Wallace, D. C., and Chalkia, D. (2013). Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease. Cold Spring Harb. Perspect. Biol. 5, a021220. doi:10.1101/cshperspect.a021220
Wu, H. M., Li, T., Wang, Z. F., Huang, S. S., Shao, Z. Q., Wang, K., et al. (2018). Mitochondrial DNA variants modulate genetic susceptibility to Parkinson's disease in Han Chinese. Neurobiol. Dis. 114, 17–23. doi:10.1016/j.nbd.2018.02.015
Xu, L., Yang, K., Fan, Q., Zhao, D., Pang, C., and Ren, S. (2021). Whole mitochondrial genome analysis in Chinese patients with keratoconus. Mol. Vis. 27, 270–282.
Yang, K., Xu, L., Fan, Q., Gu, Y., Zhang, B., Meng, F., et al. (2021). A hospital-based study on clinical data, demographic data and visual function of keratoconus patients in Central China. Sci. Rep. 11, 7559. doi:10.1038/s41598-021-87291-y
Keywords: keratoconus, mitochondrial DNA, heteroplasmy, variant, m.16180_16181delAA
Citation: Xu L, Yang K, Fan Q, Gu Y and Ren S (2023) Mitochondrial DNA heteroplasmy analysis in keratoconus patients from China. Front. Genet. 14:1251951. doi: 10.3389/fgene.2023.1251951
Received: 03 July 2023; Accepted: 05 September 2023;
Published: 18 September 2023.
Edited by:
Jared C. Roach, Institute for Systems Biology (ISB), United StatesReviewed by:
Hao Liu, China Agricultural University, ChinaCharlotte Cosemans, University of Hasselt, Belgium
Siddhesh Aras, Wayne State University, United States
Copyright © 2023 Xu, Yang, Fan, Gu and Ren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shengwei Ren, ysgzz2018@163.com