American Journal of Clinical Nutrition, Vol. 87, No. 4, 1019-1027,
April 2008
© 2008 American Society for Nutrition
ORIGINAL RESEARCH COMMUNICATION |
Functional effects of a common single-nucleotide polymorphism (GPX4c718t) in the glutathione peroxidase 4 gene: interaction with sex1,2,3
Catherine Méplan,
Lynne K Crosley,
Fergus Nicol,
Graham W Horgan,
John C Mathers,
John R Arthur and
John E Hesketh
1 From the Institute for Cell and Molecular Biosciences (CM and JEH), the Human Nutrition Research Centre (CM, JCM, and JEH), and the School of Clinical Medical Sciences (JCM), The Medical School, Newcastle University, Newcastle, United Kingdom; BioSS, Aberdeen, United Kingdom (GWH); and the Rowett Research Institute, Aberdeen, United Kingdom (LKC, FN, and JRA)
2 Supported by grant no. N05041 from the Food Standard Agency and by the Scottish Executive Environment and Rural Affairs Department (to JRA).
3 Reprints not available. Address correspondence to JE Hesketh, Institute for Cell and Molecular Biosciences, The Medical School, Newcastle University, Framlington Place, Newcastle-upon-Tyne, NE1 4HH, United Kingdom. E-mail: j.e.hesketh{at}ncl.ac.uk.
 |
ABSTRACT
|
|---|
Background: Selenium is essential for health in humans. Selenium is present as selenocysteine in selenoproteins such as the glutathione peroxidases (GPx). Selenocysteine incorporation requires specific structures in the 3'untranslated region (3'UTR) of selenoprotein mRNAs.
Objective: This study investigated the functional significance of the single-nucleotide polymorphism (SNP) GPx4c718t within the 3'UTR of the GPx4 gene.
Design: A selenium supplementation trial was carried out with prospectively genotyped individuals of both homozygote genotypes for this SNP. Blood samples were analyzed at baseline, after a 6-wk supplementation with 100 µg Se as sodium selenite/d, and during a 6-wk washout period. RNA-protein binding studies were carried out in vitro.
Results: Both lymphocyte GPx1 protein concentrations and plasma GPx3 activity increased significantly after selenium supplementation in CC but not TT participants. After selenium withdrawal, there was a significant fall in both lymphocyte GPx4 protein concentrations and GPx4 activity in TT but not in CC participants; this effect was modulated by sex. RNA-protein binding assays showed that both T and C variants of transcripts corresponding to the GPx4 3'UTR formed complexes in vitro and that the C variant bound more strongly than did either the T variant or the GPx1 3'UTR.
Conclusions: The GPX4c718t SNP both alters protein binding to the 3'UTR in vitro and influences the concentration of lymphocyte GPx4 and other selenoproteins in vivo. The latter is consistent with competition for selenium in selenoprotein synthesis, and, at low selenium intake, the SNP thus may influence susceptibility to disease.
 |
INTRODUCTION
|
|---|
The micronutrient selenium is essential for human health. Relatively high intakes can be associated with lower risks of cancer, cardiovascular disease (1-3), and viral infection (4). Conversely, lower selenium status is associated with a greater risk of some diseases (1-3) and greater mortality (5). The biological functions of selenium are exerted predominantly by means of the amino acid selenocysteine that is present in selenoproteins involved in many aspects of metabolism (6).
The incorporation of selenocysteine into selenoproteins involves the recoding of a UGA codon and the formation of a complex involving an RNA stem-loop structure, the selenocysteine-insertion sequence (SECIS), which is located in the 3'untranslated region (3'UTR) of selenoprotein mRNA (7) and a number of RNA-binding proteins including SBP2 protein (SECIS-binding protein 2) and ribosomal protein L30, a specific elongation factor (EFsec) and the tRNA(ser)sec (8-14).
The phospholipid hydroperoxide glutathione peroxidase (GPx4) is an intracellular antioxidant selenoprotein that was identified first by its ability to reduce lipid hydroperoxides, as reviewed by Brigelius-Flohé (15). GPx4 is expressed in most tissues but is found at particularly high concentrations in the testes, where it also functions in sperm maturation. The protein has also been proposed to have roles in leukotriene biosynthesis and the regulation of cytokine signaling pathways (15). GPX4-null mice die in utero between stages E7.5 and E8.5, which confirms its physiologic importance (16, 17) The lethal phenotype can be rescued by overexpression of a human GPX4 transgene that protects the animals from oxidative damage (18, 19).
A limited number of studies have investigated genetic variations within the GPX4 gene that could potentially alter its function (20-24). Although the level of GPx4 expression is correlated with sperm function, no link has been found between a series of genetic variations and fertility (20-22). In whites, there is a C
T single-nucleotide polymorphism (SNP) located within the GPx4 gene region corresponding to the 3'UTR, near the SECIS element at position 718 (23). The C variant of this SNP has a higher frequency in patients with colorectal cancer (24), which is suggestive of a link between the SNP and disease. This GPX4c718t polymorphism (rs713041) could modulate the synthesis of GPx4 by altering the affinity of the selenocysteine insertion machinery for its SECIS element. However, to date there is no direct evidence of its functional effects on GPx4 synthesis or selenoprotein metabolism. The aim of the present work was to test the hypothesis that the GPX4c718t SNP has functional consequences. To do this, we assessed lymphocyte GPx4 and other selenoprotein biomarkers in prospectively genotyped healthy volunteers before and after selenium supplementation and assessed the effects of the single-nucleotide change in vitro in RNA-protein binding assays.
 |
