Search

Many studies have focused on investigating deviations from additive interaction of two dichotomous risk factors on a binary outcome. There is, however, a gap in the literature with respect to interactions on the additive scale of >2 risk factors. In this paper, we present an approach for examining deviations from additive interaction among three or more binary exposures. The relative excess risk due to interaction (RERI) is used as measure of additive interaction. First, we concentrate on three risk factors–we propose to decompose the total RERI to: the RERI owned to the joint presence of all three risk factors and the RERI of any two risk factors, given that the third is absent. We then extend this approach, to >3 binary risk factors. For illustration, we use a sample from data from the Greek EPIC cohort and we investigate the association with overall mortality of Mediterranean diet, body mass index, and smoking. Our formulae enable better interpretability of any evidence for deviations from additivity owned to more than two risk factors and provide simple ways of communicating such results from a public health perspective by attributing any excess relative risk to specific combinations of these factors. Abbreviations: BMI: Body Mass Index; ERR: excess relative risk; EPIC: European Prospective Investigation into Cancer and nutrition; MD: Mediterranean diet; RERI: relative excess risk due to interaction; RR: relative risk; TotRERI: total relative excess risk due to interaction. © 2020 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.

Background: Although lifestyle factors have been studied in relation to individual non-communicable diseases (NCDs), their association with development of a subsequent NCD, defined as multimorbidity, has been scarcely investigated. The aim of this study was to investigate associations between five lifestyle factors and incident multimorbidity of cancer and cardiometabolic diseases. Methods: In this prospective cohort study, 291,778 participants (64% women) from seven European countries, mostly aged 43 to 58 years and free of cancer, cardiovascular disease (CVD), and type 2 diabetes (T2D) at recruitment, were included. Incident multimorbidity of cancer and cardiometabolic diseases was defined as developing subsequently two diseases including first cancer at any site, CVD, and T2D in an individual. Multi-state modelling based on Cox regression was used to compute hazard ratios (HR) and 95% confidence intervals (95% CI) of developing cancer, CVD, or T2D, and subsequent transitions to multimorbidity, in relation to body mass index (BMI), smoking status, alcohol intake, physical activity, adherence to the Mediterranean diet, and their combination as a healthy lifestyle index (HLI) score. Cumulative incidence functions (CIFs) were estimated to compute 10-year absolute risks for transitions from healthy to cancer at any site, CVD (both fatal and non-fatal), or T2D, and to subsequent multimorbidity after each of the three NCDs. Results: During a median follow-up of 11 years, 1910 men and 1334 women developed multimorbidity of cancer and cardiometabolic diseases. A higher HLI, reflecting healthy lifestyles, was strongly inversely associated with multimorbidity, with hazard ratios per 3-unit increment of 0.75 (95% CI, 0.71 to 0.81), 0.84 (0.79 to 0.90), and 0.82 (0.77 to 0.88) after cancer, CVD, and T2D, respectively. After T2D, the 10-year absolute risks of multimorbidity were 40% and 25% for men and women, respectively, with unhealthy lifestyle, and 30% and 18% for men and women with healthy lifestyles. Conclusion: Pre-diagnostic healthy lifestyle behaviours were strongly inversely associated with the risk of cancer and cardiometabolic diseases, and with the prognosis of these diseases by reducing risk of multimorbidity. © 2020 The Author(s).

Purpose: Coffee and tea constituents have shown several anti-carcinogenic activities in cellular and animal studies, including against thyroid cancer (TC). However, epidemiological evidence is still limited and inconsistent. Therefore, we aimed to investigate this association in a large prospective study. Methods: The study was conducted in the EPIC (European Prospective Investigation into Cancer and Nutrition) cohort, which included 476,108 adult men and women. Coffee and tea intakes were assessed through validated country-specific dietary questionnaires. Results: During a mean follow-up of 14 years, 748 first incident differentiated TC cases (including 601 papillary and 109 follicular TC) were identified. Coffee consumption (per 100 mL/day) was not associated either with total differentiated TC risk (HRcalibrated 1.00, 95% CI 0.97–1.04) or with the risk of TC subtypes. Tea consumption (per 100 mL/day) was not associated with the risk of total differentiated TC (HRcalibrated 0.98, 95% CI 0.95–1.02) and papillary tumor (HRcalibrated 0.99, 95% CI 0.95–1.03), whereas an inverse association was found with follicular tumor risk (HRcalibrated 0.90, 95% CI 0.81–0.99), but this association was based on a sub-analysis with a small number of cancer cases. Conclusions: In this large prospective study, coffee and tea consumptions were not associated with TC risk. © 2018, Springer-Verlag GmbH Germany, part of Springer Nature.

