Next Article in Journal
Relationship of Daily Coffee Intake with Vascular Function in Patients with Hypertension
Next Article in Special Issue
Advances in the Study of the Mechanism by Which Selenium and Selenoproteins Boost Immunity to Prevent Food Allergies
Previous Article in Journal
Influences of Vitamin D and Iron Status on Skeletal Muscle Health: A Narrative Review
Previous Article in Special Issue
Association between Selenium Status and Chronic Kidney Disease in Middle-Aged and Older Chinese Based on CHNS Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Associations of Selenoprotein Genetic Variants with the Risks of Colorectal Adenoma and Colorectal Cancer: Case–Control Studies in Irish and Czech Populations

1
Cancer Biology and Therapeutics Group, School of Biomolecular and Biomedical Science, UCD Conway Institute, University College Dublin, D04 V1W8 Dublin, Ireland
2
Institute of Biology and Medical Genetics, First Faculty of Medicine, Charles University, 121 08 Prague, Czech Republic
3
Biomedical Centre, Faculty of Medicine in Pilsen, Charles University, 323 00 Pilsen, Czech Republic
4
Institute of Experimental Medicine ASCR, 142 20 Prague, Czech Republic
5
Department of Surgery, First Faculty of Medicine, Charles University and Thomayer Hospital, 121 08 Prague, Czech Republic
6
Trinity Translational Medicine Institute, Department of Surgery, Trinity College Dublin and St. James’s Hospital, D08 NHY1 Dublin, Ireland
7
Centre for Colorectal Disease, St. Vincent’s University Hospital, D04 T6F4 Dublin, Ireland
*
Author to whom correspondence should be addressed.
Nutrients 2022, 14(13), 2718; https://doi.org/10.3390/nu14132718
Submission received: 9 June 2022 / Accepted: 27 June 2022 / Published: 29 June 2022
(This article belongs to the Special Issue Dietary Selenium Intake and Human Health)

Abstract

:
Background: Selenium manifests its biological effects through its incorporation into selenoproteins, which play several roles in countering oxidative and inflammatory responses implicated in colorectal carcinogenesis. Selenoprotein genetic variants may contribute to colorectal cancer (CRC) development, as we previously observed for SNP variants in a large European prospective study and a Czech case–control cohort. Methods: We tested if significantly associated selenoprotein gene SNPs from these studies were also associated with CRC risk in case–control studies from Ireland (colorectal neoplasia, i.e., cancer and adenoma cases: 450, controls: 461) and the Czech Republic (CRC cases: 718, controls: 646). Genotyping of 23 SNPs (20 in the Irish and 13 in the Czechs) was performed by competitive specific allele-specific PCR (KASPar). Multivariable adjusted logistic regression was used to assess the associations with CRC development. Results: We found significant associations with an increased CRC risk for rs5859 (SELENOF) and rs2972994 (SELENOP) in the Irish cohort but only with rs4802034 (SELENOV) in the Czechs. Significant associations were observed for rs5859 (SELENOF), rs4659382 (SELENON), rs2972994 (SELENOP), rs34713741 (SELENOS), and the related Se metabolism gene variant rs2275129 (SEPHS1) with advanced colorectal neoplasia development. However, none of these findings retained significance after multiple testing corrections. Conclusions: Several SNPs previously associated with CRC risk were also associated with CRC or colorectal neoplasia development in either the Irish or Czech cohorts. Selenoprotein gene variation may modify CRC risk across diverse European populations, although the specific variants may differ.

1. Introduction

Among tumors across Europe, colorectal cancer (CRC) is the second leading type of cancer for both incidence and mortality rates [1]. These rates in Ireland and also the Czech Republic are higher than the EU average [1]. To markedly reverse these rates, it is vital to improve prevention and early detection, as most CRCs are still diagnosed at advanced stages.
The marked geographic variation worldwide in the dietary availability of the essential micronutrient selenium (Se) results in a lower Se status for many populations, including those across Europe, compared with much of North America [2,3]. The suboptimal intake of Se has been associated with an increased risk of several major diseases [2,3,4], including colorectal neoplasia (polyps and tumors) [5,6].
Se exerts its physiological functions and potential anti-carcinogenic properties through its incorporation into selenoproteins as the amino acid selenocysteine. The human selenoproteome is encoded by 25 separate genes that play key roles in cell protection from oxidative stress, redox control, and inflammatory responses [7,8]. The level of Se intake in Europe and other parts of the world is not adequate for the full expression of these protective selenoproteins [3,9], and this is thought to be the major mechanism through which Se can influence CRC development risk [5].
In all, around 50% of selenoprotein genes in populations from Asia, North America, and Europe harbor single nucleotide polymorphisms (SNPs) proposed to be associated with colorectal adenoma (CRA) development and CRC risk and recurrence risk, as well as impacting survival outcomes [10,11,12,13,14,15,16]. Among these, Méplan and colleagues (2010) first reported that genetic variants known to affect the functional activity of selenoproteins in the GPX4, SELENOP, and SELENOS genes were associated with CRC risk in a case–control study from the Czech Republic [11], a European country of known low Se status [17]. The largest study assessing common variation in 154 Se-related genes was conducted by Fedirko and colleagues [10] in a multi-center, European cohort study of 1420 cases and 1420 controls. Here, it was reported that several SNPs in selenoprotein and Se metabolic pathway genes may affect CRC development risk [10], alone or in combination with suboptimal Se status, ascertained by levels of Se and the selenoprotein P (SELENOP) Se transport protein. However, other studies in generally replete Se status populations have reported more equivocal findings. A prospective study from the Women’s Health Initiative (WHI) in the USA reported no major association with CRC risk for SNPs in five selenoproteins (GPX1-4 and SELENOP) [18]. However, selenoprotein variants (in several pathway genes) were associated with CRC risk and survival outcomes in a large study of a USA population with generally adequate dietary Se intakes [16]. Thus, although no individual selenoprotein SNP variant (to our knowledge) has been confirmed in CRC–GWAS (reviewed by [14,15]), this probably reflects the complex interaction between Se intake and the genetic and environmental factors implicated in the etiology of CRC [19].
In this replication study, we analyzed the associations with CRA and CRC risk of several SNPs in selenoprotein genes (including known functional variants and other tagging SNPs) previously reported to be associated with CRC in European populations in studies led by two of the current study’s coauthors [10,11]. This study was conducted in two studies from the small and genetically homogeneous Irish (colorectal neoplasia cases: 450, controls: 461) and Czech (CRC cases: 718, controls: 646) populations, with the latter partially overlapping with a previous report [11].

2. Materials and Methods

2.1. Subjects and Case–Control Designations

