A risk variant for Barrett's esophagus and esophageal adenocarcinoma at chr8p23.1 affects enhancer activity and implicates multiple gene targets.

Nineteen genetic susceptibility loci for esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) have been identified through genome-wide association studies (GWAS). Clinical translation of such discoveries, however, has been hindered by the slow pace of discovery of functional/causal variants and gene targets at these loci. We previously developed a systematic informatics pipeline to prioritize candidate functional variants using functional potential scores, applied the pipeline to select high-scoring BE/EAC risk loci and validated a functional variant at chr19p13.11 (rs10423674). Here, we selected two additional prioritized loci for experimental interrogation: chr3p13/rs1522552 and chr8p23.1/rs55896564. Candidate enhancer regions encompassing these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two regions tested exhibited allele-specific enhancer activity - 8p23.1/rs55896564. CRISPR-mediated deletion of the putative enhancer in EAC cell lines correlated with reduced expression of three candidate gene targets: B lymphocyte kinase (BLK), nei like DNA glycosylase 2 (NEIL2) and cathepsin B (CTSB). Expression quantitative trait locus (eQTL) mapping in normal esophagus and stomach revealed strong associations between the BE/EAC risk allele at rs55896564 (G) and lower expression of CTSB, a protease gene implicated in epithelial wound repair. These results further support the utility of functional potential scores for GWAS variant prioritization, and provide the first experimental evidence of a functional variant and risk enhancer at the 8p23.1 GWAS locus. Identification of CTSB, BLK and NEIL2 as candidate gene targets suggests that altered expression of these genes may underlie the genetic risk association at 8p23.1 with BE/EAC.

Prioritization and functional analysis of GWAS risk loci for Barrett's esophagus and esophageal adenocarcinoma.

Genome-wide association studies (GWAS) have identified ~20 genetic susceptibility loci for esophageal adenocarcinoma (EAC), and its precursor, Barrett's esophagus (BE). Despite such advances, functional/causal variants and gene targets at these loci remain undefined, hindering clinical translation. A key challenge is that most causal variants map to non-coding regulatory regions such as enhancers, and typically, numerous potential candidate variants at GWAS loci require testing. We developed a systematic informatics pipeline for prioritizing candidate functional variants via integrative functional potential scores (FPS) consolidated from multi-omics annotations, and used this pipeline to identify two high-scoring variants for experimental interrogation: chr9q22.32/rs11789015 and chr19p13.11/rs10423674. Minimal candidate enhancer regions spanning these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two variants tested (rs10423674) exhibited allele-specific enhancer activity. CRISPR-mediated deletion of the putative enhancer region in EAC cell lines correlated with reduced expression of two genes-CREB-regulated transcription coactivator 1 (CRTC1) and Cartilage oligomeric matrix protein (COMP); expression of five other genes remained unchanged (CRLF1, KLHL26, TMEM59L, UBA52, RFXANK). Expression quantitative trait locus mapping indicated that rs10423674 genotype correlated with CRTC1 and COMP expression in normal esophagus. This study represents the first experimental effort to bridge GWAS associations to biology in BE/EAC and supports the utility of FPS to guide variant prioritization. Our findings reveal a functional variant and candidate risk enhancer at chr19p13.11 and implicate CRTC1 and COMP as putative gene targets, suggesting that altered expression of these genes may underlie the BE/EAC risk association.

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.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses.

BACKGROUND: Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis. METHODS: We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test). RESULTS: Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09). DISCUSSION: These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Genetic variants associated with circulating C-reactive protein levels and colorectal cancer survival: Sex-specific and lifestyle factors specific associations.

