Commentary: Special Issue - Current Understanding of Colorectal and Pancreatic Cancers.

The Commentary on Special Issue- Current Understanding of Colorectal and Pancreatic Cancers provides the reasoning for the selection of the contributions on pancreatic and colorectal cancer, and summarizes in brief the individual topics and comments upon the main outcomes. The current knowledge, contribution of the individual articles within this Special Issue, and arising priorities in the research on pancreatic ductal adenocarcinoma and colorectal cancer are highlighted.

Genetic And Environmental Associations Of Nonspecific Chromosomal Aberrations.

Nonspecific structural chromosomal aberrations (CAs) are found in around 1% of circulating lymphocytes from healthy individuals but the frequency may be higher after exposure to carcinogenic chemicals or radiation. CAs have been used in the monitoring of persons exposed to genotoxic agents and radiation. Previous studies on occupationally exposed individuals have shown associations between the frequency of CAs in peripheral blood lymphocytes and subsequent cancer risk. The cause for CA formation are believed to be unrepaired or insufficiently repaired DNA double-strand breaks or other DNA damage, and additionally telomere shortening. CAs include chromosome (CSAs) and chromatid type aberrations (CTAs). In the present review, we first describe the types of CAs, the conventional techniques used for their detection and some aspects of interpreting the results. We then focus on germline genetic variation in the frequency and type of CAs measured in a genome-wide association study (GWAS) in healthy individuals in relation to occupational and smoking-related exposure compared to non-exposed referents. The associations (at p<10-5) on 1473 healthy individuals were broadly classified in candidate genes from functional pathways related to DNA damage response/repair, including PSMA1, UBR5, RRM2B, PMS2P4, STAG3L4, BOD1, COPRS and FTO; another group included genes related to apoptosis, cell proliferation, angiogenesis and tumorigenesis, COPB1, NR2C1, COPRS, RHOT1, ITGB3, SYK, and SEMA6A; a third small group mapped to genes KLF7, SEMA5A and ITGB3 which were related to autistic traits, known to manifest frequent CAs. Dedicated studies on 153 DNA repair genes showed associations for some 30 genes, expression of which could be modified by the implicated variants. We finally point out that monitoring of CAs is so far the only method of assessing cancer risk in healthy human populations, and the use of the technology should be made more attractive by developing automated performance steps and incorporating artificial intelligence methods into the scoring.

Polymorphisms in transcription factor binding sites and enhancer regions and pancreatic ductal adenocarcinoma risk.

Genome-wide association studies (GWAS) are a powerful tool for detecting variants associated with complex traits and can help risk stratification and prevention strategies against pancreatic ductal adenocarcinoma (PDAC). However, the strict significance threshold commonly used makes it likely that many true risk loci are missed. Functional annotation of GWAS polymorphisms is a proven strategy to identify additional risk loci. We aimed to investigate single-nucleotide polymorphisms (SNP) in regulatory regions [transcription factor binding sites (TFBSs) and enhancers] that could change the expression profile of multiple genes they act upon and thereby modify PDAC risk. We analyzed a total of 12,636 PDAC cases and 43,443 controls from PanScan/PanC4 and the East Asian GWAS (discovery populations), and the PANDoRA consortium (replication population). We identified four associations that reached study-wide statistical significance in the overall meta-analysis: rs2472632(A) (enhancer variant, OR 1.10, 95%CI 1.06,1.13, p = 5.5 × 10-8), rs17358295(G) (enhancer variant, OR 1.16, 95%CI 1.10,1.22, p = 6.1 × 10-7), rs2232079(T) (TFBS variant, OR 0.88, 95%CI 0.83,0.93, p = 6.4 × 10-6) and rs10025845(A) (TFBS variant, OR 1.88, 95%CI 1.50,1.12, p = 1.32 × 10-5). The SNP with the most significant association, rs2472632, is located in an enhancer predicted to target the coiled-coil domain containing 34 oncogene. Our results provide new insights into genetic risk factors for PDAC by a focused analysis of polymorphisms in regulatory regions and demonstrating the usefulness of functional prioritization to identify loci associated with PDAC risk.

