Original research | Open | Published:
Common mitochondrial polymorphisms as risk factor for endometrial cancer
International Archives of Medicinevolume 2, Article number: 33 (2009)
Endometrial carcinoma is the most commonly diagnosed gynaecological cancer in developed countries. Although the molecular genetics of this disease has been in the focus of many research laboratories for the last 20 years, relevant prognostic and diagnostic markers are still missing. At the same time mitochondrial DNA mutations have been reported in many types of cancer during the last two decades. It is therefore very likely that the mitochondrial genotype is one of the cancer susceptibility factors. To investigate the presence of mtDNA somatic mutations and distribution of inherited polymorphisms in endometrial adenocarcinoma patients we analyzed the D-loop sequence of cancer samples and their corresponding normal tissues and moreover performed mitochondrial haplogroup analysis. We detected 2 somatic mutation and increased incidence of mtDNA polymorphisms, in particular 16223C (80% patients, p = 0.005), 16126C (23%, p = 0.025) and 207A (19%, p = 0.027). Subsequent statistical analysis revealed that endometrial carcinoma population haplogroup distribution differs from the Polish population and that haplogroup H (with its defining polymorphism - C7028T) is strongly underrepresented (p = 0.003), therefore might be a cancer-protective factor. Our report supports the notion that mtDNA polymorphisms establish a specific genetic background for endometrial adenocarcinoma development and that mtDNA analysis may result in the development of new molecular tool for cancer detection.
Endometrial cancer (EC) is the most frequently occurring invasive neoplasm of the female genital tract worldwide [1, 2]. In 2007 approximately 39,000 new cases presented in the United States and 149,300 in Europe making it the fourth most common cancer among women. At the same time approximately 7400 of women were expected to die from this cancer in USA and 46 600 in Europe annually. These data enable to calculate that altogether 2.45% of women born today will be diagnosed with EC at some time during their lifetime [3, 4].
Since 1988, the Gynecologic Oncology Committee of the International Federation of Gynecology and Obstetrics (FIGO) has recommended surgical staging of EC based on exploratory laparotomy, total abdominal hysterectomy, bilateral salpingo-oophorectomy, peritoneal cytology, and pelvic and para-aortic lymphadenectomy with the pathologic stage adding extra information. Unfortunately staging analysis still does not provide medical doctors, nor the patient with a relevant prognosis . Multiple accessory prognostic factors have been defined for endometrial cancer, also including some molecular markers [2, 6]. Nevertheless still the diagnosis is frequently uncertain because of false-positive rates of up to 25% and false-negative rates of up to 10% in cervical invasion evaluation with MRI (Magnetic Resonance Imaging) and sensitivity of invasion detection ranging from 66% to 100% (mean, 86%), and specificity from 92% to 100% (mean, 97%) . Furthermore, the overall quality of surgical staging may be poor and very different from case to case as it is related to both the year that the surgeon passed the license examination and also to specialist status and experience . It the face of presented data the need for new medical formation programmes and also novel diagnostic and prognostic markers is evident. Although endometrial carcinoma is associated with a good prognosis because patients tend to present with early disease, high-risk populations may benefit from screening, but no prospective studies have demonstrated a benefit in any population untill now. Therefore it is interesting to develop new screening tools that may enable to select populations at high EC risk and support the process of prevention and early diagnosis . Until now PTEN (phosphatase and tensin homolog), K-ras (V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog), TP53 (tumor protein 53), β-catenin, MSH2 (MutS homolog 2, colon cancer, nonpolyposis type 2), MSH6 and Her2/neu (Human Epidermal growth factor Receptor 2); and mitochondrial gene mutations and protein signaling pathways have been implicated in the process of endometrial carcinogenesis [6, 10, 11]. On the basis of recent reports, it seems possible that a molecular mtDNA-analysis-based approach may be used in clinics in the future [6, 12–19].
The first interest in mitochondrial function in carcinogenesis was reported as early as in the 1920s, when Otto Warburg discovered that cancer cells have a high glycolytic rate and produce increased levels of lactate in the presence of oxygen. Since then for more than two-thirds of the last century we have known that a common biochemical signature of many tumours, particularly those that are poorly differentiated and proliferate rapidly, is their propensity to utilize glucose at high rates [20, 21]. This cancer characteristic has opened a new field of research today referred as to "mitochondrial medicine", since mitochondria are the metabolic organelles of the cells [22, 23]. At this point the mitochondrial genome (mtDNA) came into the focus of multiple projects. mtDNA somatic mutations were described to arise in the cells of various types of human cancers including bladder, brain, breast, colon, head and neck, lung, ovarian, prostate, or thyroid [6, 12, 15, 16, 18, 24–26]. At the same time inherited polymorphism have been pointed out as contributing factors in cancer development [24, 27, 28]. Nevertheless, the difficult task of correlating mtDNA polymorphisms and somatic mutations with neoplastic phenotype is not solved yet. The key role of the mitochondria in cell apoptotic pathways and the close link of tumour - suppressor proteins with mitochondria suggest some of the mechanisms of mitochondria dependent - tumourigenesis [6, 18, 19, 25, 26]. In particular, experiments on nude mice have shown that alterations in mtDNA (and a subsequent increase of ROS production) may contribute to cancer formation and development [16, 29].
Our current research has been inspired by a number of studies indicating that mtDNA analysis may be more powerful in detecting tumour cells in bodily fluids and cytological specimens than nuclear DNA analysis [17, 30]. Fliss and colleagues have reported facile analysis of mtDNA sequence in diagnostic samples . In the present study, we examined the distribution of mtDNA inherited polymorphisms in the D-loop region of mtDNA in ECs population. Moreover we have also analyzed the presence of somatic mutations in primary ECs samples. In addition we have also investigated the distribution of mitochondrial haplogroups (haplogroup specific polymorphisms) in the patient population. Our choice of complementation of polymorphism analysis by haplogroup investigation was based on the rationale that the same mechanisms which may operate to create variation in evolution can also operate in clonal evolution in tumours . Finally, in order to test for association with cancer susceptibility we compared mtDNA-data of EC patients with those for general Polish  and European populations [33, 34].
Altogether, our analysis on 26 patients has shown that seven (27%) belong to haplogroup U, five (19%) to haplogroup J, four (15%) to haplogroup K, three (12%) to haplogroup H, two (8%) to haplogroup T and two to haplogroup W (8%). No patient of haplogroup I, V or X was identified. In three (11%) cases, no haplogroup could be assessed due to unspecific polymorphisms found in haplogroup specific positions. Statistical analysis revealed that EC population haplogroup distribution is different from distribution found in the general Polish population. Most striking is the under-representation of haplogroup H individuals among EC patients (12 vs. 38%; p = 0.003). Other haplogroups seem to be represented in the cancer population at a frequency similar to the Polish population (Figure 1.).
As haplogroup H is characterized by T7028C polymorphism presence we have analyzed its statistical parameters. Patients with endometrial adenocarcinoma had the T7028C 12% of the time while in the general Polish population  carries the T7028C SNP 42% of the time. This difference is statistically significant p = 0.003 with Fisher's exact test, which was confirmed both by Yates corrected chi square χ2 = 7.382 and p = 0.007 and 'N-1' Chi square χ2 = 8.58 and p = 0.003, and these tests have been used as expected to give relatively low Type I error. The statistics was performed as previously suggested . Moreover in cancer patients T7028C is found at RR = 0.28 (95% CI: 0.093<R.R.<0.812), which is 72% fewer than in non-cancer cohort. This gives OR = 0.18 and an inverse OR = 5.53, which means that one has a 5.53 higher chance to be T7028C carrier and not develop cancer in comparison to people who developed cancer, at specificity 0.775 (0.755 - 0.815). 23% more of 7028C do not develop endometrial adenocarcinoma in comparison to 7028T carriers. Negative Predictive Value (NPV) of 7028C test result is 0.581 (0.567 - 0.611) with Relative Risk Reduction (RRR) = 0.725 and NND = 5.69. 7028T carriers have RR = 4.51 of cancer development if compared with 7028C, with OR = 5.5316. The 7028T inheritance correlates with cancer development with Yules-Q of 0.694, and therefore indicates a moderately strong positive relationship. C7028T test has sensitivity of 0.885, which guarantees the recognition of the majority of cancer predisposition as such, thus a negative result may be used to rule out the disease susceptibility. High negative predictive value of this test at 0.950 (0.881 - 0.982 at 95% CI) also seems to confirms the utility of this test. At the same time a relatively low positive predictive value of 0.223 seems typical for a multi-factorial disease such as cancer. At the same time this test is cost-effective, NND is as low as 3.29. Inheritance of 7028T results in RRI of EC development as high as 350% in comparison to 7025C. At the same time ARI is 17.5%.
In face of those calculations 7028C is to have protective value in EC development with RRR = 78% and ARR = 18%. This percentage of 7028C corresponds to negative likelihood ratio of 0.225 (0.185 - 0.245 at 95% CI) and is considered to show a small decrease in the likelihood of disease. Diagnostic Odds Ratio of C7028T is 5.53 and as shown before DOR of 3 - 6 is hardly a useful test for diagnostic purposes, but it may provide a strong clue in investigations of etiology . Nevertheless it has to be underlined that this research must be confirmed on large sample numbers to verify its relevance for clinical application. For evidence based medicine genotyping of a large number of patients is needed before any application can be suggested [17, 37].
Mutations in the D-loop mtDNA region
In the current study, the mtDNA D-loop sequence was analyzed by sequencing mtDNA from 26 cancerous and normal tissues. Both tissues were obtained from the patients with clinically diagnosed and pathologically confirmed EC during surgery. We were able to detect 2 somatic mutations (G16153A and A16188C). These genetic changes are point mutations and A16188C may be classified as mtMSI (mitochondrial satellite instability). As anticipated these mutations are located in mtDNA hypervariable region I (HV1 - 16024-16383). G16153A is a transition typically characteristic for mtDNA mutagenesis possibly induced by oxidative stress. In contrast, mtMSI might be generated during erroneous replication. These mutations seem to be relatively infrequent [12, 15] and therefore possibly are more specific for the EC patient population.
Germ-line polymorphisms in the D-loop mtDNA region
Our population of patients differs by 80 germ-line polymorphisms (Additional file 1; Table S1., Figure 2 and 3) from haplogroup H [38, 39], the most common haplogroup in Poland, and the haplogroup we have just shown to be at lower risk for EC development [32, 34]. Specifically only 22 out of 80 (28%) polymorphisms were generally common , and A16T polymorphism has not been reported previously . We believe that polymorphisms uncommon in the global population, but abundant in the endometrial carcinoma cohort are good candidates for cancer-susceptibility loci. As predicted, polymorphisms were predominantly located in mtDNA hypervariable regions HV1 (16024-16383, Figure 3) and HV2 (57-333, Figure 2) - 34 (43%) and 27 (34%) polymorphisms, respectively.
Our special attention was aimed at 73, 143, 189, 199, 203, 204, 207, 250, 295, 16051, 16069, 16126, 16129, 16172, 16192, 16223, 16224, 16256, 16270, 16278, 16292, 16311, 16356, 16391 mtDNA position polymorphisms, as these are D-loop located haplogroup specific loci [40–42]. This analysis revealed an overabundance of specific polymorphisms in our cancer study population, in particular 16126C - specific for J or T haplogroups (p = 0.025); and 207A - typical for haplogroup W (p = 0.027). At the same time 16223T, specific for haplogroups W, X and I (p = 0.005) was underrepresented. In summary 16223C was present in 80% cases (21/26), 207A was found in 19% (5/26) and 16126C in 23% (6/26), so there might be a correlation and/or association between 16223C/207A genotype and EC development.
A number of patients harboured other unspecific polymorphisms, including 73A in patient 402 belonging to haplogroup T (G specific for haplogroup T), 16356C in patient 423 belonging to haplogroup U (T specific for haplogroup U), 16311C in patient 425 belonging to haplogroup W (C specific for haplogroups I, K, U, J and H) and 207A in patient 427 belonging to haplogroup K (G specific for haplogroup H) (data not shown). These polymorphisms would introduced a bias in haplogroup analysis if no data from RFLP and multiplex-PCR methods were available, which proves the need for whole-genome based analysis of haplogroup assignment.
Polymorphisms in the coding region
In addition to D-loop polymorphism, our study reveals we were also able to detect 11 additional germ-line polymorphisms (Additional file 2; Table S2): four in 12S rRNA, one in ND1, ND2, ND3, CO3, tRNA-Arg, tRNA-Pro, which differ in our population from rCRS (haplogroup H). In this analysis polymorphism 4216C affirms the assignment of patients numbered 401, 403, 409, 410, 418, 428 to haplogroups J and T. Remarkably, position 10463, reported here as polymorphic, has been previously found subject to somatic mutation in EC by our group and A15960T has not been reported previously .
The use of mitochondrial DNA mutation and/or polymorphism patterns as a biomarker is rapidly expanding in disciplines ranging from rare metabolic diseases, aging, to cancer and the tracing of human migration patterns, population characterization, and human identification in forensic sciences [12, 14, 15]. In the case of endometrial adenocarcinoma the presence of mtDNA mutations has been analyzed by three independent research groups from China [44–46], and one from Italy . All those studies have been rather limited. They included the analysis of: position 310 (6 cases) , nucleotide positions 514, 16189, and 16519 (51 cases) , D-loop region, the 12S and 16S rRNA genes (50 cases) , 12 microsatellite markers starting at nucleotide position 303, 514, 956, and 16184 (17 cases) , hypervariable segments of the control region, parts of 16S rRNA, tRNA(Leu) and ND1 gene (49 cases) . In summary, all those papers reported one or more point mutations in 48.4%  to 56% of cases  and indicated mononucleotide repeat (D310) as a mutation 'hot-spot' in primary tumours . In those reports no correlation between the pattern of mitochondrial abnormalities and clinical or histo-pathological features was found . The mtDNA analysis of polymorphism distribution in endometrial adenocarcinoma was not as thorough and only the 16189T>C polymorphism was associated with EC susceptibility .
EC research [43, 51] performed previously in our laboratory on a different patient cohort employed SSCP and eventually subsequent sequencing and covered ~10% of mtDNA, including nucleotides 135-433, 2986-3301, 4981-5500, 10390-10700 and 12005-12386. In this previously tested patient cohort, we found homoplasmic somatic mtDNA mutations in four out of 48 (8%) EC analyzed. These included alterations in ND2 (G5231A), ND4L (G10550A, T10640C, T10551C) and tRNAs (A12308G tRNA-Leu, C12258G tRNA-Ser, T10463C tRNA-Arg) [43, 51]. Our current research on the D-loop reveals two mtDNA mutations (again only 8% patients) and shows multiple polymorphisms, including microsatellite instability in the C-tract (between nucleotides 303 and 315) in 22 (85%) of EC patients. The frequency of mtDNA mutations found in our previous analysis and in our new patient populations are similar, which strengthens the probability of low-mutation load in EC patients. Also in the report by Schwartz et al. mutations in the mitochondrial noncoding polycytidine (C)(n) repeat (polyC) analyzed in endometrium tumours (n = 53) were present in 11% of patients . The frequency of D-loop mtDNA mutations reported by us and Schwartz et al. is much lower than that reported in other studies where the percentage was as high as 32.2% to 56% [49, 50]. The frequency of mtDNA mutations in our study also differs from that reported in other human adenocarcinoma cases, including mutations found in 37% cases of gastric tumours , 40% - in Barrett's carcinomas , and 52% - in prostate cancers . We believe that low number of mutations we detected is due to the stringent quality control we used [17, 56] and reflects real number of mtDNA mutations in EC patients. In the light of that we believe that a great role in cell transformation must be attributed to inherited mtDNA polymorphisms.
A similar result of high polymorphism rate and low mutation load within cancer cells was reported in prostate  and thyroid cancer . Also the risk of oral cancer was increased in patients bearing some of the common mitochondrial polymorphisms, not the mtDNA mutations . As many germline mtDNA polymorphisms are found in the human population  it is therefore worthwhile to consider their differential occurrence in cancer-susceptible population with emphasis on their utility as potential cancer-predisposition markers [58, 59]. In support of that hypothesis recent large review by the Wallace group revealed that a major part (72%) of previously reported tumour-specific somatic mtDNA mutations are actually mtDNA sequence variants found in the general population. The authors claim that 52% of the tumour somatic mRNA missense mutations, 83% of the tRNA, 38% of the rRNA, and 85% of the D-loop mutations are actually common sequence variants . The report by Wallace and co-workers seems to favor our results of polymorphisms analysis.
Our sequencing analysis has been focused on the D-loop region of mtDNA. The mtDNA D-loop region is highly polymorphic, contains two hypervariable regions: HV1 (16024-16383) and HV2 (57-333) and was reported as somatic mutation "hot spot" in many types of cancer . This directed our attention on this mtDNA region. We need to stress that D-loop contains crucial elements for mtDNA replication and transcription, and therefore sequence differences in this region might alter the rate of DNA replication by modifying the binding affinity of important trans-activating factors . In fact, both down- and up-regulation of mtDNA replication was found in several tumour types [62, 63]. Moreover it was shown that in some cases reduced copy number of mtDNA correlates with clinical outcome  and might be potentially used to predict prognosis . In the above-mentioned papers the D-loop sequence variability was one of key contributing factors leading to decreased mtDNA level in breast tumours. In particular, the 303-315 stretch highly polymorphic in our patient cohort localizes in the conserved sequence block II (CBSII) and is the site of replication primer binding. Although it is not clear what the long-term impact of mtMSI in this locus is, it is known that replication efficacy is altered depending on the sequence length (poly-C polymorphisms), as a result of CSB II involvement in replication priming. CBSII heavy-strand contributes to formation of a persistent RNA-DNA hybrid that serves to prime mtDNA replication. The formation of RNA-DNA hybrids is dependent on the GC-rich element. Interestingly, efficient hybrid formation is also influenced by sequences 5' to the hybrid, including the CSBIII element . Furthermore mtDNA transcription required to generate RNA primers used in the initiation of heavy strand DNA synthesis is critically dependent on the exact CSB II sequence . Premature transcription termination occurs if particular mtMSI are found in positions 300 to 282 of the mtDNA sequence and may be completely abolished in the 319-289 mutants. In contrast, 304-300 mutants show a drastic decrease in transcription termination . All these data support our hypothesis that mtDNA polymorphisms, in particular D-loop polymorphisms influence cell physiology and may result in a pro-carcinogenic phenotype of the carrier. In addition we believe that the physiological influence of D-loop (CA)(n) polymorphism is significant in a clinical perspective, not only in cell biology. This phenomenon was already identified in breast cancer population, where patients with multiple alleles of the mtDNA D-loop (CA)(n) polymorphism (state called heteroplasmy) had significantly poorer disease-free survival than those with one allele of the mtDNA D-loop (CA)(n) polymorphism. This result suggests that the mtDNA D-loop polymorphisms may be associated with cancer survival and we believe this might also be true for EC patients and should be investigated in nearest future .
The other position that was highly polymorphic in our research - 16189 has previously been in the focus of oncological research. The carriers of germ-line T to C polymorphism at 16189 are to be more susceptible to breast cancer and ganglioma development in the light of the high frequency of 16189C detected in cancer patients and low number of 16189C in healthy individuals . Interestingly, in the report by Liu et al. the T16189C polymorphism was found in 14% of ECs. In addition this polymorphism was also linked with an increased risk of type II diabetes mellitus, which is well known to be a risk factor in the aetiology of endometrial and breast cancers [48, 70]. In the screen of female cancers, deletions or insertions were detected in the poly C tract of tumour mtDNA and were linked in all cases to a germline T to C transition at 16189. The elimination of the T generates a long poly-C tract that is generally anticipated to exhibit higher instability. It was subsequently empirically indicated that this position represents a strong hotspot of mtMSI with a (C)7-14 pattern of variation . Furthermore, nucleotides 16184 - 16193 are located on the 3'-end of a termination-associated sequence (TAS) and at the 7S DNA binding site which are thought to be involved in the regulation of mtDNA synthesis . Again if polymorphisms at 16189 are in the focus, our hypothesis of mitochondrial polymorphism involvement in cancer development is again supported in many respects. The second part of our research focused on mtDNA mutation screen and we have detected A16188C which is also localized within 16184-16193 tract involved in the regulation of mtDNA synthesis. This locus was previously reported as mutated in breast cancer (16188 C→CC)  and listed as must-be-excluded in haplogroup analysis and usage of D-loop mutations as molecular clock .
The CCCCCTCCCC sequences located in the Hypervariable Regions I and II of the D-loop, but also the 12S rRNA gene (Additional file 2; Table S2.) emerge as instability 'hot spot' regions in endometrial carcinomas. Considering other mtDNA polymorphisms reported by us it is worth noting that some polymorphic sequences within mitochondrial tRNA genes (Additional file 2; Table S2.) may serve as unusual replication origins with pathogenic implications. Today we know that particular variants of tRNA genes may acquire secondary structures resembling mitochondrial origins of light strand replication, particularly structures that might invoke bi-directional replication. In consequence, this excessive replication may cause abundant mutations in genome regions not adapted to tolerate them . At the same time context analysis for 303 polymorphisms revealed a complex influence of neighbouring bases on mutagenesis in the HVS I region. Analysis suggested that a transient misalignment dislocation mutagenesis operates in monotonous runs of nucleotides and plays an important role in generating base substitutions in mitochondrial DNA . In our opinion all discussed mutation/polymorphism generation mechanisms seem to be in important in EC and have potential impact on cell physiology.
Finally we need to emphesize that as previously pointed by Salas and Bandelt  mtDNA sequencing and analysis techniques contain inherent problems, particularly with regard to the generation of authentic and useful data. In fact we understand that for the medical field it is important to employ standardized procedures based on scientific grounds, in order to have mtDNA-based evidence that may be accepted in clinics. Therefore we have put much effort in refinement of our amplification and sequencing strategies, as well as a posteriori quality control of mtDNA sequencing [17, 77]. We believe that such sequencing efforts resulted in low number of mutations, as we have eliminated phantom mutations . To further refine the authenticity of our sequencing results a posteriori, we utilized the strategy of focused database comparisons, as this method has been proven to be effective and successful in the case of modern mtDNA data [56, 76, 78]. Considering the critical analysis of mtDNA sequencing data [76, 78] we would like to stress that our summary (additional file 2; Table S2 and additional file 3; Table S3.) should be treated and interpreted with caution, as probaly some previously reported cancer mutations we cite are actually misreported nucleotide variants [78, 79].
The second part of our project covers haplogroup related polymorphism analysis (Additional file 4; Table S4.). The frequencies of mtDNA haplogroups vary between ethnic groups worldwide. Our strictly European population is exclusively distributed among the nine haplogroups designated: H, I, J, K, T, U, V, W, and X, whereas haplogroups A, B, C, D, and E that are characteristic for Asian population or haplogroups L1, L2, and L3 specific for the African population were not analyzed . The research was inspired by the fact that some haplogroup specific polymorphisms are significantly correlated with increased or decreased risk of specific human disorders, for example - LHON and Parkinson's disease with haplogroup J [80, 81], ALS with haplogroup I , AD in males with haplogroup U , and a significant increase in breast cancer development with haplogroup K . Concurrently correlation between longevity and certain haplogroups (J and U) was reported  (Additional file 4; Table S4.). Those facts prompt us to think that haplogroup mtDNA variants could be specifically connected with cancer, in particular EC. We have chosen to thoroughly analyze haplogroups in our EC cohort, as polymorphism A5178C - characteristic for haplogroup D (typical for South-East Asia) was correlated with the development of EC. The woman carrying 5178C seems to be at a significantly higher risk of EC development . Nevertheless the Chinese population has a different mitochondrial genetic background than the European or white American population, and as a result A5178C polymorphism - EC correlation is not exploitable or applicable as a clinical marker outside the region where it was observed. Until now there was no reported research on the haplogroup distribution in the Caucasian EC population (Medline© search database). As we expected, the EC cohort haplogroup distribution was not similar to the healthy Polish population with underrepresented haplogroup H individuals among EC patients. We believe that this might suggest a protective role of polymorphisms typical for haplogroup H in cancer development. Our hypothesis seems plausible as mtDNA haplogroup H is very common in Caucasoids, reaching frequencies of ~50%, which certainly suggests that it may confer some advantage. Recent studies have suggested haplogroup H3 is highly protective against AIDS progression . Moreover, mtDNA haplogroup H is a strong independent predictor of outcome during severe sepsis, conferring a over 2-fold increased chance of survival at 180 days compared with individuals with different haplogroups .
We believe that inherited mtDNA polymorphisms, both in the D-loop, and in the coding region (including haplo-group specific polymorphisms) may cause subtle differences in the encoded protein structure and function; and thus subtle changes in OXPHOS activity and free-radical production. It is therefore likely that mtDNA polymorphisms in mitochondrial genes involved in electron transport chain and oxidative phosphorylation result in increased oxidative stress and hypermutagenesis of mitochondrial as well as nuclear DNA. This predisposes an individual or population sharing the same mtDNA genotype to an earlier onset of degenerative cellular processes, such as the accumulation of somatic mtDNA variation, decline in OXPHOS capacity or faster cancer progression, as shown in a cybrid model .
In conclusion we suggest that mitochondrial research will enable to establish bio-markers helping to identify individuals at high risk for developing specific cancer types and to develop screening approaches for early diagnosis of cancer [14, 24, 84, 88]. Molecular assessment of mitochondrial abnormalities of cancer cells could represent a promising tool not only for prognosis and early diagnosis of neoplasia, but possibly also during the diagnosis or follow-up of cancer patients [15, 16, 19, 89, 90]. A number of studies have proven that the mitochondrial genome is more useful in detecting tumour cells in bodily fluids and cytological specimens than mutations in nuclear DNA  in bodily fluids  and we believe this can be applied to EC screening. Thus, we suggest that application of mtDNA polymorphism pattern analysis may be useful to select populations at increased risk of cancer development and in cancer patients render final diagnosis in those cases where the morphological abnormalities are unspecific or neoplastic cells are at difficult to detect. In those cases, the application of PCR-coupled with gel electrophoresis or DNA sequencing may lend itself to rapid analysis of multiple samples. In the clinical context, the high frequency of mitochondrial genome instability, in combination with PCR-based assays of high sensitivity, may be of potential clinical usefulness [22, 56].
Altogether we investigated twenty-seven EC patients who had undergone surgery (total abdominal hysterectomy and bilateral oophorectomy) at the IInd Department of Gynaecology, Medical University at Lublin, Lublin, Poland, between 2003 and 2006. Pelvic and para-aortic lymph node dissections were performed if material obtained at dilatation and curettage was diagnosed as non-endometrioid or poorly (G3) differentiated cancer, or when the neoplasm invaded over one-half of the myometrial thickness of the uterus at surgery. In selected cases, the patients underwent additional surgical procedures, including omentectomy, appendectomy, tumour cytoreduction or dissection of distant metastases. Non-neoplastic reference material was available in all cases. The patients received no chemotherapy, radiotherapy, or hormonal therapy before surgery. The clinical stage of the disease was classified according to the staging system of the FIGO . The material was assessed histologically at the Department of Pathology, Medical University in Lublin, Lublin, Poland, based on the WHO classification . The pathological findings were assessed with regard to histological type and grade, depth of myometrial infiltration, pattern of ovarian involvement, vascular space invasion, presence of the neoplastic material in the fallopian tube, and the presence of metastases in the ovaries and/or lymph nodes (Additional file 5; Table S5). summarizes the clinico-pathological features of EC patients enrolled in this research. The clinical records were reviewed in order to obtain information regarding patient age, menopausal status, presenting symptoms and surgical management.
Prior to the surgery, all patients had provided informed consent to the use of their post-surgical material. The project was approved by the local Ethics Committee at the Medical University of Warsaw, Warsaw, Poland (KB-0/6/2007 to AMC).
Tissue collection and DNA isolation
At surgery, after the uterus was removed, the uterine corpus was gently cut and the neoplastic material was scraped to a sterile eppendorf tube. During this procedure, material collected from the uterine cavity was not cross-contaminated by cervical cells. Tissue obtained at surgery was sub-divided into two parts. One portion was fixed in buffered formalin (pH 7.4) for routine histopathological assessment, while the rest was immediately frozen in liquid nitrogen and stored at -80°C. Non-neoplastic material was also collected at surgery and stored at -80°C until assayed.
DNA was isolated by standard proteinase K treatment followed by phenol/chloroform/isoamyl alcohol extraction. DNA was precipitated with 0.3 M sodium acetate in 70% ethanol at -20°C overnight and resuspended in Tris-EDTA (TE) buffer (pH 8.0). DNA quantification was performed with measurements at an absorbance of 260 nm.
Mutation and polymorphism analysis
A germline (inherited) polymorphism is described as a difference between normal tissue variant and the revised Cambridge Reference Sequence (rCRS) - present both in normal and tumor tissues. Whenever the change between sequences from tumor (investigated) sample and normal tissue occurs, it is defined as a somatic mtDNA mutation. Obtained alteration frequencies were compared with mtDB database data . In all positive samples (mutation), sequencing reactions were validated by a new independent amplification and sequencing.
Haplogroup analysis by RLFP
Common polymorphisms in mtDNA determining classes of related genotypes, referred to as haplogroups, were detected by restriction fragment length polymorphism (RFLP) analysis. Haplogroup RLFP analysis, restriction enzymes and primers used are summarized in Additional file 6; Table S6.
Haplogroup analysis by multiplex-PCR/sequencing
To verify haplogroups established by RFLP (Additional file 7; Table S7 and Additional file 4; Table S4.), multiplex-PCR/sequencing was performed as described previously . Moreover, haplogroups were also assigned based on specific D-loop polymorphisms according to published data [40–42]. If any unspecific RFLP or multiplex-PCR/sequencing variants were found, sequencing with appropriate primers listed in Additional file 8; Table S8. was performed in order to verify the mutation/polymorphism.
PCR amplification of D-loop segment of mtDNA
mtDNA fragment (15587 - 964) containing the D-loop region (spanning nucleotide 16024 to 576) was amplified using eight pairs of primers. The primer pairs used and the sizes of the amplified products are shown in Additional file 7; Table S7.
Fifty-microlitre reactions contained 10 ng DNA and 0.5 μM primers, 0.2 mM each of deoxynucleotide triphosphate (dNTP), 1U of FIREPol® DNA Polymerase (Solis BioDyne, Estonia) or Pfu DNA Polymerase (Fermentas AB, Lithuania) and 2.5 mM MgCl2. DNA was subjected to the following cycling conditions: initial denaturing at 95°C for 3 min followed by 94°C for 1 min, 55°C for 30 sec, and 72°C for 1 min for 40 cycles and final extension step at 72°C for 7 min. Two microlitres of PCR products was analysed on an ethidium bromide-stained, 3% agarose gel (40 min at 70 V) to demonstrate the presence of the amplification product and for its quantification.
mtDNA sequence analysis
Sequence analysis was performed by: FinchTV Version 1.4.0 (Geospiza Inc., USA) and BioEdit version 184.108.40.206 (Copyright Tom Hall 1999-2007), contig assembly was performed with Sequencher 4.1.4 (Gene Codes Corporation, Ann Arbor, MI USA) and multiple sequence alignment was performed with Clustal W . Normal and cancer tissue mtDNA sequences were compared with the revised Cambridge Reference Sequence (CRS) and sequence variants were recorded [38, 39].
Two tailed non-directional Fisher-Irwin (Fisher's exact test) was used for statistical analysis . Statistical analysis was performed with PAST - PAlaeontological STatistics, ver. 1.34 (Øyvind Hammer, D.A.T. Harper and P.D. Ryan, 2005) and Analyse-it for Microsoft Excel General & Clinical Laboratory modules Version 1.73 (Analyse-it Software, Ltd. Copyright © 1997-2005). The difference was considered statistically significant if p < 0.05. In selected cases, to confirm the result Fisher's test, Yates's chi and un-corrected chi squared test ('N - 1' chi squared test) have been used as expected to give relatively low Type I error in case of a small research cohort. The statistics was performed as previously suggested . To further understand the significance of T7028C as a factor for favorable outcomes (odds ratio, relative risk, difference in proportions, absolute and relative reduction in risk) and of the effectiveness of a diagnostic criteria (number needed to diagnose, specificity, positive & negative predictive values, positive & negative likelihood ratios, diagnostic and error odds ratios) was performed. The parameters, as well as the confidence intervals for the estimated parameters are computed by a general method [97, 98].
conserved sequence block 2
Hypervariable region I
Reactive Oxygen Species
Restriction Fragment Length Polymorphism
mitochondrial satellite instability
Relative Risk, Risk Ratio
Diagnostic Odds Ratio
Positive Predictive Value
Negative Predictive Value
Relative Risk Reduction
Absolute Risk Reduction
Relative Risk Increase
Absolute Risk Increase
Number Needed to Diagnose.
Canavan TP, Doshi NR: Endometrial cancer. Am Fam Physician 1999,59(11):3069–3077.
Ronnett BM, Zaino RJ, Ellenson LH, Kurman RJ: Endometrial carcinoma. Blaustein's Pathology of the Female Genital Tract (Edited by: Kurman RJ). New York, Berlin, Heidelberg: Springer 2002, 501–559.
Ferlay J, Autier P, Boniol M, Heanue M, Colombet M, Boyle P: Estimates of the cancer incidence and mortality in Europe in 2006. Ann Oncol 2007,18(3):581–592.
Ries L, Melbert D, Krapcho M, Mariotto A, Miller B, Feuer E, Clegg L, Horner M, Howlader N, Eisner M, et al.: SEER Cancer Statistics Review. Bethesda: National Cancer Institute 2007.
Hacker NF: Uterine cancer. Practical Gynecologic Oncology (Edited by: Berek JS, Hacker NF). Philadelphia: Lippincott Williams & Wilkins 2004, 397–442.
Klemba A, Kukwa W, Bartnik E, Krawczyk T, Scinska A, Golik P, Czarnecka AM: [Molecular biology of endometrial carcinoma]. Postepy Hig Med Dosw (Online) 2008, 62:420–432.
Kinkel K, Kaji Y, Yu KK, Segal MR, Lu Y, Powell CB, Hricak H: Radiologic staging in patients with endometrial cancer: a meta-analysis. Radiology 1999,212(3):711–718.
Crawford SC, De Caestecker L, Gillis CR, Hole D, Davis JA, Penney G, Siddiqui NA: Staging quality is related to the survival of women with endometrial cancer: a Scottish population based study. Deficient surgical staging and omission of adjuvant radiotherapy is associated with poorer survival of women diagnosed with endometrial cancer in Scotland during 1996 and 1997. Br J Cancer 2002,86(12):1837–1842.
Whitcomb BP: Gynecologic malignancies. Surg Clin North Am 2008,88(2):301–317.
Doll A, Abal M, Rigau M, Monge M, Gonzalez M, Demajo S, Colas E, Llaurado M, Alazzouzi H, Planaguma J, et al.: Novel molecular profiles of endometrial cancer-new light through old windows. J Steroid Biochem Mol Biol 2008,108(3–5):221–229.
Hecht JL, Mutter GL: Molecular and pathologic aspects of endometrial carcinogenesis. J Clin Oncol 2006,24(29):4783–4791.
Brandon M, Baldi P, Wallace DC: Mitochondrial mutations in cancer. Oncogene 2006,25(34):4647–4662.
Hendrickson SL, Hutcheson HB, Ruiz-Pesini E, Poole JC, Lautenberger J, Sezgin E, Kingsley L, Goedert JJ, Vlahov D, Donfield S, et al.: Mitochondrial DNA haplogroups influence AIDS progression. AIDS 2008,22(18):2429–2439.
Wallace DC: A mitochondrial paradigm of metabolic and degenerative diseases, aging, and cancer: a dawn for evolutionary medicine. Annu Rev Genet 2005, 39:359–407.
Czarnecka A, Golik P, Bartnik E: Mitochondrial DNA mutations in human neoplasia. J Appl Genet 2006,47(1):67–78.
Czarnecka AM, Bartnik E: Mitochondrial DNA Mutations in Tumors. Cellular Respiration and Carcinogenesis (Edited by: Apte SP, Sarangarajan R). New York City: Humana Press 2009, 1–12.
Czarnecka AM, Krawczyk T, Czarnecki JS, Kukwa W, Ścińska A, Ribbene A, Lo Verde R, Sunseri A, Peri G: Methodology For Mitochondrial DNA Research In Oncology: Goals And Pitfalls. ARS Medica Tomitana 2008,XIV(1):48–64.
Czarnecka AM, Marino Gammazza A, Di Felice V, Zummo G, Cappello F: Cancer as a "Mitochondriopathy". J Cancer Mol 2007, 3:71–79.
Plak K, Czarnecka AM, Krawczyk T, Golik P, Bartnik E: Breast cancer as a mitochondrial disorder (Review). Oncol Rep 2009,21(4):845–851.
Shaw RJ: Glucose metabolism and cancer. Curr Opin Cell Biol 2006,18(6):598–608.
Warburg O: On respiratory impairment in cancer cells. Science 1956,124(3215):269–270.
DiMauro S: Mitochondrial DNA medicine. Biosci Rep 2007,27(1–3):5–9.
Smeitink JA, Zeviani M, Turnbull DM, Jacobs HT: Mitochondrial medicine: a metabolic perspective on the pathology of oxidative phosphorylation disorders. Cell Metab 2006,3(1):9–13.
Czarnecka AM, Krawczyk T, Zdrozny M, Lubinski J, Arnold RS, Kukwa W, Scinska A, Golik P, Bartnik E, Petros JA: Mitochondrial NADH-dehydrogenase subunit 3 (ND3) polymorphism (A10398G) and sporadic breast cancer in Poland. Breast Cancer Res Treat 2009.
Pietka G, Kukwa W, Bartnik E, Scinska A, Czarnecka AM: [Mitochondrial DNA mutations in the pathogenesis in the head and neck squamous cell carcinoma]. Otolaryngol Pol 2008,62(2):158–164.
Plak K, Kukwa W, Bartnik E, Golik P, Scinska A, Krawczyk T, Czarnecka AM: [The impact of mtDNA mutations on proteins structure in selected types of cancer]. Postepy Biochem 2008,54(2):151–160.
Datta S, Majumder M, Biswas NK, Sikdar N, Roy B: Increased risk of oral cancer in relation to common Indian mitochondrial polymorphisms and Autosomal GSTP1 locus. Cancer 2007,110(9):1991–1999.
Singh KK, Kulawiec M: Mitochondrial DNA polymorphism and risk of cancer. Methods Mol Biol 2009, 471:291–303.
Petros JA, Baumann AK, Ruiz-Pesini E, Amin MB, Sun CQ, Hall J, Lim S, Issa MM, Flanders WD, Hosseini SH, et al.: mtDNA mutations increase tumorigenicity in prostate cancer. Proc Natl Acad Sci USA 2005,102(3):719–724.
Ha PK, Tong BC, Westra WH, Sanchez-Cespedes M, Parrella P, Zahurak M, Sidransky D, Califano JA: Mitochondrial C-tract alteration in premalignant lesions of the head and neck: a marker for progression and clonal proliferation. Clin Cancer Res 2002,8(7):2260–2265.
Fliss MS, Usadel H, Caballero OL, Wu L, Buta MR, Eleff SM, Jen J, Sidransky D: Facile detection of mitochondrial DNA mutations in tumors and bodily fluids. Science 2000,287(5460):2017–2019.
Piechota J, Tonska K, Nowak M, Kabzinska D, Lorenc A, Bartnik E: Comparison between the Polish population and European populations on the basis of mitochondrial morphs and haplogroups. Acta Biochim Pol 2004,51(4):883–895.
Ingman M, Gyllensten U: mtDB: Human Mitochondrial Genome Database, a resource for population genetics and medical sciences. Nucleic Acids Res 2006, (34 Database):D749–751.
Ruiz-Pesini E, Lott MT, Procaccio V, Poole JC, Brandon MC, Mishmar D, Yi C, Kreuziger J, Baldi P, Wallace DC: An enhanced MITOMAP with a global mtDNA mutational phylogeny. Nucleic Acids Res 2007, (35 Database):D823–828.
Campbell I: Chi-squared and Fisher-Irwin tests of two-by-two tables with small sample recommendations. Stat Med 2007,26(19):3661–3675.
Pepe MS, Janes H, Longton G, Leisenring W, Newcomb P: Limitations of the odds ratio in gauging the performance of a diagnostic, prognostic, or screening marker. Am J Epidemiol 2004,159(9):882–890.
Czarnecka AM, Krawczyk T, Czarnecki JS, Kukwa W, Ścińska A, Ribbene A, Lo Verde R, Sunseri A, Peri G: Methodology For Mitochondrial DNA Research In Oncology: Goals And Pitfalls. ARS Medica Tomitana 2008,XIV(1 (52)):48–64.
Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, et al.: Sequence and organization of the human mitochondrial genome. Nature 1981,290(5806):457–465.
Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N: Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet 1999,23(2):147.
Finnila S, Lehtonen MS, Majamaa K: Phylogenetic network for European mtDNA. Am J Hum Genet 2001,68(6):1475–1484.
Macaulay V, Richards M, Hickey E, Vega E, Cruciani F, Guida V, Scozzari R, Bonne-Tamir B, Sykes B, Torroni A: The emerging tree of West Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs. Am J Hum Genet 1999,64(1):232–249.
Torroni A, Huoponen K, Francalacci P, Petrozzi M, Morelli L, Scozzari R, Obinu D, Savontaus ML, Wallace DC: Classification of European mtDNAs from an analysis of three European populations. Genetics 1996,144(4):1835–1850.
Lorenc A, Bryk J, Golik P, Kupryjanczyk J, Ostrowski J, Pronicki M, Semczuk A, Szolkowska M, Bartnik E: Homoplasmic MELAS A3243G mtDNA mutation in a colon cancer sample. Mitochondrion 2003,3(2):119–124.
Li HX, Zhong S, Li CH: [Study on the mitochondrion DNA mutation in tumor tissues of gynecologic oncology patients]. Zhonghua Fu Chan Ke Za Zhi 2003,38(5):290–293.
Liu VW, Shi HH, Cheung AN, Chiu PM, Leung TW, Nagley P, Wong LC, Ngan HY: High incidence of somatic mitochondrial DNA mutations in human ovarian carcinomas. Cancer Res 2001,61(16):5998–6001.
Xu L, Hu Y, Chen B, Tang W, Han X, Yu H, Xiao C: Mitochondrial polymorphisms as risk factors for endometrial cancer in southwest China. Int J Gynecol Cancer 2006,16(4):1661–1667.
Parrella P, Seripa D, Matera MG, Rabitti C, Rinaldi M, Mazzarelli P, Gravina C, Gallucci M, Altomare V, Flammia G, et al.: Mutations of the D310 mitochondrial mononucleotide repeat in primary tumors and cytological specimens. Cancer Lett 2003,190(1):73–77.
Liu VW, Wang Y, Yang HJ, Tsang PC, Ng TY, Wong LC, Nagley P, Ngan HY: Mitochondrial DNA variant 16189T>C is associated with susceptibility to endometrial cancer. Hum Mutat 2003,22(2):173–174.
Liu VW, Yang HJ, Wang Y, Tsang PC, Cheung AN, Chiu PM, Ng TY, Wong LC, Nagley P, Ngan HY: High frequency of mitochondrial genome instability in human endometrial carcinomas. Br J Cancer 2003,89(4):697–701.
Wang Y, Liu VW, Tsang PC, Chiu PM, Cheung AN, Khoo US, Nagley P, Ngan HY: Microsatellite instability in mitochondrial genome of common female cancers. Int J Gynecol Cancer 2006,16(Suppl 1):259–266.
Semczuk A, Lorenc A, Putowski L, Futyma K, Bryk J, Miotla P, Bartnik E: Clinicoprognostical features of endometrial cancer patients with somatic mtDNA mutations. Oncol Rep 2006,16(5):1041–1045.
Schwartz S Jr, Alazzouzi H, Perucho M: Mutational dynamics in human tumors confirm the neutral intrinsic instability of the mitochondrial D-loop poly-cytidine repeat. Genes Chromosomes Cancer 2006,45(8):770–780.
Alonso A, Martin P, Albarran C, Aquilera B, Garcia O, Guzman A, Oliva H, Sancho M: Detection of somatic mutations in the mitochondrial DNA control region of colorectal and gastric tumors by heteroduplex and single-strand conformation analysis. Electrophoresis 1997,18(5):682–685.
Miyazono F, Schneider PM, Metzger R, Warnecke-Eberz U, Baldus SE, Dienes HP, Aikou T, Hoelscher AH: Mutations in the mitochondrial DNA D-Loop region occur frequently in adenocarcinoma in Barrett's esophagus. Oncogene 2002,21(23):3780–3783.
Gomez-Zaera M, Abril J, Gonzalez L, Aguilo F, Condom E, Nadal M, Nunes V: Identification of somatic and germline mitochondrial DNA sequence variants in prostate cancer patients. Mutat Res 2006,595(1–2):42–51.
Salas A, Carracedo A, Macaulay V, Richards M, Bandelt HJ: A practical guide to mitochondrial DNA error prevention in clinical, forensic, and population genetics. Biochem Biophys Res Commun 2005,335(3):891–899.
Maximo V, Soares P, Lima J, Cameselle-Teijeiro J, Sobrinho-Simoes M: Mitochondrial DNA somatic mutations (point mutations and large deletions) and mitochondrial DNA variants in human thyroid pathology: a study with emphasis on Hurthle cell tumors. Am J Pathol 2002,160(5):1857–1865.
Bandelt HJ, Salas A, Bravi CM: What is a 'novel' mtDNA mutation--and does 'novelty' really matter? J Hum Genet 2006,51(12):1073–1082.
Raule N, Sevini F, Santoro A, Altilia S, Franceschi C: Association studies on human mitochondrial DNA: methodological aspects and results in the most common age-related diseases. Mitochondrion 2007,7(1–2):29–38.
Hibi K, Nakayama H, Yamazaki T, Takase T, Taguchi M, Kasai Y, Ito K, Akiyama S, Nakao A: Detection of mitochondrial DNA alterations in primary tumors and corresponding serum of colorectal cancer patients. Int J Cancer 2001,94(3):429–431.
Pejovic T, Ladner D, Intengan M, Zheng K, Fairchild T, Dillon D, Easley S, Marchetti D, Schwartz P, Lele S, et al.: Somatic D-loop mitochondrial DNA mutations are frequent in uterine serous carcinoma. Eur J Cancer 2004,40(16):2519–2524.
Lee HC, Li SH, Lin JC, Wu CC, Yeh DC, Wei YH: Somatic mutations in the D-loop and decrease in the copy number of mitochondrial DNA in human hepatocellular carcinoma. Mutat Res 2004,547(1–2):71–78.
Liang BC, Hays L: Mitochondrial DNA copy number changes in human gliomas. Cancer Lett 1996,105(2):167–173.
Simonnet H, Alazard N, Pfeiffer K, Gallou C, Beroud C, Demont J, Bouvier R, Schagger H, Godinot C: Low mitochondrial respiratory chain content correlates with tumor aggressiveness in renal cell carcinoma. Carcinogenesis 2002,23(5):759–768.
Yu M, Zhou Y, Shi Y, Ning L, Yang Y, Wei X, Zhang N, Hao X, Niu R: Reduced mitochondrial DNA copy number is correlated with tumor progression and prognosis in Chinese breast cancer patients. IUBMB Life 2007,59(7):450–457.
Xu B, Clayton DA: RNA-DNA hybrid formation at the human mitochondrial heavy-strand origin ceases at replication start sites: an implication for RNA-DNA hybrids serving as primers. EMBO J 1996,15(12):3135–3143.
Wang Y, Liu VW, Ngan HY, Nagley P: Frequent occurrence of mitochondrial microsatellite instability in the D-loop region of human cancers. Ann N Y Acad Sci 2005, 1042:123–129.
Pham XH, Farge G, Shi Y, Gaspari M, Gustafsson CM, Falkenberg M: Conserved sequence box II directs transcription termination and primer formation in mitochondria. J Biol Chem 2006,281(34):24647–24652.
Ye C, Gao YT, Wen W, Breyer JP, Shu XO, Smith JR, Zheng W, Cai Q: Association of mitochondrial DNA displacement loop (CA)n dinucleotide repeat polymorphism with breast cancer risk and survival among Chinese women. Cancer Epidemiol Biomarkers Prev 2008,17(8):2117–2122.
Chinnery PF, Elliott HR, Patel S, Lambert C, Keers SM, Durham SE, McCarthy MI, Hitman GA, Hattersley AT, Walker M: Role of the mitochondrial DNA 16184–16193 poly-C tract in type 2 diabetes. Lancet 2005,366(9497):1650–1651.
Fernandez-Silva P, Enriquez JA, Montoya J: Replication and transcription of mammalian mitochondrial DNA. Exp Physiol 2003,88(1):41–56.
Rosson D, Keshgegian AA: Frequent mutations in the mitochondrial control region DNA in breast tissue. Cancer Lett 2004,215(1):89–94.
Ingman M, Gyllensten U: Analysis of the complete human mtDNA genome: methodology and inferences for human evolution. J Hered 2001,92(6):454–461.
Seligmann H, Krishnan NM, Rao BJ: Mitochondrial tRNA sequences as unusual replication origins: pathogenic implications for Homo sapiens. J Theor Biol 2006,243(3):375–385.
Malyarchuk BA, Rogozin IB, Berikov VB, Derenko MV: Analysis of phylogenetically reconstructed mutational spectra in human mitochondrial DNA control region. Hum Genet 2002,111(1):46–53.
Salas A, Yao YG, Macaulay V, Vega A, Carracedo A, Bandelt HJ: A critical reassessment of the role of mitochondria in tumorigenesis. PLoS Med 2005,2(11):e296.
Parson W, Bandelt HJ: Extended guidelines for mtDNA typing of population data in forensic science. Forensic Sci Int Genet 2007,1(1):13–19.
Bandelt HJ, Yao YG, Salas A, Kivisild T, Bravi CM: High penetrance of sequencing errors and interpretative shortcomings in mtDNA sequence analysis of LHON patients. Biochem Biophys Res Commun 2007,352(2):283–291.
Bandelt HJ, Quintana-Murci L, Salas A, Macaulay V: The fingerprint of phantom mutations in mitochondrial DNA data. Am J Hum Genet 2002,71(5):1150–1160.
Walt JM, Nicodemus KK, Martin ER, Scott WK, Nance MA, Watts RL, Hubble JP, Haines JL, Koller WC, Lyons K, et al.: Mitochondrial polymorphisms significantly reduce the risk of Parkinson disease. Am J Hum Genet 2003,72(4):804–811.
Zanssen S, Schon EA: Mitochondrial DNA mutations in cancer. PLoS Med 2005,2(11):e401.
Mancuso M, Conforti FL, Rocchi A, Tessitore A, Muglia M, Tedeschi G, Panza D, Monsurro M, Sola P, Mandrioli J, et al.: Could mitochondrial haplogroups play a role in sporadic amyotrophic lateral sclerosis? Neurosci Lett 2004,371(2–3):158–162.
Walt JM, Dementieva YA, Martin ER, Scott WK, Nicodemus KK, Kroner CC, Welsh-Bohmer KA, Saunders AM, Roses AD, Small GW, et al.: Analysis of European mitochondrial haplogroups with Alzheimer disease risk. Neurosci Lett 2004,365(1):28–32.
Bai RK, Leal SM, Covarrubias D, Liu A, Wong LJ: Mitochondrial genetic background modifies breast cancer risk. Cancer Res 2007,67(10):4687–4694.
Coskun PE, Ruiz-Pesini E, Wallace DC: Control region mtDNA variants: longevity, climatic adaptation, and a forensic conundrum. Proc Natl Acad Sci USA 2003,100(5):2174–2176.
Baudouin SV, Saunders D, Tiangyou W, Elson JL, Poynter J, Pyle A, Keers S, Turnbull DM, Howell N, Chinnery PF: Mitochondrial DNA and survival after sepsis: a prospective study. Lancet 2005,366(9503):2118–2121.
Munakata K, Tanaka M, Mori K, Washizuka S, Yoneda M, Tajima O, Akiyama T, Nanko S, Kunugi H, Tadokoro K, et al.: Mitochondrial DNA 3644T-->C mutation associated with bipolar disorder. Genomics 2004,84(6):1041–1050.
Brown MD, Starikovskaya E, Derbeneva O, Hosseini S, Allen JC, Mikhailovskaya IE, Sukernik RI, Wallace DC: The role of mtDNA background in disease expression: a new primary LHON mutation associated with Western Eurasian haplogroup J. Hum Genet 2002,110(2):130–138.
Jakupciak JP, Dakubo GD, Maragh S, Parr RL: Analysis of potential cancer biomarkers in mitochondrial DNA. Curr Opin Mol Ther 2006,8(6):500–506.
Raj GV, Moreno JG, Gomella LG: Utilization of polymerase chain reaction technology in the detection of solid tumors. Cancer 1998,82(8):1419–1442.
Mikuta JJ: Preoperative evaluation and staging of endometrial cancer. Cancer 1995,76(10 Suppl):2041–2043.
Tavassoli FA, Devilee P: World Health Organisation Classification of Tumours: Pathology and genetics of tumours of the breast and female genital organs. IARCPress-WHO 2003.
Malyarchuk BA, Grzybowski T, Derenko MV, Czarny J, Wozniak M, Miscicka-Sliwka D: Mitochondrial DNA variability in Poles and Russians. Ann Hum Genet 2002,66(Pt 4):261–283.
Wiesbauer M, Meierhofer D, Mayr JA, Sperl W, Paulweber B, Kofler B: Multiplex primer extension analysis for rapid detection of major European mitochondrial haplogroups. Electrophoresis 2006,27(19):3864–3868.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994,22(22):4673–4680.
Sheskin DJ: Handbook of parametric and nonparametric statistical procedures. 4 Edition Boca Raton, FL: Chapman & Hall 2007.
Fleiss JL, Levin B, Cho Paik M, Fleiss J: Statistical Methods for Rates & Proportions. 3 Edition San Francisco, CA: Wiley-Interscience Inc 2003.
Rosner B: Fundamentals of Biostatistics. 6 Edition Florence, KY: Duxbury Press 2005.
This work was supported by Polish Mitochondrial Network MitoNet.pl, Ministry of Science and Higher Education of The Republic of Poland Grant No. N N401 2327 33 to EB and AMC; Oligo.pl Minigrant G11 to EB, PG, AMC and KP, Polish Genetics Society Grant 2006/07 to EB, AMC, KP, and AK. AMC was supported by FEBS Collaborative Experimental Scholarship for Central & Eastern Europe, Fulbright Junior Research Grant and The Kosciuszko Foundation Scholarship.
The project realisation by its authors wouldn't have been possible without the support of Prof. Piotr Weglenski, PhD (Institute of Genetics and Biotechnology, Faculty of Biology, University of Warsaw). The authors would like to thank for Jerzy S. Czarnecki, PhD (University of Lodz, Lodz, Poland) and Przemyslaw Tomalski, PhD (Centre for Brain and Cognitive Development, School of Psychology, Birkbeck College, UK) for critical reading of the manuscript and fruitful discussions.
The authors declare that they have no competing interests.
AMC, AS, BK, TK, PG, EB - have made substantial contributions to the conception and design of the research, AMC, EB, AS, BK, PG - have been involved in drafting the manuscript, AMC, AK, KP, BM - performed the research, AS - collected the patients. All authors read and approved the final manuscript.