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.

CFH and ARMS2 genetic risk determines progression to neovascular age-related macular degeneration after antioxidant and zinc supplementation.

We evaluated the influence of an antioxidant and zinc nutritional supplement [the Age-Related Eye Disease Study (AREDS) formulation] on delaying or preventing progression to neovascular AMD (NV) in persons with age-related macular degeneration (AMD). AREDS subjects (n = 802) with category 3 or 4 AMD at baseline who had been treated with placebo or the AREDS formulation were evaluated for differences in the risk of progression to NV as a function of complement factor H (CFH) and age-related maculopathy susceptibility 2 (ARMS2) genotype groups. We used published genetic grouping: a two-SNP haplotype risk-calling algorithm to assess CFH, and either the single SNP rs10490924 or 372_815del443ins54 to mark ARMS2 risk. Progression risk was determined using the Cox proportional hazard model. Genetics-treatment interaction on NV risk was assessed using a multiiterative bootstrap validation analysis. We identified strong interaction of genetics with AREDS formulation treatment on the development of NV. Individuals with high CFH and no ARMS2 risk alleles and taking the AREDS formulation had increased progression to NV compared with placebo. Those with low CFH risk and high ARMS2 risk had decreased progression risk. Analysis of CFH and ARMS2 genotype groups from a validation dataset reinforces this conclusion. Bootstrapping analysis confirms the presence of a genetics-treatment interaction and suggests that individual treatment response to the AREDS formulation is largely determined by genetics. The AREDS formulation modifies the risk of progression to NV based on individual genetics. Its use should be based on patient-specific genotype.

Genome-wide diet-gene interaction analyses for risk of colorectal cancer.

Dietary factors, including meat, fruits, vegetables and fiber, are associated with colorectal cancer; however, there is limited information as to whether these dietary factors interact with genetic variants to modify risk of colorectal cancer. We tested interactions between these dietary factors and approximately 2.7 million genetic variants for colorectal cancer risk among 9,287 cases and 9,117 controls from ten studies. We used logistic regression to investigate multiplicative gene-diet interactions, as well as our recently developed Cocktail method that involves a screening step based on marginal associations and gene-diet correlations and a testing step for multiplicative interactions, while correcting for multiple testing using weighted hypothesis testing. Per quartile increment in the intake of red and processed meat were associated with statistically significant increased risks of colorectal cancer and vegetable, fruit and fiber intake with lower risks. From the case-control analysis, we detected a significant interaction between rs4143094 (10p14/near GATA3) and processed meat consumption (OR = 1.17; p = 8.7E-09), which was consistently observed across studies (p heterogeneity = 0.78). The risk of colorectal cancer associated with processed meat was increased among individuals with the rs4143094-TG and -TT genotypes (OR = 1.20 and OR = 1.39, respectively) and null among those with the GG genotype (OR = 1.03). Our results identify a novel gene-diet interaction with processed meat for colorectal cancer, highlighting that diet may modify the effect of genetic variants on disease risk, which may have important implications for prevention.

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.

The effect of 5-fluorouracil/leucovorin chemotherapy on CpG methylation, or the confounding role of leukocyte heterogeneity: An illustration.

Blood-based epigenome-wide association studies that aim at comparing CpG methylation between colorectal cancer (CRC) patients and controls can lead to the discovery of diagnostic or prognostic biomarkers. Numerous confounders can lead to spurious associations. We aimed to see if 5-fluorouracil (5-FU)/leucovorin chemotherapy administered to cases prior to the collection of their blood has an effect on methylation. 304 patients who received treatment and 273 who did not were profiled on the HumanMethylation450 array. Association tests were adjusted for confounders, including proxies for leukocyte cell counts. There were substantial methylation differences between these two groups that vanished once the leukocyte heterogeneity was accounted for. We observed a significant decrease of T cells in the treatment group (CD4+: p=10(-6); CD8+: p=0.036) and significant increase of NK cells (p=0.05) and monocytes (p=0.0006). 5-FU/leucovorin has no effect on global and local blood-based methylation profiles, other than through differences in the leukocyte compositions that the treatment induced.

Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24.

Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage 1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 x 10(-8)) and 1.14 (trend; P = 1.32 x 10(-5)), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 x 10(-4)). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 x 10(-11)). This locus has also been implicated in prostate cancer.

Treatment response to antioxidants and zinc based on CFH and ARMS2 genetic risk allele number in the Age-Related Eye Disease Study.

