Association of angiogenesis and inflammation‑related gene functional polymorphisms with early‑stage breast cancer prognosis
- Authors:
- Published online on: April 7, 2020 https://doi.org/10.3892/ol.2020.11521
- Pages: 3687-3700
-
Copyright: © Korobeinikova et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Breast cancer (BC) is one of the most commonly diagnosed cancers and the leading cause of cancer death among women worldwide (1). Improved diagnostic capabilities has led to an increased rate of BC identified at an early stage (2). However, despite early diagnosis and treatment, the rate of recurrence and metastasis following radical treatment remains disappointingly high, and survival varies considerably between patients with closely matching tumour characteristics.
Inflammation and angiogenesis are the main drivers of cancer. These processes are tightly interconnected in the sense that many pro-inflammatory proteins possess proangiogenic properties and vice versa. The formation of new blood vessels in malignant tumours ensures the supply of nutrients and oxygen, hence promoting the spread of tumour cells (3). Microvessel density is a pivotal risk factor for metastasis and a predictor of poor BC prognosis (4).
The vascular endothelial growth factor A (VEGFA, or VEGF) is the best-known proangiogenic molecule. VEGFA mediates the growth of new blood vessels by binding to the endothelial cell surface receptors. It promotes endothelial cell proliferation, migration and the formation of tubular structures (5). Another mechanism of tumour neovascularization is the so-called inflammatory angiogenesis. Such pro-inflammatory cytokines as interleukin-1β (IL-1β), interlukin-1α (IL-1α) and interleukin-6 (IL-6) through a variety of signalling pathways promote endothelial cell migration and proliferation, contributing to tumour angiogenesis that facilitates the survival of cancer cells (6,7). Apart from stimulating angiogenesis, VEGFA, IL-1β1, IL-1α and IL-6 are also involved in inflammatory processes. These proteins may prevent apoptosis and promote cancer cell proliferation, differentiation, migration and metastasis.
The association between the above-mentioned proteins and BC prognosis was demonstrated by several authors. The elevated serum level of VEGFA in metastatic BC patients is linked to worse progression-free survival (PFS) and overall survival (OS) (8). Another cytokine, IL-6, induces epithelial-mesenchymal phenotype and therapeutic resistance in BC cells (9). Higher circulating levels of IL6 were observed in more advanced stages of the disease (10). It was also found that high tumour co-expression of the VEGF and IL-6 family cytokines significantly lowers the human epidermal growth factor receptor 2 (HER2) negative BC survival (11). Moreover, higher expression of pro-inflammatory IL-1β cytokine is correlated with higher BC stage and significantly worse survival (12). In addition, IL-1α acts as a pro-inflammatory molecule itself and also promotes the activity of IL-1β, resulting in an increased growth of BC cells and tumour progression (13).
Common polymorphisms in proinflammatory and proangiogenic cytokine genes may influence their coded protein production and play a role in the course of BC. The polymorphisms with proved functional activity are VEGFA (rs699947, rs833061, rs25648, rs1005230), IL-1β (rs1143634, rs16944), IL-1α (rs1800587) and IL-6 (rs1800795, rs1800797) (14–21). Candidate gene studies as well as moderate-sized genome-wide association studies (GWAS) highlight the important role of these polymorphisms in BC risk and aggressiveness (22–32), although substantial heterogeneity across studies exists. The currently available results are inconsistent in terms of different ethnicity and cancer stage; therefore, it is necessary to further investigate the role of VEGFA, IL-1β, IL-1α and IL-6 gene polymorphisms in breast carcinogenesis and cancer progression.
This paper describes a cohort study that aimed to examine the contribution of VEGFA, IL-1β, IL-1α and IL-6 gene polymorphisms to the clinicopathologic features and survival in a homogeneous group of Eastern European (specifically, Lithuanian) early-stage BC patients.
Materials and methods
Patients
The study consisted of 202 adult Lithuanian women with primary I–II stage BC. All patients were treated in the Hospital of Lithuanian University of Health Sciences Kaunas Clinics. The exclusion criteria were other malignancies, significant comorbidities and/or incomplete medical documentation. Surgery and adjuvant therapy were chosen by the clinicians, based on pathomorphologic characteristics and validated prognostic factors. The patients were followed until 30 April, 2019 (censoring date).
Candidate polymorphisms
The genes and polymorphisms known to modulate inflammation and angiogenesis were selected. The selection criteria included: i) functional single nucleotide polymorphisms (SNPs) in the VEGFA, IL-1β, IL-1α and IL-6 genes predicting alterations in the protein level; ii) SNP relevant to outcomes in other settings; and iii) SNP with a minor allele frequency greater than 15% in the study population. We selected nine SNPs: The VEGFA gene rs699947, rs833061, rs25648, and rs1005230; the IL-1β gene rs1143634 and rs16944; the IL-1α gene rs1800587; and the IL-6 gene rs1800795 and rs1800797.
Assay methods
Peripheral blood samples from the study population were collected in 2009–2017. For genomic DNA extraction from peripheral blood leukocytes, a commercially available DNA extraction kit (i.e., Thermo Fisher Scientific Baltics, Lithuania) was used. The DNA was stored at −20°C prior to usage.
Genotyping of the selected polymorphisms was performed at the Dr. K. Janusauskas Laboratory of Genetics of the Institute of Biology Systems and Genetic Research of Lithuanian University of Health Sciences Kaunas Clinics. The SNPs of the target genes were estimated by using TaqMan SNP Genotyping Assays (C_8311602_10, C_1647381_10, C_791476_10, C_8311612_10, C_1839697_20, C_1839695_20, C_9546517_10, C_1839943_10, C_1839943_10; Applied Biosystems; Thermo Fisher Scientific, Inc.). The polymerase chain reaction was performed in a reaction volume of 25 µl containing template DNA (2 ng), 2X TaqMan Universal Master Mix II, no UNG (Applied Biosystems; Thermo Fisher Scientific, Inc.)-12,5 µl, 20X TaqMan SNP Genotyping Assay stock (initial stock of 40X or 80X was diluted to get 20X working stock)-1,25 µl. The final volume of 25 µl was adjusted by adding nuclease free ddH2O. Finally, the 2 µl of DNA was added from each sample. For negative control, nuclease free ddH2O was used instead of patient DNA, while for positive control, the DNA of the known genotype was used. Each sample genotyping was repeated twice for accuracy.
The Applied Biosystems 7900HT Real-Time Polymerase Chain Reaction System (Applied Biosystems; Thermo Fisher Scientific, Inc.) was used for SNP detection. The cycling program started from heating up to 95°C for 10 min followed by 40 cycles (at 95°C for 15 sec and at 60°C for 1 min). Finally, allelic discrimination was done by using the SDS 2.3 software provided by Applied Biosystems; Thermo Fisher Scientific, Inc.
Study design
A prospective cohort study was conducted at the Hospital of Lithuanian University of Health Sciences. For the case selection, information on primarily BC patients was retrieved from the hospital's Pathology Department. The patients who fulfilled the inclusion criteria and signed the informed consent document (approved by the Kaunas Regional Ethics Committee for Biomedical research; Protocol number BE-2-10) were enrolled in the study and their peripheral blood samples were obtained. The characteristics of clinical and pathological features and the course of the disease were obtained for all study subjects. The date of histological BC verification was time zero in the survival analysis. The endpoints of interest were disease-free survival (DFS), metastasis-free survival (MFS) and OS. We checked for associations of SNPs with the known BC prognostic factors and survival endpoints. The guidelines for the reporting of tumour marker prognostic studies were applied while conducting the study (33,34).
Statistical analysis
The allele frequency distributions of the investigated SNPs were compared with the European population data from the 1000 Genomes project phase 3 database (35). For each SNP a Hardy-Weinberg equilibrium was assessed by using Pearson's chi square and Fisher's exact tests. The Haploview v4.1 software was used to check for the linkage disequilibrium between SNPs (36). The VEGFA and IL-6 haplotypes were inferred from the tested SNPs by Bayesian methods as implemented in the Phase software (v2.1; Department of Statistics, University of Washington, Seattle, WA, USA) (37,38). The SNPs were analysed under genotype, allelic and haplotype (for IL-6 and VEGFA SNPs) models. The associations of polymorphisms with clinicopathologic variables were evaluated by Pearson's Chi-square or Fisher's exact test. The Bonferroni-corrected α level was used in the association analysis for multiple comparisons. Moreover, the Cox proportional hazards model was used to estimate the prognostic factors for DFS, MFS and OS. In addition, multivariate analysis was used to determine the interdependency of genotypes and other known prognostic factors, such as age, tumour differentiation grade, tumour size, lymph node status, oestrogen receptor status, progesterone receptor status and human epidermal growth factor receptor 2 (HER2) status. Hazard ratios (HR) and their 95% confidence intervals (CI) were recorded for each tested marker. Finally, the Kaplan-Meier analysis with the log-rank test was applied to compare the survival of the patients with different genotypes. Statistical significance was set at 5% (P<0.05). Statistical analysis was performed using SPSS for Windows v20.0 (Released 2011; IBM Corp).
Results
Sample characteristics
A total of 202 Lithuanian women with early-stage BC were included in the current analysis. The frequency distributions of clinical and tumour biological factors are shown in Table I. For all study participants, primary treatments included surgery (100%), chemotherapy (77%), hormone therapy (71%), trastuzumab (19%) and radiation therapy (97%).
All the patients were genotyped for a panel of nine SNPs: The VEGFA gene rs699947, rs833061, rs25648 and rs1005230; the IL-1β gene rs1143634 and rs16944; the IL-1α gene rs1800587; the IL-6 gene rs1800795 and rs1800797. The genotypes were found to be in Hardy-Weinberg equilibrium in all the nine SNPs. A strong linkage disequilibrium between four VEGFA and two IL-6 polymorphisms was confirmed (Fig. 1). Our cohort had similar allele distribution to that of the 1000 Genomes project phase 3 for European population. The genotype and allele frequency data is presented in Table II.
Table II.Allele and genotype frequencies of analysed single nucleotide polymorphisms in the study population and in European population data from the 1000 Genomes Project Phase 3 database. |
Inferential analysis
The data on associations between the analysed polymorphisms and clinicopathologic tumour features is shown in Tables SI–SVI. In the single-locus analysis, the genotype model revealed a significant link between the IL-6 rs1800797 genotypes and the oestrogen receptor (P=0.005) and the progesterone receptor (P=0.007) status (Bonferoni adjusted; significant p value <0.008). The allelic model showed that the A allele of this SNP is associated with positive oestrogen receptors (OR 2.23; 95% CI: 1.19–4.08; P=0.014). Patients carrying the IL-6 rs 1800797 A allele were also predisposed to higher rates of HER2 negative BC (OR 2.21; 95% CI: 1.07–4.57; P=0.042). Additionally, another IL-6 polymorphism, rs1800795, in both genotype and allelic models was linked to oestrogen receptor positive BC. Specifically, 72% of patients carrying the IL-6 rs1800795 C allele had oestrogen receptor positive disease, compared to 58% of non-carriers (OR 1.92; 95% CI: 1.04–3.57; P=0.04).
Linkage disequilibrium analysis showed a high degree of disequilibrium between the two IL-6 SNPs (r2=0.89), meaning that the associations were not independent. As such, haplotype analysis was performed to further explore the relationship between IL-6 variations and other prognostic factors. Phasing revealed three possible IL-6 (rs1800797, rs1800795) haplotypes: GG (51.7%), AC (45.5%) and GC (2.8%). It was further demonstrated that the AC haplotype was positively associated with oestrogen receptor positive BC (HR 2.13; 95% CI 1.14–3.98; P=0.018; Table SII).
By analysing the associations between the IL-1α SNP and clinicopathologic factors (Table SIII), we found a link between the IL-1α rs1800587 SNP C allele and larger primary tumour size (comparing tumours sized <2 vs. 2–5 cm) (OR 4.91; 95% CI: 1.09–22.01; P=0.022).
Other analysed polymorphisms of IL-1β and VEGFA revealed no associations with the evaluated prognostic factors in neither genotype nor allelic models (Tables SIV–SVI). Due to the strong linkage disequilibrium, four VEGFA polymorphisms were also analysed in the haplotype model. Three main haplotypes were identified: CCTC (45.5%), ATCC (33.4%) and ATCT (17.6%), as well as several rare variants. None of the VEGFA haplotypes was associated with clinicopathologic prognostic variables.
Survival analysis
In the mean follow-up time of 67 months (range 28–202), progression of the disease was observed for 33 patients. Of those who progressed, 28 had distant metastases. Twenty-two patients with progressive disease died, all due to cancer. The Kaplan-Meier survival analysis showed that the IL-1α rs1800587 SNP is associated with early-stage breast cancer DFS, MFS and OS. In particular, patients homozygous for the C allele (CC) had worse survival rates than patients homozygous and heterozygous for the A allele (CT and TT) (Fig. 2A-C).
The univariate Cox regression analysis (presented in Tables III–V) revealed that L-1α rs1800587 CC genotype carriers had 2.48 times higher risk of disease recurrence (95% CI: 1.19–5.11; P=0.014), 3.12 times higher risk of metastasis (95% CI: 1.37–7.10; P=0.007) and 2.63 times higher risk of death (95% CI: 1.07–6.46; P=0.035) than the carriers of CT and TT genotypes (Table III).
By analysing the SNPs of the IL-6 gene, we found that the rs1800797 GG genotype was a negative prognostic factor for DFS and MFS. Furthermore, IL-6 rs1800797 GG and IL-6 rs1800795 GG genotype carriers displayed a shorter OS (Table IV). In a haplotype model, the patients who inherited the GG/GG diplotype (AC haplotype non-carriers) had a higher risk of disease recurrence, metastasis and death (Table IV; Fig. 3A and B).
In a multivariate Cox regression model including age at diagnosis, tumour size, lymph node status, tumour differentiation grade, oestrogen receptor, progesterone receptor and HER2 status, the IL-1α 1800587 CC genotype remained a significant predictor of poor DFS, MFS and OS. Furthermore, the IL-6 GG/GG diplotype (non-carrying the AC haplotype) was an independent negative prognostic factor for MFS and OS (Table VI). Finally, other polymorphisms and VEGFA haplotypes were not associated with any of the survival endpoints.
Discussion
In the present study, we investigated the associations between nine functional SNPs in four cytokine genes (i.e., VEGFA, IL-1α, IL-1β and IL-6) and the clinicopathologic profiles and survival rates in a group of Lithuanian women with early-stage BC. We found an association between the IL-6 SNPs and the oestrogen receptor positive, progesterone receptor positive and HER2 negative status, and a link between IL-1α SNP and larger primary tumour size. Furthermore, we confirmed a negative prognostic value of the IL-6 rs1800797-rs1800795 GG/GG diplotype on MFS and OS and of the IL-1α rs1800587 CC genotype on DFS, MFS and OS in a highly homogeneous group of patients.
Several authors have demonstrated that carrying the IL-6 rs1800795 G or IL-6 rs1800795 G alleles is associated with higher IL-6 protein production (19–21). Specifically, IL-6 is a cytokine which is considered a prognostic marker as well as a potential therapeutic target for BC patients. This cytokine acts through several pathways, regulating the proliferation, apoptosis and metabolism of BC cells. However, the most important role of IL-6 in breast carcinogenesis is its potential to induce breast metastasis by enhancing angiogenesis and tumour cell migration (39). According to previous studies, IL-6 rs1800797 and rs1800795 polymorphisms appear to be biologically important; however, the data on their clinical importance is still heterogeneous.
The association of IL-6 polymorphisms with BC survival was analysed in several studies (Table VII). Snoussi et al (32) and DeMichele et al (27) analysed patients with non-metastatic stage I–III BC. In both studies the IL-6 rs1800795 GG genotype was associated with decreased DFS and OS. In DeMichele et al (27) the presence of at least one copy of the haplotype rs1800797G-rs1800796G-373(10A/11T)-rs1800795G was associated with worse DFS. Snoussi et al (32) also found a significant link between IL-6 rs1800797 GG polymorphism and decreased DFS and OS. However, in DeMichele et al (28), where only locally advanced stage III patients with more than 10 positive lymph nodes were studied, the results were not replicated. Specifically, there was no impact of IL-6 polymorphisms on DFS and OS in the study population. In the same study, only an unplanned sub-analysis showed an association of the rs1800795 GG and rs1800797 GG genotypes with lower DFS but not OS in the oestrogen receptor positive patient subgroup. Our study is the first to confirm the prognostic value of IL-6 polymorphisms (namely, the IL-6 GG/GG diplotype) on MFS and OS in early-stage BC. Therefore, IL-6 GG/GG may potentially be used for the selection of patients who need intensified adjuvant BC treatment.
Table VII.Studies evaluating the assocaitions of IL-6 rs1800795, IL-6 rs1800797 and IL-1α rs1800587 single nucleotide polymorphisms with breast cancer survival. |
The above-mentioned studies by Snoussi et al (32) and DeMichele et al (27,28) did not analyse MFS as a survival endpoint. To the best of our knowledge, the only study that evaluated the association between IL-6 polymorphisms and BC metastasis was a case-control study conducted by Abana et al (26). The authors demonstrated the significant link between the IL-6 rs1800795 GG genotype and the development of BC metastases with an OR of 1.52. However, due to case-control design, we have no data on MFS differences. In our study IL-6 GG/GG predicted MFS and OS, but not DFS. Taking into account the critical role of IL-6 in the development of BC metastasis and our study results, we propose the hypothesis that IL-6 polymorphisms may predict the development of metastasis rather than loco-regional relapse.
It is important to mention that, in contrast with the findings of the above-mentioned studies, Markkula et al (29), who analysed a Swedish cohort of patients with any primary stage non-metastatic BC, demonstrated that the carriers of the C allele in rs1800795 SNP had a higher risk of early BC events. The authors, however, conducted no DFS, MFS or OS analyses. Furthermore, none of the GWAS studies showed a statistically significant association between IL-6 polymorphisms and BC survival (23,24). The contrasting findings in IL-6 polymorphism studies may be due to sampling errors or differences in patient ethnicity. Nevertheless, most of the experimental data and the possible biological pathway support our findings.
IL-6 also acts as a regulator of oestrogen synthesis and may modulate the tumour cell growth related to the hormone receptor status (40). Hormone-sensitive cells exhibit a higher response to IL-6, while ER-negative cells are suppressed by IL-6 (41). However, there is no clear mechanism of association between IL-6 SNPs and the oestrogen receptor positive, progesterone receptor positive and HER2 negative status. Therefore, further studies are needed.
As far as IL-1α SNP is concerned, Um et al (18) demonstrated that the IL-1α rs1800587 CC genotype is associated with higher transcriptional activity of the IL−1α gene. Overexpression of the IL-1α promotes tumour invasiveness and metastasis by inducing the expression of angiogenic genes and growth factors (12). A large meta-analysis performed by Xia et al (30) showed that the IL-1α rs1800587 C allele and CC genotype are associated with increased cancer risk in general. The studies which investigated the associations between IL-1α rs1800587 polymorphism and BC survival are presented in Table VII. Grimm et al (31) analysed 262 Caucasian BC patients and found that the IL-1α rs1800758 C allele in the univariate model is associated with OS; however, the multivariate model failed to repeat the association. In our study, the multivariate survival analysis confirmed the statistically significant impact of the CC genotype on DFS, MFS and OS. Contrasting data is provided by Snoussi et al (32), who analysed North African BC patients and found that the rs1800587 homozygous TT genotype showed a significant association with reduced DFS and OS rate. However, the allele and genotype frequency of rs1800587 SNP in African population differs substantially from that in the population with European ancestry. This difference may be the cause of the observed inconsistency between these studies. We demonstrated that in an Eastern-European population the IL-1α rs1800758 C allele is associated with more aggressive local disease (i.e., larger tumour size) and for the first time we proposed that the IL-1α rs1800758 CC genotype is an independent negative prognostic biomarker for early-stage BC. As the allele frequencies of the IL-1α rs1800758 SNP in our study correspond to those of the 1000 Genomes project phase 3 database for European population, these findings could potentially be replicated for European population, but larger confirmatory studies are warranted.
Although experimental data suggests that VEGFA and IL-1β play a role in BC, our findings, along with the results from several other studies (22–24), do not suggest that several common polymorphisms in these genes are associated with BC clinical and morphological variables and BC survival rates.
Potential limitations of our study include a limited sample size, lack of access to tumour and/or tumour stromal tissue (which would have allowed for the assessment of the presence of SNPs in those tissues), the risk of other confounders, possible gene-gene and gene-environment interactions, and non-random sampling. We also acknowledge that VEGFA, IL-1α, IL-1β and IL-6 measurements were not available in the current study. However, our study supports the relevance of functional IL-6 and IL-1α germline polymorphisms to BC prognosis. Further investigations, preferable on larger cohorts with different ethnic origins, are needed to confirm the results of the current study.
In conclusion, we found an association between the IL−1α rs1800587 C allele and larger primary tumour size. The IL-6 rs1800797 A allele, GA genotype, IL-6 rs1800795 C allele and IL-6 (rs1800797-re1800795) AC diplotype were linked to hormonal receptor positive BC. Additionally, IL-6 rs1800797 A was associated with HER2 negative status. The multivariate IL-1α rs1800587 CC genotype was confirmed as an independent negative prognostic factor for DFS, MFS and OS, and the IL6 GG/GG diplotype for MFS and OS in early-stage BC patients. Our findings confirm the hypothesis that functional SNPs in angiogenesis- and inflammation-associated genes are associated with early-stage BC prognosis in Lithuanian population.
Supplementary Material
Supporting Data
Acknowledgements
Not applicable.
Funding
The present study was partially supported by the Lithuanian University of Health Sciences Science Foundation.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Authors' contributions
EK, RU and EJu designed the research study. EK, RU, RI, LP, VR, EJa and EJu performed the research. EK and RU analyzed the data. RU and RI contributed essential reagents or tools. EK wrote the manuscript.
Ethics approval and consent to participate
The present study was approved by Kaunas Regional Ethics Committee for Biomedical Research (protocol no. BE-2-10). Informed consent for participation in the study and use of tissue was obtained from all participants.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Glossary
Abbreviations
Abbreviations:
BC |
breast cancer |
SNP |
single nucleotide polymorphism |
VEGFA |
vascular endothelial growth factor A |
IL-6 |
interleukin-6 |
IL-1α |
interleukin-1α |
IL-1β |
interleukin-1β |
A |
adenine |
G |
guanine |
T |
thymine |
C |
cytosine |
DFS |
disease-free survival |
MFS |
metastasis-free survival |
OS |
overall survival |
HR |
hazard ratio |
OR |
odds ratio |
CI |
confidence interval |
HER2 |
human epidermal growth factor receptor 2 |
GWAS |
genome-wide association study |
References
IARC: World Cancer Report 2014. Stewart BW and Wild CP: IARC; Lyon: 2014 | |
Cancer Trends Progress Report. National Cancer Institute, NIH, DHHS; Bethesda, MD: 2019, https://progressreport.cancer.gov | |
Zhao Y and Adjei AA: Targeting angiogenesis in cancer therapy: Moving beyond vascular endothelial growth factor. Oncologist. 20:660–673. 2015. View Article : Google Scholar : PubMed/NCBI | |
Bielenberg DR and Zetter BR: The contribution of angiogenesis to the process of metastasis. Cancer J. 21:267–273. 2015. View Article : Google Scholar : PubMed/NCBI | |
Cartland SP, Genner SW, Zahoor A and Kayurma MM: Comparative evaluation of TRAIL, FGF-2 and VEGF-A-induced angiogenesis in vitro and in vivo. Int J Mol Sci. 17:E20252016. View Article : Google Scholar : PubMed/NCBI | |
De Palma M, Biziato D and Petrova TV: Microenvironmental regulation of tumour angiogenesis. Nat Rev Cancer. 17:457–474. 2017. View Article : Google Scholar : PubMed/NCBI | |
Marech I, Leporini C, Ammendola M, Porcelli M, Gadaleta CD, Russo E, De Sarro G and Ranieri G: Classical and non-classical proangiogenic factors as a target of antiangiogenic yherapy in tumor microenvironment. Cancer Lett. 380:216–226. 2016. View Article : Google Scholar : PubMed/NCBI | |
Banys-Paluchoothtwski M, Witzel I, Riethdorf S, Pantel K, Rack B, Janni W, Fasching PA, Aktas B, Kasimir-Bauer S, Hartkopf A, et al: The clinical relevance of serum vascular endothelial growth factor (VEGF) in correlation to circulating tumor cells and other serum biomarkers in patients with metastatic breast cancer. Breast Cancer Res Treat. 172:93–104. 2018. View Article : Google Scholar : PubMed/NCBI | |
Gyamfi J, Eom M, Koo JS and Choi J: Multifaceted roles of interleukin-6 in adipocyte-breast cancer cell interaction. Transl Oncol. 11:275–285. 2018. View Article : Google Scholar : PubMed/NCBI | |
Knüpfer H and Preiss R: Significance of interleukin-6 (IL-6) in breast cancer (Review). Breast Cancer Res Treat. 102:129–135. 2007. View Article : Google Scholar : PubMed/NCBI | |
Tawara K, Scott H, Emathinger J, Ide A, Fox R, Greiner D, LaJoie D, Hedeen D, Nandakumar M, Oler AJ, et al: Co-expression of VEGF and IL-6 family cytokines is associated with decreased survival in HER2 negative breast cancer patients: Subtype-specific IL-6 family cytokine-mediated VEGF secretion. Transl Oncol. 12:245–255. 2019. View Article : Google Scholar : PubMed/NCBI | |
Litmanovich A, Khazim K and Cohen I: The role of interleukin-1 in the pathogenesis of cancer and its potential as a therapeutic target in clinical practice. Oncol Ther. 6:109–127. 2018. View Article : Google Scholar | |
Balasubramanian SP, Azmy IA, Higham SE, Wilson AG, Cross SS, Cox A, Brown NJ and Reed MW: Interleukin gene polymorphisms and breast cancer: A case control study and systematic literature review. BMC Cancer. 6:1882006. View Article : Google Scholar : PubMed/NCBI | |
Koukourakis MI, Papazoglou D, Giatromanolaki A, Bougioukas G, Maltezos E and Siviridis E: VEGF gene sequence variation defines VEGF gene expression status and angiogenic activity in non-small cell lung cancer. Lung Cancer. 46:293–298. 2004. View Article : Google Scholar : PubMed/NCBI | |
Watson CJ, Webb NJ, Bottomley MJ and Brenchley PE: Identification of polymorphisms within the vascular endothelial growth factor (VEGF) gene: Correlation with variation in VEGF protein production. Cytokine. 12:1232–1235. 2000. View Article : Google Scholar : PubMed/NCBI | |
Ray D, Mishra M, Ralph S, Read I, Davies R and Brenchley P: Association of the VEGF gene with proliferative diabetic retinopathy but not proteinuria in diabetes. Diabetes. 53:861–864. 2004. View Article : Google Scholar : PubMed/NCBI | |
Hall SK, Perregaux DG, Gabel CA, Woodworth T, Durham LK, Huizinga TW, Breedveld FC and Seymour AB: Correlation of polymorphic variation in the promoter region of the interleukin-1 beta gene with secretion of interleukin-1 beta protein. Arthritis Rheum. 50:1976–1983. 2004. View Article : Google Scholar : PubMed/NCBI | |
Um JY, Rim HK, Kim SJ, Kim HL and Hong SH: Functional polymorphism of IL-1 alpha and its potential role in obesity in humans and mice. PLoS One. 6:e295242011. View Article : Google Scholar : PubMed/NCBI | |
Rana BK, Flatt SW, Health DD, Pakiz B, Quintana EL, Natarajan L and Rock CL: The IL-6 gene promoter SNP and plasma IL-6 in response to diet intervention. Nutrients. 9:E5522017. View Article : Google Scholar : PubMed/NCBI | |
Boeta-Lopez K, Duran J, Elizondo D, Gonzales E, Rentfro A, Schwarzbach AE and Nair S: Association of interleukin-6 polymorphisms with obesity or metabolic traits in young Mexican-Americans. Obes Sci Pract. 4:85–96. 2017. View Article : Google Scholar : PubMed/NCBI | |
Sampaio AM, Balseiro SC, Silva MR, Alarcão A, d'Aguiar MJ, Ferreira T and Carvalho L: Association between IL-4 and IL-6 expression variants and gastric cancer among portuguese population. GE Port J Gastroenterol. 22:143–152. 2015. View Article : Google Scholar : PubMed/NCBI | |
Langsenlehner U, Hofmann G, Renner W, Gerger A, Krenn-Pilko S, Thurner EM, Krippl P and Langsenlehner T: Association of vascular endothelial growth factor--a gene polymorphisms and haplotypes with breast cancer metastases. Acta Oncol. 54:368–376. 2015. View Article : Google Scholar : PubMed/NCBI | |
Escala-Garcia M, Guo Q, Dörk T, Canisius S, Keeman R, Dennis J, Beesley J, Lecarpentier J, Bolla MK, Wang Q, et al: Genome-wide association study of germline variants and breast cancer-specific mortality. Br J Cancer. 120:647–657. 2019. View Article : Google Scholar : PubMed/NCBI | |
Khan S, Fagerholm R, Kadalayil L, Tapper W, Aittomäki K, Liu J, Blomqvist C, Eccles D and Nevanlinna H: Meta-analysis of three genome-wide association studies identifies two loci that predict survival and treatment outcome in breast cancer. Oncotarget. 9:4249–4257. 2017.PubMed/NCBI | |
Peng X, Shi J, Sun W, Ruan X, Guo Y, Zhao L, Wang J and Li B: Genetic polymorphisms of IL-6 promoter in cancer susceptibility and prognosis: A meta-analysis. Oncotarget. 9:12351–12364. 2018. View Article : Google Scholar : PubMed/NCBI | |
Abana CO, Bingham BS, Cho JH, Graves AJ, Koyama T, Pilarski RT, Chakravarthy AB and Xia F: IL-6 variant is associated with metastasis in breast cancer patients. PLoS One. 12:e01817252017. View Article : Google Scholar : PubMed/NCBI | |
DeMichele A, Martin AM, Mick R, Gor P, Wray L, Klein-Cabral M, Athanasiadis G, Colligan T, Stadtmauer E and Weber B: Interleukin-6 −174G->C polymorphism is associated with improved outcome in high-risk breast cancer. Cancer Res. 63:8051–8056. 2003.PubMed/NCBI | |
DeMichele A, Gray R, Horn M, Chen J, Aplenc R, Vaughan WP and Tallman MS: Host genetic variants in the interleukin-6 promoter predict poor outcome in patients with estrogen receptor-positive, node-positive breast cancer. Cancer Res. 69:4184–4191. 2009. View Article : Google Scholar : PubMed/NCBI | |
Markkula A, Simonsson M, Ingvar C, Rose C and Jernström H: IL6 genotype, tumour ER-status, and treatment predicted disease-free survival in a prospective breast cancer cohort. BMC Cancer. 14:7592014. View Article : Google Scholar : PubMed/NCBI | |
Xia H, Chen Y, Meng J and Liang C: Effect of polymorphisms on IL1A to cancer susceptibility: Evidence based on 34,016 subjects. Artif Cells Nanomed Biotechnol. 47:3138–3152. 2019. View Article : Google Scholar : PubMed/NCBI | |
Grimm C, Kantelhardt E, Heinze G, Polterauer S, Zeillinger R, Kölbl H, Reinthaller A and Hefler L: The prognostic value of four interleukin-1 gene polymorphisms in Caucasian women with breast cancer: A multicenter study. BMC Cancer. 9:782009. View Article : Google Scholar : PubMed/NCBI | |
Snoussi K, Strosberg AD, Bouaouina N, Ben Ahmed S and Chouchane L: Genetic variation in pro-inflammatory cytokines (interleukin-1beta, interleukin-1alpha and interleukin-6) associated with the aggressive forms, survival, and relapse prediction of breast carcinoma. Eur Cytokine Netw. 16:253–260. 2005.PubMed/NCBI | |
Hayes DF, Ethier S and Mippman ME: New guidelines for reporting of tumor marker studies in breast cancer research and treatment: REMARK. Breast Cancer Res Treat. 100:237–238. 2006. View Article : Google Scholar : PubMed/NCBI | |
McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M and Clark GM; Statistics Subcommittee of the NCI-EORTC Working Group on Cancer Diagnostics, : Reporting recommendations for tumor marker prognostic studies (REMARK). J Natl Cancer Inst. 97:1180–1184. 2005. View Article : Google Scholar : PubMed/NCBI | |
1000 Genomes Project Consortium, ; Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA and Abecasis GR: A global reference for human genetic variation. Nature. 526:68–74. 2015. View Article : Google Scholar : PubMed/NCBI | |
Barrett JC, Fry B, Maller J and Daly MJ: Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics. 21:263–265. 2005. View Article : Google Scholar : PubMed/NCBI | |
Stephens M, Smith NJ and Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 68:978–989. 2001. View Article : Google Scholar : PubMed/NCBI | |
Stephens M and Scheet P: Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am J Hum Genet. 76:449–462. 2005. View Article : Google Scholar : PubMed/NCBI | |
Kumari N, Dwarakanath BS, Das A and Bhatt AN: Role of interleukin-6 in cancer progression and therapeutic resistance. Tumour Biol. 37:11553–11572. 2016. View Article : Google Scholar : PubMed/NCBI | |
Dethlefsen C, Højfeldt G and Hojman P: The role of intratumoral and systemic IL-6 in breast cancer. Breast Cancer Res Treat. 138:657–664. 2013. View Article : Google Scholar : PubMed/NCBI | |
Masjedi A, Hashemi V, Hojjat-Farsangi M, Ghalamfarsa G, Azizi G, Yousefi M and Jadidi-Niaragh F: The significant role of IL-6 and its signaling pathway in the immunopathogenesis and treatment of breast cancer. Biomed Pharmacother. 108:1415–1424. 2018. View Article : Google Scholar : PubMed/NCBI |