A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Formyl peptide receptor-1 (FPR1) represses intestinal oncogenesis

Formyl peptide receptor-1 (FPR1) is a pattern recognition receptor that is mostly expressed by myeloid cells. In patients with colorectal cancer (CRC), a loss-of-function polymorphism (rs867228) in the gene coding for FPR1 has been associated with reduced responses to chemotherapy or chemoradiotherapy. Moreover, rs867228 is associated with accelerated esophageal and colorectal carcinogenesis. Here, we show that dendritic cells from Fpr1-/- mice exhibit reduced migration in response to chemotherapy-treated CRC cells. Moreover, Fpr1-/- mice are particularly susceptible to chronic ulcerative colitis and colorectal oncogenesis induced by the mutagen azoxymethane followed by oral dextran sodium sulfate, a detergent that induces colitis. These experiments were performed after initial co-housing of Fpr1-/- mice and wild-type controls, precluding major Fpr1-driven differences in the microbiota. Pharmacological inhibition of Fpr1 by cyclosporin H also tended to increase intestinal oncogenesis in mice bearing the ApcMin mutation, and this effect was reversed by the anti-inflammatory drug sulindac. We conclude that defective FPR1 signaling favors intestinal tumorigenesis through the modulation of the innate inflammatory/immune response.

Predictive value of vascular endothelial growth factor polymorphisms for maintenance bevacizumab efficacy in metastatic colorectal cancer: an ancillary study of the PRODIGE 9 phase III trial

Background

Several studies have reported the impact of single nucleotide polymorphisms (SNPs) in vascular endothelial growth factor (VEGF) pathway genes on the efficacy of bevacizumab in metastatic colorectal cancer (mCRC), but results are still inconsistent. The PRODIGE 9 phase III study compared bevacizumab maintenance versus observation alone after induction chemotherapy with FOLFIRI plus bevacizumab.

Objective

We evaluated the impact of SNPs of VEGF-A, VEGF receptors (VEGFR-1, VEGFR-2), and hypoxia inducible factor-1α (HIF-1α) on tumor control duration (TCD), overall survival (OS), progression-free survival (PFS), and duration of first chemotherapy free-intervals (CFI).

Patients and methods

We included 314/491 patients from PRODIGE 9 with a DNA blood sample available. Nine SNPs were genotyped on germline DNA using real-time Polymerase Chain Reaction TaqMan TM (Thermo Fisher Scientific, Waltham, MA , USA 02451).

Results

In the bevacizumab arm, patients with the VEGFR-1 rs9582036 CC genotype (n = 14) had significantly longer TCD [22.4 months (95% confidence interval (CI): 14.75-not reached)] than patients with the AA or CA genotype [14.4 months (95% CI: 11.7-17.1)] (p = 0.036), whereas there was no significant difference in the observation arm. In the bevacizumab arm, no significant difference was found between the CC, and AA or CA genotype for OS [28.2 (95% CI: 18.1-42.8) versus 22.5 (95% CI: 18.6-24.6) months, p = 0.5], PFS [9.4 (95% CI: 7.2-11.3) versus 9.2 (95% CI: 8.71-10.1)], and duration of the first CFI [4.6 (95% CI: 1.6-13.3) versus 4.14 (95% CI: 0.5-29.0) months, p = 0.3].

Conclusion

Among mCRC patients treated with bevacizumab maintenance, those with the VEGFR-1 rs9582036 CC genotype experienced longer TCD. The presence of this genotype may thus predict a benefit of bevacizumab maintenance in mCRC.

Postoperative circulating tumor DNA detection is associated with the risk of recurrence in patients resected for a stage II colorectal cancer

Circulating tumor DNA (ctDNA) is reported to be promising in localized colorectal cancer (CRC). The present study aimed to retrospectively evaluate the impact of ctDNA in patients with a resected stage II CRC from the PROGIGE 13 trial with available paired tumor and blood samples. A group of recurrent patients were matched one-to-one with nonrecurrent patients according to sex, tumor location, treatment sequence, and blood collection timing. CtDNA was analyzed by digital PCR according to NGS of tumors. Disease-free survival (DFS) and overall survival (OS) were analyzed based on ctDNA, and the risks of recurrence and death were determined. A total of 134 patients were included, with 67 patients in each group. At least one alteration was identified in 115/134 tumors. Postoperative ctDNA was detected in 10/111 (9.0%) informative samples and was detected more frequently in the recurrent group (16.7% versus 1.8%; p = 0.02). The median DFS of ctDNA+ versus ctDNA- patients was 16.8 versus 54 months (p = 0.002), respectively, and the median OS was 51.3 versus 69.5 months (p = 0.03), respectively. CtDNA was associated with recurrence (ORa = 11.13, p = 0.03) and death (HRa = 3.15, p = 0.01). In conclusion, the presence of postoperative ctDNA is associated with both recurrence and survival in stage II CRC.

