A pan-cancer analysis of the microbiome in metastatic cancer.

Microbial communities are resident to multiple niches of the human body and are important modulators of the host immune system and responses to anticancer therapies. Recent studies have shown that complex microbial communities are present within primary tumors. To investigate the presence and relevance of the microbiome in metastases, we integrated mapping and assembly-based metagenomics, genomics, transcriptomics, and clinical data of 4,160 metastatic tumor biopsies. We identified organ-specific tropisms of microbes, enrichments of anaerobic bacteria in hypoxic tumors, associations between microbial diversity and tumor-infiltrating neutrophils, and the association of Fusobacterium with resistance to immune checkpoint blockade (ICB) in lung cancer. Furthermore, longitudinal tumor sampling revealed temporal evolution of the microbial communities and identified bacteria depleted upon ICB. Together, we generated a pan-cancer resource of the metastatic tumor microbiome that may contribute to advancing treatment strategies.

Long-term outcomes of young, node-negative, chemotherapy-naïve, triple-negative breast cancer patients according to BRCA1 status.

BACKGROUND: Due to the abundant usage of chemotherapy in young triple-negative breast cancer (TNBC) patients, the unbiased prognostic value of BRCA1-related biomarkers in this population remains unclear. In addition, whether BRCA1-related biomarkers modify the well-established prognostic value of stromal tumor-infiltrating lymphocytes (sTILs) is unknown. This study aimed to compare the outcomes of young, node-negative, chemotherapy-naïve TNBC patients according to BRCA1 status, taking sTILs into account. METHODS: We included 485 Dutch women diagnosed with node-negative TNBC under age 40 between 1989 and 2000. During this period, these women were considered low-risk and did not receive chemotherapy. BRCA1 status, including pathogenic germline BRCA1 mutation (gBRCA1m), somatic BRCA1 mutation (sBRCA1m), and tumor BRCA1 promoter methylation (BRCA1-PM), was assessed using DNA from formalin-fixed paraffin-embedded tissue. sTILs were assessed according to the international guideline. Patients' outcomes were compared using Cox regression and competing risk models. RESULTS: Among the 399 patients with BRCA1 status, 26.3% had a gBRCA1m, 5.3% had a sBRCA1m, 36.6% had tumor BRCA1-PM, and 31.8% had BRCA1-non-altered tumors. Compared to BRCA1-non-alteration, gBRCA1m was associated with worse overall survival (OS) from the fourth year after diagnosis (adjusted HR, 2.11; 95% CI, 1.18-3.75), and this association attenuated after adjustment for second primary tumors. Every 10% sTIL increment was associated with 16% higher OS (adjusted HR, 0.84; 95% CI, 0.78-0.90) in gBRCA1m, sBRCA1m, or BRCA1-non-altered patients and 31% higher OS in tumor BRCA1-PM patients. Among the 66 patients with tumor BRCA1-PM and ≥ 50% sTILs, we observed excellent 15-year OS (97.0%; 95% CI, 92.9-100%). Conversely, among the 61 patients with gBRCA1m and < 50% sTILs, we observed poor 15-year OS (50.8%; 95% CI, 39.7-65.0%). Furthermore, gBRCA1m was associated with higher (adjusted subdistribution HR, 4.04; 95% CI, 2.29-7.13) and tumor BRCA1-PM with lower (adjusted subdistribution HR, 0.42; 95% CI, 0.19-0.95) incidence of second primary tumors, compared to BRCA1-non-alteration. CONCLUSIONS: Although both gBRCA1m and tumor BRCA1-PM alter BRCA1 gene transcription, they are associated with different outcomes in young, node-negative, chemotherapy-naïve TNBC patients. By combining sTILs and BRCA1 status for risk classification, we were able to identify potential subgroups in this population to intensify and optimize adjuvant treatment.

B-cells and regulatory T-cells in the microenvironment of HER2+ breast cancer are associated with decreased survival: a real-world analysis of women with HER2+ metastatic breast cancer.

