Abstract
Background:
Currently, one of the most used strategies for the treatment of newly diagnosed patients with breast cancer is neoadjuvant chemotherapy based on the application of taxanes and anthracyclines. However, despite the high number of patients who develop a complete pathological clinical response, resistance and relapse following this therapy continue to be a clinical challenge. As a component of the innate immune system, the cytotoxic function of Natural Killer (NK) cells plays an important role in the elimination of tumor cells. However, the role of NK cells in resistance to systemic therapy in breast cancer remains unclear. The present project aims to evaluate the gene expression profile of human NK cells in breast cancer tissue resistant to treatment with taxanes–anthracyclines.
Methods:
Biopsies from tumor tissues were obtained from patients with breast cancer without prior treatment. Histopathological analysis and
Results:
A significant decrease in gene expression of cell-surface receptors related to NK cells was observed in tumor samples resistant to antineoplastic treatment compared with those that were sensitive to treatment.
Conclusion:
A decrease in NK cell infiltration into tumor tissue might be a predictive marker for failure of chemotherapeutic treatment in breast cancer.
Keywords
Background
Locally advanced breast cancer is a heterogeneous disease characterized by increased local recurrence and decreased survival. It is classified into stages IIB to IIIC by the American Joint Committee on Cancer, 1 -3 and the prognosis of patients with this disease remains highly variable. Multimodal therapy, including chemotherapy, surgery, and radiation, is the treatment of choice that has shown the highest percentage survival in such patients with breast cancer. 4,5 Particularly, chemotherapeutic schemes based on anthracycline and taxanes are widely accepted as first alternatives in the treatment of breast cancer 6 -8 ; these therapeutic strategies provide a benefit in recurrence-free survival as well as overall survival. 9,10
Management in clinical practice, using a neoadjuvant chemotherapy regimen in either a concomitant or sequential manner, is preferred in locally advanced breast cancer. The main objectives of the neoadjuvant chemotherapy regimen are to eliminate distant micrometastases that exist at the time of diagnosis, decrease tumor size allowing conservative surgeries, and to assess the
Predictive factors associated with greater probability of achieving a pCR using neoadjuvant chemotherapy regimen are tumor size, histological type (ductal–lobular carcinoma), intrinsic subtype tumor (basal–luminal), hormone receptor status (negative–positive estrogen receptor), expression of human epidermal growth factor receptor 2 (HER2), Ki-67, and the Scarff-Bloom-Richardson grade. 11,12 In this sense, large numbers of studies have focused on specific characteristics of phenotype, molecular patterns, and growth rates of tumor cells after chemotherapy. For instance, there is now sufficient evidence that chemotherapy regimens employing anthracyclines and taxanes lead to higher pCR rates in triple-negative tumors compared to estrogen and progesterone receptor-positive subtypes of breast cancer. Regrettably, only a fraction of the patients categorized in a specific subtype of breast cancer respond to neoadjuvant chemotherapy and achieve a pCR showing a better prognosis. In fact, despite the benefits of using modern chemotherapy regimens, multidrug resistance continues to be a clinical challenge, and the need for biological markers that can predict the response to chemotherapy is evident.
Distinguishing responsive from nonresponder patients can significantly improve therapeutic decisions. Therefore, much effort has been focused on the identification of specific features of the tumor microenvironment including biological and molecular markers that could help to tailor or develop new therapies. Recently, the presence of tumor infiltrating lymphocytes (TILs) has been correlated with clinical outcomes in many types of cancer. 13,14 Patients with breast cancer with prominent lymphocyte infiltration, before any treatment, show stronger responses to neoadjuvant chemotherapy. 15,16 Moreover, high lymphocyte infiltration correlates with higher pCR rates as well as with better patient prognosis. Because of this, it has been suggested that infiltration of tumor-associated lymphocytes may represent a new independent predictive factor of response to neoadjuvant chemotherapy in breast cancer. 17,18
The composition of tumor infiltrating immune cells can vary according to cancer type. Moreover, different types of infiltrating immune cells involved in both innate and adaptive immunity have diverse effects on tumor behavior with either pro-tumoral or antitumoral consequences. 14,16 In breast cancer, it has been shown that diverse patterns of gene expression in the tumoral stroma are related to good prognosis. Interestingly, overexpression of a group of genes related to the immune response, including T cells and NK (Natural Killer) cell markers, indicating a TH1-like response (granzyme A, CD52, CD247 and CD8A), are indicators of good prognosis. 19 In addition, it has been shown that tumor infiltration by T and B cells is associated with high response rates to neoadjuvant chemotherapy as well as a high rate of pCR. 18,20,21
To date, there is no clear evidence to establish an association between NK cell infiltration and clinical outcome in patients with breast cancer. Therefore, the present study aims to evaluate the expression of NK cell surface receptors in breast cancer tissues and their association with the pathological response to neoadjuvant chemotherapy.
Methods
Obtention of Tumor Explants
Infiltrating carcinoma specimens were collected from patients with breast cancer during surgery at the Hospital of Gynecology and Obstetrics of the Mexican Institute of Social Security (IMSS). Immediately after surgery, to confirm its tumorous nature and to avoid contamination, the specimen was dissected and evaluated by a pathologist. One small piece of tissue was collected in cold serum-free Dulbecco's Modified Eagle Medium: Nutrient Mixture F-12 (DMEM/F12) medium (Invitrogen, Grand Island, New York) and transported at 4°C to the organotypic culture laboratory for its immediate processing; another piece of tissue was stored in RNAlater RNA Stabilization Reagent (Qiagen, Hilden, Germany) at −80°C.
Evaluation of TILs in formalin-fixed, paraffin-embedded Specimens
Clinical and histopathological information were obtained from patients’ files archived at the Gynecology Hospital and Specialty Hospital of the Western National Medical Center (CMNO) of the IMSS, and the percentage of TILs was evaluated from formalin-fixed, paraffin-embedded (FFPE) specimens of invasive breast ductal carcinoma (n = 90) that were retrospectively obtained.
Tumor Infiltrating Lymphocytes
The distribution pattern of tumor cells and the presence of inflammation that accompanies the tumor was determined by observation under light microscopy using the 5× objective. Following this, at 10× magnification, the proportion of tumor to inflammatory cells was assessed according to Supplemental Figure 1. Finally, after analysis of different representative areas of the tissue, a mean percentage was calculated.
Treatment of Tumor Explants With Antineoplastic Drugs
From representative tumor samples, small tumor explants (4-5 mm in diameter and 250-300 μm in thickness) were prepared using the Krumdieck slicing system (Alabama Research and Development Corporation, Birmingham, Alabama) as described previously. 22 Tumor explants were placed in 6-well microplates containing DMEM/F12 culture medium–supplemented with 10% fetal bovine serum, 5 μg/mL bovine insulin, 100 μg/mL gentamicin, insulin–transferrin–selenium, and 25 mM glucose (DMEM/F12-supplemented medium). Plates were preincubated for 1 hour at 37°C, 5% CO2, 95% air, and agitated at 35 rpm. Time from harvest to slicing was kept to an absolute minimum (<2 hours). The entire process was performed under aseptic conditions. The treatment of tumor explants with antineoplastic compounds was performed after 1 h of preincubation. The tumor explants were transferred to 24-well microplates and treated with 20 μg/mL paclitaxol (TX). The control group (100% viability) consisted of untreated explants, which were incubated only with culture medium.
Alamar Blue Viability Assay
The effect of treatments on the viability of the tumor explants was assessed by the alamar Blue (AB) assay. After 48 hours of incubation with compounds, as well as with cell culture medium (control), the explants were incubated for an additional 4 hours with 10% AB in 500 μL DMEM/F12-supplemented medium at 37°C in the conditions described above. Afterward, 100 μL was collected from each sample and transferred to a 96-well microplate. Fluorescence values were read using a multimode microplate reader (Synergy BioTek HT, Winooski, Vermont) at 530 nm excitation/590 nm emission wavelengths.
Histological Analysis
Explant tissue sections prior to treatment were fixed in 10% neutral formalin and then embedded in paraffin using the conventional histological technique. Formalin-fixed, paraffin-embedded specimens of invasive breast ductal carcinoma were obtained from the Pathology Department of the CMNO, IMSS. From all samples, tissue sections of 4 μm were prepared on a microtome and mounted on glass slides, sections were deparaffinized and hematoxylin-eosin (H&E) stained. The preparations were then evaluated by a pathologist.
RNA Isolation and Gene Expression Microarrays
Biopsies stored at −80°C in RNAlater were used for the isolation of total RNA with the RNeasy Plus Mini kit (Qiagen) according to the manufacturer’s instructions. RNA was quantified by absorbance at 260/280 nm, and quality was evaluated using an Agilent 2100 Bioanalyzer and the RNA 6000 NanoChip kit (Agilent Technologies, Santa Clara, California). Double-stranded complementary DNA (ds-cDNA) was generated from 5 µg of RNA using the cDNA Synthesis Kit System (Roche Applied Science, Indianapolis, Indiana) and purified with the GenElute polymerase chain reaction (PCR) Clean-up Kit (Sigma-Aldrich Quimica, Toluca, Mexico). Afterward, the ds-cDNA was labeled with Cy3 and hybridized in a Human Gene Expression Array 12 × 135 k (Roche Applied Science). Finally, the microarrays were scanned on an MS200 Scanner (Roche Applied Science), and the data obtained were processed using the DEVA 1.2 software (Roche Applied Science). Fluorescence intensities were normalized using the RMA algorithm, and all data obtained were subsequently analyzed using the CLC Main WorkBench version 7.0.3 software (Qiagen). Intensities of sensitive versus resistant samples were compared and the genes were considered differentially expressed when the intensities showed a difference of >1.5 and a
In Silico Analysis
To further explore the molecular pathways in which differentially expressed genes are involved, we performed an analysis using the WEB-based-Gene-Set-Analysis toolkit (WebGestalt, http://www.webgestalt.org/). Additionally, genes associated with NK cell-mediated cytotoxicity were further introduced into the cBio Cancer Genomics Portal (cBioportal, http://www.cbioportal.org/); the breast invasive carcinoma (TCGA, provisional) was selected and the messenger RNA (mRNA) expression data with a
Real-Time Polymerase Chain Reaction and PCR Evaluation
Complementary DNA synthesis was performed from 5 μg of total RNA using the Transcriptor First Strand cDNA Synthesis Kit (Roche Applied Science) primed with Oligo dT. The cDNA obtained was employed in quantitative polymerase chain reactions to evaluate gene expression levels using the LightCycler FastStart DNA Master PLUS SYBR Green I Kit with 2.0 LightCycler technology (Roche Applied Science) under conditions specified by the manufacturer. Primer pairs were designed using Oligo Primer Analysis version 7.0 software (Molecular Biology Insights, Inc, Colorado Springs, Colorado) from sequences obtained at the Entrez Nucleotide Database of the National Center of Biotechnology Information (Table 1). Polymerase chain reaction products were resolved by agarose (2%) gel electrophoresis.
Primer Pairs Designed for Cell Surface Receptors Mainly Related to Natural Killer Cells.
Quantification of PCR bands
Densitometric analysis of PCR bands was carried out by using the ImageJ software version 1.51j8 (https://imagej.nih.gov/ij/). Densitometric values of RPL32 in each sample were used as reference.
Statistical Analysis
Statistical analysis was performed with SPSS version 22.0 software (IBM Corporation, Armonk, New York). Quantitative data were expressed as mean and standard error of the mean. Differences between treatments of tumor slices were analyzed by Student
Results
Response of Tumor Slices to Chemotherapeutic Treatment
Once confirmed that breast cancer slices remained viable and actively proliferating during the

Percentage viability of tumor slice samples after incubation with paclitaxol (TX). Samples with viability <25% were designated as sensitive. Data are expressed as mean ± standard of the mean (SEM) *
Clinical and Histopathological Data of Patients With Breast Cancer From Whom Tumor Slice Samples Were Obtained.
Abbreviations: ER, estrogen receptor; HER, human epidermal growth factor receptor; M, metastasis; N, node; PR, progesterone receptor; T, tumor.
Tumor Infiltrating Lymphocytes Are Decreased in Breast Cancer Tissue Resistant to Chemotherapy
Once tumor slice samples were grouped into sensitive and resistant, histopathological analysis was performed not only on breast slice cultures incubated with TX but also in tumor slices without any treatment. Hematoxylin-eosin staining showed a differential proportion of TILs between tumor slice samples cataloged as sensitive (9A, 11B, 14C, and 19D) compared to slice samples cataloged as resistant (4E, 6F, 23G, and 30H). High-level lymphocyte infiltration was evident in tumor slices sensitive to TX treatment (Figure 2).

Photomicrographs showing tumor samples from mammary gland samples stained with hematoxylin and eosin (H&E) to qualitatively evaluate the presence of tumor infiltrating lymphocytes. Sensitive tumor slices (A, B, C, and D) and resistant tumor slices (E, F, G, and H) without prior exposure to therapeutic drugs.
Additionally, TILs were evaluated in FFPE specimens of invasive breast ductal carcinoma (n = 90), obtained retrospectively within the last 5 to 10 years.
Evolution and follow-up information of patients with breast cancer was obtained from clinical records, and the specimens were classified as either good or bad prognosis, according to clinical parameters such as histological type, clinical stage, hormonal receptors, and molecular classification, but above all response to chemotherapy (pCR) and the presence of recurrent disease.
Initially, an analysis of the percentage of TILs was performed in all patients with breast cancer (n = 90) divided into patients with better prognosis (n = 42) and bad prognosis (n = 48). An increased tendency of tumor infiltrate was found in patients with breast cancer with better prognosis (28.9% ± 3.2%) compared to patients with breast cancer with bad prognosis (21.4% ± 2.4%; Figure 3A).

Percentage of tumor infiltrating lymphocytes in tumor samples of patients with breast cancer. (A) Percentage of tumor infiltrate according to patient prognosis and (B) in relation to the molecular subtype classification.
Since molecular classification in breast tumors has been employed to inform decisions regarding therapeutic treatment, the percentage of TILs was analyzed dividing patients with breast cancer into luminal A (n = 14), luminal B (n = 26), HER2 (n = 13), and triple negative (n = 37) tumors. Overall, an increased tendency for presence of TILs was also seen in patients with better prognosis compared to patients with bad prognosis classified into luminal A (15.0% ± 5.4% vs 12.5% ± 3.5%), luminal B (30% ± 5.2% vs 22.7% ± 4.4%), and triple negative tumors (34.5% ± 6.3% vs 22.56% ± 3.5%; Figure 3B).
Decreased Gene Expression of NK Cell Surface Receptors in Breast Tumors Resistant to Chemotherapy
Once a tumor slice was identified as sensitive or resistant, and the histopathologic analysis by H&E showed a decrease in TILs in either tumor slice samples or FFPE specimens of invasive breast ductal carcinoma, expression microarrays were performed to characterize the molecular profile of sensitive tumor slices compared to resistant tumor slices. Sensitive tumor slices were taken as reference, and variations in the gene expression levels were evaluated with an absolute difference value >1.5 and
Genes Regulated Significantly in Sensitive and Resistant Tumor Slices Associated With KEGG Main Routes.a
Abbreviation: FCM, fragment crystallizable M.
aSensitive tumor slices were taken as reference.
NK cell–mediated cytotoxicity was found to be one of the most significant routes, associated with down-modulation of 30 genes in resistant tumor slices compared to sensitive tumor slices, with an absolute difference value of >1.5 (
Gene Expression (Displayed as Difference) of Cell Surface Receptors Mainly Related to Natural Killer Cells. Samples of Sensitive Tumor Slices Were Taken as Reference.a
aAll genes down-regulated with an absolute difference value >1.5 (
To validate the observations found in the microarray experiments, all differentially expressed genes related to NK cells were evaluated by PCR in both sensitive and resistant tumor slice samples. For these assays, RPL32 was used as reference gene and a peripheral blood mononuclear cell sample was included as a positive control. Down-modulated expression in resistant tumor slice samples of all analyzed genes described in Table 3 was confirmed by PCR. Final PCR products were resolved by agarose (2%) gel electrophoresis to confirm the size and specificity of the product. Image processing and quantitative analysis was performed by ImageJ software; representative results are shown in Figure 4. KLRC1 to 4 were the most significantly down-modulated genes in resistant tumor slices compared to sensitive counterparts. However, all selected genes described in Table 3 showed a downward trend in resistant tumor slices compared to sensitive samples.

Agarose DNA electrophoresis of polymerase chain reaction (PCR) amplification products from complementary DNA (Cdna) of tumor slices resistant or sensitive to antineoplastic treatment (upper panels), arrows are indicating the expected products. Bands outer the expected size were considered PCR artifacts, and the 1 kb plus DNA ladder was used. Periferal Blood Mononuclear Cells (PBMCs) amplification products were included as positive controls and RPL32 as reference gene. Densitometric comparisons are shown in the lower panels. *Indicates significant differences between groups (
The genes of interest related to NK cells that were found to be significantly modulated in the present study (expression microarrays) were introduced in cBioportal and a breast invasive carcinoma cancer study was selected, which included 1105 samples, and the mRNA expression performed by RNAseq was analyzed. Finally, the differential gene expression patterns of query genes were compared with clinical tracks (overall survival and disease-free status). Messenger RNA upregulation of effector molecules (GZMB and PRF1), natural cytotoxicity receptors (NCR1), KIRs, and KLRCs were associated with overall survival (living) and disease-free status in patients with breast cancer (Figure 5).

Expression pattern of selected genes of the natural killer cell–mediated cytotoxicity pathway in patients with invasive breast cancer classified according to overall survival status and disease-free status.
Discussion
Identification of new and reliable factors that predict responses to neoadjuvant chemotherapy are urgently needed in clinical practice. Treatment decisions in breast cancer, particularly involving chemotherapy, are principally based on tumor cell phenotype; however, the percentage of patients who achieve pCR is still relatively low. 23 It has now been shown that the host immune system not only helps to elucidate cancer prognosis, it can also determine the response to antineoplastic therapies. 24 In breast cancer, the clinical relevance of TILs has been clearly established by their correlation with pCR as well as with long-term positive outcomes. 15,25 Currently, only limited data regarding particular TILs are available, and their clinical assessment has not yet been fully standardized.
Now, however, measuring the overall extent of TILs is not enough because it does not distinguish between the diversity of innate and adaptive immune cells, thereby limiting the knowledge of their specific contributions in the tumor tissue. 26 In this way, the strong association between TILs and the response to neoadjuvant chemotherapy suggests that characterization of the nature and location of the immune infiltrate should be considered.
The composition of the immune infiltrate is influenced by the type of cancer; 15 in breast cancer, there is evidence that neoadjuvant chemotherapy is particularly efficient if the patient shows signs of an antitumor immune response both preexisting and induced by therapy. 25,27 Regarding neoadjuvant anthracycline-based chemotherapy, a study demonstrated that only tumor size and TILs were independent predictors of anthracycline response in breast cancer. 28
In this work, we observed that tumor slice samples sensitive to TX treatment showed high levels of lymphocyte infiltration compared to resistant tumor slices. Furthermore, an analysis of FFPE specimens (biopsies) of invasive breast ductal carcinoma from patients who received taxane and anthracycline-based chemotherapy treatment demonstrated an increased tendency for the presence of TILs in samples classified as good prognosis compared to bad prognosis counterparts. The same behavior was found when dividing breast cancer specimens into molecular subtypes (luminal A, luminal B, and triple negative), demonstrating that TILs may be useful as an independent predictive factor.
Different chemotherapeutic agents have been described to stimulate antitumor immune responses 29 ; among them, cyclophosphamide, anthracyclines, and taxanes have been demonstrated to stimulate T cell proliferation, activate NK cells, deplete circulating regulatory T cells (Treg), and inhibit Treg infiltration into tumor tissues. 30 Furthermore, immunogenic cell death induced by chemotherapy releases tumor-associated antigens that can be recognized by antigen-presenting cells, leading to lysis of tumor cells by activated CD8+ cytotoxic T cells. 25,31,32 Because interactions between TILs and tumor cells are relevant for clinical outcomes, it is extremely important to understand the preexisting inflammatory infiltrate in the tumor before subjecting it to chemotherapeutic treatment. Evidence is accumulating that high levels of TILs in breast cancer compromising T (CD3+ CD8+ FOXP3+) and B cell populations are significantly associated with pCR and increased disease-free survival, especially if the tumors also lack immunosuppressive CD68+ macrophages. 18,25,28,33 -35
Due to the variability regarding the characterization of TILs, which may partly reside in the sensitivity of the detection methods used, the controversy continues as to what will be the best method to evaluate specific populations of TILs. This method must be easy to use and to standardize for application in the clinic. Most studies have used immunohistochemical methods to investigate the composition of immune populations infiltrating breast tumors, but a few have identified immune transcripts and gene sets relevant to predicting pCR in breast cancer subtypes. 36 -38 Chemokine (C-X-C motif) ligand 13 is a metagene signature reflecting not only T-cell infiltration but also intratumoral presence of interferon-γ-producing T cells and has been demonstrated as metagene predicting pCR after neoadjuvant breast cancer chemotherapy. 25
Our gene expression analysis, used to compare molecular profiles of sensitive tumor slices with resistant tumor slices, showed that highly significant modulated genes were related to NK cells. Thirty genes were found to be down-modulated in resistant tumor slices, among them Nkp46, KIRs (activator and inhibitor receptors), KLRCs (activator and inhibitor receptors), and effector molecules (granzyme B and perforin 1). Results obtained by gene expression microarray were validated by PCR, and so these results might indicate more than just a regulation in their state of activation, but a deficiency in the infiltration of NK cells in tumor tissue resistant to chemotherapy prior to any treatment.
NK cells are a vital component of the innate immune system that plays an important role in tumor immune surveillance and in the prevention of progressive tumor growth as well as in defense against metastatic progress. 39,40 It has been demonstrated that most human tumors show low levels of NK cell infiltration; however, human tumors that present more substantial infiltration of NK cells have been associated with improved prognosis and reduction in tumor recurrence. 41 -43
A previous report by Verma et al 44 showed a significant inhibition of a specific phenotype of circulating NK cells and NK cell cytotoxicity in women with breast cancer. A significant reduction in cytotoxicity mediated by circulating NK cells in tumors that responded poorly to neoadjuvant chemotherapy was found. On the other hand, their study showed a significant increase in NK cells in the peritumoral environment in patients who achieved a good pathological response. However, intratumoral NK cells did not show the same correlation. 44 In our study, RNA obtained from the whole tumor sample showed significantly down-modulated expression of NK cell receptor genes (KLRC1 to 4) in resistant tumor slices compared to sensitive counterparts.
Finally, it has been demonstrated that molecular signatures associated with NK cells are predictive markers of relapse-free survival in patients with breast cancer.
45
These results are in accordance with our
However, aspects related to the migration and presence of NK cells in tumor tissue, and their correlation with achievement of a pCR following chemotherapeutic treatment in breast cancer, remain poorly documented in the literature.
Conclusion
A signature of NK cell-related genes in breast cancer might be useful for identification of a group of tumors highly sensitive to neoadjuvant chemotherapy. A diminished inflammatory infiltrate, as well as decreased expression of NK cells receptor genes, may be predictive markers for a poor chemotherapeutic response.
Supplemental Material
Supplementary Figure - Expression of NK Cell Surface Receptors in Breast Cancer Tissue as Predictors of Resistance to Antineoplastic Treatment
Supplementary Figure for Expression of NK Cell Surface Receptors in Breast Cancer Tissue as Predictors of Resistance to Antineoplastic Treatment by Garcia-Chagollan Mariel, Carranza-Torres Irma Edith, Carranza-Rosales Pilar, Guzmán-Delgado Nancy Elena, Ramírez-Montoya Humberto, Martínez-Silva María Guadalupe, Mariscal-Ramirez Ignacio, Barrón-Gallardo Carlos Alfredo, Pereira-Suárez Ana Laura, Aguilar-Lemarroy Adriana, and Jave-Suárez Luis Felipe in Technology in Cancer Research & Treatment
Footnotes
Abbreviations
Authors’ Note
The data that support the findings of this study are available at the Gen Expression Omnibus (GEO,
) Series record: GSE99225.
This protocol was approved by the Ethical and Research National Committee of IMSS, with registration numbers: R-2014-785-022 (patient recruitment to test ex vivo sensitivity to taxanes) and R-2013-785-061 (obtaining formalin-fixed, paraffin-embedded [FFPE] biopsies and their clinical information to assess the infiltration of immune cells). Informed consent was obtained from all patients involved in this study. Clinical and personal information from retrospective samples (FFPE biopsies) collected from the Department of Pathology of the CMNO of IMSS were handled confidentially. All procedures performed in this study were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Authors’ Contributions
G.-C.M. performed microarray, PCR, and qRT-PCR experiments; she was also involved in the recruitment of patient samples and drafting the manuscript. C.-T.I. E. and C.-R.P. recruited biopsies of patients with breast cancer and evaluated its
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was financially supported by Fondo de Investigación en Salud, IMSS (FIS/IMSS/PROT/PRIO/14/030, to J-SLF).
Supplemental Material
Supplementary material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
