Genome-Wide Interaction Analyses between Genetic Variants and Alcohol Consumption and Smoking for Risk of Colorectal Cancer.

Genome-wide association studies (GWAS) have identified many genetic susceptibility loci for colorectal cancer (CRC). However, variants in these loci explain only a small proportion of familial aggregation, and there are likely additional variants that are associated with CRC susceptibility. Genome-wide studies of gene-environment interactions may identify variants that are not detected in GWAS of marginal gene effects. To study this, we conducted a genome-wide analysis for interaction between genetic variants and alcohol consumption and cigarette smoking using data from the Colon Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Interactions were tested using logistic regression. We identified interaction between CRC risk and alcohol consumption and variants in the 9q22.32/HIATL1 (Pinteraction = 1.76×10-8; permuted p-value 3.51x10-8) region. Compared to non-/occasional drinking light to moderate alcohol consumption was associated with a lower risk of colorectal cancer among individuals with rs9409565 CT genotype (OR, 0.82 [95% CI, 0.74-0.91]; P = 2.1×10-4) and TT genotypes (OR,0.62 [95% CI, 0.51-0.75]; P = 1.3×10-6) but not associated among those with the CC genotype (p = 0.059). No genome-wide statistically significant interactions were observed for smoking. If replicated our suggestive finding of a genome-wide significant interaction between genetic variants and alcohol consumption might contribute to understanding colorectal cancer etiology and identifying subpopulations with differential susceptibility to the effect of alcohol on CRC risk.

Etiology of hormone receptor positive breast cancer differs by levels of histologic grade and proliferation.

Limited epidemiological evidence suggests that the etiology of hormone receptor positive (HR+) breast cancer may differ by levels of histologic grade and proliferation. We pooled risk factor and pathology data on 5,905 HR+ breast cancer cases and 26,281 controls from 11 epidemiological studies. Proliferation was determined by centralized automated measures of KI67 in tissue microarrays. Odds ratios (OR), 95% confidence intervals (CI) and p-values for case-case and case-control comparisons for risk factors in relation to levels of grade and quartiles (Q1-Q4) of KI67 were estimated using polytomous logistic regression models. Case-case comparisons showed associations between nulliparity and high KI67 [OR (95% CI) for Q4 vs. Q1 = 1.54 (1.22, 1.95)]; obesity and high grade [grade 3 vs. 1 = 1.68 (1.31, 2.16)] and current use of combined hormone therapy (HT) and low grade [grade 3 vs. 1 = 0.27 (0.16, 0.44)] tumors. In case-control comparisons, nulliparity was associated with elevated risk of tumors with high but not low levels of proliferation [1.43 (1.14, 1.81) for KI67 Q4 vs. 0.83 (0.60, 1.14) for KI67 Q1]; obesity among women ≥50 years with high but not low grade tumors [1.55 (1.17, 2.06) for grade 3 vs. 0.88 (0.66, 1.16) for grade 1] and HT with low but not high grade tumors [3.07 (2.22, 4.23) for grade 1 vs. 0.85 (0.55, 1.30) for grade 3]. Menarcheal age and family history were similarly associated with HR+ tumors of different grade or KI67 levels. These findings provide insights into the etiologic heterogeneity of HR+ tumors.

Combined genetic and splicing analysis of BRCA1 c.[594-2A>C; 641A>G] highlights the relevance of naturally occurring in-frame transcripts for developing disease gene variant classification algorithms.

A recent analysis using family history weighting and co-observation classification modeling indicated that BRCA1 c.594-2A > C (IVS9-2A > C), previously described to cause exon 10 skipping (a truncating alteration), displays characteristics inconsistent with those of a high risk pathogenic BRCA1 variant. We used large-scale genetic and clinical resources from the ENIGMA, CIMBA and BCAC consortia to assess pathogenicity of c.594-2A > C. The combined odds for causality considering case-control, segregation and breast tumor pathology information was 3.23 × 10-8 Our data indicate that c.594-2A > C is always in cis with c.641A > G. The spliceogenic effect of c.[594-2A > C;641A > G] was characterized using RNA analysis of human samples and splicing minigenes. As expected, c.[594-2A > C; 641A > G] caused exon 10 skipping, albeit not due to c.594-2A > C impairing the acceptor site but rather by c.641A > G modifying exon 10 splicing regulatory element(s). Multiple blood-based RNA assays indicated that the variant allele did not produce detectable levels of full-length transcripts, with a per allele BRCA1 expression profile composed of ≈70-80% truncating transcripts, and ≈20-30% of in-frame Δ9,10 transcripts predicted to encode a BRCA1 protein with tumor suppression function.We confirm that BRCA1c.[594-2A > C;641A > G] should not be considered a high-risk pathogenic variant. Importantly, results from our detailed mRNA analysis suggest that BRCA-associated cancer risk is likely not markedly increased for individuals who carry a truncating variant in BRCA1 exons 9 or 10, or any other BRCA1 allele that permits 20-30% of tumor suppressor function. More generally, our findings highlight the importance of assessing naturally occurring alternative splicing for clinical evaluation of variants in disease-causing genes.