Elevated blood levels of C-reactive protein (CRP) have been linked to colorectal cancer (CRC) survival. We evaluated genetic variants associated with CRP levels and their interactions with sex and lifestyle factors in association with CRC-specific mortality. Our study included 16 142 CRC cases from the International Survival Analysis in Colorectal Cancer Consortium. We identified 618 common single nucleotide polymorphisms (SNPs) associated with CRP levels from the NHGRI-EBI GWAS Catalog. Cox proportional hazards regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for associations between SNPs and CRC-specific mortality adjusting for age, sex, genotyping platform/study and principal components. We investigated their interactions with sex and lifestyle factors using likelihood ratio tests. Of 5472 (33.9%) deaths accrued over up to 10 years of follow-up, 3547 (64.8%) were due to CRC. No variants were associated with CRC-specific mortality after multiple comparison correction. We observed strong evidence of interaction between variant rs1933736 at FRK gene and sex in relation to CRC-specific mortality (corrected Pinteraction  = .0004); women had higher CRC-specific mortality associated with the minor allele (HR = 1.11, 95% CI = 1.04-1.19) whereas an inverse association was observed for men (HR = 0.88, 95% CI = 0.82-0.94). There was no evidence of interactions between CRP-associated SNPs and alcohol, obesity or smoking. Our study observed a significant interaction between sex and a CRP-associated variant in relation to CRC-specific mortality. Future replication of this association and functional annotation of the variant are needed.

Genetically proxied therapeutic inhibition of antihypertensive drug targets and risk of common cancers: A mendelian randomization analysis.

BACKGROUND: Epidemiological studies have reported conflicting findings on the potential adverse effects of long-term antihypertensive medication use on cancer risk. Naturally occurring variation in genes encoding antihypertensive drug targets can be used as proxies for these targets to examine the effect of their long-term therapeutic inhibition on disease outcomes. METHODS AND FINDINGS: We performed a mendelian randomization analysis to examine the association between genetically proxied inhibition of 3 antihypertensive drug targets and risk of 4 common cancers (breast, colorectal, lung, and prostate). Single-nucleotide polymorphisms (SNPs) in ACE, ADRB1, and SLC12A3 associated (P < 5.0 × 10-8) with systolic blood pressure (SBP) in genome-wide association studies (GWAS) were used to proxy inhibition of angiotensin-converting enzyme (ACE), ß-1 adrenergic receptor (ADRB1), and sodium-chloride symporter (NCC), respectively. Summary genetic association estimates for these SNPs were obtained from GWAS consortia for the following cancers: breast (122,977 cases, 105,974 controls), colorectal (58,221 cases, 67,694 controls), lung (29,266 cases, 56,450 controls), and prostate (79,148 cases, 61,106 controls). Replication analyses were performed in the FinnGen consortium (1,573 colorectal cancer cases, 120,006 controls). Cancer GWAS and FinnGen consortia data were restricted to individuals of European ancestry. Inverse-variance weighted random-effects models were used to examine associations between genetically proxied inhibition of these drug targets and risk of cancer. Multivariable mendelian randomization and colocalization analyses were employed to examine robustness of findings to violations of mendelian randomization assumptions. Genetically proxied ACE inhibition equivalent to a 1-mm Hg reduction in SBP was associated with increased odds of colorectal cancer (odds ratio (OR) 1.13, 95% CI 1.06 to 1.22; P = 3.6 × 10-4). This finding was replicated in the FinnGen consortium (OR 1.40, 95% CI 1.02 to 1.92; P = 0.035). There was little evidence of association of genetically proxied ACE inhibition with risk of breast cancer (OR 0.98, 95% CI 0.94 to 1.02, P = 0.35), lung cancer (OR 1.01, 95% CI 0.92 to 1.10; P = 0.93), or prostate cancer (OR 1.06, 95% CI 0.99 to 1.13; P = 0.08). Genetically proxied inhibition of ADRB1 and NCC were not associated with risk of these cancers. The primary limitations of this analysis include the modest statistical power for analyses of drug targets in relation to some less common histological subtypes of cancers examined and the restriction of the majority of analyses to participants of European ancestry. CONCLUSIONS: In this study, we observed that genetically proxied long-term ACE inhibition was associated with an increased risk of colorectal cancer, warranting comprehensive evaluation of the safety profiles of ACE inhibitors in clinical trials with adequate follow-up. There was little evidence to support associations across other drug target-cancer risk analyses, consistent with findings from short-term randomized controlled trials for these medications.

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.

A Functional Variant on 9p21.3 Related to Glioma Risk Affects Enhancer Activity and Modulates Expression of CDKN2B-AS1.