Exploring the Neandertal legacy of pancreatic ductal adenocarcinoma risk in Eurasians.

BACKGROUND: The genomes of present-day non-Africans are composed of 1-3% of Neandertal-derived DNA as a consequence of admixture events between Neandertals and anatomically modern humans about 50-60 thousand years ago. Neandertal-introgressed single nucleotide polymorphisms (aSNPs) have been associated with modern human disease-related traits, which are risk factors for pancreatic ductal adenocarcinoma (PDAC), such as obesity, type 2 diabetes, and inflammation. In this study, we aimed at investigating the role of aSNPs in PDAC in three Eurasian populations. RESULTS: The high-coverage Vindija Neandertal genome was used to select aSNPs in non-African populations from 1000 Genomes project phase 3 data. Then, the association between aSNPs and PDAC risk was tested independently in Europeans and East Asians, using existing GWAS data on more than 200 000 individuals. We did not find any significant associations between aSNPs and PDAC in samples of European descent, whereas, in East Asians, we observed that the Chr10p12.1-rs117585753-T allele (MAF = 10%) increased the risk to develop PDAC (OR = 1.35, 95%CI 1.19-1.54, P = 3.59 × 10-6), with a P-value close to a threshold that takes into account multiple testing. CONCLUSIONS: Our results show only a minimal contribution of Neandertal SNPs to PDAC risk.

The dynamics of telomere length in primary and metastatic colorectal cancer lesions.

Telomeric sequences, the structures comprised of hexanucleotide repeats and associated proteins, play a pivotal role in chromosome end protection and preservation of genomic stability. Herein we address telomere length (TL) dynamics in primary colorectal cancer (CRC) tumour tissues and corresponding liver metastases. TL was measured by multiplex monochrome real-time qPCR in paired samples of primary tumours and liver metastases along with non-cancerous reference tissues obtained from 51 patients diagnosed with metastatic CRC. Telomere shortening was observed in the majority of primary tumour tissues compared to non-cancerous mucosa (84.1%, p < 0.0001). Tumours located within the proximal colon had shorter TL than those in the rectum (p < 0.05). TL in liver metastases was not significantly different from that in primary tumours (p = 0.41). TL in metastatic tissue was shorter in the patients diagnosed with metachronous liver metastases than in those diagnosed with synchronous liver metastases (p = 0.03). The metastatic liver lesions size correlated with the TL in metastases (p < 0.05). Following the neoadjuvant treatment, the patients with rectal cancer had shortened telomeres in tumour tissue than prior to the therapy (p = 0.01). Patients with a TL ratio between tumour tissue and the adjacent non-cancerous mucosa of ≥ 0.387 were associated with increased overall survival (p = 0.01). This study provides insights into TL dynamics during progression of the disease. The results show TL differences in metastatic lesions and may help in clinical practice to predict the patient's prognosis.

A Fecal MicroRNA Signature by Small RNA Sequencing Accurately Distinguishes Colorectal Cancers: Results From a Multicenter Study.

BACKGROUND & AIMS: Fecal tests currently used for colorectal cancer (CRC) screening show limited accuracy in detecting early tumors or precancerous lesions. In this respect, we comprehensively evaluated stool microRNA (miRNA) profiles as biomarkers for noninvasive CRC diagnosis. METHODS: A total of 1273 small RNA sequencing experiments were performed in multiple biospecimens. In a cross-sectional study, miRNA profiles were investigated in fecal samples from an Italian and a Czech cohort (155 CRCs, 87 adenomas, 96 other intestinal diseases, 141 colonoscopy-negative controls). A predictive miRNA signature for cancer detection was defined by a machine learning strategy and tested in additional fecal samples from 141 CRC patients and 80 healthy volunteers. miRNA profiles were compared with those of 132 tumors/adenomas paired with adjacent mucosa, 210 plasma extracellular vesicle samples, and 185 fecal immunochemical test leftover samples. RESULTS: Twenty-five miRNAs showed altered levels in the stool of CRC patients in both cohorts (adjusted P < .05). A 5-miRNA signature, including miR-149-3p, miR-607-5p, miR-1246, miR-4488, and miR-6777-5p, distinguished patients from control individuals (area under the curve [AUC], 0.86; 95% confidence interval [CI], 0.79-0.94) and was validated in an independent cohort (AUC, 0.96; 95% CI, 0.92-1.00). The signature classified control individuals from patients with low-/high-stage tumors and advanced adenomas (AUC, 0.82; 95% CI, 0.71-0.97). Tissue miRNA profiles mirrored those of stool samples, and fecal profiles of different gastrointestinal diseases highlighted miRNAs specifically dysregulated in CRC. miRNA profiles in fecal immunochemical test leftover samples showed good correlation with those of stool collected in preservative buffer, and their alterations could be detected in adenoma or CRC patients. CONCLUSIONS: Our comprehensive fecal miRNome analysis identified a signature accurately discriminating cancer aimed at improving noninvasive diagnosis and screening strategies.

