Involvement of SLC39A6 in gastric adenocarcinoma and correlation of the SLC39A6 polymorphism rs1050631 with clinical outcomes after resection

doi:10.21203/rs.2.11309/v3

Download PDF

Research article

Involvement of SLC39A6 in gastric adenocarcinoma and correlation of the SLC39A6 polymorphism rs1050631 with clinical outcomes after resection

https://doi.org/10.21203/rs.2.11309/v3

This work is licensed under a CC BY 4.0 License

Journal Publication

published 08 Nov, 2019

Read the published version in BMC Cancer →

You are reading this older preprint version

Read the latest preprint version →

Background: The single-nucleotide polymorphism SLC39A6 rs1050631 is strongly implicated in esophageal squamous cell carcinoma, leading us to question whether it may also play a role in gastric adenocarcima (GA). Methods: We genotyped the SLC39A6 rs1050631 in 512 patients who underwent GA resection. All study subjects lived in an area of China with high GA incidence. Genotypes were examined for possible correlation with survival and recurrence. The potential involvement of SLC39A6 in gastric cancer was explored in clinical samples and cell culture studies. Results: Multivariable analysis showed that patients with the CT+TT genotype at SLC39A6 rs1050631 were at greater risk of recurrence (hazard ratio, HR 1.387, p = 0.004) and death (HR 1.429, p =0.002) than patients with CC genotype. Median recurrence-free and overall survival were significantly shorter in patients with the CT+TT genotype (20, 27 months) than in patients with the CC genotype (36, 43 months, p = 0.001, p < 0.001). Patients with the CT+TT genotype who were male or ≥60 years, or who had a tumor ≥ 5 cm or a moderately differentiated tumor were at significantly higher risk of recurrence and death. SLC39A6 was overexpressed in tissues from GA patients and in GA cell lines, and SLC39A6 knockdown in GA cell lines inhibited their proliferation, migration and invasion. Conclusion: SLC39A6 rs1050631 correlates with post-resection prognosis of GA patients and SLC39A6 may participate in GA onset or progression. Keywords: SLC39A6 rs1050631, Gastric adenocarcinoma, Prognostic biomarker, High-occurrence area, Ki67, TOPOII

Gastrointestinal Surgery

Epigenetics & Genomics

Health Economics & Outcomes Research

SLC39A6 rs1050631

Gastric adenocarcinoma

Prognostic biomarker

High-occurrence area

Ki67

TOPOII

Gastric cancer is one of the most common causes of cancer-related deaths worldwide [1]. Most gastric cancer cases occur in Asia, particularly in China [2,3]. The incidence of gastric cancer, its progression and patient prognosis differ across geographic regions and ethnic groups, and the reasons for these variations are poorly understood. A high-salt diet may stimulate gastric mucosa excessively, leading to chronic gastric inflammation and gastric carcinogenesis [4]. Consumption of N-nitroso-containing foods may contribute to gastric cancer, reflecting the fact that the N-nitroso group is a well-known carcinogen [5]. Genetic factors also appear to contribute to gastric cancer development and progression [6-8]. To clarify the factors contributing to gastric cancer, it may be beneficial to study populations living in areas with a high incidence of this disease, since such groups may have unique genetic backgrounds linked to the molecular mechanisms behind the illness. Therefore we examined patients who underwent surgical resection to treat gastric adenocarcinoma (GA) in the Chinese province of Fujian, one of the areas with the highest incidence of gastric cancer in China [9,10].

We obtained clinical samples from patients in Fujian province following resection to treat GA, and examined potential association between SLC39A6 r1050631 with clinical outcomes. We also examined the relationship between SLC39A6 expression and SLC39A6 r1050631. Potential effects of knocking down SLC39A6 expression were examined in representative GA cell lines.

Patients. This retrospective study included 512 Han Chinese patients living in Fujian, China. Briefly, we examined whether polymorphism in the gene encoding solute carrier family 39 (zinc transporter) member 6, often referred to as SLC39A6 or LIV-1, is associated with GA. This gene is known to promote the development and metastasis of several human cancers [11,12]. Studies involving patients from different parts of China have generated strong evidence linking SLC39A6 overexpression with risk of esophageal squamous cell carcinoma (ESCC) and poor survival [13,14], and linking the single-nucleotide polymorphism SLC39A6 rs1050631 with survival [14]. The esophagus is connected physically and functionally to the stomach, yet we are unaware of studies exploring a potential link between SLC39A6 rs1050631 and gastric cancer. Therefore we decided to focus on this polymorphism, although other polymorphisms may also be important in gastric cancer.