OBJECTIVE: To evaluate the impact of complement factor H (CFH) and age-related maculopathy susceptibility 2 (ARMS2) risk alleles on the observed response to components of the Age-Related Eye Disease Study (AREDS) formulation. DESIGN: Genetic and statistical subgroup analysis of a randomized, prospective clinical trial. PARTICIPANTS: White patients from the AREDS with category 3 or 4 age-related macular degeneration (AMD) with available DNA (n = 989). METHODS: Four genotype groups based on CFH and ARMS2 risk allele number were defined. Progression to advanced AMD was analyzed by genotype and treatment using Cox proportionate hazards estimates and 7-year events. MAIN OUTCOME MEASURES: The effect of predefined genotype group on treatment-specific progression to advanced AMD. RESULTS: Patients with 2 CFH risk alleles and no ARMS2 risk alleles progressed more with zinc-containing treatment compared with placebo, with a hazard ratio (HR) of 3.07 (P = 0.0196) for zinc and 2.73 (P = 0.0418) for AREDS formulation (AF). Seven-year treatment-specific progression rates were: placebo, 17.0%; zinc, 43.2% (P = 0.023); and AF, 40.2% (P = 0.039). Patients with 0 or 1 CFH risk alleles and 1 or 2 ARMS2 risk alleles benefited from zinc-containing treatment compared with placebo, with an HR of 0.514 for zinc (P = 0.012) and 0.569 for AF (P = 0.0254). Seven-year treatment-specific AMD progression rates were as follows: placebo, 43.3%; zinc, 25.2% (P = 0.020); and AF, 27.3% (P = 0.011). Zinc and AF treatment each interacted statistically with these 2 genotype groups under a Cox model, with P values of 0.000999 and 0.00366, respectively. For patients with 0 or 1 CFH risk alleles and no ARMS2 risk alleles, neither zinc-containing treatment altered progression compared with placebo, but treatment with antioxidants decreased progression (HR, 0.380; P = 0.034). Seven-year progression with placebo was 22.6% and with antioxidants was 9.17% (P = 0.033). For patients with 2 CFH risk alleles and 1 or 2 ARMS2 risk alleles, no treatment was better than placebo (48.4%). CONCLUSIONS: The benefit of the AREDS formulation seems the result of a favorable response by patients in only 1 genotype group, balanced by neutral or unfavorable responses in 3 genotype groups.

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.

Pleiotropic effects of genetic risk variants for other cancers on colorectal cancer risk: PAGE, GECCO and CCFR consortia.

OBJECTIVE: Genome-wide association studies have identified a large number of single nucleotide polymorphisms (SNPs) associated with a wide array of cancer sites. Several of these variants demonstrate associations with multiple cancers, suggesting pleiotropic effects and shared biological mechanisms across some cancers. We hypothesised that SNPs previously associated with other cancers may additionally be associated with colorectal cancer. In a large-scale study, we examined 171 SNPs previously associated with 18 different cancers for their associations with colorectal cancer. DESIGN: We examined 13 338 colorectal cancer cases and 40 967 controls from three consortia: Population Architecture using Genomics and Epidemiology (PAGE), Genetic Epidemiology of Colorectal Cancer (GECCO), and the Colon Cancer Family Registry (CCFR). Study-specific logistic regression results, adjusted for age, sex, principal components of genetic ancestry, and/or study specific factors (as relevant) were combined using fixed-effect meta-analyses to evaluate the association between each SNP and colorectal cancer risk. A Bonferroni-corrected p value of 2.92×10(-4) was used to determine statistical significance of the associations. RESULTS: Two correlated SNPs--rs10090154 and rs4242382--in Region 1 of chromosome 8q24, a prostate cancer susceptibility region, demonstrated statistically significant associations with colorectal cancer risk. The most significant association was observed with rs4242382 (meta-analysis OR=1.12; 95% CI 1.07 to 1.18; p=1.74×10(-5)), which also demonstrated similar associations across racial/ethnic populations and anatomical sub-sites. CONCLUSIONS: This is the first study to clearly demonstrate Region 1 of chromosome 8q24 as a susceptibility locus for colorectal cancer; thus, adding colorectal cancer to the list of cancer sites linked to this particular multicancer risk region at 8q24.

Identification of genes expressed by immune cells of the colon that are regulated by colorectal cancer-associated variants.

