Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer

Gremlin1: a BMP antagonist with therapeutic potential in Oncology

Gremlins, originating from early 20th-century Western folklore, are mythical creatures known for causing mechanical malfunctions and electronic failures, aptly dubbed "little devils". Analogously, GREM1 acts like a horde of these mischievous entities by antagonizing the bone morphogenetic protein (BMP signaling) pathway or through other non-BMP dependent mechanisms (such as binding to Fibroblast Growth Factor Receptor 1and Epidermal Growth Factor Receptor) contributing to the malignant progression of various cancers. The overexpression of GREM1 promotes tumor cell growth and survival, enhances angiogenesis within the tumor microenvironment, and creates favorable conditions for tumor development and dissemination. Consequently, inhibiting the activity of GREM1 or blocking its interaction with BMP presents a promising strategy for suppressing tumor growth and metastasis. However, the role of GREM1 in cancer remains a subject of debate, with evidence suggesting both oncogenic and tumor-suppressive functions. Currently, several pharmaceutical companies are researching the GREM1 target, with some advancing to Phase I/II clinical trials. This article will provide a detailed overview of the GREM1 target and explore its potential role in cancer therapy.

Genetic risk impacts the association of menopausal hormone therapy with colorectal cancer risk

BACKGROUND

Menopausal hormone therapy (MHT), a common treatment to relieve symptoms of menopause, is associated with a lower risk of colorectal cancer (CRC). To inform CRC risk prediction and MHT risk-benefit assessment, we aimed to evaluate the joint association of a polygenic risk score (PRS) for CRC and MHT on CRC risk.

METHODS

We used data from 28,486 postmenopausal women (11,519 cases and 16,967 controls) of European descent. A PRS based on 141 CRC-associated genetic variants was modeled as a categorical variable in quartiles. Multiplicative interaction between PRS and MHT use was evaluated using logistic regression. Additive interaction was measured using the relative excess risk due to interaction (RERI). 30-year cumulative risks of CRC for 50-year-old women according to MHT use and PRS were calculated.

RESULTS

The reduction in odds ratios by MHT use was larger in women within the highest quartile of PRS compared to that in women within the lowest quartile of PRS (p-value = 2.7 × 10-8). At the highest quartile of PRS, the 30-year CRC risk was statistically significantly lower for women taking any MHT than for women not taking any MHT, 3.7% (3.3%-4.0%) vs 6.1% (5.7%-6.5%) (difference 2.4%, P-value = 1.83 × 10-14); these differences were also statistically significant but smaller in magnitude in the lowest PRS quartile, 1.6% (1.4%-1.8%) vs 2.2% (1.9%-2.4%) (difference 0.6%, P-value = 1.01 × 10-3), indicating 4 times greater reduction in absolute risk associated with any MHT use in the highest compared to the lowest quartile of genetic CRC risk.

CONCLUSIONS

MHT use has a greater impact on the reduction of CRC risk for women at higher genetic risk. These findings have implications for the development of risk prediction models for CRC and potentially for the consideration of genetic information in the risk-benefit assessment of MHT use.

Prioritization of risk genes in colorectal cancer by integrative analysis of multi-omics data and gene networks

Genome-wide association studies (GWASs) have identified over 140 colorectal cancer (CRC)-associated loci; however, target genes at the majority of loci and underlying molecular mechanisms are poorly understood. Here, we utilized a Bayesian approach, integrative risk gene selector (iRIGS), to prioritize risk genes at CRC GWAS loci by integrating multi-omics data. As a result, a total of 105 high-confidence risk genes (HRGs) were identified, which exhibited strong gene dependencies for CRC and enrichment in the biological processes implicated in CRC. Among the 105 HRGs, CEBPB, located at the 20q13.13 locus, acted as a transcription factor playing critical roles in cancer. Our subsequent assays indicated the tumor promoter function of CEBPB that facilitated CRC cell proliferation by regulating multiple oncogenic pathways such as MAPK, PI3K-Akt, and Ras signaling. Next, by integrating a fine-mapping analysis and three independent case-control studies in Chinese populations consisting of 8,039 cases and 12,775 controls, we elucidated that rs1810503, a putative functional variant regulating CEBPB, was associated with CRC risk (OR=0.90, 95%CI=0.86-0.93, P=1.07×10-7). The association between rs1810503 and CRC risk was further validated in three additional multi-ancestry populations consisting of 24,254 cases and 58,741 controls. Mechanistically, the rs1810503 A to T allele change weakened the enhancer activity in an allele-specific manner to decrease CEBPB expression via long-range promoter-enhancer interactions, mediated by the transcription factor, REST, and thus decreased CRC risk. In summary, our study provides a genetic resource and a generalizable strategy for CRC etiology investigation, and highlights the biological implications of CEBPB in CRC tumorigenesis, shedding new light on the etiology of CRC.

BMPR2 as a Novel Predisposition Gene for Hereditary Colorectal Polyposis

BACKGROUND & AIMS

Colorectal cancer (CRC) is one of the most prevalent tumors worldwide, with incidence quickly increasing (particularly in the context of early-onset cases), despite important prevention efforts, mainly in the form of population-wide screening programs. Although many cases present a clear familial component, the current list of hereditary CRC genes leaves a considerable proportion of the cases unexplained.

METHODS

In this work, we used whole-exome sequencing approaches on 19 unrelated patients with unexplained colonic polyposis to identify candidate CRC predisposition genes. The candidate genes were then validated in an additional series of 365 patients. CRISPR-Cas9 models were used to validate BMPR2 as a potential candidate for CRC risk.

RESULTS

We found 8 individuals carrying 6 different variants in the BMPR2 gene (approximately 2% of our cohort of patients with unexplained colonic polyposis). CRISPR-Cas9 models of 3 of these variants showed that the p.(Asn442Thrfs∗32) truncating variant completely abrogated BMP pathway function in a similar way to the BMPR2 knockout. Missense variants p.(Asn565Ser), p.(Ser967Pro) had varying effects on cell proliferation levels, with the former impairing cell control inhibition via noncanonical pathways.

CONCLUSIONS

Collectively, these results support loss-of-function BMPR2 variants as candidates to be involved in CRC germline predisposition.

Insights into radiation carcinogenesis based on dose-rate effects in tissue stem cells

PURPOSE

Increasing epidemiological and biological evidence suggests that radiation exposure enhances cancer risk in a dose-dependent manner. This can be attributed to the 'dose-rate effect,' where the biological effect of low dose-rate radiation is lower than that of the same dose at a high dose-rate. This effect has been reported in epidemiological studies and experimental biology, although the underlying biological mechanisms are not completely understood. In this review, we aim to propose a suitable model for radiation carcinogenesis based on the dose-rate effect in tissue stem cells.

METHODS

We surveyed and summarized the latest studies on the mechanisms of carcinogenesis. Next, we summarized the radiosensitivity of intestinal stem cells and the role of dose-rate in the modulation of stem-cell dynamics after irradiation.

RESULTS

Consistently, driver mutations can be detected in most cancers from past to present, supporting the hypothesis that cancer progression is initiated by the accumulation of driver mutations. Recent reports demonstrated that driver mutations can be observed even in normal tissues, which suggests that the accumulation of mutations is a necessary condition for cancer progression. In addition, driver mutations in tissue stem cells can cause tumors, whereas they are not sufficient when they occur in non-stem cells. For non-stem cells, tissue remodeling induced by marked inflammation after the loss of tissue cells is important in addition to the accumulation of mutations. Therefore, the mechanism of carcinogenesis differs according to the cell type and magnitude of stress. In addition, our results indicated that non-irradiated stem cells tend to be eliminated from three-dimensional cultures of intestinal stem cells (organoids) composed of irradiated and non-irradiated stem cells, supporting the stem-cell competition.

CONCLUSIONS

We propose a unique scheme in which the dose-rate dependent response of intestinal stem cells incorporates the concept of the threshold of stem-cell competition and context-dependent target shift from stem cells to whole tissue. The concept highlights four key issues that should be considered in radiation carcinogenesis: i.e. accumulation of mutations; tissue reconstitution; stem-cell competition; and environmental factors like epigenetic modifications.

circRNA-ZCCHC14 affects the chondrogenic differentiation ability of peripheral blood-derived mesenchymal stem cells by regulating GREM1 through miR-181a

circRNAs play an important role in the progression of osteoarthritis (OA). Therefore, we aimed to reveal the mechanism of action of circRNA-ZCCHC14 in OA. OA animal and cell models were constructed, and clinical samples were collected. The expression of circRNA-ZCCHC14 and miR-181a was detected by RT‒qPCR. The chondrogenic differentiation ability of peripheral blood-derived mesenchymal stem cells (PBMSCs) was detected by Alcian blue staining. The expression of chondrogenic differentiation-related proteins was detected by Western blotting. Double fluorescein experiments verified the targeting relationship of miR-181a with circRNA-ZCCHC14 and GREM1. Upregulation of circRNA-ZCCHC14 was observed in blood, in BMP-2- and TGF-β3-treated PBMSCs from OA patients and in animal models. Knockdown of circRNA-ZCCHC14 promoted the chondrogenic differentiation ability of PBMSCs. circRNA-ZCCHC14 was found to bind to miR-181a and negatively regulate miR-181a expression. Inhibition of miR-181a reversed the promoting effect of circRNA-ZCCHC14 knockdown on the chondrogenic differentiation ability of PBMSCs. GREM1 was identified as a target of miR-181a. Overexpression and knockdown of GREM1 regulated the expression of BMP2, which in turn affected the chondrogenic differentiation ability of PBMSCs, indicating that GREM1 and BMP2 have antagonistic effects and that they jointly regulate the chondrogenic differentiation of PBMSCs. circRNA-ZCCHC14 may promote the chondrogenic differentiation ability of PBMSCs by regulating miR-181a and inhibiting the expression of GREM1.