The PANcreatic Disease ReseArch (PANDoRA) consortium: Ten years' experience of association studies to understand the genetic architecture of pancreatic cancer.

Pancreatic cancer has an incidence that almost matches its mortality. Only a small number of risk factors and 33 susceptibility loci have been identified. so Moreover, the relative rarity of pancreatic cancer poses significant hurdles for research aimed at increasing our knowledge of the genetic mechanisms contributing to the disease. Additionally, the inability to adequately power research questions prevents small monocentric studies from being successful. Several consortia have been established to pursue a better understanding of the genetic architecture of pancreatic cancers. The Pancreatic disease research (PANDoRA) consortium is the largest in Europe. PANDoRA is spread across 12 European countries, Brazil and Japan, bringing together 29 basic and clinical research groups. In the last ten years, PANDoRA has contributed to the discovery of 25 susceptibility loci, a feat that will be instrumental in stratifying the population by risk and optimizing preventive strategies.

Differences in genome, transcriptome, miRNAome, and methylome in synchronous and metachronous liver metastasis of colorectal cancer.

Despite distant metastases being the critical factor affecting patients' survival, they remain poorly understood. Our study thus aimed to molecularly characterize colorectal cancer liver metastases (CRCLMs) and explore whether molecular profiles differ between Synchronous (SmCRC) and Metachronous (MmCRC) colorectal cancer. This characterization was performed by whole exome sequencing, whole transcriptome, whole methylome, and miRNAome. The most frequent somatic mutations were in APC, SYNE1, TP53, and TTN genes. Among the differently methylated and expressed genes were those involved in cell adhesion, extracellular matrix organization and degradation, neuroactive ligand-receptor interaction. The top up-regulated microRNAs were hsa-miR-135b-3p and -5p, and the hsa-miR-200-family while the hsa-miR-548-family belonged to the top down-regulated. MmCRC patients evinced higher tumor mutational burden, a wider median of duplications and deletions, and a heterogeneous mutational signature than SmCRC. Regarding chronicity, a significant down-regulation of SMOC2 and PPP1R9A genes in SmCRC compared to MmCRC was observed. Two miRNAs were deregulated between SmCRC and MmCRC, hsa-miR-625-3p and has-miR-1269-3p. The combined data identified the IPO5 gene. Regardless of miRNA expression levels, the combined analysis resulted in 107 deregulated genes related to relaxin, estrogen, PI3K-Akt, WNT signaling pathways, and intracellular second messenger signaling. The intersection between our and validation sets confirmed the validity of our results. We have identified genes and pathways that may be considered as actionable targets in CRCLMs. Our data also provide a valuable resource for understanding molecular distinctions between SmCRC and MmCRC. They have the potential to enhance the diagnosis, prognostication, and management of CRCLMs by a molecularly targeted approach.

MUC13-miRNA-4647 axis in colorectal cancer: Prospects to identifications of risk factors and clinical outcomes.