Genome-wide interaction analysis of menopausal hormone therapy use and breast cancer risk among 62,370 women

Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values < 5 × 10-8 as genome-wide significant, and p-values < 1 × 10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants. None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values < 1 × 105. The strongest evidence was found for rs4674019 (p-value = 2.27 × 10-7), which showed genome-wide significant interaction (p-value = 3.8 × 10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT-breast cancer risk association.

Prognostic value of circulating tumour DNA in metastatic pancreatic cancer patients: post-hoc analyses of two clinical trials

OBJECTIVE

The prognostication of metastatic pancreatic adenocarcinoma (mPDAC) patients remains uncertain, mainly based on carbohydrate antigen 19-9 (CA19-9), with limited utility. Circulating tumour DNA (ctDNA) has been suggested as a prognostic factor, but its added value has been poorly explored. The objective was to determine whether ctDNA is an independent factor for the prognostication of mPDAC.

DESIGN

Translational study based on two prospective collections of plasma samples of mPDAC patients naïve for chemotherapy. One used as a test series and the other as validation series coming from two randomised trials (Prodige 35 and Prodige 37). CtDNA was assessed by digital droplet PCR targeting two methylated markers (HOXD8 and POU4F1) according to a newly developed and validated method. Univariate and multivariate analyses were performed according to ctDNA status.

RESULTS

Of 372 plasma samples available, 354 patients were analyzed for survival. In the validation series, 145 of 255 patients were found ctDNA positive (56.8%), Median PFS and OS were 5.3 and 8.2 months in ctDNA-positive and 6.2 and 12.6 months in ctDNA-negative patients, respectively. ctDNA positivity was more often associated with young age, high CA19-9 level and neutrophils lymphocytes ratio. In multivariate analysis including these previous markers, ctDNA was confirmed as an independent prognostic marker for PFS (adjusted hazard ratio (HR) 1.5, CI 95% [1.03-2.18], p = 0.034) and OS (HR 1.62, CI 95% [1.05-2.5], p = 0.029).

CONCLUSIONS

In this first ctDNA assessment in a large series of mPDAC derived from clinical trials, ctDNA was detectable in 56.8% of patients and confirmed as an independent prognostic marker.

Circulating tumor DNA is a prognostic marker of tumor recurrence in stage II and III colorectal cancer: multicentric, prospective cohort study (ALGECOLS)

BACKGROUND

In non-metastatic colorectal cancer (CRC), we evaluated prospectively the pertinence of longitudinal detection and quantification of circulating tumor DNA (ctDNA) as a prognostic marker of recurrence.

METHOD

The presence of ctDNA was assessed from plasma collected before and after surgery for 184 patients classified as stage II or III and at each visit during 3-4 years of follow-up. The ctDNA analysis was performed by droplet-based digital polymerase chain reaction, targeting mutation and methylation markers, blindly from the clinical outcomes. Multivariate analyses were adjusted on age, gender, stage, and adjuvant chemotherapy.

RESULTS

Before surgery, 27.5% of patients were positive for ctDNA detection. The rate of recurrence was 32.7% and 11.6% in patients with or without detectable ctDNA respectively (P = 0.001). Time to recurrence (TTR) was significantly shorter in patients with detectable ctDNA before (adjusted hazard ratio [HR] = 3.58, 95% confidence interval [CI] 1.71-7.47) or immediately after surgery (adjusted HR = 3.22, 95% CI 1.32-7.89). The TTR was significantly shorter in patients with detectable ctDNA during the early postoperative follow-up (1-6 months) (adjusted HR = 5, 95% CI 1.9-12.9). Beyond this period, ctDNA remained a prognostic marker with a median anticipated diagnosis of recurrence of 13.1 weeks (interquartile range 28 weeks) when compared to imaging follow-up. The rate of ctDNA+ might be underestimated knowing that consensus pre-analytical conditions were not described at initiation of the study.

CONCLUSION

This prospective study confirms the relevance of ctDNA as a recurrence risk factor in stage II and III CRC before surgery and as a marker of minimal residual disease after surgery that may predict recurrence several months before imaging techniques.

Prognostic value of the PrPC-ILK-IDO1 axis in the mesenchymal colorectal cancer subtype

The CMS4 mesenchymal subtype of colorectal cancer (CRC) is associated with poor prognosis and resistance to treatment. The cellular prion protein PrP^C is overexpressed in CMS4 tumors and controls the expression of a panel of CMS4-specific genes in CRC cell lines. Here, we sought to investigate PrP^C downstream pathways that may underlie its role in CMS4 CRC. By combining gene set enrichment analyses and gain and loss of function approaches in CRC cell lines, we identify the integrin-linked kinase ILK as a proximal effector of PrP^C that mediates its control on the CMS4 phenotype. We further leveraged three independent large CRC cohorts to assess correlations in gene expression pattern with patient outcomes and found that ILK is overexpressed in CMS4 mesenchymal tumors and confers a poor prognosis, especially when combined with high expression of the PrP^C encoding gene PRNP. Of note, we discovered that the PrP^C-ILK signaling axis controls the expression and activity of the tryptophan metabolizing enzyme indoleamine 2,3 dioxygenase IDO1, a key player in immune tolerance. In addition, we monitored alterations in the levels of tryptophan and its metabolites of the kynurenine pathway in the plasma of metastatic CRC patients (n = 325) and we highlight their prognostic value in combination with plasma PrP^C levels. Thus, the PrP^C-ILK-IDO1 axis plays a key role in the mesenchymal subtype of CRC. PrP^C and IDO1-targeted strategies may represent new avenues for patient stratification and treatment in CRC.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 and 16q23.2, which were associated with thyroid-stimulating hormone levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Prognostic value and relation with adjuvant treatment duration of ctDNA in stage III colon cancer: a post-hoc analysis of the PRODIGE-GERCOR IDEA-France trial

PURPOSE

Circulating tumor DNA (ctDNA) has been suggested as a major prognostic factor in resected stage III colon cancer. We analyzed ctDNA of patients randomized in the phase III IDEA-France trial.

EXPERIMENTAL DESIGN

ctDNA was tested for WIF1 and NPY by droplet digital PCR with method developed and validated for colorectal cancer. Disease-free survival (DFS) and overall survival (OS) were analyzed via multivariable analysis in patients with ctDNA samples and in sub-groups according to treatment duration (3/6 months) and disease stage (high/low-risk stage III).

RESULTS

Of 2010 randomized patients, 1345 had available ctDNA samples (1017 collected both post-surgery and pre-chemotherapy). More ECOG PS of 0 (78% versus 69%) and T4 and/or N2 (40% versus 36%) were observed in patients studied (n=1017) versus not analyzed (n=993). There were 877 ctDNA-negative (86.2%) and 140 ctDNA-positive (13.8%) patients; their baseline characteristics were similar. With a median follow-up of 6.6 years, the 3-year DFS rate was 66.39% for ctDNA-positive patients and 76.71% for ctDNA-negative patients (P=0.015). ctDNA was confirmed as an independent prognostic marker for DFS (adjusted HR=1.55, 95% CI 1.13-2.12, P=0.006) and OS (HR=1.65, 95% CI 1.12-2.43, P=0.011). ctDNA was prognostic in patients treated for 3 months and with T4 and/or N2 tumors, but not in those treated for 6 months and with T1-3/N1 tumors.

CONCLUSIONS

In this first ctDNA assessment of a large series of stage III colon cancer patients enrolled in phase III trial, post-surgery ctDNA was found in 13.8% of them and was confirmed as an independent prognostic marker.

Role of DNA Repair Variants and Diagnostic Radiology Exams in Differentiated Thyroid Cancer Risk: A Pooled Analysis of Two Case-Control Studies

BACKGROUND

Given the increased use and diversity of diagnostic procedures, it is important to understand genetic susceptibility to radiation-induced thyroid cancer.

METHODS

On the basis of self-declared diagnostic radiology examination records in addition to existing literature, we estimated the radiation dose delivered to the thyroid gland from diagnostic procedures during childhood and adulthood in two case-control studies conducted in France. A total of 1,071 differentiated thyroid cancer (DTC) cases and 1,188 controls from the combined studies were genotyped using a custom-made Illumina OncoArray DNA chip. We focused our analysis on variants in genes involved in DNA damage response and repair pathways, representing a total of 5,817 SNPs in 571 genes. We estimated the OR per milli-Gray (OR/mGy) of the radiation dose delivered to the thyroid gland using conditional logistic regression. We then used an unconditional logistic regression model to assess the association between DNA repair gene variants and DTC risk. We performed a meta-analysis of the two studies.

RESULTS

The OR/mGy was 1.02 (95% confidence interval, 1.00-1.03). We found significant associations between DTC and rs7164173 in CHD2 (P = 5.79 × 10-5), rs6067822 in NFATc2 (P = 9.26 × 10-5), rs1059394 and rs699517 both in ENOSF1/THYS, rs12702628 in RPA3, and an interaction between rs7068306 in MGMT and thyroid radiation doses (P = 3.40 × 10-4).

CONCLUSIONS

Our results suggest a role for variants in CDH2, NFATc2, ENOSF1/THYS, RPA3, and MGMT in DTC risk.

IMPACT

CDH2, NFATc2, ENOSF1/THYS, and RPA3 have not previously been shown to be associated with DTC risk.

Characterization of Plasma Cell-Free DNA Integrity Using Droplet-Based Digital PCR: Toward the Development of Circulating Tumor DNA-Dedicated Assays

Background: Cellular-cell free-DNA (ccfDNA) is being explored as a diagnostic and prognostic tool for various diseases including cancer. Beyond the evaluation of the ccfDNA mutational status, its fragmentation has been investigated as a potential cancer biomarker in several studies. However, probably due to a lack of standardized procedures dedicated to preanalytical and analytical processing of plasma samples, contradictory results have been published.

Methods: ddPCR assays allowing the detection of KRAS wild-type and mutated sequences (KRAS p.G12V, pG12D, and pG13D) were designed to target different fragments sizes. Once validated on fragmented and non-fragmented DNA extracted from cancer cell lines, these assays were used to investigate the influence of the extraction methods on the non-mutated and mutated ccfDNA integrity reflected by the DNA integrity index (DII). The DII was then analyzed in two prospective cohorts of metastatic colorectal cancer patients (RASANC study n = 34; PLACOL study n = 12) and healthy subjects (n = 49).

Results and Discussion: Our results demonstrate that ccfDNA is highly fragmented in mCRC patients compared with healthy individuals. These results strongly suggest that the characterization of ccfDNA integrity hold great promise toward the development of a universal biomarker for the follow-up of mCRC patients. Furthermore, they support the importance of standardization of sample handling and processing in such analysis.

Gene network and biological pathways associated with susceptibility to differentiated thyroid carcinoma

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW_dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 2q35 and 8p12 in Europeans, Melanesians and Polynesians

Differentiated thyroid carcinoma (DTC) incidence is characterized by wide ethnic and geographic variations, with high incidence rates observed in Oceanian populations. Genome-wide association studies (GWAS) identified mainly four DTC susceptibility loci at 9q22.33, 14q13.3, 2q35 and 8p12. Here we performed fine-mapping of the 2q35 and 8p12 loci in the population of the EPITHYR consortium that includes Europeans, Melanesians and Polynesians to identify likely causal variants for DTC risk. We conducted a colocalization analysis using eQTLs data to determine the SNPs with the highest probability of causality. At 2q35, we highlighted rs16857609 located in DIRC3. This SNP has a high probability of causality in the three populations, and a significant association in Europeans (OR = 1.4, p = 1.9 x 10-10). It is also associated with expression of DIRC3 and of the nearby gene IGFBP5 in thyroid tumour cells. At 8p12, we identified rs7844425 which was significantly associated with DTC in Europeans (OR = 1.32, p = 7.6 x 10-8) and rs2439304, which was highlighted by the colocalization analysis but only moderately associated with DTC in our dataset (OR = 1.2, p = 0.001). These SNPs are linked to the expression of NRG1 in thyroid tissue. Hence, our study identified novel variants at 2q35 and 8p12 to be prioritized for further functional studies.

Circadian genes polymorphisms, night work and prostate cancer risk: Findings from the EPICAP study

Over the past two decades, several studies have attempted to understand the hypothesis that disrupting the circadian rhythm may promote the development of cancer. Some have suggested that night work and some circadian genes polymorphisms are associated with cancer, including prostate cancer. Our study aims to test the hypothesis that prostate cancer risk among night workers may be modulated by genetic polymorphisms in the circadian pathway genes based on data from the EPICAP study, a population-based case-control study including 1511 men (732 cases/779 controls) with genotyped data. We estimated odds ratio (ORs) and P values of the association between prostate cancer and circadian gene variants using logistic regression models. We tested the interaction between circadian genes variants and night work indicators that were significantly associated with prostate cancer at pathway, gene and SNP levels. Analyses were also stratified by each of these night work indicators and by cancer aggressiveness. The circadian pathway was significantly associated with aggressive prostate cancer among night workers (P = .004), particularly for men who worked at night for <20 years (P = .0002) and those who performed long night shift (>10 hours, P = .001). At the gene level, we observed among night workers significant associations between aggressive prostate cancer and ARNTL, NPAS2 and RORA. At the SNP-level, no significant association was observed. Our findings provide some clues of a potential modulating effect of circadian genes in the relationship between night work and prostate cancer. Further investigation is warranted to confirm these findings and to better elucidate the biological pathways involved.

A TLR3 Ligand Reestablishes Chemotherapeutic Responses in the Context of FPR1 Deficiency

For anthracycline-based chemotherapy to be immunogenic, dying cancer cells must release annexin A1 (ANXA1) that subsequently interacts with the pattern recognition receptor, formyl peptide receptor 1 (FPR1), on the surface of dendritic cells (DC). Approximately 30% of individuals bear loss-of-function alleles of FPR1, calling for strategies to ameliorate their anticancer immune response. Here, we show that immunotherapy with a ligand of Toll-like receptor-3, polyinosinic:polycytidylic acid (pIC), restores the deficient response to chemotherapy of tumors lacking ANXA1 developing in immunocompetent mice or those of normal cancers growing in FPR1-deficient mice. This effect was accompanied by improved DC- and T-lymphocyte-mediated anticancer immunity. Of note, carcinogen-induced breast cancers precociously developed in FPR1-deficient mice as compared with wild-type controls. A similar tendency for earlier cancer development was found in patients carrying the loss-of-function allele of FPR1. These findings have potential implications for the clinical management of FPR1-deficient patients. SIGNIFICANCE: The loss-of-function variant rs867228 in FPR1, harbored by approximately 30% of the world population, is associated with the precocious manifestation of breast, colorectal, esophageal, and head and neck carcinomas. pIC restores deficient chemotherapeutic responses in mice lacking Fpr1, suggesting a personalized strategy for compensating for the FPR1 defect.This article is highlighted in the In This Issue feature, p. 211.

Biochemically Tracked Variability of Blood Plasma Thawed-State Exposure Times in a Multisite Collection Study

The integrity of blood plasma/serum (P/S) specimens can be impacted by preanalytical handling and storage conditions that result in thawed-state exposures (> -30°C). We recently reported a simple dilute-and-shoot, intact-protein liquid chromatography/mass spectrometry (LC/MS) assay called ΔS-Cys-Albumin that quantifies cumulative exposure of P/S to thawed conditions based on the change in relative abundance of the oxidized (S-cysteinylated) proteoform of albumin (S-Cys-Albumin) in the native sample to that of an aliquot of the sample intentionally driven to its maximum oxidation state. Herein, we evaluated the effect of prestorage delay and initial storage temperature on sample integrity by applying the ΔS-Cys-Albumin assay to a set of plasma samples (n = 413) collected under a single clinical study but from 12 different collection sites. Major differences (p < 0.0001) were observed between different groups of samples with modestly inconsistent initial handling conditions (i.e., initial processing of whole blood to plasma and placement at -80°C completed in under 3 hours, 3-13 hours, and over 17 hours). ΔS-Cys-Albumin was significantly inversely correlated with delay time at 4°C before centrifugation and total delay before final storage at -80°C (p < 0.0001). Samples from two collection sites had much lower ΔS-Cys-Albumin values relative to samples from other sites, in accordance with the fact that they were stored at -20°C for an average of 7.6 months before shipment to the central repository for final storage at -80°C. Based on the rate law for S-Cys-Albumin formation in plasma ex vivo, the average time that each plasma specimen had been exposed to the equivalent of room temperature (23°C) was back calculated from the measured ΔS-Cys-Albumin values. A survey of clinical analytes in P/S whose measured concentrations are sensitive to the initial handling/storage conditions documented in this study is provided and the ramifications of the plasma integrity findings from this multisite clinical study are discussed.

Circulating Tumor DNA is Prognostic and Potentially Predictive of Eryaspase Efficacy in Second-line in Patients with Advanced Pancreatic Adenocarcinoma

PURPOSE

Eryaspase is composed of l-asparaginase encapsulated in erythrocytes and has demonstrated significant efficacy in a randomized phase II trial. We assessed the prognostic and predictive value of circulating tumor DNA (ctDNA) in patients, plasma included in this trial.

EXPERIMENTAL DESIGN

Samples prospectively collected pretreatment were centrally analyzed by next-generation sequencing. Prognostic values of baseline ctDNA and ctDNA early changes between day 0 and 28 were assessed in both arms combined on objective response rate (ORR), progression-free survival (PFS), and overall survival (OS); three groups were defined: negative ctDNA (Neg), ctDNA responders (Resp), and ctDNA nonresponders (NResp). Predictive value of ctDNA for eryaspase efficacy was investigated.

RESULTS

ctDNA was positive at baseline in 77 patients of the 113 tested patients (68%). Detectable ctDNA was an independent negative prognostic factor for OS (4.6 vs. 8.8 months; P = 0.0025) and PFS (1.6 vs. 3.3 months; P = 0.00043). Early change in ctDNA levels was correlated with ORR (20%, 26%, 0%; P < 0.04), PFS (3.7, 3.4, 1.6 months; P < 0.0001), and OS (11.7, 6.5, 4.3 months; P < 0.0001) according to the three defined groups (Neg, Res, NResp, respectively). In patients with ctDNA detectable at baseline, eryaspase was associated with better PFS [HR = 0.53; 95% confidence interval (CI): 0.3-0.94)] and OS (HR = 0.52; 95% CI: 0.29-0.91).