SUBJECTS AND METHODS
|
|---|
Selenium supplementation trial
Forty unrelated volunteers (male and female) aged 20–60 y who were nonsmokers and generally healthy were recruited from the general population in and around Newcastle (United Kingdom). Exclusion criteria included cardiovascular, hepatic, gastrointestinal and thyroid disorders; cancer; excessive alcohol consumption (>30 units/wk); and chronic intake of antiinflammatory drugs (25). Persons already taking selenium or multivitamin or vitamin E supplements also were excluded. For all participants, weight and height were measured, and basal metabolic index (BMI) was calculated. Peripheral blood samples were drawn from the antecubital vein into 10-mL EDTA-tubes (BD Vacutainer; Becton-Dickinson), processed on the same day, typically within 6 h of sample collection, and stored at –80 °C until assays were performed.
An initial 10-mL blood sample was collected for GPx4 genotype assessment. Volunteers were asked to give an additional 30-mL blood sample (baseline) and then take a daily supplement of 100 µg Se as sodium selenite (Cardinal Health) for 6 wk. At the end of the supplementation period, another 30-mL blood sample was taken, and then 3 more blood samples were taken at 2-wk intervals during a 6-wk washout period (Figure 1
A). Overall compliance with the selenium supplementation was estimated by counting the number of supplement capsules returned after the trial.
Written informed consent was obtained from all participants. The study protocol was approved by the Sunderland Local Research Ethics Committee.
DNA extraction and genotyping
DNA extraction from fresh blood was performed on buffy coat by using the Qiagen QIAmp DNA blood mini-kit according to the manufacturer's specifications. Polymerase chain reaction (PCR) was performed with 100 ng template genomic DNA, using forward (GACCTG CCCCACTATTTCTA) and reverse (GTCTGTTTATTCCCACAAGG) primers and an Expand High-Fidelity PCR system (Roche) in a ThermoHybaid Px2 thermocycler under the following conditions: an initial denaturating step at 94 °C for 4 min, which was followed by 30 cycles of denaturation (94 °C for 30 s), annealing (53.5 °C for 30 s), and extension (72 °C for 1 min). The PCR was completed by a final extension cycle at 72 °C for 7 min. An aliquot of each PCR product was subjected to electrophoresis in a 1x Tris-acetate-EDTA-2%-agarose gel and visualized with ethidium bromide. Genotype was determined by restriction fragment analysis, digesting the PCR products with StyI (Promega) overnight at 37 °C to generate either 2 fragments of 159 and 62 bp (C allele) or 3 fragments of 97, 62, and 62 bp (T allele). Subsequently digested products were resolved on a 4%-Tris-borate-EDTA polyacrylamide gel and visualized with ethidium bromide. Genotype was confirmed by direct sequencing of PCR product for 10% of samples. For this, PCR products were purified by using a polyethylene glycol (PEG) mix (26.2% wt:vol PEG 8000, 6.6 mmol MgCl2/L, 0.6 mol sodium acetate/L; pH 5.2), washed with 70% (by vol) ethanol, air-dried, and sent for sequencing (MWG-Biotech).
Plasma and erythrocyte preparation, lymphocyte isolation, and measurement of selenoprotein markers
After fractionation of the blood by centrifugation (950 x g, 4 °C, 15 min), plasma was further centrifuged at 730 x g for 12 min at 4 °C to obtain platelet-free plasma. Erythrocytes were frozen immediately. Lymphocytes were reconstituted in an equal volume of RPMI 1640/Hepes medium (GibcoBRL, Invitrogen) and isolated onto a Histopaque H1077 gradient (Sigma-Aldrich) as described elsewhere (25). After isolation, lymphocyte pellets were stored at –80 °C. Lymphocyte GPx1 activity and plasma GPx3 activity were calculated by the method of Paglia and Valentine (26), as modified by Brown et al (27), using hydrogen peroxide as a substrate. One unit of GPx1 or GPx3 activity is defined as that which oxidizes 1 µmol NADPH/min. Lymphocyte GPx4 activity was measured by the method of Weitzel et al (28) using phosphatidyl choline hydroperoxide as a substrate. One unit of GPx4 activity is defined as that which oxidizes 1 µmol NADPH/min. Lymphocyte GPx1 and GPx4 protein concentrations were measured by competitive enzyme-linked immunosorbent assay [ELISA (29)] using polyclonal antibodies and recombinant proteins GPx1 and GPx4 obtained from Labfrontier.
Cell culture and protein extraction from colon adenocarcinoma cells
Human colon adenocarcinoma cells (Caco-2) were grown in Dulbecco's Modified Essential Medium (Invitrogen) supplemented with 10% fetal calf serum (Sigma), 1% penicillin and streptomycin (Invitrogen), and 1% nonessential amino acid (Invitrogen). Cells were grown to
90% confluence, and cytoplasmic protein extracts were prepared by a modification of the method of Mahtani et al (30). Cells were cooled on ice for 5 min, rinsed once in ice-cold phosphate-buffered saline (PBS), and harvested by scraping. Cells were pelleted by centrifugation at 600 x g for 10 min at 4 °C and washed once in ice-cold PBS. The pellet was then lysed by the addition of lysis buffer [10 mmol HEPES/L, pH 7.6; 3 mmol MgCl2/L, 40 mmol KCl/L, 2 mmol dithiothreitol/L, 5% glycerol, 0.5% NP-40, and mini-EDTA–free protease inhibitor cocktail (Roche)] using 100 µL buffer/20 x 106 cells. Cells were lysed on ice for 10 min and disrupted by being passed 10 times through a 21-gauge needle. Nuclei were removed by centrifugation at 600 x g for 10 min at 4 °C. The supernatant fluid containing cytoplasmic proteins was aliquoted and snap-frozen. Protein content was measured by the Bradford method (Sigma).
In vitro transcription
Templates for transcription of human GPX4 T and C variants and GPX1 sequences corresponding to the regions of the 3'UTRs containing the SECIS structure were generated by PCR from genomic DNA from participants who were either TT or CC for GPX4c718t. PCR was performed with 100 ng template genomic DNA by using forward primers that contained the T7 promoter sequence (underlined) followed by a sequence specific for GPX4 (TAATACGACTCACTATAGGGACCTGCCCCACTATTTCTA) or GPX1 (TAATACGACTCACTATAGGCGAAGCCCTGCTGTCTCAAG), and reverse primers specific for GPX4 (GTCTGTTTATTCCCACAAGG) or GPX1 (CTGACACCCGGCACTTTATTAGTG). Annealing temperatures were 53.5 °C for GPX4 and 64 °C for GPX1. PCR products were purified by using a QIAquick PCR purification kit (Qiagen) and were of 239 (T7-GPX4) and 253 (T7-GPX1) base pairs, respectively. In vitro transcription was performed, according to the manufacturer's protocol, by using either a digoxigenin (DIG)-RNA-labeling kit (Roche) for the synthesis of labeled GPX4 probes corresponding to C or T variants or a MEGAshortscript kit (Ambion) to generate competitor unlabeled probes for GPX1 and the 2 GPX4 variants. Probes were subsequently purified by a phenol and chloroform extraction and ethanol precipitation. The precipitated RNA probes were finally resuspended in diethylpyrocarbonate-treated water in the presence of RNasin (Promega).
Electromobility-shift assay
Electromobility-shift assay (EMSA) reactions were carried out on ice for 25 min. Protein extract (20 µg) was incubated with 100 pmol digoxigenin-labeled RNA probe (previously heated to 70 °C and cooled to 40 °C at a rate of 1.3 °C/min to refold slowly in the correct conformation) in binding buffer (20 mmol HEPES/L, pH 7.6; 3 mmol MgCl2/L; 40 mmol KCl/L; 2 mmol dithiothreitol/L; protease inhibitor cocktail; and 5% glycerol) in a final 20 µL-reaction volume. For competition experiments, unlabeled transcripts were added simultaneously with the labeled transcripts. After the binding reaction, 40 units of RNase T1 were added and incubation was continued for 20 min on ice. Loading buffer (3 µL; 90% glycerol and 0.025% bromophenol blue) was then added, and complexes were separated by electrophoresis at 4 °C for 2.5 h at 20 V/cm by using 4% (wt:vol) nondenaturing PEGs in 0.5x Tris-borate EDTA. RNA-protein complexes were subsequently transferred onto a Hybond-N+ nylon membrane (Amersham Biosciences) for 1 h at 10V and UV-crosslinked to the membrane for 2 min. The membrane was then incubated in a blocking solution (2% wt:vol bovine serum albumin and 0.5 ng salmon sperm DNA in PBS) for 1 h, which was next incubated for 1 h with anti-DIG-alkaline phosphatase antibody (1:5000 dilution, Roche). After 3 washes in PBS, the membrane was rinsed in milliQ water and incubated with BM purple reagent (Roche) for 10 min to show the digoxigenin-labeled RNA-protein complexes. Band intensity was quantified by using UVIPHOTOMW software (version 11.01; Uvitec).
Statistical analysis
Biochemical variables were analyzed by using factorial analysis of variance (ANOVA). In the case of time differences, the responses to selenium supplementation and withdrawal were analyzed by using a paired t test of the values at different pairs of timepoints, and the interactions between time and other factors were examined by ANOVA of these differences. Results of the paired t test are reported only when the interaction between time and other factors was significant at P < 0.05. When a significant interaction between 2 factors was observed by ANOVA, secondary analysis was carried out to determine the main effects of each factor within subgroups. One-factor ANOVA was used to compare the binding capacity of competitor probes in the EMSA experiment.
 |
