Genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) that are associated with risk of colorectal cancer. Prior research has evaluated the presence of ...gene-environment interaction involving the first 10 identified susceptibility loci, but little work has been conducted on interaction involving SNPs at recently identified susceptibility loci, including: rs10911251, rs6691170, rs6687758, rs11903757, rs10936599, rs647161, rs1321311, rs719725, rs1665650, rs3824999, rs7136702, rs11169552, rs59336, rs3217810, rs4925386, and rs2423279.
Data on 9,160 cases and 9,280 controls from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) and Colon Cancer Family Registry (CCFR) were used to evaluate the presence of interaction involving the above-listed SNPs and sex, body mass index (BMI), alcohol consumption, smoking, aspirin use, postmenopausal hormone (PMH) use, as well as intake of dietary calcium, dietary fiber, dietary folate, red meat, processed meat, fruit, and vegetables. Interaction was evaluated using a fixed effects meta-analysis of an efficient Empirical Bayes estimator, and permutation was used to account for multiple comparisons.
None of the permutation-adjusted P values reached statistical significance.
The associations between recently identified genetic susceptibility loci and colorectal cancer are not strongly modified by sex, BMI, alcohol, smoking, aspirin, PMH use, and various dietary factors.
Results suggest no evidence of strong gene-environment interactions involving the recently identified 16 susceptibility loci for colorectal cancer taken one at a time.
Regular use of nonsteroidal anti-inflammatory drugs (NSAID) is associated with lower risk of colorectal cancer. Genome-wide interaction analysis on single variants (G × E) has identified several SNPs ...that may interact with NSAIDs to confer colorectal cancer risk, but variations in gene expression levels may also modify the effect of NSAID use. Therefore, we tested interactions between NSAID use and predicted gene expression levels in relation to colorectal cancer risk.
Genetically predicted gene expressions were tested for interaction with NSAID use on colorectal cancer risk among 19,258 colorectal cancer cases and 18,597 controls from 21 observational studies. A Mixed Score Test for Interactions (MiSTi) approach was used to jointly assess G × E effects which are modeled via fixed interaction effects of the weighted burden within each gene set (burden) and residual G × E effects (variance). A false discovery rate (FDR) at 0.2 was applied to correct for multiple testing.
Among the 4,840 genes tested, genetically predicted expression levels of four genes modified the effect of any NSAID use on colorectal cancer risk, including
(P
= 1.96 × 10
),
(P
= 2.3 × 10
),
(P
= 9.38 × 10
), and
(P
= 1.44 × 10
). There was a significant interaction between expression level of
and regular use of aspirin only on colorectal cancer risk (P
= 3.23 × 10
). No interactions were observed between predicted gene expression and nonaspirin NSAID use at FDR < 0.2.
By incorporating functional information, we discovered several novel genes that interacted with NSAID use.
These findings provide preliminary support that could help understand the chemopreventive mechanisms of NSAIDs on colorectal cancer.
Although genome-wide association studies (GWAS) have separately identified many genetic susceptibility loci for ulcerative colitis (UC), Crohn's disease (CD) and colorectal cancer (CRC), there has ...been no large-scale examination for pleiotropy, or shared genetic susceptibility, for these conditions. We used logistic regression modeling to examine the associations of 181 UC and CD susceptibility variants previously identified by GWAS with risk of CRC using data from the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. We also examined associations of significant variants with clinical and molecular characteristics in a subset of the studies. Among 11794 CRC cases and 14190 controls, rs11676348, the susceptibility single nucleotide polymorphism (SNP) for UC, was significantly associated with reduced risk of CRC (P = 7E-05). The multivariate-adjusted odds ratio of CRC with each copy of the T allele was 0.93 (95% CI 0.89-0.96). The association of the SNP with risk of CRC differed according to mucinous histological features (P heterogeneity = 0.008). In addition, the (T) allele was associated with lower risk of tumors with Crohn's-like reaction but not tumors without such immune infiltrate (P heterogeneity = 0.02) and microsatellite instability-high (MSI-high) but not microsatellite stable or MSI-low tumors (P heterogeneity = 0.03). The minor allele (T) in SNP rs11676348, located downstream from CXCR2 that has been implicated in CRC progression, is associated with a lower risk of CRC, particularly tumors with a mucinous component, Crohn's-like reaction and MSI-high. Our findings offer the promise of risk stratification of inflammatory bowel disease patients for complications such as CRC.
Over 50 loci associated with colorectal cancer (CRC) have been uncovered by genome-wide association studies (GWAS). Identifying additional loci has the potential to help elucidate aspects of the ...underlying biological processes leading to better understanding of the pathogenesis of the disease. We re-evaluated a GWAS by excluding controls that have family history of CRC or personal history of colorectal polyps, as we hypothesized that their inclusion reduces power to detect associations. This is supported empirically and through simulations. Two-phase GWAS analysis was performed in a total of 16,517 cases and 14,487 controls. We identified rs17094983, a SNP associated with risk of CRC
p
= 2.5 × 10
−10
; odds ratio estimated by re-including all controls (OR) = 0.87, 95 % confidence interval (CI) 0.83–0.91; minor allele frequency (MAF) = 13 %. Results were replicated in samples of African descent (1894 cases and 4703 controls;
p
= 0.01; OR = 0.86, 95 % CI 0.77–0.97; MAF = 16 %). Gene expression data in 195 colon adenocarcinomas and 59 normal colon tissues from two different studies revealed that this locus has genotypes that are associated with RTN1 (Reticulon 1) expression (
p
= 0.001), a protein-coding gene involved in survival and proliferation of cancer cells which is highly expressed in normal colon tissues but has significantly reduced expression in tumor cells (
p
= 1.3 × 10
−8
).
Circulating adiponectin has been associated with lower risk of colorectal cancer (CRC). Genome‐wide association studies have identified several single‐nucleotide polymorphisms (SNPs) associated with ...adiponectin levels. However, it is unclear whether these SNPs are associated with CRC risk. In addition, previous data on SNPs in the adiponectin pathway and their associations with CRC are inconsistent. Therefore, we examined 19 SNPs in genes related to adiponectin or its receptors and their associations with CRC using logistic regression among 7,020 cases and 7,631 controls drawn from ten studies included in the Genetics and Epidemiology of Colorectal Cancer Consortium. Using data from a subset of two large cohort studies, we also assessed the contribution of individual SNPs and an adiponectin genetic score to plasma adiponectin after accounting for lifestyle factors among 2,217 women and 619 men. We did not find any statistically significant association between the 19 adiponectin‐associated SNPs and CRC risk (multivariable‐adjusted odds ratios ranged from 0.89 to 1.05, all p > 0.05). Each SNP explained less than 2.50% of the variance of plasma adiponectin, and the genetic score collectively accounted for 2.95 and 1.42% of the variability of adiponectin in women and men, respectively, after adjustment for age, body mass index, physical activity, smoking, alcohol consumption, regular use of aspirin or nonsteroidal anti‐inflammatory drug and postmenopausal hormone use. In conclusion, our findings do not support an association between known adiponectin‐related common SNPs and CRC incidence. However, known common SNPs account for only a limited proportion of the interindividual variance in circulating adiponectin. Further work is warranted to investigate the relationship between adiponectin and CRC while accounting for other components in the pathway.
What's new?
Obesity is a well‐established risk factor for colorectal cancer, and mounting evidence suggests that the adipocyte‐derived hormone adiponectin is behind that link. In the present investigation, common genetic variants in the adiponectin (ADIPOQ) gene were analyzed for potential relationships with colorectal cancer risk. However, no statistically significant associations were identified. Because common genetic determinants account for only a limited proportion of the variance in circulating adiponectin concentrations, factors that influence adiponectin receptors and other components of the adiponectin pathway should be investigated.
The evaluation of less frequent genetic variants and their effect on complex disease pose new challenges for genomic research. To investigate whether epigenetic data can be used to inform aggregate ...rare-variant association methods (RVAM), we assessed whether variants more significantly associated with colorectal cancer (CRC) were preferentially located in non-coding regulatory regions, and whether enrichment was specific to colorectal tissues.
Active regulatory elements (ARE) were mapped using data from 127 tissues and cell-types from NIH Roadmap Epigenomics and Encyclopedia of DNA Elements (ENCODE) projects. We investigated whether CRC association p-values were more significant for common variants inside versus outside AREs, or 2) inside colorectal (CR) AREs versus AREs of other tissues and cell-types. We employed an integrative epigenomic RVAM for variants with allele frequency <1%. Gene sets were defined as ARE variants within 200 kilobases of a transcription start site (TSS) using either CR ARE or ARE from non-digestive tissues. CRC-set association p-values were used to evaluate enrichment of less frequent variant associations in CR ARE versus non-digestive ARE.
ARE from 126/127 tissues and cell-types were significantly enriched for stronger CRC-variant associations. Strongest enrichment was observed for digestive tissues and immune cell types. CR-specific ARE were also enriched for stronger CRC-variant associations compared to ARE combined across non-digestive tissues (p-value = 9.6 × 10-4). Additionally, we found enrichment of stronger CRC association p-values for rare variant sets of CR ARE compared to non-digestive ARE (p-value = 0.029).
Integrative epigenomic RVAM may enable discovery of less frequent variants associated with CRC, and ARE of digestive and immune tissues are most informative. Although distance-based aggregation of less frequent variants in CR ARE surrounding TSS showed modest enrichment, future association studies would likely benefit from joint analysis of transcriptomes and epigenomes to better link regulatory variation with target genes.
Menopausal hormone therapy (MHT) use has been consistently associated with a decreased risk of colorectal cancer (CRC) in women. Our aim was to use a genome-wide gene-environment interaction analysis ...to identify genetic modifiers of CRC risk associated with use of MHT.
We included 10 835 postmenopausal women (5419 cases and 5416 controls) from 10 studies. We evaluated use of any MHT, oestrogen-only (E-only) and combined oestrogen-progestogen (E+P) hormone preparations. To test for multiplicative interactions, we applied the empirical Bayes (EB) test as well as the Wald test in conventional case-control logistic regression as primary tests. The Cocktail test was used as secondary test.
The EB test identified a significant interaction between rs964293 at 20q13.2/CYP24A1 and E+P (interaction OR (95% CIs)=0.61 (0.52-0.72), P=4.8 × 10(-9)). The secondary analysis also identified this interaction (Cocktail test OR=0.64 (0.52-0.78), P=1.2 × 10(-5) (alpha threshold=3.1 × 10(-4)). The ORs for association between E+P and CRC risk by rs964293 genotype were as follows: C/C, 0.96 (0.61-1.50); A/C, 0.61 (0.39-0.95) and A/A, 0.40 (0.22-0.73), respectively.
Our results indicate that rs964293 modifies the association between E+P and CRC risk. The variant is located near CYP24A1, which encodes an enzyme involved in vitamin D metabolism. This novel finding offers additional insight into downstream pathways of CRC etiopathogenesis.
Gene expression data are often collected from tissue samples that are composed of multiple cell types. Studies of cell type composition based on gene expression data from tissue samples have recently ...attracted increasing research interest and led to new method development for cell type composition estimation. This new information on cell type composition can be associated with individual characteristics (e.g., genetic variants) or clinical outcomes (e.g., survival time). Such association analysis can be conducted for each cell type separately followed by multiple testing correction. An alternative approach is to evaluate this association using the composition of all the cell types, thus aggregating association signals across cell types. A key challenge of this approach is to account for the dependence across cell types. We propose a new method to quantify the distances between cell types while accounting for their dependencies, and use this information for association analysis. We demonstrate our method in two applied examples: to assess the association between immune cell type composition in tumor samples of colorectal cancer patients versus survival time and SNP genotypes. We found immune cell composition has prognostic value, and our distance metric leads to more accurate survival time prediction than other distance metrics that ignore cell type dependencies. In addition, survival time-associated SNPs are enriched among the SNPs associated with immune cell composition.
Smoking is associated with colorectal cancer (CRC) risk. Previous studies suggested this association may be restricted to certain molecular subtypes of CRC, but large-scale comprehensive analysis is ...lacking.
A total of 9789 CRC cases and 11 231 controls of European ancestry from 11 observational studies were included. We harmonized smoking variables across studies and derived sex study-specific quartiles of pack-years of smoking for analysis. Four somatic colorectal tumor markers were assessed individually and in combination, including
mutation,
mutation, CpG island methylator phenotype (CIMP), and microsatellite instability (MSI) status. A multinomial logistic regression analysis was used to assess the association between smoking and risk of CRC subtypes by molecular characteristics, adjusting for age, sex, and study. All statistical tests were 2-sided and adjusted for Bonferroni correction.
Heavier smoking was associated with higher risk of CRC overall and stratified by individual markers (
< .001). The associations differed statistically significantly between all molecular subtypes, which was the most statistically significant for CIMP and
. Compared with never-smokers, smokers in the fourth quartile of pack-years had a 90% higher risk of CIMP-positive CRC (odds ratio = 1.90, 95% confidence interval = 1.60 to 2.26) but only 35% higher risk for CIMP-negative CRC (odds ratio = 1.35, 95% confidence interval = 1.22 to 1.49;
= 2.1 x 10
). The association was also stronger in tumors that were
positive, MSI high, or
wild type when combined (
< .001).
Smoking was associated with differential risk of CRC subtypes defined by molecular characteristics. Heavier smokers had particularly higher risk of CRC subtypes that were CIMP positive and MSI high in combination, suggesting that smoking may be involved in the development of colorectal tumors via the serrated pathway.
Genome-wide association studies (GWAS) have successfully identified a number of single-nucleotide polymorphisms (SNPs) associated with colorectal cancer (CRC) risk. However, these susceptibility loci ...known today explain only a small fraction of the genetic risk. Gene-gene interaction (GxG) is considered to be one source of the missing heritability. To address this, we performed a genome-wide search for pair-wise GxG associated with CRC risk using 8,380 cases and 10,558 controls in the discovery phase and 2,527 cases and 2,658 controls in the replication phase. We developed a simple, but powerful method for testing interaction, which we term the Average Risk Due to Interaction (ARDI). With this method, we conducted a genome-wide search to identify SNPs showing evidence for GxG with previously identified CRC susceptibility loci from 14 independent regions. We also conducted a genome-wide search for GxG using the marginal association screening and examining interaction among SNPs that pass the screening threshold (p<10(-4)). For the known locus rs10795668 (10p14), we found an interacting SNP rs367615 (5q21) with replication p = 0.01 and combined p = 4.19×10(-8). Among the top marginal SNPs after LD pruning (n = 163), we identified an interaction between rs1571218 (20p12.3) and rs10879357 (12q21.1) (nominal combined p = 2.51×10(-6); Bonferroni adjusted p = 0.03). Our study represents the first comprehensive search for GxG in CRC, and our results may provide new insight into the genetic etiology of CRC.