Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Systematic investigation of genetically determined plasma and urinary metabolites to discover potential interventional targets for colorectal cancer.

BACKGROUND: We aimed to identify plasma and urinary metabolites related to colorectal cancer (CRC) risk and elucidate their mediator role in the associations between modifiable risk factors and CRC. METHODS: Metabolite quantitative trait loci were derived from two published metabolomics genome-wide association studies (GWASs), and summary-level data were extracted for 651 plasma metabolites and 208 urinary metabolites. Genetic associations with CRC were obtained from a large-scale GWAS meta-analysis (100,204 cases; 154,587 controls) and the FinnGen cohort (4,957 cases; 304,197 controls). Mendelian randomization (MR) and colocalization analyses were performed to evaluate the causal roles of metabolites in CRC. Druggability evaluation was employed to prioritize potential therapeutic targets. Multivariable MR and mediation estimation were conducted to elucidate the mediating effects of metabolites on the associations between modifiable risk factors and CRC. RESULTS: The study identified 30 plasma metabolites and four urinary metabolites for CRC. Plasma sphingomyelin and urinary lactose, which were positively associated with CRC risk, could be modulated by drug interventions (ie, Olipudase alfa, Tilactase). Thirteen modifiable risk factors were associated with nine metabolites and eight of these modifiable risk factors were associated with CRC risk. These nine metabolites mediated the effect of modifiable risk factors (Actinobacteria, BMI, waist-hip ratio, fasting insulin, smoking initiation) on CRC. CONCLUSION: This study identified key metabolite biomarkers associated with CRC and elucidated their mediator roles in the associations between modifiable risk factors and CRC. These findings provide new insights into the etiology and potential therapeutic targets for CRC and the etiological pathways of modifiable environmental factors with CRC.

Circulating 25-hydroxyvitamin D and survival outcomes of colorectal cancer: evidence from population-based prospective cohorts and Mendelian randomisation.

BACKGROUND: To investigate the association between circulating 25-hydroxyvitamin D (25-OHD) and colorectal cancer (CRC) survival outcomes. METHODS: We conducted analyses among the Study of Colorectal Cancer in Scotland (SOCCS) and the UK Biobank (UKBB). Both cancer-specific survival (CSS) and overall survival (OS) outcomes were examined. The 25-OHD levels were categorised into three groups, and multi-variable Cox-proportional hazard models were applied to estimate hazard ratios (HRs). We performed individual-level Mendelian randomisation (MR) through the generated polygenic risk scores (PRS) of 25-OHD and summary-level MR using the inverse-variance weighted (IVW) method. RESULTS: We observed significantly poorer CSS (HR = 0.65,95%CI = 0.55-0.76,P = 1.03 × 10-7) and OS (HR = 0.66,95%CI = 0.58-0.75,P = 8.15 × 10-11) in patients with the lowest compared to those with the highest 25-OHD after adjusting for covariates. These associations remained across patients with varied tumour sites and stages. However, we found no significant association between 25-OHD PRS and either CSS (HR = 0.98,95%CI = 0.80-1.19,P = 0.83) or OS (HR = 1.07,95%CI = 0.91-1.25,P = 0.42). Furthermore, we found no evidence for causal effects by conducting summary-level MR analysis for either CSS (IVW:HR = 1.04,95%CI = 0.85-1.28,P = 0.70) or OS (IVW:HR = 1.10,95%CI = 0.93-1.31,P = 0.25). CONCLUSION: This study supports the observed association between lower circulating 25-OHD and poorer survival outcomes for CRC patients. Whilst the genotype-specific association between better outcomes and higher 25-OHD is intriguing, we found no support for causality using MR approaches.

Gene-specific ACMG/AMP classification criteria for germline APC variants: Recommendations from the ClinGen InSiGHT Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel.