RESULTS
|
|---|
Selenium supplementation study
To investigate the functionality of GPX4c718t, a selenium supplementation trial (the SELGEN study) was carried out in 40 prospectively genotyped healthy human volunteers, of whom 22 were CC and 18 were TT for the GPX4c718t SNP. As part of the SELGEN study, volunteers also were genotyped for 2 SNPs in the selenoprotein P (SEPP) gene. Because statistical analysis of the genotypes frequencies showed that GPX4c718t was independent of the 2 SEPP SNPS, for clarity, the data are presented here only in terms of GPX4c718t. The study design is shown in Figure 1
. Volunteer characteristics and compliance are summarized in Table 1
. Volunteers were asked to take a single daily selenium supplement as sodium selenite (100 µg Se/d) for 6 wk. Blood samples were drawn before supplementation and on the last day of supplementation and at 2-wk intervals during a 6-wk washout period. Baseline (presupplementation) mean plasma selenium concentrations were 1.15 ± 0.03 µmol/L, and the concentration was significantly (P < 0.001) increased, by 20%, to 1.36 ± 0.02 µmol/L by selenium supplementation (Table 1
). In addition, supplementation resulted in a significant increase in SEPP concentrations (25) and plasma GPx3 activity (see below). The increases in plasma selenium and SEPP concentrations and in GPx3 activity indicate both the subjects' compliance and the effectiveness of the supplementation in raising selenium status. Selenium withdrawal led to a rapid decrease in plasma selenium by 2 wk and a further fall at 4 wk, after which time the concentration was similar to baseline (25). Most of the decline in plasma selenium occurred by 2 wk.
Effects of sex and GPX4c718t on GPx4 activity and protein concentrations
Lymphocyte GPx4 activity and protein concentrations were measured in samples taken before and after supplementation and during the washout period. Sex had a significant (P = 0.028) effect on GPx4 protein concentration at baseline: females had higher concentrations than did males (Figure 1B
), but no effect of sex was observed after supplementation or during the washout period. In contrast, males had a higher mean GPx4 activity than did females at baseline, but this difference was not significant (Figure 1C
). As a result, a sex effect was observed on the ratio of GPx4 activity to GPx4 protein at baseline (P = 0.027): males had a higher ratio than did females (Figure 1D
). This difference may reflect differences between the sexes in the regulation of enzyme activity.
Presupplementation GPx4 protein concentrations and activity levels were similar in participants who were TT or CC for GPX4c718t (Figure 2
A and B). Supplementation had little effect on either activity or protein, but, after the withdrawal of selenium, there was a fall in GPx4 protein concentrations and activity in TT but not in CC persons. Significant genotype x time, sex x time, and genotype x sex x time interactions were observed for GPx4 protein concentrations, particularly between before supplementation and after supplementation (P = 0.047) and 2 wk after withdrawal of supplementation (P = 0.017). ANOVA confirmed that the difference in GPx4 activity between genotypes was significant at 2 (P = 0.022) and 4 (P = 0.034) wk after selenium withdrawal, with lower activity in TT participants. A similar effect of the GPX4c718t genotype was observed on GPx4 protein after a 2-wk withdrawal (P = 0.022). Moreover a significant sex x GPX4c718t interaction was observed for GPx4 protein 4 wk after withdrawal (P = 0.044). Most of this genotype effect was attributable to females—lower GPx4 protein concentrations were seen in TT than in CC females after supplementation and at 2 (P = 0.018) and 4 (P = 0.031) wk of the washout—and the difference was less evident in males. Overall, the data indicate that, during washout of the selenium supplement, GPx4 activity and protein concentration are influenced by genotype for GPX4c718t and, in addition, GPx4 protein concentrations are influenced by sex.