CONCLUSIONS

We confirm from a prospective randomized trial that: (i) the presence of ctDNA at baseline is a major prognostic factor, (ii) the early change of ctDNA correlates with treatment outcome, and (iii) the ctDNA could be a predictive biomarker of eryaspase efficacy.

Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 × 10-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

Circulating tumor DNA is prognostic and potentially predictive of eryaspase efficacy in patients with advanced pancreatic adenocarcinoma

4617

Background: Eryaspase is composed of L-asparaginase encapsulated in erythrocytes. It has demonstrated significant efficacy in combination with chemotherapy in a randomized phase 2 trial in second-line in patients with advanced pancreatic adenocarcinoma. We assessed, in this study, the prognostic and predictive value of circulating tumor DNA (ctDNA) in plasma samples of patients included in the eryaspase phase 2 trial. Methods: Samples prospectively collected pre-treatment at each 28-day cycle were centrally analyzed by next-generation sequencing (BPER method). Prognostic values of baseline ctDNA and ctDNA early changes between day 0 and 28 were assessed in both arms combined on objective response rate (ORR), progression free survival (PFS) and overall survival (OS). We conducted interaction test between ctDNA positivity and treatment arm, and the predictive value of ctDNA for eryaspase efficacy was investigated. Results: Patients with at least one available plasma sample have been included (n = 122/141). The presence of ctDNA at baseline was identified in 68% (77/113) of patients and was an independent negative prognostic factor for OS (4.6 vs 8.8 months; p = 0.0025) and PFS (1.6 vs 3.3 months; p = 0.00043). Early change in ctDNA levels was assessed by separating patients into three categories (one without detectable ctDNA, and two according to radio median value between day 0 and day 28) that were significantly correlated with ORR, PFS and OS. A significant interaction was observed between the presence of ctDNA and eryaspase efficacy. In patients with ctDNA detectable at baseline, eryaspase was associated with better PFS (HR = 0.53; 95% CI: 0.3-0.94) and OS (HR = 0.52; 95% CI: 0.29-0.91). Conclusions: We confirm from a prospective randomized trial that 1/ the presence of ctDNA at baseline is a major prognostic factor, 2/ the early change of ctDNA correlates with treatment outcome and 3/ the ctDNA could be a predictive biomarker of eryaspase efficacy. Clinical trial information: NCT02195180 .