TCF7L2, CASC8, and GREM1 polymorphism and colorectal cancer in south-eastern Romanian population

Colorectal cancer (CRC) is a heterogeneous disease with an increasing trend and with multiple epigenetic alterations and different molecular features, a major cause of mortality and morbidity. The Wnt/β-Catenin pathway is involved in multiple aspects of cell dynamics, architecture of developing gastrointestinal tissues, and intestinal tissue homeostasis in adults, but its aberrant activity plays an important role in every aspect of colorectal carcinogenesis. The aim of our study was to investigate the association of the TCF7L2 rs7903146, CASC8 rs6983267, and Gremlin1 (GREM1) rs16969681 polymorphism in patients with CRC without other pathologies. A case-control study conducted on 31 patients diagnosed with CRC and 30 healthy controls age and sex-matched with the patients. Real time PCR was used to determine the genotypes of rs7903146, rs698267, rs1696981. We observed no association between rs6983267 and rs16969681 polymorphism and risk of CRC and low association between TCF7L2, rs7903146, polymorphism and risk of CRC. The recessive model of the TCF7L2 rs7903146 had an OR of 1.6 (95% CI 0.058-4.414, P < .05) which means that TT genotype increased the risk and possibility of development of CRC. Our study did not confirm a significant association between TCF7L2 rs7903146, CASC8 rs6983267, and GREM1 rs16969681 with CRC, but emphasizes the possibility of existence of a high risk of CRC development in patients with TT genotype of rs7903146.

Identification of shared genetic architecture between non-alcoholic fatty liver disease and type 2 diabetes: A genome-wide analysis

BACKGROUND

The incidence of complications of non-alcoholic fatty liver disease (NAFLD) and type 2 diabetes (T2D) has been increasing.

METHOD

In order to identify the shared genetic architecture of the two disease phenotypes of NAFLD and T2D, a European population-based GWAS summary and a cross-trait meta-analysis was used to identify significant shared genes for NAFLD and T2D. The enrichment of shared genes was then determined through the use of functional enrichment analysis to investigate the relationship between genes and phenotypes. Additionally, differential gene expression analysis was performed, significant differentially expressed genes in NAFLD and T2D were identified, genes that overlapped between those that were differentially expressed and cross-trait results were reported, and enrichment analysis was performed on the core genes that had been obtained in this way. Finally, the application of a bidirectional Mendelian randomization (MR) approach determined the causal link between NAFLD and T2D.

RESULT

A total of 115 genes were discovered to be shared between NAFLD and T2D in the GWAS analysis. The enrichment analysis of these genes showed that some were involved in the processes such as the decomposition and metabolism of lipids, phospholipids, and glycerophospholipids. Additionally, through the use of differential gene expression analysis, 15 core genes were confirmed to be linked to both T2D and NAFLD. They were correlated with carcinoma cells and inflammation. Furthermore, the bidirectional MR identified a positive causal relationship between NAFLD and T2D.

CONCLUSION

Our study determined the genetic structure shared between NAFLD and T2D, offering a new reference for the genetic pathogenesis and mechanism of NAFLD and T2D comorbidities.

Targeting Stress Erythropoiesis Pathways in Cancer

Cancer-related anemia (CRA) is a common multifactorial disorder that adversely affects the quality of life and overall prognosis in patients with cancer. Safety concerns associated with the most common CRA treatment options, including intravenous iron therapy and erythropoietic-stimulating agents, have often resulted in no or suboptimal anemia management for many cancer patients. Chronic anemia creates a vital need to restore normal erythropoietic output and therefore activates the mechanisms of stress erythropoiesis (SE). A growing body of evidence demonstrates that bone morphogenetic protein 4 (BMP4) signaling, along with glucocorticoids, erythropoietin, stem cell factor, growth differentiation factor 15 (GDF15) and hypoxia-inducible factors, plays a pivotal role in SE. Nevertheless, a chronic state of SE may lead to ineffective erythropoiesis, characterized by the expansion of erythroid progenitor pool, that largely fails to differentiate and give rise to mature red blood cells, further aggravating CRA. In this review, we summarize the current state of knowledge on the emerging roles for stress erythroid progenitors and activated SE pathways in tumor progression, highlighting the urgent need to suppress ineffective erythropoiesis in cancer patients and develop an optimal treatment strategy as well as a personalized approach to CRA management.

Risk Stratification for Early-Onset Colorectal Cancer Using a Combination of Genetic and Environmental Risk Scores: An International Multi-Center Study

BACKGROUND

The incidence of colorectal cancer (CRC) among individuals aged younger than 50 years has been increasing. As screening guidelines lower the recommended age of screening initiation, concerns including the burden on screening capacity and costs have been recognized, suggesting that an individualized approach may be warranted. We developed risk prediction models for early-onset CRC that incorporate an environmental risk score (ERS), including 16 lifestyle and environmental factors, and a polygenic risk score (PRS) of 141 variants.

METHODS

Relying on risk score weights for ERS and PRS derived from studies of CRC at all ages, we evaluated risks for early-onset CRC in 3486 cases and 3890 controls aged younger than 50 years. Relative and absolute risks for early-onset CRC were assessed according to values of the ERS and PRS. The discriminatory performance of these scores was estimated using the covariate-adjusted area under the receiver operating characteristic curve.

RESULTS

Increasing values of ERS and PRS were associated with increasing relative risks for early-onset CRC (odds ratio per SD of ERS = 1.14, 95% confidence interval [CI] = 1.08 to 1.20; odds ratio per SD of PRS = 1.59, 95% CI = 1.51 to 1.68), both contributing to case-control discrimination (area under the curve = 0.631, 95% CI = 0.615 to 0.647). Based on absolute risks, we can expect 26 excess cases per 10 000 men and 21 per 10 000 women among those scoring at the 90th percentile for both risk scores.

CONCLUSIONS

Personal risk scores have the potential to identify individuals at differential relative and absolute risk for early-onset CRC. Improved discrimination may aid in targeted CRC screening of younger, high-risk individuals, potentially improving outcomes.

F-box only protein 9 and its role in cancer

The F-box proteins (FBP), substrate recognition subunit of the SCF (Skp1-Cullin1-F-box protein complex) E3 ligase, play important roles in the ubiquitylation and subsequent degradation of the target proteins from several cellular processes. Disorders of F-box protein-mediated proteolysis lead to human malignancies. FBP plays an important role in many cellular processes, including cell proliferation, cell cycle, apoptosis, migration, invasion, and metastasis, suggesting that it can be associated with tumorigenesis, cancer development and progression. However, the expression and function of FBXO9 (F-box only protein 9) differ in various types of human cancer. Due to the ability to regulate the stability and activity of oncogenes and tumor-suppressor genes, and the physiological functions of many of the F-box proteins remain subtle, further genetic and mechanistic studies will elaborate and help define FBXO9's role. Targeting F-box protein or F-box protein signaling pathways could be an effective strategy for preventing or treating human cancer. This review is presented to summarize the part of FBXO9 in different types of human cancer and its regulation mechanism, and to pave the way to design FBXO9-targeting anticancer therapies.

Morphogen regulation of stem cell plasticity in intestinal regeneration and carcinogenesis

The intestinal epithelium is a tissue with high cell turnover, supported by adult intestinal stem cells. Intestinal homeostasis is underpinned by crypt basal columnar stem cells, marked by expression of the LGR5 gene. However, recent research has demonstrated considerable stem cell plasticity following injury, with dedifferentiation of a range of other intestinal cell populations, induced by a permissive microenvironment in the regenerating mucosa. The regulation of this profound adaptive cell reprogramming response is the subject of current research. There is a demonstrable contribution from disruption of key homeostatic signaling pathways such as wingless-related integration site and bone morphogenetic protein, and an emerging signaling hub role for the mechanoreceptor transducers Yes-associated protein 1/transcriptional coactivator with PDZ-binding motif, negatively regulated by the Hippo pathway. However, a number of outstanding questions remain, including a need to understand how tissues sense damage, and how pathways intersect to mediate dynamic changes in the stem cell population. Better understanding of these pathways, associated functional redundancies, and how they may be both enhanced for recovery of inflammatory diseases, and co-opted in neoplasia development, may have significant clinical implications, and could lead to development of more targeted molecular therapies which target individual stem or stem-like cell populations.

Alterations to Genome Organisation in Stem Cells, Their Differentiation and Associated Diseases

The organisation of the genome in its home, the cell nucleus, is reliant on a number of different aspects to establish, maintain and alter its functional non-random positioning. The genome is dispersed throughout a cell nucleus in specific chromosome territories which are further divided into topologically associated domains (TADs), where regions of the genome from different and the same chromosomes come together. This organisation is both controlled by DNA and chromatin epigenetic modification and the association of the genome with nuclear structures such as the nuclear lamina, the nucleolus and nuclear bodies and speckles. Indeed, sequences that are associated with the first two structures mentioned are termed lamina-associated domains (LADs) and nucleolar-associated domains (NADs), respectively. The modifications and nuclear structures that regulate genome function are altered through a cell's life from stem cell to differentiated cell through to reversible quiescence and irreversible senescence, and hence impacting on genome organisation, altering it to silence specific genes and permit others to be expressed in a controlled way in different cell types and cell cycle statuses. The structures and enzymes and thus the organisation of the genome can also be deleteriously affected, leading to disease and/or premature ageing.

