Abstract Background IRAK1 has been shown to be abnormally expressed in a set of tumors leading to tumorigenesis and progression. IRAK1 is a therapeutic target that drives breast cancer metastasis and resistance to paclitaxel. However, the exact role of IRAK1 in neoadjuvant chemotherapy (NCT) for breast cancer remains unclear. The aim of this study was to investigate the relationship between the expression of IRAK1 and the clinicopathological parameters and survival prognosis of breast cancer patients treated with NCT. Patients and methods Based on the clinical data and mRNA microarray data from 1,085 breast cancer patients in The Cancer Genome Atlas, the correlation between IRAK1 expression and clinicopathological parameters of breast cancer was analyzed. Immunohistochemistry was performed to evaluate the expression of IRAK1. The Human Protein Atlas and the String database were used to analyze the expression of IRAK1 protein and its interaction with altered neighboring proteins in breast cancer. IRAK1 alteration was analyzed in cBioPortal database. GEO enrichment of IRAK1 was performed using WEB-based Gene SeT AnaLysis Toolkit. Results The expression of IRAK1 was significantly downregulated following NCT. The decreased expression of IRAK1 following NCT was positively correlated with reduced tumor size. Finally, survival analysis confirmed that a shorter survival period was correlated to higher expression of IRAK1 both before and after NCT. Conclusion These findings advanced our understanding about the expression pattern of IRAK1 in breast cancer, especially in those patients who were treated with NCT, suggesting that IRAK1 could be used as a prognostic indicator, as well as a potential indicator for evaluating NCT efficacy for breast cancer patients. Keywords: IRAK1, breast cancer, neoadjuvant chemotherapy, bioinformatics, survival Introduction Breast cancer has the highest incidence among Chinese women, accounting for 15% of new cancers in women.[32]^1 Breast cancer is the leading cause of death in women younger than 45 years old.[33]^1 According to different expression status of receptors including estrogen receptor (ER), progesterone receptor (PR), and Her-2, breast cancer can be divided into four molecular subtypes: luminal A type, luminal B type, triple-negative type, and Her-2 overexpression type. Among them, triple-negative breast cancer has the worst prognosis and is still lacking effective prognostic indicators.[34]^2 Neoadjuvant therapy (NAT) refers to anti-tumor treatment for patients before surgery. For breast cancer, NAT includes neoadjuvant chemotherapy (NCT), neoadjuvant hormonotherapy, and neoadjuvant Her-2 targeted therapy.[35]^3 NCT has become the standard treatment for patients with locally advanced breast cancer. After NCT, some inoperable breast cancers can be treated with radical surgery. NCT makes breast-conserving surgery possible by reducing the size of primary tumors. In patients who initially need axillary lymph node dissections, NCT also increases the possibility of a separate sentinel lymph node biopsy.[36]^4^–[37]^7 Because endocrine therapy is ineffective, chemotherapy is extremely important for triple-negative breast cancer. Unfortunately, there is still a lack of effective predictors for NCT. There are four members in the IRAK family, including IRAK1, IRAK2, IRAKM, and IRAK4. The IRAK1 gene consists of 14 exons and is located on the X chromosome.[38]^8 IRAK1 is the first member of this kinase family and is a key component of the IL-1R signaling pathway. In addition, IRAK1 is a downstream effector molecule of the TLR signaling pathway. The abnormal activation of TLR signaling has an important impact on oncogenesis, allergies, septicemia, and autoimmunity.[39]^9 In recent years, reports on the relationship between IRAK1 and lung cancer, liver cancer, thyroid cancer, melanoma, and lymphoma have been frequently published.[40]^10^–[41]^13 By using bioinformatics methods, Pilarsky et al found that IRAK1 is one of the overexpressed genes commonly found in solid tumors.[42]^14 Wee et al confirmed that IRAK1 is highly expressed in triple-negative breast cancer compared to other types, and inhibition of IRAK1 expression can effectively suppress the growth and metastasis of triple-negative breast cancer.[43]^15 However, the exact role of IRAK1 in NCT for breast cancer remains unclear. In this study, we aimed to study the expression pattern of IRAK1 protein in breast cancer, and explore the relationship between IRAK1 expression and clinicopathological characteristics in breast cancer patients treated with NCT. Patients and methods Patients and tissue samples A total of 75 primary female breast cancer patients who received NCT treatment at First Hospital of China Medical University from April 2007 to February 2011 were selected for the study. According to WHO standards, all cancers were confirmed as invasive breast cancer by histopathology.[44]^16 Formalin-fixed, paraffin-embedded core biopsy specimens and resection specimens were stored in the Department of Pathology. ER, PR, HER-2, and Ki67 status were routinely examined and could be obtained from the electronic medical records. Male patients, patients with other types of breast cancer or with history of any malignancy were excluded from our study. The clinicopathological parameters of the 75 patients were summarized in [45]Table 1. Follow-up was updated by review of records and telephone calls. The median follow-up was 49 months. The date of death and the date of relapse were used to calculate disease-free survival (DFS). Formalin-fixed, paraffin-embedded breast tissues were collected from core biopsies (pre-NCT) and matched with resection tissues (post-NCT). Six patients achieved pathological complete response (pCR). This study was approved by the Institutional Ethics Committees of the First Affiliated Hospital of China Medical University. All patients received an explanation of the study aims and signed written informed consent. The study was conducted in accordance with the Declaration of Helsinki. Table 1. Clinicopathological characteristics of the patients N (%) Age  <45 42 (56)  >45 33 (44) Pre-NAC stage (based on ultrasound)  T2 38 (50.7)  T3 21 (28)  T4 16 (21.3) Post-NAC stage (based on resection pathology)  T0 5 (6.7)  T1 21 (28)  T2 31 (41.3)  T3 10 (13.3)  T4 8 (10.7) Tumor size change  Increase 4 (5.3)  Decrease 66 (88)  pCR 5 (6.7) Chemotherapy regimens  Platinum + platinum (TP) 28 (37.3)  Platinum + cyclophosphamide (TC) 15 (20)  5-fluorouracil + epirubicin + cyclophosphamide (CEF) 32 (42.7) Axillary lymph node metastasis 52 (69.3) ER positive 31 (41.3) Her-2 positive 44 (58.7) Surgery  Breast conserving 3 (4)  Mastectomy 72 (96) Follow-up  Median: 49 Range 32–130 Recurrence 29 (38.6) Death 21 (28.0) [46]Open in a new tab Abbreviations: ER, estrogen receptor; NAC, neoadjuvant chemotherapy; pCR, pathological complete response. The breast cancer mRNA expression dataset was downloaded and processed from the The Cancer Genome Atlas (TCGA) database ([47]https://cancergenome.nih.gov). Online survival analysis was performed using K-M plotter ([48]http://kmplot.com/analysis/). We only included breast cancer patients who had complete RNAseq data and complete clinical information including tumor size, lymph node status, TNM stage, ER, PR, HER-2, and overall survival (OS) data. Finally, 1,085 cases of breast cancer patients were included in our study. The RNAseq data of all 291 normal samples were used for the comparison of mRNA expression difference between tumor and normal tissues. Relevant clinical and pathological parameters were described in [49]Table 2. Table 2. Correlation between IRAK1 expression and clinicopathological parameters based on The Cancer Genome Atlas database Characteristics N IRAK1 P-value Positive Negative Molecular type <0.001  Luminal A 419 58 361  Luminal B 180 28 152  Her-2 62 15 47  Basal 140 22 118 Pathology M stage 0.048  M0 869 120 749  M1 21 17 4 Pathology N stage 0.284  N0 488 44 444  N1 355 59 296  N2 112 22 90  N3 76 15 61 Pathology T stage 0.809  T1 277 33 244  T2 603 77 526  T3 131 25 106  T4 37 15 22 [50]Open in a new tab Immunohistochemical staining Immunohistochemistry was performed with 4 µm formalin-fixed, paraffin-embedded sections. After deparaffinization and rehydration, 3% H[2]O[2] solution was used to block endogenous peroxidase at room temperature. Antigen retrieval was performed by using sodium citrate buffer, and normal serum was applied to block non-specific antibody binding (reagent B; UltraSensitive™ SP immunohistochemistry kit; Maxim Biotech Inc., Rockville, MD, USA). Slides were then incubated overnight at 4°C with the primary antibody against IRAK1 (F-4, 1:50 dilution) obtained from Santa Cruz Biotechnology Inc. ([sc-5288] Dallas, TX, USA). After overnight incubation, samples were then treated with reagent C and reagent D (UltraSensitive™ SP immunohistochemistry kit; Maxim Biotech Inc.), respectively. Finally, DAB staining was performed and the staining intensity was graded. The intensities of IRAK1 staining were scored as 0 (no staining), 1 (weak staining), 2 (medium staining), and 4 (strong staining). The results were evaluated and scored independently by at least two pathologists in a blinded manner using a semi-quantitative scoring system considering the extent and intensity. The mean percentage of positive cells was scored as 0 (0%), 1 (1%–30%), 2 (31%–50%), 3 (51%–80%), and 4 (81%–100%). The two scores were multiplied, resulting in a final immunoreactivity score ranging between 0 and 12. The six patients who achieved pCR were regarded as negative post-neoadjuvant chemotherapy (NAC). Statistical analysis All the statistical analyses were carried out using SPSS version 23.0 software and Graphpad 6.0. Wilcoxon signed-rank test was used to analyze the difference in IRAK1 expression before and after NCT. Mann–Whitney U test was carried out to assess the difference in expression change between group A and group B (divided according to the decrease in tumor size). Spearman’s correlation was used to reveal the correlation between the expression level of IRAK1 and the clinicopathological characteristics. ROC curves were drawn to find the cutoff values (giving the highest combined sensitivity and specificity) to dichotomize pre- and post-NAC expression scores of IRAK1 for the end point of DFS. DFS was estimated using Kaplan–Meier analyses. Additionally, multivariate Cox linear analysis was carried out between the expression of IRAK1 and clinicopathological features, including tumor grade, receptor status, pre-NAC tumor stage, and axillary metastasis. All P-values presented are two-sided, and P≤0.05 was considered to be statistically significant. Results The expression pattern of IRAK1 in breast cancer We first investigated the difference in mRNA expression for IRAK1 using the mRNA HiSeq expression data from TCGA. As shown in [51]Figure 1A, IRAK1 expression was upregulated in breast cancer tissues compared with the normal breast tissues. The Human Protein Atlas is dedicated to providing information on the organization and cell distribution of all 24,000 human proteins and provides free access to the public. The Swedish Knut & Alice Wallenberg Foundation, which founded the database, uses immunohistochemistry with specific antibodies to examine each protein in 48 human normal tissues, 20 tumor tissues, 47 cell lines, and 12 blood cells. The distribution and expression of the results were expressed in at least 576 immunohistochemical stains and read and indexed by a professional. The examined tissues were from 144 different individuals and 216 tumor tissues, ensuring that the staining results were fully representative.[52]^17 Data in [53]Figure 1B and C showed IRAK1 expressed in breast cancer cells (MCF-7) and breast cancer tissues, respectively. To further explore the expression pattern of IRAK1 in patients treated with NCT, immunohistochemistry was carried out in 75 patients. IRAK1 was mainly expressed in the cytoplasm. IRAK1 expression level was reduced obviously following NAC treatment in 58/75 patients (77.3%). Representative staining images of IRAK1 in matched breast tumor tissues ([54]Figure 2A and C) pre- and ([55]Figure 2B and D) post-NAC chemotherapy were shown in [56]Figure 2, Wilcoxon signed rank test revealed a significant difference in IRAK1 expression between pre-NAC and post-NAC (IRAK1, P<0.001), [57]Figure 3A and B. Spearman’s correlation analysis was performed between IRAK1 expression (both pre- and post-NAK) and clinicopathological characteristics (including the age when the patient was diagnosed, molecular subtypes, histological grade, TNM stage, ER, PR, HER-2 levels, and lymph node status), and it was shown that the expression of IRAK1 both pre-NAC and post-NAC was not related to the patient’s clinicopathological characteristics. Figure 1. [58]Figure 1 [59]Open in a new tab IRAK1 expression in The Cancer Genome Atlas (TCGA) database and The Human Protein Altas. Notes: (A) IRAK1 was highly expressed in breast cancer tissues compared with normal breast tissues based on TCGA database. (B) IRAK1 expressed in MCF-7 cell line. (C) IRAK1 expressed in breast cancer tissue. Figure 2. [60]Figure 2 [61]Open in a new tab Immunohistochemistry images of IRAK1 expression in matched breast tumor tissues pre- and post-neoadjuvant chemotherapy (NCT). Notes: (A–D) Representative staining images of IRAK1 in matched breast tumor tissues pre- and post-NCT. (A) Strongly stained. (B) Moderately stained. (C) Moderately stained. (D) Negatively stained. Magnification ×400. Figure 3. [62]Figure 3 [63]Open in a new tab The expression of IRAK1 before and after neoadjuvant chemotherapy. Notes: (A) Difference in expression before and after neoadjuvant chemotherapy. (B) Distribution of the values before and after chemotherapy, respectively. (C) Patients whose tumor sizes diminished by <2 cm, showed no change or increased, were defined as group A, while patients whose tumor sizes diminished by >2 cm or those who achieved pathological complete response were defined as group B. The Mann–Whitney U test was used (P=0.019). Decrease in IRAK1 expression following NAC was associated with reduction of tumor size The relationship between the reduction of tumor size and the expression of IRAK1 pre- and post-NAC was investigated. Patients were divided into two groups according to the reduction of tumor size. For patients with tumors that had diminished by <2 cm, and the patients with an unchanged or increased tumor size were defined as group A, while patients whose tumor sizes had diminished by >2 cm, or patients who had achieved pCR were defined as another group (group B). Our data showed that there was a significant difference in the decrease of IRAK1 expression between the two groups (Mann–Whitney U test, P=0.019, [64]Figure 3C). The correlation between DFS and the expression of IRAK1 pre- and post-NAC treatment First, we analyzed the correlation between IRAK1 expression and survival online ([65]www.kmplot.com). In order to assess the differential survival with respect to the expression of IRAK1 pre- and post-NAC, as well as the expression change of IRAK1 following NAC, Kaplan–Meier survival analyses were performed. ROC curve analyses were used to dichotomize the expression scores into high and low expression groups. The cutoff values were obtained from the highest combined sensitivity and specificity at the end point of DFS, and were shown as follows: pre-NAC:7; post-NAC:5. A high expression of IRAK1 both in pre-NAC and post-NAC was found to be significantly correlated with poor DFS. High level of IRAK1 at post-NAC was more valuable in predicting poor DFS compared with that at pre-NAC ([66]Figure 4). Figure 4. [67]Figure 4 [68]Open in a new tab ROC curve and survival curve of IRAK1 expression both pre-neoadjuvant chemotherapy (NCT) and post-NCT. Notes: (A) ROC curve of IRAK1 pre-NCT, an AUC value of 0.681 indicated (P=0.008). (B) Kaplan–Meier survival curve of IRAK1 pre-NCT (log rank: P=0.001). (C) ROC curve of IRAK1 post-NCT. An AUC value of 0.794 (P<0.001). (D) Kaplan–Meier survival curve of IRAK1 post-NCT (log rank: P<0.001). IRAK1 alteration analyzed in cBioPortal database The cBio Cancer Genomics Portal ([69]http://cbioportal.org) is an open platform, providing data of more than 5,000 tumor samples from 20 cancer studies, which can help cancer researchers to study in depth.[70]^18^–[71]^20 We investigated the genetic alteration of IRAK1 with a cohort of 816 cases of invasive breast cancer available in the cBioPortal database. The result showed that 6.9% (56/816) of cases of breast cancer in this cohort exhibited IRAK1 alteration ([72]Figure 5A). Moreover, a network of the interaction between IRAK1 and several frequently altered neighboring genes was established ([73]Figure 5B). Survival analysis (OS and DFS) was performed to compare the prognosis of the patient groups with and without IRAK1 alteration. However, the results did not show statistical significance ([74]Figure 5C and D). Figure 5. [75]Figure 5 [76]Open in a new tab The alteration of IRAK1 and its interaction in altered neighboring genes from the cBioPortal database. Notes: (A) IRAK1 gene alteration. In this cohort, 7% (56/816) of invasive breast cancer cases exhibited IRAK1 gene alteration. (B) Network of IRAK1 and altered neighboring genes. Analysis of the overall survival (C) and disease-free survival (D) of invasive breast cancer patients with and without IRAK1 gene alteration. No statistically significant difference was identified. IRAK1 protein expression in the String database The String database is a database that is used to search the potential interactions between known proteins and predicted proteins. This database can be applied to 2,031 species and contains interactions between 9.6 million proteins and 13.8 million proteins.[77]^21^,[78]^22 An analysis was performed to investigate the correlations between IRAK1 expression and its potential targets. The results revealed that the expressions of several target proteins were significantly correlated with IRAK1 expression according to the GEPIA online database.[79]^23 (MYD88: r=0.19, IKBKG: r=−0.15, IL1R1: r=−0.073 TAB1: r=0.048, TRAF6: r=−0.035, IRAK4: r=−0.055, NFKBIA: r=0.01, IRF7: r=0.12, TAB2: r=−0.021, TOLLIP: r=0.092, all P<0.05 besides TRAF6). Remarkably, MYD88 and IRF7 showed a relatively good correlation with IRAK1 ([80]Figure 6). Both MYD88 and IRF7 play a major role in the TLR4 signaling pathway, which is consistent with Wee et al’s study of IRAK1 in breast cancer. This also suggests that IRAK1 may mediate its role in NCT for breast cancer via TLR signaling pathway. Figure 6. [81]Figure 6 [82]Open in a new tab IRAK1 protein expression in the String database. Notes: (A) Corresponding proteins of IRAK1 protein in network visualization. (B,C) Pearson correlation analysis revealed that IRAK1 expression was significantly correlated with the expression of IRF7 and MYD88, respectively. Gene set enrichment analysis of IRAK1 in WebGestalt 2017 WebGestalt (Web-based Gene SeT AnaLysis Toolkit, [83]http://www.webgestalt.org) was one of the first web tools for functional enrichment analysis.[84]^24^,[85]^25 The GOSlim classification plot is used to show the distribution of IRAK1 across the major branches of the GO biological process, cellular component and molecular function ontologies. The enrichment results were shown in [86]Figure 7. The results showed that IRAK1 is mainly localized in the cell membrane in breast cancer, mainly enriched in biological regulation, metabolic process, and protein binding, which also paves the way for further research on the mechanism of IRAK1 in NCT for breast cancer in the future. Figure 7. [87]Figure 7 [88]Open in a new tab The enrichment results of IRAK1 using web-based Gene SeT AnaLysis Toolkit. Notes: (A) IRAK1 enrichment in biological process, cellular component, and molecular function. (B) Positively related genes of IRAK1. (C) Negatively related genes of IRAK1. (D) Volcano plot of IRAK1 in breast cancer. Discussion IRAK1 has been shown to play an important role in the IL-1/TLR signaling pathway. TLRs are key regulators of the innate immune response and have been observed in many malignancies,[89]^26 including breast cancer. IRAK1, a downstream effector of TLR2, inhibits IRAK1 expression and can also inhibit the growth of breast cancer in vitro and in vivo, indicating that IRAK1 may be a potential therapeutic and prophylactic target for breast cancer through the TLR2 pathway.[90]^27 TLR4 overexpression in human breast tumors is often associated with paclitaxel chemoresistance.[91]^28 In recent years, some studies have confirmed that IRAK1 can act as a proto-oncogene to promote tumor progression, metastasis, and chemotherapy resistance. IRAK1 is overexpressed and over-activated in myelodysplastic syndrome (MDS). Treatment with small-molecule IRAK1 inhibitors (IRAK1/4-Inh) resulted in increased proliferation and apoptosis of clones of MDS, which was consistent with the inhibition of TRAF6/ NF-κB. These findings suggest that IRAK1 can be used as a therapeutic target for MDS.[92]^29 The expression level of IRAK1 protein may be used as a parameter to identify whether patients are at high risk for lung cancer and may serve as a potential target for lung cancer diagnosis and gene therapy.[93]^13 IRAK1 is an oncogene in the development of hepatocellular carcinoma and is overexpressed in hepatocellular carcinoma. In addition, by analyzing the immunohistochemistry and RNA sequencing data in cancer genomic maps, IRAK1 is expected to become a diagnostic and therapeutic target for hepatocellular carcinoma.[94]^30 However, to the best of our knowledge, there has been no study investigating the potential role of IRAK1 in NCT for breast cancer. In this study, we, for the first time, used immunohistochemistry to detect the expression of IRAK1 protein in breast cancer tissues before and after NCT in breast cancer patients undergoing NCT. Neoadjuvant systemic chemotherapy is a possible method for the treatment of locally advanced, operable, primary non-surgical or inflammatory breast cancer. There is a clear correlation between response to chemotherapy drugs and DFS and OS in breast cancer patients. However, the great challenge for NAC is chemotherapy resistance, and there is also a lack of excellent predictors of NAC response. Although many different prognostic indicators are being developed, few molecular markers are used to predict NAC responses. In this study, we found that changes in IRAK1 expression were correlated with changes in tumor size before and after chemotherapy. Although this finding was not sufficient to predict chemotherapy sensitivity prior to NAC, it can be used as a marker for evaluating the efficacy of NCT and as a reference for postoperative chemotherapy regimens. In this study, we first studied the correlation between IRAK1 expression and the clinicopathological parameters of breast cancer and its impact on prognosis based on the clinical data and mRNA microarray data of 1,085 breast in NCT patients was analyzed online with km-plotter. Immunohistochemical method was used to detect the expression pattern of IRAK1 in NCT breast cancer patients. The results showed that the expression of IRAK1 after NCT was significantly lower than before NCT, and the decreased expression of IRAK1 was correlated with change of tumor size. IRAK1 may be a potential prognostic indicator of the efficacy of NCT. The Human Protein Atlas validated IRAK1 protein expression in breast cancer tissue. The IRAK1 mutations were analyzed online using cBioPortal database to understand IRAK1 mutations in breast cancer and their association with survival. The IRAK1 protein interaction network was obtained through the String database, and its upstream and downstream related proteins were predicted and verified in TCGA database. Gene enrichment analysis of IRAK1 was performed using the web-based Gene SeT AnaLysis Toolkit to understand its biological process. These bioinformatics analyses partially revealed the mechanism of IRAK1 in breast cancer. However, the molecular mechanism of IRAK1 in NCT for breast cancer remains unclear. Bioinformatics analysis showed that the mutation rate of IRAK1 gene in breast cancer was as high as 7%, and most of the mutant forms were amplified. Significantly altered neighboring genes within the network are APK2, MAP2K6, MAP2K4, MAP3K3, AKT3, MTOR, PRKCA, and TLR5. Therefore, IRAK1 may also promote the development of breast cancer through the MAP3K–MAP2K–MAPK signaling pathway and the AKT– mTOR signaling pathway in addition to the TLR signaling pathway. Enrichment analysis suggested that IRAK1 plays a major role in biological regulation, metabolic processes, and protein binding. These findings will guide our further study in exploring the mechanism of action of IRAK1 in NCT. Footnotes Disclosure We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work. There is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of the manuscript. The authors report no other conflicts of interest in this work. References