Polymorphisms in the selenoprotein S and 15-kDa selenoprotein genes are associated with altered susceptibility to colorectal cancer

Abstract

Selenium (Se), a dietary trace metal essential for human health, is incorporated into ~25 selenoproteins including selenoprotein S (SelS) and the 15-kDa selenoprotein (Sep15) both of which have functions in the endoplasmic reticulum protein unfolding response. The aim of this study was to investigate whether genetic variants in such selenoprotein genes are associated with altered risk of colorectal cancer (CRC). A Korean population of 827 patients with CRC and 733 healthy controls was genotyped for 7 SNPs in selenoprotein genes and one SNP in the gene encoding manganese superoxide dismutase using Sequenom technology. Multivariate logistic regression analysis showed that after adjustment for lifestyle factors three SNP variants were associated with altered disease risk. There was a mean odds ratio of 2.25 [95% CI 1.13,4.48] in females homozygous TT for rs34713741 in SELS with the T variant being associated with higher risk of rectal cancer, and odds ratios of 2.47 and 2.51, respectively, for rs5845 and rs5859 in SEP15 with the minor A and T alleles being associated with increased risk of male rectal cancer. The data indicate that the minor alleles for rs5845, rs5859 and rs34713741 are associated with increased rectal cancer risk and that the effects of the three SNPs are dependent on gender. The results highlight potential links between Se, the function of two selenoproteins involved in the protein unfolding response and CRC risk. Further studies are required to investigate whether the effects of the variants on CRC risk are also modulated by dietary Se intake.

Introduction

Selenium (Se) is a dietary micronutrient that is essential for human health [23, 29]. Cancer mortality is inversely correlated with estimated Se intake [8, 29]. Low Se intake or low plasma Se concentrations are associated with increased risk of colorectal cancer (CRC) [7, 25] whereas in contrast, higher Se status is associated with lower risk of a recurrence of colonic tumours [15]. In the United States, daily supplementation with 200 μg Se has been found to lead to reduced mortality from CRC [6].

The biological functions of Se are thought to be brought about through its presence as the amino acid selenocysteine in ~25 selenoproteins [2, 19]. Selenocysteine is incorporated into selenoproteins during their synthesis with a specific stem-loop structure in the 3′untranslated region (3′UTR) of the mRNAs required to recode a UGA codon from “stop” to selenocysteine [2]. These selenoproteins include the family of glutathione peroxidises [4], selenoprotein P (SePP) which has a Se transport role [5], selenoprotein S which is an endoplasmic reticulum protein involved in removing unfolded proteins and the 15-kDa selenoprotein which is also an endoplasmic reticulum protein involved in the unfolded protein response [11, 18]. Because the selenoproteins have important roles in cell stress responses and redox control [2, 4], it is possible that genetic variation due to single nucleotide polymorphisms (SNP) in the genes encoding the selenoproteins may influence susceptibility to cancer and Se requirements for optimal health.

Several SNPs in selenoprotein genes have been reported to have functional consequences. For example, rs1050450 causes a Pro-Leu amino acid change in glutathione peroxidise 1 and affects protein function [12], rs713041 causes a T-C substitution in a region of the glutathione peroxidise 4 (GPX4) gene corresponding to the 3′UTR of the mRNA and alters the protein binding to the 3′UTR and reporter gene activity [3, 21], rs5859 and rs5845 in the 3′UTR of the 15-kDa selenoprotein affect reporter gene activity [13] and rs7579 and rs3877899 in the SEPP1 gene affect biomarkers of selenoprotein status in vivo [20]. Variants in the promoter of SELS are functional in that they modulate SelS expression and plasma levels of inflammatory cytokines [10].

Only a very limited number of studies have explored whether allelic variants of SNPs in selenoprotein genes are associated with altered disease susceptibility, and only three studies have investigated association with CRC or adenoma risk. Variants of a SNP in the promoter of the selenoprotein P gene (SEPP1) were not associated with altered risk of CRC [1], but a more recent study indicated that three variants in SEPP1 and one in the thioredoxin reductase 1 gene TXNRD1 showed an association with risk of adenoma [22]. In addition, a relatively small case–control study in the United Kingdom showed an association of the C variant of rs713041 in GPX4 with increased risk of CRC [3]. The aim of the present work was to extend these studies by investigating the association between a series of SNPs and CRC risk, focusing on SNPs present in selenoprotein genes and known to be functional. To do this, we analysed DNA from a case–control study of a Korean population.

Subjects and methods

Subjects