Colorectal cancer risk susceptibility loci in a Swedish population

To search for colorectal cancer (CRC) risk loci, Swedish samples were used for a genome-wide haplotype analysis. A logistic regression model was employed in 2663 CRC cases and 1642 controls in the discovery analysis. Three analyses were done, on all, familial-, and nonfamilial CRC samples and only results with odds ratio (OR) > 1 were analyzed. single nucleotide polymorphism (SNP) analysis did not generate any statistically significant results. Haplotype analysis suggested novel loci, on chromosome 2q36.1 (OR = 1.71, p value = 5.6924 × 10^-8 ) in all CRC samples, chromosome 1q43 (OR = 4.04 p value = 3.24 × 10^-8 ) in familial CRC samples, and two hits in nonfamilial CRC samples, chromosomes 2q36.1 (OR = 1.71 p value = 5.69 × 10^-8 ) and 3p24.3 (OR = 1.62 p value = 6.21 × 10^-9 ). Moreover, one locus on chromosome 20q13.33 was suggested in analyses of all samples, and five more novel loci were suggested on chromosomes 10q25.3, 15q,22.31, 17p11.2, 1p34.2, and 3q24. The haplotypes from the analysis of all samples were replicated in a second study of CRC cases and controls from the same part of Sweden. In summary, using haplotype analysis in Swedish CRC samples, the best hits were novel loci and the locus on chromosomes 2q36.1 and 20q13.33 suggested in the analysis of all samples were confirmed in a second cohort. The ORs were often higher than ORs from published genome-wide association study (GWAS). The study suggested it was possible that a risk locus could involve more than one gene, and that haplotypes could give information on the gene or genes possibly involved in the risk at specific locus.

Repurposing non-invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome

In our previous work, genomic data generated through non-invasive prenatal testing (NIPT) based on low-coverage massively parallel whole-genome sequencing of total plasma DNA of pregnant women in Slovakia was described as a valuable source of population specific data. In the present study, these data were used to determine the population allele frequency of common risk variants located in genes associated with colorectal cancer (CRC) and Lynch syndrome (LS). Allele frequencies of identified variants were compared with six world populations to detect significant differences between populations. Finally, variants were interpreted, functional consequences were searched for and clinical significance of variants was investigated using publicly available databases. Although the present study did not identify any pathogenic variants associated with CRC or LS in the Slovak population using NIPT data, significant differences were observed in the allelic frequency of risk CRC variants previously reported in genome-wide association studies and common variants located in genes associated with LS. As Slovakia is one of the leading countries with the highest incidence of CRC among male patients in the world, there is a need for studies dedicated to investigating the cause of such a high incidence of CRC in Slovakia. The present study also assumed that extensive cross-country data aggregation of NIPT results would represent an unprecedented source of information concerning human genome variation in cancer research.

Non-Lynch Familial and Early-Onset Colorectal Cancer Explained by Accumulation of Low-Risk Genetic Variants

A large proportion of familial and/or early-onset cancer patients do not carry pathogenic variants in known cancer predisposing genes. We aimed to assess the contribution of previously validated low-risk colorectal cancer (CRC) alleles to familial/early-onset CRC (fCRC) and to serrated polyposis. We estimated the association of CRC with a 92-variant-based weighted polygenic risk score (wPRS) using 417 fCRC patients, 80 serrated polyposis patients, 1077 hospital-based incident CRC patients, and 1642 controls. The mean wPRS was significantly higher in fCRC than in controls or sporadic CRC patients. fCRC patients in the highest (20th) wPRS quantile were at four-fold greater CRC risk than those in the middle quantile (10th). Compared to low-wPRS fCRC, a higher number of high-wPRS fCRC patients had developed multiple primary CRCs, had CRC family history, and were diagnosed at age ≥50. No association with wPRS was observed for serrated polyposis. In conclusion, a relevant proportion of mismatch repair (MMR)-proficient fCRC cases might be explained by the accumulation of low-risk CRC alleles. Validation in independent cohorts and development of predictive models that include polygenic risk score (PRS) data and other CRC predisposing factors will determine the implementation of PRS into genetic testing and counselling in familial and early-onset CRC.

Individual and Joint Associations of Genetic Risk and Healthy Lifestyle Score with Colorectal Neoplasms Among Participants of Screening Colonoscopy

Genetic and lifestyle factors contribute to colorectal cancer risk. We investigated their individual and joint associations with various stages of colorectal carcinogenesis. We assessed associations of a polygenic risk score (PRS) and a healthy lifestyle score (HLS) with presence of nonadvanced adenomas and advanced neoplasms among 2,585 participants of screening colonoscopy from Germany. The PRS and HLS individually showed only weak associations with presence of nonadvanced adenomas; stronger associations were observed with advanced neoplasms (ORs, 95% CI, for highest vs. lowest risk tertile: PRS 2.27, 1.78-2.88; HLS 1.96, 1.53-2.51). The PRS was associated with higher odds of advanced neoplasms among carriers of any neoplasms (1.65, 1.23-2.22). Subjects in the highest risk tertile (vs. lowest tertile) of both scores had higher risks for nonadvanced adenomas (1.77, 1.09-2.86), for advanced neoplasms (3.95, 2.53-6.16) and, among carriers of any neoplasms, for advanced versus nonadvanced neoplasms (2.26, 1.31-3.92). Both scores were individually associated with increased risk of nonadvanced adenomas and, much more pronounced, advanced neoplasms. The similarly strong association in relative terms across all levels of genetic risk implies that a healthy lifestyle may be particularly beneficial in those at highest genetic risk, given that the same relative risk reduction in this group would imply a stronger absolute risk reduction. Genetic factors may be of particular relevance for the transition of nonadvanced to advanced adenomas. PREVENTION RELEVANCE: Genetic factors have strong impact on the risk of colorectal neoplasms, which may be reduced by healthy lifestyle. Similarly strong associations in relative terms across all levels of genetic risk imply that a healthy lifestyle may be beneficial due to higher absolute risk reduction in those at highest genetic risk.

Identifying Novel Susceptibility Genes for Colorectal Cancer Risk From a Transcriptome-Wide Association Study of 125,478 Subjects

BACKGROUND AND AIMS

Susceptibility genes and the underlying mechanisms for the majority of risk loci identified by genome-wide association studies (GWAS) for colorectal cancer (CRC) risk remain largely unknown. We conducted a transcriptome-wide association study (TWAS) to identify putative susceptibility genes.

METHODS

Gene-expression prediction models were built using transcriptome and genetic data from the 284 normal transverse colon tissues of European descendants from the Genotype-Tissue Expression (GTEx), and model performance was evaluated using data from The Cancer Genome Atlas (n = 355). We applied the gene-expression prediction models and GWAS data to evaluate associations of genetically predicted gene-expression with CRC risk in 58,131 CRC cases and 67,347 controls of European ancestry. Dual-luciferase reporter assays and knockdown experiments in CRC cells and tumor xenografts were conducted.

RESULTS

We identified 25 genes associated with CRC risk at a Bonferroni-corrected threshold of P < 9.1 × 10-6, including genes in 4 novel loci, PYGL (14q22.1), RPL28 (19q13.42), CAPN12 (19q13.2), MYH7B (20q11.22), and MAP1L3CA (20q11.22). In 9 known GWAS-identified loci, we uncovered 9 genes that have not been reported previously, whereas 4 genes remained statistically significant after adjusting for the lead risk variant of the locus. Through colocalization analysis in GWAS loci, we additionally identified 12 putative susceptibility genes that were supported by TWAS analysis at P < .01. We showed that risk allele of the lead risk variant rs1741640 affected the promoter activity of CABLES2. Knockdown experiments confirmed that CABLES2 plays a vital role in colorectal carcinogenesis.

CONCLUSIONS

Our study reveals new putative susceptibility genes and provides new insight into the biological mechanisms underlying CRC development.

Subversion of Niche-Signalling Pathways in Colorectal Cancer: What Makes and Breaks the Intestinal Stem Cell

The intestinal epithelium fulfils pleiotropic functions in nutrient uptake, waste elimination, and immune surveillance while also forming a barrier against luminal toxins and gut-resident microbiota. Incessantly barraged by extraneous stresses, the intestine must continuously replenish its epithelial lining and regenerate the full gamut of specialized cell types that underpin its functions. Homeostatic remodelling is orchestrated by the intestinal stem cell (ISC) niche: a convergence of epithelial- and stromal-derived cues, which maintains ISCs in a multipotent state. Following demise of homeostatic ISCs post injury, plasticity is pervasive among multiple populations of reserve stem-like cells, lineage-committed progenitors, and/or fully differentiated cell types, all of which can contribute to regeneration and repair. Failure to restore the epithelial barrier risks seepage of toxic luminal contents, resulting in inflammation and likely predisposing to tumour formation. Here, we explore how homeostatic niche-signalling pathways are subverted in tumorigenesis, enabling ISCs to gain autonomy from niche restraints ("ISC emancipation") and transform into cancer stem cells capable of driving tumour initiation, progression, and therapy resistance. We further consider the implications of the pervasive plasticity of the intestinal epithelium for the trajectory of colorectal cancer, the emergence of distinct molecular subtypes, the propensity to metastasize, and the development of effective therapeutic strategies.

