Open Access

Single nucleotide polymorphisms in genes MACC1, RAD18, MMP7 and SDF-1a as prognostic factors in resectable colorectal cancer


Cite

Introduction

Colorectal cancer (CRC) represents the third most common malignancy worldwide in men and second most common malignancy in women, accounting for approximately 10% of all tumour types worldwide and 8% of cancer related mortality.1 In Slovenia according to Cancer registry of Slovenia yearly reports, CRC is the second most common cancer in men and women. The incidence was steadily increasing in the last decades. From 2001 to 2011 it has risen by 35%, from 1110 in 2001 to 1709 in 2010 when it reached its peak. The incidence is now declining with 1530 in 2012. There is also a relative proximal shift of tumour location with the incidence of colon cancer increasing faster than the incidence of rectal cancer.25

The survival of CRC is improving as a consequence of screening programs, new chemotherapy regimens and targeted treatments as well as improved surgical treatment of metastatic disease. According to Surveillance, Epidemiology, and End Results (SEER) register and a worldwide CONCORD-2 study, the 5-year CRC survival in developed countries is more than 60%.6, 7

CRC is a heterogeneous disease and its cancero-genesis a multistep process. CRC develops in 75% sporadically because of mutations acquired during a person’s lifetime and in 25% as a combination of hereditary syndromes, a higher risk because of CRC familial burden without criteria for a hereditary syndrome or as a consequence of inflammatory bowel syndrome.8 It evolves through distinct genetic pathways: chromosomal instability, microsatellite instability (MSI-H) and the CpG island methylator phenotype (CIMP). Next to »classical adenoma-carcinoma« sequence proposed by Fearon and Voglestein there is also a newly recognized »serrated neoplasia pathway«, where CRC evolves through different precancerous lesions, as for instance serrated adenoma.9 The carcinogenesis might also differ in regard to CRC arising either in right or left hemicolon.10

The current treatment for resectable CRC of stage I, II and III is surgical resection. For patients of stage I, surgical resection is the only recommended treatment without adjuvant chemotherapy. For patients of stage III adjuvant chemotherapy is recommended in all patients. In contrast for patients of stage II adjuvant chemotherapy is not recommended for unselected patients, but for those with clinical or pathological risk factors.1113

Potential clinical and pathological risk factors for recurrence of stage II CRC have been investigated and incorporated in different guidelines, but a definite consensus has not yet been reached. According to European and American guidelines (The European Society for Medical Oncology (ESMO), The American Society of Clinical Oncology (ASCO) and The National Comprehensive Cancer Network (NCCN)), negative prognostic risk factors according to all three sets of guidelines are: T4 tumours, bowel perforation, extension of surgical lymphadenectomy, inadequate pathological sampling of lymph nodes and poor tumour differentiation grade. Further negative prognostic markers included in one or two sets of guidelines are: bowel obstruction, lymphovascular invasion and/or perineural invasion and indeterminate or positive margins. Consensus on them has not been reached yet. There is no clear message regarding adjuvant chemotherapy patient selection in stage II CRC.1113

Although stage I and early stage II CRC are prognostically very favourable, with a small burden of disease, a proportion of these tumours have certain characteristics, making them clinically more malignant and therefore predisposing them to disease recurrence or metachronous colon cancer.9 Up to 30% of patients with stage I and up to 50% of patients with stage II of CRC are going to relapse.14, 15 Considering these facts it is clear that classical TNM staging system has its limitations, so it is necessary to determine molecular or immuno-logical prognostic and predictive markers to implement in routine clinical practice.16, 17

Single nucleotide polymorphisms (SNPs) are molecular factors that might be useful as prognostic markers in CRC. Preliminary genome wide association (GWA) study in non-caucasian population has indicated a role of SNPs in resectable CRC.18 We hypothesized that SNPs participating in genetic risk for CRC and metastasis might prove as a prognostic factor in resectable CRC. In our study, we have selected SNPs with higher frequency in patients with either local lymph node involvement or systemic dissemination in genes participating in CRC carcinogenesis and disease dissemination: SDF-1a (stromal derived factor 1 alpha) located on chromosome 10, MMP7 located on chromosome 11, RAD18 located on chromosome 3, and MACC1 (metastasis associated in colorectal cancer 1) located on chromosome 7.1922 The aim of our study was to evaluate the role of SNPs in selected genes as prognostic markers in resectable CRC.