View larger version (16K):
[in this window]
[in a new window]
|
FIGURE 2.. Mean (±SEM) effects of GPX4c718t genotype on lymphocyte GPx4 activity and protein concentration. GPx4 activity was measured by spectrophotometric assay, and protein concentrations were measured by using an enzyme-linked immunosorbent assay. Presupp, before supplementation (baseline); postsupp, after supplementation. GPx4 activity (B) and protein concentrations (A, C, and D) are shown for groups of persons of CC (n = 22) and TT (n = 18) genotype for GPX4c718t. Data are shown for the different genotypes in all volunteers (B) and when separated into females (C) and males (D). *Both GPx4 protein and activity differed significantly between CC and TT participants at 2 or 4 wk after withdrawal of supplementation, P < 0.05 (ANOVA).
|
|
Effect of polymorphism on other selenoprotein biomarkers
There is a hierarchy in selenoprotein synthesis that results from a competition for available selenium and for components of the selenoprotein synthetic machinery (31-33). Therefore, in addition to GPx4 protein and activity, plasma GPx3 activity, lymphocyte GPx1 activity, and protein and plasma TR1 concentrations were measured. GPX4c718t genotype affected the response to supplementation of the GPx1 protein concentration and its decrease 2 and 6 wk after selenium withdrawal (Figure 3
A): persons of CC genotype had significant differences in GPx1 protein (P = 0.019, 0.019, and 0.023, respectively; paired t test) (Figure 3A
). Further statistical analysis showed that the increase in GPx1 protein concentrations after supplementation was significant (P = 0.017, paired t test) in CC females (Figure 3B
) but not in males (Figure 3C
). Both the difference between sex and the sex x time interaction were significant (post- to presupplementation differences—genotype x sex interaction, P = 0.033; sex main effect, P = 0.006; week 6 washout-postsupplementation differences—genotype x sex interaction, P = 0.053; sex main effect, P = 0.007). In addition, sex x time x genotype interactions were significant. Finally, TT females showed a rapid drop in GPx1 activity by 2 wk after selenium withdrawal (P = 0.016), but TT males showed no significant decrease in GPx1 activity during the washout period (data not shown).

View larger version (7K):
[in this window]
[in a new window]
|
FIGURE 3.. Mean (±SEM) effects of GPX4c718t genotype on lymphocyte GPx1 protein concentrations, measured by using an enzyme-linked immunosorbent assay. Presupp, before supplementation (baseline); postsupp, after supplementation. Data are shown from 22 CC (11 female, 11 male) and 18 TT (13 female, 5 male) participants on the basis of GPX4c7187t genotype (A) and of genotype for females (B) and males (C). *Significant responses in GPx1 activity to selenium supplementation and withdrawal were assessed by paired t test. *, **Significant differences were indicated by 2 symbols joined by a line: *P < 0.05, **P < 0.01.
|
|
In addition, the GPX4c718t genotype influenced plasma concentrations of the selenoproteins GPx3 and TR1. GPX4c718t genotype affected GPx3 activity after supplementation (P = 0.029): CC participants had a higher GPx3 activity than did TT participants (Figure 4
A). Moreover a significant response to selenium supplementation occurred only in CC volunteers (P = 0.005, paired t test), and this difference in response was significant between genotypes (P = 0.045, ANOVA). Similarly, there was a significant sex x GPX4c718t genotype interaction that determines plasma TR1 concentrations at all time-points (P
0.045 for all; Figure 4B
). TT males and CC females showed higher TR1 concentrations than did CC males and TT females after supplementation and throughout the washout period (Figure 4B
).