Polygenic risk score is a predictor of adenomatous polyps at screening colonoscopy

BACKGROUND

Single nucleotide polymorphism (SNP)-based polygenic risk scoring is predictive of colorectal cancer (CRC) risk. However, few studies have investigated the association of genetic risk score (GRS) with detection of adenomatous polyps at screening colonoscopy.

METHODS

We randomly selected 1769 Caucasian subjects who underwent screening colonoscopy from the Genomic Health Initiative (GHI), a biobank of NorthShore University HealthSystem. Outcomes from initial screening colonoscopy were recorded. Twenty-two CRC risk-associated SNPs were obtained from the Affymetrix™ SNP array and used to calculate an odds ratio (OR)-weighted and population-standardized GRS. Subjects with GRS of < 0.5, 0.5-1.5, and > 1.5 were categorized as low, average and elevated risk.

RESULTS

Among 1,769 subjects, 520 (29%) had 1 or more adenomatous polyps. GRS was significantly higher in subjects with adenomatous polyps than those without; mean (95% confidence interval) was 1.02 (1.00-1.05) and 0.97 (0.95-0.99), respectively, p < 0.001. The association remained significant after adjusting for age, gender, body mass index, and family history, p < 0.001. The detection rate of adenomatous polyps was 10.8%, 29.0% and 39.7% in subjects with low, average and elevated GRS, respectively, p-trend < 0.001. Higher GRS was also associated with early age diagnosis of adenomatous polyps, p < 0.001. In contrast, positive family history was not associated with risk and age of adenomatous polyps.

CONCLUSIONS

GRS was significantly associated with adenomatous polyps in subjects undergoing screening colonoscopy. This result may help in stratifying average risk patients and facilitating personalized colonoscopy screening strategies.

SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Several genome wide association studies of colorectal cancer (CRC) have identified single nucleotide polymorphisms (SNPs) on chromosome 15q13.3 associated with CRC risk. To identify functional variant(s) underlying this association, we investigated SNPs in linkage disequilibrium with the risk‐associated SNP rs4779584 that overlapped regulatory regions/enhancer elements characterized in colon‐related tissues and cells. We identified several SNP‐containing regulatory regions that exhibited enhancer activity in vitro, including one SNP (rs1406389) that correlated with allele‐specific effects on enhancer activity. Deletion of either this enhancer or another enhancer that had previously been reported in this region correlated with decreased expression of GREM1 following CRISPR/Cas9 genome editing. That GREM1 is one target of these enhancers was further supported by an expression quantitative trait loci correlation between rs1406389 and GREM1 expression in the transverse but not sigmoid colon in the Genotype‐Tissue Expression dataset. Taken together, we conclude that the 15q13.3 region contains at least two functional variants that map to distinct enhancers and impact CRC risk through modulation of GREM1 expression.

Multi-omics analysis to identify susceptibility genes for colorectal cancer

Most genetic variants for colorectal cancer (CRC) identified in genome-wide association studies (GWAS) are located in intergenic regions, implying pathogenic dysregulations of gene expression. However, comprehensive assessments of target genes in CRC remain to be explored. We conducted a multi-omics analysis using transcriptome and/or DNA methylation data from the Genotype-Tissue Expression, The Cancer Genome Atlas and the Colonomics projects. We identified 116 putative target genes for 45 GWAS-identified variants. Using summary-data-based Mendelian randomization approach (SMR), we demonstrated that the CRC susceptibility for 29 out of the 45 CRC variants may be mediated by cis-effects on gene regulation. At a cutoff of the Bonferroni-corrected PSMR < 0.05, we determined 66 putative susceptibility genes, including 39 genes that have not been previously reported. We further performed in vitro assays for two selected genes, DIP2B and SFMBT1, and provide functional evidence that they play a vital role in colorectal carcinogenesis via disrupting cell behavior, including migration, invasion and epithelial-mesenchymal transition. Our study reveals a large number of putative novel susceptibility genes and provides additional insight into the underlying mechanisms for CRC genetic risk loci.

Use of Polygenic Risk Scores to Select Screening Intervals After Negative Findings From Colonoscopy

BACKGROUND & AIMS

Polygenic risk scores (PRSs) could help to define starting ages for colorectal cancer (CRC) screening. However, the role of PRS in determining the length of screening interval after negative findings from colonoscopies is unclear. We aimed to evaluate CRC risk according to PRS and time since last negative colonoscopy.

METHODS

We collected data from 3827 cases and 2641 CRC-free controls in a population-based case-control study in Germany. We constructed a polygenic risk scoring system, based on 90 single-nucleotide polymorphisms, associated with risk of CRC in people of European descent. Participants were classified as having low, medium, or high genetic risk according to tertiles of PRSs among controls. Multiple logistic regression models were used to assess CRC risk according to PRS and time since last negative colonoscopy.

RESULTS

Compared to individuals without colonoscopy in the low PRS category, a 42%-85% lower risk of CRC was observed for individuals who had a negative finding from colonoscopy within 10 years. Beyond 10 years after a negative finding from colonoscopy, significantly lower risk only persisted for the low and medium PRS groups, but not for the high PRS group. Adjusted odds ratios were 0.44 (95% CI, 0.29-0.68), 0.51 (95% CI, 0.34-0.77), and 0.85 (95% CI, 0.58-1.23) in the low, medium, and high PRS group, respectively. Within any time interval, risks were lower for distal than for proximal CRCs.

CONCLUSIONS

Based on findings from a population-based case-control study, the recommended 10-year screening interval for colonoscopy may not need to be shortened among people with high PRSs, but could potentially be prolonged for people with low and medium PRSs. Studies are needed to address personalized time intervals for repeat colonoscopies in average-risk screening cohorts.

The Impact of a Comprehensive Risk Prediction Model for Colorectal Cancer on a Population Screening Program

BACKGROUND

In many countries, population colorectal cancer (CRC) screening is based on age and family history, though more precise risk prediction could better target screening. We examined the impact of a CRC risk prediction model (incorporating age, sex, lifestyle, genomic, and family history factors) to target screening under several feasible screening scenarios.

METHODS

We estimated the model's predicted CRC risk distribution in the Australian population. Predicted CRC risks were categorized into screening recommendations under 3 proposed scenarios to compare with current recommendations: 1) highly tailored, 2) 3 risk categories, and 3) 4 sex-specific risk categories. Under each scenario, for 35- to 74-year-olds, we calculated the number of CRC screens by immunochemical fecal occult blood testing (iFOBT) and colonoscopy and the proportion of predicted CRCs over 10 years in each screening group.

RESULTS

Currently, 1.1% of 35- to 74-year-olds are recommended screening colonoscopy and 56.2% iFOBT, and 5.7% and 83.2% of CRCs over 10 years were predicted to occur in these groups, respectively. For the scenarios, 1) colonoscopy was recommended to 8.1% and iFOBT to 37.5%, with 36.1% and 50.1% of CRCs in each group; 2) colonoscopy was recommended to 2.4% and iFOBT to 56.0%, with 13.2% and 76.9% of cancers in each group; and 3) colonoscopy was recommended to 5.0% and iFOBT to 54.2%, with 24.5% and 66.5% of cancers in each group.

CONCLUSIONS

A highly tailored CRC screening scenario results in many fewer screens but more cancers in those unscreened. Category-based scenarios may provide a good balance between number of screens and cancers detected and are simpler to implement.

Secreted BMP antagonists and their role in cancer and bone metastases

Bone morphogenetic proteins (BMPs) are multifunctional secreted cytokines that act in a highly context-dependent manner. BMP action extends beyond the induction of cartilage and bone formation, to encompass pivotal roles in controlling tissue and organ homeostasis during development and adulthood. BMPs signal via plasma membrane type I and type II serine/threonine kinase receptors and intracellular SMAD transcriptional effectors. Exquisite temporospatial control of BMP/SMAD signalling and crosstalk with other cellular cues is achieved by a series of positive and negative regulators at each step in the BMP/SMAD pathway. The interaction of BMP ligand with its receptors is carefully controlled by a diverse set of secreted antagonists that bind BMPs and block their interaction with their cognate BMP receptors. Perturbations in this BMP/BMP antagonist balance are implicated in a range of developmental disorders and diseases, including cancer. Here, we provide an overview of the structure and function of secreted BMP antagonists, and summarize recent novel insights into their role in cancer progression and bone metastasis. Gremlin1 (GREM1) is a highly studied BMP antagonist, and we will focus on this molecule in particular and its role in cancer. The therapeutic potential of pharmacological inhibitors for secreted BMP antagonists for cancer and other human diseases will also be discussed.

Estimation of Absolute Risk of Colorectal Cancer Based on Healthy Lifestyle, Genetic Risk, and Colonoscopy Status in a Population-Based Study

BACKGROUND & AIMS

Estimates of absolute risk of colorectal cancer (CRC) are needed to facilitate communication and better inform the public about the potentials and limits of cancer prevention.

METHODS

Using data from a large population-based case-control study in Germany (Darmkrebs: Chancen der Verhütung durch Screening [DACHS] study, which began in 2003) and population registry data, we calculated 30-year absolute risk estimates for development of CRC based on a healthy lifestyle score (derived from 5 modifiable lifestyle factors: smoking, alcohol consumption, diet, physical activity, and body fatness), a polygenic risk score (based on 90 single-nucleotide polymorphisms), and colonoscopy history.

RESULTS