Patients and methods

We have conducted a study, regarding the role of selected SNPs in resectable CRC. In total, 163 consecutive patients treated surgically at University Medical Centre in Maribor in years 2007 and 2008 have been investigated. The inclusion criterion was colorectal adenocarcinoma of stages I, II or III. The exclusion criteria were: history of inflammatory bowel disease, preoperative chemoradiotherapy/radiotherapy, perioperative mortality within 30 days and confirmed familial CRC (FAP, HNPCC or other familial syndromes). Mean age of the patients at diagnosis was 67 years +/- 11.4 years (range 26–88 years). Clinical and pathohistological characteristics examined were: age, gender, TNM stage, differentiation grade, perineural invasion and lymphovascular invasion. Time to progression was defined as time from diagnosis to progression. Patients were followed-up on average for 69 months. Data regarding their vital status was acquired from Cancer registry of Slovenia. Clinical and pathohistological characteristics of patients and tumours are summarized in Table 1.

Clinical and pathohistological characteristics of patients included in study

Clinical and histopathological characteristicsCRC patients (N = 163)
Sex
Male/female, N(%)92/71 (56.4/43.6)
Age at diagnosis
Mean +/- SD67.28 +/- 11.44
Stage of disease, N (%)
I29 (17.8)
II81 (49.7)
III53 (32.5)
TNM staging
Tumour, N (%)
T18 (4.9)
T224 (14.7)
T3120 (73.6)
T411 (6.7)
Lymph nodes, N (%)
N0110 (67.5)
N138 (23.3)
N215 (9.2)
Vital status (5.10.2015), N (%)
dead65 (39.9) (42 due to CRC
progression, 23 other cause)
alive98 (60.1)
Disease progression, N (%)
yes46 (28.2) (42 dead, 4 alive)
no117 (71.8)
Disease progression according to stage, N (%)
stage I5 (17.2)
stage II18 (22.2)
stage III23 (43.4)
Clinical characteristics, N (%)
more than 12 lymph nodes resected70 (42.9)
adjuvant therapy41 (25.2)
Differentiation grade, N (%)
I66 (40.5)
II74 (45.4)
III21 (12.9)
no data2 (1.2)
Perineural invasion, N (%)
yes13 (8.0)
no150 (92.0)
Lymphovascular invasion, N (%)
yes27 (16.6)
no136 (83.4)

N = number of patients

DNA isolation and SNP genotyping

DNA was extracted from the formalin-fixed paraffin-embedded (FFPE) CRC tissues of 163 patients from the Department of Pathology, University Medical Centre in Maribor. FFPE tissues were prepared as follow: at macroscopic examination of resected specimen, the path olog ist sampled representative tissue samples (2 x 2 x 0.4 cm) from tumour, bowel wall outside the tumour and all lymph nodes from the resected pericolic or mesorectal fat. For determination of resection margin status, the representative tissue samples were taken also from the proximal and distal intestinal resection margin, from circumferential resection margin and dyed with the indian-ink. All representative tissue samples and all obtained lymph nodes were put into labeled histo-casettes and standardly processed in automated histoprocessors (dehydrated and paraffined). Paraffined tissue samples were embedded in paraffin blocks. For the extraction of DNA, 12μm thick tissue slices were cut with microtome from tumour tissue blocks.

DNA was isolated from FFPE tissues using BiOstic FFPE Tissue DNA Isolation Kit ®(MO BIO Laboratories, Inc.) according to manufacturer recommendations. SNPs of genes SDF-1a (rs1801157), MMP7 (rs11568818), RAD18 (rs373572) in MACC1 (rs1990172) were genotyped using polymerase chain reaction (PCR) followed by high resolution melting (HRM) or restriction fragment length polymorphism (RFLP) techniques. Forward and reverse primer sequences, size of product after PCR, primer concentrations, annealing temperatures and genotyping method are shown in Table 2. HRM genotyping was performed using real time PCR LC480 instrument (Roche, Germany). PCRHRM was carried out using LC480 HRM Master Mix (Roche, Germany). Conditions were as follow: initial denaturation at 95°C for 10 min, followed by 45 cycles of 95°C for 10 s, 57 or 60°C (primer pair dependent) for 15 s and 72°C for 10 s, followed by HRM step of 95°C for 1 min, 40°C for 1 min and 60 – 90°C at 0.02°C/s. PCR-RFLP conditions were as follow: initial denaturation at 95°C for 5 min, followed by 35 cycles of 95°C for 30 s, 60 or 63°C (primer pair dependent) for 30 s and 72°C for 30 s. After PCR, products were incubated with restriction enzymes shown in Table 2 at 37°C overnight. Digested products were resolved in 2% agarose gel.

