CYP24A1 variant modifies the association between use of oestrogen plus progestogen therapy and colorectal cancer risk

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

Association of aspirin and NSAID use with risk of colorectal cancer according to genetic variants

IMPORTANCE

Use of aspirin and other nonsteroidal anti-inflammatory drugs (NSAIDs) is associated with lower risk of colorectal cancer.

OBJECTIVE

To identify common genetic markers that may confer differential benefit from aspirin or NSAID chemoprevention, we tested gene × environment interactions between regular use of aspirin and/or NSAIDs and single-nucleotide polymorphisms (SNPs) in relation to risk of colorectal cancer.

DESIGN, SETTING, AND PARTICIPANTS

Case-control study using data from 5 case-control and 5 cohort studies initiated between 1976 and 2003 across the United States, Canada, Australia, and Germany and including colorectal cancer cases (n=8634) and matched controls (n=8553) ascertained between 1976 and 2011. Participants were all of European descent.

EXPOSURES

Genome-wide SNP data and information on regular use of aspirin and/or NSAIDs and other risk factors.

MAIN OUTCOMES AND MEASURES

Colorectal cancer.

RESULTS

Regular use of aspirin and/or NSAIDs was associated with lower risk of colorectal cancer (prevalence, 28% vs 38%; odds ratio [OR], 0.69 [95% CI, 0.64-0.74]; P = 6.2 × 10(-28)) compared with nonregular use. In the conventional logistic regression analysis, the SNP rs2965667 at chromosome 12p12.3 near the MGST1 gene showed a genome-wide significant interaction with aspirin and/or NSAID use (P = 4.6 × 10(-9) for interaction). Aspirin and/or NSAID use was associated with a lower risk of colorectal cancer among individuals with rs2965667-TT genotype (prevalence, 28% vs 38%; OR, 0.66 [95% CI, 0.61-0.70]; P = 7.7 × 10(-33)) but with a higher risk among those with rare (4%) TA or AA genotypes (prevalence, 35% vs 29%; OR, 1.89 [95% CI, 1.27-2.81]; P = .002). In case-only interaction analysis, the SNP rs16973225 at chromosome 15q25.2 near the IL16 gene showed a genome-wide significant interaction with use of aspirin and/or NSAIDs (P = 8.2 × 10(-9) for interaction). Regular use was associated with a lower risk of colorectal cancer among individuals with rs16973225-AA genotype (prevalence, 28% vs 38%; OR, 0.66 [95% CI, 0.62-0.71]; P = 1.9 × 10(-30)) but was not associated with risk of colorectal cancer among those with less common (9%) AC or CC genotypes (prevalence, 36% vs 39%; OR, 0.97 [95% CI, 0.78-1.20]; P = .76).

CONCLUSIONS AND RELEVANCE

In this genome-wide investigation of gene × environment interactions, use of aspirin and/or NSAIDs was associated with lower risk of colorectal cancer, and this association differed according to genetic variation at 2 SNPs at chromosomes 12 and 15. Validation of these findings in additional populations may facilitate targeted colorectal cancer prevention strategies.

Functional annotation of putative regulatory elements at cancer susceptibility Loci

Most cancer-associated genetic variants identified from genome-wide association studies (GWAS) do not obviously change protein structure, leading to the hypothesis that the associations are attributable to regulatory polymorphisms. Translating genetic associations into mechanistic insights can be facilitated by knowledge of the causal regulatory variant (or variants) responsible for the statistical signal. Experimental validation of candidate functional variants is onerous, making bioinformatic approaches necessary to prioritize candidates for laboratory analysis. Thus, a systematic approach for recognizing functional (and, therefore, likely causal) variants in noncoding regions is an important step toward interpreting cancer risk loci. This review provides a detailed introduction to current regulatory variant annotations, followed by an overview of how to leverage these resources to prioritize candidate functional polymorphisms in regulatory regions.

Whole-exome imputation of sequence variants identified two novel alleles associated with adult body height in African Americans