This study examined two colorectal neoplasia study cohorts, from the Czech Republic and from Ireland. Selected baseline characteristics of both study populations are summarized in Table 1. The Czech cohort comprised 718 CRC cases and 646 controls collected at several oncology departments in the Czech Republic between September 2003 and December 2016. Controls were those requiring colonoscopy for various gastrointestinal complaints but who showed no evidence of malignant or idiopathic bowel disease and had no history of chronic illness requiring repeated hospitalization [20]. The Irish cohort comprised 450 cases of 251 CRC and 209 advanced adenomas (i.e., adenomas with high-grade dysplasia (HGD), adenomas with at least 20% tubular villous or villous features, all adenomas greater than 10 mm, and the presence of three or more adenomas) [21,22]. The control group consisted of blood donor controls (>50 years at blood collection) and subjects with a positive immunochemical fecal occult blood test (FIT) upon CRC screening but either with no detectable pathology after colonoscopy (n = 461) or with small (1–2 mm) hyperplastic polyps (n = 110). Cases and FIT screening samples were collected from Tallaght University Hospital from 2008 to 2011 from study volunteers within the departments of Gastroenterology and Surgery, and most were recruited patients with a positive FIT test from a CRC screening pilot program [23], in addition to further consecutive CRC patients from the surgery department. Blood samples (for DNA extraction) from an additional 198 incident CRC cases came from the CRC biobank of the UCD Centre for Colorectal Disease, St Vincent University Hospital, Dublin.
Colorectal neoplasia histology was verified by experienced pathologists at each study center. The clinical characteristics of the study cohorts, including age, sex, tumor localization and histology, and pathological tumor–node–metastasis (pTNM) stage (tumor stage, regional lymph node involvement, and distant metastasis) were taken from patients’ medical records. There were limited epidemiological data for the Irish study (restricted to age, sex, and family history). However, all Czech cohort study participants completed demographic/lifestyle questionnaires to provide information regarding their education, lifestyle habits, body mass index, diabetes, family/personal history of cancer, and long-term (at least six consecutive months) drug use. The data for these Czech samples were retrieved and acquired as previously described [24].
Cancer cases were coded using the 10th Revision of the International Classification of Diseases (ICD-10) and the second revision of the International Classification of Disease for Oncology (ICDO-2). Tumors that were overlapping or unspecified (C18.8 and C18.9) were not included in the assignment of proximal and distal colon tumors. Tumors for the proximal and distal colon as well as overlapping or unspecified tumors were combined to define tumors from the whole colon. Cancers of the rectum were defined as tumors occurring at the recto-sigmoid junction (C19) or the rectum (C20). Subjects with anal canal tumors (C21) were excluded from this study.
All patients provided informed consent according to the Helsinki Declaration (2013 update). The study received approval from the Ethical Committees of the St. James’s Hospital and Federated Dublin Voluntary Hospitals (Dublin, Ireland), and St Vincent’s Hospital Healthcare Group (Dublin, Ireland) and from the Ethical Committee of the Medical Faculty and Teaching Hospital in Pilsen, Czech Republic. All biosamples and patient data were coded to protect anonymity; blood samples and the extracted DNA were stored at −80 °C.

2.2. Blood Collection and DNA Extraction

For the Czech cohort, germline DNA was isolated from peripheral leucocytes within 4 weeks of blood sample collection using the standard proteinase K digestion, phenol/chloroform extraction, and ethanol precipitation and stored at −80 °C. Concentrations were then determined using a Nanodrop 2000c Spectrophotometer. Blood samples were taken from the Irish patients from Tallaght University Hospital 1 day before surgery/colonoscopy in EDTA-coated 6 mL VACUTAINER® tubes, (Cruinn Diagnostics, Dublin, Ireland). The buffy coat layer was separated by centrifuging the samples at 2000× g for 10 min (within 4 h of collection), aliquoted into cryovials, and stored at −80 °C. For the CRC samples from St Vincent’s, the DNA was extracted from peripheral blood mononuclear cells that had been stored frozen at −80 °C. All DNA was extracted using the automated Autopure Instrumentation (Qiagen, Hilden, Germany) at St James University Hospital. A NanoDrop 2000c Spectrophotometer (ThermoScientific, Asheville, NC, USA) was used to quantify the extracted DNA.

2.3. Single Nucleotide Polymorphism Genotyping

Participants were genotyped for a total of 23 SNPs in 17 genes, of which 17 were in 13 selenoprotein genes, while the other 6 were in 6 wider Se-related pathway genes (antioxidant and redox reactions, apoptosis, and Wnt signaling). Ten of these (rs2074451 in GPX4, rs445870 in GPX5, rs11705137 in SELENOM, rs11247710 and rs4659382 in SELENON, rs2275129 in SEPHS1, rs11111979 in TXNRD1, rs4645887 in BAX, rs3813498 in FOXO3, and rs17265803 in FRZB) were genotyped in both the Irish and the Czech cohorts. Following previously published data in Czech populations [11], these 10 SNPs were genotyped only in the Irish samples (rs1050450 in GPX1, rs8177447 in GPX3, rs713041 in GPX4, rs5859 in SELENOF, rs2972994 and rs7579 in SELENOP, rs8177426 and rs34713741 in SELENOS, rs9605031 in TXNRD2, and rs4880 in SOD2). Finally, a further 3 SNPs (rs9818758 in GPX1, rs4802034 in SELENOV, and rs7953266 in TXNRD1) were genotyped only in the Czechs, as this CRC case–control cohort was primarily used to validate the results from our 2019 study in the European Prospective Investigation of Cancer and Nutrition (EPIC) [10]. Table S1 lists the allele frequencies for these SNPs observed in both the Czech and Irish study populations. The genotyping was carried out using a robust, competitive allele-specific PCR system (KASPar) by LGC Genomics (Hoddesdon, Hertfordshire, EN11 OEX, UK). Genotyping analysis was blinded to case–control status, and a random 2% of samples as duplicate quality controls were validated to ensure complete concordance. Samples with inconclusive or failed genotype calls were excluded from the analysis; the genotyping success rate was 94.3%. The genotype correlation between the duplicate samples was 100%. Details of allele probe sequences are available upon request.
For cases and controls combined, all SNP genotypes were in Hardy–Weinberg equilibrium (HWE) in both populations, except for GPX4 rs713041, TXNRD1 rs11111979, and FRZB rs17265803 in the Irish cohort and SELENOV rs4802034 in the Czech cohort (Table S1). When separated by case–control status, all studied SNPs were in HWE in the Czech controls, while SELENOV rs4802034 was not in HWE in the Czech cases. Finally, rs713041 and rs11111979 were not in HWE in the Irish cases (with polyps either included or excluded), and SELENOS rs17265803 was not in HWE in the cases with polyps excluded.

2.4. Statistical Analyses