Primer sequences used in PCR reaction, expected sizes of products, annealing temperatures, primer concentrations and genotyping method of selected SNPs; restriction enzymes and sizes of fragments after restriction for genotyping of SNPs rs1801157 (CXCL12) and rs1990172 (MACC1)

GeneSNP IDForward and reverse primerProduct size [bp]Annealing temperature [°C]Primer concentration [nM]Genotyping method
CXCL12rs1801157 A/GGTGGGATGGGATGGTGGAG CCTCAGCTCAGGGTAGCC10960650RFLP
MACC1rs1990172 G/TCAGGGAAAGAAATGGTTATTGCA GGAAAAGGAGGGAAGCATGTG11563300RFLP
MMP7rs11568818 A/GTGGAGTCAATTTATGCAGCAG CGAGGAAGTATTACATCGTTATTGG9357250HRM
RAD18rs373572 A/GTGTGATTAACCTAGTGGTTATTTTCTT GCATCCTAGTCTTCTCTATATTTTCG8560300HRM
GeneSNP IDRestriction enzymeSize of fragments after restriction [bp]
CXCL12rs1801157MspIAA: 109, AG: 109+62+47, GG: 62+47
MACC1rs1990172BseGIGG: 80+35, GT: 105+80+35, TT: 105

Our study was approved by the National Ethics committee of Slovenia (clinical trial registration number: 65/02/13) and is listed at University Medical Centre Maribor as research project: IRP-2014/01-21.

Statistical analysis

The clinical endpoint of our trial was evaluating the role of selected SNPs as prognostic factors by determining disease-free survival from the date of the surgery. Kaplan-Meier survival curves were constructed and compared using the log-rank test. Multivariate analysis was carried out using a Cox proportional hazard model. Group distribution for each clinicopathological characteristic was compared using two-tailed Fischer exact test. Data are expressed as the mean +/- standard deviation. Statistical significance was defined as p < 0.05. All analyses were performed using SPSS.

Results

Over an average follow-up period of 69 months, 65 deaths were recorded (39.9%). Out of those 65 patients, 42 died because of disease progression and 23 patients died of another cause.

Selected SNPs were genotyped in 163 CRC patients with well-defined clinical and histopathological characteristics. Genotype and allele frequencies were calculated for all patients and are shown in Table 3. When comparing selected SNPs with clinical and pathohistological characteristics of patients, higher frequency of genotype GG of MMP7 rs11568818 SNP was found in patients with T3/T4 staging (29.6%) compared to patients with T1/T2 staging (7.1%, p = 0.014), in patients with N1/N2 staging (36.2%) compared to N0 staging (19.8%, p = 0.041) and in patients with lymphovascular invasion (45.8%) compared to patients without lymphovascular invasion (21.0%, p = 0.018). For MACC1 rs1990172 SNP, we found higher frequency of genotype TT in patients with T3/T4 staging (62.5%) compared to patients with T1/T2 staging (36.0%, p = 0.024). All correlations between selected SNPs and clinical and histopathological characteristics are presented in Table 3.

Associations between selected SNPs and clinico-histopathological characteristics of patients