All patients were diagnosed with primary GA. Surgical resection of the primary gastric tumors was performed between July 2003 and December 2009 at 900 Hospital of the Joint Logistics Team (Fujian, China). Pathologists confirmed the diagnosis of GA following histopathological examination of the tumor tissues. All patients had complete medical records, including detailed clinical pathological features. Recurrence was defined based on our previously described method [15]. Survival was defined as the interval from the date of surgery to the date of death or the last follow-up (November 2014). Survival information was obtained primarily through telephone interview and the Social Security Death Index system. None of the patients included into this study had received preoperative chemotherapy. Of the 512 patients, 329 received postoperative chemotherapy with epirubicin, cisplatin, fluorouracil, or one or two of these three drugs plus the remaining one or two drugs. The following data were extracted from medical records in the hospital database: age, sex, tumor differentiation grade, tumor size, tumor-node-metastasis (TNM) stage, lymph node metastasis, distant metastasis, chemotherapy status, and other clinicopathological information. TNM staging and histologic classification were performed by experienced pathologists as described [16].

Immunohistochemical detection. SLC39A6 expression was examined in a subset of 198 randomly selected GA tissue blocks and 83 non-cancerous gastric tissues using standard immunohistochemical method. The anti-SLC39A6 antibody was from Abcam (Cambridge, MA). Immunostaining was assessed as described [16,17]. Tissues showing scores of ≥1+ for SLC39A6 staining were defined as positive; scores of ≥2+ were defined as high expression and <2+ as low expression.

SNP Selection and Genotyping. SLC39A6 rs1050631 was selected as the focus of the present study because of the strong evidence linking this gene to proliferation and invasion of ESCC cells, and this polymorphism to survival outcomes in ESCC patients, based on analysis of different groups of individuals from different parts of China [13,14]. The esophagus and stomach are physically and functionally connected in the digestive tract, and SLC39A6 rs1050631 has never been investigated in GA. Genomic DNA was extracted from 512 GA tissue samples using a QIAamp DNA FFPE Tissue Kit (Qiagen GmbH). The tissue samples had been formalin-fixed and parafﬁn-embedded immediately after surgical resection. Genotyping and analysis of the single-nucleotide polymorphism were performed as described [15-17]. The assay involved PCR to amplify the DNA, PCR product extension using a single primer, and product identification using MassARRAY SpectroCHIP and matrix-assisted laser desorption/ionization–time-of-flight mass spectrometry (Sequenom). Data were analyzed using TYPER software (Sequenom) [15-17].

Cell culture and siRNA transfection. One normal gastric cell line (GES-1) and four GA cell lines (BGC-823, SGC-7901, AGS and MGC-803) were used in this investigation. The lines GES-1, BGC-823, and AGS were obtained and cultured as described. The lines SGC-7901 and MGC-803 were purchased from the Cell Bank of the Shanghai Institute of Biochemistry and Cell Biology (Chinese Academy of Sciences, Shanghai, China) and cultured in the same way as BGC-823. Two siRNAs targeting different sites in the SLC39A6 gene were designed and synthesized by GenePharma (Shanghai, China). siRNA1 targeted nt 1213-1231 in the SLC39A6 cDNA (sense: 5′ –GCGAGGAAGUUAUCUGUAATT–3′; anti-sense: 5′–UUACAGAUAACUUCCUCGCTT–3′). siRNA2 targeted nt 658-676 in the SLC39A6 cDNA (sense: 5′ –GCAAUUUCCACACGGCAAUTT–3′; anti-sense: 5′–AUUGCCGUGUGGAAAUUGCTT–3′). A scrambled siRNA (sense: 5′–UUCUCCGAACGUGUCACGUTT–3′; anti-sense: 5′–ACGUGACACGUUCGGAGAATT–3′) was used as the negative control (NC). siRNAs were transfected into BGC-823 and SGC-7901 cells at a final concentration of 100 nM using Lipofectamine 2000.

Quantitative RT-PCR assay. Relative levels of SLC39A6 mRNA in GES-1 cells and four GA cell lines (AGS, BGC-823, SGC-7901 and MGC-803) were determined by quantitative RT-PCR as described for lncRNA HNF-AS1 [18], with the following primers: SLC39A6, 5′-GCCTGCAGTCTTGGAAGAAG-3′ and 5′-GCCGAGTGTATCGTGGAAAT-3′; GAPDH, 5′-GGGAGCCAAAAGGGTCA-3′ and 5′-GAGTCCTTCCACGATACCAA-3 ′.

Cell proliferation. BGC-823 and SGC-7901 cells were transfected for 48 h with anti-SLC39A6 siRNA, then 2,000 cells in 100 μl medium were seeded into each well of a 96-well plate. At 2 h (0 d), 1 d, 3 d and 5 d after seeding, CCK-8 solution (10 µl per well; Boster, Wuhan, China) was placed into each well and cells were incubated at 37 ℃ for 1 h. Absorbance at 450 nm was measured using a SpectraMax microplate reader (Molecular Devices, Thermo Fisher, Waltham, MA, USA). Culture dish containing no cells and only medium was used as the blank control. BGC-823 and SGC-7901 proliferation were plotted with time based on absorbance at 450 nm.

Cell migration and invasion. BGC-823 and SGC-7901 cell migration and invasion were examined almost as described [19], except that we performed measurements at 36 h and 48 h, respectively.

