Open Access

Single nucleotide polymorphisms in the mitochondrial displacement loop and outcome of esophageal squamous cell carcinoma

  • Ruixing Zhang1,
  • Rui Wang2,
  • Fengbin Zhang1,
  • Chensi Wu1,
  • Haiyan Fan1,
  • Yan Li3,
  • Cuiju Wang4 and
  • Zhanjun Guo1, 5Email author
Contributed equally
Journal of Experimental & Clinical Cancer Research201029:155

https://doi.org/10.1186/1756-9966-29-155

Received: 25 October 2010

Accepted: 26 November 2010

Published: 26 November 2010

Abstract

Backgroud

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.

Methods

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.

Results

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).

Conclusion

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.

Background

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 [1]. Cumulative mortality attributed to esophageal cancer is approximately 20% for women and 25% for men [2]. 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 [3].

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 [4]. 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 [710]. 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 [1114]. 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.

Methods

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 [15]. 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.

Statistical analysis

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.

Results

A total of 66 patients were enrolled in this study. Six of these patients were lost to follow-up. A review was conducted every six months over a five-year period. Those patients lost to follow-up during this time period were as follows: 1 patient in Year 2; 1 patient in Year 3; 3 patients in Year 4; and, 1 patient in Year 5. Sixty patients shared the same performance status (ECOG Score: Zero). Those analyzed included: 12 patients at Stage I; 33 patients at Stage II, and 15 patients at Stage III. Thirty-six patients died during follow-up. None of these patients had received any adjuvant chemotherapy or radiation therapy after ESCC resection. Data for the 5 year follow-up period were analysed with clinical characteristics using the Kaplan-Meier method and were compared by the log-rank test. Sex, age and local lymphatic metastasis were not statistically significant predictors of the length of post-operational survival, but TNM stage was correlated with survival in these patients (Table 1). As expected, patients at different stages had different 5 year survival rates: stage I, 75%, stage II, 36.4% and stage III, 20%. The survival length distribution between any two stages was significantly different (p < 0.05) by the log-rank test. These data demonstrated that TNM stage is a good predictor of ESCC outcome.
Table 1

Univariate analysis of clinical characteristics associated with post-operational survival in ESCC patients

Characteristics

No. cases

5 years survival rate (%)

p value

Gender

   

0.129

 

Male

37

35.10

 
 

Female

23

47.80

 

Age (years)

  

0.282

 

≤ 55

17

23.50

 
 

> 55

43

46.50

 

TNM classificationa

  

0.012

 

I

12

75

 
 

II

33

36.40

 
 

III

15

20

 

Lymphatic metastases

  

0.418

 

Yes

12

33.30

 
 

No

48

41.70

 

aThe survival in each stage was compared as I versus II, I versus III and II versus III

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.

The relationship between mtDNA genotype and survival was compared subsequently, the ESCC patients were divided into two groups on the basis of their genotype at each SNP site, the post-operational survival curve was plotted using the Kaplan-Meier method for all ESCC patients at these sites. A dramatic difference in survival rate appeared at 16274, 16278 (refers to rs41458645 in NCBI SNP database, http://www.ncbi.nlm.nih.gov/snp/) and 16399 alleles by the log-rank test (Figure 1). The 3 SNPs were previously identified in mitochondria database (http://www.mitomap.org). The frequent allele 16274G, and the rare alleles 16278T and 16399G were associated with a shorter period of survival, with p = 0.0431, 0.0064 and 0.0028, respectively (Figure 1A, B and 1C). We performed multivariate analysis with Cox proportional hazards model including the factors of three SNPs and TNM stage. We coded G as 0 and A as 1 at 16274 position, C as 0 and T as 1 at 16278, A as 0 and G as 1 at 16399. We also coded I of TNM stage as 0, II as 1 and III as 2. As shown in Table 2, the 16278 and 16399 alleles were identified as independent predictors for ESCC outcome. The length of survival for patients with the rare allele 16278T genotype was significantly shorter than that for patients with the frequent allele 16278C (relative risk, 3.001; 95% CI, 1.029 - 8.756; p = 0.044) at the 16278 site. The same was seen for the rare allele 16399G genotype when compared with matched alleles 16399A at the 16399 site in ESCC patients (relative risk, 3.483; 95% CI, 1.068 - 11.359; p = 0.039) (Table 2). These data demonstrated the strong prediction power of 16278C/T and 16399A/G on outcome for ESCC patients.
Figure 1

Survival curve according to the nucleotide at position (A) 16274, (B) 16278 and (C) 16399 in D-loop of ESCC patients.

Table 2

Multivariate analysis of prognostic factors associated with post-operational survival in ESCC patients with Cox proportional hazards model

Factors

Relative risk

95% C.I.

p value

Stage of tumor

1.328

0.955-1.848

0.092

16274(G/A)

0

0

0.975

16278(C/T)

3.001

1.029-8.756

0.044

16399((A/G)

3.483

1.068-11.359

0.039

Discussion

Selected SNPs in the D-loop region have been examined for the ability to predict cancer risk in other types of tumour [1114]. 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 [16], 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 [1719]. These three SNPs may altered transcription of mitochondrial genome, and that the production of ROS is enhanced when the mitochondrial transcription is altered [20], 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.

Notes

Declarations

Acknowledgements

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.

Authors’ Affiliations

(1)
Department of Gastroenterology and Hepatology, The Fourth Hospital of Hebei Medical University
(2)
Department of Thoracic Sugery, The Fourth Hospital of Hebei Medical University
(3)
Department of Molecular Biology, The Fourth Hospital of Hebei Medical University
(4)
Department of Gynecology Ultrasound, The Fourth Hospital of Hebei Medical University
(5)
Hebei Key Lab of Laboratory Animal Science, Hebei Medical University

References

  1. 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
  2. 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
  3. Blanchard P, Quero L, Hennequin C: Prognostic and predictive factor of oesophageal carcinoma. Bull Cancer. 2009, 96: 379-389.Google Scholar
  4. 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
  5. 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
  6. Beal MF: Mitochondia, free radicals, and neurodegeneration. Curr Opin Neurobiol. 1996, 6: 661-666. 10.1016/S0959-4388(96)80100-0.View ArticleGoogle Scholar
  7. 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
  8. 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
  9. 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
  10. 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
  11. 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
  12. 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
  13. 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
  14. 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
  15. 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
  16. 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
  17. 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
  18. 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
  19. 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
  20. 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

Copyright

© Zhang et al; licensee BioMed Central Ltd. 2010

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.

Advertisement