To perform the case–control comparisons in the Irish colorectal neoplasia cohort, we used two different case–control allocation models dependent on the inclusion of small hyperplastic polyps (<2 mm) in the control or pathology (case) group. This is because most of these small hyperplastic polyps have generally been regarded as benign and are found in 20–40% of people over the age of 50 [25], but there is some evidence to suggest that these polyps may be precancerous [26].
The association between individual SNPs and CRC risk was estimated by unconditional logistic regression to calculate the odds ratio (OR) and 95% confidence interval (CI), adjusting for age (as a continuous variable) and sex. Associations in the Czech cohort were also adjusted for the individual’s smoking habits, body mass index (BMI), and alcohol consumption. ORs are presented using the most frequent homozygous genotype as a reference. Three standard genetic analysis models were tested for disease penetrance: common dominant, common recessive, and additive models [27]. Subgroup analyses by sex and anatomical subsite of the colorectum (colon and rectum) were conducted. A two-sided p-value below 0.05 was considered statistically significant. Multiple testing corrections considering the false discovery rate (FDR) were performed using the Benjamin–Hochberg procedure [28]. Analyses were conducted using the R 4.1.2 and R packages (R Foundation for Statistical Computing, Vienna, Austria; http://www.R-project.org/ (accessed on 1 November 2021)) and PLINK (version 1.9) (http://pngu.mgh.harvard.edu/purcell/plink/ (accessed on 10 October 2021)) [29].

3. Results

In total, 23 SNPs in 17 genes were assessed in either the Czech or the Irish cohorts or both. The ORs and 95% CIs for their associations with CRC risk are presented in Table 2 and Table S2, respectively, for each study.
In the larger Czech cohort, a significant association with an increased CRC risk was only observed for SNP rs4802034 in the SELENOV gene (recessive model OR = 2.14; 95% CI: 1.23, 3.72, p = 0.007; additive model OR = 1.44; 95% CI: 1.08, 1.91, p = 0.01; Table 2).
In the Irish cohort, significant associations with higher CRC risks were observed for variants in two selenoprotein genes: rs5859 in SELENOF (dominant model OR = 1.54; 95% CI: 1.10, 2.15, p = 0.012; additive model OR = 1.49; 95% CI: 1.12, 1.98, p = 0.006) and rs2972994 in SELENOP (dominant model OR = 1.44; 95% CI: 1.01, 2.05, p = 0.044), with similar risk estimates for advanced colorectal neoplasia (i.e., tumors and adenomas combined). An association with a decreased disease risk in the dominant genetic model was observed for rs34713741 in SELENOS with CRC (OR = 0.73; 95% CI: 0.53, 0.99, p = 0.049), advanced colorectal neoplasia (p = 0.011), and CRA (i.e., advanced adenomas including those with HGD, p = 0.046), although this variant was borderline for HWE in cases (p = 0.93). Significant associations were also observed with an increased risk of advanced colorectal neoplasia and CRA for rs4659382 in SELENON and rs2275129 in SEPHS1. It was also observed that rs3813498 in FOXO3 was associated (p = 0.011) with a decreased risk of colorectal neoplasia in the dominant genetic model with OR = 0.69 and CI: 0.52, 0.92. The significant findings for colorectal neoplasia categories are shown in Table 3, while all the analyses for CRC are presented in Table S2.
None of the SNP associations in either cohort retained significance following multiple test corrections. Tables S3 and S4 contain the ORs, 95% CIs, p-values, and FDR corrections for the full stratified analyses of the Irish and Czech cohorts, respectively. In the Irish cohort, these analyses were stratified based on neoplasia subtype (i.e., polyp, adenoma, HGD, cancer) with adjustment for age and sex; the Czech cohort was stratified based on sex and cancer subsite (i.e., colon, rectum) with adjustment for age and sex.
There were no major findings in the Czech stratified analyses by major anatomical subsite. Only rs2275129 (SEPSH1) was nominally associated with increased rectal cancer risk, but the raw p-value was just under the significance threshold (p = 0.044), and this may well have been a chance finding due to the multiple comparisons (Table S4).

4. Discussion

In this study, we observed limited associations for colorectal neoplasia risk with polymorphisms in selenoprotein (and related Se metabolism) genes for the studied populations from the Czech Republic and Ireland. Nominally significantly increased CRC risks were observed for rs4802034-SELENOV in the Czech cohort (p = 0.007) and rs5859-SELENOF (p = 0.006) and rs2972994-SELENOP (p = 0.044) in the Irish cohort, while there was a decreased CRC risk observed for rs34713741 in SELENOS in the latter group.
In the stratified analyses by cancer anatomic location and neoplasia type, only rs2275129 in the SEPSH1 gene was associated with an increased risk of rectal cancer in the Czechs. For neoplasia stratification in the Irish cohort, several variants in selenoprotein and Se metabolism genes, i.e., rs5859-SELENOF, rs4659382-SELENON, rs2972994-SELENOP, and rs2275129-SEPHS1, were associated with an increased risk of advanced colorectal neoplasia, whereas SELENOS rs34713741 was associated with a decreased risk of advanced colorectal neoplasia.
While none of the findings retained significance after FDR multiple-testing adjustment, such corrections could be considered over-stringent given that all our analyses were planned a priori, based on our stated hypothesis, and with a modest number of related Se pathway genetic variants previously implicated in CRC development.
The association of rs4802034 in the SELENOV gene with increased CRC risk in the Czechs (p = 0.007) was also observed for the recessive genetic model (p = 0.021) in the EPIC study [10]. Its protein product has a poorly characterized function, particularly regarding its role in influencing colorectal function and the development of CRC [30]. As recently shown in murine models, SELENOV gene variation may alter the response to the endoplasmic reticulum (ER) and oxidative stresses [31], and these processes are also implicated in colorectal carcinogenesis [10].
The current study demonstrated that SELENOF rs5859 was associated with an increased risk of CRC and all advanced colorectal neoplasia in the Irish cohort. SELENOF is thought to play a role in ER protein folding and quality control [32] and has been implicated in prostate cancer and survival after CRC diagnosis in previous studies [33,34]. Two-loci interactions of rs5859 with several different SNPs in SELENOP showed significant associations with altered CRC risk in a previous Czech study by Méplan et al. [11], although the sample sizes for these analyses were small, and the variant was not associated with CRC risk alone. Our observation that rs2972994-SELENOP was associated with an increased risk of CRC and all colorectal neoplasia in the Irish cohort (p = 0.044) aligned with previous studies conducted by Peters et al. [35,36] on advanced CRA. The SELENOP protein is synthesized in the liver and is primarily responsible for the transport of Se to distal tissue, although it also has antioxidant properties in the extracellular space [37]. Genetic variation in the SELENOP gene may contribute to colorectal carcinogenesis by reducing the protein’s antioxidant properties and reducing the colonic tissue bioavailability of Se [38]. Méplan et al. [11] observed a sex-specific differential risk direction for rs2972994-SELENOP, in that rs2972994-SELENOP was not significantly associated with CRC risk for men and women combined, but the authors observed that men with one T allele variant were significantly associated with a higher risk of CRC whereas, in women, the risk of CRC was significantly decreased. This may point to complex interactions between the hierarchy of Se metabolism and sex, contributing to varying effects of this SNP on CRC development risk. The rs2972994 is located on the 44321bp 3′ of STP in proximity to the selenocysteine insertion sequence in the 3′ promoter region of the SELENOP gene [36]. Although the range of functions of the SELENOP gene is not fully understood, polymorphic antisense transcripts influence gene expression by regulating post-transcriptional modifications [28]. The rs2972994-SELENOP variant is expected to interrupt the STAF transcription factor binding site of the Matlnspector-based SELENOP gene [39].
SELENON encodes an ER-localized glycoprotein that is involved in the oxidative protection and regulation of calcium homeostasis; it has thus been implicated in muscle development and function [40]. Hughes et al. [33] found that this protein was significantly (p = 0.001) downregulated in CRC tumor tissues. SELENON rs4659382 was nominally associated (p = 0.035) in the additive model with an increased risk of CRA and with all advanced neoplasia in the Irish cohort, as observed for CRC (dominant model: p = 0.005) in the EPIC study [10].
SEPHS1 encodes selenophosphate synthetase 1. While this is not itself a selenoprotein, it plays a crucial role in selenoprotein synthesis by using selenide and ATP to synthesize selenophosphate, the key Se donor in the production of selenocysteine [41]. SNP variants may decrease the efficiency of this enzyme and reduce selenoprotein production, thus contributing to an increased risk of CRC development. We observed that rs2275129 in SEPHS1 was associated with increased colorectal neoplasia risk in the Irish population, although not CRC risk alone, as was shown in the previous EPIC study [10].
SELENOS plays a vital role in ER stress response and inflammation control. Interestingly, SELENOS-rs34713741 (a functional promoter variant) showed nominal significance for a decreased association with CRA risk in the Irish cohort, with or without small polyps in the control group. Méplan et al. [11] and Li et al. [42] both found this SNP to be associated with an increased risk of CRC. It also appears in our study that this SNP is significantly associated with a decreased risk of advanced colorectal neoplasia in the dominant model only (p = 0.011), whereas the point estimates for the recessive and additive model yielded non-significant inverse associations with CRA risk. This may simply have resulted from type 1 errors as the raw p values were close to the significance cut-off. However, it may also reflect unknown genotyping errors for this SNP or a possible chromosomal crossover event in the tagged variant, which may also explain why the SNP was not in HWE in the Irish cases (with polyps excluded). Alternatively, although this is unlikely, there may be an interaction between this selenoprotein gene variant and the hierarchy of Se metabolism that led to differences in risk directions. SELENOS encodes a protein involved in the unfolded protein response, shuttling misfolded proteins from the ER lumen to the cytosol where they are degraded by ubiquitin-dependent mechanisms. It also influences inflammatory signaling pathways, which reduces ER stress and may also affect colonic cell growth [7,43]. Curan et al. [43] conducted a functional analysis of rs34713741 and provided evidence that the A allele significantly altered the expression of SELENOS after exposure to ER stress agents (p = 6 × 10−5). Furthermore, the downregulation of the SELENOS gene in macrophage cells by short interfering RNA enhanced the IL-6 and TNF-a (inflammatory cytokines) release, depicting its importance as an inflammation mediator. Hence, as for SELENOS, functional variants in this gene may also impact ER and inflammatory stress responses implicated in CRC development.
Among the wider Se pathway genes, FOXO3 plays a vital role as a central transcription factor that mediates various pathological and physiological processes including progression of the cell cycle, apoptosis, and survival [44,45]. We demonstrated that rs3813498 in FOXO3 was associated with a decreased risk of all colorectal neoplasia (OR dominant model: 0.69; CI: 0.52, 0.92, p = 0.011) in the Irish cohort, as previously observed in EPIC (OR: 0.79; CI: 0.67, 0.92, p = 0.003).
Several studies in different regions of the world suggest that genetic variation in different selenoproteins can influence CRC risk [11,12,41] and that this can be modified by Se status [10,14,46]. This association with selenoprotein variants and CRC risk is thought to be more apparent in areas of suboptimal Se availability; one possible explanation for these differences is that a higher intake of Se may counteract the decreased functional efficiency of selenoprotein synthesis (and possible subsequent modified CRC risk) caused by selenoprotein gene variants and the hierarchical regulation of selenoprotein biosynthesis [11,44].
A limitation of our study was the absence of comprehensive Se status data (not available for the Czechs and only available for a subset of the Irish study). Thus, we lacked the power to perform a rigorous SNPxSe status interaction analysis.
The Irish study was analyzed using different case–control allocations, with generally no major impact on the results when small non-advanced polyps were excluded or included from the case or control groups. For significant associations, there seemed to be more consistent associations observed with polyps in the control group or omitted from the analysis. It is not possible to make a definitive conclusion on the effects of polyps on cancer associations since none of these associations were significant following multiple testing corrections. Many CRC case–control studies either are unable to account for the presence of these small polyps, or deliberately include them in the control group due to their common occurrence in controls of an appropriately similar age range to the cases (as for the Czech cohort in this study). A strength of this study is its consideration of these issues and the inference from the results that it is reasonable to include small polyps in control groups, though it may be more robust to exclude them from analyses or consider them a distinct pathology group.

5. Conclusions

Overall, the results of this study add to the finding from our previous studies [10,11] that the risk of CRC may be modified by selenoprotein genotypes, although the context of the contribution of this genetic variation appears complicated and divergent in different study settings. Individuals with those selenoprotein genotypes associated with increased colorectal neoplasia risk may benefit from increasing their Se status.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu14132718/s1, Table S1: Allele frequencies for SNPs observed in both the Czech and Irish study populations. Table S2: Associations of gene variants with colorectal cancer risk in the Irish cohort; Table S3: OR’s, 95% CIs, p-values and FDR corrections for the full stratified analyses of the Irish cohorts.; Table S4: OR’s, 95% CIs, p-values and FDR corrections for the full stratified analyses of the Czech Republic cohort.

Author Contributions

Conceptualization: D.J.H. and P.V. Sample collection: D.J.H., P.V., M.L., V.L., L.V., V.V., J.O., J.B., P.B., B.N., G.D. and K.S. Experimental procedures: D.J.H. Statistical analysis: M.M. and D.J.H. Funding acquisition: D.J.H., P.V., J.O., G.D. and K.S. Writing original draft: D.J.H., M.M. and N.A. Writing, review, and editing: P.V., L.V., V.V., M.L., V.L., J.B., P.B., J.O., B.N., G.D. and K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from the Health Research Board of Ireland (HRA-POR/2013/397) and the International Hundred K+ Cohorts Consortium/Global Genomic Medicine Collaborative (IHCC/G2MC) No.31312020, to D.J.H. Support for this work was also provided by COST Action CA17118 supported by COST (European Cooperation in Science and Technology, www.cost.EU) to D.J.H, V.V., and P.V. P.V. acknowledges support from Czech Science Foundation GACR 20-03997S and VV 22-05942S and from the Czech Health Research Council of the Ministry of Health of the Czech Republic (AZV NV18-03-00199).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki. The study approvals were from the respective Ethical Committees of the St. James’s Hospital and Federated Dublin Voluntary Hospitals (Dublin, Ireland, reference 2007-37-17, May 2010), the St Vincent’s Hospital Healthcare Group (Dublin, Ireland, reference ‘The identification of genetic and molecular factors associated with bowel diseases, treatment sensitivity & resistance’, June 2006), and the Medical Faculty and Teaching Hospital in Pilsen (the Czech Republic, reference NT12025, April 2011).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Niall Swan for the pathology designation of the Irish colorectal neoplasia samples (Departments of Gastroenterology and Surgery, The Adelaide & Meath Hospital).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ferlay, J.; Colombet, M.; Soerjomataram, I.; Parkin, D.M.; Piñeros, M.; Znaor, A.; Bray, F. Cancer statistics for the year 2020: An overview. Int. J. Cancer 2021, 149, 778–789. [Google Scholar] [CrossRef] [PubMed]
  2. Rayman, M.P. Food-chain selenium and human health: Emphasis on intake. Br. J. Nutr. 2008, 100, 254–268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rayman, M.P. Selenoproteins and human health: Insights from epidemiological data. Biochim. Biophys. Acta Gen. Subj. 2009, 1790, 1533–1540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fairweather-Tait, S.J.; Bao, Y.; Broadley, M.; Collings, R.; Ford, D.; Hesketh, J.E.; Hurst, R. Selenium in Human Health and Disease. Antioxid. Redox Signal. 2011, 14, 1337–1383. [Google Scholar] [CrossRef]
  5. Hughes, D.J.; Fedirko, V.; Jenab, M.; Schomburg, L.; Méplan, C.; Freisling, H.; Bueno-De-Mesquita, H.; Hybsier, S.; Becker, N.-P.; Czuban, M.; et al. Selenium status is associated with colorectal cancer risk in the European prospective investigation of cancer and nutrition cohort. Int. J. Cancer 2014, 136, 1149–1161. [Google Scholar] [CrossRef]
  6. Méplan, C. Selenium and Chronic Diseases: A Nutritional Genomics Perspective. Nutrients 2015, 7, 3621–3651. [Google Scholar] [CrossRef]
  7. Bellinger, F.P.; Raman, A.V.; Reeves, M.A.; Berry, M.J. Regulation and function of selenoproteins in human disease. Biochem. J. 2009, 422, 11–22. [Google Scholar] [CrossRef] [Green Version]
  8. Hatfield, D.L.; Tsuji, P.A.; Carlson, B.A.; Gladyshev, V.N. Selenium and selenocysteine: Roles in cancer, health, and development. Trends Biochem. Sci. 2014, 39, 112–120. [Google Scholar] [CrossRef] [Green Version]
  9. Hesketh, J. Nutrigenomics and Selenium: Gene Expression Patterns, Physiological Targets, and Genetics. Annu. Rev. Nutr. 2008, 28, 157–177. [Google Scholar] [CrossRef] [Green Version]
  10. Fedirko, V.; Jenab, M.; Méplan, C.; Jones, J.S.; Zhu, W.; Schomburg, L.; Siddiq, A.; Hybsier, S.; Overvad, K.; Tjønneland, A.; et al. Association of Selenoprotein and Selenium Pathway Genotypes with Risk of Colorectal Cancer and Interaction with Selenium Status. Nutrients 2019, 11, 935. [Google Scholar] [CrossRef] [Green Version]
  11. Meplan, C.; Hughes, D.J.; Pardini, B.; Naccarati, A.; Soucek, P.; Vodickova, L.; Hlavata, I.; Vrana, D.; Vodicka, P.; Hesketh, J.E. Genetic variants in selenoprotein genes increase risk of colorectal cancer. Carcinogenesis 2010, 31, 1074–1079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Sutherland, A.; Kim, D.-H.; Relton, C.; Ahn, Y.-O.; Hesketh, J. Polymorphisms in the selenoprotein S and 15-kDa selenoprotein genes are associated with altered susceptibility to colorectal cancer. Genes Nutr. 2010, 5, 215–223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Haug, U.; Poole, E.; Xiao, L.; Curtin, K.; Duggan, D.; Hsu, L.; Makar, K.W.; Peters, U.; Kulmacz, R.J.; Potter, J.; et al. Glutathione peroxidase tagSNPs: Associations with rectal cancer but not with colon cancer. Genes Chromosom. Cancer 2012, 51, 598–605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Méplan, C.; Hesketh, J. The influence of selenium and selenoprotein gene variants on colorectal cancer risk. Mutagenesis 2012, 27, 177–186. [Google Scholar] [CrossRef] [Green Version]
  15. Méplan, C.; Hesketh, J. Selenium and Cancer: A Story that Should not be Forgotten-Insights from Genomics. Cancer Treat Res. 2013, 159, 145–166. [Google Scholar] [CrossRef]
  16. Slattery, M.L.; Lundgreen, A.; Welbourn, B.; Corcoran, C.; Wolff, R.K. Genetic Variation in Selenoprotein Genes, Lifestyle, and Risk of Colon and Rectal Cancer. PLoS ONE 2012, 7, e37312. [Google Scholar] [CrossRef]
  17. Batáriová, A.; Černá, M.; Spěváčková, V.; Čejchanová, M.; Beneš, B.; Šmíd, J. Whole blood selenium content in healthy adults in the Czech Republic. Sci. Total Environ. 2005, 338, 183–188. [Google Scholar] [CrossRef]
  18. Takata, Y.; Kristal, A.R.; King, I.B.; Song, X.; Diamond, A.M.; Foster, C.B.; Hutter, C.M.; Hsu, L.; Duggan, D.J.; Langer, R.D.; et al. Serum Selenium, Genetic Variation in Selenoenzymes, and Risk of Colorectal Cancer: Primary Analysis from the Women’s Health Initiative Observational Study and Meta-analysis. Cancer Epidemiol. Biomark. Prev. 2011, 20, 1822–1830. [Google Scholar] [CrossRef] [Green Version]
  19. Murphy, N.; Moreno, V.; Hughes, D.J.; Vodicka, L.; Vodicka, P.; Aglago, E.K.; Gunter, M.J.; Jenab, M. Lifestyle and dietary environmental factors in colorectal cancer susceptibility. Mol. Asp. Med. 2019, 69, 2–9. [Google Scholar] [CrossRef]
  20. Landi, D.; Gemignani, F.; Naccarati, A.; Pardini, B.; Vodicka, P.; Vodickova, L.; Novotny, J.; Försti, A.; Hemminki, K.; Canzian, F.; et al. Polymorphisms within micro-RNA-binding sites and risk of sporadic colorectal cancer. Carcinogenesis 2007, 29, 579–584. [Google Scholar] [CrossRef] [Green Version]
  21. Kim, D.H.; Pickhardt, P.J.; Taylor, A.J. Characteristics of Advanced Adenomas Detected at CT Colonographic Screening: Implications for Appropriate Polyp Size Thresholds for Polypectomy Versus Surveillance. Am. J. Roentgenol. 2007, 188, 940–944. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Corley, D.A.; Jensen, C.D.; Marks, A.R.; Zhao, W.K.; Lee, J.K.; Doubeni, C.A.; Zauber, A.G.; de Boer, J.; Fireman, B.H.; Schottinger, J.E.; et al. Adenoma Detection Rate and Risk of Colorectal Cancer and Death. N. Engl. J. Med. 2014, 370, 1298–1306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kelley, L.; Swan, N.; Hughes, D.J. An analysis of the duplicate testing strategy of an Irish immunochemical faecal occult blood test colorectal cancer screening programme. Color. Dis. 2013, 15, e512–e521. [Google Scholar] [CrossRef]
  24. Bendova, P.; Pardini, B.; Susova, S.; Rosendorf, J.; Levy, M.; Skrobanek, P.; Buchler, T.; Kral, J.; Liska, V.; Vodickova, L.; et al. Genetic variations in microRNA-binding sites of solute carrier transporter genes as predictors of clinical outcome in colorectal cancer. Carcinogenesis 2020, 42, 378–394. [Google Scholar] [CrossRef]
  25. Liljegren, A.; Lindblom, A.; Rotstein, S.; Nilsson, B.; Rubio, C.; Jaramillo, E. Prevalence and incidence of hyperplastic polyps and adenomas in familial colorectal cancer: Correlation between the two types of colon polyps. Gut 2003, 52, 1140–1147. [Google Scholar] [CrossRef]
  26. Jass, J.R. Hyperplastic polyps and colorectal cancer: Is there a link? Clin. Gastroenterol. Hepatol. 2004, 2, 1–8. [Google Scholar] [CrossRef]
  27. Clarke, G.M.; Anderson, C.; Pettersson, F.H.; Cardon, L.R.; Morris, A.P.; Zondervan, K. Basic statistical analysis in genetic case-control studies. Nat. Protoc. 2011, 6, 121–133. [Google Scholar] [CrossRef] [Green Version]
  28. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  29. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  30. Peters, K.M.; Carlson, B.A.; Gladyshev, V.N.; Tsuji, P.A. Selenoproteins in colon cancer. Free. Radic. Biol. Med. 2018, 127, 14–25. [Google Scholar] [CrossRef]
  31. Zhang, X.; Xiong, W.; Chen, L.-L.; Huang, J.-Q.; Lei, X.G. Selenoprotein V protects against endoplasmic reticulum stress and oxidative injury induced by pro-oxidants. Free Radic. Biol. Med. 2020, 160, 670–679. [Google Scholar] [CrossRef] [PubMed]
  32. Takeda, Y.; Seko, A.; Hachisu, M.; Daikoku, S.; Izumi, M.; Koizumi, A.; Fujikawa, K.; Kajihara, Y.; Ito, Y. Both isoforms of human UDP-glucose:glycoprotein glucosyltransferase are enzymatically active. Glycobiology 2014, 24, 344–350. [Google Scholar] [CrossRef] [PubMed]
  33. Hughes, D.J.; Kunická, T.; Schomburg, L.; Liška, V.; Swan, N.; Souček, P. Expression of Selenoprotein Genes and Association with Selenium Status in Colorectal Adenoma and Colorectal Cancer. Nutrients 2018, 10, 1812. [Google Scholar] [CrossRef] [Green Version]
  34. Ekoue, D.N.; Ansong, E.; Liu, L.; Macias, V.; Deaton, R.; Lacher, C.; Picklo, M.; Nonn, L.; Gann, P.H.; Kajdacsy-Balla, A.; et al. Correlations of SELENOF and SELENOP genotypes with serum selenium levels and prostate cancer. Prostate 2018, 78, 279–288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Peters, U.; Chatterjee, N.; Hayes, R.B.; Schoen, R.E.; Wang, Y.; Chanock, S.J.; Foster, C.B. Variation in the selenoenzyme genes and risk of advanced distal colorectal adenoma. Cancer Epidemiol. Prev. Biomark. 2008, 17, 1144–1154. [Google Scholar] [CrossRef] [Green Version]
  36. Peters, U.; Chatterjee, N.; Hayes, R.B.; Schoen, R.E.; Chanock, S.J.; Foster, C.B. Genetic variation in the selenoenzymes and risk of advanced colorectal adenoma. Cancer Res. 2007, 67, 2590. [Google Scholar]
  37. Burk, R.F.; Hill, K.E. Selenoprotein P: An Extracellular Protein with Unique Physical Characteristics and a Role in Selenium Homeostasis. Annu. Rev. Nutr. 2005, 25, 215–235. [Google Scholar] [CrossRef]
  38. Calvo, A.; Xiao, N.; Kang, J.; Best, C.J.; Leiva, I.; Emmert-Buck, M.R.; Green, J.E. Alterations in gene expression profiles during prostate cancer progression: Functional correlations to tumorigenicity and down-regulation of selenoprotein-P in mouse and human tumors. Cancer Res. 2002, 62, 5325–5335. [Google Scholar]
  39. Quandt, K.; Frech, K.; Karas, H.; Wingender, E.; Werner, T. Matlnd and Matlnspector: New fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res. 1995, 23, 4878–4884. [Google Scholar] [CrossRef]
  40. Petit, N.; Lescure, A.; Rederstorff, M.; Krol, A.; Moghadaszadeh, B.; Wewer, U.M.; Guicheney, P. Selenoprotein N: An endoplasmic reticulum glycoprotein with an early developmental expression pattern. Hum. Mol. Genet. 2003, 12, 1045–1053. [Google Scholar] [CrossRef] [Green Version]
  41. Tobe, R.; Carlson, B.A.; Huh, J.H.; Castro, N.P.; Xu, X.-M.; Tsuji, P.A.; Lee, S.-G.; Bang, J.; Na, J.-W.; Kong, Y.-Y.; et al. Selenophosphate synthetase 1 is an essential protein with roles in regulation of redox homoeostasis in mammals. Biochem. J. 2016, 473, 2141–2154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Li, J.; Zhu, Y.; Zhou, Y.; Jiang, H.; Chen, Z.; Lu, B.; Shen, X. The SELS rs34713741 Polymorphism Is Associated with Susceptibility to Colorectal Cancer and Gastric Cancer: A Meta-Analysis. Genet. Test. Mol. Biomark. 2020, 24, 835–844. [Google Scholar] [CrossRef] [PubMed]
  43. Curran, J.E.; Jowett, J.; Elliott, K.S.; Gao, Y.; Gluschenko, K.; Wang, J.; Azim, D.M.A.; Cai, G.; Mahaney, M.C.; Comuzzie, A.G.; et al. Genetic variation in selenoprotein S influences inflammatory response. Nat. Genet. 2005, 37, 1234–1241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Schomburg, L.; Schweizer, U. Hierarchical regulation of selenoprotein expression and sex-specific effects of selenium. Biochim. Biophys. Acta Gen. Subj. 2009, 1790, 1453–1462. [Google Scholar] [CrossRef] [PubMed]
  45. Combs, J.G.F. Biomarkers of Selenium Status. Nutrients 2015, 7, 2209–2236. [Google Scholar] [CrossRef] [Green Version]
  46. Bermano, G.; Pagmantidis, V.; Holloway, N.; Kadri, S.; Mowat, N.A.G.; Shiel, R.S.; Arthur, J.R.; Mathers, J.C.; Daly, A.K.; Broom, J.; et al. Evidence that a polymorphism within the 3′UTR of glutathione peroxidase 4 is functional and is associated with susceptibility to colorectal cancer. Genes Nutr. 2007, 2, 225–232. [Google Scholar] [CrossRef] [Green Version]
Table 1. Selected baseline characteristics of the colorectal neoplasia cases and the controls in the Czech and Irish studies.
Table 1. Selected baseline characteristics of the colorectal neoplasia cases and the controls in the Czech and Irish studies.
Study CharacteristicsCzech RepublicRepublic of Ireland
CRC CasesControlsCRC CasesControls
N (%)718 (52.64)646 (47.36)241 (26.5)461 (50.6)
Age64.13 ± 12.0150.73 ± 13.6265.75 ± 11.5052.05 (±11.31)
Sex (males, %)432 (60.18)365 (56.50)156 (64.7)210 (45.6)
Smoking status
Non-smokers (%)414 (57.66)395 (61.15)
Ex-smokers (%)186 (25.91)111 (17.18)
Current smokers (%)118 (16.43)140 (21.67)
Mean BMI (kg/m2) ± SD27.17 ± 5.0226.01 ± 4.62
Alcohol Consumption (%)27.17 ± 5.0226.01 ± 4.62
Yes (%)182 (25.3)202 (31.3)
No (%)128 (17.8)129 (19.9)
Positive Family History of Cancer
Yes (%)275 (38.3)284 (43.9)
No (%)359 (50.0)308 (47.7)
Colorectal cancer
Colon cancer (%)487 (67.83) 129 (53.53)
Rectal cancer (%)231 (32.17) 1 114 (47.30)
T staging n (T1/T2/T3/T4)26/116/406/91 10/9/20/7
N staging n (N0/N1/N2)353/188/87 37/5/4
M staging n (M0/M1)513/117 9/5
Stage (I/II/III/IV)92/399/103/3 13/20/9/5
Other neoplasia (%) 209 (22.9)
Polyps 2 (%) 110 (52.6)
Adenoma 3 (%) 75 (35.9)
Adenomas with high grade dysplasia 4 (%) 24 (11.5)
1 Percentages do not add up to 100% due to missing values. 2 Polyps were generally hyperplastic and less than 2 mm. 3 Adenomas were all >10 mm and/or with villous or tubulovillous components. 4 Adenomas with high-grade dysplasia were advanced adenomas with at least 5% HGD. Abbreviations: CRC, colorectal cancer; TNM staging,: tumor stage, regional lymph node involvement, and distant metastasis; BMI, body mass index.
Table 2. Associations of the selenoprotein gene variants with colorectal cancer risk in the Czech cohort.
Table 2. Associations of the selenoprotein gene variants with colorectal cancer risk in the Czech cohort.
GeneSNPGenotypeCasesControlsOR95% CIp *FDR
GPX1rs9818758 GG/GA/AA497/201/19437/171/19
Dominant (GG vs. GA + AA)7176271.090.84, 1.440.510.86
Recessive (AA vs. GG + GA)7176270.880.41, 1.910.740.94
Additive7176270.950.65, 1.400.810.94
GPX4rs2074451 GG/GT/TT221/350/141191/324/108
Dominant (GG vs. GT + TT)7126231.090.83, 1.430.540.86
Recessive (TT vs. GG + GT)7126231.250.91, 1.730.170.86
Additive7126231.130.94, 1.360.190.86
GPX5rs445870 AA/AG/GG388/271/48336/236/55
Dominant (AA vs. AG + GG)7076271.000.78, 1.290.990.99
Recessive (GG vs. AA + AG)7076270.960.59, 1.540.860.94
Additive7076270.980.77, 1.250.870.94
SELENOMrs11705137 TT/TC/CC224/350/135190/290/137
Dominant (TT vs. TC + CC)7096170.910.69, 1.190.490.86
Recessive (CC vs. TT + TC)7096170.860.63, 1.170.330.86
Additive7096170.910.76, 1.090.290.86
SELENONrs11247710 CC/CG/GG229/370/115209/313/104
Dominant (CC vs. CG + GG)7146261.090.84, 1.430.510.86
Recessive (GG vs. CC + CG)7146260.870.62, 1.220.420.86
Additive*7146260.970.80, 1.180.450.86
SELENONrs4659382 CC/CG/GG447/240/30386/214/31
Dominant (CC vs. CG + GG)7176311.020.79, 1.320.860.94
Recessive (GG vs. CC + CG)7176310.770.42, 1.410.410.86
Additive7176310.890.65, 1.210.450.86
SELENOVrs4802034 CC/CT/TT386/277/48334/262/27
Dominant (CC vs. CT + TT)7116231.040.81, 1.330.790.94
Recessive (TT vs. CC + CT)7116232.141.23, 3.720.0070.19
Additive7116231.441.08, 1.910.010.19
TXNRD1rs11111979 GG/GC/CC203/328/174171/316/139
Dominant (GG vs. GC + CC)7056260.930.70, 1.230.590.86
Recessive (CC vs. GG + GC)7056261.080.81, 1.450.610.86
Additive7056261.000.84, 1.190.980.99
TXNRD1rs7953266 TT/TC/CC256/323/133210/312/105
Dominant (TT vs. TC + CC)7126270.860.66, 1.110.240.86
Recessive (CC vs. TT + TC)7126271.010.72, 1.390.970.99
Additive7126270.950.79, 1.140.610.86
BAXrs4645887TT/TA/AA261/334/119208/322/94
Dominant (TT vs. TA + AA)7146240.890.69, 1.160.390.86
Recessive (AA vs. TT + TA)7146241.290.91, 1.830.160.86
Additive7146241.070.88, 1.300.480.86
FRZBrs17265803 TT/TC/CC521/179/14477/140/10
Dominant (TT vs. TC + CC)7146271.220.92, 1.630.180.86
Recessive (CC vs. TT + TC)7146271.540.61, 3.870.360.86
Additive7146271.270.79, 2.010.310.86
FOXO3rs3813498 TT/TC/CC477/207/28421/180/26
Dominant (TT vs. TC + CC)7126271.080.82, 1.410.590.86
Recessive (CC vs. TT + TC)7126271.080.56, 2.070.830.94
Additive7126271.050.75, 1.460.780.94
SEPHS1rs2275129 GG/GC/CC224/339/150190/326/112
Dominant (GG vs. GC + CC)7136280.930.71, 1.230.620.86
Recessive (CC vs. GG + GC)7136281.360.99, 1.880.060.78
Additive7136281.110.92, 1.330.280.86
SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval. * p-value obtained from logistic regression model (additive, dominant, recessive), adjusted by age, sex, smoking, BMI, and alcohol consumption; p-value considered significant at < 0.05 and noted in bold.
Table 3. The significant associations between SNPs and colorectal neoplasia risk for different subsets of the stratified analysis in the Irish cohort.
Table 3. The significant associations between SNPs and colorectal neoplasia risk for different subsets of the stratified analysis in the Irish cohort.
GeneSNPCase/ControlGenotypeCasesControlsOR95%CIp *FDR
SELENOFrs5859Cancer/ControlCC/CT/TT129/93/14239/118/11
Dominant (CC vs. CT + TT)2363681.541.10, 2.150.0120.50
Recessive (TT vs. CT + CC)2363682.050.91, 4.590.0820.74
Additive2363681.491.12, 1.980.0060.50
SELENOFrs5859Cancer/Control + PolypCC/CT/TT129/93/14305/156/16
Dominant (CC vs. CT + TT)2364771.471.07, 2.020.0170.50
Recessive (TT vs. CT + CC)2364771.820.87, 3.790.1110.79
Additive2364771.421.09, 1.860.0100.50
SELENOFrs5859Adenoma + HGD + Cancer/ControlCC/CT/TT185/129/20239/118/11
Dominant (CC vs. CT + TT)3343681.421.07, 1.880.0160.50
Recessive (TT vs. CT + CC)3343681.940.94, 4.000.0720.71
Additive3343681.391.09, 1.780.0080.50
SELENOFrs5859Adenoma + HGD + Cancer/Control + PolypCC/CT/TT185/129/20305/156/16
Dominant (CC vs. CT + TT)3344771.431.07, 1.900.0140.50
Recessive (TT vs. CT + CC)3344771.840.94, 3.590.0080.50
Additive3344771.391.09, 1.770.0070.50
SELENONrs4659382Adenoma + HGD + Cancer/ControlTT/TA/AA248/164/30133/76/5
Dominant (TT vs. TA +AA)4422141.280.92, 1.790.1420.68
Recessive (AA vs. TA + TT)4422143.041.16, 7.960.0230.80
Additive4422141.361.02, 1.800.0350.57
SELENONrs4659382Adenoma + HGD/ControlTT/TA/AA49/41/7133/76/5
Dominant (TT vs. TA +AA)972141.610.99, 2.610.0540.69
Recessive (AA vs. TA + TT)972143.251.01, 10.50.0490.69
Additive972141.641.08, 2.480.0190.50
SELENONrs4659382Polyp + Adenoma + HGD/ControlTT/TA/AA110/80/16133/76/5
Dominant (TT vs. TA +AA)2062141.430.97, 2.110.0690.50
Recessive (AA vs. TA + TT)2062143.521.26, 9.790.0160.72
Additive2062141.511.08, 2.100.0150.50
SELENOPrs2972994Cancer/Control + PolypGG/GA/AA58/123/81148/221/102
Dominant (GG vs. GA + AA)2404711.441.01, 2.050.0440.69
Recessive (AA vs. GA + GG2404711.180.82, 1.700.3780.89
Additive2404711.220.98, 1.511.2270.69
SELENOPrs2972994Adenoma + HGD + Cancer/Control + PolypGG/GA/AA85/136/87148/221/102
Dominant (GG vs. GA + AA)3084711.371.00, 1.870.0490.69
Recessive (AA vs. GA + GG3084711.250.89, 1.730.1820.80
Additive3084711.221.00, 1.490.0430.69
SELENOSrs34713741Cancer/Control + PolypCC/CT/TT136/77/22236/207/29
Dominant (CC vs. CT + TT)2354720.730.53, 0.990.0490.69
Recessive (TT vs. CT + CC)2354721.580.88, 2.810.1220.80
Additive2354720.890.69, 1.140.3520.89
SELENOSrs34713741Adenoma+HGD/ControlCC/CT/TT59/31/5183/159/20
Dominant (CC vs. CT + TT)953620.620.39, 0.990.0460.75
Recessive (TT vs. CT + CC)953620.950.35, 2.600.9210.99
Additive953620.710.48, 1.050.0880.69
SELENOSrs34713741Adenoma + HGD/Control + PolypCC/CT/TT59/31/5236/207/29
Dominant (CC vs. CT + TT)954720.610.39, 0.960.0320.71
Recessive (TT vs. CT + CC)954720.850.32, 2.250.7410.94
Additive954720.690.47, 1.010.0600.68
SEPSH1rs2275129Adenoma + HGD + Cancer/ControlGG/GC/CC78/136/8756/116/44
Dominant (GG vs. GC + CC)3012160.960.66, 1.400.8310.85
Recessive (CC vs. GC + GG)3012161.521.02, 2.250.0390.99
Additive3012161.140.91, 1.430.2630.69
SEPSH1rs2275129Adenoma + HGD + Cancer/Control + PolypGG/GC/CC78/136/8787/165/71
Dominant (GG vs. GC + CC)3013231.050.74, 1.510.7720.80
Recessive (CC vs. GC + GG)3013231.441.00, 2.070.0470.95
Additive3013231.170.94, 1.450.1690.69
SEPSH1rs2275129Polyp + Adenoma + HGD/ControlGG/GC/CC51/91/5956/116/44
Dominant (GG vs. GC + CC)2012161.030.66, 1.590.8970.80
Recessive (CC vs. GC + GG)2012161.621.04, 2.550.0340.99
Additive2012161.210.92, 1.580.1710.68
SEPSH1rs2275129Adenoma + HGD/ControlGG/GC/CC20/42/3256/116/44
Dominant (GG vs. GC + CC)942161.290.72, 2.310.3820.68
Recessive (CC vs. GC + GG)942162.021.17, 3.460.0110.89
Additive942161.461.03, 2.080.0350.50
SEPSH1rs2275129Adenoma + HGD/Control + PolypGG/GC/CC20/42/3287/165/71
Dominant (GG vs. GC + CC)943231.360.79, 2.370.2700.68
Recessive (CC vs. GC + GG)943231.831.11, 3.020.0180.85
Additive943231.431.03, 1.980.0340.50
FOXO3rs3813498Adenoma + HGD + Cancer/Control + PolypTT/TC/CC221/103/14212/102/8
Dominant (TT vs. TC + CC)3383220.690.52, 0.920.0110.50
Recessive (CC vs. TC + TT)3383221.700.70, 4.090.2410.83
Additive3383221.070.81, 1.420.6210.92
Polyps were generally hyperplastic and less than 2 mm. Adenomas were all >10 mm and/or with villous or tubulovillous components. Adenomas with high-grade dysplasia were advanced adenomas with at least 5% HGD. SNP, single nucleotide polymorphism; OR, odd ratio; CI, confidence interval; Ca/Co, case/control; HGD, high-grade dysplasia. * p-value obtained from logistic regression model (additive, dominant, recessive), adjusted by age; p-value considered significant at <0.05 and noted in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mukhtar, M.; Ashfield, N.; Vodickova, L.; Vymetalkova, V.; Levy, M.; Liska, V.; Bruha, J.; Bendova, P.; O’Sullivan, J.; Doherty, G.; et al. The Associations of Selenoprotein Genetic Variants with the Risks of Colorectal Adenoma and Colorectal Cancer: Case–Control Studies in Irish and Czech Populations. Nutrients 2022, 14, 2718. https://doi.org/10.3390/nu14132718

AMA Style

Mukhtar M, Ashfield N, Vodickova L, Vymetalkova V, Levy M, Liska V, Bruha J, Bendova P, O’Sullivan J, Doherty G, et al. The Associations of Selenoprotein Genetic Variants with the Risks of Colorectal Adenoma and Colorectal Cancer: Case–Control Studies in Irish and Czech Populations. Nutrients. 2022; 14(13):2718. https://doi.org/10.3390/nu14132718

Chicago/Turabian Style

Mukhtar, Maryam, Niall Ashfield, Ludmila Vodickova, Veronika Vymetalkova, Miroslav Levy, Václav Liska, Jan Bruha, Petra Bendova, Jacintha O’Sullivan, Glen Doherty, and et al. 2022. "The Associations of Selenoprotein Genetic Variants with the Risks of Colorectal Adenoma and Colorectal Cancer: Case–Control Studies in Irish and Czech Populations" Nutrients 14, no. 13: 2718. https://doi.org/10.3390/nu14132718

APA Style

Mukhtar, M., Ashfield, N., Vodickova, L., Vymetalkova, V., Levy, M., Liska, V., Bruha, J., Bendova, P., O’Sullivan, J., Doherty, G., Sheahan, K., Nolan, B., Vodicka, P., & Hughes, D. J. (2022). The Associations of Selenoprotein Genetic Variants with the Risks of Colorectal Adenoma and Colorectal Cancer: Case–Control Studies in Irish and Czech Populations. Nutrients, 14(13), 2718. https://doi.org/10.3390/nu14132718

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop