Genetic Predictors of Circulating 25-Hydroxyvitamin D and Prognosis after Colorectal Cancer

BACKGROUND

Low serum 25-hydroxyvitamin D [25(OH)D] concentrations in patients with colorectal cancer have been consistently associated with higher mortality in observational studies. It is unclear whether low 25(OH)D levels directly influence colorectal cancer mortality. To minimize bias, we use genetic variants associated with vitamin D levels to evaluate the association with overall and colorectal cancer-specific survival.

METHODS

Six genetic variants have been robustly identified to be associated with 25(OH)D levels in genome-wide association studies. On the basis of data from the International Survival Analysis in Colorectal Cancer Consortium, the individual genetic variants and a weighted genetic risk score were tested for association with overall and colorectal cancer-specific survival using Cox proportional hazards models in 7,657 patients with stage I to IV colorectal cancer, of whom 2,438 died from any cause and 1,648 died from colorectal cancer.

RESULTS

The 25(OH)D decreasing allele of SNP rs2282679 (GC gene, encodes group-specific component/vitamin D-binding protein) was associated with poorer colorectal cancer-specific survival, although not significant after multiple-testing correction. None of the other five SNPs showed an association. The genetic risk score showed nonsignificant associations with increased overall [HR = 1.54; confidence interval (CI), 0.86-2.78] and colorectal cancer-specific mortality (HR = 1.76; 95% CI, 0.86-3.58). A significant increased risk of overall mortality was observed in women (HR = 3.26; 95% CI, 1.45-7.33; P _{heterogeneity} = 0.01) and normal-weight individuals (HR = 4.14; 95% CI, 1.50-11.43, P _{heterogeneity} = 0.02).

CONCLUSIONS

Our results provided little evidence for an association of genetic predisposition of lower vitamin D levels with increased overall or colorectal cancer-specific survival, although power might have been an issue.

IMPACT

Further studies are warranted to investigate the association in specific subgroups.

DNA methylation profiling to explore colorectal tumor differences according to menopausal hormone therapy use in women

Aim: Use of menopausal hormone therapy (MHT) has been associated with a reduced risk for colorectal cancer, but mechanisms underlying this relationship are not well understood. In the colon, MHT appears to act through estrogen receptor β (ERβ) which may influence DNA methylation by binding to DNA. Using genome-wide methylation profiling data, we aimed to identify genes that may be differentially methylated according to MHT use. Materials & methods: DNA methylation was measured using Illumina HumanMethylation450k arrays in two independent tumor sample sets of colorectal cancer patients. Differential methylation was determined using R/limma. Results: In the discovery analysis, two CpG sites showed differential DNA methylation according to MHT use, both were not replicated. In stratified analyses, 342 CpG sites were associated with current MHT use only in ERβ-positive tumors. Conclusion: The suggestive findings of differential methylation according to current MHT use in ERβ-positive tumors warrant further investigation.

Mendelian randomization analysis of C-reactive protein on colorectal cancer risk

BACKGROUND

Chronic inflammation is a risk factor for colorectal cancer (CRC). Circulating C-reactive protein (CRP) is also moderately associated with CRC risk. However, observational studies are susceptible to unmeasured confounding or reverse causality. Using genetic risk variants as instrumental variables, we investigated the causal relationship between genetically elevated CRP concentration and CRC risk, using a Mendelian randomization approach.

METHODS

Individual-level data from 30 480 CRC cases and 22 844 controls from 33 participating studies in three international consortia were used: the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colorectal Transdisciplinary Study (CORECT) and the Colon Cancer Family Registry (CCFR). As instrumental variables, we included 19 single nucleotide polymorphisms (SNPs) previously associated with CRP concentration. The SNP-CRC associations were estimated using a logistic regression model adjusted for age, sex, principal components and genotyping phases. An inverse-variance weighted method was applied to estimate the causal effect of CRP on CRC risk.

RESULTS

Among the 19 CRP-associated SNPs, rs1260326 and rs6734238 were significantly associated with CRC risk (P = 7.5 × 10-4, and P = 0.003, respectively). A genetically predicted one-unit increase in the log-transformed CRP concentrations (mg/l) was not associated with increased risk of CRC [odds ratio (OR) = 1.04; 95% confidence interval (CI): 0.97, 1.12; P = 0.256). No evidence of association was observed in subgroup analyses stratified by other risk factors.

CONCLUSIONS

In spite of adequate statistical power to detect moderate association, we found genetically elevated CRP concentration was not associated with increased risk of CRC among individuals of European ancestry. Our findings suggested that circulating CRP is unlikely to be a causal factor in CRC development.

Association analyses identify 31 new risk loci for colorectal cancer susceptibility

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

Genome-wide DNA methylation differences according to oestrogen receptor beta status in colorectal cancer

Involvement of sex hormones in colorectal cancer (CRC) development has been linked to oestrogen receptor β (ERβ). Expression of ERβ is found reduced in tumour tissue and inversely related to mortality. However, mechanisms are not well understood. Our study aimed to detect differentially methylated genes associated with ERβ expression, which could point to mechanisms by which ERβ could influence risk and prognosis of CRC. Epigenome-wide DNA methylation profiling was performed using Illumina HumanMethylation450k BeadChip arrays in two independent tumour sample sets of CRC patients recruited in 2003-2010 by the German DACHS study (discovery cohort n = 917, replication cohort n = 907). ERβ expression was measured using immunohistochemistry and scored as negative, moderate and high. Differentially methylated CpG sites and genomic regions were determined using limma in the R-package RnBeads. For the comparison of tumours with moderate/high ERβ versus negative expression, differentially methylated CpG sites were identified but not confirmed by replication. Comparing tumours of high with tumours of negative ERβ expression revealed 2,904 differentially methylated CpG sites of which 403 were replicated (FDR adjusted p-value<0.05). Replicated CpGs were annotated to genes such as CD36, HK1 or LRP5. A survival analysis indicates that 30 of the replicated CpGs are also associated with overall survival (FDR-adjusted p-value<0.05). The regional analysis identified 60 differentially methylated promotor regions. The epigenome-wide analysis identified both novel genes as well as genes already implicated in CRC. Follow-up mechanistic studies to better understand the regulatory role of ERβ could inform potential targets for improving treatment or prevention of CRC.

Novel Common Genetic Susceptibility Loci for Colorectal Cancer

BACKGROUND

Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10-8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk.

METHODS

We conducted a GWAS in European descent CRC cases and control subjects using a discovery-replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10-8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided.

RESULTS

The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10-8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0.

CONCLUSIONS