Adult body height is a quantitative trait for which genome-wide association studies (GWAS) have identified numerous loci, primarily in European populations. These loci, comprising common variants, explain <10% of the phenotypic variance in height. We searched for novel associations between height and common (minor allele frequency, MAF ≥5%) or infrequent (0.5% < MAF < 5%) variants across the exome in African Americans. Using a reference panel of 1692 African Americans and 471 Europeans from the National Heart, Lung, and Blood Institute's (NHLBI) Exome Sequencing Project (ESP), we imputed whole-exome sequence data into 13 719 African Americans with existing array-based GWAS data (discovery). Variants achieving a height-association threshold of P < 5E-06 in the imputed dataset were followed up in an independent sample of 1989 African Americans with whole-exome sequence data (replication). We used P < 2.5E-07 (=0.05/196 779 variants) to define statistically significant associations in meta-analyses combining the discovery and replication sets (N = 15 708). We discovered and replicated three independent loci for association: 5p13.3/C5orf22/rs17410035 (MAF = 0.10, β = 0.64 cm, P = 8.3E-08), 13q14.2/SPRYD7/rs114089985 (MAF = 0.03, β = 1.46 cm, P = 4.8E-10) and 17q23.3/GH2/rs2006123 (MAF = 0.30; β = 0.47 cm; P = 4.7E-09). Conditional analyses suggested 5p13.3 (C5orf22/rs17410035) and 13q14.2 (SPRYD7/rs114089985) may harbor novel height alleles independent of previous GWAS-identified variants (r(2) with GWAS loci <0.01); whereas 17q23.3/GH2/rs2006123 was correlated with GWAS-identified variants in European and African populations. Notably, 13q14.2/rs114089985 is infrequent in African Americans (MAF = 3%), extremely rare in European Americans (MAF = 0.03%), and monomorphic in Asian populations, suggesting it may be an African-American-specific height allele. Our findings demonstrate that whole-exome imputation of sequence variants can identify low-frequency variants and discover novel variants in non-European populations.

Identification of Genetic Susceptibility Loci for Colorectal Tumors in a Genome-Wide Meta-analysis

BACKGROUND & AIMS

Heritable factors contribute to the development of colorectal cancer. Identifying the genetic loci associated with colorectal tumor formation could elucidate the mechanisms of pathogenesis.

METHODS

We conducted a genome-wide association study that included 14 studies, 12,696 cases of colorectal tumors (11,870 cancer, 826 adenoma), and 15,113 controls of European descent. The 10 most statistically significant, previously unreported findings were followed up in 6 studies; these included 3056 colorectal tumor cases (2098 cancer, 958 adenoma) and 6658 controls of European and Asian descent.

RESULTS

Based on the combined analysis, we identified a locus that reached the conventional genome-wide significance level at less than 5.0 × 10(-8): an intergenic region on chromosome 2q32.3, close to nucleic acid binding protein 1 (most significant single nucleotide polymorphism: rs11903757; odds ratio [OR], 1.15 per risk allele; P = 3.7 × 10(-8)). We also found evidence for 3 additional loci with P values less than 5.0 × 10(-7): a locus within the laminin gamma 1 gene on chromosome 1q25.3 (rs10911251; OR, 1.10 per risk allele; P = 9.5 × 10(-8)), a locus within the cyclin D2 gene on chromosome 12p13.32 (rs3217810 per risk allele; OR, 0.84; P = 5.9 × 10(-8)), and a locus in the T-box 3 gene on chromosome 12q24.21 (rs59336; OR, 0.91 per risk allele; P = 3.7 × 10(-7)).

CONCLUSIONS

In a large genome-wide association study, we associated polymorphisms close to nucleic acid binding protein 1 (which encodes a DNA-binding protein involved in DNA repair) with colorectal tumor risk. We also provided evidence for an association between colorectal tumor risk and polymorphisms in laminin gamma 1 (this is the second gene in the laminin family to be associated with colorectal cancers), cyclin D2 (which encodes for cyclin D2), and T-box 3 (which encodes a T-box transcription factor and is a target of Wnt signaling to β-catenin). The roles of these genes and their products in cancer pathogenesis warrant further investigation.