Background: Coffee consumption in relation to female hormone-related cancers has been investigated but metaanalyses regarding breast and ovarian cancer include studies published up to 2012 with inconsistent results for ovarian cancer. Methods: We conducted two updated meta-analyses of studies published up to June 2016 to quantify the association of coffee intake with breast and ovarian cancer risk with random effects models. We used the dataset developed by the International Agency for Research on Cancer Working Group for Monograph 116 meeting (May 2016). We additionally performed a PubMed search in June 2016. Results: Summary relative risks (RRs) (95% confidence intervals (CI)) for the study-specific highest vs. lowest coffee consumption were for breast and ovarian cancer respectively: 0.97 (0.93–1.00, Ι2 5.5%, 40 studies, 76,728 cases) and 1.03 (0.93–1.14, Ι2 31.9%, 31 studies, 13,111 cases). For decaffeinated coffee the corresponding RRs were: 1.00 (0.93-1.08, I2 32.2%, 13 studies) and 0.83 (0.71-0.96, I2 about 0%, 9 studies). The association of coffee with ovarian cancer risk was higher among publications before (RR=1.37, 1.12–1.69) compared to after 2000 (RR=0.96, 0.86-1.06). Conclusion: Our meta-analyses provide strong, quantitative evidence that coffee consumption is not related to breast cancer risk and appears to be unrelated to ovarian cancer risk. © 2019, Indian Journal of Public Health Research and Development. All rights reserved.

Background & Aims: To date, evidence on the association between physical activity and risk of hepatobiliary cancers has been inconclusive. We examined this association in the European Prospective Investigation into Cancer and Nutrition cohort (EPIC). Methods: We identified 275 hepatocellular carcinoma (HCC) cases, 93 intrahepatic bile duct cancers (IHBCs), and 164 non-gallbladder extrahepatic bile duct cancers (NGBCs) among 467,336 EPIC participants (median follow-up 14.9 years). We estimated cause-specific hazard ratios (HRs) for total physical activity and vigorous physical activity and performed mediation analysis and secondary analyses to assess robustness to confounding (e.g. due to hepatitis virus infection). Results: In the EPIC cohort, the multivariable-adjusted HR of HCC was 0.55 (95% CI 0.38–0.80) comparing active and inactive individuals. Regarding vigorous physical activity, for those reporting >2 hours/week compared to those with no vigorous activity, the HR for HCC was 0.50 (95% CI 0.33–0.76). Estimates were similar in sensitivity analyses for confounding. Total and vigorous physical activity were unrelated to IHBC and NGBC. In mediation analysis, waist circumference explained about 40% and body mass index 30% of the overall association of total physical activity and HCC. Conclusions: These findings suggest an inverse association between physical activity and risk of HCC, which is potentially mediated by obesity. Lay summary: In a pan-European study of 467,336 men and women, we found that physical activity is associated with a reduced risk of developing liver cancers over the next decade. This risk was independent of other liver cancer risk factors, and did not vary by age, gender, smoking status, body weight, and alcohol consumption. © 2019 European Association for the Study of the Liver

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