Statistical analysis. Statistical analysis was performed using SPSS 17.0 (IBM, Chicago, IL, USA). All statistical tests were two-sided, and p < 0.05 was considered significant. The chi-square test was used to identify associations of SLC39A6 rs1050631 genotypes with clinico-demographic characteristics, tumor pathology and outcomes (recurrence and survival). The variables used in uni- and multivariable Cox rgression were age, sex, tumor size, differentiation grade, gross findings, lymph node metastasis, distant metastasis, and chemotherapy regimen. Only sex, age and common variants or those with p < 0.2 in chi-square testing were adjusted in subsequent multivariable Cox models. Risks are presented as hazard ratios (HRs) and 95% conﬁdence intervals (CIs). Chi-square partitioning was used to assess the significance of relationships between SLC39A6 rs1050631 genotypes and survival or recurrence. Log-rank testing of Kaplan–Meier curves was used to assess associations of genotypes with recurrence-free and overall survival. Relationships among SLC39A6 rs1050631 genotypes and recurrence/survival were assessed using multivariable Cox models that adjusted for age, sex, and lymph node metastasis status. Stratification was also performed based on certain clinico-demographic characteristics in an effort to identify factors that strongly influenced recurrence or survival; Cox models in stratified analyses did not adjust for the variable that was being tested. Two-way ANOVA was used to evaluate differences in cell proliferation among different treatment groups. One-way ANOVA and Dunnett's multiple comparison test was used to assess the significance of differences in cell expression level, cell migration and invasion among different groups.

Associations between SLC39A6 rs1050631 genotype and clinico-demographic characteristics. Among the 512 tissue samples from patients resected for GA, 3.9% (20 cases), 27.1% (139 cases), and 69.0% (353 cases) had the TT, CT, or CC genotype, respectively (Table 1). Genotype frequencies did not differ significantly (based on the chi-squared test) with age (<60 vs. ≥60 years), degree of differentiation, gross findings (apophysis vs. invasion), Lauren classification (intestinal type, diffuse type, uncertain type), tumor location (upper third, middle third, lower third), tumor size (≥5 vs. <5 cm), distant metastasis status or chemotherapy status (Table 1). Potential association was observed between genotype and lymph node metastasis status (p <0.2). Genotype showed a significant association with post-resection survival (p = 0.002) and recurrence (p = 0.004) (Table 1).

Associations between SLC39A6 rs1050631 genotype and GA recurrence. Of all 512 patients, 334 (65.2%) experienced recurrence. Of the 353 patients with the CC genotype, 214 (60.6%) experienced recurrence (Tables 1 and 2). A significantly higher proportion of patients with the CT (104 of 139, 74.8%) or TT genotype experienced recurrence (16 of 20, 80.0%).

Univariable analysis showed risk of recurrence to be markedly higher in patients with the CT genotype (HR 1.813, p = 0.022) or the CT+TT genotype (HR 1.440, p = 0.001) (Table 2). Similar results were obtained with multivariable analysis: CT genotype, HR 1.379, p = 0.007; CT+TT genotype, HR 1.387, p = 0.004. In order to isolate factors that may influence post-resection outcomes in patients with GA, we stratified our patients by sex, age, tumor size, histologic grade,postoperative chemotherapy status and lymph node metastasis.Multivariable analysis showed that patients with the CT genotype were at increased recurrence risk if they were male (HR 1.564, p = 0.001), aged ≥60 years (HR 1.512, p = 0.009), with a tumor ≥5 cm (HR 1.807, p = 0.001) or with a moderately differentiated tumor (HR 1.830, p < 0.001) than those patients with CC genotype had (Table S1). Similar results were obtained for patients with the or CT+TT genotype: male, HR 1.545, p = 0.001; age, HR 1.529, p = 0.005; tumor ≥ 5 cm, HR 1.789, p = 0.001; and moderately differentiated tumor, HR 1.780, p < 0.001 (Table S1). In terms of the associations among the genotypes and recurrence risk involved chemotherapy status, Multivariable analysis exhibited that patients with CT, TT or CT+TT had higher recurrence risk only appeared in the group of patients without performing postoperative chemotherapy (Table S1).

Kaplan–Meier and log-rank analyses showed that median recurrence-free survival time was only 20 months in patients with the CT+TT genotype, significantly shorter than the 36 months in patients with the CC genotype (p = 0.001, Fig. 1A). Stratification analyses based on sex, age, differentiation grade, tumor size, histologic grade and postoperative chemotherapy status showed that patients with the CT+TT genotype who were male, aged ≥60 years or who had a tumor ≥ 5 cm or had a moderately differentiated tumor or no matter patients who whether performed postoperative chemotherapy had significantly shorter median recurrence-free survival time than patients with the CC genotype had (sex: 18 vs. 38 months, p < 0.001; age, 13 vs. 35 months, p = 0.001; tumor size, 9 vs. 35 months, p < 0.001; moderately differentiated tumor, 22 vs. 56 months, p < 0.001; had performed postoperative chemotherapy: 23 vs. 49months, p=0.031; or had not perform postoperative chemotherapy: 15 vs. 24 months, p=0.004. Fig. 1B, 1C and Fig. 2).