PURPOSE: The Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel (VCEP) was established by the International Society for Gastrointestinal Hereditary Tumours and the Clinical Genome Resource, who set out to develop recommendations for the interpretation of germline APC variants underlying Familial Adenomatous Polyposis, the most frequent hereditary polyposis syndrome. METHODS: Through a rigorous process of database analysis, literature review, and expert elicitation, the APC VCEP derived gene-specific modifications to the ACMG/AMP (American College of Medical Genetics and Genomics and Association for Molecular Pathology) variant classification guidelines and validated such criteria through the pilot classification of 58 variants. RESULTS: The APC-specific criteria represented gene- and disease-informed specifications, including a quantitative approach to allele frequency thresholds, a stepwise decision tool for truncating variants, and semiquantitative evaluations of experimental and clinical data. Using the APC-specific criteria, 47% (27/58) of pilot variants were reclassified including 14 previous variants of uncertain significance (VUS). CONCLUSION: The APC-specific ACMG/AMP criteria preserved the classification of well-characterized variants on ClinVar while substantially reducing the number of VUS by 56% (14/25). Moving forward, the APC VCEP will continue to interpret prioritized lists of VUS, the results of which will represent the most authoritative variant classification for widespread clinical use.

Identification of novel protein biomarkers and drug targets for colorectal cancer by integrating human plasma proteome with genome.

BACKGROUND: The proteome is a major source of therapeutic targets. We conducted a proteome-wide Mendelian randomization (MR) study to identify candidate protein markers and therapeutic targets for colorectal cancer (CRC). METHODS: Protein quantitative trait loci (pQTLs) were derived from seven published genome-wide association studies (GWASs) on plasma proteome, and summary-level data were extracted for 4853 circulating protein markers. Genetic associations with CRC were obtained from a large-scale GWAS meta-analysis (16,871 cases and 26,328 controls), the FinnGen cohort (4957 cases and 304,197 controls), and the UK Biobank (9276 cases and 477,069 controls). Colocalization and summary-data-based MR (SMR) analyses were performed sequentially to verify the causal role of candidate proteins. Single cell-type expression analysis, protein-protein interaction (PPI), and druggability evaluation were further conducted to detect the specific cell type with enrichment expression and prioritize potential therapeutic targets. RESULTS: Collectively, genetically predicted levels of 13 proteins were associated with CRC risk. Elevated levels of two proteins (GREM1, CHRDL2) and decreased levels of 11 proteins were associated with an increased risk of CRC, among which four (GREM1, CLSTN3, CSF2RA, CD86) were prioritized with the most convincing evidence. These protein-coding genes are mainly expressed in tissue stem cells, epithelial cells, and monocytes in colon tumor tissue. Two interactive pairs of proteins (GREM1 and CHRDL2; MMP2 and TIMP2) were identified to be involved in osteoclast differentiation and tumorigenesis pathways; four proteins (POLR2F, CSF2RA, CD86, MMP2) have been targeted for drug development on autoimmune diseases and other cancers, with the potentials of being repurposed as therapeutic targets for CRC. CONCLUSIONS: This study identified several protein biomarkers to be associated with CRC risk and provided new insights into the etiology and promising targets for the development of screening biomarkers and therapeutic drugs for CRC.

Dissecting the pathogenic effects of smoking and its hallmarks in blood DNA methylation on colorectal cancer risk.

BACKGROUND: Tobacco smoking is suggested as a risk factor for colorectal cancer (CRC), but the complex relationship and the potential pathway are not fully understood. METHODS: We performed two-sample Mendelian randomisation (MR) analyses with genetic instruments for smoking behaviours and related DNA methylation in blood and summary-level GWAS data of colorectal cancer to disentangle the relationship. Colocalization analyses and prospective gene-environment interaction analyses were also conducted as replication. RESULTS: Convincing evidence was identified for the pathogenic effect of smoking initiation on CRC risk and suggestive evidence was observed for the protective effect of smoking cessation in the univariable MR analyses. Multivariable MR analysis revealed that these associations were independent of other smoking phenotypes and alcohol drinking. Genetically predicted methylation at CpG site cg17823346 [ZMIZ1] were identified to decrease CRC risk; while genetically predicted methylation at cg02149899 would increase CRC risk. Colocalization and gene-environment interaction analyses added further evidence to the relationship between epigenetic modification at cg17823346 [ZMIZ1] as well as cg02149899 and CRC risk. DISCUSSION: Our study confirms the significant association between tobacco smoking, DNA methylation and CRC risk and yields a novel insight into the pathogenic effect of tobacco smoking on CRC risk.