Gallstones, a common gastrointestinal condition, can lead to several digestive complications and can result in inflammation. Risk factors for gallstones include obesity, diabetes, smoking and physical inactivity, all of which are known risk factors for colorectal cancer (CRC), as is inflammation. However, it is unclear whether gallstones are a risk factor for CRC. We examined the association between history of gallstones and CRC in the European Prospective Investigation into Cancer and Nutrition (EPIC) study, a prospective cohort of over half a million participants from ten European countries. History of gallstones was assessed at baseline using a self-reported questionnaire. The analytic cohort included 334,986 participants; a history of gallstones was reported by 3,917 men and 19,836 women, and incident CRC was diagnosed among 1,832 men and 2,178 women (mean follow-up: 13.6 years). Hazard ratios (HR) and 95% confidence intervals (CI) for the association between gallstones and CRC were estimated using Cox proportional hazards regression models, stratified by sex, study centre and age at recruitment. The models were adjusted for body mass index, diabetes, alcohol intake and physical activity. A positive, marginally significant association was detected between gallstones and CRC among women in multivariable analyses (HR = 1.14, 95%CI 0.99–1.31, p = 0.077). The relationship between gallstones and CRC among men was inverse but not significant (HR = 0.81, 95%CI 0.63–1.04, p = 0.10). Additional adjustment for details of reproductive history or waist circumference yielded minimal changes to the observed associations. Further research is required to confirm the nature of the association between gallstones and CRC by sex. © 2018 UICC

Background: The impact of cisplatin use on long-term survival of unselected patients with advanced urinary tract cancer (aUTC) has not been adequately investigated. We used a multinational database to study long-term survival and the impact of treatment type in unselected patients with aUTC. Materials and Methods: A total of 1,333 patients with aUTC (cT4bN0M0, cTanyN+M0, cTanyNanyM+), transitional-cell, squamous, or adenocarcinoma histology who received systemic chemotherapy and had available survival data were selected. Long-term survival was defined as alive at 3 years following initiation of first-line chemotherapy. Conditional overall survival (COS) analysis was employed to study change in prognosis given time survived from initiation of first-line chemotherapy. Results: Median follow-up was 31.7 months. The combination of cisplatin use and cisplatin eligibility accurately predicted long-term survival. Eligible patients treated with cisplatin conferred a 31.6% probability of 3-year survival (95% confidence interval [CI]: 25.1–38.3), and 2-year COS for patients surviving 3 years after initiation of cisplatin-based chemotherapy was 83% (95% CI: 59.7–93.5). The respective probabilities for patients who were ineligible for cisplatin or not treated with cisplatin despite eligibility were 14% (95% CI: 10.8–17.6) and 49.3% (95% CI: 28.2–67.4). Two-year COS remained significantly different between these two groups up to 3 years after chemotherapy initiation. Conclusion: Cisplatin-based therapy was associated with the highest likelihood of long-term survival in patients with aUTC and should be used in patients who fulfill the established eligibility criteria. Novel therapies are necessary to increase long-term survival in cisplatin-ineligible patients. Implications for Practice: Long-term, disease-free survival is possible in one in four eligible-for-cisplatin patients with advanced urinary tract cancer (aUTC) treated with cisplatin-based combination chemotherapy. Therefore, deviations from eligibility criteria should be avoided. Consolidation surgery should be considered in responders. These data provide benchmarks for the study of novel therapies in aUTC. © AlphaMed Press 2019