MUC13, a transmembrane mucin glycoprotein, is overexpressed in colorectal cancer (CRC), however, its regulation and functions are not fully understood. It has been shown that MUC13 protects colonic epithelial cells from apoptosis. Therefore, studying MUC13 and MUC13-regulated pathways may reveal promising therapeutic approaches for CRC treatment. Growing evidence suggests that microRNAs (miRs) are involved in the development and progression of CRC. In the present study, the MUC13-miR-4647 axis was addressed in association with survival of patients. miR-4647 is predicted in silico to bind to the MUC13 gene and was analyzed by RT-qPCR in 187 tumors and their adjacent non-malignant mucosa of patients with CRC. The impact of previously mentioned genes on survival and migration abilities of cancer cells was validated in vitro. Significantly upregulated MUC13 (P=0.02) in was observed tumor tissues compared with non-malignant adjacent mucosa, while miR-4647 (P=0.05) showed an opposite trend. Higher expression levels of MUC13 (log-rank P=0.05) were associated with worse patient's survival. The ectopic overexpression of studied miR resulted in decreased migratory abilities and worse survival of cells. Attenuated MUC13 expression levels confirmed the suppression of colony forming of CRC cells. In summary, the present data suggested the essential role of MUC13-miR-4647 in patients' survival, and this axis may serve as a novel therapeutic target. It is anticipated MUC13 may hold significant potential in the screening, diagnosis and treatment of CRC.

An investigation of DNA damage and DNA repair in chemical carcinogenesis triggered by small-molecule xenobiotics and in cancer: Thirty years with the comet assay.

In the present review we addressed the determination of DNA damage induced by small-molecule carcinogens, considered their persistence in DNA and mutagenicity in in vitro and in vivo systems over a period of 30 years. The review spans from the investigation of the role of DNA damage in the cascade of chemical carcinogenesis. In the nineties, this concept evolved into the biomonitoring studies comprising multiple biomarkers that not only reflected DNA/chromosomal damage, but also the potential of the organism for biotransformation/elimination of various xenobiotics. Since first years of the new millennium, dynamic system of DNA repair and host susceptibility factors started to appear in studies and a considerable knowledge has been accumulated on carcinogens and their role in carcinogenesis. It was understood that the final biological links bridging the arising DNA damage and cancer onset remain to be elucidated. In further years the community of scientists learnt that cancer is a multifactorial disease evolving over several decades of individual´s life. Moreover, DNA damage and DNA repair are inseparable players also in treatment of malignant diseases, but affect substantially other processes, such as degeneration. Functional monitoring of DNA repair pathways and DNA damage response may cast some light on above aspects. Very little is currently known about the relationship between telomere homeostasis and DNA damage formation and repair. DNA damage/repair in genomic and mitochondrial DNA and crosstalk between these two entities emerge as a new interesting topic.

Genetic Predictors for Fecal Propionate and Butyrate-Producing Microbiome Pathway Are Not Associated with Colorectal Cancer Risk: A Mendelian Randomization Analysis.

BACKGROUND: Mechanistic data indicate the benefit of short-chain fatty acids (SCFA) produced by gut microbial fermentation of fiber on colorectal cancer, but direct epidemiologic evidence is limited. A recent study identified SNPs for two SCFA traits (fecal propionate and butyrate-producing microbiome pathway PWY-5022) in Europeans and showed metabolic benefits. METHODS: We conducted a two-sample Mendelian randomization analysis of the genetic instruments for the two SCFA traits (three SNPs for fecal propionate and nine for PWY-5022) in relation to colorectal cancer risk in three large European genetic consortia of 58,131 colorectal cancer cases and 67,347 controls. We estimated the risk of overall colorectal cancer and conducted subgroup analyses by sex, age, and anatomic subsites of colorectal cancer. RESULTS: We did not observe strong evidence for an association of the genetic predictors for fecal propionate levels and the abundance of PWY-5022 with the risk of overall colorectal cancer, colorectal cancer by sex, or early-onset colorectal cancer (diagnosed at <50 years), with no evidence of heterogeneity or pleiotropy. When assessed by tumor subsites, we found weak evidence for an association between PWY-5022 and risk of rectal cancer (OR per 1-SD, 0.95; 95% confidence intervals, 0.91-0.99; P = 0.03) but it did not surpass multiple testing of subgroup analysis. CONCLUSIONS: Genetic instruments for fecal propionate levels and the abundance of PWY-5022 were not associated with colorectal cancer risk. IMPACT: Fecal propionate and PWY-5022 may not have a substantial influence on colorectal cancer risk. Future research is warranted to comprehensively investigate the effects of SCFA-producing bacteria and SCFAs on colorectal cancer risk.