Development and Validation of Risk Prediction Models for Colorectal Cancer in Patients with Symptoms.

We aimed to develop and validate prediction models incorporating demographics, clinical features, and a weighted genetic risk score (wGRS) for individual prediction of colorectal cancer (CRC) risk in patients with gastroenterological symptoms. Prediction models were developed with internal validation [CRC Cases: n = 1686/Controls: n = 963]. Candidate predictors included age, sex, BMI, wGRS, family history, and symptoms (changes in bowel habits, rectal bleeding, weight loss, anaemia, abdominal pain). The baseline model included all the non-genetic predictors. Models A (baseline model + wGRS) and B (baseline model) were developed based on LASSO regression to select predictors. Models C (baseline model + wGRS) and D (baseline model) were built using all variables. Models' calibration and discrimination were evaluated through the Hosmer-Lemeshow test (calibration curves were plotted) and C-statistics (corrected based on 1000 bootstrapping). The models' prediction performance was: model A (corrected C-statistic = 0.765); model B (corrected C-statistic = 0.753); model C (corrected C-statistic = 0.764); and model D (corrected C-statistic = 0.752). Models A and C, that integrated wGRS with demographic and clinical predictors, had a statistically significant improved prediction performance. Our findings suggest that future application of genetic predictors holds significant promise, which could enhance CRC risk prediction. Therefore, further investigation through model external validation and clinical impact is merited.

Exploring the cross-cancer effect of smoking and its fingerprints in blood DNA methylation on multiple cancers: A Mendelian randomization study.

Aberrant smoking-related DNA methylation has been widely investigated as a carcinogenesis mechanism, but whether the cross-cancer epigenetic pathways exist remains unclear. We conducted two-sample Mendelian randomization (MR) analyses respectively on smoking behaviors (age of smoking initiation, smoking initiation, smoking cessation, and lifetime smoking index [LSI]) and smoking-related DNA methylation to investigate their effect on 15 site-specific cancers, based on a genome-wide association study (GWAS) of 1.2 million European individuals and an epigenome-WAS (EWAS) of 5907 blood samples of Europeans for smoking and 15 GWASs of European ancestry for multiple site-specific cancers. Significantly identified CpG sites were further used for colocalization analysis, and those with cross-cancer effect were validated by overlapping with tissue-specific eQTLs. In the genomic MR, smoking measurements of smoking initiation, smoking cessation and LSI were suggested to be casually associated with risk of seven types of site-specific cancers, among which cancers at lung, cervix and colorectum were provided with strong evidence. In the epigenetic MR, methylation at 75 CpG sites were reported to be significantly associated with increased risks of multiple cancers. Eight out of 75 CpG sites were observed with cross-cancer effect, among which cg06639488 (EFNA1), cg12101586 (CYP1A1) and cg14142171 (HLA-L) were validated by eQTLs at specific cancer sites, and cg07932199 (ATXN2) had strong evidence to be associated with cancers of lung (coefficient, 0.65, 95% confidence interval [CI], 0.31-1.00), colorectum (0.90 [0.61, 1.18]), breast (0.31 [0.20, 0.43]) and endometrium (0.98 [0.68, 1.27]). These findings highlight the potential practices targeting DNA methylation-involved cross-cancer pathways.

Exploring the Complex Relationship between Gut Microbiota and Risk of Colorectal Neoplasia Using Bidirectional Mendelian Randomization Analysis.

