Single nucleotide polymorphisms in the mitochondrial displacement loop and outcome of esophageal squamous cell carcinoma
- Ruixing Zhang†1,
- Rui Wang†2,
- Fengbin Zhang1,
- Chensi Wu1,
- Haiyan Fan1,
- Yan Li3,
- Cuiju Wang4 and
- Zhanjun Guo1, 5Email author
© Zhang et al; licensee BioMed Central Ltd. 2010
Received: 25 October 2010
Accepted: 26 November 2010
Published: 26 November 2010
Accumulation of single nucleotide polymorphisms (SNPs) in the displacement loop (D-loop) of mitochondrial DNA (mtDNA) has been described for different types of cancers and might be associated with cancer risk and disease outcome. We used a population-based series of esophageal squamous cell carcinoma (ESCC) patients for investigating the prediction power of SNPs in mitochondrial D-loop.
The D-loop region of mtDNA was sequenced for 60 ESCC patients recorded in the Fourth Hospital of Hebei Medical University between 2003 and 2004. The 5 year survival curve were calculated with the Kaplan-Meier method and compared by the log-rank test at each SNP site, a multivariate survival analysis was also performed with the Cox proportional hazards method.
The SNP sites of nucleotides 16274G/A, 16278C/T and 16399A/G were identified for prediction of post-operational survival by the log-rank test. In an overall multivariate analysis, the 16278 and 16399 alleles were identified as independent predictors of ESCC outcome. The length of survival of patients with the minor allele 16278T genotype was significantly shorter than that of patients with 16278C at the 16278 site (relative risk, 3.001; 95% CI, 1.029 - 8.756; p = 0.044). The length of survival of patients with the minor allele 16399G genotype was significantly shorter than that of patients with the more frequent allele 16399A at the 16399 site in ESCC patients (relative risk, 3.483; 95% CI, 1.068 - 11.359; p = 0.039).
Genetic polymorphisms in the D-loop are independent prognostic markers for patients with ESCC. Accordingly, the analysis of genetic polymorphisms in the mitochondrial D-loop can help identify patient subgroups at high risk of a poor disease outcome.
Esophageal cancer is one of the commonest cancers in the population of northern central China with an age-standardized annual incidence rate > 125/100,000 . Cumulative mortality attributed to esophageal cancer is approximately 20% for women and 25% for men . The prognosis of esophageal cancer remains poor, despite improved diagnosis and therapeutic strategies, mostly because of its aggressive nature. The performance status, the TNM stage, and lymph node metastases seem to be the predictive factors of esophageal cancer; some molecular factors, such as p53 mutaion and NF-kappaB expression level, also show predictive power for esophageal cancer outcome .
The human mitochondrial genome is 16 kb in length and is a closed-circular duplex molecule that contains 37 genes, including 2 ribosomal RNAs and a complete set of 22 tRNAs . mtDNA is believed to be more susceptible to DNA damage and acquires mutations at a higher rate than nuclear DNA, because of the high levels of reactive oxygen species (ROS), the lack of protective histones, and limited capacity for DNA repair in the mitochondria [5, 6].
In cancers patients, sequence changes accumulated extensively in the mitochondrial D-loop region, which is important for regulating both replication and expression of the mitochondrial genome, because it contains the leading-strand origin of replication and the main promoter for transcription [7–10]. Only a few germline single nucleotide polymorphisms (SNPs) in the D-loop have been shown to be prognostic of cancer risk and outcome, but their predictive values have not been fully determined [11–14]. The D-loop contains a length of 1122 bps (nucleotide 16024-16569 and 1-576) refers to mitochondria database (http://www.mitomap.org). In this study, we sequenced a region of about 1 kb in the mitochondrial D-loop from 60 patients with esophageal squamous cell carcinoma (ESCC), the main pathological type of esophageal cancer in China, to assess the relationships between germline SNPs of the D-loop with cancer risk and gemline SNPs with disease outcome in ESCC patients.
Tissue specimens and DNA extraction
Blood samples were collected at the Fourth Hospital of Hebei University from 66 ESCC patients who underwent esophageal cancer resection in the Department of Thoracic Surgery between 2003 and 2004. The patients were selected when they received endoscopy examination and specimen were confirmed as ESCC by pathologist. All the patients comes from the Hebei Province of China a high risk area of ESCC. The tumor-free controls as determined per endoscopy, radiograph, and blood examination, were randomly selected from the same area. Both patients and controls contain 42 males and 24 females with the mean age of 59.78 ± 8.32 in ESCC patients and 60.84 ± 8.77 in controls. Genomic DNA was extracted immediately with a Wizard Genomic DNA extraction kit (Promega, Madison, WI) from blood samples. The study was approved by the Human Tissue Research Committee of the Fourth Hospital of Hebei Medical University. All patients provided written informed consent for the collection of samples and subsequent analysis.
PCR amplification and sequence analysis
The forward primer 5'-CCCCATGCTTACAAGCAAGT-3' (nucleotide 16190-16209) and reverse primer 5'-GCTTTGAGGAGGTAAGCTAC-3' (nucleotide 602-583) were used for amplification of a 982 bp product from mtDNA D-Loop region as described previously . PCR was performed according to the protocol of PCR Master Mix Kit (Promega, Madison, WI) and purified prior to sequencing. Cycle sequencing was carried out with the Dye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystem, Foster City, CA) and the products were then separated on the ABIPRISM Genetic Analyzer 3100 (Applied Biosystem). Polymorphisms were confirmed by repeated analyses from both strands. SNPs were identified directly from blood mitochondria.
The χ2 test was used to analyze dichotomous values, such as the presence or absence of an individual SNP between ESCC patients and healthy controls. The survival curve was calculated using the Kaplan-Meier method, and compared with the log-rank test. Multivariate survival analysis was performed using a Cox proportional hazards model. All of the statistical analysis was done with the SPSS 11.5 software package (SPSS Company, Chicago, IL). A p value of < 0.05 was considered statistically significant.
Univariate analysis of clinical characteristics associated with post-operational survival in ESCC patients
5 years survival rate (%)
SNPs in reference to GenBank accession AC_000021 were detected in 88 sites of the 982-bp mitochondria D-Loop region from blood samples [see Additional file 1], The sequence chromatograms show a clear single peak at each nucleotide position, indicating that mitochondria in ESCC individuals were homoplasmic. At first, we compared the distribution of germline SNPs at each site between ESCC and control patients to identify any link between an SNP and cancer risk; no association with ESCC cancer risk was detected in any SNP in the D-loop at p < 0.05 levels. We assessed the relationships between these SNPs and post-operational survival of these ESCC patients.
Multivariate analysis of prognostic factors associated with post-operational survival in ESCC patients with Cox proportional hazards model
Stage of tumor
Selected SNPs in the D-loop region have been examined for the ability to predict cancer risk in other types of tumour [11–14]. The present study has extended those analyses to determine the cancer risk and the post-operational survival-associated germline SNPs in a continuous sequence of mtDNA between nucleotides 16190 and 583 in ESCC patients. Three SNPs, 16274G/A, 16278C/T and 16399A/G, were identified for their association with post-operational survival at statistically significant levels by the log-rank test. Multivariate survival analysis identified 16278C/T and 16399A/G to be independent prediction markers for ESCC outcome.
We suggest for the first time that SNPs in the D-loop is a prognostic factor in ESCC patients. The relative risk (RR) of death in patients was significantly higher (16278C versus 16278T, RR, 3.001; 95% CI, 1.029 - 8.756; p = 0.044. 16399A versus 16399G, RR, 3.483; 95% CI, 1.068 - 11.359; p = 0.039). Nucleotides 16278 and 16399 are located in hypervariable segment 1 (HV1), which is associated with high rates of mutation , but the functional significance of these SNPs in HV1 is not known. Minor alleles of 16278T and 16399G are associated with dramatically shorter period of postoperative survival; the survival curve decreased rapidly in patients carrying these alleles (Figure 1). We compared the distribution frequency of these two SNPs between ESCC patients and normal controls; among 60 age-matched controls, only one carried the 16278T allele and none carried the 16399G allele. These two minor alleles appear to be correlated with ESCC risk, but a much larger sample size will be needed to obtain a statistically significant difference of distribution frequency for these alleles between ESCC patients and controls. It is very interesting if a nucleotide substitution can render an individual susceptible to a tumour but subsequently modulates prognosis of that tumour. The 5 year survival rate displays a trend for distinction on the basis of sex (35.1% male versus 47.8% female), age (46.5% older than 55 versus 23.5% younger than 55) and lymphatic metastasis (33.3% yes versus 41.7% no). The prediction value of these factors need to be investigated by further study.
The D-loop region of mtDNA is important for regulation of mitochondrial genome replication and expression. SNPs in this region might affect mtDNA replication and lead to alteration of the electron transport chain, which is responsible for the release of highly reactive oxygen species (ROS) and could contribute to nuclear genome damage as well as cancer initiation and promotion [17–19]. These three SNPs may altered transcription of mitochondrial genome, and that the production of ROS is enhanced when the mitochondrial transcription is altered , these ROS-mediated mechanism may accelerate the tumor development.
In conclusion, SNPs in the D-loop were found to be independent prognostic markers for ESCC outcome. The analysis of genetic polymorphisms in the D-loop might help to identify patient subgroups at high risk for a disease outcome, thereby helping to refine therapeutic decisions in ESCC cancers.
This work was supported by National Natural Science Foundation of PR China No. 30801384. The research was supported in part by Natural Science Foundation of Hebei Province No. C2008000958.
- Blot WJ, Li JY: Some considerations in the design of a nutrition intervention trial in Linxian, People's republic of China. Natl Cancer Inst Monogr. 1985, 69: 29-34.Google Scholar
- Abnet CC, Huppi K, Carrera A, Armistead D, McKenney K, Hu N, Tang ZZ, Taylor PR, Dawsey SM: Control region mutations and the 'common deletion' are frequent in the mitochondrial DNA of patients with esophageal squamous cell carcinoma. BMC Cancer. 2004, 4: 30-10.1186/1471-2407-4-30.View ArticleGoogle Scholar
- Blanchard P, Quero L, Hennequin C: Prognostic and predictive factor of oesophageal carcinoma. Bull Cancer. 2009, 96: 379-389.Google Scholar
- Shadel GS, Clayton DA: Mitochondrial DNA maintenance in vertebrates. Annu Rev Biochem. 1997, 66: 409-435. 10.1146/annurev.biochem.66.1.409.View ArticleGoogle Scholar
- DiMauro S, Schon EA: Mitochondrial DNA mutations in human disease. Am J Med Genet. 2001, 106: 18-26. 10.1002/ajmg.1392.View ArticleGoogle Scholar
- Beal MF: Mitochondia, free radicals, and neurodegeneration. Curr Opin Neurobiol. 1996, 6: 661-666. 10.1016/S0959-4388(96)80100-0.View ArticleGoogle Scholar
- Yoneyama H, Hara T, Kato Y, Yamori T, Matsuura ET, Koike K: Nucleotide sequence variation is frequently in the mitochondrial DNA displacement loop region of individual human tumor cells. Mol Cancer Res. 2005, 3: 14-20.Google Scholar
- Nashikawa M, Nishiguchi S, Shiomi S, Tamori A, Koh N, Takeda T, Kubo S, Hirohashi K, Kinoshita H, Sato E, Inoue M: Somatic mutation of mitochondrial DNA in cancerous and noncancerous liver tissue in individuals with hepatocellular carcinoma. Cancer Res. 2001, 61: 1843-1845.Google Scholar
- Sanchez-Cespedes M, Parrella P, Nomoto S, Cohen D, Xiao Y, Esteller M, Jeronimo C, Jordan RC, Nicol T, Koch WM, Schoenberg M, Mazzarelli P, Fazio VM, Sidransky D: Identification of a mononucleotide repaet as a major target for mitochondrial DNA alterations in human tumors. Cancer Res. 2001, 61: 7015-7019.Google Scholar
- Taanman JW: The mitochondrial genome: structure, transcription, translation and replication. Biochim Biophys Acta. 1999, 1410: 103-123. 10.1016/S0005-2728(98)00161-3.View ArticleGoogle Scholar
- Navaglia F, Basso D, Fogar P, Sperti C, Greco E, Zambon CF, Stranges A, Falda A, Pizzi S, Parenti A, Pedrazzoli S, Plebani M: Mitochondrial DNA D-loop in pancreatic cancer: somatic mutations are epiphenomena while the germline 16519 T variant worsens metabolism and outcome. Am J Clin Pathol. 2006, 126: 593-601. 10.1309/GQFCCJMH5KHNVX73.View ArticleGoogle Scholar
- Wang L, Bamlet WR, de Andrade M, Boardman LA, Cunningham JM, Thibodeau SN, Petersen GM: Mitochondrial genetic polymorphisms and pancreatic cancer risk. Cancer Epidemiol Biomarkers Prev. 2007, 16: 1455-1459. 10.1158/1055-9965.EPI-06-0367.View ArticleGoogle Scholar
- Wang L, McDonnell SK, Hebbring SJ, Cunningham JM, St Sauver J, Cerhan JR, Isaya G, Schaid DJ, Thibodeau SN: Polymorphisms in mitochondrial genes and prostate cancer risk. Cancer Epidemiol Biomarkers Prev. 2008, 17: 3558-3566. 10.1158/1055-9965.EPI-08-0434.View ArticleGoogle Scholar
- Bai RK, Leal SM, Covarrubias D, Liu A, Wong LJ: Mitochondrial genetic background modifies breast cancer risk. Cancer Res. 2007, 67: 4687-4694. 10.1158/0008-5472.CAN-06-3554.View ArticleGoogle Scholar
- 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. Mutant Res. 2004, 547: 71-78.View ArticleGoogle Scholar
- Stoneking M: Hypervariable sites in the mtDNA control region are mutational hotspots. Am J Hum Genet. 2000, 67: 1029-1032. 10.1086/303092.View ArticleGoogle Scholar
- Bandy B, Davision AJ: Mitochondrial mutations may increase oxidtaive stress: implications for carcinogenesis and aging?. Free Radic Biol Med. 1990, 8: 523-539. 10.1016/0891-5849(90)90152-9.View ArticleGoogle Scholar
- Gille JJ, Joenje H: Cell culture models for oxidative stress: Superoxide and hydrogen peroxidative versus normobaric heperoxia. Mutat Res. 1992, 275: 405-414.View ArticleGoogle Scholar
- Shigenaga MK, Hagen TM, Ames BN: Oxidative damage and mitochondrial decay in aging. Proc Natl Acad Sci USA. 1994, 91: 10771-10778. 10.1073/pnas.91.23.10771.View ArticleGoogle Scholar
- Dement GA, Maloney SC, Reeves R: Nuclear HMGA1 nonhistone chromatin proteins directly influence mitochondrial transcription, maintenance, and function. Exp Cell Res. 2007, 313: 77-87. 10.1016/j.yexcr.2006.09.014.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.