Role of GSTM1 and GSTT1 genotypes in differentiated thyroid cancer and interaction with lifestyle factors: Results from case-control studies in France and New Caledonia

Background

GSTM1 and GSTT1 are involved in detoxification of xenobiotics, products of oxidative stress and in steroid hormones metabolism. We investigated whether GSTM1 and GSTT1 gene deletion was associated with DTC risk and explored interaction with non-genetic risk factors of DTC.

Methods

The study included 661 DTC cases and 736 controls from two case-control studies conducted in France and New Caledonia. Odds ratios (OR) and their confidence interval (CI) for DTC associated with GST genotypes, alcohol drinking, tobacco smoking, body mass index and hormonal factors were calculated using logistic regression models.

Results

Results are presented for Europeans and Melanesians combined, as no heterogeneity between groups was detected. We found that DTC risk increased with obesity and decrease with alcohol drinking. After stratification by gene deletion status, the OR for obesity was 5.75, (95%CI 2.25–14.7) among individuals with GSTT1 and GSTM1-deleted genotype, and 1.26, (95%CI 0.89–1.77) in carriers of both genes (p-interaction = 0.02). The OR for drinking ≥1 glass/week was 0.33 (95%CI 0.15–0.74) in GSTT1-null individuals while it was 1.01 (95%CI 0.67–1.52) in non-null carriers of the gene (p-interaction = 0.01). No interaction between GST genotypes and other non-genetic risk factors was detected.