Associations between SLC39A6 rs1050631 genotype and overall survival. In the entire cohort of 512 patients, 330 (64.5%) died, similar to the proportion of patients with the CC genotype who died (210 of 353, 59.5%; p = 0.125). A significantly higher rate of death occurred among patients with the CT (104 of 139, 74.8%) or TT genotype (16 of 20, 80%), based on the chi-squared test and chi-square partitioning (Table 1).

Univariable Cox analysis revealed markedly increased risk of death in patients with the CT genotype (HR 1.440, p = 0.002) and also in patients with the CT+TT genotype (HR 1.492, p < 0.001). Similarly, multivariable Cox analysis revealed markedly increased risk of death in patients with the CT genotype (HR 1.416, p = 0.004) or the CT+TT genotype (HR 1.429, p=0.002), after adjusting for age, sex, and lymph node metastasis status (Table 2). Stratification by sex showed that risk of death was significantly higher among male patients with the CT genotype (HR 1.601, p = 0.001) or CT+TT genotype (HR 1.586, p < 0.001) than among the patients with CC genotype (Table S2). Similarly, stratification by age, tumor size, differentiation grade, postoperative chemotherapy status or lymph node metastasis revealed significantly increased risk of death among patients with the CT or CT+TT genotype when they were aged ≥60 years (HR 1.531, p = 0.007; HR 1.542, p = 0.004) or had a tumor ≥ 5 cm (HR 1.928, p < 0.001; HR 1.904, p < 0.001) or a moderately differentiated GA (HR 1.735, p = 0.001; HR 1.706, p = 0.001) (Table S2). In terms of the associations among the genotypes and death risk involved chemotherapy status, no matter patients had or had not performed postoperative chemotherapy, who carried CT or CT+TT genotypes showed higher death risk than those carried CC genotype (Table S2).

Patients with the CC genotype showed median overall survival time of 43 months, compared to 27 months for patients with the CT+TT genotype (p < 0.001; Fig. 3A). Stratification analyses showed that survival time was significantly shorter in patients with the CT+TT genotype if they had any of the following characteristics: male (26 vs. 45 months, p < 0.001), aged ≥60 years (20 vs. 42.0 months, p = 0.001), had a tumor ≥ 5 cm (16 vs. 42 months, p < 0.001), had moderately differentiated GA (31 vs. 59 months, p < 0.001) or no matter patients who whether performed postoperative chemotherapy: had performed postoperative chemotherapy: 31 vs. 57 months, p=0.05, which almost equal to have statistically significant difference ; or had not perform postoperative chemotherapy: 22 vs. 28 months, p=0.012.(Fig. 3B, 3C and Fig. 4).

SLC39A6 overexpression in GA. The results of the above mentioned experiments suggest an association between SLC39A6 rs1050631 and post-resection outcomes in patients with GA. This raised the question of whether SLC39A6 expression might be associated with GA. We used quantitative PCR and immunohistochemistry to assess expression levels in four GA cell lines (AGS, BGC-823, SGC-7901 and MGC-803) and our cohort of GA tissues. We found that the protein was overexpressed in all four cell lines relative to the normal gastric cell line GES-1, and that the protein was present (immunopositivity scores of ≥1+) in 150 of 198 (75.76%) GA tissues, compared to only 48 of 83 (57.83%) non-cancerous gastric tissues (p = 0.003; Fig.5A and 5B).

SLC39A6 knockdown inhibits proliferation, migration and invasion of GA cells. We designed two short interfering siRNAs to inhibit SLC39A6 expression in the GA cell lines BGC-823 and SGC-7901 (Fig. 6A and 6B), and this knockdown led to significantly less proliferation in both lines than in untransfected cells (Mock) or cells transfected with scrambled control siRNA (Fig. 6C and 6D). Knockdown also significantly reduced migration (Fig. 7A and 7B) and invasion by both cell lines (Fig. 7C and 7D). These results suggest that SLC39A6 may contribute to GA by functioning as a typical oncogene.

Association of different genotype at SLC39A6 rs1050631 with SLC39A6 expression in GA. Immunohistochemistry of 198 GA tissues from our patient cohort detected high SLC39A6 protein expression (≥2+) in 86 tissues (43.43%; Fig. 8); high expression was detected in 6 of 14 (42.86%) patients with TT, 50 of 92 (54.35%) patients with CT, 56 of 106 (52.83%) patients with CT+TT, and 31 of 92 (33.70%) patients with CC genotype. Though the high positive expression rate among the three different genotypes groups showed no statistically significantly different from that in the whole group (p>0.05) , the rate of high SLC39A6 expression was significantly greater among patients with CT genotype than among those with CC genotype (p=0.005) based on chi-square partitioning (Table 3). We speculate that the SLC39A6 rs1053631 genotype is associated with post-resection prognosis of patients with GA because the CT genotype leads to higher SLC39A6 expression.

Occurrence and prognosis of gastric cancer are influenced by genetics, lifestyle and the environment [20,21]. Single-nucleotide polymorphisms in certain genes may alter the levels or activity of proteins in a way that leads to cancer [22]. Here we analyzed the potential association of a particular polymorphism with GA and post-resection recurrence and death in a population living in a part of China with one of the highest incidences of gastric cancer in this country. This population may contain genetic elements that predispose them to GA, offering unique opportunities for studying genetic risk factors and prognostic markers, especially since most genetic risk studies are not conducted in high-incidence areas. Our data may also strengthen our understanding of risk factors in a high-risk environment.

Although several polymorphisms have been linked to gastric cancer, we focused on polymorphism in the SLC39A6 (LIV-1) gene because it is a downstream target of the STAT3 oncoprotein and promotes cancer progression by increasing cell growth and differentiation [23]. SLC39A6 can influence the epithelial–mesenchymal transition in pancreatic cancer cells, making their phenotype more aggressive [24]; and treatments that target SLC39A6 can reduce metastasis and predict prognosis of patients with hepatocellular carcinoma [25], and it can inhibit progression of metastatic breast cancer [26]. In particular, we focused on the SLC39A6 rs1050631 polymorphism. This polymorphism shows a strong association with death risk of ESCC in extensive studies involving more than a thousand of Chinese from at least two different parts of the country [14]. We are unaware of reports exploring a possible association between this polymorphism and GA.

Our study suggests that GA patients with the CT+TT genotype at rs1050631 are at significantly higher risk of recurrence and death than patients with CC genotype. Male gender, age older than 60 years, having a tumor ≥5 cm or having a moderately differentiated tumor are associated with worse prognosis in patients with the CT+TT genotype than in those with the CC genotype. It is still unknown the reason why SLC39A6 rs1050631 polymorphism had different effect on affecting the recurrence risk and death risk in GA patients with and without performing postoperative chemotherapy. Patients with CT or CT+TT genotypes were only associated with increased recurrence risk in the group of patients no carrying out postoperative chemotherapy. Whereas, patients with CT or CT+TT genotypes were associated with increased death risk both in the groups of patients carry and no carrying out postoperative chemotherapy. However, patients with CT+TT genotypes had both significantly decreased recurrence-free survival time and overall survival time than those with CC genotype. About these, it needs further investigation to reveal the reasons.

This observation of an association between SLC39A6 rs1050631 and post-resection outcomes in patients with GA prompted to ask whether SLC39A6 expression is associated with GA. Therefore we compared expression levels in a panel of GA cell lines as well as in our cohort of GA patients. The results show an association between SLC39A6 up-regulation and GA. We further found that the CT+TT genotype at SLC39A6 rs1050631, which was associated with worse clinical outcomes in our cohort. Consistent with a role of SLC39A6 in GA, we found that knocking down the protein in GA cell lines inhibited their proliferation, migration and invasion. These results illustrating SLC39A6 overexpression in GA tissues and GA cell lines as well as anti-cancer effects of SLC39A6 down-regulation mirror results published for SLC39A6 in ESCC [14]. These findings indicate that, as in ESCC and several other cancers, SLC39A6 appears to participate in GA, such that polymorphisms affecting its expression may be useful prognostic markers. Our results further suggest a model in which SLC39A6 rs1050631 alters the expression of SLC39A6 and thereby influences post-resection outcomes.

We identified several clinico-demographic characteristics associated with worse post-resection outcomes, including male gender, older age and having larger, moderately differentiated tumors or chemotherapy status.

In summary, our study provides preliminary evidence that SLC39A6 is involved in GA, and that genotype at SLC39A6 rs1050631 can predict post-resection prognosis of GA patients, at least in a population in an area with high GA incidence. Although we focused here on SLC39A6 rs1050631 because of the strong evidence in favor of a linkage with ESCC [14], future research is needed to explore other single-nucleotide polymorphisms that may be involved in gastric cancer. Due to this study has its limitation: all the clinical samples were obtained from single center in Fujian Province, further studies are also needed to examine whether our findings in this Chinese population can be generalized to other ethnic groups.

GA:gastric adenocarcima; ESCC:esophageal squamous cell carcinoma TNM:tumor-node-metastasis ; NC:negative control; HR:hazard ratios;CI:conﬁdence intervals.

Ethics approval and consent to participate

The present study was approved by the Ethics Committee of the Fuzhou General Hospital(Dongfang Hospital), and informed written consent was obtained from all patients. Tumor tissues of GA patients were obtained based on the Declaration of Helsinki

Consent to publish

Not applicable.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare that they have no competing interests.

Funding

This study was supported by the Innovation Team Foundation of Fuzhou General Hospital(grant no. 2014CXTD04), the National Natural Science Foundation of Fujian Province, China (grant no. 2017J01215) and the Large Data Application and Demonstration Project of Fujian Development and Reform Commission.The funding source had no role in the study design, data collection or analyses, the decision to publish, or the preparation of this manuscript.

Authors' Contributions

JG and WJR performed the experiments including immunohistochemical staining, cell culture, cell transfection, and cell function assays. JG,WJR and CHX conducted the statistical analysis. PCW,HW,XHFand ZZZ collected the data. YD and LYL analyzed and interpreted the data. QJH wrote the manuscript. LW and SQC revised the manuscript. LYL and SQC designed the project. All authors have read and approved the fnal manuscript.