Prognostic impact of surgery for early-stage invasive breast cancer on breast cancer-specific survival, overall survival, and recurrence risk: a population-based analysis.

PURPOSE: Recent cohort studies demonstrated better overall survival (OS) or breast cancer-specific survival (BCS) for breast-conserving therapy (BCT) followed by radiation (RT) compared to mastectomy alone (MT). This is the first observational study in which adjustments for a comprehensive set of prognostic factors, adjuvant therapies, mode of detection, and comorbidities were possible to investigate OS, BCS, as well as recurrence risk of patients undergoing BCT + RT, MT + RT, or MT. METHODS: Women aged 50-74 years at diagnosis of early-stage invasive breast cancer (I-IIIa) between 2001 and 2005 at the German population-based case-control study (MARIE study) were recruited and followed prospectively as a case cohort until 2015. Kaplan-Meier estimates and stepwise adjusted multivariable Cox models were used to estimate hazard ratios (HRs) and 95% confidence intervals (95% CI). RESULTS: The 2762 patients included were followed up for a median of 11.9 years (95% CI 11.8-12.0). 74.2% of patients underwent BCT + RT; 10.3% MT + RT and 15.6% MT alone. Compared to patients treated with MT alone, patients treated with BCT + RT showed non-statistically significant improved OS (HR 0.79, 95% CI 0.61-1.02), BCS (HR 0.79, 95% CI 0.55-1.12), and no difference in recurrence risks (HR 1.01, 95% CI 0.74-1.37). For patients treated with MT + RT, there were no differences in OS (HR 1.06, 95% CI 0.75-1.50), BCS (HR 1.17, 95% CI 0.75-1.82), or recurrence risk (HR 1.33, 95% CI 0.89-1.97). CONCLUSIONS: Among patients with early-stage breast cancer, clinical outcomes more than 10 years after diagnosis did not differ between the primary treatment options BCT + RT, MT + RT versus MT alone after full adjustment.

Reproductive profiles and risk of breast cancer subtypes: a multi-center case-only study.

BACKGROUND: Previous studies have shown that reproductive factors are differentially associated with breast cancer (BC) risk by subtypes. The aim of this study was to investigate associations between reproductive factors and BC subtypes, and whether these vary by age at diagnosis. METHODS: We used pooled data on tumor markers (estrogen and progesterone receptor, human epidermal growth factor receptor-2 (HER2)) and reproductive risk factors (parity, age at first full-time pregnancy (FFTP) and age at menarche) from 28,095 patients with invasive BC from 34 studies participating in the Breast Cancer Association Consortium (BCAC). In a case-only analysis, we used logistic regression to assess associations between reproductive factors and BC subtype compared to luminal A tumors as a reference. The interaction between age and parity in BC subtype risk was also tested, across all ages and, because age was modeled non-linearly, specifically at ages 35, 55 and 75 years. RESULTS: Parous women were more likely to be diagnosed with triple negative BC (TNBC) than with luminal A BC, irrespective of age (OR for parity = 1.38, 95% CI 1.16-1.65, p = 0.0004; p for interaction with age = 0.076). Parous women were also more likely to be diagnosed with luminal and non-luminal HER2-like BCs and this effect was slightly more pronounced at an early age (p for interaction with age = 0.037 and 0.030, respectively). For instance, women diagnosed at age 35 were 1.48 (CI 1.01-2.16) more likely to have luminal HER2-like BC than luminal A BC, while this association was not significant at age 75 (OR = 0.72, CI 0.45-1.14). While age at menarche was not significantly associated with BC subtype, increasing age at FFTP was non-linearly associated with TNBC relative to luminal A BC. An age at FFTP of 25 versus 20 years lowered the risk for TNBC (OR = 0.78, CI 0.70-0.88, p < 0.0001), but this effect was not apparent at a later FFTP. CONCLUSIONS: Our main findings suggest that parity is associated with TNBC across all ages at BC diagnosis, whereas the association with luminal HER2-like BC was present only for early onset BC.