View larger version (9K):
[in this window]
[in a new window]
|
FIGURE 4.. Mean (±SEM) effects of GPX4 GPX4c718t genotype on plasma GPx3 activity and plasma thioredoxin reductase 1 (TR1) protein concentrations. Presupp, before supplementation (baseline); postsupp, after supplementation. GPx3 activity (A) was measured by spectrophotometric assay, and TR1 concentrations (B) were measured by radioimmunoassay. A: Data are presented for persons of CC (n = 22) and TT (n = 18) genotype. *Significant difference between CC and TT participants after supplementation, P < 0.05 (ANOVA). *, **Significant responses to supplementation and wash-out at 2 wk (paired t test): *P < 0.05, **P < 0.01. B: Data are shown for groups of males (5 TT, 11 CC) and females (13 TT, 11 CC) as a function of genotype for GPX4c718t. *, **, ***Significant differences between CC and TT participants at all time-points (ANOVA): *P < 0.05, **P < 0.01, ***P < 0.001).
|
|
Effect of GPx4 polymorphism on in vitro RNA-protein binding
EMSA experiments were carried out to determine whether GPX4c718t SNP modulated the ability of RNA transcripts containing the GPX4 3'UTR (including the SECIS) to form specific RNA-protein complexes and compete with GPx1 mRNA for protein binding. Both T and C variants of digoxigenin-labeled GPx4 transcripts formed complexes with proteins extracted from Caco-2 cells, as judged by lower transcript mobility (Figure 5
A, B). Relative binding capacity of the 2 variants was assessed by the incubation of labeled C or T variant with an increasing molar excess of unlabeled competitor transcripts corresponding to either itself, the alternative variant, or the GPX1 3'UTR.

View larger version (23K):
[in this window]
[in a new window]
|
FIGURE 5.. Effects of GPX4c718t genotype on RNA-protein complex formation by electromobility-shift assay (EMSA). Transcripts of either C or T variant of the GPX4 3'untranslated region were labeled with digoxigenin (DIG)-UTP and incubated with protein extracts from colon adenocarcinoma cells. Comp, competitor; DIG, digoxigenin; FP, free probe. EMSAs were carried out using either DIG-T–(A) or DIG-C–(B) labeled transcripts in the absence or presence of competitors. A: Complex formation with DIG-T transcripts and effective competition with unlabeled T transcripts at 10–15x molar excess and with unlabeled C competitor transcripts at 5x molar excess. B: Complex formation with DIG-C transcripts and effective competition with unlabeled C competitor transcripts at 5x molar excess and with unlabeled GPx1 transcripts at 10–15x molar excess. The extent of complex formation was quantified by densitometry in 4–5 independent experiments, and the mean (±SEM) amount of complex formation relative to the absence of competitor (100%) in the presence of increasing molar excess of unlabeled RNA competitors is shown in C and D. The dotted line rectangle in A indicates the RNA-protein complexes used for the quantification. The data from experiments with labeled T transcripts (C) showed significant differences in the extent of complex formation between T and C competitors at 10x and 15x molar excess (*) and a significant difference between C and GPx1 competitors at 15x molar excess (#) (ANOVA). In experiments with labeled C transcripts, significant differences were found between C and both T and GPx1 competitors at 5x, 10x, and 15x molar excess (* for T competitors; # for GPx1 competitors) (ANOVA). *, **, ***Significant differences: *P < 0.05, **P < 0.01, ***P < 0.001. #, ##, ###Significant differences: #P < 0.05, ##P < 0.01, ###P < 0.001.
|
|
Both the unlabeled T and C transcripts could compete with their own digoxigenin-labeled equivalent so that complex formation was no longer evident with a 5–10x molar excess (Figure 5A, B
). Quantification of band intensities showed that unlabeled C transcripts competed strongly with labeled T transcripts at 5x molar excess so that complex formation was disrupted by 45%, whereas competition of unlabeled T transcripts was significantly weaker, and complex formation was affected by only 35% at 10–15x molar excess (P = 0.0013, one-factor ANOVA; Figure 5
C). In the complementary experiments, the C variant competed more strongly with itself than with an unlabeled T competitor (Figure 5D
). These competition experiments suggest that the C variant transcripts bind protein more strongly than do the T variant transcripts.
Unlabeled GPX1 transcripts competed with labeled T GPx4 transcripts for complex formation (Figure 5B
), and quantification of band intensities showed that this competition required a higher molar excess (10–15x) of GPX1 competitor compared with the C competitor, which indicated that the C variant has stronger protein binding than does either the T variant or the GPX1 SECIS element (Figure 5C,D
). GPX1 competitor transcripts lowered RNA-protein complex formation with the T variant by 60% and that with the C variant by only 30% (Figure 5C, D
). Competition experiments with unlabeled rat metallothionein-I 3'UTR transcripts (a negative control lacking a SECIS element) showed no competition for binding to GPX4 and therefore confirmed the specificity of protein binding to GPX4 and GPX1 probes (data not shown). Overall, the data indicate that the C variant GPX4 transcripts bind proteins more strongly than do the T variant and GPX1 transcripts, which is consistent with the in vivo observation that GPX4c718t influences both GPX4 and GPX1 protein expression (Figure 3
).
Effect of GPx4 polymorphism on relative activity of GPx1 and GPx4
To assess the relative effects of the SNP on GPx1 and GPx4 in vivo, the ratio of GPx1 activity to GPx4 activity was calculated from the lymphocyte enzyme activities determined in the supplementation trial (Figure 6
). An effect of GPX4c718t was evident 2 wk after selenium withdrawal (P = 0.047). For CC participants, the ratio of the 2 enzyme activities returned to presupplementation values after 2 wk of selenium withdrawal, whereas the ratio stayed significantly higher in TT participants during the washout period at 2 and 4 wk. This observation suggests that CC persons maintain GPx4 synthesis preferentially over GPx1 when selenium intake falls, which is consistent with the stronger competition of the C variant transcripts in the RNA-binding assays.