BACKGROUND: Human gut microbiome has complex relationships with the host, contributing to metabolism, immunity, and carcinogenesis. METHODS: Summary-level data for gut microbiota and metabolites were obtained from MiBioGen, FINRISK and human metabolome consortia. Summary-level data for colorectal cancer were derived from a genome-wide association study meta-analysis. In forward Mendelian randomization (MR), we employed genetic instrumental variables (IV) for 24 gut microbiota taxa and six bacterial metabolites to examine their causal relationship with colorectal cancer. We also used a lenient threshold for nine apriori gut microbiota taxa as secondary analyses. In reverse MR, we explored association between genetic liability to colorectal neoplasia and abundance of microbiota studied above using 95, 19, and 7 IVs for colorectal cancer, adenoma, and polyps, respectively. RESULTS: Forward MR did not find evidence indicating causal relationship between any of the gut microbiota taxa or six bacterial metabolites tested and colorectal cancer risk. However, reverse MR supported genetic liability to colorectal adenomas was causally related with increased abundance of two taxa: Gammaproteobacteria (ß = 0.027, which represents a 0.027 increase in log-transformed relative abundance values of Gammaproteobacteria for per one-unit increase in log OR of adenoma risk; P = 7.06×10-8), Enterobacteriaceae (ß = 0.023, P = 1.29×10-5). CONCLUSIONS: We find genetic liability to colorectal neoplasia may be associated with abundance of certain microbiota taxa. It is more likely that subset of colorectal cancer genetic liability variants changes gut biology by influencing both gut microbiota and colorectal cancer risk. IMPACT: This study highlights the need of future complementary studies to explore causal mechanisms linking both host genetic variation with gut microbiome and colorectal cancer susceptibility.

Cross-cancer pleiotropic analysis identifies three novel genetic risk loci for colorectal cancer.

BACKGROUND: To understand the shared genetic basis between colorectal cancer (CRC) and other cancers and identify potential pleiotropic loci for compensating the missing genetic heritability of CRC. METHODS: We conducted a systematic genome-wide pleiotropy scan to appraise associations between cancer-related genetic variants and CRC risk among European populations. Single nucleotide polymorphism (SNP)-set analysis was performed using data from the UK Biobank and the Study of Colorectal Cancer in Scotland (10 039 CRC cases and 30 277 controls) to evaluate the overlapped genetic regions for susceptibility of CRC and other cancers. The variant-level pleiotropic associations between CRC and other cancers were examined by CRC genome-wide association study meta-analysis and the pleiotropic analysis under composite null hypothesis (PLACO) pleiotropy test. Gene-based, co-expression and pathway enrichment analyses were performed to explore potential shared biological pathways. The interaction between novel genetic variants and common environmental factors was further examined for their effects on CRC. RESULTS: Genome-wide pleiotropic analysis identified three novel SNPs (rs2230469, rs9277378 and rs143190905) and three mapped genes (PIP4K2A, HLA-DPB1 and RTEL1) to be associated with CRC. These genetic variants were significant expressions quantitative trait loci in colon tissue, influencing the expression of their mapped genes. Significant interactions of PIP4K2A and HLA-DPB1 with environmental factors, including smoking and alcohol drinking, were observed. All mapped genes and their co-expressed genes were significantly enriched in pathways involved in carcinogenesis. CONCLUSION: Our findings provide an important insight into the shared genetic basis between CRC and other cancers. We revealed several novel CRC susceptibility loci to help understand the genetic architecture of CRC.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

Contrast enhanced magneto-motive ultrasound in lymph nodes - modelling and pre-clinical imaging using magnetic microbubbles.