Genome-wide association studies have identified SNPs associated with glioma risk on 9p21.3, but biological mechanisms underlying this association are unknown. We tested the hypothesis that a functional SNP on 9p21.3 affects activity of an enhancer, causing altered expression of nearby genes. We considered all SNPs in linkage disequilibrium with the 9p21.3 sentinel SNP rs634537 that mapped to putative enhancers. An enhancer containing rs1537372 exhibited allele-specific effects on luciferase activity. Deletion of this enhancer in GBM cell lines correlated with decreased expression of CDKN2B-AS1. Expression quantitative trait loci analysis using non-diseased brain samples showed rs1537372 to be a consistently significant eQTL for CDKN2B-AS1. Additionally, our analysis of Hi-C data generated in neural progenitor cells showed that the bait region containing rs1537372 interacted with the CDKN2B-AS1 promoter. These data suggest rs1537372, a SNP at the 9p21.3 risk locus, is a functional variant that modulates expression of CDKN2B-AS1.

A functional variant on 20q13.33 related to glioma risk alters enhancer activity and modulates expression of multiple genes.

Genome-wide association studies (GWAS) have identified single-nucleotide polymorphisms (SNPs) associated with glioma risk on 20q13.33, but the biological mechanisms underlying this association are unknown. We tested the hypothesis that a functional SNP on 20q13.33 impacted the activity of an enhancer, leading to an altered expression of nearby genes. To identify candidate functional SNPs, we identified all SNPs in linkage disequilibrium with the risk-associated SNP rs2297440 that mapped to putative enhancers. Putative enhancers containing candidate functional SNPs were tested for allele-specific effects in luciferase enhancer activity assays against glioblastoma multiforme (GBM) cell lines. An enhancer containing SNP rs3761124 exhibited allele-specific effects on activity. Deletion of this enhancer by CRISPR-Cas9 editing in GBM cell lines correlated with an altered expression of multiple genes, including STMN3, RTEL1, RTEL1-TNFRSF6B, GMEB2, and SRMS. Expression quantitative trait loci (eQTL) analyses using nondiseased brain samples, isocitrate dehydrogenase 1 (IDH1) wild-type glioma, and neurodevelopmental tissues showed STMN3 to be a consistent significant eQTL with rs3761124. RTEL1 and GMEB2 were also significant eQTLs in the context of early CNS development and/or in IDH1 wild-type glioma. We provide evidence that rs3761124 is a functional variant on 20q13.33 related to glioma/GBM risk that modulates the expression of STMN3 and potentially other genes across diverse cellular contexts.

Lack of an association between gallstone disease and bilirubin levels with risk of colorectal cancer: a Mendelian randomisation analysis.

BACKGROUND: Epidemiological studies of the relationship between gallstone disease and circulating levels of bilirubin with risk of developing colorectal cancer (CRC) have been inconsistent. To address possible confounding and reverse causation, we examine the relationship between these potential risk factors and CRC using Mendelian randomisation (MR). METHODS: We used two-sample MR to examine the relationship between genetic liability to gallstone disease and circulating levels of bilirubin with CRC in 26,397 patients and 41,481 controls. We calculated the odds ratio per genetically predicted SD unit increase in log bilirubin levels (ORSD) for CRC and tested for a non-zero causal effect of gallstones on CRC. Sensitivity analysis was applied to identify violations of estimator assumptions. RESULTS: No association between either gallstone disease (P value = 0.60) or circulating levels of bilirubin (ORSD = 1.00, 95% confidence interval (CI) = 0.96-1.03, P value = 0.90) with CRC was shown. CONCLUSIONS: Despite the large scale of this study, we found no evidence for a causal relationship between either circulating levels of bilirubin or gallstone disease with risk of developing CRC. While the magnitude of effect suggested by some observational studies can confidently be excluded, we cannot exclude the possibility of smaller effect sizes and non-linear relationships.

A Mixed-Effects Model for Powerful Association Tests in Integrative Functional Genomics.