Complex molecular profile of DNA repair genes in epithelial ovarian carcinoma patients with different sensitivity to platinum-based therapy.

Epithelial ovarian carcinoma (EOC) is known for high mortality due to diagnosis at advanced stages and frequent therapy resistance. Previous findings suggested that the DNA repair system is involved in the therapeutic response of cancer patients and DNA repair genes are promising targets for novel therapies. This study aimed to address complex inter-relations among gene expression levels, methylation profiles, and somatic mutations in DNA repair genes and EOC prognosis and therapy resistance status. We found significant associations of DUT expression with the presence of peritoneal metastases in EOC patients. The high-grade serous EOC subtype was enriched with TP53 mutations compared to other subtypes. Furthermore, somatic mutations in XPC and PRKDC were significantly associated with worse overall survival of EOC patients, and higher FAAP20 expression in platinum-resistant than platinum-sensitive patients was observed. We found higher methylation of RAD50 in platinum-resistant than in platinum-sensitive patients. Somatic mutations in BRCA1 and RAD9A were significantly associated with higher RBBP8 methylation in platinum-sensitive compared to platinum-resistant EOC patients. In conclusion, we discovered associations of several candidate genes from the DNA repair pathway with the prognosis and platinum resistance status of EOC patients, which deserve further validation as potential predictive biomarkers.

MiR-140 leads to MRE11 downregulation and ameliorates oxaliplatin treatment and therapy response in colorectal cancer patients.