The cases comprised of a consecutive series of patients with histologically confirmed incident colorectal adenocarcinomas, who were admitted to two university hospitals in Seoul, Korea between 1998 and 2004. At diagnosis, distinction was made between colon and rectal cancers. Tumours in the sigmoid colon were classified as colonic and those at rectosigmoid junction as rectal. A total of 827 patients were recruited. Seven hundred and thirty-three controls were selected from patients admitted to orthopaedics, general surgery, or otorhinolaryngology wards, for a wide spectrum of non-neoplastic conditions, including acute appendicitis, acute otitis media, inguinal hernia, and non-traumatic orthopaedic disorder, during the same period as the case admissions; those with a prior history of malignant neoplasms were excluded. The gender distribution of the two groups was similar with ~57% being men in both case and control groups. The demographic and lifestyle characteristics of the cases and controls are shown in Table 1. Venous blood was taken at the time of interview with written informed consent.

Table 1 Demographic and lifestyle characteristics of subjects

SNP Genotyping

After centrifugation of blood, DNA was extracted from the buffy coat using a commercial kit (QIAGEN) as described previously [16]. All participants were genotyped for 7 SNPs in selenoprotein genes: rs1050450 in GPX1, rs713041 in GPX4, rs5859 and rs5845 in the 15-kDa selenoprotein gene (SEP15), rs34713741 in SELS and rs7579 and rs3877899 in SEPP1. In addition, since an association between rs1050450 in GPX1 and breast cancer has been found to depend on genotype for the SNP in manganese superoxide dismutase gene that causes a Val/Ala change [9], the participants were also genotyped for SNP rs4880 in this gene. All SNPs with the exception of rs1050450 in GPX1 were genotyped using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry of primer extension products (Sequenom iPlex system, Hamburg, Germany), especially suited to high throughput genotyping. The primer sequences were determined with Sequenom SpectroDESIGNER software, and amplification primers and extension primers are shown in Table 2. All reactions were carried out in a final volume of 5 μl according to manufacturer’s conditions. Genotyping of 10% of the samples was repeated blind to confirm assay fidelity. The genotyping success rate was >93%.

Table 2 PCR primer sequences used for amplification of DNA for Sequenom analysis

Statistical analysis

The case–control association studies were analysed using χ² tests on 2 × 2 and 2 × 3 contingency tables for allele and genotype frequencies, respectively. No significant deviation from Hardy–Weinberg equilibrium was observed for any of the SNPs in this study, at the 5% level. Regression analysis was performed using the software program STATA (version 10.0). Initially, univariate logistic regression was used to assess the relationship between each variant and CRC risk. Subsequently, data were stratified according to disease subsite, and multivariate logistic analysis was carried out to calculate odds ratios, 95% confidence intervals (CI) and P for trend values after adjustment for age and gender only or for age, sex, intake of aspirin, intake of multivitamins, family history of colorectal cancer, smoking, alcohol, intake of red meat, vegetables and fruit, and physical activity.

Results

All the SNPs tested were in Hardy–Weinberg equilibrium. As indicated in Table 1, the control group was significantly younger than the CRC group, smoked less, had less family history of CRC and reported significant differences in alcohol, fruit and vegetable intake. Therefore, the genotype data were adjusted for age, intake of aspirin, intake of multivitamin, family history of colorectal cancer, smoking, alcohol, intakes of red meat, vegetables and fruits, and physical activity.

As shown in Table 3, after adjustment for age, sex and lifestyle factors, logistic regression analysis showed a mean odds ratio of 1.25 (95% CI 0.99,1.57, P for trend = 0.26) for carriage of one variant T allele of rs34713741 in SELS. Stratification of the data according to disease subsite (rectal or colon cancer) showed a statistically significant increase in risk of rectal cancer for carriage of one T allele of this SNP (OR 1.4 [95% CI 1.04,1.89], P for trend = 0.018) and a borderline increase (OR 1.57 [95% CI 0.98,2.51]) for the homozygous TT individuals (Table 4). This variant did not appear to alter risk of colon cancer. When the data were analysed according to gender, the association of the T variant allele for rs34713741 with increased risk of CRC was not observed in men but was associated with increased risk of rectal cancer in women (Table 5) with an odds ratios of 2.25 (95% CI 1.13,4.48, P for trend = 0.017) for TT individuals.

Table 3 Genotype frequencies in Korean patients with colorectal cancer and controls

Table 4 Genotype frequencies in Korean patients with colorectal cancer and controls stratified according to disease subsite

Table 5 Genotype frequencies in female Korean patients with colorectal cancer and controls stratified according to disease subsite

As shown in Table 6, analysis of the data according to gender also showed that in men there was an increased risk of rectal cancer was associated with carriage of one variant allele of either rs5845 or rs5859 in the SEP15 gene and this was of borderline statistical significance (OR 2.47; [95% CI 0.99,6.19], and OR 2.51; [95% CI 1.00,6.28], P = 0.052 and P = 0.049, respectively). There was, however, no increased risk in women or in patients with colon cancer. There was no difference in allele frequencies between patients with either rectal cancer or colon cancer and controls for rs1050450, rs713041, rs4880 and rs3877899.