We investigated the predictive and prognostic significance of tumoral messenger RNA levels of vascular endothelial growth factor A (VEGFS) splice variants in metastatic colorectal cancer (mCRC) patients treated with bevacizumab. VEGFA145b had negative predictive significance predominantly in those patients with right-sided primary tumors. All VEGFAxxxb variants were negative prognosticators for patients with right-sided mCRC, whereas VEGFA165b was of favorable prognostic significance in patients with left-sided tumors. © 2019 Elsevier Inc. Background: Alternative splicing of vascular endothelial growth factor A (VEGFA) results in VEGFAxxxb antiangiogenic isoforms that fail to activate angiogenesis. Bevacizumab, widely used in patients with metastatic colorectal cancer (CRC), binds both VEGFA and VEGFAxxxb isoforms. Patients and Methods: Formalin-fixed, paraffin-embedded primary tumors from metastatic CRC patients treated with first-line FOLFIRI (leucovorin, 5-fluorouracil, irinotecan, and oxaliplatin) + bevacizumab (n = 285) or FOLFIRI only (n = 75) were collected. The relative expression of VEGFA121a, 121b, 145a, 145b, 165a, and 165b was assessed with custom TaqMan-MGB assays and quantitative PCR. Results: At a median follow-up of 101.5 months, left-sided primary CRC was a favorable prognosticator (median survival, 29.2 vs. 18.2 months; P = .015). Positive high VEGFA145b was an unfavorable factor for progression-free survival (PFS; hazard ratio [HR] = 1.66; 95% confidence interval [CI], 1.13-2.44; P = .009) in patients who received FOLFIRI + bevacizumab, without prognostic significance in FOLFIRI-only patients (HR = 0.70; 95% CI, 0.34-1.44; P = .33). The adverse effect on PFS of 145b was more pronounced in patients with right-sided colon cancer (HR = 2.62; 95% CI, 1.35-5.12; P = .005), especially in those who received bevacizumab (HR = 2.85; 95% CI, 1.31-6.21; P = .008). In patients with right-sided colon primary tumors, isoform 121b correlated with inferior PFS (HR = 1.73; 95% CI, 0.94-3.18; P = .076) and overall survival (OS; HR = 2.0; 95% CI, 1.08-3.72; P = .028). In patients with left-sided primary tumors, positive high 165b correlated with superior PFS (HR = 0.76; 95% CI, 0.59-0.99; P = .044) and OS (HR = 0.68; 95% CI, 0.52-0.90; P = .006). At multivariate analysis, right-sided primary tumor was associated with inferior PFS (HR = 1.28; 95% CI, 1.00-1.64), while 145b consistently retained predictive significance for lack of benefit in PFS with bevacizumab (HR = 1.71; 95% CI, 1.16-2.53). Multivariate analysis for OS showed that VEGFA165b expression was favorable in patients with left-sided but unfavorable in patients with right-sided primary tumors (Pinteraction < .001). Conclusion: The antiangiogenic isoform VEGFA145b messenger RNA may predict resistance to bevacizumab. Differences in biological relevance and prognostic significance of various VEGFA isoforms were found for right- versus left-sided primary tumors. © 2019 Elsevier Inc.

Aim: To further quantify the association between physical activity (PA) after breast cancer diagnosis and all-cause mortality, breast cancer mortality and/or breast cancer recurrence. Methods and results: PubMed was searched until November 2017 for observational studies investigating any type of PA in association with total mortality, breast cancer mortality and/or breast cancer recurrence among women with breast cancer diagnosis. Pooled hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated using random-effects models for highest versus lowest categories of PA. Ten studies were included in the meta-analysis. During an average follow-up ranging from 3.5 to 12.7 years there were 23,041 breast cancer survivors, 2,522 deaths from all causes, 841 deaths from breast cancer and 1,398 recurrences/remissions. Compared to women in the lowest recreational PA level (lowest quintile/quartile), women in the highest level had a lower risk of all-cause mortality (HR = 0.58, 95% CIs: 0.45–0.75; 8 studies), of death from breast cancer (HR = 0.60, 95% CIs 0.36–0.99; 5 studies) and a lower, albeit non-significant, risk of recurrence (HR = 0.79, 95% CIs 0.60–1.05; 5 studies). There was evidence of heterogeneity between the studies evaluating recreational PA and total mortality (Ι2 = 52.4%) and even higher for breast cancer mortality (Ι2 = 77.7%) or recurrence (Ι2 = 66.4%). Conclusion: Highest recreational PA after breast cancer diagnosis was associated with lower all-cause and breast cancer mortality. This finding probably reflects the favorable impact of PA on cardiovascular mortality, and a possible favorable role on breast cancer survival, though reverse causation cannot be excluded. © 2019 Elsevier Ltd

There is a growing evidence of the protective role of the Mediterranean diet (MD) on cancer. However, no prospective study has yet investigated its influence on lymphoma. We evaluated the association between adherence to the MD and risk of lymphoma and its subtypes in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. The analysis included 476,160 participants, recruited from 10 European countries between 1991 and 2001. Adherence to the MD was estimated through the adapted relative MD (arMED) score excluding alcohol. Cox proportional hazards regression models were used while adjusting for potential confounders. During an average follow-up of 13.9 years, 3,136 lymphomas (135 Hodgkin lymphoma [HL], 2,606 non-HL and 395 lymphoma not otherwise specified) were identified. Overall, a 1-unit increase in the arMED score was associated with a 2% lower risk of lymphoma (95% CI: 0.97; 1.00, p-trend = 0.03) while a statistically nonsignificant inverse association between a high versus low arMED score and risk of lymphoma was observed (hazard ratio [HR]: 0.91 [95% CI 0.80; 1.03], p-trend = 0.12). Analyses by lymphoma subtype did not reveal any statistically significant associations. Albeit with small numbers of cases (N = 135), a suggestive inverse association was found for HL (HR 1-unit increase = 0.93 [95% CI: 0.86; 1.01], p-trend = 0.07). However, the study may have lacked statistical power to detect small effect sizes for lymphoma subtype. Our findings suggest that an increasing arMED score was inversely related to the risk of overall lymphoma in EPIC but not by subtypes. Further large prospective studies are warranted to confirm these findings. © 2018 UICC