We analyzed data from 4220 patients with CRC and 3338 individuals without CRC. Adherence to a healthy lifestyle and colonoscopy in the preceding 10 years were associated with a reduced relative risk of CRC in men and women. We observed a higher CRC risk in participants with high or intermediate genetic risk scores. For 50-year-old men and women without a colonoscopy, the absolute risk of CRC varied according to the polygenic risk score and the healthy lifestyle score (men, 3.5%-13.4%; women, 2.5%-10.6%). For 50-year-old men and women with a colonoscopy, the absolute risk of developing CRC was much lower but still varied according to the polygenic risk score and the healthy lifestyle score (men, 1.2%-4.8%; women, 0.9%-4.2%). Among all risk factor profiles, the 30-year absolute risk estimates consistently decreased with adherence to a healthy lifestyle.

CONCLUSIONS

In a population-based study, we found that a colonoscopy can drastically reduce the absolute risk of CRC and that the genetically predetermined risk of CRC can be further reduced by adherence to a healthy lifestyle. Our results show the magnitude of CRC prevention possible through colonoscopy and lifestyle at a predefined genetic risk. This observational study has been registered in the German Clinical Trials Register (DRKS00011793), which is a primary registry in the World Health Organization Registry Network.

Colorectal cancer genetic variants are also associated with serrated polyposis syndrome susceptibility

Background

Serrated polyposis syndrome (SPS) is a clinical entity characterised by large and/ormultiple serrated polyps throughout the colon and increased risk for colorectal cancer (CRC). The basis for SPS genetic predisposition is largely unknown. Common, low-penetrance genetic variants have been consistently associated with CRC susceptibility, however, their role in SPS genetic predisposition has not been yet explored.

Objective

The aim of this study was to evaluate if common, low-penetrance genetic variants for CRC risk are also implicated in SPS genetic susceptibility.

Methods

A case-control study was performed in 219 SPS patients and 548 asymptomatic controls analysing 65 CRC susceptibility variants. A risk prediction model for SPS predisposition was developed.

Results

Statistically significant associations with SPS were found for seven genetic variants (rs4779584-GREM1, rs16892766-EIF3H, rs3217810-CCND2, rs992157-PNKD1/TMBIM1, rs704017-ZMIZ1, rs11196172-TCF7L2, rs6061231-LAMA5). The GREM1 risk allele was remarkably over-represented in SPS cases compared with controls (OR=1.573, 1.21–2.04, p value=0.0006). A fourfold increase in SPS risk was observed when comparing subjects within the highest decile of variants (≥65) with those in the first decile (≤50).

Conclusions

Genetic variants for CRC risk are also involved in SPS susceptibility, being the most relevant ones rs4779584-GREM1, rs16892766-EIF3H and rs3217810-CCND2.

Risk SNP-Mediated Enhancer-Promoter Interaction Drives Colorectal Cancer through Both FADS2 and AP002754.2

Although genome-wide association studies (GWAS) have identified more than 100 colorectal cancer risk loci, most of the biological mechanisms associated with these loci remain unclear. Here we first performed a comprehensive expression quantitative trait loci analysis in colorectal cancer tissues adjusted for multiple confounders to test the determinants of germline variants in established GWAS susceptibility loci on mRNA and long noncoding RNA (lncRNA) expression. Combining integrative functional genomic/epigenomic analyses and a large-scale population study consisting of 6,024 cases and 10,022 controls, we then prioritized rs174575 with a C>G change as a potential causal candidate for colorectal cancer at 11q12.2, as its G allele was associated with an increased risk of colorectal cancer (OR = 1.26; 95% confidence interval = 1.17-1.36; P = 2.57 × 10^-9). rs174575 acted as an allele-specific enhancer to distally facilitate expression of both FADS2 and lncRNA AP002754.2 via long-range enhancer-promoter interaction loops, which were mediated by E2F1. AP002754.2 further activated a transcriptional activator that upregulated FADS2 expression. FADS2, in turn, was overexpressed in colorectal cancer tumor tissues and functioned as a potential oncogene that facilitated colorectal cancer cell proliferation and xenograft growth in vitro and in vivo by increasing the metabolism of PGE2, an oncogenic molecule involved in colorectal cancer tumorigenesis. Our findings represent a novel mechanism by which a noncoding variant can facilitate long-range genome interactions to modulate the expression of multiple genes including not only mRNA, but also lncRNA, which provides new insights into the understanding of colorectal cancer etiology. SIGNIFICANCE: This study provides an oncogenic regulatory circuit among several oncogenes including E2F1, FADS2, and AP002754.2 underlying the association of rs174575 with colorectal cancer risk, which is driven by long-range enhancer-promoter interaction loops. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/9/1804/F1.large.jpg.

Identification of Novel Loci and New Risk Variant in Known Loci for Colorectal Cancer Risk in East Asians

BACKGROUND

Risk variants identified so far for colorectal cancer explain only a small proportion of familial risk of this cancer, particularly in Asians.

METHODS

We performed a genome-wide association study (GWAS) of colorectal cancer in East Asians, including 23,572 colorectal cancer cases and 48,700 controls. To identify novel risk loci, we selected 60 promising risk variants for replication using data from 58,131 colorectal cancer cases and 67,347 controls of European descent. To identify additional risk variants in known colorectal cancer loci, we performed conditional analyses in East Asians.

RESULTS

An indel variant, rs67052019 at 1p13.3, was found to be associated with colorectal cancer risk at P = 3.9 × 10-8 in Asians (OR per allele deletion = 1.13, 95% confidence interval = 1.08-1.18). This association was replicated in European descendants using a variant (rs2938616) in complete linkage disequilibrium with rs67052019 (P = 7.7 × 10-3). Of the remaining 59 variants, 12 showed an association at P < 0.05 in the European-ancestry study, including rs11108175 and rs9634162 at P < 5 × 10-8 and two variants with an association near the genome-wide significance level (rs60911071, P = 5.8 × 10-8; rs62558833, P = 7.5 × 10-8) in the combined analyses of Asian- and European-ancestry data. In addition, using data from East Asians, we identified 13 new risk variants at 11 loci reported from previous GWAS.

CONCLUSIONS

In this large GWAS, we identified three novel risk loci and two highly suggestive loci for colorectal cancer risk and provided evidence for potential roles of multiple genes and pathways in the etiology of colorectal cancer. In addition, we showed that additional risk variants exist in many colorectal cancer risk loci identified previously.

IMPACT

Our study provides novel data to improve the understanding of the genetic basis for colorectal cancer risk.

Performance of individual and joint risk stratification by an environmental risk score and a genetic risk score in a colorectal cancer screening setting

Early detection of colorectal neoplasms can reduce the disease burden of colorectal cancer by timely intervention of individuals at high risk. Our aim was to evaluate a joint environmental-genetic risk score as a risk stratification tool for early detection of advanced colorectal neoplasm (ACRN). Known environmental risk factors and high-risk genetic loci were summarized into risk scores for ACRN in 1014 eligible participants of a screening study. The performances of single and joint environmental-genetic scores were evaluated with estimates and 95% confidence intervals (CI) of the absolute risk, relative risk and predictive ability using the area under the curve (AUC). Individuals with higher environmental risk scores showed increasing ACRN risk, with 3.1-fold for intermediate risk and 4.8-fold for very high risk, compared to the very low environmental risk group. Similarly, individuals with higher genetic risk scores showed increasing ACRN risk, with 2.2-fold for intermediate risk and 3.5-fold for very high risk, compared to the lowest genetic risk group. Moreover, the joint environmental-genetic score improved the ACRN risk stratification and showed higher predictive values (AUC = 0.64; 95%CI = 0.60-0.67) with substantial difference (p = 0.0002) compared to the single environmental score (0.58; 0.55-0.62). The integration of environmental and genetic factors looks promising for improving targeting individuals at high-risk of colorectal neoplasm. Applications in practical screening programs require optimization with additional genetic and other biomarkers involved in colorectal carcinogenesis.

Colorectal cancer susceptibility variants and risk of conventional adenomas and serrated polyps: results from three cohort studies

BACKGROUND

Increasing evidence suggests that conventional adenomas (CAs) and serrated polyps (SPs) represent two distinct groups of precursor lesions for colorectal cancer (CRC). The influence of common genetic variants on risk of CAs and SPs remain largely unknown.

METHODS

Among 27 426 participants within three prospective cohort studies, we created a weighted genetic risk score (GRS) based on 40 CRC-related single nucleotide polymorphisms (SNPs) identified in previous genome-wide association studies; and we examined the association of GRS (per one standard deviation increment) with risk of CAs, SPs and synchronous CAs and SPs, by multivariable logistic regression. We also analysed individual variants in the secondary analysis.

RESULTS

During 18-20 years of follow-up, we documented 2952 CAs, 1585 SPs and 794 synchronous CAs and SPs. Higher GRS was associated with increased risk of CAs [odds ratio (OR) = 1.17, 95% confidence interval (CI): 1.12-1.21] and SPs (OR = 1.09, 95% CI: 1.03-1.14), with a stronger association for CAs than SPs (Pheterogeneity=0.01). An even stronger association was found for patients with synchronous CAs and SPs (OR = 1.32), advanced CAs (OR = 1.22) and multiple CAs (OR = 1.25). Different sets of variants were associated with CAs and SPs, with a Spearman correlation coefficient of 0.02 between the ORs associating the 40 SNPs with the two lesions. After correcting for multiple testing, three variants were associated with CAs (rs3802842, rs6983267 and rs7136702) and two with SPs (rs16892766 and rs4779584).