BACKGROUND: Despite major improvements in treatment of HER2-positive metastatic breast cancer (MBC), only few patients achieve complete remission and remain progression free for a prolonged time. The tumor immune microenvironment plays an important role in the response to treatment in HER2-positive breast cancer and could contain valuable prognostic information. Detailed information on the cancer-immune cell interactions in HER2-positive MBC is however still lacking. By characterizing the tumor immune microenvironment in patients with HER2-positive MBC, we aimed to get a better understanding why overall survival (OS) differs so widely and which alternative treatment approaches may improve outcome. METHODS: We included all patients with HER2-positive MBC who were treated with trastuzumab-based palliative therapy in the Netherlands Cancer Institute between 2000 and 2014 and for whom pre-treatment tissue from the primary tumor or from metastases was available. Infiltrating immune cells and their spatial relationships to one another and to tumor cells were characterized by immunohistochemistry and multiplex immunofluorescence. We also evaluated immune signatures and other key pathways using next-generation RNA-sequencing data. With nine years median follow-up from initial diagnosis of MBC, we investigated the association between tumor and immune characteristics and outcome. RESULTS: A total of 124 patients with 147 samples were included and evaluated. The different technologies showed high correlations between each other. T-cells were less prevalent in metastases compared to primary tumors, whereas B-cells and regulatory T-cells (Tregs) were comparable between primary tumors and metastases. Stromal tumor-infiltrating lymphocytes in general were not associated with OS. The infiltration of B-cells and Tregs in the primary tumor was associated with unfavorable OS. Four signatures classifying the extracellular matrix of primary tumors showed differential survival in the population as a whole. CONCLUSIONS: In a real-world cohort of 124 patients with HER2-positive MBC, B-cells, and Tregs in primary tumors are associated with unfavorable survival. With this paper, we provide a comprehensive insight in the tumor immune microenvironment that could guide further research into development of novel immunomodulatory strategies.

Characterization of the Tumor Microenvironment of De Novo Oligometastatic Breast Cancer in a Nationwide Cohort.

PURPOSE: Oligometastatic breast cancer (OMBC) has a more favorable outcome than widespread metastatic breast cancer. Some patients with OMBC achieve long-term remission if treated with multimodality therapy, including systemic and locally ablative therapies. However, not all patients with OMBC benefit from such treatment, while all experience toxicity. To explore biomarkers identifying patients with OMBC and potential long-term survival, we compared tumor-immune characteristics of patients with OMBC and long-term versus shorter-term survival. MATERIALS AND METHODS: We collected tumor tissue of 97 patients with de novo OMBC (≤5 metastases) via the Dutch nationwide cancer and pathology registries using a case-control design. Long-term survivors (LTS) were defined as patients alive ≥10 years since OMBC diagnosis. Fifty-five LTS and 42 shorter-term survivors (STS) were included. Median follow-up was 15 years (IQR, 14-16). Tumor characteristics and infiltrating immune cells were assessed by immunohistochemistry and next-generation RNA-sequencing. Association of the resulting 52 biomarkers with long-term survival was assessed using logistic regression. Associations with survival within LTS were assessed using Cox-proportional hazards modeling. P values were adjusted for multiple hypothesis testing. RESULTS: Most patients had estrogen receptor (ER)-positive OMBC (n = 86; 89%) and 23 (24%) had human epidermal growth factor receptor 2-positive disease. ER positivity in primary tumors distinguished LTS from STS. In addition, extracellular matrix (ECM)2-low and ECM4-high distinguished between long-term and shorter-term survival. Immune levels in the primary tumor did not associate with LTS. However, within the LTS subset, higher immune levels associated with improved progression-free survival. CONCLUSION: We identified tumor and ECM features in the primary tumor of patients with de novo OMBC that were associated with long-term survival. Our data should be validated in other patients with OMBC before they can be used in clinical practice.

Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

Breast cancer genomes from CHEK2 c.1100delC mutation carriers lack somatic TP53 mutations and display a unique structural variant size distribution profile.