Gene-environment interactions involving functional variants: Results from the Breast Cancer Association Consortium.

Investigating the most likely causal variants identified by fine-mapping analyses may improve the power to detect gene-environment interactions. We assessed the interplay between 70 single nucleotide polymorphisms identified by genetic fine-scale mapping of susceptibility loci and 11 epidemiological breast cancer risk factors in relation to breast cancer. Analyses were conducted on up to 58,573 subjects (26,968 cases and 31,605 controls) from the Breast Cancer Association Consortium, in one of the largest studies of its kind. Analyses were carried out separately for estrogen receptor (ER) positive (ER+) and ER negative (ER-) disease. The Bayesian False Discovery Probability (BFDP) was computed to assess the noteworthiness of the results. Four potential gene-environment interactions were identified as noteworthy (BFDP < 0.80) when assuming a true prior interaction probability of 0.01. The strongest interaction result in relation to overall breast cancer risk was found between CFLAR-rs7558475 and current smoking (ORint = 0.77, 95% CI: 0.67-0.88, pint = 1.8 × 10-4 ). The interaction with the strongest statistical evidence was found between 5q14-rs7707921 and alcohol consumption (ORint =1.36, 95% CI: 1.16-1.59, pint = 1.9 × 10-5 ) in relation to ER- disease risk. The remaining two gene-environment interactions were also identified in relation to ER- breast cancer risk and were found between 3p21-rs6796502 and age at menarche (ORint = 1.26, 95% CI: 1.12-1.43, pint =1.8 × 10-4 ) and between 8q23-rs13267382 and age at first full-term pregnancy (ORint = 0.89, 95% CI: 0.83-0.95, pint = 5.2 × 10-4 ). While these results do not suggest any strong gene-environment interactions, our results may still be useful to inform experimental studies. These may in turn, shed light on the potential interactions observed.

Shared genetics underlying epidemiological association between endometriosis and ovarian cancer.

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci.

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.

Circulating miRNAs with prognostic value in metastatic breast cancer and for early detection of metastasis.

Metastasis is the principal cause of high morbidity and mortality among breast cancer (BC) patients. Identification of markers that can be routinely monitored to predict onset of metastasis in BC patients and prognosis of metastatic breast cancer (MBC) patients would increase their median survival. In this study, plasma miRNAs of 40 MBC patients were profiled by TaqMan low density arrays and miRNAs with prognostic capacity were identified. The candidates were validated initially in the samples of 237 MBC patients and subsequently in 335 samples from an independent study cohort of BC patients. Sixteen miRNAs were established to be significantly associated with overall survival, and were termed as prognostic miRNA panel template (PROMPT). These included miR-141, miR-144, miR-193b, miR-200a, miR-200b, miR-200c, miR-203, miR-210, miR-215, miR-365, miR-375, miR-429, miR-486-5p, miR-801, miR-1260 and miR-1274a. Additionally, 11 of these miRNAs were also associated with progression-free survival. Their prognostic significance was further confirmed in samples from a second study cohort of BC patients. In addition, miR-200a, miR-200b, miR-200c, miR-210, miR-215 and miR-486-5p were found to be significantly associated with onset of metastasis up to 2 years prior to clinical diagnosis in BC patients. We have thus identified panels of miRNAs, which include metastasis promoting miR-200 family and miR-203, as well as oncogenic and tumor-suppressive miRNAs, that can serve as prognostic markers for MBC, and early detection markers of metastasis in BC.

Prognostic value of automated KI67 scoring in breast cancer: a centralised evaluation of 8088 patients from 10 study groups.