CONCLUSIONS

Common genetic variants play a potential role in the conventional and serrated pathways of CRC. Different sets of variants are identified for the two pathways, further supporting the aetiological heterogeneity of CRC.

Establishing a valid approach for estimating familial risk of cancer explained by common genetic variants

We critically examined existing approaches for the estimation of the excess familial risk of cancer that can be attributed to identified common genetic risk variants and propose an alternative, more straightforward approach for calculating this proportion using well-established epidemiological methodology. We applied the underlying equations of the traditional approaches and the new epidemiological approach for colorectal cancer (CRC) in a large population-based case-control study in Germany with 4,447 cases and 3,480 controls, who were recruited from 2003 to 2016 and for whom interview, medical and genomic data were available. Having a family history of CRC (FH) was associated with a 1.77-fold risk increase in our study population (95% CI 1.52-2.07). Traditional approaches yielded estimates of the FH-associated risk explained by 97 common genetics variants from 9.6% to 23.1%, depending on various assumptions. Our alternative approach resulted in smaller and more consistent estimates of this proportion, ranging from 5.4% to 14.3%. Commonly employed methods may lead to strongly divergent and possibly exaggerated estimates of excess familial risk of cancer explained by associated known common genetic variants. Our results suggest that familial risk and risk associated with known common genetic variants might reflect two complementary major sources of risk.

Single nucleotide polymorphisms associated with susceptibility for development of colorectal cancer: Case-control study in a Basque population

Given the significant population diversity in genetic variation, we aimed to investigate whether single nucleotide polymorphisms (SNPs) previously identified in studies of colorectal cancer (CRC) susceptibility were also relevant to the population of the Basque Country (North of Spain). We genotyped 230 CRC cases and 230 healthy controls for 48 previously reported CRC-susceptibility SNPs. Only the rs6687758 in DUPS10 exhibited a statistically significant association with CRC risk based on the crude analysis. The rs6687758 AG genotype conferred about 2.13-fold increased risk for CRC compared to the AA genotype. Moreover, we found significant associations in cases between smoking status, physical activity, and the rs6687758 SNP. The results of a Genetic Risk Score (GRS) showed that the risk alleles were more frequent in cases than controls and the score was associated with CRC in crude analysis. In conclusion, we have confirmed a CRC susceptibility locus and the existence of associations between modifiable factors and the rs6687758 SNP; moreover, the GRS was associated with CRC. However, further experimental validations are needed to establish the role of this SNP, the function of the gene identified, as well as the contribution of the interaction between environmental factors and this locusto the risk of CRC.

Combined effect of modifiable and non-modifiable risk factors for colorectal cancer risk in a pooled analysis of 11 population-based studies

OBJECTIVE

'Environmental' factors associated with colorectal cancer (CRC) risk include modifiable and non-modifiable variables. Whether those with different non-modifiable baseline risks will benefit similarly from reducing their modifiable CRC risks remains unclear.

DESIGN

Using 7945 cases and 8893 controls from 11 population-based studies, we combined 17 risk factors to characterise the overall environmental predisposition to CRC (environmental risk score (E-score)). We estimated the absolute risks (ARs) of CRC of 10 and 30 years across E-score using incidence-rate data from the Surveillance, Epidemiology, and End Results programme. We then combined the modifiable risk factors and estimated ARs across the modifiable risk score, stratified by non-modifiable risk profile based on genetic predisposition, family history and height.

RESULTS

Higher E-score was associated with increased CRC risk (ORquartile, 1.33; 95% CI 1.30 to 1.37). Across E-scores, 30-year ARs of CRC increased from 2.5% in the lowest quartile (Q1) to 5.9% in the highest (Q4) quartile for men, and from 2.1% to 4.5% for women. The modifiable risk score had a stronger association in those with high non-modifiable risk (relative excess risk due to interaction=1.2, 95% CI 0.5 to 1.9). For those in Q4 of non-modifiable risk, a decrease in modifiable risk reduced 30-year ARs from 8.9% to 3.4% for men and from 6.0% to 3.2% for women, a level lower or comparable to the average population risk.

CONCLUSIONS

Changes in modifiable risk factors may result in a substantial decline in CRC risk in both sexes. Those with high inherited risk may reap greater benefit from lifestyle modifications. Our results suggested comprehensive evaluation of environmental factors may facilitate CRC risk stratification.

Ability of known susceptibility SNPs to predict colorectal cancer risk for persons with and without a family history

Before SNP-based risk can be incorporated in colorectal cancer (CRC) screening, the ability of these SNPs to estimate CRC risk for persons with and without a family history of CRC, and the screening implications need to be determined. We estimated the association with CRC of a 45 SNP-based risk using 1181 cases and 999 controls, and its correlation with CRC risk predicted from detailed family history. We estimated the predicted change in the distribution across predefined risk categories, and implications for recommended screening commencement age, from adding SNP-based risk to family history. The inter-quintile risk ratio for colorectal cancer risk of the SNP-based risk was 3.28 (95% CI 2.54-4.22). SNP-based and family history-based risks were not correlated (r = 0.02). For persons with no first-degree relatives with CRC, screening could commence 4 years earlier for women (5 years for men) in the highest quintile of SNP-based risk. For persons with two first-degree relatives with CRC, screening could commence 16 years earlier for men and women in the highest quintile, and 7 years earlier for the lowest quintile. This 45 SNP panel in conjunction with family history, can identify people who could benefit from earlier screening. Risk reclassification by 45 SNPs could inform targeted screening for CRC prevention, particularly in clinical genetics settings when mutations in high-risk genes cannot be identified. Yet to be determined is cost-effectiveness, resources requirements, community, patient and clinician acceptance, and feasibility with potentially ethical, legal and insurance implications.

Exome-wide search and functional annotation of genes associated in patients with severe tick-borne encephalitis in a Russian population

Background

Tick-borne encephalitis (TBE) is a viral infectious disease caused by tick-borne encephalitis virus (TBEV). TBEV infection is responsible for a variety of clinical manifestations ranging from mild fever to severe neurological illness. Genetic factors involved in the host response to TBEV that may potentially play a role in the severity of the disease are still poorly understood. In this study, using whole-exome sequencing, we aimed to identify genetic variants and genes associated with severe forms of TBE as well as biological pathways through which the identified variants may influence the severity of the disease.

Results

Whole-exome sequencing data analysis was performed on 22 Russian patients with severe forms of TBE and 17 Russian individuals from the control group. We identified 2407 candidate genes harboring rare, potentially pathogenic variants in exomes of patients with TBE and not containing any rare, potentially pathogenic variants in exomes of individuals from the control group. According to DAVID tool, this set of 2407 genes was enriched with genes involved in extracellular matrix proteoglycans pathway and genes encoding proteins located at the cell periphery. A total of 154 genes/proteins from these functional groups have been shown to be involved in protein-protein interactions (PPIs) with the known candidate genes/proteins extracted from TBEVHostDB database. By ranking these genes according to the number of rare harmful minor alleles, we identified two genes (MSR1 and LMO7), harboring five minor alleles, and three genes (FLNA, PALLD, PKD1) harboring four minor alleles.

When considering genes harboring genetic variants associated with severe forms of TBE at the suggestive P-value < 0.01, 46 genes containing harmful variants were identified. Out of these 46 genes, eight (MAP4, WDFY4, ACTRT2, KLHL25, MAP2K3, MBD1, OR10J1, and OR2T34) were additionally found among genes containing rare pathogenic variants identified in patients with TBE; and five genes (WDFY4,ALK, MAP4, BNIPL, EPPK1) were found to encode proteins that are involved in PPIs with proteins encoded by genes from TBEVHostDB. Three genes out of five (MAP4, EPPK1, ALK) were found to encode proteins located at cell periphery.

Conclusions

Whole-exome sequencing followed by systems biology approach enabled to identify eight candidate genes (MAP4, WDFY4, ACTRT2, KLHL25, MAP2K3, MBD1, OR10J1, and OR2T34) that can potentially determine predisposition to severe forms of TBE. Analyses of the genetic risk factors for severe forms of TBE revealed a significant enrichment with genes controlling extracellular matrix proteoglycans pathway as well as genes encoding components of cell periphery.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0503-x) contains supplementary material, which is available to authorized users.

Association analyses identify 31 new risk loci for colorectal cancer susceptibility

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

Systematic evaluation of cancer-specific genetic risk score for 11 types of cancer in The Cancer Genome Atlas and Electronic Medical Records and Genomics cohorts

BACKGROUND

Genetic risk score (GRS) is an odds ratio (OR)-weighted and population-standardized method for measuring cumulative effect of multiple risk-associated single nucleotide polymorphisms (SNPs). We hypothesize that GRS is a valid tool for risk assessment of most common cancers.

METHODS

Utilizing genotype and phenotype data from The Cancer Genome Atlas (TCGA) and Electronic Medical Records and Genomics (eMERGE), we tested 11 cancer-specific GRSs (bladder, breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, prostate, renal, and thyroid cancer) for association with the respective cancer type. Cancer-specific GRSs were calculated, for the first time in these cohorts, based on previously published risk-associated SNPs using the Caucasian subjects in these two cohorts.

RESULTS