BACKGROUND: CHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still unclear how exactly CHEK2 c.1100delC promotes tumorigenesis. Since the mutational landscape of a tumor reflects the processes that have operated on its development, the aim of this study was to uncover the somatic genomic landscape of CHEK2-associated BC. METHODS: We sequenced primary BC (pBC) and normal genomes of 20 CHEK2 c.1100delC mutation carriers as well as their pBC transcriptomes. Including pre-existing cohorts, we exhaustively compared CHEK2 pBC genomes to those from BRCA1/2 mutation carriers, those that displayed homologous recombination deficiency (HRD) and ER- and ER+ pBCs, totaling to 574 pBC genomes. Findings were validated in 517 metastatic BC genomes subdivided into the same subgroups. Transcriptome data from 168 ER+ pBCs were used to derive a TP53-mutant gene expression signature and perform cluster analysis with CHEK2 BC transcriptomes. Finally, clinical outcome of CHEK2 c.1100delC carriers was compared with BC patients displaying somatic TP53 mutations in two well-described retrospective cohorts totaling to 942 independent pBC cases. RESULTS: BC genomes from CHEK2 mutation carriers were most similar to ER+ BC genomes and least similar to those of BRCA1/2 mutation carriers in terms of tumor mutational burden as well as mutational signatures. Moreover, CHEK2 BC genomes did not show any evidence of HRD. Somatic TP53 mutation frequency and the size distribution of structural variants (SVs), however, were different compared to ER+ BC. Interestingly, BC genomes with bi-allelic CHEK2 inactivation lacked somatic TP53 mutations and transcriptomic analysis indicated a shared biology with TP53 mutant BC. Moreover, CHEK2 BC genomes had an increased frequency of > 1 Mb deletions, inversions and tandem duplications with peaks at specific sizes. The high chromothripsis frequency among CHEK2 BC genomes appeared, however, not associated with this unique SV size distribution profile. CONCLUSIONS: CHEK2 BC genomes are most similar to ER+ BC genomes, but display unique features that may further unravel CHEK2-driven tumorigenesis. Increased insight into this mechanism could explain the shorter survival of CHEK2 mutation carriers that is likely driven by intrinsic tumor aggressiveness rather than endocrine resistance.

Immune landscape of breast tumors with low and intermediate estrogen receptor expression.

Immune checkpoint blockade (ICB) is currently approved for patients with triple-negative breast cancer (TNBC), whereas responses to ICB are also observed in a small subgroup of Estrogen Receptor (ER)-positive breast cancer. The cut-off for ER-positivity (≥1%) is based on likelihood of endocrine treatment response, but ER-positive breast cancer represents a very heterogeneous group. This raises the question whether selection based on ER-negativity should be revisited to select patients for ICB treatment in the context of clinical trials. Stromal tumor-infiltrating lymphocytes (sTILs) and other immune parameters are higher in TNBC compared to ER-positive breast cancer, but it is unknown whether lower ER levels are associated with more inflamed tumor microenvironments (TME). We collected a consecutive series of primary tumors from 173 HER2-negative breast cancer patients, enriched for tumors with ER expression between 1 and 99% and found levels of stromal TILs, CD8 + T cells, and PD-L1 positivity in breast tumors with ER 1-9% and ER 10-50% to be comparable to tumors with ER 0%. Expression of immune-related gene signatures in tumors with ER 1-9% and ER 10-50% was comparable to ER 0%, and higher than in tumors with ER 51-99% and ER 100%. Our results suggest that the immune landscape of ER low tumors (1-9%) and ER intermediate tumors (10-50%) mimic that of primary TNBC.

IFN-γ signature enables selection of neoadjuvant treatment in patients with stage III melanoma.

Neoadjuvant ipilimumab + nivolumab has demonstrated high pathologic response rates in stage III melanoma. Patients with low intra-tumoral interferon-γ (IFN-γ) signatures are less likely to benefit. We show that domatinostat (a class I histone deacetylase inhibitor) addition to anti-PD-1 + anti-CTLA-4 increased the IFN-γ response and reduced tumor growth in our murine melanoma model, rationalizing evaluation in patients. To stratify patients into IFN-γ high and low cohorts, we developed a baseline IFN-γ signature expression algorithm, which was prospectively tested in the DONIMI trial. Patients with stage III melanoma and high intra-tumoral IFN-γ scores were randomized to neoadjuvant nivolumab or nivolumab + domatinostat, while patients with low IFN-γ scores received nivolumab + domatinostat or ipilimumab + nivolumab + domatinostat. Domatinostat addition to neoadjuvant nivolumab ± ipilimumab did not delay surgery but induced unexpected severe skin toxicity, hampering domatinostat dose escalation. At studied dose levels, domatinostat addition did not increase treatment efficacy. The baseline IFN-γ score adequately differentiated patients who were likely to benefit from nivolumab alone versus patients who require other therapies.

Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry.

BACKGROUND: Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes. METHODS: We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry. RESULTS: In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C. CONCLUSIONS: Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

Uncovering the Contribution of Moderate-Penetrance Susceptibility Genes to Breast Cancer by Whole-Exome Sequencing and Targeted Enrichment Sequencing of Candidate Genes in Women of European Ancestry.