Conclusion

GSTM1 and GSTT1 genotypes may modulate the DTC risk associated with BMI and alcohol consumption.

The cellular prion protein controls the mesenchymal-like molecular subtype and predicts disease outcome in colorectal cancer

Background

Comprehensive transcriptomic analyses have shown that colorectal cancer (CRC) is heterogeneous and have led to the definition of molecular subtypes among which the stem-cell, mesenchymal-like group is associated with poor prognosis. The molecular pathways orchestrating the emergence of this subtype are incompletely understood. In line with the contribution of the cellular prion protein PrP^C to stemness, we hypothesize that deregulation of this protein could lead to a stem-cell, mesenchymal-like phenotype in CRC.

Methods

We assessed the distribution of the PrP^C-encoding PRNP mRNA in two large CRC cohorts according to molecular classification and its association with patient survival. We developed cell-based assays to explore the impact of gain and loss of PrP^C function on markers of the mesenchymal subtype and to delineate the signalling pathways recruited by PrP^C. We measured soluble PrP^C in the plasmas of 325 patients with metastatic CRC and probed associations with disease outcome.

Findings

We found that PRNP gene expression is enriched in tumours of the mesenchymal subtype and is associated with poor survival. Our in vitro analyses revealed that PrP^C controls the expression of genes that specify the mesenchymal subtype through the recruitment of the Hippo pathway effectors YAP and TAZ and the TGFß pathway. We showed that plasma levels of PrP^C are elevated in metastatic CRC and are associated with poor disease control.

Interpretation

Our findings define PrP^C as a candidate driver of the poor-prognosis mesenchymal subtype of CRC. They suggest that PrP^C may serve as a potential biomarker for patient stratification in CRC.

Funding

Grant support was provided by the following: Cancéropôle Ile de France (grant number 2016-1-EMERG-36-UP 5-1), Association pour la Recherche sur le Cancer (grant number PJA 20171206220), SATT Ile de France Innov (grant number 415) as well as INSERM.

Association of breast cancer risk with polymorphisms in genes involved in the metabolism of xenobiotics and interaction with tobacco smoking: A gene-set analysis

Single nucleotide polymorphisms (SNPs) in genes involved in xenobiotics metabolism (XM) are suspected to play a role in breast cancer risk. However, previous findings based on a SNP by SNP approach need to be replicated taking into account the combined effects of multiple SNPs. We used a gene-set analysis method to study the association between breast cancer risk and genetic variation in XM genes (seen as a set of SNPs) and in the XM pathway (seen as a set of genes). We also studied the interaction between variants in XM genes and tobacco smoking. The analysis was conducted in a case-control study of 1,125 cases and 1,172 controls. Using a dedicated chip, genotyping data of 585 SNPs in 68 XM genes were available. Genetic variation in the whole XM pathway was significantly associated with premenopausal breast cancer risk (p = 0.008). This association was mainly driven by genetic variation in NAT2, CYP2C18, CYP2C19, AKR1C2 and ALDH1A3. The association between the XM gene pathway and breast cancer was observed among current and previous smokers, but not among never smokers (p = 0.013 for interaction between XM genes and tobacco smoking status). The association with breast cancer risk indicates that XM genes variants may play a role in breast carcinogenesis through their detoxification function of environmental pollutants, such as those contained in tobacco smoke.

Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Circadian genes and risk of prostate cancer: Findings from the EPICAP study

Circadian rhythms regulate several physiological functions and genes controlling the circadian rhythm were found to regulate cell proliferation, cell cycle and apoptosis. Few studies have investigated the role of those circadian genes in prostate cancer occurrence. We aim to investigate the relationship between circadian genes polymorphisms and prostate cancer risk based on data from the EPICAP study, a population-based case-control study including 1,515 men (732 cases / 783 controls) with genotyped data. Odds Ratios (ORs) for association between prostate cancer and circadian gene variants were estimated for each of the 872 single nucleotide polymorphisms (SNPs) in 31 circadian clock genes. We also used a gene-based and pathway-based approach with a focus on the pathway including 9 core circadian genes. Separate analyses were conducted by prostate cancer aggressiveness. The core-circadian pathway (p = 0.0006) was significantly associated to prostate cancer, for either low (p = 0.002) or high (p = 0.01) grade tumor. At the gene level, we observed significant associations between all prostate cancer and NPAS2 and PER1 after correcting for multiple testing, while only RORA was significant for aggressive tumors. At the SNP-level, no significant association was observed. Our findings provide additional evidence of a potential link between genetic variants in circadian genes and prostate cancer risk. Further investigation is warranted to confirm these findings and to better understand the biological pathways involved.

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

Association of Prognostic Value of Primary Tumor Location in Stage III Colon Cancer With RAS and BRAF Mutational Status

Importance

We know of no data on the prognostic value of primary tumor location (PTL) according to BRAF, RAS, and microsatellite instability (MSI) status in patients who have undergone resection for colon cancer (CC) and have been treated with current standard adjuvant chemotherapy.

Objective

To determine the prognostic and predictive value of PTL according to BRAF, RAS, and MSI status in patients with stage III CC receiving adjuvant treatment with FOLFOX (folinic acid [leucovorin calcium], fluorouracil, and oxaliplatin) with or without cetuximab.

Design, Setting, and Participants

This post hoc analysis included patients with available tumor blocks of resected stage III colon adenocarcinoma who participated in the Pan-European Trials in Alimentary Tract Cancer (PETACC)-8 phase 3 randomized trial. Among the 2559 patients who underwent randomization, 1900 were screened by next-generation sequencing, which showed that 1869 had full information concerning PTL. We categorized primary tumor site as located proximal (right) or distal (left) to the splenic flexure.

Main Outcomes and Measures

The associations between PTL (right- vs left-sided) and disease-free survival (DFS), survival after relapse (SAR), and overall survival (OS) were assessed by Cox models and adjusted for clinical and pathological features, treatment, and MSI, BRAF, and RAS status.

Results