Table 6 Genotype frequencies in male Korean patients with colorectal cancer and controls stratified according to disease subsite

Discussion

The results presented here show that for three SNPs in selenoprotein genes, rs5845 and rs5859 in the SEP15 gene and rs34713741 in the SELS gene, there were differences in allele frequency between patients with rectal cancer and healthy controls in a Korean population; the results suggest that the minor allele for all three SNPs is associated with increased disease risk. These SNPs are found in regulatory gene regions, one in the SELS promoter and two in the region of the SEP15 gene corresponding to the 3′UTR of the mRNA. Regulatory elements within the 3′UTR are essential for selenium incorporation into selenoproteins [2] and therefore, SNPs in gene regions corresponding to the 3′UTR of selenoprotein mRNAs have the potential to influence selenoprotein expression. Indeed, minor allelic variants of rs5845 and rs5859 in SEP15 have been shown to have functional consequences [13]. There are no available data on the functional effects of rs34713741, but it is situated very close (~100 bp) to rs28665122, a SNP that has been shown to modulate SelS expression and plasma levels of inflammatory cytokines [10]. This is the first report of an association of these SNPs with risk of CRC, but replication in a second and larger cohort is necessary to rule out the possibility that the findings could be explained by chance.

In the case of all three SNPs, the association of the minor allele with increased disease risk was modulated by gender. The effect of rs5845 and rs5859 in SEP15 on rectal cancer risk was only significant in men, and the effect of rs34713741 in the SELS gene was only found in women (Tables 5, 6). The precise reason for this is not clear but it most likely reflects gender-specific differences in selenoprotein metabolism since such differences have been observed previously [24, 26]. In addition, all three SNPs affected risk of rectal but not colon cancer. The reason for this is not clear, and nothing is known of selenoprotein expression in the rectum when compared to the colon.

Selenoprotein S is an endoplasmic reticulum protein involved in removing misfolded proteins, and SNP rs28665122 in the SELS promoter has been reported to be associated with altered levels of pro-inflammatory cytokines [10]. The observed difference in allele frequency of rs34713741 between patients with rectal cancer and controls (Tables 3, 4, 5) may reflect the close proximity (~100 bp) and likely linkage of rs28665122 and rs3471374 and changes in inflammatory events in the rectum, especially since altered inflammatory processes have been associated with increased CRC risk [17]. The previous failure to find any association between genotype for rs28665122 and susceptibility to inflammatory bowel disease [27] may reflect the small size of the population examined in that study. The present results suggest investigation of SELS association with gut inflammation using either mechanistic approaches or larger population groups would be worthwhile.

Sep15, a member of a novel family of selenoproteins, is also involved in protein folding mechanisms within the endoplasmic reticulum [11, 18]. Rs5845 and rs5859 are SNPs in a region of the SEP15 gene corresponding to the 3′UTR; they are closely linked, and their modulation of the ability of the 3′UTR to promote reporter gene activity [13] indicates that they have functional consequences. Furthermore, the SNPs have been reported to affect breast cancer risk [13] and risk of lung cancer in smokers [14]. The present results indicate that the SNPs influence rectal cancer risk in men.

Glutathione peroxidise 4 has functions in detoxification of reactive oxygen species, especially lipid hydroperoxides [4] and has been implicated in regulation of inflammatory and apoptotic pathways [28]. A previous study of rs713041, a T-C substitution in GPX4, indicated that in a cohort of Caucasians carriage of the T allele was associated with lower risk of with CRC [3]. The present data fail to confirm this observation, and this may reflect the larger size of the present cohort, ethnic or gender distribution differences in the populations studied (Caucasians in United Kingdom versus Koreans), or dietary and/or other environmental factors that differed in the two study groups.

In conclusion, the present study of a Korean population group shows that allele frequencies for three SNPs in selenoprotein genes, rs5845 and 5859 in the SEP15 gene and rs34713741 in the SELS gene differ between patients with rectal cancer and healthy controls. The results suggest that the minor alleles for all three SNPs are associated with increased rectal cancer risk and that the effects of the SNPs are dependent on gender. It is notable that these SNPs are in genes that code for two selenoproteins involved in the metabolism of unfolded or misfolded proteins. Thus, the present results suggest that Se and these two selenoproteins have an important role in the unfolded protein response in cells of the large bowel and that genetic variation in these mechanisms may underlie rectal cancer susceptibility. Further studies are required to define the functional effects of the variants, especially rs34713741, to confirm whether the same effects on disease risk are seen in other ethnic populations and to investigate whether the effects of the variants are influenced by dietary Se intake.