TNM stagingGrade of differentiationPerineural invasionLymphovascular invasion
Gene/ SNP IDFrequencyT1+T2T3+T4N0N1+N21+23NoYesNoYes
CXCL12 rs1801157AA4.760.05.97.30.04.85.35.20.05.80.0
AG29.2531.028.826.035.327.842.129.923.127.338.5
GG65.9969.065.366.764.767.552.664.976.966.961.5
Statistical analysis0.3460.0961.0001.0000.354p-valueAA vs AG+GG
1.2611.5730.9001.1021.228OR
1.159–1.3731.388–1.7830.102–7.9161.045–1.1621.135–1.32995% CI
0.8281.0000.3030.5430.656p-valueAA+AG vs GG
0.8140.9600.5561.8600.820OR
0.340–1.9480.471–1.9580.210–1.4710.488–7.0890.342–1.96795% CI
MACC1 rs1990172GG5.8412.04.57.82.15.25.35.67.77.10.0
GT36.5052.033.037.834.037.926.338.715.435.441.7
TT57.6636.062.554.463.856.968.455.676.957.558.3
Statistical analysis0.0240.3631.0000.5600.350p-valueGG vs GT+TT
2.9631.4770.9820.7181.229OR
1.202–7.3010.715–3.0500.112–8.6410.081–6.3381.131–1.33495% CI
0.2630.1600.4530.2371.000p-valueGG+GT vs TT
3.8802.9181.6412.6571.034OR
0.463–32.5160.649–13.1190.583–4.6200.697–10.1270.523–2.52595% CI
MMP7 rs11568818AA27.2728.627.028.125.527.926.326.733.327.725.0
AG47.5564.343.552.138.347.552.648.933.351.329.2
GG25.177.129.619.836.224.621.124.433.321.045.8
Statistical analysis1.0000.8431.0000.7361.000p-valueAA vs AG+GG
1.0841.1411.0820.7291.151OR
0.433–2.7130.517–2.5210.362–3.2340.207–2.5730.420–3.15295% CI
0.0140.0411.0000.4970.018p-valueAA+AG vs GG
5.4572.2960.8181.5473.182OR
1.226–24.2841.054–5.0020.252–2.6540.437–5.4791.273–7.95295% CI
RAD18 rs373572AA51.5946.752.848.658.050.455.050.761.551.950.0
AG38.2233.339.441.132.039.335.039.623.138.238.5
GG10.1920.07.910.310.010.410.09.715.49.911.5
Statistical analysis0.6850.3060.8120.5671.000p-valueAA vs AG+GG
0.7840.6850.8300.6431.079OR
0.353–1.7390.348–1.3480.323–2.1330.201–2.0580.465–2.50495% CI
0.0851.0001.0000.6240.731p-valueAA+AG vs GG
0.3420.9700.9601.6881.184OR
0.113–1.0300.318–2.9570.201–4.5800.339–8.3990.312–4.48895% CI

We have further evaluated the correlation of genotype frequencies in patients with progression compared to patients without disease progression. We found association between tumour TNM staging and SNP rs373572 in gene RAD18. Higher frequency of genotype GG was found in patients with T1/T2 staging with disease progression (60.0%) compared to patients with T1/T2 staging without disease progression (12.0%, p = 0.041). No statistically significant differences were discovered in rs1990172 in gene MACC1, rs1801157 in gene SDF-1a or rs11568818 in gene MMP7.

The results of survival analysis showed association with SNP rs1990172 in gene MACC1 and with SNP rs373572 in gene RAD18. We found worse disease-free survival (DFS) for patients with TT genotype of SNP rs1990172 in gene MACC1 compared to patients with GT or GG genotype (p = 0.029, Figure 1). One year, 3 years and 5 years DFS were in patients with TT genotype 94.8%, 67.3% and 62.9%, respectively compared to patients with GT/GG genotype, where 1 year, 3 years and 5 years DFS were 91.1%, 85.7% and 78.3%, respectively. Next, we found worse survival for patients with GG genotype compared to patients with AG or AA genotype for SNP rs373572 in gene RAD18 (p = 0.020, Figure 2). One year, 3 years and 5 years DFS were in patients with GG genotype 86.7%, 53.3% and 45.7%, respectively, compared to patients with AG/AA genotype, where 1 year, 3 years and 5 years DFS were 94.9%, 78.9% and 74.2%, respectively. The survival analysis for SNP rs11568818 in gene MMP7 and rs1801157 in gene SDF-1a did not show statistically significant differences.

Survival analysis curves for different genotype groups of CRC patients according to SNP rs1990172 in gene MACC1.

Survival analysis curves for different genotype groups of CRC patients according to SNP rs373572 in gene RAD18.