Among the 1869 patients (1056 [57%] male; mean [SD] age, 59.4 [9.5] years) with full molecular data analyzed, 755 (40%) had a right-sided tumor, 164 (10%) had MSI, 942 (50%) had RAS mutations, and 212 (11%) had BRAF mutations. Right-sided tumor location was not prognostic for DFS in the whole population but was associated with a shorter SAR (hazard ratio [HR], 1.54; 95% CI, 1.23-1.93; P = .001) and OS (HR, 1.25; 95% CI, 1.02-1.54; P = .03). When looking at DFS in the different molecular subgroups, we found similar results for microsatellite-stable tumors and tumors with MSI; a better DFS in right-sided vs left-sided tumors in patients with RAS mutations (HR, 0.80; 95% CI, 0.64-1.00; P = .046); and a worse DFS in right-sided vs left-sided tumors in patients with RAS and BRAF double wild type (HR, 1.39; 95% CI, 1.01-1.92; P = .04). These results were found independently of the treatment received, and no beneficial effect of cetuximab on DFS or OS was observed in left-sided tumors.

Conclusions and Relevance

Although right-sided tumor location is associated with poor survival in patients with metastatic CC as previously reported, the association with disease recurrence appears to vary for patients with stage III CC and RAS or BRAF mutations vs those with double wild type.

Prognostic Value of Methylator Phenotype in Stage III Colon Cancer Treated with Oxaliplatin-based Adjuvant Chemotherapy

Purpose: There are conflicting results concerning the prognostic value of the CpG island methylator phenotype (CIMP) in patients with nonmetastatic colon cancer. We studied this phenotype in stage III colon cancer characterized for mismatch repair (MMR), RAS, and BRAF status, and treated with adjuvant FOLFOX-based regimen.Experimental Design: Tumor samples of 1,907 patients enrolled in the PETACC-8 adjuvant phase III trial were analyzed. The method used was methylation-specific PCR, where CIMP⁺ status was defined by methylation of at least 3 of 5 following genes: IGF2, CACNA1G, NEUROG1, SOCS1, and RUNX3 Association between CIMP status and overall survival (OS), disease-free survival (DFS), and survival after recurrence (SAR), was assessed by Cox model adjusted for prognostic factors and treatment arm (FOLFOX4 ± cetuximab).Results: CIMP status was successfully determined in 1,867 patients (97.9%): 275 (14.7%) tumors were CIMP⁺ Compared with CIMP^- patients, CIMP⁺ patients were more frequently older (P = 0.002), females (P = 0.04), with right-sided (P < 0.0001), grade 3-4 (P < 0.0001), pN2 (P = 0.001), dMMR (P < 0.0001), BRAF mutated (P < 0.0001), and RAS wild-type (P < 0.0001) tumors. In multivariate analysis, CIMP⁺ status was associated with shorter OS [HR, 1.46; 95% confidence interval (CI), 1.02-1.94; P = 0.04] and SAR [HR, 1.76; 95% CI, 1.20-2.56; P < 0.0004]; but not DFS [HR, 1.15; 95% CI, 0.86-1.54; P = 0.34]. A nonsignificant trend of detrimental effect of cetuximab was observed in patients with CIMP⁺ tumors for OS, DFS, and SAR.Conclusions: In a large cohort of well-defined patients with stage III colon cancer, CIMP⁺ phenotype is associated with a shorter OS and SAR but not to DFS. Clin Cancer Res; 24(19); 4745-53. ©2018 AACR.

Smoking and Parkinson disease: Evidence for gene-by-smoking interactions

OBJECTIVE

To investigate whether cigarette smoking interacts with genes involved in individual susceptibility to xenobiotics for the risk of Parkinson disease (PD).

METHODS

Two French population-based case-control studies (513 patients, 1,147 controls) were included as a discovery sample to examine gene-smoking interactions based on 3,179 single nucleotide polymorphisms (SNPs) in 289 genes involved in individual susceptibility to xenobiotics. SNP-by-cigarette smoking interactions were tested in the discovery sample through an empirical Bayes (EB) approach. Nine SNPs were selected for replication in a population-based case-control study from California (410 patients, 845 controls) with standard logistic regression and the EB approach. For SNPs that replicated, we performed pooled analyses including the discovery and replication datasets and computed pooled odds ratios and confidence intervals (CIs) using random-effects meta-analysis.

RESULTS

Nine SNPs interacted with smoking in the discovery dataset and were selected for replication. Interactions of smoking with rs4240705 in the RXRA gene and rs1900586 in the SLC17A6 gene were replicated. In pooled analyses (logistic regression), the interactions between smoking and rs4240705-G and rs1900586-G were 1.66 (95% CI 1.28-2.14, p = 1.1 × 10^-4, p for heterogeneity = 0.366) and 1.61 (95% CI 1.17-2.21, p = 0.003, p for heterogeneity = 0.616), respectively. For both SNPs, while smoking was significantly less frequent in patients than controls in AA homozygotes, this inverse association disappeared in G allele carriers.

CONCLUSIONS

We identified and replicated suggestive gene-by-smoking interactions in PD. The inverse association of smoking with PD was less pronounced in carriers of minor alleles of both RXRA-rs4240705 and SLC17A6-rs1900586. These findings may help identify biological pathways involved in the inverse association between smoking and PD.

A “chamber of errors” adaptation to assess pharmaceutical assistants’ knowledge in chemotherapy preparation

French preparation guidelines state that pharmacy staff who manipulate cytotoxic drugs have to follow specific training. In order to assess the pharmaceutical assistants’ skills and knowledge, we developed a “Cytotoxic Preparation Centralized Unit (CPCU) of errors,” derived from the Canadian concept of “Chamber of horrors.” A table listing 20 mistakes to track down was created and each pharmaceutical assistant spent 20 min in the “CPCU of errors” with the pharmacist, who wrote down the spotted mistakes in real time. Among the 21 trained pharmaceutical assistants, 15 were evaluated. On average, 11.9 mistakes on 20 were detected. The lowest score was 7 spotted errors on 20 and the highest was 16 on 20. Those results should be qualified depending on pharmaceutical assistants’ years of experience in the preparation of chemotherapy. Those results may be explained by the way the role-playing was conducted. The simulation was not conducted during an actual preparation using the usual equipment. One of the major obstacles was the difficulty to clear some time for this project because its realization required a full-time pharmacist and the referring pharmaceutical assistant in addition to the evaluated pharmaceutical assistants. Overall, the staff feedback was positive and the role-playing led to a reminder of theoretical knowledge and the good use of some devices. It would be interesting to develop this type of project through a regional oncology network to create a medium that can be used by other hospitals.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

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.

Pooled analysis of the HLA-DRB1 by smoking interaction in Parkinson disease

OBJECTIVE

Inflammatory response plays an important role in Parkinson disease (PD). Previous studies have reported an association between human leukocyte antigen (HLA)-DRB1 and the risk of PD. There has also been growing interest in investigating whether inflammation-related genes interact with environmental factors such as smoking to influence PD risk. We performed a pooled analysis of the interaction between HLA-DRB1 and smoking in PD in 3 population-based case-control studies from Denmark and France.

METHODS