Genome-wide search for gene-gene interactions in colorectal cancer

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⁻⁴). For the known locus rs10795668 (10p14), we found an interacting SNP rs367615 (5q21) with replication p = 0.01 and combined p = 4.19×10⁻⁸. 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⁻⁶; 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.

Abstract PR3: Framework for post-GWAS functional annotation of regulatory regions associated with susceptibility loci for colorectal cancer

Genome-wide association (GWA) studies of colorectal cancer (CRC) have successfully replicated several susceptibility loci in populations of European ancestry. Most of the statistically associated tag SNPs do not correlate with any coding variants suggesting that the underlying functional variants responsible for the associations likely impact regulatory regions. Advancing basic research and translating GWAS findings into clinical benefit will be facilitated by identifying causal variants and their associated biological mechanism.

We developed a ‘post-GWAS’ strategy to functionally dissect susceptibility regions that combines in silico prediction of functional alleles with a novel allelic reporter assay to quantify effects on gene expression. Functional predictions were made by compiling a list of correlated SNPs for each published CRC susceptibility loci and annotating each list for potential regulatory evidence consistent with enhancers, promoters, insulators and silencers. Using the UCSC Genome Browser, we aligned each variant with a combined browser view of several ENCODE tracks and compared allelic regions in JASPAR and ConSite to potentially identify variants that alter transcription factor binding sites (TFBS). A transparent overlay of various cell lines in the Histone Modifications Track was used to characterize functional regulation at low resolution. The DNase hypersensitivity track provided a more precise demarcation of the region for plasmid design. Evidence of altered binding sites through ENCODE's ChIP-Seq TFBS track, JASPAR and ConSite helped prioritize variants for further functional analysis. Use of 46-way PhastCons track in the Genome Browser was used as secondary evidence for a regulatory region, but lack of conservation did not rule out a candidate. Priority for in vitro analysis was based on strength of the functional annotation and ability to identify CRC associated genes that could reasonably be regulated by these predicted regulatory elements. Our in silico analysis identified 2 strong candidate SNPs falling in the promoter regions of GREM1 and CDH1, as well as, 4 probable candidates in the putative enhancer regions of CDH1 and BMP2.

In order to test our in silico functional candidates, we designed allelic promoter and enhancer constructs for each candidate variant and inserted them into custom pCI-neo vectors to drive expression of the ZSgreen reporter gene. These plasmids were semi-stably co-transfected into a colon cancer cell line, and after 14 days of growth cells with integrated plasmid were Fluorescence-Activated Cell Sorted into 4 quartiles of ZSgreen expression. To quantify the relative activity of our allelic plasmids, each pool was sequenced with Illumina's high through-put sequencing technology. We tested for statistically significant log fold changes between each cotransfected plasmid. Due to interdependence of signal between alleles we tested differences using a Poisson Generalized Linear Model in a step-wise manner. Blinded modeling of this method revealed that this approach could be used to test at least 24 allelic plasmids simultaneously. Although this study focused on susceptibility loci for CRC, this framework could be extended to rapidly test many functional hypotheses in other complex genetic diseases providing a high throughput approach to gain substantial new insight into the functional effects of GWAS findings.

This proffered talk is also presented as Poster 63.

Citation Format: Stephanie A. Rosse, Christina T. Chen, Orsalem J. Kahsai, Ahmad S. Zebari, Paul W. Livermore Auer, Graham Casey, Ulrike Peters, Christopher S. Carlson. Framework for post-GWAS functional annotation of regulatory regions associated with susceptibility loci for colorectal cancer. [abstract]. In: Proceedings of the AACR Special Conference on Post-GWAS Horizons in Molecular Epidemiology: Digging Deeper into the Environment; 2012 Nov 11-14; Hollywood, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2012;21(11 Suppl):Abstract nr PR3.