BACKGROUND: The value of KI67 in breast cancer prognostication has been questioned due to concerns on the analytical validity of visual KI67 assessment and methodological limitations of published studies. Here, we investigate the prognostic value of automated KI67 scoring in a large, multicentre study, and compare this with pathologists' visual scores available in a subset of patients. METHODS: We utilised 143 tissue microarrays containing 15,313 tumour tissue cores from 8088 breast cancer patients in 10 collaborating studies. A total of 1401 deaths occurred during a median follow-up of 7.5 years. Centralised KI67 assessment was performed using an automated scoring protocol. The relationship of KI67 levels with 10-year breast cancer specific survival (BCSS) was investigated using Kaplan-Meier survival curves and Cox proportional hazard regression models adjusted for known prognostic factors. RESULTS: Patients in the highest quartile of KI67 (>12 % positive KI67 cells) had a worse 10-year BCSS than patients in the lower three quartiles. This association was statistically significant for ER-positive patients (hazard ratio (HR) (95 % CI) at baseline = 1.96 (1.31-2.93); P = 0.001) but not for ER-negative patients (1.23 (0.86-1.77); P = 0.248) (P-heterogeneity = 0.064). In spite of differences in characteristics of the study populations, the estimates of HR were consistent across all studies (P-heterogeneity = 0.941 for ER-positive and P-heterogeneity = 0.866 for ER-negative). Among ER-positive cancers, KI67 was associated with worse prognosis in both node-negative (2.47 (1.16-5.27)) and node-positive (1.74 (1.05-2.86)) tumours (P-heterogeneity = 0.671). Further classification according to ER, PR and HER2 showed statistically significant associations with prognosis among hormone receptor-positive patients regardless of HER2 status (P-heterogeneity = 0.270) and among triple-negative patients (1.70 (1.02-2.84)). Model fit parameters were similar for visual and automated measures of KI67 in a subset of 2440 patients with information from both sources. CONCLUSIONS: Findings from this large-scale multicentre analysis with centrally generated automated KI67 scores show strong evidence in support of a prognostic value for automated KI67 scoring in breast cancer. Given the advantages of automated scoring in terms of its potential for standardisation, reproducibility and throughput, automated methods appear to be promising alternatives to visual scoring for KI67 assessment.

Patient survival and tumor characteristics associated with CHEK2:p.I157T - findings from the Breast Cancer Association Consortium.

BACKGROUND: P.I157T is a CHEK2 missense mutation associated with a modest increase in breast cancer risk. Previously, another CHEK2 mutation, the protein truncating c.1100delC has been associated with poor prognosis of breast cancer patients. Here, we have investigated patient survival and characteristics of breast tumors of germ line p.I157T carriers. METHODS: We included in the analyses 26,801 European female breast cancer patients from 15 studies participating in the Breast Cancer Association Consortium. We analyzed the association between p.I157T and the clinico-pathological breast cancer characteristics by comparing the p.I157T carrier tumors to non-carrier and c.1100delC carrier tumors. Similarly, we investigated the p.I157T associated risk of early death, breast cancer-associated death, distant metastasis, locoregional relapse and second breast cancer using Cox proportional hazards models. Additionally, we explored the p.I157T-associated genomic gene expression profile using data from breast tumors of 183 Finnish female breast cancer patients (ten p.I157T carriers) (GEO: GSE24450). Differential gene expression analysis was performed using a moderated t test. Functional enrichment was investigated using the DAVID functional annotation tool and gene set enrichment analysis (GSEA). The tumors were classified into molecular subtypes according to the St Gallen 2013 criteria and the PAM50 gene expression signature. RESULTS: P.I157T was not associated with increased risk of early death, breast cancer-associated death or distant metastasis relapse, and there was a significant difference in prognosis associated with the two CHEK2 mutations, p.I157T and c.1100delC. Furthermore, p.I157T was associated with lobular histological type and clinico-pathological markers of good prognosis, such as ER and PR expression, low TP53 expression and low grade. Gene expression analysis suggested luminal A to be the most common subtype for p.I157T carriers and CDH1 (cadherin 1) target genes to be significantly enriched among genes, whose expression differed between p.I157T and non-carrier tumors. CONCLUSIONS: Our analyses suggest that there are fundamental differences in breast tumors of CHEK2:p.I157T and c.1100delC carriers. The poor prognosis associated with c.1100delC cannot be generalized to other CHEK2 mutations.

Relationship between menopausal hormone therapy and mortality after breast cancer The MARIEplus study, a prospective case cohort.