Rare variants in at least 10 genes, including BRCA1, BRCA2, PALB2, ATM, and CHEK2, are associated with increased risk of breast cancer; however, these variants, in combination with common variants identified through genome-wide association studies, explain only a fraction of the familial aggregation of the disease. To identify further susceptibility genes, we performed a two-stage whole-exome sequencing study. In the discovery stage, samples from 1528 breast cancer cases enriched for breast cancer susceptibility and 3733 geographically matched unaffected controls were sequenced. Using five different filtering and gene prioritization strategies, 198 genes were selected for further validation. These genes, and a panel of 32 known or suspected breast cancer susceptibility genes, were assessed in a validation set of 6211 cases and 6019 controls for their association with risk of breast cancer overall, and by estrogen receptor (ER) disease subtypes, using gene burden tests applied to loss-of-function and rare missense variants. Twenty genes showed nominal evidence of association (p-value < 0.05) with either overall or subtype-specific breast cancer. Our study had the statistical power to detect susceptibility genes with effect sizes similar to ATM, CHEK2, and PALB2, however, it was underpowered to identify genes in which susceptibility variants are rarer or confer smaller effect sizes. Larger sample sizes would be required in order to identify such genes.

Prognostic Value of Stromal Tumor-Infiltrating Lymphocytes in Young, Node-Negative, Triple-Negative Breast Cancer Patients Who Did Not Receive (neo)Adjuvant Systemic Therapy.

PURPOSE: Triple-negative breast cancer (TNBC) is considered aggressive, and therefore, virtually all young patients with TNBC receive (neo)adjuvant chemotherapy. Increased stromal tumor-infiltrating lymphocytes (sTILs) have been associated with a favorable prognosis in TNBC. However, whether this association holds for patients who are node-negative (N0), young (< 40 years), and chemotherapy-naïve, and thus can be used for chemotherapy de-escalation strategies, is unknown. METHODS: We selected all patients with N0 TNBC diagnosed between 1989 and 2000 from a Dutch population-based registry. Patients were age < 40 years at diagnosis and had not received (neo)adjuvant systemic therapy, as was standard practice at the time. Formalin-fixed paraffin-embedded blocks were retrieved (PALGA: Dutch Pathology Registry), and a pathology review including sTILs was performed. Patients were categorized according to sTILs (< 30%, 30%-75%, and ≥ 75%). Multivariable Cox regression was performed for overall survival, with or without sTILs as a covariate. Cumulative incidence of distant metastasis or death was analyzed in a competing risk model, with second primary tumors as competing risk. RESULTS: sTILs were scored for 441 patients. High sTILs (≥ 75%; 21%) translated into an excellent prognosis with a 15-year cumulative incidence of a distant metastasis or death of only 2.1% (95% CI, 0 to 5.0), whereas low sTILs (< 30%; 52%) had an unfavorable prognosis with a 15-year cumulative incidence of a distant metastasis or death of 38.4% (32.1 to 44.6). In addition, every 10% increment of sTILs decreased the risk of death by 19% (adjusted hazard ratio: 0.81; 95% CI, 0.76 to 0.87), which are an independent predictor adding prognostic information to standard clinicopathologic variables (χ2 = 46.7, P < .001). CONCLUSION: Chemotherapy-naïve, young patients with N0 TNBC with high sTILs (≥ 75%) have an excellent long-term prognosis. Therefore, sTILs should be considered for prospective clinical trials investigating (neo)adjuvant chemotherapy de-escalation strategies.

Circadian PERformance in breast cancer: a germline and somatic genetic study of PER3VNTR polymorphisms and gene co-expression.

Polymorphisms in the PER3 gene have been associated with several human disease phenotypes, including sleep disorders and cancer. In particular, the long allele of a variable number of tandem repeat (VNTR) polymorphism has been previously linked to an increased risk of breast cancer. Here we carried out a combined germline and somatic genetic analysis of the role of the PER3VNRT polymorphism in breast cancer. The combined data from 8284 individuals showed a non-significant trend towards increased breast cancer risk in the 5-repeat allele homozygous carriers (OR = 1.17, 95% CI: 0.97-1.42). We observed allelic imbalance at the PER3 locus in matched blood and tumor DNA samples, showing a significant retention of the long variant (risk) allele in tumor samples, and a preferential loss of the short repetition allele (p = 0.0005). Gene co-expression analysis in healthy and tumoral breast tissue samples uncovered significant associations between PER3 expression levels with those from genes which belong to several cancer-associated pathways. Finally, relapse-free survival (RFS) analysis showed that low expression levels of PER3 were linked to a significant lower RSF in luminal A (p = 3 × 10-12) but not in the rest of breast cancer subtypes.

Absence of PD-L1 expression on tumor cells in the context of an activated immune infiltrate may indicate impaired IFNγ signaling in non-small cell lung cancer.