Recent evidence suggested a weak relationship between alcohol consumption and pancreatic cancer (PC) risk. In our study, the association between lifetime and baseline alcohol intakes and the risk of PC was evaluated, including the type of alcoholic beverages and potential interaction with smoking. Within the European Prospective Investigation into Cancer and Nutrition (EPIC) study, 1,283 incident PC (57% women) were diagnosed from 476,106 cancer-free participants, followed up for 14 years. Amounts of lifetime and baseline alcohol were estimated through lifestyle and dietary questionnaires, respectively. Cox proportional hazard models with age as primary time variable were used to estimate PC hazard ratios (HR) and their 95% confidence interval (CI). Alcohol intake was positively associated with PC risk in men. Associations were mainly driven by extreme alcohol levels, with HRs comparing heavy drinkers (>60 g/day) to the reference category (0.1–4.9 g/day) equal to 1.77 (95% CI: 1.06, 2.95) and 1.63 (95% CI: 1.16, 2.29) for lifetime and baseline alcohol, respectively. Baseline alcohol intakes from beer (>40 g/day) and spirits/liquors (>10 g/day) showed HRs equal to 1.58 (95% CI: 1.07, 2.34) and 1.41 (95% CI: 1.03, 1.94), respectively, compared to the reference category (0.1–2.9 g/day). In women, HR estimates did not reach statistically significance. The alcohol and PC risk association was not modified by smoking status. Findings from a large prospective study suggest that baseline and lifetime alcohol intakes were positively associated with PC risk, with more apparent risk estimates for beer and spirits/liquors than wine intake. © 2018 IARC/WHO

Background: Few published breast cancer (BC) risk prediction models consider the heterogeneity of predictor variables between estrogen-receptor positive (ER+) and negative (ER-) tumors. Using data from two large cohorts, we examined whether modeling this heterogeneity could improve prediction. Methods: We built two models, for ER+ (ModelER+) and ER- tumors (ModelER-), respectively, in 281,330 women (51% postmenopausal at recruitment) from the European Prospective Investigation into Cancer and Nutrition cohort. Discrimination (C-statistic) and calibration (the agreement between predicted and observed tumor risks) were assessed both internally and externally in 82,319 postmenopausal women from the Women's Health Initiative study. We performed decision curve analysis to compare ModelER+ and the Gail model (ModelGail) regarding their applicability in risk assessment for chemoprevention. Results: Parity, number of full-term pregnancies, age at first full-term pregnancy and body height were only associated with ER+ tumors. Menopausal status, age at menarche and at menopause, hormone replacement therapy, postmenopausal body mass index, and alcohol intake were homogeneously associated with ER+ and ER- tumors. Internal validation yielded a C-statistic of 0.64 for ModelER+ and 0.59 for ModelER-. External validation reduced the C-statistic of ModelER+ (0.59) and ModelGail (0.57). In external evaluation of calibration, ModelER+ outperformed the ModelGail: the former led to a 9% overestimation of the risk of ER+ tumors, while the latter yielded a 22% underestimation of the overall BC risk. Compared with the treat-all strategy, ModelER+ produced equal or higher net benefits irrespective of the benefit-to-harm ratio of chemoprevention, while ModelGail did not produce higher net benefits unless the benefit-to-harm ratio was below 50. The clinical applicability, i.e. the area defined by the net benefit curve and the treat-all and treat-none strategies, was 12.7 × 10- 6 for ModelER+ and 3.0 × 10- 6 for ModelGail. Conclusions: Modeling heterogeneous epidemiological risk factors might yield little improvement in BC risk prediction. Nevertheless, a model specifically predictive of ER+ tumor risk could be more applicable than an omnibus model in risk assessment for chemoprevention. © 2018 The Author(s).