View larger version (14K):
[in this window]
[in a new window]
|
FIGURE 6.. Mean (±SEM) effects of GPX4c718t genotype on the ratio of lymphocyte GPx1 activity to lymphocyte GPx4 activity, measured by spectrophotometric assay. Values shown are for groups of participants of CC (n = 22) and TT (n = 18) genotype. Presupp, before supplementation (baseline); postsupp, after supplementation. *Significant differences between CC and TT participants at 2 and 4 wk after withdrawal of the selenium supplement, P < 0.05 (ANOVA).
|
|
 |
DISCUSSION
|
|---|
SNPs in selenoprotein genes have been described, but detailed data on their functionality are less comprehensive (25, 34-37). A T/C variation has been described in the GPX4 gene at position 718 with the 3'UTR (23) close to the predicted SECIS structure, which is essential for selenoprotein synthesis (38), and this SNP has recently been confirmed (36). All 3 genotypic variants at this locus have been reported to occur at relatively high frequency in whites, Chinese, and South Asian ethnic groups (results not shown), which indicates that, if functional, this SNP could affect a significant number of people. However, to date, there is only indirect evidence—namely, measurement of 5'lipoxygenase metabolite concentrations (23)—that the SNP may be functional. The present work provides direct evidence, both in vivo and in vitro, that this SNP has functional consequences related to selenium metabolism.
First, the fall in GPx4 activity and protein concentrations in lymphocytes after withdrawal of a daily selenium supplement is affected by GPX4c718t genotype (Figure 2
). The data suggest that CC persons maintain GPx4 concentrations better than do TT persons when selenium intake falls. The effects of the SNP were observed after withdrawal of the selenium supplement but not at baseline. This may be because initial interindividual variations in selenium intake mask the influence of the polymorphism at baseline. Alternatively, revelation of the genotype effect after selenium supplementation is stopped may be due to the sudden decrease in selenium intake that leads to a transient negative selenium balance and the need to prioritize the available selenium.
Second, in vivo, GPX4c718t influences the concentrations of other selenoprotein biomarkers. After the supplementation, plasma GPx3 increased in CC but not in TT participants, and plasma TR1 was greater in CC than in TT females. In addition, presupplementation concentrations of GPx1 were higher in TT than in CC participants. At first, it may be surprising that an SNP in one selenoprotein gene influences the expression of other selenoproteins. However, there is a strong body of evidence of competition among selenoproteins for available selenium, in the form of selenocysteine that is required for synthesis of the selenoproteins, and for the selenium incorporation machinery (31-33, 39, 40). Furthermore, the mechanism for this competition may involve, at least in part, differences in 3'UTR sequences (32, 33, 41-43). GPX4 variants may influence the position of GPx4 in the hierarchy of demands on the components of the selenoprotein synthesis apparatus, with consequential effects on the synthesis of other selenoproteins.
Third, in vitro binding experiments showed that GPX4 transcripts with the C variant compete more strongly with either the T variant transcripts or GPX1 3'UTR transcripts for protein binding, which suggests that the C and T variants differ in their ability to form RNA-protein complexes. Our hypothesis is that the T and C allelic variants differ in their capacity to promote selenocysteine incorporation into GPx4; this possibility is supported by recent data indicating that the C variant of the GPX4 3'UTR is stronger than the T variant at driving synthesis of a selenoprotein reporter (24). In addition, the hierarchy in selenoprotein synthesis (31-33, 44) would predict that GPX4c718t CC participants could maintain their GPx4 concentrations as selenium intake falls, whereas the abundances of some other selenoproteins would be more sensitive to selenium supply; this prediction was confirmed in the present experiments for GPx1, TR1, and GPx3. Because the selenoproteins and their mRNAs are present in different abundances [eg, there is 4 times as much GPx1 as GPx4 in many tissues (44)]; the effects of the T/C variation in the 3'UTR most likely reflect not only a competition for available selenocysteine but also the ability of different transcripts to compete for the selenoprotein synthetic machinery. Thus, the stronger in vitro protein binding of C variant than of the GPX1 3'UTR is paralleled by a greater response of GPx1 protein to supplementation in persons of CC genotype than in those of TT genotype.
The effect of GPX4c718t on GPx4 expression was dependent on sex, and it was evident in females but not males. This could result from differing, perhaps less stringent, demands for GPx4 in females that would make synthesis of the protein susceptible to the SNP. Additional regulation of GPx4 expression in males may apply because of the role of GPx4 in male fertility or of sex differences in immune function. It it interesting that, although the effects of the SNP on GPx4 activity and protein concentrations after withdrawal of the selenium supplement occurred in parallel, sex effects on activity and protein concentrations at baseline did not. As a result, GPx4 activity:GPx4 protein was lower in females than in males, which indicates that regulatory processes may be affecting GPx4 expression in males, who maintain lymphocyte GPx4 at a higher concentration than do females. The mechanism underlying this effect is not known, but it may well involve differences in the stability of GPx4 protein. To our knowledge, this report is the first time an effect of sex on human lymphocyte GPx4 expression has been described. However, in rats, there is a similar sexual dimorphism in the synthesis of the selenoprotein Dio1 in the liver and kidney (45). This sex difference was abolished by ovariectomy or orchidectomy, which suggests a role of sex steroids in posttranscriptional regulation of selenoproteins.
In conclusion, the present data provide evidence that an SNP in the GPX4 gene, a T/C variation at position 718 close to the SECIS, has functional consequences. Pilot studies suggest that this SNP may affect susceptibility to colon cancer (24). Such observations, taken with the present data on functionality, indicate the potential importance of this SNP in the GPX4 gene and could be especially relevant when selenium intake is suboptimal. We propose that the SNP alters protein binding to the 3'UTR, which in turn modulates the ability of the GPX4 3'UTR both to promote GPx4 synthesis and to compete for components of the selenoprotein synthetic machinery—and thus influence the synthesis of a range of selenoproteins
 |