Discussion

Our study is the first report of the association between SNP rs373572 in RAD18 gene and SNP rs1990172 in MACC1 gene with DFS in resectable CRC. We also identified the association of SNP rs373572 in RAD18 gene in patients with stage I CRC and disease relapse.

Association between SNP rs373572 in RAD18 gene and DFS has been found, where 5-year DFS was significantly shorter for patients with GG genotype compared to patients with AA or AG genotype. Multivariate analysis showed that GG genotype could be independent negative prognostic factor, however association was at a borderline of statistical significance. In patients with stage I CRC, the GG genotype was more frequent in patients who relapsed, making it a potential negative prognostic factor for early stage CRC. RAD18 gene combines two distinct pathways maintaining genome stability.23 So far, association between SNP rs373572 and increased risk for CRC has been confirmed and significant association between SNP and clinicopathological features, specifically in differentiated grade and lymph node metastasis has been found in Japanese population.21 Association between rs373572 and CRC (colon and rectum cancer) has been also found in Han Chinese and it was shown that SNP is significantly related with increased risk of metastasis in CRC.24 Interestingly, this study showed in contrary better prognosis for patients with GG genotype. Further, statistically significant association has been found between MACC1 intron rs1990172 SNP genotype GT or TT and higher T stage of TNM. Patients with GT or GG genotype also had higher 5-years DFS compared to patients with TT genotype. Multivariate analysis showed that TT genotype was an independent negative prognostic factor. This may partially explain the poorer prognosis of patients with TT genotype. MACC1 is a regulator of the HGF/ Met signaling pathway which plays a key role in regulating many biological processes including cellular proliferation, cell metastasis, cell invasiveness, angiogenesis, epithelial mesenchymal transition, inducing increased invasiveness, tumorigenesis and also chemoresistance.25 The expression levels of MACC1 in colon cancer without distant metastases was significantly higher in primary tumours that later developed distant metastases, compared to those that did not metastasize within a 10-year-follow-up period.26 Several SNPs have been discovered in human MACC1 gene and Lang et al.27 conducted a trial researching six SNPs in MACC1 gene and they report a positive association of the MACC1 tagging SNP rs1990172 with reduced overall survival in patients with CRC. The study showed, in contrary to our results, better prognosis for patients with TT genotype. There are some differences in the design and clinical endpoint comparing both studies. Lang et al. also included patients with metastatic disease making the patient population more heterogenous. Their clinical endpoint was also overall survival (OS) not DFS. The authors concluded that rs1990172 was significantly associated with an increased risk for any death, not just CRC related. As we see in our study, out of 65 patients who died in course of follow up, 42 died of CRC recurrence and in 23 patients’cause of death was different. Considering this, DFS is more reliable to evaluate the role of a prognostic factor in adjuvant setting than OS.

Association between MACC1 rs1990172 SNP and other cancers has been found. For example, in patients with HER2-positive breast cancer, increased risk for progression or death for carriers of the G allele of SNP rs1990172 has been found.28

In our study we also found association between genotype GG of MMP7 gene promoter SNP rs11568818 and higher tumour (T) TNM stage, with regional lymph node dissemination and lymphovascular invasion. Our results suggest, that GG genotype of MMP7 rs11568818 SNP is negative prognostic biomarker for CRC. Consistent with our results, genotype GG of SNP rs11568818 in MMP7 has been recently associated with patients with CRC and aggressive course of the disease and with higher regional lymph node involvement.29 In this study relationship between GG genotype of MMP7 rs11568818 SNP and advanced tumour infiltration has also been reported. Higher prevalence of MMP7 rs11568818 GG genotype was also found among bladder cancer patients compared to controls, however, SNP was not associated with the tumour grade or stage.30 In meta-analysis of MMP2, MMP7 and MMP9 promoter polymorphisms, AG and GG genotype carriers of MMP7 rs11568818 SNP had an increased gastric cancer risk, suggesting it may play allele-specific role in cancer development.31 In breast cancer, statistically significant association with disease-free survival (DFS) was found for MMP7 rs11568818 SNP, where patients homozygous for G allele had significantly worse prognosis.32 In our study, statistically significant association with 5-year DFS for MMP7 rs11568818 SNP has not been found. Altogether, MMP7 seems very promising candidate for predicting tumour progression and metastasis also in CRC patients, particularly since it was also found, that the over-expression of MMP7 has considerable metastatic potential and correlates with unfavourable clinico-pathological characteristics.33