Acknowledgements

Not applicable.

Authors information

¹Department of General Surgery, Dongfang Hospital (900 Hospital of the Joint Logistics Team), 156 North Xi-er Huan Road, Fuzhou, Fujian 350025, China

²China Clinical Institute of Fuzhou General Hospital (900 Hospital of the Joint Logistics Team), Fujian Medical University, 156 North Xi-er Huan Road, Fuzhou, Fujian 350025, China

³Department of Experimental Medicine, Dongfang Hospital (900 Hospital of the Joint Logistics Team), 156 North Xi-er Huan Road, Fuzhou, Fujian 350025, China

⁴Fujian Meiya Aijiankang Health Management Co., Ltd. 4602#, Building 1, Shimao International Center, 108 Guangda Road, Fuzhou, Fujian 350025, China

⁵Union Medical College, Fujian Medical University, 29 XinQuan Road, Fuzhou, Fujian 350025, China

* These authors contributed equally to the article and should be considered co-first authors

Correspondence to: Professor Shaoquan Chen and Mr. Liying Lin, Department of General Surgery, Dongfang Hospital (900 Hospital of the Joint Logistics Team), 156 North Xi-er Huan Road, Fuzhou, Fujian 350025, China.

E-mail: [email protected] for Shaoquan Chen and [email protected] for Liying Lin.

Xie XS, Lin JX, Li P, Xie JW, Wang JB, Lu J, Chen QY, Cao LL, Lin M, Tu RH, et al: A risk prediction system of postoperative hemorrhage following laparoscopy-assisted radical gastrectomy with D2 lymphadenectomy for primary gastric cancer. Oncotarget 8(46):81511-81519, 2017.
Pan QX, Su ZJ, Zhang JH, Wang CR and Ke SY: A comparison of the prognostic value of preoperative inflammation-based scores and TNM stage in patients with gastric cancer. Onco Targets Ther 8:1375-1385, 2015.
Zhuo C, Ying M, Lin R, Wu X, Guan S and Yang C: Negative lymph node count is a significant prognostic factor in patient with stage IV gastric cancer after palliative gastrectomy. Oncotarget 8:71197-71205, 2017.
Lin SH, Li YH, Leung K, Huang CY and Wang XR: Salt processed food and gastric cancer in a Chinese population. Asian Pac J Cancer Prev 15: 5293-5298, 2014.
Yang Q, Xu E, Dai J, Wu J, Zhang S, Peng B and Jiang Y: miR-21 regulates N-methyl-N-nitro-N'-nitrosoguanidine-induced gastric tumorigenesis by targeting FASLG and BTG2. Toxicol Lett 228: 147-156, 2014.
Zhu H, Li X, Zhang X, Chen D, Li D, Ren J, Gu H, Shu Y and Wang D: Polymorphisms in mismatch repair genes are associated with risk and microsatellite instability of gastric cancer, and interact with life exposures. Gene 579(1): 52-57, 2016.
Li XW, Shi BY, Yang QL, Wu J, Wu HM, Wang YF, Wu ZJ, Fan YM and Wang YP: Epigenetic regulation of CDH1 exon 8 alternative splicing in gastric cancer. BMC Cancer 15: 954, 2015.
Li-Chang HH, Kasaian K, Ng Y, Lum A, Kong E, Lim H, Jones SJ, Huntsman DG, Schaeffer DF and Yip S: Retrospective review using targeted deep sequencing reveals mutational differences between gastroesophageal junction and gastric carcinomas. BMC Cancer 15: 32, 2015.
Ang TL, FocK KM: Clinical epidemiology of gastric cancer. Singapore Med J 55(12): 621–628, 2014.
Cai L, Yu SZ, Ye WM and Yi YN: Fish sauce and gastric cancer: an ecological study in Fujian Province, China. World J Gastroenterol6(5): 671–675, 2000.
Lue HW, Yang X, Wang R, Qian W, Xu RZ, Lyles R, Osunkoya AO, Zhou BP, Vessella RL, Zayzafoon M, et al: LIV-1 promotes prostate cancer epithelial-to-mesenchymal transition and metastasis through HB-EGF shedding and EGFR-mediated ERK signaling. PLoS One 6: e27720, 2011.
Zhao L, Chen W, Taylor KM, Cai B and Li X: LIV-1 suppression inhibits HeLa cell invasion by targeting ERK1/2-Snail/Slug pathway. Biochem Biophys Res Commun 363: 82-88, 2007.
Cui XB, Shen YY, Jin TT, Li S, Li TT, Zhang SM, Peng H, Liu CX, Li SG, Yang L, et al: SLC39A6: a potential target for diagnosis and therapy of esophageal carcinoma. J Transl Med 13: 321, 2015.
Wu C, Li D, Jia W, Hu Z, Zhou Y, Yu D, Tong T, Wang M, Lin D, Qiao Y, et al: Genome-wide association study identifies common variants in SLC39A6 associated with length of survival in esophageal squamous-cell carcinoma. Nat Genet 45(6): 632-638, 2013.
Wang X, Lin Y, Lan F, Yu Y, Ouyang X, Wang X, Huang Q, Wang L, Tan J and Zheng F: A GG allele of 3'-side AKT1 SNP is associated with decreased AKT1 activation and better prognosis of gastric cancer. J Cancer Res Clin Oncol 140: 1399-1411, 2014.
Liu W, Zhong S, Chen J, Yu Y: HER-2/neu overexpression is an independent prognostic factor for intestinal-type and early-stage gastric cancer patients. J Clin Gastroenterol 46: e31-37, 2012.
Wang X, Lin Y, Lan F, Yu Y, Ouyang X, Liu W, Xie F, Wang X and Huang Q: BAX and CDKN1A polymorphisms correlated with clinical outcomes of gastric cancer patients treated with postoperative chemotherapy. Med Oncol 31: 249, 2014.
Dang Y, Lan F, Ouyang X, Wang K, Lin Y, Yu Y, Wang L, Wang Y and Huang Q: Expression and clinical significance of long non-coding RNA HNF1A-AS1 in human gastric cancer. World J Surg Oncol 13: 302, 2015.
Ma L, Lan F, Zheng Z, Xie F, Wang L, Liu W, Han J, Zheng F, Xie Y and Huang Q: Epidermal growth factor (EGF) and interleukin (IL)-1β synergistically promote ERK1/2-mediated invasive breast ductal cancer cell migration and invasion. Mol Cancer 11: 79, 2012.
Lynch HT, Grady W, Suriano G, Huntsman D: Gastric cancer: new genetic developments. J Surg Oncol 90(3): 114-133, 2005.
Yang L: Incidence and mortality of gastric cancer in China.World J Gastroenterol 12: 17-20, 2006.
Belyavskaya VA, Prudnikova TY, Domanitskaya NV, Litviakov NV, Maksimov VN, Cherdyntseva NV and Grigorieva EV: GLCE rs3865014 (Val597Ile) polymorphism is associated with breast cancer susceptibility and triple-negative breast cancer in Siberian population. Gene 628: 224-229, 2017.
Hogstrand C, Kille P, Ackland ML, Hiscox S and Taylor KM: A mechanism for epithelial-mesenchymal transition and anoikis resistance in breast cancer triggered by zinc channel ZIP6 and STAT3 (signal transducer and activator of transcription 3). Biochem J 455(2): 229-237, 2013.
Unno J, Satoh K, Hirota M, Kanno A, Hamada S, Ito H, Masamune A, Tsukamoto N, Motoi F, Egawa S, et al: LIV-1 enhances the aggressive phenotype through the induction of epithelial to mesenchymal transition in human pancreatic carcinoma cells. Int J Oncol 35(4): 813-821, 2009.
Lian J, Jing Y, Dong Q, Huan L, Chen D, Bao C, Wang Q, Zhao F, Li J, Yao M, et al: miR-192, a prognostic indicator, targets the SLC39A6/SNAIL pathway to reduce tumor metastasis in human hepatocellular carcinoma. Oncotarget 7: 2672-2683, 2016.
Sussman D, Smith LM, Anderson ME, Duniho S, Hunter JH, Kostner H, Miyamoto JB, Nesterova A, Westendorf L, Van Epps HA, et al: SGN-LIV1A: a novel antibody-drug conjugate targeting LIV-1 for the treatment of metastatic breast cancer. Mol Cancer Ther 13: 2991-3000, 2014.