Cancer therapy failure is a fundamental challenge in cancer treatment. One of the most common reasons for therapy failure is the development of acquired resistance of cancer cells. DNA-damaging agents are frequently used in first-line chemotherapy regimens and DNA damage response, and DNA repair pathways are significantly involved in the mechanisms of chemoresistance. MRE11, a part of the MRN complex involved in double-strand break (DSB) repair, is connected to colorectal cancer (CRC) patients' prognosis. Our previous results showed that single-nucleotide polymorphisms (SNPs) in the 3' untranslated region (3'UTR) microRNA (miRNA) binding sites of MRE11 gene are associated with decreased cancer risk but with shorter survival of CRC patients, which implies the role of miRNA regulation in CRC. The therapy of colorectal cancer utilizes oxaliplatin (oxalato(trans-l-1,2-diaminocyclohexane)platinum), which is often compromised by chemoresistance development. There is, therefore, a crucial clinical need to understand the cellular processes associated with drug resistance and improve treatment responses by applying efficient combination therapies. The main aim of this study was to investigate the effect of miRNAs on the oxaliplatin therapy response of CRC patients. By the in silico analysis, miR-140 was predicted to target MRE11 and modulate CRC prognosis. The lower expression of miR-140 was associated with the metastatic phenotype (p < 0.05) and poor progression-free survival (odds ratio (OR) = 0.4, p < 0.05). In the in vitro analysis, we used miRNA mimics to increase the level of miR-140 in the CRC cell line. This resulted in decreased proliferation of CRC cells (p < 0.05). Increased levels of miR-140 also led to increased sensitivity of cancer cells to oxaliplatin (p < 0.05) and to the accumulation of DNA damage. Our results, both in vitro and in vivo, suggest that miR-140 may act as a tumor suppressor and plays an important role in DSB DNA repair and, consequently, CRC therapy response.

Monitoring of telomere dynamics in peripheral blood leukocytes in relation to colorectal cancer patients' outcomes.

We investigated the possible associations between leukocyte telomere length, therapy outcomes, and clinicopathological features in patients with colorectal cancer. Additionally, telomerase reverse transcriptase (TERT) expression was evaluated. Telomere length was measured using singleplex qPCR in 478 consecutive leukocyte DNA samples from 198 patients. Blood was drawn at diagnosis prior to any therapy and then at 6-month intervals for 18 months. Following diagnosis, the telomeres gradually shortened during the course of the treatment regardless of the patient's age. The most pronounced decrease was observed 12 months after the diagnosis (p < 0.0001). Based on tumor localization, the decrease in telomere length one year after the diagnosis followed different trajectories (p = 0.03). In patients treated with adjuvant therapy, telomere length correlated with the time elapsed after completion of therapy (p = 0.03). TERT expression did not correlate with the telomere length; however, it was higher in women than men (1.35-fold, 95% CI 1.11-1.65, p = 0.003) and in smokers than non-smokers (1.27-fold, 95% CI 1.01-1.61, p = 0.04). Leukocyte telomere length declines naturally during aging, but the accelerated shortening observed in our patients was age-independent. Telomere length manifestly reflected chemotherapy impact and could be linked to therapy toxicity.

Common variability in oestrogen-related genes and pancreatic ductal adenocarcinoma risk in women.

The incidence of pancreatic ductal adenocarcinoma (PDAC) is different among males and females. This disparity cannot be fully explained by the difference in terms of exposure to known risk factors; therefore, the lower incidence in women could be attributed to sex-specific hormones. A two-phase association study was conducted in 12,387 female subjects (5436 PDAC cases and 6951 controls) to assess the effect on risk of developing PDAC of single nucleotide polymorphisms (SNPs) in 208 genes involved in oestrogen and pregnenolone biosynthesis and oestrogen-mediated signalling. In the discovery phase 14 polymorphisms showed a statistically significant association (P < 0.05). In the replication none of the findings were validated. In addition, a gene-based analysis was performed on the 208 selected genes. Four genes (NR5A2, MED1, NCOA2 and RUNX1) were associated with PDAC risk, but only NR5A2 showed an association (P = 4.08 × 10-5) below the Bonferroni-corrected threshold of statistical significance. In conclusion, despite differences in incidence between males and females, our study did not identify an effect of common polymorphisms in the oestrogen and pregnenolone pathways in relation to PDAC susceptibility. However, we validated the previously reported association between NR5A2 gene variants and PDAC risk.

TElomeric repeat-containing RNA (TERRA): Physiological functions and relevance in cancer.

Telomeres are complex protective structures located at the ends of linear eukaryotic chromosomes. Their purpose is to prevent genomic instability. Research progress in telomere biology during the past decades has identified a network of telomeric transcripts of which the best-studied is TElomeric Repeat-containing RNA (TERRA). TERRA was shown to be important not only for the preservation of telomere homeostasis and genomic stability but also for the expression of hundreds of genes across the human genome. These findings added a new level of complexity to telomere biology. Herein we provide insights on the telomere transcriptome, its relevance for proper telomere function, and its implications in human pathology. We also discuss possible clinical opportunities of exosomal telomere transcripts detection as a biomarker in cancer precision medicine.

Discovery of Long Non-Coding RNA MALAT1 Amplification in Precancerous Colorectal Lesions.

A colorectal adenoma, an aberrantly growing tissue, arises from the intestinal epithelium and is considered as precursor of colorectal cancer (CRC). In this study, we investigated structural and numerical chromosomal aberrations in adenomas, hypothesizing that chromosomal instability (CIN) occurs early in adenomas. We applied array comparative genomic hybridization (aCGH) to fresh frozen colorectal adenomas and their adjacent mucosa from 16 patients who underwent colonoscopy examination. In our study, histologically similar colorectal adenomas showed wide variability in chromosomal instability. Based on the obtained results, we further stratified patients into four distinct groups. The first group showed the gain of MALAT1 and TALAM1, long non-coding RNAs (lncRNAs). The second group involved patients with numerous microdeletions. The third group consisted of patients with a disrupted karyotype. The fourth group of patients did not show any CIN in adenomas. Overall, we identified frequent losses in genes, such as TSC2, COL1A1, NOTCH1, MIR4673, and GNAS, and gene gain containing MALAT1 and TALAM1. Since long non-coding RNA MALAT1 is associated with cancer cell metastasis and migration, its gene amplification represents an important event for adenoma development.

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

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.