Abstract
To identify the role of mitochondrial DNA (mtDNA) mutations in uterine fibroids patients, genomic DNA isolated from paired myometrium and fibroid tissues was screened for mutations. The present study represents the first investigation to report that 10.4% of uterine fibroids cases had either mtDNA mutations or polymorphisms or both. Among the 14 mitochondrial sequence variants identified, seven are somatic mutations (A3327C, G3352A, G3376A, G3380A, G3421A, T15312G, and C15493G) and the remaining (G3316A, C3342A, C3442T, T10205A, A10188G, A10229C, and A10301T) are gene polymorphisms. Somatic mutations were both homo- and heteroplasmic in nature. Of the seven somatic mutations located in the MTND1 and MTCYB genes, five (71.42%) are nonsynonymous in nature, whereas four (57.14%) of the polymorphisms located in MTND1 and MTND3 genes are found to be nonsynonymous. Sequence variants such as G3380A, G3421A, T15312G, G3376A, and G3316A have been earlier described in different human pathologies, but the remaining are novel ones. Mitochondrial somatic mutations and polymorphisms may predispose women to an earlier onset of degenerative cellular processes, which impair oxidative phosphorylation capacity and thereby promote tumorigenesis in uterine smooth muscle cells. Detection of mtDNA sequence variations in fibroid patients raises the need for larger case-control studies to screen the whole mitochondrial genome and evaluate as a future diagnostic biomarker in fibroid patients.
Introduction
It is apparent that estrogen acts as a promoter of fibroid growth and mitochondria constitute an important target of action (Felty and Roy, 2005). Estrogen has been reported to alter the reactive oxygen species (ROS) levels inside the mitochondrial matrix region (Liehr and Roy, 1998; Felty et al., 2005). Mitochondrial DNA (mtDNA) is more vulnerable to oxidative damage and undergoes a higher rate of mutation than does the nuclear genome (Jeon et al., 2005); hence, it is conceivable that estrogen-stimulated ROS levels could lead to mitochondrial genetic instability in tumor cells. Several studies have identified the existence of mtDNA alterations, such as missense mutations, frame shift mutations, and multiple gene insertions/deletions in breast, ovarian, cervical, and uterine endometrial neoplasms and cell lines (Copeland et al., 2002; Chatterjee et al., 2006; Czarnecka et al., 2006). Mutations in mtDNA could in turn elevate ROS production, leading to additional mutations and oxidative stress; this could provide replication advantage to transformed cells to evade apoptosis and drive cellular proliferation (Rustin, 2002). Given the critical role of estrogen in the tumorigenesis of fibroids and in mitochondria, the role of mtDNA mutations could be hypothesized in the pathogenesis of uterine fibroids. Therefore, mutation screening was done at different gene regions of mtDNA isolated from fibroid tumors and matched myometrium tissues.
Materials and Methods
Enrollment of subjects and clinical specimen collection
Sixty-nine uterine fibroid patients from government maternity hospitals in Hyderabad were enrolled in this study. Case diagnosis was based on the gynecological and ultrasonography examinations. Clinical details including age at menarche, associated symptoms, and length and duration of menstrual cycle were recorded in case sheets. A total of 138 tissue samples, comprised of uterine fibroids and adjacent normal myometrium of the same patient, were collected after vaginal or abdominal hysterectomy. Histological characterization of tissues was performed by an experienced pathologist. The present study protocol was approved by the hospital ethics committee and informed consent was obtained from all the patients.
DNA isolation, polymerase chain reaction, and mutation detection
Genomic DNA was extracted from tissue samples by following the method described by Alluri et al. (2005). Individual polymerase chain reactions (PCRs) were performed to amplify the 2993-3493, 5909-6301, 7471-8566, 8211-8411, 10076-10357, and 14975-15744 gene regions using the sequence-specific primers (Table 1). Specific conditions for PCR including the annealing temperature and its duration, in addition to the number of cycles, are mentioned in Table 2. Single-stranded conformational polymorphism (SSCP) and restriction enzyme digestion methods were followed to screen the mutations, based on the size of PCR amplicons (Dong and Zhu, 2005). SSCP with PCR amplicons of the 2993-3493 region was carried out after initial digestion with HaeII enzyme. PCR amplicons that showed mobility shift on SSCP analysis or altered restriction pattern by restriction enzyme digestion analysis were initially purified using the QIAquick (Qiagen) gel extraction kit, followed by a bidirectional sequencing reaction using an ABI PRISM 310 Genetic Analyzer (Perkin-Elmer).
Analysis of the mtDNA sequences
The results of the mtDNA sequence analysis were compared with the published NCBI sequences, using the BLAST search tool for mutation identification (www.ncbi.com/blast). Further, sequence variants were compared against mitochondrial databanks, such as the MITOMAP and mtDBase banks to identify whether those variants were reported earlier (www.mitomap.org/cgibin/tbl14gen.pl and www.genpat.uu.se/mtDB/index.html). Somatic mtDNA mutations are defined as those present in the fibroids but not in the myometrium tissues samples from the patient concerned, whereas sequence variants, which are present in both the fibroids and matched myometrium tissues, were scored as germline polymorphisms. Somatic mutations or polymorphisms that are not found in MITOMAP database are considered as novel.
Results
Clinical analysis
Of the cases included in the study, 19 were diagnosed to have single and 50 to have multiple fibroids, which varied in size between 1.0×1.0 and 9.0×6.0 cm, with a mean area of 20.0 cm. The mean age at fibroids detection was 40±1.02 years and the age at menarche was 13.8±0.72 years. Based on three-generation pedigree analysis, 9 of 69 cases (13.04%) reported positive family history of fibroids in their first-degree relatives. Nulliparity was reported in 3 of 69 (4.34%) fibroid cases.
Somatic mutations and polymorphisms identification
On sequencing analysis of PCR amplicons, which showed altered mobility on SSCP gels (six cases) or altered banding pattern on agarose gels (one case), a total of 14 mutations were detected in seven fibroid cases (10.14%). Among these, seven are mitochondrial somatic variations (A3327C, G3352A, G3376A, G3380A, G3421A, T15312G, and C15493G) and the remaining seven (G3316A, C3342A, C3442T, T10205A, A10188G, A10229C, and A10301T) are germline polymorphisms (Table 2). Except T15312G and C15493G, all other sequence variations were seen in homoplasmic form. Of the 14 different sequence variants identified, 8 mutations were located in the MTND1 gene and the remaining 6 mutations were in the MTND3 and MTCYB genes. Of the seven somatic mutations of the MTND1 and MTCYB genes, five (71.42%) are nonsynonymous, that is, resulting in an amino acid change, and the remaining two (28.57%) are synonymous, whereas among the seven gene polymorphisms, four (57.14%) were found to be nonsynonymous and three (42.85%) are synonymous. Nine of the 14 sequence variants, A3327C, G3352A, C15493G, A10188G, T10204A, A10229C, A10301T, C3342A, and C3442T, were not recorded in the mtDBase database and MITOMAP databases; hence, these are novel mtDNA sequence variations.
Discussion
The application of mDNA mutation and/or polymorphism patterns as a biomarker is rapidly expanding in disciplines ranging from rare metabolic diseases and aging to tumorigenesis. Mitochondrial mutations are associated with breast (15.6%-93%), ovarian (26%), and cervical (27%-95%) neoplasms, which indicates the rising significance of mitochondria in estrogen-sensitive pathologies (Bianchi et al., 1995; Parrella et al., 2001; Tan et al., 2002; Zhu et al., 2004; Sharma et al., 2005; Tseng et al., 2006; Wang et al., 2007; Chen and Zhan, 2009; Shen et al., 2010). The present study represents the first investigation reporting that 10.4% of uterine fibroids cases had either mtDNA mutations or polymorphisms or both.
Upon comparison of all the identified somatic mitochondrial mutations in mitochondrial genome databases, G3380A was previously found to be detected in colon tumors (Habano et al., 1999). Interestingly, this same mutation was described in a heteroplasmic fashion in muscle fibers of mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) patients as well (Horvath et al., 2008). On the other hand, the G3421A mutation was reported in pancreatic cancer patients (Jones et al., 2001) and also in a family with maternally inherited diabetes and deafness but in a homoplasmic fashion (Chen et al., 2006). The G3376A mutation has been implicated as pathogenic for Leber's hereditary optic neuropathy (LHON) and MELAS patients (Blakely et al., 2005). The nonsynonymous nature of all these transitions could affect the functionality of the highly conserved region of mitochondrial NADH dehydrogenase subunit 1 (MTND1) protein. Another nonsynonymous mutation, T15312G, located in mitochondrial cytochrome b (MTCYB), was previously reported in thyroid tumors (Chatterjee et al., 2006). The higher percentage of homoplasmic somatic mutations or polymorphisms observed in fibroid cases could have been driven by random genetic drift mechanism, by which preexisting heteroplasmic sequence variants were outgrown in successive proliferation of fibroid progenitor cells (Coller et al., 2001; Jones et al., 2001). On the basis of these findings, it is suggested that the observed somatic mutations are pathogenic and are functionally selected during the neoplastic initiation and growth. None of the sequence variations showed specific clinopathological correlation with tumor types, nulliparity, and family history. Differential occurrence of mtDNA mutations in various tumors is partly due to the prominence of different repair mechanisms. A majority of fibroid cases are asymptomatic and go undetected and this could be attributed to the low frequency of family history among the mutant cases.
Among the gene polymorphisms identified, G3316A, which is seen in one fibroid case, was previously reported in patients with gestational diabetes and LHON (Odawara et al., 1996; Yu et al., 2004; Phasukkijwatana et al., 2006; Tang et al., 2006). Interestingly, the C3342A transversion, seen in the same patient, was earlier reported as C to T in a healthy Caucasian population (Snejina et al., 2006). Polymorphisms, such as G3316A, T10204A, A10301T, A10188G, and A10229C, can cause subtle differences in the encoded protein structure and function. No functional significance is attributed to sequence variations such as A3327C, C3342A, C3442T, A10301T, and C15493G, as they are synonymous. Gene polymorphisms of mtDNA predispose individuals to an earlier onset of degenerative cellular processes and decline in OXPHOS capacity (Trounce and Pinkert, 2007). These finding indicates that mitochondrial polymorphisms are indeed frequent in fibroid patients and raises the need to conduct larger case-control studies to investigate them as risk markers to fibroid susceptibility.
MTND1, MTND3, and MTCYB genes encode integral polypeptide subunits of mitochondrial respiratory chain complexes I, II, and III, respectively. Nonsynonymous sequence variations in these genes could alter the sequence in the encoded proteins, which in turn affects the stability and functioning of mitochondrial respiration, resulting in a change in mitochondrial ROS levels (Ugalde et al., 2003; Hinttala et al., 2006; Malfatti et al., 2007). Elevated ROS levels damage mtDNA and also result in membrane fluidity loss and defective apoptotic effects. The mechanisms of the differences in phenotype associated with mtDNA mutations and polymorphisms are not fully understood, but differences in organ-specific thresholds for energy demand and/or the presence of additional nuclear gene mutations may also modulate mitochondrial pathologies. In view of the high energy demand that exists in uterine smooth muscle tissues, disruption in cellular energy production could promote fibroid tumorigenesis.
Although none of the sequence variations are observed in the MTCO1, MTCO2, MTRNR2, MTTL1, and MTATP8 gene regions, the possibility of mutations at other sequences cannot be ruled out. Hence, comprehensive screening of the whole mtDNA genome in fibroid tumors could identify more mitochondrial sequence variants. We were unable to observe any trend between mtDNA mutations and clinical and pathological characteristics, or family history, possibly because of our relatively small sample size and detection bias.
In conclusion, the present study demonstrates the existence of mtDNA somatic mutations and polymorphisms in uterine fibroid patients. Additional evaluation of the molecular consequences of somatic mutations in fibroids may elucidate their role in energy metabolism of fibroid tumors. In a clinical context, application of PCR-coupled restriction enzyme and/or SSCP methods for rapid multisampling analysis of fibroid patients may potentiate the detection of mtDNA sequence variants with minimum cost.
Footnotes
Acknowledgments
The authors acknowledge DBT, Government of India (BT/PR/5516/MED/14/646/2004), for their financial support. The authors greatly thank the management and staff of Bhagwan Mahaveer Medical Research Centre, Hyderabad, Govt. Maternity Hospitals (Nayapool and Kothi), Hyderabad, and Bhagwan Mahaveer Medical Research Centre, Hyderabad, for providing clinical specimens and data. The authors sincerely thank Miss. Maria Arland for helping in sample collection and S. Movva for technical help.
Disclosure Statement
No competing financial interests exist.