ACKNOWLEDGMENTS
|
|---|
We thank GJ Beckett and AF Howie (University of Edinburgh) for TR1 measurements and advice, B Burtle and W Bal (Newcastle University) for technical support, F McArdle (University of Liverpool) for ICPMS measurements, and the volunteers for their participation in the SELGEN Study.
The authors' responsibilities were as follows—CM: the supplementation trial, genotyping, RNA protein binding; and data analysis; LKC and FN: measurements of selenoprotein biomarkers; GWH: helped with the statistics; JEH, JCM, and JRA: designed and supervised the project; and JEH and CM: wrote the manuscript with contributions from JCM and JRA. None of the authors had a personal or financial conflict of interest.
 |
REFERENCES
|
|---|
- Combs GF Jr. Current evidence and research needs to support a health claim for selenium and cancer prevention. J Nutr 2005;135:343–7.[Abstract/Free Full Text]
- Rayman MP. Selenium in cancer prevention: a review of the evidence and mechanism of action. Proc Nutr Soc 2005;64:527–42.[Medline]
- Whanger PD. Selenium and its relationship to cancer: an update dagger. Br J Nutr 2004;91:11–28.[Medline]
- Beck MA, Nelson HK, Shi Q, et al. Selenium deficiency increases the pathology of an influenza virus infection. FASEB J 2001;15:1481–3.[Abstract/Free Full Text]
- Akbaraly NT, Arnaud J, Hininger-Favier I, Gourlet V, Roussel AM, Berr C. Selenium and mortality in the elderly: results from the EVA study. Clin Chem 2005;51:2117–23.[Abstract/Free Full Text]
- Méplan C, Pagmantidis V, Hesketh JE. Advances in selenoprotein expression—patterns and individual variations. In: Brigelius-Flohe R, Joost H-G, eds. Nutritional genomics, nutrients, genes and genetic variation in health and disease. Weinheim, Germany: Wiley-VCH, 2006:132–58.
- Berry MJ. Knowing when not to stop. Nat Struct Mol Biol 2005;12:389–90.[Medline]
- Copeland PR, Fletcher JE, Carlson BA, Hatfield DL, Driscoll DM. A novel RNA binding protein, SBP2, is required for the translation of mammalian selenoprotein mRNAs EMBO J 2000;19:306–14.[Medline]
- Fagegaltier D, Hubert N, Yamada K, Mizutani T, Carbon P, Krol A. Characterization of mSelB, a novel mammalian elongation factor for selenoprotein translation. EMBO J 2000;19:4796–805.[Medline]
- Fletcher JE, Copeland PR, Driscoll DM, Krol A. The selenocysteine incorporation machinery: interactions between the SECIS RNA and the SECIS-binding protein SBP2. RNA 2001;7:1442–53.[Abstract]
- Tujebajeva RM, Copeland PR, Xu XM, et al. Decoding apparatus for eukaryotic selenocysteine insertion. EMBO Rep 2000;1:158–63.[Medline]
- Zavacki AM, Mansell JB, Chung M, Klimovitsky B, Harney JW, Berry MJ. Coupled tRNA(Sec)-dependent assembly of the selenocysteine decoding apparatus. Mol Cell 2003;11:773–81.[Medline]
- Chavatte L, Brown BA, Driscoll DM. Ribosomal protein L30 is a component of the UGA-selenocysteine recoding machinery in eukaryotes. Nat Struct Mol Biol 2005;12:408–16.[Medline]
- Small-Howard A, Morozova N, Stoytcheva Z, et al. Supramolecular complexes mediate selenocysteine incorporation in vivo. Mol Cell Biol 2006;26:2337–46.[Abstract/Free Full Text]
- Brigelius-Flohe R. Tissue-specific functions of individual glutathione peroxidases. Free Radic Biol Med 1999;27:951–65.[Medline]
- Imai H, Hirao F, Sakamoto T, et al. Early embryonic lethality caused by targeted disruption of the mouse PHGPx gene. Biochem Biophys Res Commun 2003;305:278–86.[Medline]
- Yant LJ, Ran Q, Rao L, et al. The selenoprotein GPX4 is essential for mouse development and protects from radiation and oxidative damage insults. Free Radic Biol Med 2003;34:496–502.[Medline]
- Ran Q, Gu M, Van Remmen H, Strong R, Roberts JL, Richardson A. Glutathione peroxidase 4 protects cortical neurons from oxidative injury and amyloid toxicity. J Neurosci Res 2006;84:202–8.[Medline]
- Ran Q, Liang H, Gu M, et al. Transgenic mice overexpressing glutathione peroxidase 4 are protected against oxidative stress-induced apoptosis. J Biol Chem 2004;279:55137–46.[Abstract/Free Full Text]
- Foresta C, Flohe L, Garolla A, Roveri A, Ursini F, Maiorino M. Male fertility is linked to the selenoprotein phospholipid hydroperoxide glutathione peroxidase. Biol Reprod 2002;67:967–71.[Abstract/Free Full Text]
- Imai H, Suzuki K, Ishizaka K, et al. Failure of the expression of phospholipid hydroperoxide glutathione peroxidase in the spermatozoa of human infertile males. Biol Reprod 2001;64:674–83.[Abstract/Free Full Text]
- Diaconu M, Tangat Y, Bohm D, et al. Failure of phospholipid hydroperoxide glutathione peroxidase expression in oligoasthenozoospermia and mutations in the PHGPx gene. Andrologia 2006;38:152–7.[Medline]
- Villette S, Kyle JA, Brown KM, et al. A novel single nucleotide polymorphism in the 3' untranslated region of human glutathione peroxidase 4 influences lipoxygenase metabolism. Blood Cells Mol Dis 2002;29:174–8.[Medline]
- Bermano G, Pagmantidis V, Holloway N, et al. Evidence that a polymorphism within the 3'UTR of glutathione peroxidase 4 is functional and associated with susceptibility to colorectal cancer. Genes Nutr 2007;2:225–32.
- Meplan C, Crosley LK, Nicol F, et al. Genetic polymorphisms in the human selenoprotein P gene determine the response of selenoprotein markers to selenium supplementation in a gender-specific manner (the SELGEN study). FASEB J 2007;21:3063–74 (Epub 2007 May 29).[Abstract/Free Full Text]
- Paglia DE, Valentine WN. Studies on the quantitative and qualitative characterization of erythrocyte glutathione peroxidase. J Lab Clin Med 1967;70:158–69.[Medline]
- Brown KM, Pickard K, Nicol F, Beckett GJ, Duthie GG, Arthur JR. Effects of organic and inorganic selenium supplementation on selenoenzyme activity in blood lymphoctyes, granulocytes, platelets and erythrocytes. Clin Sci (Lond) 2000;98:593–9.[Medline]
- Weitzel F, Ursini F, Wendel A. Phospholipid hydroperoxide glutathione peroxidase in various mouse organs during selenium deficiency and repletion. Biochim Biophys Acta 1990;1036:88–94.[Medline]
- Campbell L, Howie F, Arthur JR, Nicol F, Beckett G. Selenium and sulforaphane modify the expression of selenoenzymes in the human endothelial cell line EAhy926 and protect cells from oxidative damage. Nutrition 2007;23:138–44.[Medline]
- Mahtani KR, Brook M, Dean JL, Sully G, Saklatvala J, Clark AR. Mitogen-activated protein kinase p38 controls the expression and posttranslational modification of tristetraprolin, a regulator of tumor necrosis factor alpha mRNA stability. Mol Cell Biol 2001;21:6461–9.[Abstract/Free Full Text]
- Bermano G, Nicol F, Dyer JA, et al. Tissue-specific regulation of selenoenzyme gene expression during selenium deficiency in rats. Biochem J 1995;311(pt 2):425–30.[Medline]
- Berry MJ. Insights into the hierarchy of selenium incorporation. Nat Genet 2005;37:1162–3.[Medline]
- Brigelius-Flohe R, Muller C, Menard J, Florian S, Schmehl K, Wingler K. Functions of GI-GPx: lessons from selenium-dependent expression and intracellular localization. Biofactors 2001;14:101–6.[Medline]
- Curran JE, Jowett JB, Elliott KS, et al. Genetic variation in selenoprotein S influences inflammatory response. Nat Genet 2005;37:1234–41.[Medline]
- Forsberg L, de Faire U, Marklund SL, Andersson PM, Stegmayr B, Morgenstern R. Phenotype determination of a common Pro-Leu polymorphism in human glutathione peroxidase 1. Blood Cells Mol Dis 2000;26:423–6.[Medline]
- Foster CB, Aswath K, Chanock SJ, McKay HF, Peters U. Polymorphism analysis of six selenoprotein genes: support for a selective sweep at the glutathione peroxidase 1 locus (3p21) in Asian populations. BMC Genet 2006;7:56.[Medline]
- Hu YJ, Korotkov KV, Mehta R, et al. Distribution and functional consequences of nucleotide polymorphisms in the 3'-untranslated region of the human Sep15 gene. Cancer Res 2001;61:2307–10.[Abstract/Free Full Text]
- Low SC, Berry MJ. Knowing when not to stop: selenocysteine incorporation in eukaryotes. Trends Biochem Sci 1996;21:203–8.[Medline]
- Wingler K, Bocher M, Flohe L, Kollmus H, Brigelius-Flohe R. mRNA stability and selenocysteine insertion sequence efficiency rank gastrointestinal glutathione peroxidase high in the hierarchy of selenoproteins. Eur J Biochem 1999;259:149–57.[Medline]
- Xu XM, Mix H, Carlson BA, et al. Evidence for direct roles of two additional factors, SECp43 and soluble liver antigen, in the selenoprotein synthesis machinery. J Biol Chem 2005;280:41568–75.[Abstract/Free Full Text]
- Bermano G, Arthur JR, Hesketh JE. Role of the 3' untranslated region in the regulation of cytosolic glutathione peroxidase and phospholipid-hydroperoxide glutathione peroxidase gene expression by selenium supply. Biochem J 1996;320(pt 3):891–5.[Medline]
- Low SC, Grundner-Culemann E, Harney JW, Berry MJ. SECIS-SBP2 interactions dictate selenocysteine incorporation efficiency and selenoprotein hierarchy. Embo J 2000;19:6882–90.[Medline]
- Weiss Sachdev S, Sunde RA. Selenium regulation of transcript abundance and translational efficiency of glutathione peroxidase-1 and -4 in rat liver. Biochem J 2001;357:851–8.[Medline]
- Behne D, Hilmert H, Scheid S, Gessner H, Elger W. Evidence for specific selenium target tissues and new biologically important selenoproteins. Biochim Biophys Acta 1988;966:12–21.[Medline]
- Riese C, Michaelis M, Mentrup B, et al. Selenium-dependent pre- and posttranscriptional mechanisms are responsible for sexual dimorphic expression of selenoproteins in murine tissues. Endocrinology 2006;147:5883–92.[Abstract/Free Full Text]
Received for publication August 22, 2007.
Accepted for publication November 8, 2007.
This article has been cited by other articles:

|
 |

|
 |
 
K. Ashton, L. Hooper, L. J Harvey, R. Hurst, A. Casgrain, and S. J Fairweather-Tait
Methods of assessment of selenium status in humans: a systematic review
Am. J. Clinical Nutrition,
June 1, 2009;
89(6):
2025S - 2039S.
[Abstract]
[Full Text]
[PDF]
|
 |
|