We included 2,056 cases and 2,723 controls from 3 PD studies (Denmark, France) that obtained information on smoking through interviews. Genotyping of the rs660895 polymorphism in the HLA-DRB1 region was based on saliva or blood DNA samples. To assess interactions, we used logistic regression with product terms between rs660895 and smoking. We performed random-effects meta-analysis of marginal associations and interactions.

RESULTS

Both carrying rs660895-G (AG vs AA: odds ratio [OR] = 0.81; GG vs AA: OR = 0.56; p-trend = 0.003) and ever smoking (OR = 0.56, p < 0.001) were inversely associated with PD. A multiplicative interaction was observed between rs660895 and smoking using codominant, additive (interaction parameter = 1.37, p = 0.005), and dominant (interaction parameter = 1.54, p = 0.001) genetic models without any heterogeneity (I² = 0.0%); the inverse association of rs660895-(AG+GG) with PD seen in never smokers (OR = 0.64, p < 0.001) disappeared among ever smokers (OR = 1.00, p = 0.99). Similar interactions were observed when we investigated light and heavy smokers separately.

INTERPRETATION

Our study provides the first evidence that smoking modifies the previously reported inverse association of rs660895-G with PD, and suggests that smoking and HLA-DRB1 are involved in common pathways, possibly related to neuroinflammation. Ann Neurol 2017;82:655-664.

Early Evaluation of Circulating Tumor DNA as Marker of Therapeutic Efficacy in Metastatic Colorectal Cancer Patients (PLACOL Study)

Purpose: Markers of chemotherapy efficacy in metastatic colorectal cancer (mCRC) are essential for optimization of treatment strategies. We evaluated the applicability of early changes in circulating tumor DNA (ctDNA) as a marker of therapeutic efficacy.Experimental Design: This prospective study enrolled consecutive patients with mCRC receiving a first- or second-line chemotherapy. CtDNA was assessed in plasma collected before the first (C₀), second (C₁) and/or third (C₂) chemotherapy cycle, using picodroplet-digital PCR assays based either on detection of gene mutation (KRAS, BRAF, TP53) or hypermethylation (WIF1, NPY). CT scans were centrally assessed using RECIST v1.1 criteria. Multivariate analyses were adjusted on age, gender, ECOG performance status (PS), metastatic synchronicity, and treatment line.Results: Eighty-two patients with mCRC treated in first- (82.9%) or second- (17.1%) line chemotherapy were included. Patients with a high (>10 ng/mL) versus low (≤0.1 ng/mL) ctDNA concentration at C₀ had a shorter overall survival (OS; 6.8 vs. 33.4 months: adjusted HR, 5.64; 95% CI, 2.5-12.6; P < 0.0001). By analyzing the evolution of the ctDNA concentration between C₀ and C₂ or C₁ (C_2or1), we classified the patients in two groups (named "good" or "bad ctDNA responders"). In multivariate analysis, patients belonging to the group called "good ctDNA responder" (n = 58) versus "bad ctDNA responder" (n = 15) had a better objective response rate (P < 0.001), and a longer median progression-free survival (8.5 vs. 2.4 months: HR, 0.19; 95% CI, 0.09-0.40; P < 0.0001) and OS (27.1 vs. 11.2 months: HR, 0.25; 95% CI, 0.11-0.57; P < 0.001).Conclusions: This study suggests that early change in ctDNA concentration is a marker of therapeutic efficacy in patients with mCRC. Clin Cancer Res; 23(18); 5416-25. ©2017 AACR.

Association of prognostic value of primary tumor location in stage III colon cancer with RAS and BRAF mutational status

3515

Background: Recent data suggest that the anatomic site of colon primary tumor may be an important factor in the interpretation of molecular markers with clinical outcome in metastatic colon cancer (CC) patients (pts). We assessed here the prognostic value of primary location in fully resected stage III CC pts and its relationship to MSI, RAS and BRAF mutational status. Methods: Pts enrolled in the PETACC-8 trial were analyzed. We categorized tumor site as located proximal (left-sided) or distal (right-sided) to the splenic flexure. The association between tumor location and disease free survival (DFS), survival after relapse (SAR) and overall survival (OS) were assessed by Cox models and adjusted for clinical and pathological features, MSI, BRAF and RASmutation status. The outcome of pts receiving FOLFOX or FOLFOX and cetuximab in the adjuvant setting were also determined according to tumor site. Results: Among the 1869 pts with full molecular data available, 755 (40%) had a right-sided tumor, 164 (10%) were MSI, 942 (50%) were mutated for RAS and 212 (11%) were mutated for BRAF. Right-sided tumor was not prognostic for DFS in the whole population but was associated to a shorter SAR (HR: 1.54 [1.23 - 1.93], p = 0.001) and OS (HR: 1.25 [1.02 - 1.54], p = 0.03). Same results were observed for MSS and for MSI pts. However, when looking at pts mutated for RAS or BRAF(MUT) and those double wild type (DWT) for those mutations, we found that right-sided tumors, when compared to left-sided tumors, was associated with a worst DFS in DWT patients (HR:1.39[1.01-1.92], p = 0.04) and a better DFS in MUT patients (HR:0.77[0.63-0.95], p = 0.01). These results were found independently of the treatment received and no beneficial effect of cetuximab on DFS or OS was observed in left-sided tumors. Conclusions: In the whole study population of stage III CC pts, though right-sided tumor location influences OS as previously reported, it does not seem to influence DFS but only SAR, when disease becomes metastatic. Interestingly, sidedness seems to influence DFS when splitting the population in MUT or DWT for RAS and BRAF, with a worst DFS for right-sided tumors in DWT and a worst DFS for left-sided tumors in RAS or BRAF mutants. Clinical trial information: 2005-003463-23.

Adjuvant FOLFOX +/- cetuximab in full RAS and BRAF wildtype stage III colon cancer patients

Background

RAS mutations have been shown to confer resistance to anti- epidermal growth factor receptor (EGFR) treatment. We analysed the results of the PETACC8 trial (cetuximab + FOLFOX vs FOLFOX) in full RAS and BRAF wildtype (WT) patients (pts) with resected stage III colon cancer.

Patients and methods

Exons 2, 3 and 4 of KRAS and NRAS, and BRAF exons 11 and 15, were sequenced using the Ampliseq colon-lung cancer panel version 2, in PETACC8 trial pts who consented to translational research. The impact of cetuximab on time to recurrence (TTR), disease-free survival (DFS) and overall survival (OS) was investigated in pts with tumours harbouring RAS and BRAF WT, and RAS mutations. The prognostic value of each individual mutation was also tested.

Results