A locus on human chromosome 11q23 tagged by marker rs3802842 was associated with colorectal cancer (CRC) in a genome-wide association study; this finding has been replicated in case-control studies worldwide. In order to identify biologic factors at this locus that are related to the etiopathology of CRC, we used microarray-based target selection methods, coupled to next-generation sequencing, to study 103 kb at the 11q23 locus. We genotyped 369 putative variants from 1,030 patients with CRC (cases) and 1,061 individuals without CRC (controls) from the Ontario Familial Colorectal Cancer Registry. Two previously uncharacterized genes, COLCA1 and COLCA2, were found to be co-regulated genes that are transcribed from opposite strands. Expression levels of COLCA1 and COLCA2 transcripts correlate with rs3802842 genotypes. In colon tissues, COLCA1 co-localizes with crystalloid granules of eosinophils and granular organelles of mast cells, neutrophils, macrophages, dendritic cells and differentiated myeloid-derived cell lines. COLCA2 is present in the cytoplasm of normal epithelial, immune and other cell lineages, as well as tumor cells. Tissue microarray analysis demonstrates the association of rs3802842 with lymphocyte density in the lamina propria (p = 0.014) and levels of COLCA1 in the lamina propria (p = 0.00016) and COLCA2 (tumor cells, p = 0.0041 and lamina propria, p = 6 × 10(-5)). In conclusion, genetic, expression and immunohistochemical data implicate COLCA1 and COLCA2 in the pathogenesis of colon cancer. Histologic analyses indicate the involvement of immune pathways.

Cumulative impact of common genetic variants and other risk factors on colorectal cancer risk in 42,103 individuals.

OBJECTIVE: Colorectal cancer (CRC) has a substantial heritable component. Common genetic variation has been shown to contribute to CRC risk. A study was conducted in a large multi-population study to assess the feasibility of CRC risk prediction using common genetic variant data combined with other risk factors. A risk prediction model was built and applied to the Scottish population using available data. DESIGN: Nine populations of European descent were studied to develop and validate CRC risk prediction models. Binary logistic regression was used to assess the combined effect of age, gender, family history (FH) and genotypes at 10 susceptibility loci that individually only modestly influence CRC risk. Risk models were generated from case-control data incorporating genotypes alone (n=39,266) and in combination with gender, age and FH (n=11,324). Model discriminatory performance was assessed using 10-fold internal cross-validation and externally using 4187 independent samples. The 10-year absolute risk was estimated by modelling genotype and FH with age- and gender-specific population risks. RESULTS: The median number of risk alleles was greater in cases than controls (10 vs 9, p<2.2 × 10(-16)), confirmed in external validation sets (Sweden p=1.2 × 10(-6), Finland p=2 × 10(-5)). The mean per-allele increase in risk was 9% (OR 1.09; 95% CI 1.05 to 1.13). Discriminative performance was poor across the risk spectrum (area under curve for genotypes alone 0.57; area under curve for genotype/age/gender/FH 0.59). However, modelling genotype data, FH, age and gender with Scottish population data shows the practicalities of identifying a subgroup with >5% predicted 10-year absolute risk. CONCLUSION: Genotype data provide additional information that complements age, gender and FH as risk factors, but individualised genetic risk prediction is not currently feasible. Nonetheless, the modelling exercise suggests public health potential since it is possible to stratify the population into CRC risk categories, thereby informing targeted prevention and surveillance.

CFH and ARMS2 genetic polymorphisms predict response to antioxidants and zinc in patients with age-related macular degeneration.

OBJECTIVE: The Age-Related Eye Disease Study (AREDS) demonstrated that antioxidant and zinc supplementation decreases progression to advanced age-related macular degeneration (AMD) in patients with moderate to severe disease. We evaluated the interaction of genetics and type of nutritional supplement on progression from moderate to advanced AMD. DESIGN: Genetic analysis of a randomized, prospective clinical trial. PARTICIPANTS: White patients with AREDS category 3 AMD in 1 eye and AREDS categories 1 through 4 AMD in the fellow eye enrolled in the AREDS with available peripheral blood-derived DNA (995). METHODS: Subjects were evaluated for known AMD genetic risk markers and treatment category. The progression rate to advanced AMD was analyzed by genotypes and AREDS treatment group using Cox regression. MAIN OUTCOME MEASURES: The effect of inherited gene polymorphisms on treatment group-specific rate of progression to advanced AMD. RESULTS: Over an average of 10.1 years, individuals with 1 or 2 complement factor H (CFH) risk alleles derived maximum benefit from antioxidants alone. In these patients, the addition of zinc negated the benefits of antioxidants. Treatment with zinc and antioxidants was associated with a risk ratio (RR) of 1.83 with 2 CFH risk alleles (P = 1.03E-02), compared with outcomes for patients without CFH risk alleles. Patients with age-related maculopathy sensitivity 2 (ARMS2) risk alleles derived maximum benefit from zinc-containing regimens, with a deleterious response to antioxidants in the presence of ARMS2 risk alleles. Treatment with antioxidants was associated with an RR of 2.58 for those with 1 ARMS2 risk allele and 3.96 for those with 2 ARMS2 risk alleles (P = 1.04E-6), compared with patients with no ARMS2 risk alleles. Individuals homozygous for CFH and ARMS2 risk alleles derived no benefit from any category of AREDS treatment. CONCLUSIONS: Individuals with moderate AMD could benefit from pharmacogenomic selection of nutritional supplements. In this analysis, patients with no CFH risk alleles and with 1 or 2 ARMS2 risk alleles derived maximum benefit from zinc-only supplementation. Patients with one or two CFH risk alleles and no ARMS2 risk alleles derived maximum benefit from antioxidant-only supplementation; treatment with zinc was associated with increased progression to advanced AMD. These recommendations could lead to improved outcomes through genotype-directed therapy.