Table 1. Patient characteristics and distribution of SLC39A6 rs1050631 genotypes

Variables	Patients N (%)	SLC39A6 rs1050631 genotype (N)					P^*	P_dominant^a	P_additive^b
Variables	Patients N (%)	TT(%)	CT(%)		CC(%)		P^*	P_dominant^a	P_additive^b
Total	512	20(3.9)	139(27.1)		353(69.0)
Age						0.26		0.25	0.14
<60	235(45.9)	6(1.2)	61(11.9)	168(32.8)
≥60	277(54.1)	14(2.7)	78(15.2)	185(36.1)
Gender						0.38		0.18	0.25
Male	383(74.8)	15(2.9)	110(21.5)	258(50.4)
Female	129(25.2)	5(1.0)	29(5.7)	95(18.5)
Grade of differentiation						0.85		0.93	0.86
Well differentiated	21(4.1)	1(0.2)	6(1.2)	14(2.7)
Moderately differentiated	263(51.4)	9(1.8)	75(14.6)	179(35.0)
Poorly differentiated	228(44.5)	10(2.0)	58(11.3)	160(31.2)
Gross findings						0.79		0.96	0.81
Apophysis	35(6.8)	2(0.4)	9(1.8)	24(4.7)
Invasion	477(93.2)	18(3.5)	130(25.4)	329(64.3)
Lauren classification						0.65		0.45	0.22
Intestinal type	420(82.1)	18(3.5)	119(23.2)	283(55.3)
Diffuse type	81(15.8)	2(0.4)	18(3.5)	61((11.9))
Uncertain type	11(2.1)	0(0)	2(0.4)	9(1.8)
Tumor location						0.32		0.38	0.42
Upper third	179(35.0)	4(0.8)	52(10.2)	123(24.0)
Middle third	78(15.2)	4(0.8)	16(3.1)	58(11.3)
Lower third	252(49.2)	12(2.3)	69(13.5)	171(33.4)
Total	3(0.6)	0(0)	2(0.4)	1(0.2)
Tumor size						0.99		0.88	0.92
≥5cm	210(41.0)	8(1.6)	58(11.3)	144(28.1)
<5cm	302(59.0)	12(2.3)	81(15.8)	209(40.8)
Survival						0.002		<0.001	<0.001
Alive	182(35.5)	4(0.4)	35(6.8)	143(27.9)
Dead	330(64.5)	16(3.1)	104(20.3)	210(41.0)
Recurrence						0.004		<0.001	<0.001
Yes	334(65.2)	16(3.1)	104(20.3)	214(41.8)
No	178(34.8)	4(0.8)	35(6.8)	139(27.1)
Distant metastasis						0.43		0.87	0.81
Yes	63(12.3)	4(0.8)	15(2.9)	44(8.6)
No	449(87.7)	16(3.1)	124(24.2)	309(60.4)
Lymph node metastasis						0.12		0.21	0.08
Yes	398(77.7)	19(3.7)	110(21.5)	269(52.5)
No	114(22.3)	1(0.2)	29(5.7)	84(16.4)
Chemotherapy
Yes	329(64.3)	13(2.5)	93(18.2)	223(43.6)		0.61		0.78	0.69
No	183(35.7)	7(1.4)	46(9.0)	130(25.4)