Mean cancer-specific GRS in the population controls of eMERGE approximated the expected value of 1.00 (between 0.98 and 1.02) for all 11 types of cancer. Mean cancer-specific GRS was consistently higher in respective cancer patients than controls for all 11 types of cancer (P < 0.05). When subjects were categorized into low-, average-, and high-risk groups based on cancer-specific GRS (<0.5, 0.5-1.5, and >1.5, respectively), significant dose-response associations of higher cancer-specific GRS with higher OR of respective type of cancer were found for nine types of cancer (P_-trend  < 0.05). More than 64% subjects in the population controls of eMERGE can be classified as high risk for at least one type of these cancers.

CONCLUSION

Validity of GRS for predicting cancer risk is demonstrated for most types of cancer. If confirmed in larger studies, cancer-specific GRS may have the potential for developing personalized cancer screening strategy.

GWAS identifies two novel colorectal cancer loci at 16q24.1 and 20q13.12

Colorectal cancer (CRC) is the fourth leading cause of cancer mortality worldwide. Genome-wide association studies (GWAS) identified more than 50 CRC loci. However, most of the previous studies were conducted in European population, and host genetic factors among Japanese population are largely remained to be identified. To identify novel loci in the Japanese population, here, we performed a large-scale GWAS using 6692 cases and 27 178 controls followed by a replication analysis using more than 11 000 case-control samples. We found the significant association of 10 loci (P < 5 × 10-8), including 2 novel loci on 16q24.1 (IRF8-FOXF1, rs847208, P = 3.15 × 10-9 and odds ratio = 1.107 with 95% confidence interval (CI) of 1.071-1.145) and 20q13.12 (TOX2, rs6065668, P = 4.47 × 10-11 and odds ratio = 0.897 with 95% CI of 0.868-0.926). Moreover, 35 previously reported single nucleotide polymorphisms (SNPs) in 24 regions were validated in the Japanese population (P < 0.05) with the same risk allele as in the previous studies. SNP rs6065668 was significantly associated with TOX2 expression in the sigmoid colon. In addition, nucleotide substitutions in the regulatory region of TOX2 were predicted to alter the binding of several transcription factors, including KLF5. Our findings elucidate the important role of genetic variations in the development of CRC in the Japanese population.

Discovery of common and rare genetic risk variants for colorectal cancer

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

Genetic susceptibility in the development of colorectal adenomas according to family history of colorectal cancer

Our study aimed to evaluate the relevance of genetic susceptibility in the development of colorectal adenomas (CRA) and its relationship with the presence of family history of colorectal cancer (CRC). Genomic DNA from 750 cases (first degree relatives of patients with CRC) and 750 controls (subjects with no family history of CRC) was genotyped for 99 single nucleotide polymorphisms (SNPs) previously associated with CRC/CRA risk by GWAS and candidate gene studies by using the MassArray™ (Sequenom) platform. Cases and controls were matched by gender, age and histological lesion. Eight hundred and fifty‐eight patients showed no neoplastic lesions, whereas 288 patients showed low‐risk adenomas, and 354 patients presented high‐risk adenomas. Two SNPs (rs10505477, rs6983267) in the CASC8 gene were associated with a reduced risk of CRA in controls (log‐additive models, OR: 0.67, 95%CI:0.54–0.83, and OR:0.66, 95%CI:0.54–0.84, respectively). Stratified analysis by histological lesion revealed the association of rs10505477 and rs6983267 variants with reduced risk of low‐ and high‐risk adenomas in controls, being this effect stronger in low‐risk adenomas (log‐additive models, OR:0.63, 95%CI:0.47–0.84 and OR:0.64, 95%CI:0.47–0.86, respectively). Moreover, 2 SNPs (rs10795668, rs11255841) in the noncoding LINC00709 gene were significantly associated with a reduced risk of low‐risk adenomas in cases (recessive models, OR:0.22, 95%CI:0.06–0.72, and OR:0.08, 95%CI:0.03–0.61) and controls (dominant models, OR:0.50, 95%CI:0.34–0.75, and OR:0.52, 95%CI:0.35–0.78, respectively). In conclusion, some variants associated with CRC risk (rs10505477, rs6983267, rs10795668 and rs11255841) are also involved in the susceptibility to CRA and specific subtypes. These associations are influenced by the presence of family history of CRC.

What's new?

While numerous candidate gene variants have been associated with colorectal cancer, little is known about the relevance of genetic susceptibility or influence of family history in the development of precancerous colorectal adenomas. In the present study, certain genetic variants previously associated with colorectal cancer risk, including two variants in the CASC8 gene and two in the lnc‐RNA LINC00709 gene, were found to be also involved in susceptibility to colorectal adenomas. The associations were modified by family history of colorectal cancer. The results could have implications for colorectal cancer screening and the identification of individuals at increased risk of colorectal adenoma.

Knockdown GREM1 suppresses cell growth, angiogenesis, and epithelial-mesenchymal transition in colon cancer

Gremlin 1 (GREM1), as a bone morphogenetic protein (BMP) antagonist and vascular endothelial growth factor receptor-2 (VEGFR2) novel agonist, has been confirmed as overexpressed in colorectal cancer (CRC) tissues but its role in carcinogenesis remains unclear. Here we reported that the GREM1 expression in mesenchymal-like colon cancer cells (SW620 and SW480) was significantly higher than that of epithelial-like colon cancer cells (Caco-2, HTC116, and HT29) and normal colon cell. Simultaneously, we analyzed two series of CRC transcriptomes from Gene Expression Omnibus (GEO) databases and found the great majority of primary CRC tissues expressed high level of GREM1 messenger RNA (mRNA) compared with adjacent normal tissues, and that the GREM1 mRNA expression is correlated with low histological grade development and stage 2 to 3 metastatic recurrence in CRC based on a data analysis of 104 different stage CRC tissue from the GEO databases. Functional studies showed that GREM1 silencing by short hairpin RNA (shRNA) significantly inhibited CRC cells proliferation, migration, the formation of vascular endothelial growth factor (VEGF)-induced capillary structure of human umbilical vein endothelial cells (HUVECs), and epithelial-mesenchymal transition in colon cancer cells by repressing phosphorylation levels of BMP downstream signal Smad1, vascular endothelial growth factor (VEGF) downstream signal matrix metallopeptidase 2 (MMP2), and metastasis-related factor C-X-C motif chemokine ligand 12 (CXCL12) expression. In addition, shGREM1 combined with VEGF inhibitor BAW2881 displayed more effective antiangiogenesis to inhibit the tube formation of HUVEC. Hence, these experiments demonstrated that GREM1 is involved in CRC development and procession and provide a new idea for CRC diagnosis, resistance therapy, and prognosis.

Genetic Risk Score Is Associated With Prevalence of Advanced Neoplasms in a Colorectal Cancer Screening Population

BACKGROUND & AIMS

The presence of specific single nucleotide polymorphisms (SNPs) can be used to calculate an individual's risk for colorectal cancer (CRC), called a genetic risk score (GRS). We investigated whether GRS can identify individuals with clinically relevant neoplasms in a screening colonoscopy population.

METHODS

We derived a GRS based on 48 SNPs associated with CRC, identified in a comprehensive literature search. We obtained genetic data from 1043 participants (50-79 years old) in a screening colonoscopy study in Germany, recruited from 2005 through 2013 (294 with advanced neoplasms, 249 with non-advanced adenoma (NAAs), and 500 without neoplasms). Each participant was assigned a GRS by aggregating their risk alleles (0, 1, or 2). Risk of advanced neoplasms and NAA according to GRS was calculated by multiple logistic regression. Risk advancement periods were calculated. We replicated our findings using data from a subset of the Tennessee Colorectal Polyp Study.

RESULTS

An increased GRS was associated with higher prevalence of advanced neoplasms, but not NAAs. Participants in the middle and upper tertiles of GRS had a 2.2-fold and 2.7-fold increase in risk, respectively, of advanced neoplasms compared to those in the lower tertile. Adjusted odds ratios (ORs) were 1.09 (95% confidence interval [CI], 0.76-1.57) for NAA in the middle tertile and 1.05 (95% CI, 0.70-1.55) for NAA in the upper tertile. The ORs were largest for proximal advanced neoplasms for participants in the middle tertile (OR, 3.55; 95% CI 1.85-6.82) and the upper tertile (OR, 3.61; 95% CI 1.84-7.10). The risk advancement period for medium vs low GRS was 13.4 years (95% CI 4.8-22.0) and for high vs low GRS was 17.5 years (95% CI, 7.8-27.3).

CONCLUSIONS

In a genetic analysis of participants in a CRC screening study in Germany, an increased GRS (based on CRC-associated SNPs) was associated with increased prevalence of advanced neoplasms. These findings might be used in defining risk-adapted screening ages.

Determining Risk of Colorectal Cancer and Starting Age of Screening Based on Lifestyle, Environmental, and Genetic Factors

BACKGROUND & AIMS

Guidelines for initiating colorectal cancer (CRC) screening are based on family history but do not consider lifestyle, environmental, or genetic risk factors. We developed models to determine risk of CRC, based on lifestyle and environmental factors and genetic variants, and to identify an optimal age to begin screening.

METHODS

We collected data from 9748 CRC cases and 10,590 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colorectal Transdisciplinary study, from 1992 through 2005. Half of the participants were used to develop the risk determination model and the other half were used to evaluate the discriminatory accuracy (validation set). Models of CRC risk were created based on family history, 19 lifestyle and environmental factors (E-score), and 63 CRC-associated single-nucleotide polymorphisms identified in genome-wide association studies (G-score). We evaluated the discriminatory accuracy of the models by calculating area under the receiver operating characteristic curve values, adjusting for study, age, and endoscopy history for the validation set. We used the models to project the 10-year absolute risk of CRC for a given risk profile and recommend ages to begin screening in comparison to CRC risk for an average individual at 50 years of age, using external population incidence rates for non-Hispanic whites from the Surveillance, Epidemiology, and End Results program registry.