Despite advances in MRI, the detection and characterisation of lymph nodes in rectal cancer remains complex, especially when assessing the response to neo-adjuvant treatment. An alternative approach is functional imaging, previously shown to aid characterization of cancer tissues. We report proof-of-concept of the novel technique Contrast-Enhanced Magneto-Motive Ultrasound (CE-MMUS) to recover information relating to local perfusion and lymphatic drainage, and interrogate tissue mechanical properties through magnetically induced tissue deformations. The feasibility of the proposed application was explored using a combination of pre-clinical ultrasound imaging and finite element analysis. First, contrast enhanced ultrasound imaging on one wild type mouse recorded lymphatic drainage of magnetic microbubbles after bolus injection. Second, preliminary CE-MMUS data were acquired as a proof of concept. Third, the magneto-mechanical interactions of a magnetic microbubble with an elastic solid were simulated using finite element software. Accumulation of magnetic microbubbles in the inguinal lymph node was verified using contrast enhanced ultrasound, with peak enhancement occurring 3.7 s post-injection. Preliminary CE-MMUS indicates the presence of magnetic contrast agent in the lymph node. The finite element analysis explores how the magnetic force is transferred to motion of the solid, which depends on elasticity and bubble radius, indicating an inverse relation with displacement. Combining magnetic microbubbles with MMUS could harness the advantages of both techniques, to provide perfusion information, robust lymph node delineation and characterisation based on mechanical properties. Clinical Relevance- Robust detection and characterisation of lymph nodes could be aided by visualising lymphatic drainage of magnetic microbubbles using contrast enhanced ultrasound imaging and magneto-motion, which is dependent on tissue mechanical properties.

Transcriptional dynamics of colorectal cancer risk associated variation at 11q23.1 correlate with tuft cell abundance and marker expression in silico.

Colorectal cancer (CRC) is characterised by heritable risk that is not well understood. Heritable, genetic variation at 11q23.1 is associated with increased colorectal cancer (CRC) risk, demonstrating eQTL effects on 3 cis- and 23 trans-eQTL targets. We sought to determine the relationship between 11q23.1 cis- and trans-eQTL target expression and test for potential cell-specificity. scRNAseq from 32,361 healthy colonic epithelial cells was aggregated and subject to weighted gene co-expression network analysis (WGCNA). One module (blue) included 19 trans-eQTL targets and was correlated with POU2AF2 expression only. Following unsupervised clustering of single cells, the expression of 19 trans-eQTL targets was greatest and most variable in cluster number 11, which transcriptionally resembled tuft cells. 14 trans-eQTL targets were found to demarcate this cluster, 11 of which were corroborated in a second dataset. Intra-cluster WGCNA and module preservation analysis then identified twelve 11q23.1 trans-eQTL targets to comprise a network that was specific to cluster 11. Finally, linear modelling and differential abundance testing showed 11q23.1 trans-eQTL target expression was predictive of cluster 11 abundance. Our findings suggest 11q23.1 trans-eQTL targets comprise a POU2AF2-related network that is likely tuft cell-specific and reduced expression of these genes correlates with reduced tuft cell abundance in silico.

Modelling of magnetic microbubbles to evaluate contrast enhanced magnetomotive ultrasound in lymph nodes - a pre-clinical study.

OBJECTIVES: Despite advances in MRI the detection and characterisation of lymph nodes in rectal cancer remains complex, especially when assessing the response to neoadjuvant treatment. An alternative approach is functional imaging, previously shown to aid characterisation of cancer tissues. We report proof of concept of the novel technique Contrast-Enhanced Magneto-Motive Ultrasound (CE-MMUS) to recover information relating to local perfusion and lymphatic drainage, and interrogate tissue mechanical properties through magnetically induced deformations. METHODS: The feasibility of the proposed application was explored using a combination of experimental animal and phantom ultrasound imaging, along with finite element analysis. First, contrast-enhanced ultrasound imaging on one wild type mouse recorded lymphatic drainage of magnetic microbubbles after bolus injection. Second, tissue phantoms were imaged using MMUS to illustrate the force- and elasticity dependence of the magnetomotion. Third, the magnetomechanical interactions of a magnetic microbubble with an elastic solid were simulated using finite element software. RESULTS: Accumulation of magnetic microbubbles in the inguinal lymph node was verified using contrast enhanced ultrasound, with peak enhancement occurring 3.7 s post-injection. The magnetic microbubble gave rise to displacements depending on force, elasticity, and bubble radius, indicating an inverse relation between displacement and the latter two. CONCLUSION: Combining magnetic microbubbles with MMUS could harness the advantages of both techniques, to provide perfusion information, robust lymph node delineation and characterisation based on mechanical properties. ADVANCES IN KNOWLEDGE: (a) Lymphatic drainage of magnetic microbubbles visualised using contrast-enhanced ultrasound imaging and (b) magnetomechanical interactions between such bubbles and surrounding tissue could both contribute to (c) robust detection and characterisation of lymph nodes.