Meta-analysis of new genome-wide association studies of colorectal cancer risk.

Colorectal cancer is the second leading cause of cancer death in developed countries. Genome-wide association studies (GWAS) have successfully identified novel susceptibility loci for colorectal cancer. To follow up on these findings, and try to identify novel colorectal cancer susceptibility loci, we present results for GWAS of colorectal cancer (2,906 cases, 3,416 controls) that have not previously published main associations. Specifically, we calculated odds ratios and 95% confidence intervals using log-additive models for each study. In order to improve our power to detect novel colorectal cancer susceptibility loci, we performed a meta-analysis combining the results across studies. We selected the most statistically significant single nucleotide polymorphisms (SNPs) for replication using ten independent studies (8,161 cases and 9,101 controls). We again used a meta-analysis to summarize results for the replication studies alone, and for a combined analysis of GWAS and replication studies. We measured ten SNPs previously identified in colorectal cancer susceptibility loci and found eight to be associated with colorectal cancer (p value range 0.02 to 1.8 × 10(-8)). When we excluded studies that have previously published on these SNPs, five SNPs remained significant at p < 0.05 in the combined analysis. No novel susceptibility loci were significant in the replication study after adjustment for multiple testing, and none reached genome-wide significance from a combined analysis of GWAS and replication. We observed marginally significant evidence for a second independent SNP in the BMP2 region at chromosomal location 20p12 (rs4813802; replication p value 0.03; combined p value 7.3 × 10(-5)). In a region on 5p33.15, which includes the coding regions of the TERT-CLPTM1L genes and has been identified in GWAS to be associated with susceptibility to at least seven other cancers, we observed a marginally significant association with rs2853668 (replication p value 0.03; combined p value 1.9 × 10(-4)). Our study suggests a complex nature of the contribution of common genetic variants to risk for colorectal cancer.

The utility and predictive value of combinations of low penetrance genes for screening and risk prediction of colorectal cancer.

Despite the fact that colorectal cancer (CRC) is a highly treatable form of cancer if detected early, a very low proportion of the eligible population undergoes screening for this form of cancer. Integrating a genomic screening profile as a component of existing screening programs for CRC could potentially improve the effectiveness of population screening by allowing the assignment of individuals to different types and intensities of screening and also by potentially increasing the uptake of existing screening programs. We evaluated the utility and predictive value of genomic profiling as applied to CRC, and as a potential component of a population-based cancer screening program. We generated simulated data representing a typical North American population including a variety of genetic profiles, with a range of relative risks and prevalences for individual risk genes. We then used these data to estimate parameters characterizing the predictive value of a logistic regression model built on genetic markers for CRC. Meta-analyses of genetic associations with CRC were used in building science to inform the simulation work, and to select genetic variants to include in logistic regression model-building using data from the ARCTIC study in Ontario, which included 1,200 CRC cases and a similar number of cancer-free population-based controls. Our simulations demonstrate that for reasonable assumptions involving modest relative risks for individual genetic variants, that substantial predictive power can be achieved when risk variants are common (e.g., prevalence > 20%) and data for enough risk variants are available (e.g., approximately 140-160). Pilot work in population data shows modest, but statistically significant predictive utility for a small collection of risk variants, smaller in effect than age and gender alone in predicting an individual's CRC risk. Further genotyping and many more samples will be required, and indeed the discovery of many more risk loci associated with CRC before the question of the potential utility of germline genomic profiling can be definitively answered.