Genome-wide association studies (GWASs) have successfully identified thousands of genetic variants for many complex diseases; however, these variants explain only a small fraction of the heritability. Recently, genetic association studies that leverage external transcriptome data have received much attention and shown promise for discovering novel variants. One such approach, PrediXcan, is to use predicted gene expression through genetic regulation. However, there are limitations in this approach. The predicted gene expression may be biased, resulting from regularized regression applied to moderately sample-sized reference studies. Further, some variants can individually influence disease risk through alternative functional mechanisms besides expression. Thus, testing only the association of predicted gene expression as proposed in PrediXcan will potentially lose power. To tackle these challenges, we consider a unified mixed effects model that formulates the association of intermediate phenotypes such as imputed gene expression through fixed effects, while allowing residual effects of individual variants to be random. We consider a set-based score testing framework, MiST (mixed effects score test), and propose two data-driven combination approaches to jointly test for the fixed and random effects. We establish the asymptotic distributions, which enable rapid calculation of p values for genome-wide analyses, and provide p values for fixed and random effects separately to enhance interpretability over GWASs. Extensive simulations demonstrate that our approaches are more powerful than existing ones. We apply our approach to a large-scale GWAS of colorectal cancer and identify two genes, POU5F1B and ATF1, which would have otherwise been missed by PrediXcan, after adjusting for all known loci.

Circulating Levels of Insulin-like Growth Factor 1 and Insulin-like Growth Factor Binding Protein 3 Associate With Risk of Colorectal Cancer Based on Serologic and Mendelian Randomization Analyses.

BACKGROUND & AIMS: Human studies examining associations between circulating levels of insulin-like growth factor 1 (IGF1) and insulin-like growth factor binding protein 3 (IGFBP3) and colorectal cancer risk have reported inconsistent results. We conducted complementary serologic and Mendelian randomization (MR) analyses to determine whether alterations in circulating levels of IGF1 or IGFBP3 are associated with colorectal cancer development. METHODS: Serum levels of IGF1 were measured in blood samples collected from 397,380 participants from the UK Biobank, from 2006 through 2010. Incident cancer cases and cancer cases recorded first in death certificates were identified through linkage to national cancer and death registries. Complete follow-up was available through March 31, 2016. For the MR analyses, we identified genetic variants associated with circulating levels of IGF1 and IGFBP3. The association of these genetic variants with colorectal cancer was examined with 2-sample MR methods using genome-wide association study consortia data (52,865 cases with colorectal cancer and 46,287 individuals without [controls]) RESULTS: After a median follow-up period of 7.1 years, 2665 cases of colorectal cancer were recorded. In a multivariable-adjusted model, circulating level of IGF1 associated with colorectal cancer risk (hazard ratio per 1 standard deviation increment of IGF1, 1.11; 95% confidence interval [CI] 1.05-1.17). Similar associations were found by sex, follow-up time, and tumor subsite. In the MR analyses, a 1 standard deviation increment in IGF1 level, predicted based on genetic factors, was associated with a higher risk of colorectal cancer risk (odds ratio 1.08; 95% CI 1.03-1.12; P = 3.3 × 10-4). Level of IGFBP3, predicted based on genetic factors, was associated with colorectal cancer risk (odds ratio per 1 standard deviation increment, 1.12; 95% CI 1.06-1.18; P = 4.2 × 10-5). Colorectal cancer risk was associated with only 1 variant in the IGFBP3 gene region (rs11977526), which also associated with anthropometric traits and circulating level of IGF2. CONCLUSIONS: In an analysis of blood samples from almost 400,000 participants in the UK Biobank, we found an association between circulating level of IGF1 and colorectal cancer. Using genetic data from 52,865 cases with colorectal cancer and 46,287 controls, a higher level of IGF1, determined by genetic factors, was associated with colorectal cancer. Further studies are needed to determine how this signaling pathway might contribute to colorectal carcinogenesis.

A latent unknown clustering integrating multi-omics data (LUCID) with phenotypic traits.