Development of Preclinical Ultrasound Imaging Techniques to Identify and Image Sentinel Lymph Nodes in a Cancerous Animal Model.

Lymph nodes (LNs) are believed to be the first organs targeted by colorectal cancer cells detached from a primary solid tumor because of their role in draining interstitial fluids. Better detection and assessment of these organs have the potential to help clinicians in stratification and designing optimal design of oncological treatments for each patient. Whilst highly valuable for the detection of primary tumors, CT and MRI remain limited for the characterization of LNs. B-mode ultrasound (US) and contrast-enhanced ultrasound (CEUS) can improve the detection of LNs and could provide critical complementary information to MRI and CT scans; however, the European Federation of Societies for Ultrasound in Medicine and Biology (EFSUMB) guidelines advise that further evidence is required before US or CEUS can be recommended for clinical use. Moreover, knowledge of the lymphatic system and LNs is relatively limited, especially in preclinical models. In this pilot study, we have created a mouse model of metastatic cancer and utilized 3D high-frequency ultrasound to assess the volume, shape, and absence of hilum, along with CEUS to assess the flow dynamics of tumor-free and tumor-bearing LNs in vivo. The aforementioned parameters were used to create a scoring system to predict the likelihood of a disease-involved LN before establishing post-mortem diagnosis with histopathology. Preliminary results suggest that a sum score of parameters may provide a more accurate diagnosis than the LN size, the single parameter currently used to predict the involvement of an LN in disease.

Bidirectional Mendelian randomisation analysis of the relationship between circulating vitamin D concentration and colorectal cancer risk.

Epidemiological evidence is consistent with a protective effect of vitamin D against colorectal cancer (CRC), but the observed strong associations are open to confounders and potential reverse causation. Previous Mendelian randomisation (MR) studies were limited by poor genetic instruments and inadequate statistical power. Moreover, whether genetically higher CRC risk can influence vitamin D level, namely the reverse causation, still remains unknown. Herein, we report the first bidirectional MR study. We employed 110 newly identified genetic variants as proxies for vitamin D to obtain unconfounded effect estimates on CRC risk in 26 397 CRC cases and 41 481 controls of European ancestry. To test for reserve causation, we estimated effects of 115 CRC-risk variants on vitamin D level among 417 580 participants from the UK Biobank. The causal association was estimated using the random-effect inverse-variance weighted (IVW) method. We found no significant causal effect of vitamin D on CRC risk [IVW estimate odds ratio: 0.97, 95% confidence interval (CI) = 0.88-1.07, P = .565]. Similarly, no significant reverse causal association was identified between genetically increased CRC risk and vitamin D levels (IVW estimate ß: -0.002, 95% CI = -0.008 to 0.004, P = .543). Stratified analysis by tumour sites did not identify significant causal associations in either direction between vitamin D and colon or rectal cancer. Despite the improved statistical power of this study, we found no evidence of causal association of either direction between circulating vitamin D and CRC risk. Significant associations reported by observational studies may be primarily driven by unidentified confounders.

Phenome-wide association study (PheWAS) of colorectal cancer risk SNP effects on health outcomes in UK Biobank.