^* The χ² test was used to calculate p values. Fisher’s exact test was used in the analysis of contingency tables when sample size was smaller than 5.

^aLogistic regression analysis for the dominant model (TT+CT vs CC) was used to calculate the p value.

^bLogistic regression analysis for the log-additive model was used to calculate the p value by comparing the effect of each additional variant allele.

Table 2.Associations between SLC39A6 rs1050631 genotypes and recurrence and survival

	SNP		Genotype		Recurrence		Unadjusted					Adjusted^a
					No(n)/Yes(n)		HR		95%CI	p		HR			95%CI			p
	Total (n=512)				178/334
	rs1050631		CC		139/214		1(reference)					1(reference)
			CT		35/104		1.813		1.09-3.01	0.022		1.379		1.090-1.744			0.007
			TT		4/16		1.400		1.11-1.77	0.005		1.441		0.864-2.403			0.161
			CT+TT		39/120		1.440		1.15-1.80	0.001		1.387		1.108-1.735			0.004
	SNP		Genotype		Survival		Unadjusted					Adjusteda
					Ali(n)/Dea(n)		HR		95%CI	p		HR		95%CI			p
	Total(n=512)				182/330
rs1050631		CC		143/210		1(reference)							1(reference)
		CT		35/104		1.440		1.138-1.823			0.002		1.416			1.119-1.792			0.004
		TT		4/16		1.947		1.170-3.239			0.010		1.523			0.912-2.542			0.107
		CT+TT		39/120		1.492		1.192-1.868			<0.001		1.429			1.141-1.790			0.002

The bold values indicate p value < 0.05.

^aData were calculated using multivariable Cox models that adjusted for gender, age and lymph node metastasis.

Abbreviations: HR, hazard ratio; CI, Confidence Interval.

Table 3 Chi-square partitioning to test the association between SLC39A6 rs1050631 and SLC39A6 expression detected by immunohistochemistry

SNP	Genotype	expression of SLC39A6		Case-total^a	Case-case^b
SNP	Genotype	Low(n)/ High(n)	High expression rate (%)	p<0.00833	p<0.0125
Total(n=198)		112/86	43.4
rs1050631	CC	61/31	33.7	0.116	1vs2:0.005
	CT	42/50	54.3	0.083	2vs3:0.422
	TT	8/6	42.9	0.966	1vs3:0.712

^aCase–total: p value for the comparison between the group with the indicated genotype and the total group of 198 samples;

^bCase–case: p value for the pairwise comparisons between groups with different genotypes.

Additional file 1:Table S1.Associations between SLC39A6 rs1050631 genotypes and recurrent after stratification by sex, age, tumor size, differentiation grade,chemotherapy status and lymph node metastasis.Table S2. Associations between SLC39A6 rs1050631 genotypes with survival, after stratification by sex, age, tumor size, differentiation grade,chemotherapy status and lymph node metastasis.

Additionaltable.docx

Download PDF

Journal Publication

published 08 Nov, 2019

Read the published version in BMC Cancer →

Editorial decision: Minor revision
24 Sep, 2019

You are reading this older preprint version

Read the latest preprint version →

Involvement of SLC39A6 in gastric adenocarcinoma and correlation of the SLC39A6 polymorphism rs1050631 with clinical outcomes after resection

Status:

Journal Publication

Version 3

Abstract

Figures

Introduction

Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

References

Tables

Additional File

Supplementary Files

Status:

Journal Publication

Version 3