MOTIVATION: Epidemiologic, clinical and translational studies are increasingly generating multiplatform omics data. Methods that can integrate across multiple high-dimensional data types while accounting for differential patterns are critical for uncovering novel associations and underlying relevant subgroups. RESULTS: We propose an integrative model to estimate latent unknown clusters (LUCID) aiming to both distinguish unique genomic, exposure and informative biomarkers/omic effects while jointly estimating subgroups relevant to the outcome of interest. Simulation studies indicate that we can obtain consistent estimates reflective of the true simulated values, accurately estimate subgroups and recapitulate subgroup-specific effects. We also demonstrate the use of the integrated model for future prediction of risk subgroups and phenotypes. We apply this approach to two real data applications to highlight the integration of genomic, exposure and metabolomic data. AVAILABILITY AND IMPLEMENTATION: The LUCID method is implemented through the LUCIDus R package available on CRAN (https://CRAN.R-project.org/package=LUCIDus). SUPPLEMENTARY INFORMATION: Supplementary materials are available at Bioinformatics online.

Genome-wide association study of circulating folate one-carbon metabolites.

Experimental, observational, and clinical trials support a critical role of folate one-carbon metabolism (FOCM) in colorectal cancer (CRC) development. In this report, we focus on understanding the relationship between common genetic variants and metabolites of FOCM. We conducted a genome-wide association study of FOCM biomarkers among 1,788 unaffected (without CRC) individuals of European ancestry from the Colon Cancer Family Registry. Twelve metabolites, including 5-methyltetrahydrofolate, vitamin B2 (flavin mononucleotide and riboflavin), vitamin B6 (4-pyridoxic acid, pyridoxal, and pyridoxamine), total homocysteine, methionine, S-adenosylmethionine, S-adenosylhomocysteine, cystathionine, and creatinine were measured from plasma using liquid chromatography-mass spectrometry (LC-MS) or LC-MS/MS. For each individual biomarker, we estimated genotype array-specific associations followed by a fixed-effect meta-analysis. We identified the variant rs35976024 (at 2p11.2 and intronic of ATOH8) associated with total homocysteine (p = 4.9 × 10-8 ). We found a group of six highly correlated variants on chromosome 15q14 associated with cystathionine (all p < 5 × 10-8 ), with the most significant variant rs28391580 (p = 2.8 × 10-8 ). Two variants (rs139435405 and rs149119426) on chromosome 14q13 showed significant (p < 5 × 10-8 ) associations with S-adenosylhomocysteine. These three biomarkers with significant associations are closely involved in homocysteine metabolism. Furthermore, when assessing the principal components (PCs) derived from seven individual biomarkers, we identified the variant rs12665366 (at 6p25.3 and intronic of EXOC2) associated with the first PC (p = 2.3 × 10-8 ). Our data suggest that common genetic variants may play an important role in FOCM, particularly in homocysteine metabolism.

DNA repair and cancer in colon and rectum: Novel players in genetic susceptibility.

Interindividual differences in DNA repair systems may play a role in modulating the individual risk of developing colorectal cancer. To better ascertain the role of DNA repair gene polymorphisms on colon and rectal cancer risk individually, we evaluated 15,419 single nucleotide polymorphisms (SNPs) within 185 DNA repair genes using GWAS data from the Colon Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), which included 8,178 colon cancer, 2,936 rectum cancer cases and 14,659 controls. Rs1800734 (in MLH1 gene) was associated with colon cancer risk (p-value = 3.5 × 10-6 ) and rs2189517 (in RAD51B) with rectal cancer risk (p-value = 5.7 × 10-6 ). The results had statistical significance close to the Bonferroni corrected p-value of 5.8 × 10-6 . Ninety-four SNPs were significantly associated with colorectal cancer risk after Binomial Sequential Goodness of Fit (BSGoF) procedure and confirmed the relevance of DNA mismatch repair (MMR) and homologous recombination pathways for colon and rectum cancer, respectively. Defects in MMR genes are known to be crucial for familial form of colorectal cancer but our findings suggest that specific genetic variations in MLH1 are important also in the individual predisposition to sporadic colon cancer. Other SNPs associated with the risk of colon cancer (e.g., rs16906252 in MGMT) were found to affect mRNA expression levels in colon transverse and therefore working as possible cis-eQTL suggesting possible mechanisms of carcinogenesis.

Large-Scale Genome-Wide Association Study of East Asians Identifies Loci Associated With Risk for Colorectal Cancer.