Cohort studies of breast cancer (BC) patients, but not of disease-free women at inclusion, have found menopausal hormone therapy (MHT) to be associated with decreased BC specific mortality (BCM). Here, the German population-based MARIEplus BC cohort was analyzed to further elucidate associations of prediagnostic MHT with BCM (and modification by tumor characteristics), recurrence, and secondarily with other cause and overall mortality. Enrolled 2002-2005, incident invasive BC cases (N = 3,321) were followed up for a median of 6.1 years. Cox proportional hazards models adjusted for tumor characteristics, mammography and lifestyle were applied. Compared with never users of MHT, current users at date of diagnosis had significantly lower BCM (Hazard ratio (HR) 0.72, 95% CI 0.53-0.97) and risk of recurrence (HR 0.61, 95% CI 0.46-0.82). The MHT related reduced BCM was confined to patients with low grade tumors (HR 0.44, 95% CI 0.28-0.70; phet = 0.01) and not modified by estrogen receptor or nodal status. BCM decreased with MHT duration in current and increased in past users (phet = 0.015). Mortality due to causes other than BC and overall mortality were also reduced in current MHT users (HR 0.51, 95% CI 0.32-0.81, HR 0.66, 95% CI 0.52-0.86, respectively). Favorable tumor characteristics and mammographic surveillance could not fully explain associations of current MHT use with BCM and recurrence risk. Thus, the study contributes to the evidence that prediagnostic MHT does not have a negative impact on prognosis after BC. The restriction of a reduced BCM to low grade tumors should be confirmed in independent studies.

Plasma hyaluronic acid level as a prognostic and monitoring marker of metastatic breast cancer.

Conventional tumor markers have limited value for prognostication and treatment monitoring in metastatic breast cancer (MBC) patients and novel circulating tumor markers therefore need to be explored. Hyaluronic acid (HA) is a major macropolysaccharide in the extracellular matrix and is reported to be associated with tumor progression. In our study, we investigated plasma HA level with respect to progression free survival (PFS) and overall survival (OS), as well as the treatment monitoring value in MBC patients. The prognostic value of plasma HA level was investigated in a discovery cohort of 212 MBC patients with 2.5-year follow-up and validated in an independent validation cohort of 334 patients with 5-year follow-up. The treatment monitoring value of plasma HA level was investigated in 61 MBC patients from discovery cohort who had been radiographically examined after first complete cycle of chemo therapy. We found a robust association between high plasma HA level and poor prognosis of MBC patients in both discovery (pPFS = 7.92 × 10(-6) and pOS = 5.27 × 10(-5)) and validation studies (pPFS = 3.66 × 10(-4) and pOS = 1.43 × 10(-4)). In the discovery cohort, the plasma HA level displayed independent prognostic value after adjusted for age and clinicopathological factors, with respect to PFS and OS. Further, the decrease of plasma HA level displayed good concordance with treatment response evaluated by radiographic examination (AUC = 0.79). Plasma HA level displays prognostic value, as well as treatment monitoring value for MBC patients.

A splicing variant of TERT identified by GWAS interacts with menopausal estrogen therapy in risk of ovarian cancer.

Menopausal estrogen-alone therapy (ET) is a well-established risk factor for serous and endometrioid ovarian cancer. Genetics also plays a role in ovarian cancer, which is partly attributable to 18 confirmed ovarian cancer susceptibility loci identified by genome-wide association studies. The interplay among these loci, ET use and ovarian cancer risk has yet to be evaluated. We analyzed data from 1,414 serous cases, 337 endometrioid cases and 4,051 controls across 10 case-control studies participating in the Ovarian Cancer Association Consortium (OCAC). Conditional logistic regression was used to determine the association between the confirmed susceptibility variants and risk of serous and endometrioid ovarian cancer among ET users and non-users separately and to test for statistical interaction. A splicing variant in TERT, rs10069690, showed a statistically significant interaction with ET use for risk of serous ovarian cancer (pint = 0.013). ET users carrying the T allele had a 51% increased risk of disease (OR = 1.51, 95% CI 1.19-1.91), which was stronger for long-term ET users of 10+ years (OR = 1.85, 95% CI 1.28-2.66, pint = 0.034). Non-users showed essentially no association (OR = 1.08, 95% CI 0.96-1.21). Two additional genomic regions harboring rs7207826 (C allele) and rs56318008 (T allele) also had significant interactions with ET use for the endometrioid histotype (pint = 0.021 and pint = 0.037, respectively). Hence, three confirmed susceptibility variants were identified whose associations with ovarian cancer risk are modified by ET exposure; follow-up is warranted given that these interactions are not adjusted for multiple comparisons. These findings, if validated, may elucidate the mechanism of action of these loci.