Among the 2559 pts analysed, 745 pts (29%) were known to have KRAS exon 2 mutations and 163 pts (6.4%) the BRAF V600E mutation. Of the remaining 1651 pts, 1054 were assessed by NGS, showing that a further 227 pts (21%) had KRAS exon 2, 3, 4 or NRAS exon 2, 3, 4 mutations, and that 46 pts (4.4%) had a newly diagnosed BRAF mutation. Cetuximab added to FOLFOX did not significantly improve TTR, DFS or OS in pts with RAS WT or RAS and BRAF WT tumours (HR 0.77-1.03, all P > 0.05). Cetuximab addition was not either significantly deleterious in RAS mutant pts or in pts with rare RAS or BRAF mutations. In the overall trial population, NRAS and KRAS codon 61 mutations were the only rare mutations with the same pejorative prognostic value as KRAS exon 2 or BRAF V600E mutations.

Conclusion

Though not significant, the clinically relevant 0.76 adjusted HR observed for DFS in favour of adding cetuximab to FOLFOX, in full RAS and BRAF WT stage III colon cancer pts, may justify a new randomized controlled trial testing EGFR inhibitors in this setting.

Clinical trial number

This is an ancillary study of the PETACC8 trial: EUDRACT 2005-003463-23.

A Study of Hypermethylated Circulating Tumor DNA as a Universal Colorectal Cancer Biomarker

BACKGROUND

Circulating tumor DNA (ctDNA) has emerged as a good candidate for tracking tumor dynamics in different cancer types, potentially avoiding repeated tumor biopsies. Many different genes can be mutated within a tumor, complicating procedures for tumor monitoring, even with highly sensitive next-generation sequencing (NGS) strategies. Droplet-based digital PCR (dPCR) is a highly sensitive and quantitative procedure, allowing detection of very low amounts of circulating tumor genetic material, but can be limited in the total number of target loci monitored.

METHODS

We analyzed hypermethylation of 3 genes, by use of droplet-based dPCR in different stages of colorectal cancer (CRC), to identify universal markers for tumor follow-up.

RESULTS

Hypermethylation of WIF1 (WNT inhibitory factor 1) and NPY (neuropeptide Y) genes was significantly higher in tumor tissue compared to normal tissue, independently of tumor stage. All tumor tissues appeared positive for one of the 2 markers. Methylated ctDNA (MetctDNA) was detected in 80% of metastatic CRC and 45% of localized CRC. For samples with detectable mutations in ctDNA, MetctDNA and mutant ctDNA (MutctDNA) fractions were correlated. During follow-up of different stage CRC patients, MetctDNA changes allowed monitoring of tumor evolution.

CONCLUSIONS

These results indicate that MetctDNA could be used as a universal surrogate marker for tumor follow-up in CRC patients, and monitoring MetctDNA by droplet-based dPCR could avoid the need for monitoring mutations.

Fine-mapping of two differentiated thyroid carcinoma susceptibility loci at 9q22.33 and 14q13.3 detects novel candidate functional SNPs in Europeans from metropolitan France and Melanesians from New Caledonia

Incidence of differentiated thyroid carcinoma varies considerably between countries and ethnic groups, with particularly high incidence rates in Melanesians of New Caledonia. Differentiated thyroid cancer (DTC) has a familial relative risk higher than other cancers, highlighting the contribution of inherited factors to the disease. Recently, genome-wide association studies (GWAS) identified several DTC susceptibility loci. The most robust associations were reported at loci 9q22 (rs965513 and rs1867277) and 14q13 (rs944289 and rs116909734). In this study, we performed a fine-mapping study of the two gene regions among Europeans and Melanesians from Metropolitan France and New Caledonia. We examined 81 single nucleotide polymorphisms (SNPs) at 9q22 and 561 SNPs at 14q13 in Europeans (625 cases/776 controls) and in Melanesians (244 cases/189 controls). The association with the four SNPs previously identified in GWAS was replicated in Europeans while only rs944289 was replicated in Melanesians. Among Europeans, we found that the two SNPs previously reported at 9q22 were not independently associated to DTC and that rs965513 was the predominant signal; at 14q13, we showed that the haplotype rs944289[C]-rs116909374[C]-rs999460[T] was significantly associated with DTC risk and that the association with rs116909374 differed by smoking status (p-interaction = 0.03). Among Melanesians, a new independent signal was observed at 14q13 for rs1755774 which is strongly correlated to rs2787423; this latter is potentially a functional variant. Significant interactions with parity (p < 0.05) and body mass index were observed for rs1755774 and rs2787423. This study contributed to a better characterization of the DTC loci 9q22 and 14q13 in Europeans and in Melanesians and has identified novel variants to be prioritized for further functional studies.

Hypermethylated circulating DNA detection using picoliter droplet-based PCR in colorectal cancer

622

Background: Circulating tumor DNA (ctDNA) is thoroughly investigated as a surrogate biomarker of tumor follow-up, in different cancer types, such as colorectal cancer (CRC). Droplet-based digital PCR (ddPCR) is a highly sensitive and also quantitative method for detection of very low amount of ctDNA. Since many different genes can be mutated within a specific tumor type and also wide mutation spectrum can occur within a specific gene, procedures for ctDNA monitoring can be time consuming and need to be improved for a routinely use. To overcome these drawbacks, we characterized the methylation status of 3 genes frequently hypermethylated in CRC to identify universal markers for tumor follow-up. Methods: The characterization of the methylated status of the WIF, NPY and PENK genes in the tumor DNA was performed in 56 CRC of different stages and 45 corresponding plasma samples using droplet-based dPCR, after DNA bisulfite conversion. A two-panels assay (with albumin as a reference) was developed. Methylation level of these 3 genes in tumor tissues was compared to corresponding normal tissues (n = 22) and plasma samples (MetctDNA). To validate, plasma samples of additional 91 patients were analyzed for the presence of ctDNA both by the characterization of KRAS, BRAF, TP53 and PIK3CA mutations (MutctDNA) and of MetctDNA, at various stages of their follow-up, and 9 of them had MetctDNA assessment during treatment follow-up. Results: All tumor samples were positive for WIF1 and/or NPY markers. Hypermethylation of these two genes was significantly higher in tumor tissue compared to normal, independently of the tumor stage (p < 0.0001). MetctDNA could be detected in 75% of metastatic CRC patients and 24% of localized CRC patients (stage 1 to 3). MetctDNA and MutctDNA fractions were strongly correlated (R2 > 0.9, p < 0.0001). During follow-up, MetctDNA levels changes allowed monitoring of tumor evolution in different stages CRC patients. Conclusions: These results indicate that determination of MetctDNA by droplet-based dPCR can reach same efficiency than MutctDNA for ctDNA assessment, using only 2 markers, and thus could be considered as a universal surrogate marker of tumor follow-up in CRC patients.

Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer

Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.

Identification and characterization of novel associations in the CASP8/ALS2CR12 region on chromosome 2 with breast cancer risk

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.