BACKGROUND & AIMS: Genome-wide association studies (GWASs) have associated approximately 50 loci with risk of colorectal cancer (CRC)-nearly one third of these loci were initially associated with CRC in studies conducted in East Asian populations. We conducted a GWAS of East Asians to identify CRC risk loci and evaluate the generalizability of findings from GWASs of European populations to Asian populations. METHODS: We analyzed genetic data from 22,775 patients with CRC (cases) and 47,731 individuals without cancer (controls) from 14 studies in the Asia Colorectal Cancer Consortium. First, we performed a meta-analysis of 7 GWASs (10,625 cases and 34,595 controls) and identified 46,554 promising risk variants for replication by adding them to the Multi-Ethnic Global Array (MEGA) for genotype analysis in 6445 cases and 7175 controls. These data were analyzed, along with data from an additional 5705 cases and 5961 controls genotyped using the OncoArray. We also obtained data from 57,976 cases and 67,242 controls of European descent. Variants at identified risk loci were functionally annotated and evaluated in correlation with gene expression levels. RESULTS: A meta-analyses of all samples from people of Asian descent identified 13 loci and 1 new variant at a known locus (10q24.2) associated with risk of CRC at the genome-wide significance level of P < 5 × 10-8. We did not perform experiments to replicate these associations in additional individuals of Asian ancestry. However, the lead risk variant in 6 of these loci was also significantly associated with risk of CRC in European descendants. A strong association (44%-75% increase in risk per allele) was found for 2 low-frequency variants: rs201395236 at 1q44 (minor allele frequency, 1.34%) and rs77969132 at 12p11.21 (minor allele frequency, 1.53%). For 8 of the 13 associated loci, the variants with the highest levels of significant association were located inside or near the protein-coding genes L1TD1, EFCAB2, PPP1R21, SLCO2A1, HLA-G, NOTCH4, DENND5B, and GNAS. For other intergenic loci, we provided evidence for the possible involvement of the genes ALDH7A1, PRICKLE1, KLF5, WWOX, and GLP2R. We replicated findings for 41 of 52 previously reported risk loci. CONCLUSIONS: We showed that most of the risk loci previously associated with CRC risk in individuals of European descent were also associated with CRC risk in East Asians. Furthermore, we identified 13 loci significantly associated with risk for CRC in Asians. Many of these loci contained genes that regulate the immune response, Wnt signaling to ß-catenin, prostaglandin E2 catabolism, and cell pluripotency and proliferation. Further analyses of these genes and their variants is warranted, particularly for the 8 loci for which the lead CRC risk variants were not replicated in persons of European descent.

Circulating bilirubin levels and risk of colorectal cancer: serological and Mendelian randomization analyses.

BACKGROUND: Bilirubin, a byproduct of hemoglobin breakdown and purported anti-oxidant, is thought to be cancer preventive. We conducted complementary serological and Mendelian randomization (MR) analyses to investigate whether alterations in circulating levels of bilirubin are associated with risk of colorectal cancer (CRC). We decided a priori to perform analyses separately in men and women based on suggestive evidence that associations may differ by sex. METHODS: In a case-control study nested in the European Prospective Investigation into Cancer and Nutrition (EPIC), pre-diagnostic unconjugated bilirubin (UCB, the main component of total bilirubin) concentrations were measured by high-performance liquid chromatography in plasma samples of 1386 CRC cases and their individually matched controls. Additionally, 115 single-nucleotide polymorphisms (SNPs) robustly associated (P < 5 × 10-8) with circulating total bilirubin were instrumented in a 2-sample MR to test for a potential causal effect of bilirubin on CRC risk in 52,775 CRC cases and 45,940 matched controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colon Cancer Family Registry (CCFR), and the Colorectal Transdisciplinary (CORECT) study. RESULTS: The associations between circulating UCB levels and CRC risk differed by sex (Pheterogeneity = 0.008). Among men, higher levels of UCB were positively associated with CRC risk (odds ratio [OR] = 1.19, 95% confidence interval [CI] = 1.04-1.36; per 1-SD increment of log-UCB). In women, an inverse association was observed (OR = 0.86 (0.76-0.97)). In the MR analysis of the main UGT1A1 SNP (rs6431625), genetically predicted higher levels of total bilirubin were associated with a 7% increase in CRC risk in men (OR = 1.07 (1.02-1.12); P = 0.006; per 1-SD increment of total bilirubin), while there was no association in women (OR = 1.01 (0.96-1.06); P = 0.73). Raised bilirubin levels, predicted by instrumental variables excluding rs6431625, were suggestive of an inverse association with CRC in men, but not in women. These differences by sex did not reach formal statistical significance (Pheterogeneity ≥ 0.2). CONCLUSIONS: Additional insight into the relationship between circulating bilirubin and CRC is needed in order to conclude on a potential causal role of bilirubin in CRC development.