Assessing diet through dietary patterns has become popular in diet-disease investigations due to the appealing simplicity of expressing several dietary exposures through one variable. Current evidence suggests that there may exist a priori and a posteriori dietary patterns that are consistently associated with reduced all-cause, as well as site-specific cancer risk across different populations. Typical examples include the Mediterranean dietary pattern, the Healthy Eating Index, a number of "healthy" a posteriori dietary patterns, and others. Despite their apparent differences, by inspecting their components it seems that multiple dietary patterns reflect core constituents of a healthy diet. Ongoing research is targeted to: (a) identify the common features underlying the beneficial-for-cancer-prevention dietary patterns, (b) decompose the estimated associations of dietary patterns with cancer risk to the individual associations of their components, and (c) construct dietary patterns for site-specific cancer risk prediction. Results of these and other initiatives could be helpful for public health policy makers and responsible bodies to (a) better formulate relevant messages for informing people about the components of cancer-preventing diets, and (b) promote small changes in eating patterns that can lower cancer risk and improve cancer outcomes. © 2018 Macmillan Publishers Limited, part of Springer Nature.

Background: Methylation measures quantified by microarray techniques can be affected by systematic variation due to the technical processing of samples, which may compromise the accuracy of the measurement process and contribute to bias the estimate of the association under investigation. The quantification of the contribution of the systematic source of variation is challenging in datasets characterized by hundreds of thousands of features. In this study, we introduce a method previously developed for the analysis of metabolomics data to evaluate the performance of existing normalizing techniques to correct for unwanted variation. Illumina Infinium HumanMethylation450K was used to acquire methylation levels in over 421,000 CpG sites for 902 study participants of a case-control study on breast cancer nested within the EPIC cohort. The principal component partial R-square (PC-PR2) analysis was used to identify and quantify the variability attributable to potential systematic sources of variation. Three correcting techniques, namely ComBat, surrogate variables analysis (SVA) and a linear regression model to compute residuals were applied. The impact of each correcting method on the association between smoking status and DNA methylation levels was evaluated, and results were compared with findings from a large meta-analysis. Results: A sizeable proportion of systematic variability due to variables expressing 'batch' and 'sample position' within 'chip' was identified, with values of the partial R2 statistics equal to 9.5 and 11.4% of total variation, respectively. After application of ComBat or the residuals' methods, the contribution was 1.3 and 0.2%, respectively. The SVA technique resulted in a reduced variability due to 'batch' (1.3%) and 'sample position' (0.6%), and in a diminished variability attributable to 'chip' within a batch (0.9%). After ComBat or the residuals' corrections, a larger number of significant sites (k = 600 and k = 427, respectively) were associated to smoking status than the SVA correction (k = 96). Conclusions: The three correction methods removed systematic variation in DNA methylation data, as assessed by the PC-PR2, which lent itself as a useful tool to explore variability in large dimension data. SVA produced more conservative findings than ComBat in the association between smoking and DNA methylation. © 2018 The Author(s).

Several studies support an important role of angiogenesis in breast cancer growth and metastasis. The main objectives of the study were to investigate the immunohistochemical expression of vascular endothelial growth factor (VEGF) family ligands (VEGF-A and VEGF-C) and receptors (VEGFR1, VEGFR2 and VEGFR3) in breast cancer and their associations with clinicopathological parameters, cancer subtypes/subgroups and patient outcome. Formalin-fixed paraffin-embedded tumor tissue samples were collected from early-stage breast cancer patients treated with anthracycline-based chemotherapy within a randomized trial. Immunohistochemistry was performed on serial 2.5 μm thick tissue sections from tissue microarray blocks. High VEGF-A, VEGF-C, VEGFR1, VEGFR2 and VEGFR3 protein expression was observed in 11.8% (N = 87), 80.8% (N = 585), 28.1% (N = 202), 64.6% (N = 359) and 71.8% (N = 517) of the cases, respectively. Significant associations were observed among all proteins (all p-values