We also investigated SNP rs1801157 in CXCL12 (SDF-1a) gene. The CXCL12/CXCR4 axis promotes metastasis in numerous cancers. CXCL12 is being produced and released from tissues as liver or lung and trigger the migration of tumour cells expressing with CXCR4 receptor thereby promoting invasion, proliferation and survival under suboptimal condition.34 SNPs in CXCL12 gene have also been studied as a factor of increased likelihood developing cancer and increased likelihood of dissemination.35, 36 In a clinical trial conducted by Chang et al.37, GA/AA genotype of SNP rs1801157 was significantly higher in patients with lymph node metastasis among T3 tumours. In addition, an investigation of the relationship between CXCL12 genotypes and different clinico-pathological prognostic factors revealed a positive association between the GA/AA genotype and lympho-vascular invasion. Both of these results indicate a predisposition to worse prognosis. In our study we haven’t confirmed statistically significant association for SNP rs1801157 CXCL12 gene in 5-year DFS nor correlation between genotype distribution and clinicopathological characteristics. However, the distribution of genotypes was in our study similar compared to other studies.37 They however discovered higher AA genotype frequency in patients with T3 CRC in regional lymph node dissemination and lymphovascular invasion. Survival analysis showed worse DFS for the AA genotype in patients with lymph metastases.37

Stage III CRC patients are those who benefit from adjuvant chemotherapy resulting in increased DFS and OS at 6 and 10 years of follow up.38, 39 Patients with stage II disease with negative prognostic factors also benefit from adjuvant chemotherapy. There is however some controversy regarding unselected patients of stage II and also to some extent of stage I, because it is not always objectively possible to make a clear conclusion regarding negative regional lymph node status and the absence of negative pathohistological prognostic factors.40

TNM staging has in proportion of patients low prognostic value. Resection of appropriate number of lymph nodes is frequently impossible. In only about 50% of patients the required resection of 12 lymph nodes is achieved.16, 41, 42 Next to that, light microscopy has its limitations regarding sensitivity regarding detection of malignant cells in regional lymph nodes.43 Also only a proportion of a pathological specimen can be examined and this can lead to false negative reports regarding pathological risk factors.17 Inadequate resection of appropriate number of lymph nodes and limitations of light microscopy may lead to false downward stage migration. Lymphovascular and perineural invasion are negative prognostic factors, but they can also be underreported.44, 45 Considering these facts we see, that a substantial proportion of patients is undertreated.

A proportion of patients can have a more malignant phenotype irrespective of the TNM stage and known pathohistological prognostic factors, predisposing them to a more aggressive course of the disease. An intrinsic molecular characteristic, like nucleotide polymorphism, might prove extremely helpful in this regard. Polymorphisms of genes participating in carcinogenesis and disease dissemination thereby represent a potential new prognostic marker.

We have investigated patients from northeastern part of Slovenia almost exclusively of Slovenian origin. Our results regarding genotype distribution and connection with pathohistological characteristics are most consistent with study conducted by Dziki et al.29 on Polish population, which is as Slovenian also of Slavic origin. Our results compared to other studies are less consistent, but other studies were mainly performed in east Asia and Austria. The differences may be attributable to population genetic differences.

There is however also a question of appropriate study design, especially clinical endpoint. OS is of course the most important endpoint in oncology research, but it may be influenced by many treatment related factors beyond the point of disease relapse. DFS is therefore more appropriate clinical endpoint in resectable disease in research of prognostic markers. Although our results indicate the role of polymorphisms, further research is needed to validate our findings.

Screening of polymorphisms in selected genes of CRC patients in our study suggested that they may have a role as a prognostic factors in resectable CRC. In conclusion, the goal is to identify patients who are going to derive most clinical benefit, without facing unnecessary side effects.

eISSN:
1581-3207
Language:
English
Publication timeframe:
4 times per year
Journal Subjects:
Medicine, Clinical Medicine, Radiology, Internal Medicine, Haematology, Oncology