BACKGROUND: In non-small cell lung cancer (NSCLC), PD-L1 expression on either tumor cells (TC) or both TC and tumor-infiltrating immune cells (IC) is currently the most used biomarker in cancer immunotherapy. However, the mechanisms involved in PD-L1 regulation are not fully understood. To provide better insight in these mechanisms, a multiangular analysis approach was used to combine protein and mRNA expression with several clinicopathological characteristics. PATIENTS AND METHODS: Archival tissues from 640 early stage, resected NSCLC patients were analyzed with immunohistochemistry for expression of PD-L1 and CD8 infiltration. In addition, mutational status and expression of a selection of immune genes involved in the PD-L1/PD-1 axis and T-cell response was determined. RESULTS: Tumors with high PD-L1 expression on TC or on IC represent two subsets of NSCLC with minimal overlap. We observed that PD-L1 expression on IC irrespective of expression on TC is a good marker for inflammation within tumors. In the tumors with the highest IC expression and absent TC expression an association with reduced IFNγ downstream signaling in tumor cells was observed. CONCLUSIONS: These results show that PD-L1 expression on TC and IC are both independent hallmarks of the inflamed phenotype in NSCLC, and TC-negative/IC-high tumors can also be categorized as inflamed. The lack of correlation between PD-L1 TC and IC expression in this subgroup may be caused by impaired IFNγ signaling in tumor cells. These findings may bring a better understanding of the tumor-immune system interaction and the clinical relevance of PD-L1 expression on IC irrespective of PD-L1 expression on TC.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer.

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Evaluation of the EGFR polymorphism R497K in two cohorts of neoadjuvantly treated breast cancer patients.

Pathological response of breast cancer to neoadjuvant chemotherapy (NAC) presents great variability, and new prognostic biomarkers are needed. Our aim was to evaluate the association of the epidermal growth factor receptor gene (EGFR) polymorphism R497K (rs2227983) with prognostic features and clinical outcomes of breast cancer, including the pathological response to NAC and the recurrence-free survival (RFS). Tumoral complete response (tCR) was defined by no remaining invasive cancer in the excised breast, whereas pathological complete response (pCR) was defined by no remaining invasive cancer both in the excised breast and lymph nodes. Two independent cohorts were analyzed: one from Brazil (INCA, n = 288) and one from The Netherlands (NKI-AVL, n = 255). In the INCA cohort, the variant (Lys-containing) genotypes were significantly associated with lower proportion of tCR (ORadj = 0.92; 95%CI = 0.85-0.99), whereas in the NKI-AVL cohort they were associated with tumor grade 3 (p = 0.035) and with triple-negative subtype (p = 0.032), but not with clinical outcomes. Such distinct prognostic associations may have arisen due to different neoadjuvant protocols (p < 0.001), or to lower age at diagnosis (p < 0.001) and higher proportion of tumor grade 3 (p = 0.018) at the NKI-AVL cohort. Moreover, NKI-AVL patients achieved better proportion of pCR (21.2% vs 8.3%, p < 0.001) and better RFS (HRadj = 0.48; 95% adjCI = 0.26-0.86) than patients from INCA. In conclusion, large scale studies comprehending different populations are needed to evaluate the impact of genome variants on breast cancer outcomes.

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.

Diabetes and Breast Cancer Subtypes.