This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Genetic polymorphisms in genes of class switch recombination and multiple myeloma risk and survival: an IMMEnSE study

Genetic variants in genes acting during the maturation process of immature B-cell to differentiated plasma cell could influence the risk of developing multiple myeloma (MM). During B-cell maturation, several programmed genetic rearrangements occur to increase the variation of the immunoglobulin chains. Class switch recombination (CSR) is one of the most important among these mechanisms. Germline polymorphisms altering even subtly this process could play a role in the etiology and outcome of MM. We performed an association study of 30 genetic variants in the key CSR genes, using 2632 MM patients and 2848 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium, the Heidelberg MM Group and the ESTHER cohort. We found an association between LIG4-rs1555902 and decreased MM risk, which approached statistical significance, as well as significant associations between AICDA-rs3794318 and better outcome. Our results add to our knowledge on the genetic component of MM risk and survival.

Discovery of common and rare genetic risk variants for colorectal cancer

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

Evidence for PTGER4, PSCA, and MBOAT7 as risk genes for gastric cancer on the genome and transcriptome level

Genetic associations between variants on chromosome 5p13 and 8q24 and gastric cancer (GC) have been previously reported in the Asian population. We aimed to replicate these findings and to characterize the associations at the genome and transcriptome level. We performed a fine‐mapping association study in 1926 GC patients and 2012 controls of European descent using high dense SNP marker sets on both chromosomal regions. Next, we performed expression quantitative trait locus (eQTL) analyses using gastric transcriptome data from 143 individuals focusing on the GC associated variants. On chromosome 5p13 the strongest association was observed at rs6872282 (P = 2.53 × 10⁻⁰⁴) and on chromosome 8q24 at rs2585176 (P = 1.09 × 10⁻⁰⁹). On chromosome 5p13 we found cis‐eQTL effects with an upregulation of PTGER4 expression in GC risk allele carrier (P = 9.27 × 10⁻¹¹). On chromosome 8q24 we observed cis‐eQTL effects with an upregulation of PSCA expression in GC risk allele carrier (P = 2.17 × 10⁻⁴⁷). In addition, we found trans‐eQTL effects for the same variants on 8q24 with a downregulation of MBOAT7 expression in GC risk allele carrier (P = 3.11 × 10⁻⁰⁹). In summary, we confirmed and refined the previously reported GC associations at both chromosomal regions. Our data point to shared etiological factors between Asians and Europeans. Furthermore, our data imply an upregulated expression of PTGER4 and PSCA as well as a downregulated expression of MBOAT7 in gastric tissue as risk‐conferring GC pathomechanisms.

Pathway analysis of genetic variants in folate-mediated one-carbon metabolism-related genes and survival in a prospectively followed cohort of colorectal cancer patients

Folate-mediated one-carbon metabolism (FOCM) is a key pathway essential for nucleotide synthesis, DNA methylation, and repair. This pathway is a critical target for 5-fluorouracil (5-FU), which is predominantly used for colorectal cancer (CRC) treatment. A comprehensive assessment of polymorphisms in FOCM-related genes and their association with prognosis has not yet been performed. Within 1,739 CRC cases aged ≥30 years diagnosed from 2003 to 2007 (DACHS study), we investigated 397 single nucleotide polymorphisms (SNPs) and 50 candidates in 48 FOCM-related genes for associations with overall- (OS) and disease-free survival (DFS) using multiple Cox regression (adjusted for age, sex, stage, grade, BMI, and alcohol). We investigated effect modification by 5-FU-based chemotherapy and assessed pathway-specific effects. Correction for multiple testing was performed using false discovery rates (FDR). After a median follow-up time of 5.0 years, 585 patients were deceased. For one candidate SNP in MTHFR and two in TYMS, we observed significant inverse associations with OS (MTHFR: rs1801133, C677T: HRhet  = 0.81, 95% CI: 0.67-0.97; TYMS: rs1001761: HRhet  = 0.82, 95% CI: 0.68-0.99 and rs2847149: HRhet  = 0.82, 95% CI: 0.68-0.99). After FDR correction, one polymorphism in paraoxonase 1 (PON1; rs3917538) was significantly associated with OS (HRhet  = 1.28, 95% CI: 1.07-1.53; HRhzv  = 2.02, 95% CI:1.46-2.80; HRlogAdd  = 1.31, pFDR  = 0.01). Adjusted pathway analyses showed significant associations for pyrimidine biosynthesis (P = 0.04) and fluorouracil drug metabolism (P < 0.01) with significant gene-chemotherapy interactions, including PON1 rs3917538. This study supports the concept that FOCM-related genes could be associated with CRC survival and may modify effects of 5-FU-based chemotherapy in genes in pyrimidine and fluorouracil metabolism, which are relevant targets for therapeutic response and prognosis in CRC. These results require confirmation in additional clinical studies.

Mendelian randomisation study of age at menarche and age at menopause and the risk of colorectal cancer

Background

Substantial evidence supports an association between use of menopausal hormone therapy and decreased colorectal cancer (CRC) risk, indicating a role of exogenous sex hormones in CRC development. However, findings on endogenous oestrogen exposure and CRC are inconsistent.

Methods

We used a Mendelian randomisation approach to test for a causal effect of age at menarche and age at menopause as surrogates for endogenous oestrogen exposure on CRC risk. Weighted genetic risk scores based on 358 single-nucleotide polymorphisms associated with age at menarche and 51 single-nucleotide polymorphisms associated with age at menopause were used to estimate the association with CRC risk using logistic regression in 12,944 women diagnosed with CRC and 10,741 women without CRC from three consortia. Sensitivity analyses were conducted to address pleiotropy and possible confounding by body mass index.

Results

Genetic risk scores for age at menarche (odds ratio per year 0.98, 95% confidence interval: 0.95–1.02) and age at menopause (odds ratio 0.98, 95% confidence interval: 0.94–1.01) were not significantly associated with CRC risk. The sensitivity analyses yielded similar results.

Conclusions

Our study does not support a causal relationship between genetic risk scores for age at menarche and age at menopause and CRC risk.

Do pancreatic cancer and chronic pancreatitis share the same genetic risk factors? A PANcreatic Disease ReseArch (PANDoRA) consortium investigation