Adiposity, metabolites, and colorectal cancer risk: Mendelian randomization study.

BACKGROUND: Higher adiposity increases the risk of colorectal cancer (CRC), but whether this relationship varies by anatomical sub-site or by sex is unclear. Further, the metabolic alterations mediating the effects of adiposity on CRC are not fully understood. METHODS: We examined sex- and site-specific associations of adiposity with CRC risk and whether adiposity-associated metabolites explain the associations of adiposity with CRC. Genetic variants from genome-wide association studies of body mass index (BMI) and waist-to-hip ratio (WHR, unadjusted for BMI; N = 806,810), and 123 metabolites from targeted nuclear magnetic resonance metabolomics (N = 24,925), were used as instruments. Sex-combined and sex-specific Mendelian randomization (MR) was conducted for BMI and WHR with CRC risk (58,221 cases and 67,694 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium, Colorectal Cancer Transdisciplinary Study, and Colon Cancer Family Registry). Sex-combined MR was conducted for BMI and WHR with metabolites, for metabolites with CRC, and for BMI and WHR with CRC adjusted for metabolite classes in multivariable models. RESULTS: In sex-specific MR analyses, higher BMI (per 4.2 kg/m2) was associated with 1.23 (95% confidence interval (CI) = 1.08, 1.38) times higher CRC odds among men (inverse-variance-weighted (IVW) model); among women, higher BMI (per 5.2 kg/m2) was associated with 1.09 (95% CI = 0.97, 1.22) times higher CRC odds. WHR (per 0.07 higher) was more strongly associated with CRC risk among women (IVW OR = 1.25, 95% CI = 1.08, 1.43) than men (IVW OR = 1.05, 95% CI = 0.81, 1.36). BMI or WHR was associated with 104/123 metabolites at false discovery rate-corrected P ≤ 0.05; several metabolites were associated with CRC, but not in directions that were consistent with the mediation of positive adiposity-CRC relations. In multivariable MR analyses, associations of BMI and WHR with CRC were not attenuated following adjustment for representative metabolite classes, e.g., the univariable IVW OR for BMI with CRC was 1.12 (95% CI = 1.00, 1.26), and this became 1.11 (95% CI = 0.99, 1.26) when adjusting for cholesterol in low-density lipoprotein particles. CONCLUSIONS: Our results suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. Adiposity was associated with numerous metabolic alterations, but none of these explained associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.

Type 2 diabetes mellitus, blood cholesterol, triglyceride and colorectal cancer risk in Lynch syndrome.

BACKGROUND: Type 2 diabetes mellitus and high total cholesterol and triglycerides are known to be associated with increased colorectal cancer risk for the general population. These associations are unknown for people with a germline DNA mismatch repair gene mutation (Lynch syndrome), who are at high risk of colorectal cancer. METHODS: This study included 2023 (56.4% female) carriers with a mismatch repair gene mutation (737 in MLH1, 928 in MSH2, 230 in MSH6, 106 in PMS2, 22 in EPCAM) recruited by the Colon Cancer Family Registry between 1998 and 2012. Weighted Cox regression was used to estimate the hazard ratios (HR) and 95% confidence intervals (CI) for the associations between self-reported type 2 diabetes, high cholesterol, triglyceride and colorectal cancer risk. RESULTS: Overall, 802 carriers were diagnosed with colorectal cancer at a median age of 42 years. A higher risk of colorectal cancer was observed in those with self-reported type-2 diabetes (HR 1.92; 95% CI, 1.03-3.58) and high cholesterol (HR 1.76; CI 1.23-2.52) compared with those without these conditions. There was no evidence of high triglyceride being associated with colorectal cancer risk. CONCLUSION: For people with Lynch syndrome, self-reported type-2 diabetes mellitus and high cholesterol were associated with increased colorectal cancer risk.