A Comprehensive Multistate Model Analyzing Associations of Various Risk Factors With the Course of Breast Cancer in a Population-Based Cohort of Breast Cancer Cases.

We employed a semi-Markov multistate model for the simultaneous analysis of various endpoints describing the course of breast cancer. Results were compared with those from standard analyses using a Cox proportional hazards model. We included 3,012 patients with invasive breast cancer newly diagnosed between 2001 and 2005 who were recruited in Germany for a population-based study, the Mamma Carcinoma Risk Factor Investigation (MARIE Study), and prospectively followed up until the end of 2009. Locoregional recurrence and distant metastasis were included as intermediate states, and deaths from breast cancer, secondary cancer, and other causes were included as competing absorbing states. Tumor characteristics were significantly associated with all breast cancer-related endpoints. Nodal involvement was significantly related to local recurrence but more strongly related to distant metastases. Smoking was significantly associated with mortality from second cancers and other causes, whereas menopausal hormone use was significantly associated with reduced distant metastasis and death from causes other than cancer. The presence of cardiovascular disease at diagnosis was solely associated with mortality from other causes. Compared with separate Cox models, multistate models allow for dissection of prognostic factors and intermediate events in the analysis of cause-specific mortality and can yield new insights into disease progression and associated pathways.

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.

Post-GWAS gene-environment interplay in breast cancer: results from the Breast and Prostate Cancer Cohort Consortium and a meta-analysis on 79,000 women.

We studied the interplay between 39 breast cancer (BC) risk SNPs and established BC risk (body mass index, height, age at menarche, parity, age at menopause, smoking, alcohol and family history of BC) and prognostic factors (TNM stage, tumor grade, tumor size, age at diagnosis, estrogen receptor status and progesterone receptor status) as joint determinants of BC risk. We used a nested case-control design within the National Cancer Institute's Breast and Prostate Cancer Cohort Consortium (BPC3), with 16 285 BC cases and 19 376 controls. We performed stratified analyses for both the risk and prognostic factors, testing for heterogeneity for the risk factors, and case-case comparisons for differential associations of polymorphisms by subgroups of the prognostic factors. We analyzed multiplicative interactions between the SNPs and the risk factors. Finally, we also performed a meta-analysis of the interaction ORs from BPC3 and the Breast Cancer Association Consortium. After correction for multiple testing, no significant interaction between the SNPs and the established risk factors in the BPC3 study was found. The meta-analysis showed a suggestive interaction between smoking status and SLC4A7-rs4973768 (Pinteraction = 8.84 × 10(-4)) which, although not significant after considering multiple comparison, has a plausible biological explanation. In conclusion, in this study of up to almost 79 000 women we can conclusively exclude any novel major interactions between genome-wide association studies hits and the epidemiologic risk factors taken into consideration, but we propose a suggestive interaction between smoking status and SLC4A7-rs4973768 that if further replicated could help our understanding in the etiology of BC.

A large-scale assessment of two-way SNP interactions in breast cancer susceptibility using 46,450 cases and 42,461 controls from the breast cancer association consortium.

Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70,917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46,450 breast cancer cases and 42,461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10(-4)) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10(-8). Results from the second analytic approach were consistent with those from the first (P > 10(-10)). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

Gene-environment interaction and risk of breast cancer.

Hereditary, genetic factors as well as lifestyle and environmental factors, for example, parity and body mass index, predict breast cancer development. Gene-environment interaction studies may help to identify subgroups of women at high-risk of breast cancer and can be leveraged to discover new genetic risk factors. A few interesting results in studies including over 30,000 breast cancer cases and healthy controls indicate that such interactions exist. Explorative gene-environment interaction studies aiming to identify new genetic or environmental factors are scarce and still underpowered. Gene-environment interactions might be stronger for rare genetic variants, but data are lacking. Ongoing initiatives to genotype larger sample sets in combination with comprehensive epidemiologic databases will provide further opportunities to study gene-environment interactions in breast cancer. However, based on the available evidence, we conclude that associations between the common genetic variants known today and breast cancer risk are only weakly modified by environmental factors, if at all.

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.