Pancreatic ductal adenocarcinoma (PDAC) is a very aggressive tumor with a five-year survival of less than 6%. Chronic pancreatitis (CP), an inflammatory process in of the pancreas, is a strong risk factor for PDAC. Several genetic polymorphisms have been discovered as susceptibility loci for both CP and PDAC. Since CP and PDAC share a consistent number of epidemiologic risk factors, the aim of this study was to investigate whether specific CP risk loci also contribute to PDAC susceptibility. We selected five common SNPs (rs11988997, rs379742, rs10273639, rs2995271 and rs12688220) that were identified as susceptibility markers for CP and analyzed them in 2,914 PDAC cases, 356 CP cases and 5,596 controls retrospectively collected in the context of the international PANDoRA consortium. We found a weak association between the minor allele of the PRSS1-PRSS2-rs10273639 and an increased risk of developing PDAC (ORhomozygous  = 1.19, 95% CI 1.02-1.38, p = 0.023). Additionally all the SNPs confirmed statistically significant associations with risk of developing CP, the strongest being PRSS1-PRSS2-rs10273639 (ORheterozygous  = 0.51, 95% CI 0.39-0.67, p = 1.10 × 10-6 ) and MORC4-rs 12837024 (ORhomozygous  = 2.07 (1.55-2.77, ptrend  = 0.7 × 10-11 ). Taken together, the results from our study do not support variants rs11988997, rs379742, rs10273639, rs2995271 and rs12688220 as strong predictors of PDAC risk, but further support the role of these SNPs in CP susceptibility. Our study suggests that CP and PDAC probably do not share genetic susceptibility, at least in terms of high frequency variants.

Association analysis identifies 65 new breast cancer risk loci

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

Leukocyte Telomere Length and All-Cause, Cardiovascular Disease, and Cancer Mortality: Results From Individual-Participant-Data Meta-Analysis of 2 Large Prospective Cohort Studies

We studied the associations of leukocyte telomere length (LTL) with all-cause, cardiovascular disease, and cancer mortality in 12,199 adults participating in 2 population-based prospective cohort studies from Europe (ESTHER) and the United States (Nurses' Health Study). Blood samples were collected in 1989-1990 (Nurses' Health Study) and 2000-2002 (ESTHER). LTL was measured by quantitative polymerase chain reaction. We calculated z scores for LTL to standardize LTL measurements across the cohorts. Cox proportional hazards regression models were used to calculate relative mortality according to continuous levels and quintiles of LTL z scores. The hazard ratios obtained from each cohort were subsequently pooled by meta-analysis. Overall, 2,882 deaths were recorded during follow-up (Nurses' Health Study, 1989-2010; ESTHER, 2000-2015). LTL was inversely associated with age in both cohorts. After adjustment for age, a significant inverse trend of LTL with all-cause mortality was observed in both cohorts. In random-effects meta-analysis, age-adjusted hazard ratios for the shortest LTL quintile compared with the longest were 1.23 (95% confidence interval (CI): 1.04, 1.46) for all-cause mortality, 1.29 (95% CI: 0.83, 2.00) for cardiovascular mortality, and 1.10 (95% CI: 0.88, 1.37) for cancer mortality. In this study population with an age range of 43-75 years, we corroborated previous evidence suggesting that LTL predicts all-cause mortality beyond its association with age.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Fine-Mapping of Common Genetic Variants Associated with Colorectal Tumor Risk Identified Potential Functional Variants

Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) associated with colorectal cancer risk. These SNPs may tag correlated variants with biological importance. Fine-mapping around GWAS loci can facilitate detection of functional candidates and additional independent risk variants. We analyzed 11,900 cases and 14,311 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. To fine-map genomic regions containing all known common risk variants, we imputed high-density genetic data from the 1000 Genomes Project. We tested single-variant associations with colorectal tumor risk for all variants spanning genomic regions 250-kb upstream or downstream of 31 GWAS-identified SNPs (index SNPs). We queried the University of California, Santa Cruz Genome Browser to examine evidence for biological function. Index SNPs did not show the strongest association signals with colorectal tumor risk in their respective genomic regions. Bioinformatics analysis of SNPs showing smaller P-values in each region revealed 21 functional candidates in 12 loci (5q31.1, 8q24, 11q13.4, 11q23, 12p13.32, 12q24.21, 14q22.2, 15q13, 18q21, 19q13.1, 20p12.3, and 20q13.33). We did not observe evidence of additional independent association signals in GWAS-identified regions. Our results support the utility of integrating data from comprehensive fine-mapping with expanding publicly available genomic databases to help clarify GWAS associations and identify functional candidates that warrant more onerous laboratory follow-up. Such efforts may aid the eventual discovery of disease-causing variant(s).

Epigenome-wide search for association of serum 25-hydroxyvitamin D concentration with leukocyte DNA methylation in a large cohort of older men

Aim: We aimed for an epigenome-wide identification of vitamin D-associated CpG sites in leukocyte DNA. Materials & methods: Infinium HumanMethylation450BeadChip measurements in 402 Caucasian older men were evaluated for significant association with 25-hydroxy-vitamin (25(OH)D) using Spearman's correlation and median regression to adjust for confounding variables. A cross-validation approach as well as a bootstrapping procedure were implemented to determine the replicability of significant associations. Multiple testing was corrected for by Benjamini–Hochberg or Bonferroni. Results: Although in the screening subcohorts significant associations of DNAm with 25(OH)D were observed in the validation cohorts these associations were not replicated after adjustment for potential confounders. At none of the 361,945 CpGs a significant association of DNAm with 25(OH)D was found in all 100 random bootstrap samples, but in comparison at 462 CpGs for the well-established association with age. Conclusion: Leukocyte DNAm was not associated with 25(OH)D levels after validation and consideration of confounders.

Genome-wide Association Studies Identify Genetic Loci Associated With Albuminuria in Diabetes

Elevated concentrations of albumin in the urine, albuminuria, are a hallmark of diabetic kidney disease and are associated with an increased risk for end-stage renal disease and cardiovascular events. To gain insight into the pathophysiological mechanisms underlying albuminuria, we conducted meta-analyses of genome-wide association studies and independent replication in up to 5,825 individuals of European ancestry with diabetes and up to 46,061 without diabetes, followed by functional studies. Known associations of variants in CUBN, encoding cubilin, with the urinary albumin-to-creatinine ratio (UACR) were confirmed in the overall sample (P = 2.4 × 10(-10)). Gene-by-diabetes interactions were detected and confirmed for variants in HS6ST1 and near RAB38/CTSC. Single nucleotide polymorphisms at these loci demonstrated a genetic effect on UACR in individuals with but not without diabetes. The change in the average UACR per minor allele was 21% for HS6ST1 (P = 6.3 × 10(-7)) and 13% for RAB38/CTSC (P = 5.8 × 10(-7)). Experiments using streptozotocin-induced diabetic Rab38 knockout and control rats showed higher urinary albumin concentrations and reduced amounts of megalin and cubilin at the proximal tubule cell surface in Rab38 knockout versus control rats. Relative expression of RAB38 was higher in tubuli of patients with diabetic kidney disease compared with control subjects. The loci identified here confirm known pathways and highlight novel pathways influencing albuminuria.

BRCA2 Polymorphic Stop Codon K3326X and the Risk of Breast, Prostate, and Ovarian Cancers

BACKGROUND

The K3326X variant in BRCA2 (BRCA2*c.9976A>T; p.Lys3326*; rs11571833) has been found to be associated with small increased risks of breast cancer. However, it is not clear to what extent linkage disequilibrium with fully pathogenic mutations might account for this association. There is scant information about the effect of K3326X in other hormone-related cancers.

METHODS

Using weighted logistic regression, we analyzed data from the large iCOGS study including 76 637 cancer case patients and 83 796 control patients to estimate odds ratios (ORw) and 95% confidence intervals (CIs) for K3326X variant carriers in relation to breast, ovarian, and prostate cancer risks, with weights defined as probability of not having a pathogenic BRCA2 variant. Using Cox proportional hazards modeling, we also examined the associations of K3326X with breast and ovarian cancer risks among 7183 BRCA1 variant carriers. All statistical tests were two-sided.

RESULTS

The K3326X variant was associated with breast (ORw = 1.28, 95% CI = 1.17 to 1.40, P = 5.9x10(-) (6)) and invasive ovarian cancer (ORw = 1.26, 95% CI = 1.10 to 1.43, P = 3.8x10(-3)). These associations were stronger for serous ovarian cancer and for estrogen receptor-negative breast cancer (ORw = 1.46, 95% CI = 1.2 to 1.70, P = 3.4x10(-5) and ORw = 1.50, 95% CI = 1.28 to 1.76, P = 4.1x10(-5), respectively). For BRCA1 mutation carriers, there was a statistically significant inverse association of the K3326X variant with risk of ovarian cancer (HR = 0.43, 95% CI = 0.22 to 0.84, P = .013) but no association with breast cancer. No association with prostate cancer was observed.

CONCLUSIONS

Our study provides evidence that the K3326X variant is associated with risk of developing breast and ovarian cancers independent of other pathogenic variants in BRCA2. Further studies are needed to determine the biological mechanism of action responsible for these associations.

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.

Type 2 diabetes and leucocyte DNA methylation: an epigenome-wide association study in over 1,500 older adults

AIMS/HYPOTHESIS

Development of type 2 diabetes depends on environmental and genetic factors. We investigated the epigenome-wide association of prevalent diabetes with DNA methylation (DNAm) in peripheral blood.

METHODS

DNAm was measured in whole blood with the Illumina Infinium HumanMethylation450 BeadChip in two subsamples of participants from the ESTHER cohort study. Cohort 1 included 988 participants, who were consecutively recruited between July and October 2000 and cohort 2 included 527 randomly selected participants. The association of DNAm with prevalent type 2 diabetes at recruitment was estimated using median regression analysis adjusting for sex, age, BMI, smoking behaviour, cell composition and batch at 361,922 CpG sites.

RESULTS

Type 2 diabetes was prevalent in 16% of the participants, and diabetes was poorly controlled in 45% of the diabetic patients. In cohort 1 (discovery) DNAm at 39 CpGs was significantly associated with prevalent diabetes after correction for multiple testing. In cohort 2 (replication) at one of these CpGs, DNAm was still significantly associated. Decreasing methylation levels at cg19693031 with increasing fasting glucose and HbA1c concentrations were observed using restricted cubic spline analysis. In diabetic patients with poorly controlled diabetes, the decrease in estimated DNAm levels was approximately 5% in comparison with participants free of diagnosed diabetes.

CONCLUSIONS/INTERPRETATION

Cg19693031, which is located within the 3'-untranslated region of TXNIP, might play a role in the pathophysiology of type 2 diabetes. This result appears biologically plausible given that thioredoxin-interacting protein is overexpressed in diabetic animals and humans and 3'-untranslated regions are known to play a regulatory role in gene expression.

Smoking-Associated DNA Methylation Biomarkers and Their Predictive Value for All-Cause and Cardiovascular Mortality

BACKGROUND

With epigenome-wide mapping of DNA methylation, a number of novel smoking-associated loci have been identified.

OBJECTIVES

We aimed to assess dose-response relationships of methylation at the top hits from the epigenome-wide methylation studies with smoking exposure as well as with total and cause-specific mortality.

METHODS

In a population-based prospective cohort study in Germany, methylation was quantified in baseline blood DNA of 1,000 older adults by the Illumina 450K assay. Deaths were recorded during a median follow-up of 10.3 years. Dose-response relationships of smoking exposure with methylation at nine CpGs were modeled by restricted cubic spline regression. Associations of individual and aggregate methylation patterns with all-cause, cardiovascular, and cancer mortality were assessed by multiple Cox regression.

RESULTS

Clear dose-response relationships with respect to current and lifetime smoking intensity were consistently observed for methylation at six of the nine CpGs. Seven of the nine CpGs were also associated with mortality outcomes to various extents. A methylation score based on the top two CpGs (cg05575921 and cg06126421) showed the strongest associations with all-cause, cardiovascular, and cancer mortality, with adjusted hazard ratios (95% CI) of 3.59 (2.10, 6.16), 7.41 (2.81, 19.54), and 2.48 (1.01, 6.08), respectively, for participants with methylation levels in the lowest quartile at both CpGs. Adding methylation at those two CpGs into a model that included the variables of the Systematic Coronary Risk Evaluation chart for fatal cardiovascular risk prediction improved the predictive discrimination.

CONCLUSION

The novel methylation biomarkers are highly informative for both smoking exposure and smoking-related mortality outcomes. In particular, these biomarkers may substantially improve cardiovascular risk prediction. Nevertheless, the findings of the present study need to be further validated in additional large longitudinal studies.

CITATION

Zhang Y, Schöttker B, Florath I, Stock C, Butterbach K, Holleczek B, Mons U, Brenner H. 2016. Smoking-associated DNA methylation biomarkers and their predictive value for all-cause and cardiovascular mortality. Environ Health Perspect 124:67-74; http://dx.doi.org/10.1289/ehp.1409020.

TERT gene harbors multiple variants associated with pancreatic cancer susceptibility

A small number of common susceptibility loci have been identified for pancreatic cancer, one of which is marked by rs401681 in the TERT-CLPTM1L gene region on chromosome 5p15.33. Because this region is characterized by low linkage disequilibrium, we sought to identify whether additional single nucleotide polymorphisms (SNPs) could be related to pancreatic cancer risk, independently of rs401681. We performed an in-depth analysis of genetic variability of the telomerase reverse transcriptase (TERT) and the telomerase RNA component (TERC) genes, in 5,550 subjects with pancreatic cancer and 7,585 controls from the PANcreatic Disease ReseArch (PANDoRA) and the PanScan consortia. We identified a significant association between a variant in TERT and pancreatic cancer risk (rs2853677, odds ratio = 0.85; 95% confidence interval = 0.80-0.90, p = 8.3 × 10(-8)). Additional analysis adjusting rs2853677 for rs401681 indicated that the two SNPs are independently associated with pancreatic cancer risk, as suggested by the low linkage disequilibrium between them (r(2) = 0.07, D' = 0.28). Three additional SNPs in TERT reached statistical significance after correction for multiple testing: rs2736100 (p = 3.0 × 10(-5) ), rs4583925 (p = 4.0 × 10(-5) ) and rs2735948 (p = 5.0 × 10(-5) ). In conclusion, we confirmed that the TERT locus is associated with pancreatic cancer risk, possibly through several independent variants.

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same region.

Prediction of individual genetic risk to prostate cancer using a polygenic score

BACKGROUND

Polygenic risk scores comprising established susceptibility variants have shown to be informative classifiers for several complex diseases including prostate cancer. For prostate cancer it is unknown if inclusion of genetic markers that have so far not been associated with prostate cancer risk at a genome-wide significant level will improve disease prediction.

METHODS

We built polygenic risk scores in a large training set comprising over 25,000 individuals. Initially 65 established prostate cancer susceptibility variants were selected. After LD pruning additional variants were prioritized based on their association with prostate cancer. Six-fold cross validation was performed to assess genetic risk scores and optimize the number of additional variants to be included. The final model was evaluated in an independent study population including 1,370 cases and 1,239 controls.

RESULTS

The polygenic risk score with 65 established susceptibility variants provided an area under the curve (AUC) of 0.67. Adding an additional 68 novel variants significantly increased the AUC to 0.68 (P = 0.0012) and the net reclassification index with 0.21 (P = 8.5E-08). All novel variants were located in genomic regions established as associated with prostate cancer risk.

CONCLUSIONS

Inclusion of additional genetic variants from established prostate cancer susceptibility regions improves disease prediction.

Smoking habits and leukocyte telomere length dynamics among older adults: Results from the ESTHER cohort

BACKGROUND & AIMS

Leukocyte telomere length (LTL) shortens with age and short LTL has been associated with increased mortality and increased risk for some age-related outcomes. This study aims to analyse the associations of smoking habits with LTL and rate of LTL change per year in older adults.

METHODS

LTL was measured by quantitative PCR at baseline in 3600 older adults, who were enrolled in a population-based cohort study in Germany. For longitudinal analyses, measurements were repeated in blood samples obtained at 8-year follow-up from 1000 participants. Terminal Restriction Fragment analysis was additionally performed in a sub-sample to obtain absolute LTL in base pairs. Multivariate linear regression models were used to estimate associations of smoking habits with baseline LTL and changes in LTL over time.

RESULTS

LTL was inversely associated with age (r=-0.090, p<0.0001). Women had longer LTL than men (p<0.0001). Smoking was inversely associated with LTL. On average, current smokers had 73 base pairs (BP) shorter LTL compared to never smokers. Smoking intensity and pack-years of smoking were also inversely associated with LTL, and a positive association was observed with years since smoking cessation. Slower LTL attrition rates were observed in ever smokers over 8years of follow-up.

CONCLUSIONS

Our cross-sectional analysis supports suggestions that smoking might contribute to shortening of LTL but this relationship could not be shown longitudinally. The overall rather small effect sizes observed for smoking-related variables suggest that LTL reflects smoking-related health hazards only to a very limited extent.

Discovery of a novel epigenetic cancer marker related to the oxidative status of human blood

Long-lasting oxidative stress exposure may lead to relatively stable epigenetic modifications of the DNA in order to activate anti-oxidative defence mechanisms. Oxidative stress related DNA methylation may therefore be associated (causally or as a by-product) with cancer. We measured derivatives of reactive oxygen metabolites (D-ROM), total thiol levels (TTL) and DNA methylation with the Illumina Infinium 450K BeadChip in three samples of German individuals aged ≥50 years: n = 1,000 ESTHER study baseline participants (DNA methylation only), n = 99 ESTHER eight-year follow-up participants and n = 142 participants of the BLITZ study. The correlation coefficient of methylation at cg10342304 and D-ROM in the ESTHER 8-year follow-up sample (r = -0.427; P = 1 × 10(-5)) was replicated with a P-value indicating statistical significance after correction for multiple testing in the BLITZ sample (r = -0.192; P = 0.022). The association was robust to adjusting for potential confounders. In the ESTHER baseline sample, the hazard ratio for cancer development in 11 years of follow-up comparing bottom and top quartile of DNA methylation at cg10342304 was 1.86 (95%-confidence-interval 1.01-3.43). In summary, this first epigenome-wide screening and replication study with oxidative status markers observed a negative correlation of D-ROM levels and DNA methylation at cg10342304 in two independent cohorts. This CpG site is located in the body region of the nucleoredoxin gene. The nucleoredoxin protein is a redox-dependent inhibitor of the Wnt/ß-catenin signaling pathway, a well-characterized cancer pathway. If the observed CpG-cancer association can be successfully replicated by other studies, this epigenetic marker could be an interesting biomarker of cancer risk.

Annexin A1 expression in a pooled breast cancer series: association with tumor subtypes and prognosis

BACKGROUND

Annexin A1 (ANXA1) is a protein related with the carcinogenesis process and metastasis formation in many tumors. However, little is known about the prognostic value of ANXA1 in breast cancer. The purpose of this study is to evaluate the association between ANXA1 expression, BRCA1/2 germline carriership, specific tumor subtypes and survival in breast cancer patients.

METHODS

Clinical-pathological information and follow-up data were collected from nine breast cancer studies from the Breast Cancer Association Consortium (BCAC) (n = 5,752) and from one study of familial breast cancer patients with BRCA1/2 mutations (n = 107). ANXA1 expression was scored based on the percentage of immunohistochemical staining in tumor cells. Survival analyses were performed using a multivariable Cox model.

RESULTS

The frequency of ANXA1 positive tumors was higher in familial breast cancer patients with BRCA1/2 mutations than in BCAC patients, with 48.6 % versus 12.4 %, respectively; P <0.0001. ANXA1 was also highly expressed in BCAC tumors that were poorly differentiated, triple negative, EGFR-CK5/6 positive or had developed in patients at a young age. In the first 5 years of follow-up, patients with ANXA1 positive tumors had a worse breast cancer-specific survival (BCSS) than ANXA1 negative (HRadj = 1.35; 95 % CI = 1.05-1.73), but the association weakened after 10 years (HRadj = 1.13; 95 % CI = 0.91-1.40). ANXA1 was a significant independent predictor of survival in HER2+ patients (10-years BCSS: HRadj = 1.70; 95 % CI = 1.17-2.45).

CONCLUSIONS

ANXA1 is overexpressed in familial breast cancer patients with BRCA1/2 mutations and correlated with poor prognosis features: triple negative and poorly differentiated tumors. ANXA1 might be a biomarker candidate for breast cancer survival prediction in high risk groups such as HER2+ cases.

Somatic MyD88 mutations and a germline IRAK2 variant contribute to lymphoma by differential effects on NF-κB (INM7P.354)

Within innate immune signaling pathways, Myeloid differentiation (MyD) 88 and Interleukin-1 receptor associated kinases (IRAKs) fulfill pivotal roles downstream of multiple Toll-like receptors (TLR). Somatic MyD88 mutations in B cells were reported to affect lymphomagenesis by constitutive hyperactivation of NF-κB. We here show that oncogenic mutants augment spontaneous TIR domain oligomerization into so-called Myddosome signaling complexes driving NF-κB activation. Blocking of MyD88 oligomerization consequently induced the death of MyD88-mutated but not unmutated lymphoma cells. In line with the involvement of the MyD88 axis in lymphomagenesis we detected a genetic association of a coding germline variant of IRAK2, a MyD88 interactor, with increased risk (OR = 1.43, 95% confidence interval 1.06-1.93, p=0.018) to develop lymphoma in a German cohort (n=387). This variant is frequently found in 15-44% of individuals of different ethnic groups. Unexpectedly, this variant was associated with reduced NF-κB activation and cytokine responses to TLR agonists in healthy donors. Mechanistically, the IRAK2 variant retained intact binding to but impaired ubiquitination of TRAF6, a vital step in signal transduction. Thus our findings confirm the involvement of MyD88/IRAK in lymphoma but suggest that MyD88 and IRAK2 may contribute differently to oncogenic signaling. Additionally we highlight IRAK2 genetic variants as a critical factor and potentially novel biomarker for human lymphomagenesis.

Common germline polymorphisms associated with breast cancer-specific survival

INTRODUCTION

Previous studies have identified common germline variants nominally associated with breast cancer survival. These associations have not been widely replicated in further studies. The purpose of this study was to evaluate the association of previously reported SNPs with breast cancer-specific survival using data from a pooled analysis of eight breast cancer survival genome-wide association studies (GWAS) from the Breast Cancer Association Consortium.

METHODS

A literature review was conducted of all previously published associations between common germline variants and three survival outcomes: breast cancer-specific survival, overall survival and disease-free survival. All associations that reached the nominal significance level of P value <0.05 were included. Single nucleotide polymorphisms that had been previously reported as nominally associated with at least one survival outcome were evaluated in the pooled analysis of over 37,000 breast cancer cases for association with breast cancer-specific survival. Previous associations were evaluated using a one-sided test based on the reported direction of effect.

RESULTS

Fifty-six variants from 45 previous publications were evaluated in the meta-analysis. Fifty-four of these were evaluated in the full set of 37,954 breast cancer cases with 2,900 events and the two additional variants were evaluated in a reduced sample size of 30,000 samples in order to ensure independence from the previously published studies. Five variants reached nominal significance (P <0.05) in the pooled GWAS data compared to 2.8 expected under the null hypothesis. Seven additional variants were associated (P <0.05) with ER-positive disease.

CONCLUSIONS

Although no variants reached genome-wide significance (P <5 x 10(-8)), these results suggest that there is some evidence of association between candidate common germline variants and breast cancer prognosis. Larger studies from multinational collaborations are necessary to increase the power to detect associations, between common variants and prognosis, at more stringent significance levels.

F2RL3 methylation, lung cancer incidence and mortality

Smoking accounts for a large share of lung cancer. F2RL3 methylation was recently identified as a biomarker closely reflecting both current and past smoking exposure. We aimed to assess the associations of F2RL3 methylation with lung cancer incidence and mortality. In a large population-based cohort study, F2RL3 methylation was measured in baseline blood samples of 4,987 participants by MassARRAY. Associations of F2RL3 methylation and smoking with lung cancer incidence/mortality during a median follow-up of 10.9 years were assessed by Cox regression, controlling for potential confounders. The ability of F2RL3 methylation to predict lung cancer was examined by Harrell's C statistics. Hypomethylation at F2RL3 was strongly associated with both lung cancer incidence and mortality, with age- and sex-adjusted hazard ratios (HR; 95% CI) of 9.99 (5.61-17.79) and 16.86 (8.53-33.34), respectively, for participants whose methylation intensity were ≤0.54 compared with whose methylation intensity were ≥0.75. Strongly elevated HRs of 2.88 (1.42-5.84) and 5.17 (2.28-11.70) persisted even after controlling for multiple covariates including smoking status and pack-years. With fully adjusted HRs of 9.92 (2.88-34.12) and 16.48 (4.10-66.15), the associations between methylation and the two outcomes were particularly strong among participants≥65 years. Combination of F2RL3 methylation and pack-years predicted lung cancer incidence with high accuracy (optimism-corrected Harrell's C statistics = 0.86 for participants≥65 years). These findings suggested that F2RL3 methylation is a very strong predictor of lung cancer risk and mortality, particularly at older age. The potential implications of F2RL3 methylation for early detection, risk stratification and prevention of lung cancer warrant further exploration.

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.

Abstract 2188: Genetic variation in angiogenesis-related genes is associated with colorectal cancer risk and prognosis

Introduction

Angiogenesis, the generation of new blood vessels, is crucial in tumor growth, progression, and metastasis. The process involves a variety of factors including signaling, adhesion and chemotactic molecules, ECM (extracellular matrix) proteins, proteinases, transcription factors, growth factors and receptors. Variation in genes encoding these factors potentially influence angiogenic processes. We thus evaluated the association between variants in angiogenesis-related genes and risk as well as prognosis of colorectal cancer (CRC).

Approach

In a candidate pathway approach, we investigated 437 variants in 36 angiogenesis-related genes for association with CRC risk and prognosis in ∼1800 patients and ∼1800 controls of the German DACHS/IMPACT study. Patients were aged 30 years or older and diagnosed between 2003 and 2007. CRC risk was estimated using conditional logistic regression based on the co-dominant inheritance model. In addition, we investigated the association between polymorphisms and overall survival using multivariable Cox regression. Correction for multiple testing was performed using false discovery rates (FDR).

Results

Risk

Several of the investigated variants in ANGTP1, ETS1, FLT4, MMP2, NOTCH4, PDGFRB, TGFB2 were associated with risk of CRC, and one variant in DLL1 (delta-like 1) remained significant after accounting for multiple testing. DLL1 rs9348307 was associated with decreased risk of CRC (ORGC 0.81, 95% CI 0.68-0.96; ORCC 0.62, 95% CI 0.36-1.06, p=0.02).

Survival

After a median follow up time of 5 years, variants in twelve genes (ANGPT1, EFNB2, ETS1, FLT4, JAG1, KDR, MMP2, MMP9, NRP1, NRP2, PDGFRB, TGFB2) were associated with overall survival in CRC patients. When accounting for multiple testing, variants in EFNB2, JAG1 and MMP2 remained significant.

Most of the EFNB2 variants were associated with poorer survival, while one variant in EFNB2 was associated with better survival (rs2391333: HRCT 0.88, 95% CI 0.74-1.05; HRTT 0.60, 95% CI 0.45-0.78, p&lt;0.001). Furthermore, two variants in JAG1 were associated with worse survival in CRC patients (rs3748480: HRTG 1.27, 95% CI 1.05-1.53; HRGG 1.79, 95% CI 1.02-3.13, p=0.02; rs6040062: HRGA 1.26, 95% CI 1.04-1.52; HRAA 1.78, 95% CI 1.02-3.12, p=0.02). Finally, three variants in MMP2 (rs11639960, rs1561219, rs17301608) were associated with poorer survival, while one MMP2 variant was associated with improved survival in CRC patients (rs243847: HRTC 0.87, 95% CI 0.73-1.04, HRCC 0.60, 95% CI 0.46-0.78, p&lt;0.001).

Conclusion

Variants in angiogenesis-related genes were associated with CRC risk and overall survival in CRC patients. Two of the identified genes (DLL1 and JAG1) are involved in NOTCH signaling, which is crucial to cell differentiation, proliferation, apoptosis and angiogenesis, underlining the importance of this oncogenic pathway in colorectal carcinogenesis.

Citation Format: Dominique Scherer, Yesilda Balavarca, Nina Habermann, Katharina Buck, Petra Seibold, Lisanne Kap, Katja Butterbach, Katrin Pfütze, Axel Benner, Michael Hoffmeister, Hermann Brenner, Barbara Burwinkel, Jenny Chang-Claude, Cornelia M. Ulrich. Genetic variation in angiogenesis-related genes is associated with colorectal cancer risk and prognosis. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 2188. doi:10.1158/1538-7445.AM2014-2188

Risk of multiple myeloma is associated with polymorphisms within telomerase genes and telomere length

Compelling biological and epidemiological evidences point to a key role of genetic variants of the TERT and TERC genes in cancer development. We analyzed the genetic variability of these two gene regions using samples of 2,267 multiple myeloma (MM) cases and 2,796 healthy controls. We found that a TERT variant, rs2242652, is associated with reduced MM susceptibility (OR = 0.81; 95% CI: 0.72-0.92; p = 0.001). In addition we measured the leukocyte telomere length (LTL) in a subgroup of 140 cases who were chemotherapy-free at the time of blood donation and 468 controls, and found that MM patients had longer telomeres compared to controls (OR = 1.19; 95% CI: 0.63-2.24; p(trend)  = 0.01 comparing the quartile with the longest LTL versus the shortest LTL). Our data suggest the hypothesis of decreased disease risk by genetic variants that reduce the efficiency of the telomerase complex. This reduced efficiency leads to shorter telomere ends, which in turn may also be a marker of decreased MM risk.

GWAS-identified common variants for obesity are not associated with the risk of developing colorectal cancer

BACKGROUND

Observational studies have consistently associated obesity with colorectal cancer risk. Because both traits are genetically determined and share some metabolic biomarkers, we hypothesized that obesity-related polymorphisms could also influence the risk of developing colorectal cancer.

METHODS

We conducted a comprehensive population-based case-control study in 1,792 German colorectal cancer cases and 1,805 controls to explore associations between 28 obesogenic variants identified through genome-wide association studies (GWAS) and colorectal cancer risk. We also evaluated interactions between polymorphisms and body mass index (BMI), type II diabetes (T2D), and gender.

RESULTS

No evidence of association between obesogenic variants and colorectal cancer risk was observed after correction for multiple testing. There was only a remarkable interaction between the LTArs1041981 polymorphism and gender, which modified the risk of colorectal cancer [Pinteraction = 0.002; males: odds ratio (OR), 1.14; 95% confidence intervals (CI), 1.00-1.30 vs. females: OR, 0.83; 95% CI, 0.71-0.97].

CONCLUSIONS

Our findings showed that obesogenic variants are not a major pathogenetic risk factor for colorectal cancer.

IMPACT

This comprehensive population-based case-control study does not provide evidence of a shared genetic component between obesity and colorectal cancer.

Estimating the heritability of colorectal cancer

A sizable fraction of colorectal cancer (CRC) is expected to be explained by heritable factors, with heritability estimates ranging from 12 to 35% twin and family studies. Genome-wide association studies (GWAS) have successfully identified a number of common single-nucleotide polymorphisms (SNPs) associated with CRC risk. Although it has been shown that these CRC susceptibility SNPs only explain a small proportion of the genetic risk, it is not clear how much of the heritability these SNPs explain and how much is left to be detected by other, yet to be identified, common SNPs. Therefore, we estimated the heritability of CRC under different scenarios using Genome-Wide Complex Trait Analysis in the Genetics and Epidemiology of Colorectal Cancer Consortium including 8025 cases and 10 814 controls. We estimated that the heritability explained by known common CRC SNPs identified in GWAS was 0.65% (95% CI:0.3-1%; P = 1.11 × 10-16), whereas the heritability explained by all common SNPs was at least 7.42% (95% CI: 4.71-10.12%; P = 8.13 × 10(-8)), suggesting that many common variants associated with CRC risk remain to be detected. Comparing the heritability explained by the common variants with that from twin and family studies, a fraction of the heritability may be explained by other genetic variants, such as rare variants. In addition, our analysis showed that the gene × smoking interaction explained a significant proportion of the CRC variance (P = 1.26 × 10(-2)). In summary, our results suggest that known CRC SNPs only explain a small proportion of the heritability and more common SNPs have yet to be identified.

Cross-sectional and longitudinal changes in DNA methylation with age: an epigenome-wide analysis revealing over 60 novel age-associated CpG sites

Understanding the role of epigenetic modifications, e.g. DNA methylation, in the process of aging requires the characterization of methylation patterns in large cohorts. We analysed >480 000 CpG sites using Infinium HumanMethylation450 BeadChip (Illumina) in whole blood DNA of 965 participants of a population-based cohort study aged between 50 and 75 years. In an exploratory analysis in 400 individuals, 200 CpG sites with the highest Spearman correlation coefficients for the association between methylation and age were identified. Of these 200 CpGs, 162 were significantly associated with age, which was verified in an independent cohort of 498 individuals using mixed linear regression models adjusted for gender, smoking behaviour, age-related diseases and random batch effect and corrected for multiple testing by Bonferroni. In another independent cohort of 67 individuals without history of major age-related diseases and with a follow-up of 8 years, we observed a gain in methylation at 96% (52%, significant) of the positively age-associated CpGs and a loss at all (89%, significant) of the negatively age-associated CpGs in each individual while getting 8 years older. A regression model for age prediction based on 17 CpGs as predicting variables explained 71% of the variance in age with an average accuracy of 2.6 years. In comparison with cord blood samples obtained from the Ulm Birth Cohort Study, we observed a more than 2-fold change in mean methylation levels from birth to older age at 86 CpGs. We were able to identify 65 novel CpG sites with significant association of methylation with age.

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.

Vitamin D receptor genotype rs731236 (Taq1) and breast cancer prognosis

Several studies have suggested that the anticancerogenous effects of vitamin D might be modulated by genetic variants in the vitamin D receptor (VDR) gene. The association of VDR polymorphisms with breast cancer-specific and all-cause mortality after a breast cancer diagnosis remains, however, largely unexplored. We assessed the association of genetic variants in VDR (rs731236, rs1989969, rs2228570, and 11568820) with breast cancer survival in a sample of 498 patients with breast cancer with a mean age at diagnosis of 61 years from Saarland, Germany, who were followed for up to 5 years with respect to total and breast cancer-specific mortality (56 and 48 events, respectively). Adjusted HRs with 95% confidence intervals (CI) were estimated by Cox regression models. We found that patients with breast cancer homozygous for the rare allele of rs731236 (15% of the women in our cohort) had a tendency toward an increased risk for breast cancer-specific mortality. The HR (95% CI) adjusted for age and breast cancer stage was 2.8 (1.1-7.2) for breast cancer-specific mortality and 2.1 (0.9-4.9) for total mortality. Additional adjustment for family history of breast cancer, radical mastectomy, and body mass index only marginally changed the estimates. No association was found for rs1989969, rs2228570, and rs11568820. Our analysis suggests that VDR polymorphism rs731236 might be associated with breast cancer-specific mortality, and if our findings are confirmed in future bigger studies rs731236 might deserve consideration as a prognostic factor in clinical care of patients with breast cancer.

A comprehensive study of polymorphisms in the ABCB1, ABCC2, ABCG2, NR1I2 genes and lymphoma risk

Owing to their role in controlling the efflux of toxic compounds, transporters are central players in the process of detoxification and elimination of xenobiotics, which in turn is related to cancer risk. Among these transporters, ATP-binding cassette B1/multidrug resistance 1 (ABCB1/MDR1), ABCC2/multidrug resistance protein 2 (MRP2) and ABCG2/breast cancer resistance protein (BCRP) affect susceptibility to many hematopoietic malignancies. The maintenance of regulated expression of these transporters is governed through the activation of intracellular "xenosensors" like the nuclear receptor 1I2/pregnane X receptor (NR1I2/PXR). SNPs in genes encoding these regulators have also been implicated in the risk of several cancers. Using a tagging approach, we tested the hypothesis that common polymorphisms in the transporter genes ABCB1, ABCC2, ABCG2 and the regulator gene NR1I2 could be implicated in lymphoma risk. We selected 68 SNPs in the four genes, and we genotyped them in 1,481 lymphoma cases and 1,491 controls of the European case-control study (EpiLymph) using the Illumina GoldenGate™ assay technology. Carriers of the SNP rs6857600 minor allele in ABCG2 was associated with a decrease in risk of B-cell lymphoma (B-NHL) overall (p < 0.001). Furthermore, a decreased risk of chronic lymphocytic leukemia (CLL) was associated with the ABCG2 rs2231142 variant (p = 0.0004), which could be replicated in an independent population. These results suggest a role for this gene in B-NHL susceptibility, especially for CLL.

Association of JAK-STAT pathway related genes with lymphoma risk: results of a European case-control study (EpiLymph)

Previous studies have suggested an important role for the Janus kinase-signal transducer and activator of transcription (JAK-STAT) signalling pathway in tumour development. Therefore, we explored genetic variants in JAK-STAT pathway associated genes with lymphoma risk. In samples of the EpiLymph case-control study we genotyped 1536 single nucleotide polymorphisms (SNPs) using GoldenGate BeadArray™ Technology (Illumina, San Diego, CA, USA). Here, we report the associations between selected SNPs and haplotypes of the JAK-STAT pathway and risk of Hodgkin lymphoma (HL), B-cell non-Hodgkin lymphoma (B-NHL) and most frequent B-NHL subtypes. Among 210 relevant JAK-STAT pathway-related SNPs, polymorphisms in nine genes (BMF, IFNG, IL12A, SOCS1, STAT1, STAT3, STAT5A, STAT6, TP63) were significantly associated with lymphoma risk. At a study-wise significance level, we obtained a risk reduction of 28% among carriers of the heterozygous genotype of the STAT3 variant (rs1053023) for B-NHL. For six other variants within the STAT3 gene we observed an inverse association with different lymphoma subtypes. A reduced risk for HL was observed for the heterozygous genotype of the STAT6 SNP (rs324011). This is an explorative investigation to examine associations between JAK-STAT signalling related genes and lymphoma risk. The results implicate a relevant role of certain pathway-related genes in lymphomagenesis, but still need to be approved by independent studies.

Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32

To identify susceptibility loci for non-Hodgkin lymphoma (NHL) subtypes, we conducted a three-stage genome-wide association study. We identified two variants associated with follicular lymphoma (FL) in 1,465 FL cases/6,958 controls at 6p21.32 (rs10484561, rs7755224, r²=1.0; combined p-values=1.12×10^-29, 2.00×10^-19), providing further support that MHC genetic variation influences FL susceptibility. Confirmatory evidence of a previously reported association was also found between chronic lymphocytic leukemia/small lymphocytic lymphoma and rs735665 (combined p-value=4.24×10^-9).