BACKGROUND: Associations between colorectal cancer (CRC) and other health outcomes have been reported, but these may be subject to biases, or due to limitations of observational studies. METHODS: We set out to determine whether genetic predisposition to CRC is also associated with the risk of other phenotypes. Under the phenome-wide association study (PheWAS) and tree-structured phenotypic model (TreeWAS), we studied 334,385 unrelated White British individuals (excluding CRC patients) from the UK Biobank cohort. We generated a polygenic risk score (PRS) from CRC genome-wide association studies as a measure of CRC risk. We performed sensitivity analyses to test the robustness of the results and searched the Danish Disease Trajectory Browser (DTB) to replicate the observed associations. RESULTS: Eight PheWAS phenotypes and 21 TreeWAS nodes were associated with CRC genetic predisposition by PheWAS and TreeWAS, respectively. The PheWAS detected associations were from neoplasms and digestive system disease group (e.g. benign neoplasm of colon, anal and rectal polyp and diverticular disease). The results from the TreeWAS corroborated the results from the PheWAS. These results were replicated in the observational data within the DTB. CONCLUSIONS: We show that benign colorectal neoplasms share genetic aetiology with CRC using PheWAS and TreeWAS methods. Additionally, CRC genetic predisposition is associated with diverticular disease.

Vitamin D treatment induces in vitro and ex vivo transcriptomic changes indicating anti-tumor effects.

Vitamin D deficiency is associated with risk of several common cancers, including colorectal cancer (CRC). Here we have utilized patient derived epithelial organoids (ex vivo) and CRC cell lines (in vitro) to show that calcitriol (1,25OHD) increased the expression of the CRC tumor suppressor gene, CDH1, at both the transcript and protein level. Whole genome expression analysis demonstrated significant differential expression of a further six genes after 1,25OHD treatment, including genes with established links to carcinogenesis GADD45, EFTUD1 and KIAA1199. Furthermore, gene ontologies relevant to carcinogenesis were enriched by 1,25OHD treatment (e.g., 'regulation of Wnt signaling pathway', 'regulation of cell death'), with common enriched processes across in vitro and ex vivo cultures including 'negative regulation of cell proliferation', 'regulation of cell migration' and 'regulation of cell differentiation'. Our results identify genes and pathways that are modifiable by calcitriol that have links to CRC tumorigenesis. Hence the findings provide potential mechanism to the epidemiological and clinical trial data indicating a causal association between vitamin D and CRC. We suggest there is strong rationale for further well-designed trials of vitamin D supplementation as a novel CRC chemopreventive and chemotherapeutic agent.

Differential genetic influences over colorectal cancer risk and gene expression in large bowel mucosa.

Site-specific variation in colorectal cancer (CRC) incidence, biology and prognosis are poorly understood. We sought to determine whether common genetic variants influencing CRC risk might exhibit topographical differences on CRC risk through regional differences in effects on gene expression in the large bowel mucosa. We conducted a site-specific genetic association study (10 630 cases, 31 331 controls) to identify whether established risk variants exert differential effects on risk of proximal, compared to distal CRC. We collected normal colorectal mucosa and blood from 481 subjects and assessed mucosal gene expression using Illumina HumanHT-12v4 arrays in relation to germline genotype. Expression quantitative trait loci (eQTLs) were explored by anatomical location of sampling. The rs3087967 genotype (chr11q23.1 risk variant) exhibited significant site-specific effects-risk of distal CRC (odds ratio [OR] = 1.20, P = 8.20 × 10-20 ) with negligible effects on proximal CRC risk (OR = 1.05, P = .10). Expression of 1261 genes differed between proximal and distal colonic mucosa (top hit PRAC gene, fold-difference = 10, P = 3.48 × 10-57 ). In eQTL studies, rs3087967 genotype was associated with expression of 8 cis- and 21 trans-genes. Four of these (AKAP14, ADH5P4, ASGR2, RP11-342M1.7) showed differential effects by site, with strongest trans-eQTL signals in proximal colonic mucosa (eg, AKAP14, beta = 0.61, P = 5.02 × 10-5 ) and opposite signals in distal mucosa (AKAP14, beta = -0.17, P = .04). In summary, genetic variation at the chr11q23.1 risk locus imparts greater risk of distal rather than proximal CRC and exhibits site-specific differences in eQTL effects in normal mucosa. Topographical differences in genomic control over gene expression relevant to CRC risk may underlie site-specific variation in CRC. Results may inform individualised CRC screening programmes.