Abstract
Because energy production involves oxidative phosphorylation, mitochondria are major sources of reactive oxygen species in the cell. Recent findings indicate that mitochondrial DNA (mtDNA) variants may play a role in the etiology of certain autoimmune and chronic inflammatory diseases. The aim of this study was to investigate the possible association between mtDNA polymorphisms and susceptibility to endometriosis. This study included 198 patients with histologically confirmed endometriosis and 167 patients without endometriosis as controls. Common variants of mtDNA at nt10398 (A/G transition), nt13708 (G/A transition), and nt16189 (T/C transition) were detected using polymerase chain reaction. An association study was performed with a chi-square test and logistic regression analysis. The prevalence of the mtDNA nt16189 variant was higher in patients with endometriosis (46.0%, 91 of 198) than in controls (34.7%, 58 of 167) (p=0.030) with odds ratio (OR) of 1.98 (95% confidence interval [CI]: 1.04–3.78). A combination of the 10398 and 16189 variants was also associated with increased risk for endometriosis (OR=1.90, 95% CI: 1.13–3.18, p=0.015). These associations remained significant even after adjusting for age and body mass index. Our data strongly suggest that the mtDNA 16189 variants and the combination of mtDNA 16189 and 10398 variants increase susceptibility to endometriosis.
Introduction
Mitochondria are organelles found in all nucleated cells. They are composed of about 1500 proteins encoded by nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) (Taylor and Turnbull, 2005). The mitochondrial genome (mtDNA) consists of closed circular DNA encoding 2 rRNA, 22 tRNA, and 13 proteins (Anderson et al., 1981). All 13 proteins encoded by mtDNA are involved in mitochondrial oxidative phosphorylation (OXPHOS) enzyme complexes that are essential for ATP production. Because energy production involves OXPHOS, mitochondria are also major sources of ROS in the cell. It is estimated that 90% of cellular oxygen is consumed in the mitochondria and that about 2%–4% of that oxygen is converted to ROS. As there is no system that efficiently protects and repairs mtDNA under much higher levels of oxidative stress, mtDNA is more vulnerable to alteration than nDNA (Kang and Hamasaki, 2005).
Although defective mitochondrial function is most commonly thought to be related to neuromuscular and neurodegenerative diseases, studies indicate that mtDNA mutations are also involved in other human diseases and malignancies (Carew and Huang, 2002; Mootha et al., 2003; Suzuki et al., 2003; Rohan et al., 2010). A number of recent studies have suggested that mtDNA variants may also play a role in the etiology of certain cancers, autoimmune diseases, and chronic inflammatory diseases (Canter et al., 2005; Yu et al., 2008, 2009; Jonsen et al., 2009; Singh and Kulawiec, 2009). However, the data concerning the association of defective mitochondrial function with endometriosis are limited (Kao et al., 2005). Given that mitochondria are involved in ROS formation, apoptosis, and energy production required for the activation and proliferation of peripheral lymphocytes and that mitochondria and autoimmunity are closely linked to each other, we hypothesized that mtDNA variants are involved in the pathogenesis of endometriosis. In this study, we investigated three common mtDNA variants at nt10398 (A→G transition), nt13708 (G→A transition), and nt16189 (T→C transition), which have been shown to be involved with ROS production and oxidative damage, and assessed whether these mtDNA polymorphisms are associated with susceptibility to endometriosis.
Materials and Methods
Study subjects
The present study included 365 premenopausal women who underwent elective surgery with subsequent pathological evaluation at the Department of Obstetrics and Gynecology, Gangnam Severance Hospital, between January 2007 and March 2010. The study sample was biologically homogeneous, as all subjects were of Korean descent. This study was approved by the Institutional Review Board of Gangnam Severance Hospital and all subjects provided written informed consent. One hundred ninety-eight women had direct visualization of endometriotic lesions and/or pathologically confirmed endometriosis at the time of surgery, and the stage of endometriosis was determined using the revised American Society of Reproductive Medicine classification (American Society for Reproductive Medicine, 1997). Thirty-seven patients were diagnosed with minimal to mild disease (26 patients with stage I and 11 patients with stage II endometriosis) and 161 patients were diagnosed with moderate to severe disease (80 patients with stage III and 81 patients with stage IV endometriosis). One hundred sixty-seven patients without evidence of endometriosis served as control subjects. Controls included 55 cases of dermoid cysts, 40 cases of serous cystadenoma, 28 cases of mucinous cystadenoma, 38 cases of myomas, and 6 cases of paratubal cysts that were otherwise healthy. Postmenopausal women, previous hormone or GnRH agonist users, and patients who had adenomyosis, endometrial cancer, endometrial hyperplasia or endometrial polyps, tubo-ovarian abscess, and/or pelvic inflammatory disease were excluded from this study.
For all patients, pretreatment serum CA-125 levels were measured with CA-125 II electrochemiluminescence immunoassay using a Roche/Hitachi Modular Analytics E170 (Roche Diagnostics, Tokyo, Japan). The presence of dysmenorrhea before surgery was recorded using a four-point scale and a Visual Analog Scale (VAS). The scores for the four-point scale were as follows: 0=no dysmenorrhea; 1=minimal (can work, somewhat uncomfortable); 2=moderate (can work, but quite uncomfortable); or 3=severe (miss work, has to be in bed). VAS consisted of a 10-cm-linear analog scale marked from 0 to 10, where 0 represented no pain and 10 represented the most severe pain. The score was recorded by marking a point on the 10-cm line. Body mass index (BMI) was calculated as weight in kilograms divided by height in meters squared.
Sample DNA extraction and mtDNA genotyping
Blood samples were collected from all patients using a standard venipuncture technique and ethylenediaminetetraacetic acid tubes. DNA was subsequently extracted from peripheral blood lymphocytes using standard methods.
The three single-nucleotide polymorphisms (SNPs) of mtDNA nt10398 (A/G), nt13708 (G/A), and nt16189 (T/C) were amplified by three polymerase chain reactions (PCRs). The primers were generally positioned approximately 50–60 bp from the boundary to allow for the detection of the SNPs (Table 1). PCR using genomic DNA was carried out in 30 μL reaction volumes containing 100 ng genomic DNA, 0.2 mM primers, 100 mM dNTP, reaction buffer (100 mM Tris [pH 8.3], 500 mM KCl, 15 mM MgCl2, 0.01% gelatin), and 5 U/μL SP-Taq polymerase (COSMO Genetech. Co. Ltd., Seoul, Korea) with the following cycling profile: denaturation at 95°C for 5 min followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 60°C for 30 s, and extension at 72°C for 30 s and a final extension step at 72°C for 30 s. All thermal cycles were run on TaKaRa PCR Thermal Cycler (TaKaRa). All PCR fragments were purified using a PCR purification kit (COSMO Genetech. Co. Ltd.) after gel electrophoresis with a 1.5% agarose gel. All PCR products were directly sequenced in both directions using a 3100 Genetic Analyzer (PE Applied Biosystems, Foster City, CA) with PRISM Dye Terminator and Dye Primer Cycle Sequencing chemistries.
F, forward; R, reverse.
Statistical analysis
Statistical analysis was conducted using SPSS version 18 (SPSS, Chicago, IL). Continuous variables were expressed as means±standard deviation and compared using Student's t-test. The statistical difference in the frequency of occurrence of the variants of mtDNA between endometriosis and control group was assessed by Pearson's χ 2 test. Logistic regression analysis was used to analyze the association of variants of mtDNA with endometriosis after correction for age and BMI. p<0.05 was considered significant.
Results
The clinical characteristics of the patients are shown in Table 2. Mean age, mean parity, and BMI were significantly lower in patients with endometriosis than in controls, whereas scores from the four-point scale and VAS for dysmenorrhea were significantly higher in the endometriosis group than in the control group. The mean serum CA-125 levels were 74.87±180.21 IU/mL for endometriosis group and 32.45±82.79 IU/mL for the control group (p=0.005).
Data are expressed as mean±standard deviation.
BMI, body mass index; VAS, Visual Analog Scale.
Among the three mtDNA variants, only the prevalence of nt16189 was significantly higher in patients with endometriosis than in controls (45.9% vs. 34.7%, p=0.030) with odds ratio (OR) of 1.598 (95% confidence interval [CI]: 1.046–2.442) (Table 3). The other two mtDNA variants were not associated with susceptibility to endometriosis in the cohort. We further analyzed the frequency of combination of the three mtDNA variants, and the prevalence of the combination of nt10398 and nt16189 variants was significantly higher in the endometriosis group than in the control group (26.8% vs. 16.2%, p=0.015, OR: 1.895, 95% CI: 1.129–3.183).
Data are expressed as n (%).
SNP, single-nucleotide polymorphism; OR, odds ratio; CI, confidence interval.
For further analysis, the endometriosis group was divided into subgroups according to severity of disease (Table 4). When the prevalence of three mtDNA variants and combination of these three mtDNA variants in patients with endometriosis were compared with controls according to severity of disease, the differences in the prevalence of nt16189 variant and the combination of nt10398 and nt16189 variants were more prominent in patients with minimal to mild disease, whereas only a marginal significance was noted in patients with moderate to severe disease. However, no significant differences in the frequency of mtDNA variants within each endometriosis subgroup were noted. We also conducted a subgroup analysis according to the degree of dysmenorrhea, but found no significant differences in the frequency of mtDNA polymorphisms between each of the endometriosis subgroups and the control group (data not shown).
Versus control.
Data are expressed as n (%).
The accumulation of mtDNA mutations has been linked to the aging process and the association of BMI with endometriosis has been previously documented (Hediger et al., 2005; Sato et al., 2007; Yi et al., 2009; Stefanatos and Sanz, 2011). As there were significant differences in age and BMI between study groups, logistic regression analysis was performed. A multivariable logistic regression analysis showed that mtDNA nt16189 variant and the combination of mtDNA nt10398 and nt16189 variants were significantly associated with susceptibility to endometriosis even after adjusting for age and BMI (Table 5).
Multivariate regression included age and BMI.
Discussion
In this study, we demonstrated that mtDNA nt16189 variant and the combination of mtDNA nt10398 and nt16189 variants are associated with an increased risk for endometriosis. To the best of our knowledge, this is the first study to date that has examined the association of mtDNA polymorphisms with endometriosis, and our results imply that mtDNA polymorphism may play an important role in the pathogenesis of endometriosis.
Although mitochondria perform multiple cellular functions including control of cell death, growth, and development and integration of signals from mitochondria to nucleus and nucleus to mitochondria, they are the major sources of energy and ROS production (Singh and Kulawiec, 2009). Although mitochondria are major sources of ROS in the cell, it is true that mitochondria are not the only major source of ROS production. However, as the mitochondrial genome is close to the site of ROS production, lacks histones and introns, and has much less efficient DNA repair mechanisms than nDNA, mtDNA is particularly susceptible to damage by ROS (Canter et al., 2005). In addition to somatic changes, mtDNA polymorphisms may have subtle effects on ROS production. However, it has been postulated that if the variants reduce efficacy of mitochondrial functioning, the accumulation of ROS may lead to cancer (Bai et al., 2007). Given the various roles of mitochondria in cellular functions, it is biologically plausible that mtDNA mutations affecting mitochondrial respiratory enzyme complexes may be involved in the pathogenesis of endometriosis by accumulation of ROS and dysfunctional mitochondria induced apoptosis.
In this study, we demonstrated that mtDNA nt16189 variant is associated with endometriosis. The mtDNA nt16189 variant has been shown to be associated with insulin resistance, metabolic syndrome, and susceptibility to type 2 diabetes (Poulton et al., 2002; Weng et al., 2005; Bhat et al., 2007). It is located in a cytosine run within the noncoding D-loop of mtDNA, close to sequences that are conserved across multiple species, which may be involved in the regulation of replication and transcription (Sbisa et al., 1997). Previous studies indicate that mtDNA nt16189 variant may impair the ability of a cell to properly respond to oxidative stress and oxidative damage (Lin et al., 2005a, 2005b). Because oxidative stress plays an important role in both immune response and inflammation, this finding may support a mechanism for the association of the nt16189 variant with susceptibility to endometriosis.
We also evaluated the association of the mtDNA nt10398 variant with susceptibility to endometriosis. The mtDNA nt10398 variation results in a nonconserved amino acid substitution of threonine (encoded by the A allele) with alanine (encoded by the G allele) within the NADH dehydrogenase 3 (ND3) subunit of Complex I. Complex I serves to dehydrogenate NADH and shuttle electrons to coenzyme Q, generating superoxide as a byproduct (Lenaz et al., 2006). There is also an increasing support for the association of the 10398 variation in the ND3 locus of mtDNA with increased rates of electron leakage and ROS production due to altered Complex I function (Ross et al., 2001; van der Walt et al., 2003). This polymorphism has been shown to be associated with type 2 diabetes, breast cancer, and secondary antiphospholipid syndrome in systemic lupus erythematosus patients (Canter et al., 2005; Bhat et al., 2007; Jonsen et al., 2009). Mechanisms of the adverse effects of the 10398 A allele in other diseases involving oxidative stress include less efficient functioning of Complex I in individuals whose mitochondria carry this allele and, alternatively, the reduced ability of cells in these individuals to constrain an oxidative challenge (van der Walt et al., 2003; Canter et al., 2005). Although our results showed no significant association between the nt10398 variant and endometriosis, the combination of the nt10398 and nt16189 variants displayed significant association with endometriosis. This finding may be explained by the fact that these two variants are closely related to ROS production and ROS clearance, which seem to be a pivotal event in the development of endometriosis. We assume that the accumulation of excessive oxidative stress caused by overproduction of ROS and the alteration of Complex I due to mitochondrial 10398 variation is further exacerbated by the diminished ability of the cell to properly respond to oxidative stress and oxidative damage due to mtDNA nt16189 variant. Thus, the combination of these variants may lead to susceptibility to endometriosis.
One interesting finding in this study is that these associations were more prominent in patients with minimal to mild disease than in patients with moderate to severe disease. Other studies have also associated genetic polymorphisms with minimal to mild endometriosis, suggesting that the genotypes associated with susceptibility to endometriosis may differ according to severity of disease (Baranova et al., 1999; Kim et al., 2005). The sample size in this study, however, is too small to draw such conclusions and therefore warrants a follow-up study using a larger number of patients in this specific disease population.
Limitations to our study include a relatively small sample size. Salanti et al. (2005) suggested that the most realistic OR between SNPs and a complex trait disease can be achieved with a sample size of at least 1000 subjects. Given that our sample size is smaller than the ideal sample size, meta-analysis with multiple studies is needed. Second, healthy women without pelvic pain or disease would have been an ideal comparison group to distinguish between endometriosis and normal patients. Choosing adequate control patients is a complex and often overlooked problem in endometriosis research. Although we only included patients surgically proven to be free of endometriosis, the majority of our control group had other benign diseases, which may have different effects on mtDNA polymorphisms. Further, correction for multiple comparisons has not been applied to the results in this hypothesis-generating study, increasing the possibility for chance findings.
In conclusion, this study demonstrated that mtDNA 16189 variant and the combination of mtDNA 16189 and 10398 variants are associated with increased risk of endometriosis. Although the mechanism underlying this association remains unclear, these findings can help to define the role of the mitochondrial genome in the pathogenesis of endometriosis. However, this is a preliminary study showing associations with several limitations. Therefore, follow-up studies with a large number of patients are warranted to confirm the possible association of mtDNA polymorphisms with susceptibility to endometriosis and to define their exact role in the pathogenesis of endometriosis.
Footnotes
Acknowledgments
This research was supported by Faculty Research Grant from Yonsei University, College of Medicine, Seoul, Korea (6-2011-0093).
Disclosure Statement
No potential conflicts of interest exist.