RESULTS

In our models, E-score and G-score each determined risk of CRC with greater accuracy than family history. A model that combined both scores and family history estimated CRC risk with an area under the receiver operating characteristic curve value of 0.63 (95% confidence interval, 0.62-0.64) for men and 0.62 (95% confidence interval, 0.61-0.63) for women; area under the receiver operating characteristic curve values based on only family history ranged from 0.53 to 0.54 and those based only E-score or G-score ranged from 0.59 to 0.60. Although screening is recommended to begin at age 50 years for individuals with no family history of CRC, starting ages calculated based on combined E-score and G-score differed by 12 years for men and 14 for women, for individuals with the highest vs the lowest 10% of risk.

CONCLUSIONS

We used data from 2 large international consortia to develop CRC risk calculation models that included genetic and environmental factors along with family history. These determine risk of CRC and starting ages for screening with greater accuracy than the family history only model, which is based on the current screening guideline. These scoring systems might serve as a first step toward developing individualized CRC prevention strategies.

Telomere length and genetics are independent colorectal tumour risk factors in an evaluation of biomarkers in normal bowel

BACKGROUND

Colorectal cancer (CRC) screening might be improved by using a measure of prior risk to modulate screening intensity or the faecal immunochemical test threshold. Intermediate molecular biomarkers could aid risk prediction by capturing both known and unknown risk factors.

METHODS

We sampled normal bowel mucosa from the proximal colon, distal colon and rectum of 317 individuals undergoing colonoscopy. We defined cases as having a personal history of colorectal polyp(s)/cancer, and controls as having no history of colorectal neoplasia. Molecular analyses were performed for: telomere length (TL); global methylation; and the expression of genes in molecular pathways associated with colorectal tumourigenesis. We also calculated a polygenic risk score (PRS) based on CRC susceptibility polymorphisms.

RESULTS

Bowel TL was significantly longer in cases than controls, but was not associated with blood TL. PRS was significantly and independently higher in cases. Hypermethylation showed a suggestive association with case:control status. No gene or pathway was differentially expressed between cases and controls. Gene expression often varied considerably between bowel locations.

CONCLUSIONS

PRS and bowel TL (but not blood TL) may be clinically-useful predictors of CRC risk. Sample collection to assess these biomarkers is feasible in clinical practice, especially where population screening uses flexible sigmoidoscopy or colonoscopy.

Integrative expression quantitative trait locus-based analysis of colorectal cancer identified a functional polymorphism regulating SLC22A5 expression

Multiple single nucleotide polymorphisms (SNPs) have been found to be highly correlated with colorectal cancer (CRC) risk. However, the variants identified thus far only explain a small part of the cases, suggesting the existence of many uncharacterised genetic determinants. In this study, using the multilevel 'omics' data provided in The Cancer Genome Atlas, we systematically performed expression quantitative trait locus (eQTL) analysis for CRC and identified nine SNPs with significant effects on mRNA expression (correlation |r| > 0.3 and FDR < 0.01). Then we conducted a two-stage case-control study consisting of 1528 cases and 1528 controls to examine the associations between candidate SNPs and CRC risk. We found that rs27437 in SLC22A5 was significantly correlated with CRC risk in both stages and the combined study (additive model, OR = 1.31, 95%CI = 1.17-1.47, P = 1.97 × 10^-6). eQTL analysis showed that rs27437 GG and GA genotypes were associated with lower expression of SLC22A5 compared with the AA genotype. Dual-luciferase reporter assays confirmed that the G risk allele could decrease the expression of luciferase. SLC22A5 was significantly decreased in CRC tumour tissues compared with adjacent normal tissues, indicating that SLC22A5 may play important roles in CRC, and rs27437 in SLC22A5 might serve as a novel biomarker for early detection and prevention of CRC.

Strongly enhanced colorectal cancer risk stratification by combining family history and genetic risk score

Background and aim

Family history (FH) and genetic risk scores (GRSs) are increasingly used for risk stratification for colorectal cancer (CRC) screening. However, they were mostly considered alternatively rather than jointly. The aim of this study was to assess the potential of individual and joint risk stratification for CRC by FH and GRS.

Patients and methods

A GRS was built based on the number of risk alleles in 53 previously identified single-nucleotide polymorphisms among 2,363 patients with a first diagnosis of CRC and 2,198 controls in DACHS [colorectal cancer: chances for prevention through screening], a population-based case-control study in Germany. Associations between GRS and FH with CRC risk were quantified by multiple logistic regression.

Results

A total of 316 cases (13.4%) and 214 controls (9.7%) had a first-degree relative (FDR) with CRC (adjusted odds ratio [aOR] 1.86, 95% CI 1.52–2.29). A GRS in the highest decile was associated with a 3.0-fold increased risk of CRC (aOR 3.00, 95% CI 2.24–4.02) compared with the lowest decile. This association was tentatively more pronounced in older age groups. FH and GRS were essentially unrelated, and their joint consideration provided more accurate risk stratification than risk stratification based on each of the variables individually. For example, risk was 6.1-fold increased in the presence of both FH in a FDR and a GRS in the highest decile (aOR 6.14, 95% CI 3.47–10.84) compared to persons without FH and a GRS in the lowest decile.

Conclusion

Both FH and the so far identified genetic variants carry essentially independent risk information and in combination provide great potential for CRC risk stratification.

Identification of a functional variant for colorectal cancer risk mapping to chromosome 5q31.1

Genome-wide association studies (GWASs) have established chromosome 5q31.1 as a risk locus for colorectal cancer (CRC). We previously identified a potentially regulatory single nucleotide polymorphism (SNP) rs17716310 within 5q31.1. Now, we extended our study with another independent Chinese population, functional assays and analyses of TCGA (The Cancer Genome Atlas) data. Significant associations between rs17716310 and CRC risk were found in Present Study including 1075 CRC cases and 1999 controls (additive model: OR = 1.149, 95% CI = 1.027–1.286, P = 0.016), and in Combined Study including 1766 cases and 2708 controls (additive model: OR = 1.145, 95% CI = 1.045–1.254, P = 0.004). Dual luciferase reporter gene assays indicated that the variant C allele obviously increased transcriptional activity. Using TCGA datasets, we indicated rs17716310 as a cis expression quantitative trait locus (eQTL) for the gene SMAD5, whose expression was significantly higher in CRC tissues. These findings suggested that the functional polymorphism rs17716310 A > C might be a genetic modifier for CRC, promoting the expression of SMAD5 that belonged to the transforming growth factor beta (TGF-β) signaling pathway.

An evolutionary perspective on field cancerization

Tumorigenesis begins long before the growth of a clinically detectable lesion and, indeed, even before any of the usual morphological correlates of pre-malignancy are recognizable. Field cancerization, which is the replacement of the normal cell population by a cancer-primed cell population that may show no morphological change, is now recognized to underlie the development of many types of cancer, including the common carcinomas of the lung, colon, skin, prostate and bladder. Field cancerization is the consequence of the evolution of somatic cells in the body that results in cells that carry some but not all phenotypes required for malignancy. Here, we review the evidence of field cancerization across organs and examine the biological mechanisms that drive the evolutionary process that results in field creation. We discuss the clinical implications, principally, how measurements of the cancerized field could improve cancer risk prediction in patients with pre-malignant disease.

Effects of interactions between common genetic variants and smoking on colorectal cancer

Background

Although genome-wide association studies (GWAS) have identified variants in approximately 40 susceptibility loci for colorectal cancer (CRC), there are few studies on the interactions between identified single-nucleotide polymorphisms (SNPs) and lifestyle risk factors. We evaluated whether smoking could modify associations between these genetic variants and CRC risk.

Methods

A total of 703 CRC patients and 1406 healthy controls were included in this case-control study from the National Cancer Center in Korea. Thirty CRC susceptibility SNPs identified in previous GWAS were genotyped. A logistic regression model was used to examine associations between the SNPs and smoking behaviors by sex. The interaction was estimated by including an additional interaction term in the model.

Results

In men, an increased CRC risk was observed for longer durations (OR_{>28 vs. ≤28years} = 1.49 (95% CI = 1.11–1.98)), greater quantities (OR_{≥20 vs. <20cigarettes/day} = 2.12 (1.61–2.79)), and longer pack-years of smoking (OR_{≥21 vs. <21pack-years} = 1.78 (1.35–2.35)). In women, longer pack-years of smoking significantly increased CRC risk (OR_{≥5 vs. <5pack-years} = 6.11 (1.10–34.00)). Moreover, there were significant interactions between smoking status and the polymorphisms rs1957636 at 14q22.3 (P_interaction = 5.5 × 10⁻⁴) and rs4813802 at 20p12.3 (P_interaction = 0.04) in men. Interactions between smoking status and the rs6687758 at 1q41 (P_interaction = 0.03), duration and the rs174537 at 11q12.2 (P_interaction = 0.05), and pack-years and the rs4813802 (P_interaction = 0.04) were also found in women.

Conclusions

Associations between susceptibility SNPs and CRC risk may be modified by smoking behaviors, supporting the existence of gene-smoking interactions.

Electronic supplementary material

The online version of this article (10.1186/s12885-017-3886-0) contains supplementary material, which is available to authorized users.