BACKGROUND: Women with diabetes have a worse survival after breast cancer diagnosis compared to women without diabetes. This may be due to a different etiological profile, leading to the development of more aggressive breast cancer subtypes. Our aim was to investigate whether insulin and non-insulin treated women with diabetes develop specific clinicopathological breast cancer subtypes compared to women without diabetes. METHODS AND FINDINGS: This cross-sectional study included randomly selected patients with invasive breast cancer diagnosed in 2000-2010. Stratified by age at breast cancer diagnosis (≤50 and >50 years), women with diabetes were 2:1 frequency-matched on year of birth and age at breast cancer diagnosis (both in 10-year categories) to women without diabetes, to select ~300 patients with tumor tissue available. Tumor MicroArrays were stained by immunohistochemistry for estrogen and progesterone receptor (ER, PR), HER2, Ki67, CK5/6, CK14, and p63. A pathologist scored all stains and revised morphology and grade. Associations between diabetes/insulin treatment and clinicopathological subtypes were analyzed using multivariable logistic regression. Morphology and grade were not significantly different between women with diabetes (n = 211) and women without diabetes (n = 101), irrespective of menopausal status. Premenopausal women with diabetes tended to have more often PR-negative (OR = 2.44(95%CI:1.07-5.55)), HER2-negative (OR = 2.84(95%CI:1.11-7.22)), and basal-like (OR = 3.14(95%CI:1.03-9.60) tumors than the women without diabetes, with non-significantly increased frequencies of ER-negative (OR = 2.48(95%CI:0.95-6.45)) and triple negative (OR = 2.60(95%CI:0.88-7.67) tumors. After adjustment for age and BMI, the associations remained similar in size but less significant. We observed no evidence for associations of clinicopathological subtypes with diabetes in postmenopausal women, or with insulin treatment in general. CONCLUSIONS: We found no compelling evidence that women with diabetes, treated with or without insulin, develop different breast cancer subtypes than women without diabetes. However, premenopausal women with diabetes tended to develop breast tumors that do not express hormonal receptors, which are typically associated with poor prognosis.

Recurrent HOXB13 mutations in the Dutch population do not associate with increased breast cancer risk.

The HOXB13 p.G84E mutation has been firmly established as a prostate cancer susceptibility allele. Although HOXB13 also plays a role in breast tumor progression, the association of HOXB13 p.G84E with breast cancer risk is less evident. Therefore, we comprehensively interrogated the entire HOXB13 coding sequence for mutations in 1,250 non-BRCA1/2 familial breast cancer cases and 800 controls. We identified two predicted deleterious missense mutations, p.G84E and p.R217C, that were recurrent among breast cancer cases and further evaluated their association with breast cancer risk in a larger study. Taken together, 4,520 familial non-BRCA1/2 breast cancer cases and 3,127 controls were genotyped including the cases and controls of the whole gene screen. The concordance rate for the genotyping assays compared with Sanger sequencing was 100%. The prostate cancer risk allele p.G84E was identified in 18 (0.56%) of 3,187 cases and 16 (0.70%) of 2,300 controls (OR = 0.81, 95% CI = 0.41-1.59, P = 0.54). Additionally, p.R217C was identified in 10 (0.31%) of 3,208 cases and 2 (0.087%) of 2,288 controls (OR = 3.57, 95% CI = 0.76-33.57, P = 0.14). These results imply that none of the recurrent HOXB13 mutations in the Dutch population are associated with breast cancer risk, although it may be worthwhile to evaluate p.R217C in a larger study.

Age- and Tumor Subtype-Specific Breast Cancer Risk Estimates for CHEK2*1100delC Carriers.

PURPOSE: CHEK2*1100delC is a well-established breast cancer risk variant that is most prevalent in European populations; however, there are limited data on risk of breast cancer by age and tumor subtype, which limits its usefulness in breast cancer risk prediction. We aimed to generate tumor subtype- and age-specific risk estimates by using data from the Breast Cancer Association Consortium, including 44,777 patients with breast cancer and 42,997 controls from 33 studies genotyped for CHEK2*1100delC. PATIENTS AND METHODS: CHEK2*1100delC genotyping was mostly done by a custom Taqman assay. Breast cancer odds ratios (ORs) for CHEK2*1100delC carriers versus noncarriers were estimated by using logistic regression and adjusted for study (categorical) and age. Main analyses included patients with invasive breast cancer from population- and hospital-based studies. RESULTS: Proportions of heterozygous CHEK2*1100delC carriers in controls, in patients with breast cancer from population- and hospital-based studies, and in patients with breast cancer from familial- and clinical genetics center-based studies were 0.5%, 1.3%, and 3.0%, respectively. The estimated OR for invasive breast cancer was 2.26 (95%CI, 1.90 to 2.69; P = 2.3 × 10(-20)). The OR was higher for estrogen receptor (ER)-positive disease (2.55 [95%CI, 2.10 to 3.10; P = 4.9 × 10(-21)]) than it was for ER-negative disease (1.32 [95%CI, 0.93 to 1.88; P = .12]; P interaction = 9.9 × 10(-4)). The OR significantly declined with attained age for breast cancer overall (P = .001) and for ER-positive tumors (P = .001). Estimated cumulative risks for development of ER-positive and ER-negative tumors by age 80 in CHEK2*1100delC carriers were 20% and 3%, respectively, compared with 9% and 2%, respectively, in the general population of the United Kingdom. CONCLUSION: These CHEK2*1100delC breast cancer risk estimates provide a basis for incorporating CHEK2*1100delC into breast cancer risk prediction models and into guidelines for intensified screening and follow-up.