Impact Factor 4.848 | CiteScore 3.5
More on impact ›

Original Research ARTICLE

Front. Oncol., 05 August 2020 | https://doi.org/10.3389/fonc.2020.01206

Biological Landscape of Triple Negative Breast Cancers Expressing CTLA-4

  • 1Departamento de Biología Molecular y Biotecnología, Instituto de Investigaciones Biomédicas, UNAM, Ciudad de Mexico, Mexico
  • 2Departamento de Medicina Genómica y Toxicología Ambiental, Instituto de Investigaciones Biomédicas, UNAM, Ciudad de Mexico, Mexico
  • 3Programa Institucional de Cáncer de Mama, Instituto de Investigaciones Biomédicas, UNAM, Ciudad de Mexico, Mexico
  • 4American British Cowdray Medical Center, Cancer Center, Ciudad de Mexico, Mexico
  • 5American British Cowdray Medical Center, Pathology Service, Ciudad de Mexico, Mexico
  • 6Biochemistry Department, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Ciudad de Mexico, Mexico
  • 7Red de Apoyo a la Investigación (RAI), Universidad Nacional Autónoma de Mexico-Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico City, Mexico
  • 8Laboratorio de Células Troncales Mesenquimales, Unidad de Investigación Médica en Enfermedades Oncológicas, Hospital de Oncología, Centro Médico Nacional Siglo XXI, Instituto Mexicano del Seguro Social (IMSS), Ciudad de Mexico, Mexico
  • 9Escuela de Medicina, Universidad Panamericana, Ciudad de Mexico, Mexico

Patients with triple-negative breast cancer (TNBC) have a poor prognosis, partly because of the absence of targeted therapies. Recognition of the key role of immune responses against cancer has allowed the advent of immunotherapy, focused on the inhibition of negative immune checkpoints, such as CTLA-4. CTLA-4 is also expressed in some cancer cells, but its activity in tumor cells is not completely understood. Thus, the aim of the present work was to determine the biological landscape and functions of CTLA-4 expressed in TNBC cells through preclinical and in silico analysis. Exploration of CTLA-4 by immunohistochemistry in 50 TNBC tumors revealed membrane and cytoplasmic expression at different intensities. Preclinical experiments, using TNBC cell lines, showed that stimulation of CTLA-4 with CD80 enhances activation of the ERK1/2 signaling pathway, while CTLA-4 blockade by Ipilimumab induces the activation of AKT and reduces cell proliferation in vitro. We then developed an analytic pipeline to define the effects of CTLA-4 in available public data that allowed us to identify four distinct tumor clusters associated with CTLA-4 activation, which are characterized by enrichment of distinctive pathways associated with cell adhesion, MAPK signaling, TGF-ß, VEGF, TNF-α, drug metabolism, ion and amino acid transport, and KRAS signaling, among others. In addition, blockade of CTLA-4 induced increased secretion of IL-2 by tumor cells, suggesting that the receptor regulates cellular functions that may impact the immune microenvironment. This is relevant because a deep characterization of immune infiltrate, conducted using public data to estimate the abundancies of immune-cell types, showed that CTLA-4-activated-like tumors present a conditional immune state similar to an escape phenotype exploited by cancer cells. Finally, by interrogating transcriptional predictors of immunotherapy response, we defined that CTLA-4 activation correlates with high immune scores related to good clinical predicted responses to anti-CTLA-4 therapy. This work sheds new light on the roles of activated CLTA-4 in the tumor compartment and suggests an important interplay between tumor CLTA-4-activated portraits and immune-infiltrating cell populations.

Introduction

Breast cancer is the most prevalent invasive malignancy, and the most common cause of cancer-related death in women worldwide (1). Triple-negative breast cancer (TNBC) is a subtype that lacks the expression of estrogen receptor (ER), progesterone receptor (PgR), and human epidermal growth factor receptor 2 (HER2). Overall, TNBC represents 10–20% of all breast cancer cases (2), being more prevalent in women of Hispanic (3) and African descent (4). TNBC tumors exhibit an unfavorable clinical evolution, given the lack of specific targeted therapies (5). Currently, the standard treatment for TNBC patients is based on chemotherapy (6), and although the initial clinical response is acceptable, the overall prognosis for these patients is poor, compared with other breast cancer subtypes. Thus, there is a great need for the development of new therapeutic options for this pathology.

Recently, the recognition of the key role of immune responses against cancer cells has allowed the advent of immunotherapy based on immune checkpoint inhibitors. Immune checkpoints are negative regulators of cytotoxic T-lymphocyte (CTL) activity. CTLs are the main effector cells that recognize tumor antigens and mediate immune anti-tumor responses; CTL activity is regulated by the interplay between co-stimulatory pathways and checkpoints. After tumor antigen presentation, CTL-associated CD28 receptor binds to B7 (CD80/CD86) on the surface of antigen-presenting cells (APCs), resulting in CTL activation. Afterwards, the cytotoxic T- lymphocyte-associated antigen 4 (CTLA-4) is translocated to the CTL membrane, where it competes for and binds to B7, with a higher affinity than CD28, inhibiting the previously activated pathways (7). CTLA-4 upregulation in cancer patients is recognized as an important mediator of immune evasion. Thus, efforts to block CTLA-4 resulted in the advent of the monoclonal anti-CTLA-4 antibody, Ipilimumab, which was approved to be used as a checkpoint inhibitor for the treatment of melanoma (8, 9). However, the proportion of patients who respond to Ipilimumab remains modest (~22% objective response rate), as tumors have multiple mechanisms of immune evasion (10).

The use of anti-CTLA-4 immunotherapy in other types of cancer has been hampered by the limited understanding of driving responses and resistance pathways. For instance, breast tumors have been considered immunologically quiescent or “cold” tumors. However, recent multi-omic profiling has provided robust characterization of the immune portrait of breast cancers, describing TNBC as a dynamic and heterogeneous immunogenic subtype (11, 12). Since immunotherapy aims at stimulating and restoring anti-tumor immune responses, it is reasonable to hypothesize that the immune microenvironment surrounding tumors is crucial for its success. Recent reports that integrate transcriptional-based deconvolution algorithms and effector-regulatory immune signatures in TNBC demonstrated the existence of three immunoclusters. Among them, an immune-activated cluster was shown to be associated with an inflammatory phenotype, high T and natural killer (NK) cell infiltration, and a favorable prognosis. Interestingly, this immunocluster was also associated with enriched expression of CTLA-4 (13). The expression of CTLA-4 in various tumor types and cancer cell lines has been previously described (1417). It has also been reported that expression of CTLA-4 occurs in around 50% of breast carcinomas, but not in normal breast tissues (18). However, evidence describing CTLA-4 signal transduction pathways in TNBC cells, and the potential effects of its modulation, is still limited.

Thus, here we characterized the biological role of CTLA-4 expression in TNBC cells through preclinical in vitro experiments and benchmark analysis. Immunohistochemistry analyses with a robust scoring system of TNBC biopsies corroborated CTLA-4 expression in different cellular compartments. We then investigated CTLA-4 functions and associated signaling pathways by activating or blocking the receptor on TNBC cell lines. On the basis of public gene expression profiles of TNBC, the transcriptional landscape of tumors over-expressing CTLA-4 with activated downstream pathways was described. Additionally, we characterized the interactions between tumor-expressed CTLA-4 and immune infiltration. Finally, an overview of the possible clinical immunotherapy responses of tumors with activated CTLA-4-associated signaling was explored through public signatures. Improving our knowledge on the activity of CTLA-4 on tumor cells will help to understand the potential effects of the receptor on the clinical response to immunotherapy.

Materials and Methods

Clinical Samples and Cell Lines

A total of 50 patients diagnosed with invasive TNBC between 2005 and 2019, at the American British Cowdray Medical Center (ABC Medical Center) (Mexico City, Mexico) were recruited. The study was approved by the institutional research and ethics committees from the ABC Medical Center. Patients were selected if: (i) they were females; (ii) had histological diagnosis, (iii) had molecular diagnosis showing negative ER, PgR, and HER2; (iv) had electronic or physical clinical record to obtain clinical information on the stage of diagnosis and treatment; and (v) had a tumor percentage >10%. Paraffin-embedded tissue sections were retrospectively collected, and immunochemistry information, including ER, PgR, and HER2 expression, and Ki67 index were collected from the Pathology Department of the ABC Medical Center. The clinicopathological characteristics of the study population are summarized in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Clinicopathological characteristics of TNBC cases.

TNBC cell lines MDA-MB-231 (ATCC® HTB-26™), HCC1937 (ATCC® CRL-2336™), DU4475 (ATCC® HTB-123™), and HCC70 (ATCC® CRL-2315™) were purchased from the American Type Culture Collection (ATCC, Rockville, MD, USA). All cell lines were maintained in Roswell Park Memorial Institute (RPMI)-1640 medium (Biowest SAS, Nuaillé, France) supplemented with 10% fetal bovine serum (FBS; Biowest SAS, Nuaillé, France), 100 mg/ml of streptomycin, and 100 U/ml of penicillin (Invitrogen, Carlsbad, CA, USA).

Immunohistochemistry (IHC) and CTLA-4 Evaluation

CTLA-4 expression was evaluated by IHC in the selected TNBC biopsies. 2-μm sections from formalin-fixed paraffin-embedded tissue were deparaffinized and rehydrated. Antigen retrieval and immunodetection were achieved by using the Mouse/Rabbit Immuno-Detector DAB HRP Brown Detection System (Bio SB, Inc. CA, USA). Tissue sections were stained with anti-CTLA-4 (Cat. F-8 sc-376016, Santa Cruz Biotechnology, Inc. TX, USA) diluted 1:30. Standard hematoxylin-eosin (H-E) staining was performed for each tumor sample and CTLA-4 staining intensity was analyzed by an expert pathologist. Assessment of CTLA-4 expression on tumor cells was defined based on the criteria previously published for scoring PD-L1 (19). This protocol defines a positive stain as the percentage of membrane staining of any intensity in tumor cells, and tumor cells are scored as the proportion of tumor area. The scores were categorized as follows: T0 (negative, <1%), T1 (low/weak, >1 and <5%), T2 (medium, >5 and <50%) and T3 (high/strong, >50%) (Roche, VENTANA PD-L1). The slides were scanned using an iScan Coreo Au (Ventana Medical Systems, Inc, AZ, USA).

Reverse Transcriptase PCR Analysis

Total RNA was isolated from the cell lines using the RNeasy Mini Kit (Qiagen, Valencia, CA, USA) following the manufacturer's protocol. Reverse transcriptase-polymerase chain reaction (RT-PCR) for the detection of CTLA-4 transcripts was carried out using the SuperScript One-Step RT-PCR System (Qiagen, Valencia, CA, USA), with the following oligonucleotides: 5′-ATGGCTTGCCTTGGATTTCAGCGGCACAAGG-3′ and 5′-TCAATTGATGGGAATAAAATAAGGCTGAAATTGC-3′, which amplify a 660-bp fragment corresponding to the entire CTLA-4 coding sequence (20). The reverse transcription reaction was carried out at 50°C for 30 min. The CTLA-4 DNA amplification was performed using the following PCR program: 95°C for 15 min, followed by 40 cycles of 94°C for 1 min, 55°C for 1 min, 72°C for 1 min, and a final extension at 72°C for 10 min. Amplification of the β-actin gene, using the forward 5′-GGG TCA GAA GGA TTC CTA TG-3′ and reverse 5′-GGT CTC AAA CAT GAT CTG GG-3′ oligonucleotides, was included as an internal control. RNA from phytohemagglutinin-L-activated human peripheral blood mononuclear cells (PBMC) was included as a positive control.

Flow Cytometry Analysis

The expression of cytoplasmic CTLA-4 was analyzed by flow cytometry using the APC-conjugated mouse monoclonal anti-human CD152 antibody (Clone BNI3; BD Biosciences. San Jose, CA, USA). Cell surface CTLA-4 was evaluated by flow cytometry using mouse monoclonal Fluorochrome-labeled anti-CTLA-4 (Cat: 369612; Biolegend San Diego, CA, USA). CD80 was analyzed using the Phycoerythrin-conjugated anti-B7-1/CD80 (Clone No. 3771; R&D Systems Inc. Minneapolis, MN, USA). CD86 was evaluated using the APC-conjugated anti-B7-2/CD86 (Clone No. 37301; R&D Systems Inc. Minneapolis, MN, USA). Phytohemagglutinin-activated PBMC were included as a positive control for the expression of CTLA-4. Human peripheral blood monocytes cultured in vitro were included as a positive control for the expression of CD80 and CD86 (21). The cells were analyzed in a FACS Canto II Flow Cytometer (BD Biosciences Co. San Jose, CA, USA), capturing 10,000 events per sample. The percentages of positive cells and average fluorescence intensities were obtained and analyzed with the FlowJo 10 software.

Cell Proliferation and Invasion Assays

To evaluate cell proliferation, cells were seeded in 96-well plates at a density of 10,000 cells per well and incubated with either recombinant human CD80 (0.025, 0.15, and 1 μg/ml; Abcam, Cambridge, MA, USA) or Ipilimumab (1, 5, and 10 μg/ml; Bristol-Myers Squibb Company). Cell viability was measured at 24, 48, and 72 h after incubation using the colorimetric 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay.

To assess cell invasion, a transwell assay using extracellular matrix-coated Boyden chambers (Sigma-Aldrich, St. Louis, MO, USA) was performed. Cells were seeded in the upper chamber in FBS-free medium and treated with Ipilimumab (10 μg/ml) or CD80 (1 μg/ml) for 24 h. FBS-supplemented medium was added to the lower chamber. Cells which had passed through the matrix-coated membrane were recovered from the lower compartment, stained with the CellTracker Red reagent (Thermo Fisher Scientific Inc., Waltham, MA, USA) and evaluated in a Synergy H4 hybrid plate reader (BioTek Instruments Inc., Winooski, VT, USA) using the Gen5 software.

Western Blot

For protein extraction, cells were resuspended in lysis buffer (50 mM Tris–HCl, pH 7.4; 150 mM NaCl; 1 mM EDTA; 1% NP40; 0.25% sodium deoxycholate), containing 100 μl/ml of protease inhibitor cocktail (Roche, Manheim, Germany) and 10 μl/ml of phosphatase inhibitors (Sigma-Aldrich, St. Louis, MO, USA). Protein concentration was quantified using the DC protein assay kit (BioRad Laboratories, Hercules, CA, USA). Proteins (30 μg) were resolved by 10% SDS-PAGE and transferred onto polyvinylidene fluoride (PVDF) membranes (Millipore, Billerica, MA, USA). The following antibodies were used for the detection of CTLA-4-associated signaling pathways: Phospho-p44/42 MAPK (ERK1/2) (Thr202/Tyr204) (Cat. 9101), phospho-Akt (Thr308) (244F9) (Cat. 4056), phospho-Akt (Ser473) (Cat. 9271), Akt (pan) (C67E7) (Cat. 4691), Snail (C15D3) (Cat. 3879), ZEB1 (D80D3) (Cat.3396), and PTEN (138G6) (Cat. 9559), all from Cell Signaling (Cell Signaling Technology, Inc., Danvers, MA, USA); ERK1 and ERK2 (ERK1/2) (Cat. GTX17942), phospho-ERK1/2 (Thr185 + Thr202 + Tyr204 + Tyr187 (Cat. GTX24819); and anti-GAPDH (Cat. GTX100118), Vimentin (Cat. GTX-100619), and Twist (GTX-127310), all from GeneTex (GeneTex Inc., Irvine, CA, USA); E-cadherin (Cat. sc-71009), Vinculin (Cat. sc-25336), β-catenin (Cat. sc-7963), Actin (Cat. sc-47778), all from Santa Cruz (Santa Cruz Biotechnology, Inc., Dallas TX, USA). Densitometric analysis was performed with ImageJ software.

IL-2 Quantitative ELISA

Cells were incubated with Ipilimumab (10 μg/ml), CD80 (1 μg/ml), or left untreated for 48 h. Conditioned medium was harvested. Soluble IL-2 was determined using the Human IL-2 Mini TMB ELISA Development kit (PeproTech, Inc. NJ, USA) according to the manufacturer's protocol. Absorbance values were measured at 450 nm in a Synergy H4 hybrid plate reader (BioTek Instruments Inc. Winooski, VT, USA).

Co-culture of CTLA-4-Expressing TNBC Cells and Peripheral Blood Mononuclear Cells

Blood samples were obtained from healthy female donors. Peripheral blood mononuclear cells (PBMC) were isolated by standard Histopaque density centrifugation (Histopaque-1077. Sigma-Aldrich, St. Louis, MO, USA). PBMC were cultured in RPMI-1640 medium supplemented with 10% FBS, 100 mg/ml of streptomycin and 100 U/ml of penicillin, and 2 mM L-glutamine (Invitrogen, Carlsbad, CA, USA). To induce the expression of CTLA-4, PBMC were stimulated with 7.5 μg/ml Phytohemagglutinin (PHA-P. Sigma-Aldrich, St Louis, MO, USA) in supplemented RPMI for 6 days. CTLA-4 expression was evaluated by flow cytometry. Stimulated PBMC (16 × 104) were washed with PBS and co-cultured with HCC1937 or MDA-MB-231 (4 × 104) in 96-well plates for 48 h. The effect of CTLA-4 blockade was evaluated by adding Ipilimumab (10 μg/ml). At the end of the experiment, non-adherent cells were harvested and the number of viable cells was evaluated by MTT assay.

Mining Data and Microarray Data Processing

Publicly available data from 10 different tumor data sets, profiled on Affymetrix platforms deposited in the Gene Expression Omnibus (GEO); RNA sequencing data available at The Cancer Genome Atlas (GDC TCGA Breast Cancer) from the Xena browser (22); and normalized data from METABRIC study (23), recovered from cBioportal (24), were downloaded, and a total number of 1,320 TNBC tumors were evaluated. ER, PgR, and HER2 immunohistochemistry status was used to define the triple-negative phenotype. Supplementary Table 1 contains relevant characteristics of the data sets analyzed. All data retrieved from GEO were downloaded as raw data (cel file) and processed as follows: Microarray data from U133A, U133 Plus 2.0 and human transcriptome arrays were independently normalized and background corrected using RMA and quantile algorithms implemented in oligo package from Bioconductor (25). Genes from each platform were annotated with biomaRt (26), and only common cross-platform genes (including data from TCGA and METABRIC cohorts) were selected for further analysis. Duplicated probes were collapsed by selecting the probe with the highest interquartile range across the evaluated samples. Normalized data were scaled by median-absolute-deviation (MAD) for each sample. TCGA data were downloaded as raw counts and processed with limma-voom in limma R package (27).

Tumor Purity

Tumor purity was inferred according to the ESTIMATE algorithm (28). We applied the formula to the normalized MAD values and for limma-voom transformed data from TCGA tumors. Because some samples presented low tumor composition, we restricted our analysis to samples with an inferred tumor purity value over 60%.

Adjusted CTLA-4 Expression

To reduce bias due to remainder immune contribution to CTLA-4 expression levels, we adjusted it through a lineal regression model:

Adjusted CTL4 expression = lm(CTLA4~Est+CD4*CD8).

GSEA and ssGSEA

Non-preranked gene set enrichment analysis (GSEA) (29) was carried out via local software and Molecular Signatures collection (MSigDB: hallmarks, KEGG, Biocarta, Reactome, and Gene Ontology Biological Process Sets). GSEA was computed on the common cross-platform genes of the studied cohorts adjusted through median centering. Significantly enriched gene sets were reported based on NES and FDR value (<0.25). Single-sample GSEA (ssGSEA) (30) was performed with previously published gene signatures or gene sets from MSigDB, using GSVA Bioconductor package (31). For activated-like phenotype, we analyzed the MSigDB sets: for ERK1/2, (1) GO positive regulation of ERK1 and ERK2 cascade, (2) Reactome ERK MAPK targets, and (3) GO ERK1 and ERK2 cascade; for AKT pathway, (1) Hallmark PI3K/AKT/mTOR signaling, (2) Reactome AKT phosphorylated targets in the cytosol, and (3) Reactome AKT phosphorylated targets in the nucleus (Supplementary Table 2). The gene list of immune signatures was summarized by Charoentong et al. (32).

Immune Infiltration Characterization

The CIBERSORT web tool (https://cibersort.stanford.edu) was used to infer the relative proportions of immune cells using the leukocyte gene signature matrix (LM22) to deconvolve our MAD normalized datasets with 1,000 permutations (33). Relative proportion values were plotted in a heatmap. Additionally, xCell tool using default parameters (http://xcell.ucsf.edu/) was used to confirm immune infiltration composition and contribution of stem cells compartment. MAD-normalized gene expression data was analyzed. Normalized cell enrichment scores (ratio between the enrichment score of each sample and the maximum value among all samples) were evaluated.

Lehmann Subtypes

TNBC subtypes described by Lehmann et al. (34) were inferred via the web-based tool “TNBC type” (http://cbc.mc.vanderbilt.edu/tnbc/) based on MAD-normalized gene expression.

Characterization by Gene Signatures

Proliferation signature was recovered from the literature (proliferation Gene Cluster UNC337 intrinsic clustering) (35), and mean expressions of the multiple genes were computed as the score value. Cytolytic activity score has been defined as the mean of GZMA and PRF1 normalized gene expression (36). Tumor Inflammation Signature (TIS) score was defined by the average mean of normalized expression of the identified genes (37). Immune-related signatures were computed by the average gene expression value of each identified signature, as described previously (38, 39). To define responder patients based on anti-CLTA-4 immune scores, a support vector machine (SVM) approach was conducted on R. To train the model, data from TCGA TNBC tumors reported in Supplementary Table 3 (TNBC, N = 148) (39) was used. After the model construction, benchmark datasets were analyzed in R.

Statistical Analysis

Clinical data were analyzed using the SPSS™ Statistics Software (19.0 IBM Co, NY, USA). Differences in smoking history and treatment features between CTLA-4-expressing and CTLA-4-negative groups were analyzed by the Fisher's exact test. Chi-square test was applied to compare the expression of CTLA-4 among groups with different tumor stage at diagnosis. Results of in vitro experiments were presented as the mean ± standard error (SEM). T-test was used to compare the in vitro treatment groups and cell lines. Confidence intervals (CI; 95%) and p-values were calculated. Two-tailed p-values < 0.05 were considered statistically significant. Data were plotted using ggplot in R environment and GraphPad Prism software (GraphPad software). The statistical differences between the different conditions from benchmark TNBC data were assessed by Wilcoxon test taking HH condition as reference (*p-value ≤ 0.05; **p-value ≤ 0.01; ***p-value ≤ 0.001; ****p-value ≤ 0.0001). Fisher's exact test and Benjamin-Hochberg adjusted p-values were computed to examine the distributions of Lehmann subtypes and responder vs. non-responder groups.

Results

Tumor Cells of TNBC Patients Express CTLA-4 in Different Cell Compartments

A total of 50 TNBC cases were included. At diagnosis, 34% of patients had a tumor in early stage, 48% were locally advanced cases, and 14% presented metastatic disease. Complete clinicopathological data of the study population is summarized in Table 1. Expression of CTLA-4 was evaluated in tumor cells and infiltrating lymphocytes of TNBC paraffin-embedded biopsies, using an antibody that was previously validated in a breast cancer study (40). The quality of all tumor tissues was examined in H-E-stained slides (Figures 1A–D). All tumors showed cytoplasmic CTLA-4 staining, whereas different staining intensities were observed at the cell membrane. Overall, 15 tissue samples were negative (<1%) for CTLA-4 membrane staining, and they were classified as TC0 (Figure 1E). Surface CTLA-4 immunostaining was detected in 70% (35/50) of TNBC biopsies; 14% (7/50) of them showed a low/weak expression, scored as TC1 (>1 and <5%) (Figure 1F); 42% (21/50) showed a medium expression, scored as TC2 (>5 and <50%) (Figure 1G); and 14% (7/50) showed a high/strong membrane expression, scored as TC3 (>50%) (Figure 1H). Immunohistochemistry analysis also showed positive CTLA-4 immunostaining in lymphocytes of samples presenting tumor-infiltrating lymphocytes (TILs) (Table 1).

FIGURE 1
www.frontiersin.org

Figure 1. Immunohistochemical staining of CTLA-4 in TNBC tissues. (A–D) Representative images of hematoxylin-eosin staining and CTLA-4 membrane staining showing different intensity scores: (E) TC0 (<1%), (F) TC1 (>1 and <5%), (G) TC2 (>5 and <50%), (H) TC3 (>50%). All images are 20 × magnification.

Categorization of samples by tumor stage showed that the proportion of samples with expression of CTLA-4 at the cell-surface did not differ among patients diagnosed with early, locally advanced, and metastatic tumors (Fisher's Exact Test, p = 0.535) (Table 2). Expression of CTLA-4 was not significantly associated with history of smoking (p = 0.546) or treatment (neoadyuvant vs. non-neoadyuvant chemotherapy, p = 0.071).

TABLE 2
www.frontiersin.org

Table 2. Association of CTLA-4 expression with tumor stage, smoking history, and chemotherapy.

CTLA-4 and Its Primary Ligands CD80/CD86 Are Expressed on TNBC Cell Lines and Tumors

Our previous results demonstrated the expression of CTLA-4 on the membrane of TNBC tissue. Thus, we next analyzed the expression of CTLA-4 in TNBC-derived cell lines (DU4475, HCC70, HCC1937, and MDA-MB-231). The intracellular expression of CTLA-4 was confirmed by flow cytometry and was observed in all evaluated cell lines. The proportion of intracellular CTLA-4 positive cells ranged from 6.11 to 16.4%. The MDA-MB-231 cell line presented the highest proportion of intracellular CTLA-4 (16.4%) (Figure 2A). Since only CTLA-4 expressed on the cell surface is functional (41), we next analyzed membrane expression of CTLA-4. Membrane CTLA-4 was observed in all cell lines. The proportion of cells expressing membrane CTLA-4 was higher than 90% in HCC70, HCC1937, and MDA-MB-231 cell lines, while only 4.3% of DU4475 cells expressed CTLA-4 in this compartment (Figure 2B).

FIGURE 2
www.frontiersin.org

Figure 2. Expression of CTLA-4, CD80/CD86 ligands, and CTLA-4-associated pathways in TNBC cell lines. Expression of CTLA-4 protein was evaluated by flow cytometry using phytohemagglutinin-L-activated human lymphocytes (PBMC) as a positive control in (A) intracellular compartment and (B) cell-surface compartment. Expression of cell-surface (C) CD80 and (D) CD86 in HCC1937 and MDA-MB-231 cells by flow cytometry, including in vitro-cultivated human peripheral blood monocytes as a positive control. (E) The expression of the CTLA-4 gene was investigated by RT-PCR in the DU4457, HCC70, HCC1937, and MDA-MB-231 cell lines. Phytohemagglutinin-L-activated human lymphocyte RNA was included as a positive control (+), and RNA from in vitro-cultured human peripheral blood monocytes as a negative control (–). Lower panel: Representation of protein expression results (F) Scatterplot showing the correlation between CTLA-4-normalized expression and ssGSEA-enriched score of immune signatures (CD4/CD8, lymphocyte infiltration) among human triple-negative tumors from public databases (N = 1,320). CTLA-4 expression was highly concordant with lymphoid infiltration and tumor purity (Spearman correlation R = 0.72; p < 0.05). Colors of each point represent tumor purity values inferred by the ESTIMATE algorithm. (G) Donut charts showing the distribution of TN tumors over-expressing CTLA-4/CD80/CD86 (mRNA level). Color shades indicate high or low CD80 or CD86.

CTLA-4 activation occurs by interaction with its primary ligands CD80/CD86 on antigen-presenting cells (4244). Thus, we next investigated the expression of CD80 and CD86 on the surface of HCC1937 and MDA-MB-231 cells. Flow cytometry analysis showed that CD80 was weakly expressed in both cell lines (<13% of positive cells), compared to CD86, which was detected in 22.9% of HCC1937, and 44.3% of MDA-MB-231 (Figures 2C,D), suggesting a potential interaction of CD80/CD86-expressing cells with CTLA-4-expressing cells that might activate the receptor.

To extend our observations, we collected a comprehensive benchmark dataset evaluating gene expression profiles of 1,320 TNBC tumors. To assess a surrogate measurement of CTLA-4 protein expression based on mRNA levels, we first analyzed the presence of CTLA-4 mRNA in four TNBC-derived cell lines by RT-PCR. CTLA-4 mRNA was detected in all TNBC cells (Figure 2E). The amplified fragments observed were of the same size as those detected in activated human leukocytes, and a good relationship between CTLA-4 protein and mRNA level was observed. This suggests that it is suitable to use CTLA-4 gene expression as a surrogate value of protein status (Figure 2E).

Since a major source of CTLA-4 in TNBC tumors is the infiltrated lymphocyte immunome, we computed CD8/CD4 lymphocyte infiltration score, determined through ssGSEA of publicly available signatures (32) and compared it to CTLA-4 gene expression levels. Indeed, we observed a highly significant correlation of CTLA-4 gene expression and the CD8/CD4 enrichment score (N = 1,320, R = 72%, p < 0.05) (Figure 2F). To guarantee robust power to detect tumor CTLA-4 expression by assuring the highest tumor cell content, we selected tumors with at least 60% tumor purity, defined by the ESTIMATE algorithm (N = 783, Supplementary Table 1). Moreover, to minimize the bias due to the contribution of infiltrating lymphocytes, the reported CTLA-4 level was adjusted by the ssGSEA CD8/CD4 infiltration score and purity estimated value, ensuring a minimal impact of the immune compartment on the evaluated CTLA-4 expression. Notably, when looking at the proportions of TNBC over-expressing CTLA-4, 73% were found up-modulated in CD86 and 68% in CD80, whereas tumors with low levels of CTLA-4 generally presented a lower frequency of CD86 (34%) and CD80 (38%) up-modulation (Figure 2G). Taken together these observations indicate that the coordinated expression of the CTLA-4/CD80/CD86 axis is a widespread feature of TNBC.

Transcriptional Landscape of CTLA-4-Expressing Tumors Is Related to Relevant Oncogenic Pathways

It is known that activation of CTLA-4 on T-cells inhibits the PI3K/AKT pathway (4547). Thus, to investigate the effect of CTLA-4 on the activation of the PI3K/AKT cascade, HCC1937 and MDA-MB-231 cells were incubated with either recombinant CD80 or Ipilimumab. Activation of AKT by phosphorylation of the Thr308 or Ser473 residues was evaluated at different time points by Western blot. Our analysis showed constitutive phosphorylation of the Ser473 residue in HCC1937 cells, which was significantly increased after incubation with Ipilimumab, but experienced no alteration after incubation with CD80.

Conversely, incubation with CD80 induced a significant reduction of Ser473 phosphorylation in MDA-MB-231 cells, whereas Ipilimumab had no effect on the phosphorylation of this residue (Figure 3A). In sharp contrast, we observed that phosphorylation of the Thr308 was absent in both cell lines. Incubation with Ipilimumab resulted in the phosphorylation of Thr308 in HCC1937 cells but not in MDA-MB-231. No changes were observed after CD80 incubation. These observations suggest that blockade of CTLA-4 induces full activation of AKT through the phosphorylation of Ser473 and Thr308 in HCC1937 cells. It has also been described that CTLA-4 is able to activate the MAP kinase ERK1/2 pathway (48). Thus, we analyzed the effect of activating or blocking CTLA-4 on the phosphorylation of ERK1/2. As shown in Figure 3A, ERK1/2 is constitutively activated in MDA-MB-231 cells and phosphorylation was not affected by any of the treatments. In contrast, incubation of HCC1937 cells with CD80 induced the phosphorylation of ERK1/2 after 3 min of treatment. ERK1/2 remained phosphorylated for up to 15 min. When the cells were incubated with Ipilimumab, a modest phosphorylation of ERK1/2 was observed. However, it was only detected in early time points, at 3 and 5 min post-incubation. We demonstrated that CTLA-4 expressed on breast cancer cells, mainly in the membrane of HCC1937 cell line, is functional and induces the activation of the MAP kinase pathway, while AKT signaling is abrogated.

FIGURE 3
www.frontiersin.org

Figure 3. Transcriptional landscape of activated-like CTLA-4 in human tumors from publicly available data. (A) Evaluation of AKT and ERK1/2 activation signaling pathways after incubation with either CD80 (0.15 μg/ml) or Ipilimumab (10 μg/ml) at the indicated time points. Activation of AKT was evaluated by Western blot using antibodies against phosphorylated Ser473 (pAKT-Ser473) and Thr308 (pAKT-Thr308) residues. Activation of ERK1/2 was analyzed by using an antibody that recognizes phosphorylated Thr185 + Thr202 + Tyr204 + Tyr187 residues (pERK1/2). Expression of GAPDH was included as loading control. Densitometric analysis indicated as percentage below each blot, taking time 0 as reference for each evaluated protein. (B) Heatmap of decoding tumor subgroups based on CTLA-4 expression and activation. First panel shows the inferred tumor purity value (%) among the complete TNBC cohort (N = 1,320). The following panels display selected tumors based on tumor purity above 60% (N = 783). Middle panel shows tumor subtype features of CTLA-4 ERK1-2 tumors, while lower panel shows features of CTLA-4 AKT tumors. Samples represented in the plot were sorted by tumor purity value in each cluster. CTLA-4 expression value was adjusted by CD8/CD4/lymphocyte infiltration and ESTIMATE value. Boxplots of CD80 and CD86 normalized expression values in (C) CTLA-4/ERK cluster and (D) CTLA-4/AKT cluster (HH: CTLA-4H/H−ERK1/2 and CTLA-4H/H−AKT; HL: CTLA-4H/L−ERK1/2 and CTLA-4H/L−AKT; LH: CTLA-4L/H−ERK1/2 and CTLA-4L/H−AKT; LL: CTLA-4L/L−ERK1/2 and CTLA-4L/L−AKT). (*p-value ≤ 0.05; ***p-value ≤ 0.001; ****p-value ≤ 0.0001).

To study the role of CTLA-4 on TNBC tumors, we explored public data with a tumor purity >60%. Tumors showing a CTLA-4 expression values over the median were assumed to be CTLA-4High (CTLA-4H). Then, to infer CTLA-4 activation by cell-membrane expression, the results of in vitro experiments were considered, and tumors were subjected to ssGSEA analysis using the enrichment scores of ERK1/2 and its target gene set, as well as PI3K/AKT/mTOR and AKT target gene sets (nucleus and cytoplasm). Four major groups were defined for each oncogenic pathway. The resulting clusters were described based on the ERK1/2 pathway as follows: tumors that are CTLA-4H/H−ERK1/2 (N = 193), representing the biological relevant group for our study where tumor cells over-express CTLA-4 and are activated-like tumor as shown by ERK1/2 enrichment pathway; CTLA-4H/L−ERK1/2 (N = 125), which probably do not express CTLA-4 at the membrane level; and two low CTLA-4 expression groups, CTLA-4L/H−ERK1/2 (N = 203) and CTLA-4L/L−ERK1/2 (N = 226) (Figure 3B). Based on the AKT pathway, tumors were categorized as: CTLA-4H/H−AKT (N = 156), which do not express membrane CTLA-4, and CTLA-4H/L−AKT, representing our study phenotype where CTLA-4 is up-regulated in tumor cells and is activated-like, as defined by the down-modulation of the AKT pathway. Finally, CTLA-4L/H−AKT and CTLA-4L/L−AKT were also detected (Figure 3B).

To explore whether expression of the CTLA-4/CD80/CD86 axis was associated with an activated-like phenotype, gene expression of CD80/CD86 was evaluated on TNBC. CTLA-4H/H−ERK1/2 tumors presented upregulation of CD80 (Figure 3C), but a more significant over-expression of CD86 compared to CTLA-4H/L−ERK1/2 samples (Figure 3D). This observation recapitulates the results of our in vitro assays and suggests the presence of active CTLA-4/CD80/CD86 paracrine signaling.

To get insight into the transcriptional landscape of the categorized CTLA-4 tumors, a GSEA analysis was performed based on common cross-platform genes evaluated by the public expression profiles of TNBC tumors. Several potentially actionable molecular alterations were enriched in CTLA-4-over-expressing activated-like tumors. We initially focused on the comparison between ERK-activated tumors (CTLA-4H/H−ERK1/2 vs. CTLA-4H/L−ERK1/2). As expected, GSEA denoted the enrichment of ERK1/2 down-stream signaling pathways such as TGF-β, multiple signal transduction processes, such as tumor necrosis factor α (TNFα), extracellular-signal-regulated kinases (ERKs), cJun kinases (JNKs), and GTPase activity in CTLA-4H/H−ERK1/2 tumors. GSEA also pointed out pathways activated by TNFα including apoptosis, NF-κB, and MAP kinases (Figure 4A and Supplementary Table 3). Other enriched processes in this tumor phenotype included KRAS, a critical regulator of oncogenic properties, cell–cell adhesion, androgen feedback, and stress responses to hypoxia and UV exposure (Figure 4A and Supplementary Table 3). Hypoxic conditions result in the transcription of various cytokines, including VEGF, which is a pro-angiogenic factor and was enriched in CTLA-4H/H−ERK1/2 tumors. Some pathways overlap with CTLA-4-activated-like tumors with an abrogated AKT program (CTLA-4H/L−AKT vs. CTLA-4H/H−AKT), such as KRAS signaling and stem-like related pathways including Notch and Hedgehog. GSEA also revealed unique significant enrichment of gene sets related to drug metabolism and regulation of ion and amino acid transport (Figure 4B).

FIGURE 4
www.frontiersin.org

Figure 4. Biological effects of CTLA-4 activation in human tumors and in vitro cell line models. Results from GSEA as bubble plot comparing (A) CTLA-4H/H−ERK1/2 vs. CTLA-4H/L−ERK1/2 and (B) CTLA-4H/L−ERK1/2 vs. CTLA-4H/H−ERK1/2, where the size of the circles gives the normalized enrichment score (NES) and the color Indicates FDR. In vitro assays to evaluate the effect of the activation and blockade of CTLA-4 on cell proliferation. (C) The effect of blocking CTLA-4 on the proliferation of HCC1937 and MDA-MB-231 cells was analyzed by incubating cells with the monoclonal antibody Ipilimumab (10.0 μg/ml), cell number was evaluated at 24 and 48 h post-incubation. (D) The effect of activation of CTLA-4 on the proliferation of HCC1937 and MDA-MB-231 cells was analyzed by incubating cells with recombinant human CD80 (0.15 μg/ml), cell number was evaluated at 24 and 48 h post-incubation. Boxplots for the expression of proliferation gene signatures in (E) CTLA-4/ERK1-2 and (F) CTLA-4/AKT cluster. (*p-value ≤ 0.05; ***p-value ≤ 0.001; ****p-value ≤ 0.0001.

As a negative regulator of T cell activity, alterations of AKT and ERK1/2 cascades induced by CTLA-4 have been associated with the inhibition of cell proliferation (49, 50). Therefore, we investigated whether incubation of CTLA-4-expressing cells with various concentrations of recombinant CD80 or with Ipilimumab might regulate proliferation rates at different time points. As shown in Figure 4C, blockade of CTLA-4 by incubation with Ipilimumab induced a dose-dependent reduction of HCC1937 and MDA-MB-231 cell numbers 48 h post-treatment (p < 0.01), but not in HCC70 cells. CD80 incubation only modified the number of cells at early time points (24 h) in MDA-MB-231, but not at 48 h (Figure 4D). Incubation with CD80 had no effect on the proliferation of HCC70 and HCC1937 cells. Further inspection of proliferation among TNBCs revealed that CTLA-4-activated-like tumors present lower expression of proliferation gene clusters in comparison with tumors that only over-expressed the immune checkpoint (Figures 4E,F). Although the GSEA analysis indicated the significant enrichment of potential cell adhesion regulators, which may impact cell motility/invasion, and could be associated with mesenchymal features, we were not able to demonstrate an effect of CTLA-4 modulation on cell invasion (Supplementary Figure 1A) or on the expression of epithelial-to-mesenchymal biomarkers, in vitro (Supplementary Figure 1B).

Together, these results revealed a unique transcriptomic landscape of CTLA-4-activated-like tumors that might impact relevant cancer programs locally, in the tumor cell, and in the tumor microenvironment.

CTLA-4-Activated Tumors Are Enriched in Immunogenic Pathways and CTLA-4 Blockade Increases IL-2 Secretion

Rationalizing that tumor expression of CTLA-4 not only regulates cellular functions, but may also modulate the immune microenvironment, we inspected the immunogenicity characteristics of cell lines and tumors. First, GSEA revealed an enrichment of immune-related pathways only in the CLTA-4-activated-like phenotype, inferred by ERK status (Figure 5A and Supplementary Table 3). This analysis showed that CTLA-4-activated tumors are highly enriched in immune cells (T cells, NK cells, neutrophils) and immune pathways (immune activation, interleukin and immune communication signaling) with relevant regulatory functions.

FIGURE 5
www.frontiersin.org

Figure 5. Immunophenotypes of TNBC with activated CTLA-4. (A) Bubble plot showing the most significant enriched immunological gene sets from GSEA in CTLA-4H/HERK1/2 compared to CTLA-4H/LERK1/2 tumors. Bubble size represents normalized enrichment score (NES) and color FDR value (<0.05). Box plots of (B) cytolytic activity (CYT) and (C) Tumor Inflammation Signature (TIS) according to gene signatures among TNBC (HH: CTLA-4H/H-ERK1/2; HL: CTLA-4H/L-ERK1/2; LH: CTLA-4L/H-ERK1/2; LL: CTLA-4L/L-ERK1/2). Heatmap of relative proportions of immune subpopulations computed with CIBERSORT in (D) CTLA-4H/H-ERK1/2 and (E) CTLA-4H/L-ERK1/2 tumors. Samples are sorted according to the relative proportion of each immune cell (Barplots in upper panel). (F) Heatmap of significant enriched ssGSEA scores (FDR < 0.05) of immune cell subpopulations according to immune signatures described by Charoentong in CTLA-4H/H-ERK1/2 TNBC compared to CTLA-4H/L-ERK1/2. Upper: adaptive immunity, middle: adaptive-innate immunity, and lower: innate immunity. Boxplots showing the xCell scores of (G) Lymphoid population scores and (H) Myeloid population scores differently distributed (FDR < 0.05) within HH (CTLA-4H/H-ERK1/2) and HL (CTLA-4H/L-ERK1/2) tumors. (I) Secretion of IL-2 by HCC1937 and MDA-MB-231 was evaluated by quantitative ELISA after incubation with CD80 (0.15 μg/ml) or Ipilimumab (10 μg/ml). (J) Effect of CTLA-4 blockade on lymphocytes proliferation (phytohemagglutinin-stimulated peripheral blood leukocytes) in co-culture with HCC1937 and MDA-MB-231 cells in the presence and absence of Ipilimumab. (*p-value ≤ 0.05; ****p-value ≤ 0.0001).

Considering that physiological anti-tumor immunity requires a cytolytic immune response, cytolytic activity scores (CYT) (36) were computed in TNBC datasets. CTLA-4H/H−ERK1/2 tumors exhibited significantly higher CYT scores compared to CTLA-4H/L−ERK1/2 (p = 0.0001) (Figure 5B). Furthermore, a pre-existing adaptive immune response within tumors was inferred via the Tumor Inflammation Signature (TIS) algorithm (37); a similar enriched profile was found for CTLA-4H/H−ERK1/2 tumors, which presented the highest TIS score (p = 0.0001) (Figure 5C).

To define the crosstalk between the activated-CTLA-4 axis in tumor cells and the intra-tumor immune infiltrate, we estimated the immune cellular composition of the analyzed TNBC tumors based on their transcriptome profiles, through dedicated algorithms. First, we estimated the relative immune cellular composition using the CIBERSORT deconvolution method (33). We observed that the frequency of the immunological sub-populations appears to be distinct for each tumor subtype (Figures 5D,E and Supplementary Table 4). We then validated the predicted immune composition based on benchmark signatures (32) by ssGSEA analysis. We identified 16 kinds of immune cell subtypes that had different enrichment scores between CTLA-4H/H−ERK1/2 and CTLA-4H/L−ERK1/2. Five adaptive immune cells were significantly associated with the CTLA-4H/H−ERK1/2 phenotype (Figure 5F), and these tumors also showed a significant up-regulation of three classes of NK cells and two dendritic cells. A high degree of consistency with immune scores estimated by xCell (51) in CTLA-4H/H−ERK1/2 tumors vs. CTLA-4H/L−ERK1/2 was observed (Figures 5G,H and Supplementary Table 4). xCell analysis also pointed out the loss of stem cell compartment, particularly the mesenchymal stem cell (MSC) population and common lymphoid progenitors (CLP) on CTLA-4-activated-like tissues (Supplementary Figure 2A).

In agreement with the above described results, ERK-activated-like tumors over-expressing CTLA-4 (CTLA-4H/H−ERK1/2) were enriched in immune (IM) and basal-like 1 (BL1) subtypes, as described by the TNBC molecular classification defined by Lehmann (34), in comparison to down-modulating CTLA-4 tumors [Fisher's exact, Benjamin–Hochberg (BH) p < 0.05] (Supplementary Figure 2B). Additionally, no differential enrichment was observed in IM or BL1 subtype distribution in CTLA-4H/L−AKT and CTLA-4H/H−AKT tissues, corroborating the results computed through the GSEA (Supplementary Figure 2C).

A major effect of CTLA-4 activation on lymphocytes is the inhibition of interleukin IL-2 production and secretion (49, 50). Thus, to address the question of whether CTLA-4 modifies IL-2 secretion by TNBC cells, the concentration of soluble IL-2 in culture medium of untreated and cells treated with either recombinant CD80 or Ipilimumab was evaluated by ELISA. Incubation with CD80 significantly reduced IL-2 secretion by HCC1937 cells, but had no effect on MDA-MB-231 cells (Figure 5I). In contrast, incubation with Ipilimumab markedly enhanced IL-2 secretion on HCC1937 cells (Figure 5I), while MDA-MB-231 cells were unaffected (Figure 5I). These observations suggest the potential effects of CTLA-4 activation on IL-2 secretion in some TNBC cells.

Since IL-2 is a critical stimulator of T cell proliferation and differentiation, we next investigated whether proliferation of lymphocytes can be modulated by TNBC cell lines secreting IL-2. Co-culture of lymphocytes with tumor cells induced a significant increment on the number of viable lymphocytes (p < 0.05), which may be explained by the natural production of IL-2 by both cell lines (Figure 5J). Further, blocking CTLA-4 with Ipilimumab increased the number of viable lymphocytes compared to non-treated cells (p < 0.05). The most significant effect was observed in lymphocytes co-cultured with HCC1937, which is consistent with the observation that Ipilimumab boosts IL-2 secretion in this cell line (Figure 5J). Together, these observations suggest that CTLA-4 blockade in tumor cells may contribute to generating a favorable microenvironment for immune responses.

Although IL-2 plays an important role in the expansion of cytotoxic T cells, recent studies have demonstrated that IL-2 is also crucial for the establishment and maintenance of T regulatory cells (Tregs) (52, 53). Thus, we explored the presence of Tregs in CTLA-4-expressing TNBC. Interestingly, application of deconvolution methods showed that tumors of the CTLA-4-activated subtype displayed a greater range of CD8/CD4 signatures, but we did not observe a differential enrichment of Tregs between the CTLA-4H/H−ERK1/2 and CTLA-4H/L−ERK1/2 tissues (Supplementary Figure 2D). Accordingly, the expression of FOXP3, a regulatory element necessary for Treg development (54), did not vary across CTLA-4 over-expressing tumors (Supplementary Figure 2E) and a similar correlation was detected between CD8/CD4 enrichment score and FOXP3 in CTLA-4H/H−ERK1/2 and CTLA-4H/L−ERK1/2 tumors (Supplementary Figure 2F).

Moreover, our results showed the presence of NK cells in CTLA-4H/H−ERK1/2 tumors. NK cell proliferation is regulated by IL-15 (55), but studies of IL-15 expression in breast cancer are scarce (56). Thus, we analyzed the level of IL-15 mRNA in tumors. The level of IL-15 mRNA in CTLA-4H/H−ERK1/2 tumors was significantly higher than the observed in CTLA-4H/L−ERK1/2 tumors (p < 0.05) (Supplementary Figure 3A). Likewise, immunogenomic analysis demonstrated a significant correlation between IL-15 mRNA expression and the NK fraction (R = 0.34, p < 0.05) (Supplementary Figure 3B). IL-15 is mainly produced by macrophages after stimulation with Colony Stimulating Factor-1 (CSF-1). Accordingly, we observed that CSF-1 is highly expressed in CTLA-4H/H−ERK1/2 tumors (Supplementary Figure 3C).

Overall, characterization of the tumor infiltrated microenvironment supports the idea that there are different patterns of tumor immune response associated with CTLA-4 expression levels and activated-like status in tumor-cells.

CTLA-4 Activated Tumors Present Highly Predicted Response Scores to Immunotherapy Based on Immune Signature

Enhancing our understanding of the interaction between cancer cells and the immune system would lead to the improvement of immunotherapy for treating cancer (5759). Thus, we aimed to assess the potential benefit of immunotherapy across TNBC benchmark data by interrogating their transcriptional profiles. We applied the Immune Score (IS) (39), which infers the response to anti-CTLA-4 immunotherapies in accordance with tumor CTLA-4 status. ISs were more elevated in CTLA-4H/H−ERK1/2 than in CTLA-4H/L−ERK1/2 tissues (Figure 6A). In accordance, the proportion of potential responders to CTLA-4 immunotherapy were highly enriched in CTLA-4H/H−ERK1/2 tumors (61% in CTLA-4H/H−ERK1/2 vs. 40% in CTLA-4H/L−ERK1/2, BH p = 0.031) (Figure 6B). To further dissect the biological role of the described immune phenotype, we assessed the tissue-resident memory T gene signature (TRM) (38), which described high enrichment of T cell populations in TNBCs and is able to distinguish patients with high CD8 expression who have a good or poor prognosis. TRM signature scores were higher in activated-CTLA-4 tissues (Figure 6C). Taken together, these results suggest that tumor cells over-expressing CTLA-4 may show distinct clinical and therapy responses as a result of the established biological signaling pathways.

FIGURE 6
www.frontiersin.org

Figure 6. Immune signatures associated with immunotherapy response and clinical behavior in TNBC with diverse CTLA-4 portrait. (A) Boxplot of CTLA-4 immune signature scores according to tumor subtypes (HH: CTLA-4H/H−ERK1/2 and HL: CTLA-4H/L−ERK1/2). (B) Proportion of potential responders and non-responders to anti-CTLA-4 immunotherapy is shown according to tumor subtypes (Fisher BH p-value < 0.05). (C) Boxplot showing the distribution of CD8 TRM signature scores among TNBC stratified by CTLA-4 subtypes (HH: CTLA-4H/H−ERK1/2 and HL: CTLA-4H/L−ERK1/2) (****p-value ≤ 0.0001). (D) Graphical representation of the biological landscape of TNBC expressing CTLA-4. Immunohistochemistry (N = 50), preclinical in vitro evaluation (MDA-MB-231 and HCC1937) and in silico analysis of publicly available gene expression data on TNBC (N = 783).

Discussion

CTLA-4 was initially described as a T cell-associated molecule. However, there is evidence showing the expression of CTLA-4 in various types of non-lymphoid cells (60). Here we describe the cytoplasmic and cell surface expression of CTLA-4 on TNBC samples and cell lines. Expression of CTLA-4 on breast cancer cells was first reported by Contardi et al. (15) in five ductal carcinoma patients with positive immunostaining at the membrane and cytoplasmic compartments, with similar intensity in all samples. Thereafter, the expression of CTLA-4 at protein and mRNA levels was demonstrated in breast cancer (60). More recently, CTLA-4 has been reported to be over-expressed in a high proportion of breast cancer samples (18). In contrast with former reports, we exclusively focused our study on the expression of CTLA-4 in TNBC. We observed that all analyzed TNBCs (N = 50) showed cytoplasmic expression of CTLA-4, as previously reported (61). However, differences in membrane expression were detected. Around 70% of tissues were positive for CTLA-4 at the membrane showing different levels of staining. Considering that only surface-expressed CTLA-4 is functional (41), these findings might have potential clinical implications. Unfortunately, the small number of patients included in the present work prevented us from correlating CTLA-4 surface expression and clinical features. Thus, evaluation of CTLA-4 expression in a larger TNBC population is necessary.

In agreement with earlier publications, we found that various TNBC-derived cell lines express CTLA-4 at the cell surface (14, 15). CTLA-4 is activated after binding to the CD80/CD86 ligands on the surface of antigen-presenting cells (62). Here, we detected the expression of CD80 and CD86 on the cell membrane of TNBC cell lines and human tumors, suggesting that this paracrine signal triggers CTLA-4 activation and its associated pathways.

Decoding the biological states of cell lines expressing CTLA-4 at the cell surface allowed us to detect a constitutive phosphorylation of the AKT-Ser473 residue in all cell lines, which was increased by Ipilimumab in HCC1937, but not in MDA-MB-231 cells. Surprisingly, we were not able to detect phosphorylation of the Thr308 residue in unstimulated cells. Nonetheless, incubation with Ipilimumab induced Thr308 phosphorylation in HCC1937 cells. In order to acquire full kinase activity, AKT must be phosphorylated in Ser473 at the C-terminal motif site and in Thr308 residue, which is located at the catalytic site (63, 64). Thus, it seems that CTLA-4 blockade may induce full activation of AKT signaling in particular cells, like HCC1937. Interestingly, AKT is a downstream target of PP2A, a crucial serine-threonine phosphatase, which associates to CTLA-4 and regulates its inhibitory activity (65), and reduced PP2A activity has been related with enhanced Thr308 phosphorylation, without altering Ser473, in acute myeloid leukemia (66). It was also reported that over-expression of CTLA-4 promotes the activation of the ERK1/2/MAPK pathway in different tumor types (67, 68) and stimulated T-cells (69). Additionally, coordinated anti-CTLA-4 and anti-PD-1 treatment suppressed the protein levels of phosphorylated ERK1/2 (70). Our results demonstrate that CD80 treatment induces the phosphorylation of ERK1/2 in HCC1937 cells at short post-treatment times.

We also systematically evaluated the transcriptional landscapes of publicly available gene expression profiles in accordance with CTLA-4 levels and their inferred activation statuses. Different contingency strategies to minimize bias due to the immune infiltrated compartment were applied to evaluate CTLA-4 mRNA portraits. We only studied samples with a predicted tumor purity higher than 60%, and CTLA-4 expression profiles were adjusted using the lymphocyte enrichment scores and tumor purity values. Our study decodes the heterogenous tumor clusters guided by CTLA-4 expression and activation and suggests the existence of diverse tumorigenic pathways in TNBC in association with CTLA-4 status. When the biological groups were stratified in accordance with their activated-CTLA-4 expression, we found that MAPK signaling, TGF-ß, VEGF, TNF-α, drug metabolism, ion and amino acid transport, and KRAS signaling were significantly associated with a CTLA-4 activated phenotype. Together, these observations indicate that CTLA-4 induces tumor driver pathways by distinct cellular mechanisms, primarily differing in ERK1/2 and AKT signaling processes (Figure 6D).

It is generally accepted that TNBC are tumors with a high rate of proliferation (71). Under this perspective, we analyzed the potential effect of CTLA-4 on proliferation rates. Molecular characterization of CLTA-4-activated-like tumors did not show differences in the expression of proliferation genes compared to tumors that only over-express CTLA-4. This could be partially explained by the reduced prevalence of Basal-like 1 subtype in the activated-like lesions, which are heavily enriched in annotations associated with proliferation (34). On the other hand, CTLA-4 blockade with Ipilimumab decreased total cell number of all treated cell lines. These results were unexpected because inhibition of CTLA-4 is traditionally associated with an increase in T lymphocyte proliferation. However, the activity of CTLA-4 in different cell types is incompletely understood. In line with our results, it was reported that treatment with Ipilimumab significantly reduced proliferation of CTLA-4-expressing CD4+ T-cells (71). Furthermore, the effect of CTLA-4 blockade seems to be associated with a decreased production of cytokines, such as IFN-γ, TNF-α, IL-13, and IL-17F, rather than with the induction of apoptosis (72). Although we did not evaluate apoptosis markers, differences in cell viability were not observed between cells incubated with Ipilimumab and untreated cells, suggesting that the reduction of proliferation might be associated with unknown signaling cascades or with extracellular signals, such as cytokines, that need to be studied to better understand the molecular mechanisms governing CTLA-4 functions in TNBC cells. Importantly, anti-proliferative effects of anti-CTLA therapies in tumor cells has novel prognostic and clinical implications.

We observed that the immune component was up-regulated in the CTLA-4 ERK1/2 activated-like subtype. GSEA screening clearly delineated enrichment of inflammation and immune infiltration signaling, as well as the antigen processing/presentation machinery in association to CTLA-4 activity. Interestingly, translocation of CTLA-4 to the T-cell membrane induces the inhibition of IL-2 production and secretion (72). The aberrant production of IL-2 has been documented in several tumors (7376), including breast cancer (77). Accordingly, incubation with CD80 decreased the concentration of soluble IL-2 in HCC1937 cells, whereas treatment with Ipilimumab significantly increased soluble IL-2, which was able to sustain the proliferation of lymphocytes already expressing CTLA-4. It is known that CTLA-4 regulates transcription of IL-2 mainly by suppressing the activity of the PI3K/AKT cascade (78). Thus, we hypothesize that the reduction of secreted IL-2 might be due to either the inhibition of AKT-Thr308 phosphorylation or the over-activation of ERK/TGF-β axis, since it has been reported that TGF-β inhibits effector/memory T lymphoblast proliferation and IL-2 production (79, 80).

Our data suggest that Ipilimumab may stimulate CTLA-4-expressing tumor cells to secrete IL-2, contributing to the establishment of a favorable immunologic microenvironment, that may support the action of Ipilimumab. In our study, CTLA-4 expression and activation were correlated with higher immune infiltration of cells typically expressing CTLA-4 as an evasion mechanism (8184), such as exhausted-like lymphoid cells, mainly CD8/CD4 T-cells and NK cells, as well as the myeloid compartment, including macrophages. Tumor microenvironment may influence immunotherapy responses by building complex interactions between tumor cells, immune cells, and secreted soluble factors. We described that CTLA-4 might regulate tumor immune microenvironment by enhancing relevant immunomodulators, such as IL-15 and CSF-1, which prime NK and M1 macrophages, two cells interconnected by their immunomodulatory functions. These data point out the existence of a fine crosstalk between immune infiltrated cells and CTLA-4 expressing tumor cells, that deserves further study.

At present, the possibility of including immune checkpoint inhibitors as part of primary TNBC treatment is being discussed (8587). However, the lack of biomarkers for the selection of patients still limits the use of immune checkpoint blockade in TNBC therapy. Here we have found that, in addition to the expected expression of CTLA-4 on infiltrating lymphocytes, tumor cells from some patients, and cell lines, are also positive for the expression of CTLA-4 at the cell-surface. Interestingly, blockade of the receptor did not produce the same effect in all CTLA-4-expressing cell lines, and it seems to depend on the potential modulation of the canonical signaling cascades (AKT, ERK1/2) and the resulting production of regulatory cytokines, such as IL-2, suggesting that checkpoint inhibitors would work best in patients with tumors over-expressing CTLA-4 at the cell membrane. In fact, the immune signature score, which is well associated with response to anti-CTLA-4 treatment, revealed a high prevalence of predicted responders to immunotherapy in CLTA-4 activated patients in the analyzed datasets. Overall, these findings provide novel insights into the biological landscape of tumors over-expressing CTLA-4 and could lead to the identification of useful biomarkers for immunotherapy, increasing the spectrum of patients to be benefited.

An important limitation of this work is the small number of samples evaluated by immunochemistry, which prevented us from making any generalization of the results. To overcome this limitation, a larger number of tissue samples are currently being studied. Another restriction is the small number of cell lines analyzed that do not completely recapitulate TNBC heterogeneity. However, taking advantage of benchmark data we reinforced our analysis including 730 triple-negative tumors fulfilling the inclusion criteria. This allowed us to identify different biological patterns associated with CTLA-expression in tumors and in cell lines. We also defined relevant mechanisms and biological portraits of human tumors over-expressing activated CTLA-4, which may predict favorable responses to anti-CTLA-4 therapy.

In conclusion, we have found that CTLA-4 is expressed at the cell membrane of a subgroup of TNBC tumors and our results provide novel mechanistic data by which CTLA-4 regulates key signal transducers in tumor cells (Figure 6D). Moreover, blockade of CTLA-4 with Ipilimumab significantly activated molecular cascades that probably cooperate to create a favorable microenvironment that enhances immune responses against tumor cells (Figure 6D). These observations are promising since they suggest that CTLA-4 expressed at the tumor-cell surface might be itself a target of immune checkpoint inhibitors, and might also be a candidate biomarker for immunotherapy, which is urgently needed to provide robust tools to select patients who might benefit from this new therapeutic strategy.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: the NCBI Gene Expression Omnibus (GSE25066, GSE2603, GSE19615, GSE21653, GSE76275, GSE102484, GSE76250, and GSE86948), the UCSC Xena Browser (https://xenabrowser.net/) (GDC TCGA Breast Cancer, BRCA), and cBioPortal (https://www.cbioportal.org/) (METABRIC).

Ethics Statement

The studies involving human participants were reviewed and approved by American British Cowdray Medical Center (ABC Medical Center), Mexico City, Mexico. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

LR-Z conceived the project and designed experiments, supervised laboratory staff, analyzed and discussed data, and wrote the manuscript. SR-C conceived and conducted all genomic data analyses, including preprocessing of raw microarrays or RNA-seq data, enrichment analysis and public gene signatures, analyzed and discussed data, and wrote the manuscript. MN-B performed cell proliferation and invasion assays, IL-2 ELISA, flow cytometry, and Western blot analysis. MC-B identified TNBC cases, created the patient database, performed immunohistochemistry, evaluated CTLA-4 staining, analyzed data, and wrote manuscript. JM-H supervised patient data revision, performed statistical analysis regarding CTLA-4 scores, and clinical features. CL-T supervised immunohistochemistry, performed pathological assessment, and CTLA-4 scores. RG-C and AZ-D coordinated study operation at the ABC Medical Center and participated in discussion and analysis of results. MI-S and JE-L supervised and performed Western blot analysis. JM and VC-M performed flow cytometry analysis. DO-P, DV-O, and ER-S participated in the selection of patients, retrieving of specimens, and reviewing of patient's data. PS-S and LT-U performed co-culture assays, analyzed data, and produced figures. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National University of Mexico (UNAM. PAPIIT. IN204419), Consejo Nacional de Ciencia y Tecnología (CONACyT. SALUD-2016-1-272793); DGAPA (Fellowship 2017 to MC-B), Programa Institucional de Cáncer de Mama, Instituto de Investigaciones Biomédicas, UNAM (Fellowship to MC-B), and Cátedra Salvador Zubirán, UNAM/INCMNSZ (to SR-C).

Conflict of Interest

LR-Z has received personal fees from Bristol-Myers Squibb and Roche.

The remaining authors declare that the research was conducted in the absence of any commercial of financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank Ma. Cecilia Aguilar Zacarías, Diana E. Aguilar León, María de los Angeles Ibarra Meneses, Laura García Alanís, and José Rafael Cervantes Roldán for excellent technical assistance and the Flow Cytometry core facility from Coodinación de Investigación en Salud at Centro Médico Nacional, Siglo XXI of IMSS for instrumentation.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.01206/full#supplementary-material

References

1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2018) 68:394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Andreopoulou E, Schweber SJ, Sparano JA, McDaid HM. Therapies for triple negative breast cancer. Expert Opin Pharmacother. (2015) 16:983–98. doi: 10.1517/14656566.2015.1032246

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Reynoso-Noveron N, Villarreal-Garza C, Soto-Perez-de-Celis E, Arce-Salinas C, Matus-Santos J, Ramirez-Ugalde MT, et al. Clinical and epidemiological profile of breast cancer in Mexico: results of the seguro popular. J Glob Oncol. (2017) 3:757–64. doi: 10.1200/JGO.2016.007377

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Stead LA, Lash TL, Sobieraj JE, Chi DD, Westrup JL, Charlot M, et al. Triple-negative breast cancers are increased in black women regardless of age or body mass index. Breast Cancer Res. (2009) 11:R18. doi: 10.1186/bcr2242

CrossRef Full Text | Google Scholar

5. Lebert JM, Lester R, Powell E, Seal M, McCarthy J. Advances in the systemic treatment of triple-negative breast cancer. Curr Oncol. (2018) 25:S142–50. doi: 10.3747/co.25.3954

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Yao H, He G, Yan S, Chen C, Song L, Rosol TJ, et al. Triple-negative breast cancer: is there a treatment on the horizon? Oncotarget. (2017) 8:1913–24. doi: 10.18632/oncotarget.12284

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Qureshi OS, Zheng Y, Nakamura K, Attridge K, Manzotti C, Schmidt EM, et al. Trans-endocytosis of CD80 and CD86: a molecular basis for the cell-extrinsic function of CTLA-4. Science. (2011) 332:600–3. doi: 10.1126/science.1202947

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Robert C, Thomas L, Bondarenko I, O'Day S, Weber J, Garbe C, et al. Ipilimumab plus dacarbazine for previously untreated metastatic melanoma. N Engl J Med. (2011) 364:2517–26. doi: 10.1056/NEJMoa1104621

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hodi FS, O'Day SJ, McDermott DF, Weber RW, Sosman JA, Haanen JB, et al. Improved survival with ipilimumab in patients with metastatic melanoma. N Engl J Med. (2010) 363:711–23. doi: 10.1056/NEJMoa1003466

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Schadendorf D, Hodi FS, Robert C, Weber JS, Margolin K, Hamid O, et al. Pooled analysis of long-term survival data from phase II and phase III trials of ipilimumab in unresectable or metastatic melanoma. J Clin Oncol. (2015) 33:1889–94. doi: 10.1200/JCO.2014.56.2736

CrossRef Full Text | Google Scholar

11. Johnson DB, Nixon MJ, Wang Y, Wang DY, Castellanos E, Estrada MV, et al. Tumor-specific MHC-II expression drives a unique pattern of resistance to immunotherapy via LAG-3/FCRL6 engagement. JCI Insight. (2018) 3:e120360. doi: 10.1172/jci.insight.120360

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of cancer. Immunity. (2018) 48:812−30 e14. doi: 10.1016/j.immuni.2018.03.023

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Romero-Cordoba S, Meneghini E, Sant M, Iorio MV, Sfondrini L, Paolini B, et al. Decoding immune heterogeneity of triple negative breast cancer and its association with systemic inflammation. Cancers. (2019) 11:911. doi: 10.3390/cancers11070911

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chen X, Shao Q, Hao S, Zhao Z, Wang Y, Guo X, et al. CTLA-4 positive breast cancer cells suppress dendritic cells maturation and function. Oncotarget. (2017) 8:13703–15. doi: 10.18632/oncotarget.14626

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Contardi E, Palmisano GL, Tazzari PL, Martelli AM, Fala F, Fabbi M, et al. CTLA-4 is constitutively expressed on tumor cells and can trigger apoptosis upon ligand interaction. Int J Cancer. (2005) 117:538–50. doi: 10.1002/ijc.21155

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Pistillo MP, Tazzari PL, Palmisano GL, Pierri I, Bolognesi A, Ferlito F, et al. CTLA-4 is not restricted to the lymphoid cell lineage and can function as a target molecule for apoptosis induction of leukemic cells. Blood. (2003) 101:202–9. doi: 10.1182/blood-2002-06-1668

CrossRef Full Text | Google Scholar

17. Zhang H, Dutta P, Liu J, Sabri N, Song Y, Li WX, et al. Tumour cell-intrinsic CTLA4 regulates PD-L1 expression in non-small cell lung cancer. J Cell Mol Med. (2019) 23:535–42. doi: 10.1111/jcmm.13956

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Kassardjian A, Shintaku PI, Moatamed NA. Expression of immune checkpoint regulators, cytotoxic T lymphocyte antigen 4 (CTLA-4) and programmed death-ligand 1 (PD-L1), in female breast carcinomas. PLoS ONE. (2018) 13:e0195958. doi: 10.1371/journal.pone.0195958

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Scheel AH, Dietel M, Heukamp LC, Johrens K, Kirchner T, Reu S, et al. Harmonized PD-L1 immunohistochemistry for pulmonary squamous-cell and adenocarcinomas. Mod Pathol. (2016) 29:1165–72. doi: 10.1038/modpathol.2016.117

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wang XB, Zheng CY, Giscombe R, Lefvert AK. Regulation of surface and intracellular expression of CTLA-4 on human peripheral T cells. Scand J Immunol. (2001) 54:453–8. doi: 10.1046/j.1365-3083.2001.00985.x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Fleischer J, Soeth E, Reiling N, Grage-Griebenow E, Flad HD, Ernst M. Differential expression and function of CD80 (B7-1) and CD86 (B7-2) on human peripheral blood monocytes. Immunology. (1996) 89:592–8. doi: 10.1046/j.1365-2567.1996.d01-785.x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Cruz UoCS. The UCSC Xena Platform for Public and Private Cancer Genomics Data Visualization and Interpretation. (2020). Available online at: https://xenabrowser.net/. (accessed January 8, 2020).

Google Scholar

23. Pereira B, Chin SF, Rueda OM, Vollan HK, Provenzano E, Bardwell HA, et al. The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nat Commun. (2016) 7:11479. doi: 10.1038/ncomms11908

CrossRef Full Text | Google Scholar

24. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. (2012) 2:401–4. doi: 10.1158/2159-8290.CD-12-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. (2010) 26:2363–7. doi: 10.1093/bioinformatics/btq431

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc. (2009) 4:1184–91. doi: 10.1038/nprot.2009.97

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. (2013) 4:2612. doi: 10.1038/ncomms3612

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. (2009) 462:108−12. doi: 10.1038/nature08460

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. (2013) 14:7. doi: 10.1186/1471-2105-14-7

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. (2017) 18:248–62. doi: 10.1016/j.celrep.2016.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Lehmann BD, Bauer JA, Chen X, Sanders ME, Chakravarthy AB, Shyr Y, et al. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest. (2011) 121:2750–67. doi: 10.1172/JCI45014

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, et al. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res. (2010) 12:R68. doi: 10.1186/bcr2635

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. (2015) 160:48–61. doi: 10.1016/j.cell.2014.12.033

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ayers M, Lunceford J, Nebozhyn M, Murphy E, Loboda A, Kaufman DR, et al. IFN-gamma-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Invest. (2017) 127:2930–40. doi: 10.1172/JCI91190

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Savas P, Virassamy B, Ye C, Salim A, Mintoff CP, Caramia F, et al. Single-cell profiling of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis. Nat Med. (2018) 24:986–93. doi: 10.1038/s41591-018-0078-7

CrossRef Full Text | Google Scholar

39. Ock CY, Hwang JE, Keam B, Kim SB, Shim JJ, Jang HJ, et al. Genomic landscape associated with potential response to anti-CTLA-4 treatment in cancers. Nat Commun. (2017) 8:1050. doi: 10.1038/s41467-017-01018-0

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Acs B, Madaras L, Tokes AM, Kovacs AK, Kovacs E, Ozsvari-Vidakovich M, et al. PD-1, PD-L1 and CTLA-4 in pregnancy-related - and in early-onset breast cancer: a comparative study. Breast. (2017) 35:69–77. doi: 10.1016/j.breast.2017.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Maszyna F, Hoff H, Kunkel D, Radbruch A, Brunner-Weinzierl MC. Diversity of clonal T cell proliferation is mediated by differential expression of CD152 (CTLA-4) on the cell surface of activated individual T lymphocytes. J Immunol. (2003) 171:3459–66. doi: 10.4049/jimmunol.171.7.3459

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Walunas TL, Lenschow DJ, Bakker CY, Linsley PS, Freeman GJ, Green JM, et al. CTLA-4 can function as a negative regulator of T cell activation. Immunity. (1994) 1:405–13. doi: 10.1016/1074-7613(94)90071-X

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Krummel MF, Allison JP. CD28 and CTLA-4 have opposing effects on the response of T cells to stimulation. J Exp Med. (1995) 182:459–65. doi: 10.1084/jem.182.2.459

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Chambers CA, Kuhns MS, Egen JG, Allison JP. CTLA-4-mediated inhibition in regulation of T cell responses: mechanisms and manipulation in tumor immunotherapy. Annu Rev Immunol. (2001) 19:565–94. doi: 10.1146/annurev.immunol.19.1.565

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Parry RV, Chemnitz JM, Frauwirth KA, Lanfranco AR, Braunstein I, Kobayashi SV, et al. CTLA-4 and PD-1 receptors inhibit T-cell activation by distinct mechanisms. Mol Cell Biol. (2005) 25:9543–53. doi: 10.1128/MCB.25.21.9543-9553.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Olsson C, Riesbeck K, Dohlsten M, Michaelsson E. CTLA-4 ligation suppresses CD28-induced NF-kappaB and AP-1 activity in mouse T cell blasts. J Biol Chem. (1999) 274:14400–5. doi: 10.1074/jbc.274.20.14400

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Kubsch S, Graulich E, Knop J, Steinbrink K. Suppressor activity of anergic T cells induced by IL-10-treated human dendritic cells: association with IL-2- and CTLA-4-dependent G1 arrest of the cell cycle regulated by p27Kip1. Eur J Immunol. (2003) 33:1988–97. doi: 10.1002/eji.200323600

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Calvo CR, Amsen D, Kruisbeek AM. Cytotoxic T lymphocyte antigen 4 (CTLA-4) interferes with extracellular signal-regulated kinase (ERK) and Jun NH2-terminal kinase (JNK) activation, but does not affect phosphorylation of T cell receptor zeta and ZAP70. J Exp Med. (1997) 186:1645–53. doi: 10.1084/jem.186.10.1645

CrossRef Full Text | Google Scholar

49. Brunner MC, Chambers CA, Chan FK, Hanke J, Winoto A, Allison JP. CTLA-4-Mediated inhibition of early events of T cell proliferation. J Immunol. (1999) 162:5813–20.

PubMed Abstract | Google Scholar

50. Blair PJ, Riley JL, Levine BL, Lee KP, Craighead N, Francomano T, et al. CTLA-4 ligation delivers a unique signal to resting human CD4 T cells that inhibits interleukin-2 secretion but allows Bcl-X(L) induction. J Immunol. (1998) 160:12–5.

PubMed Abstract | Google Scholar

51. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. (2017) 18:220. doi: 10.1186/s13059-017-1349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Nelson BH. IL-2, regulatory T cells, and tolerance. J Immunol. (2004) 172:3983–8. doi: 10.4049/jimmunol.172.7.3983

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Ahmadzadeh M, Rosenberg SA. IL-2 administration increases CD4+ CD25(hi) Foxp3+ regulatory T cells in cancer patients. Blood. (2006) 107:2409–14. doi: 10.1182/blood-2005-06-2399

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Hori S, Nomura T, Sakaguchi S. Control of regulatory T cell development by the transcription factor Foxp3. Science. (2003) 299:1057–61. doi: 10.1126/science.1079490

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Conlon KC, Lugli E, Welles HC, Rosenberg SA, Fojo AT, Morris JC, et al. Redistribution, hyperproliferation, activation of natural killer cells and CD8 T cells, and cytokine production during first-in-human clinical trial of recombinant human interleukin-15 in patients with cancer. J Clin Oncol. (2015) 33:74–82. doi: 10.1200/JCO.2014.57.3329

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, et al. A subcellular map of the human proteome. Science. (2017) 356:eaal3321. doi: 10.1126/science.aal3321

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Wang E, Monaco A, Monsurro V, Sabatino M, Pos Z, Uccellini L, et al. Antitumor vaccines, immunotherapy and the immunological constant of rejection. IDrugs. (2009) 12:297–301.

PubMed Abstract | Google Scholar

58. Leach DR, Krummel MF, Allison JP. Enhancement of antitumor immunity by CTLA-4 blockade. Science. (1996) 271:1734–6. doi: 10.1126/science.271.5256.1734

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Dong H, Strome SE, Salomao DR, Tamura H, Hirano F, Flies DB, et al. Tumor-associated B7-H1 promotes T-cell apoptosis: a potential mechanism of immune evasion. Nat Med. (2002) 8:793–800. doi: 10.1038/nm730

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Mao H, Zhang L, Yang Y, Zuo W, Bi Y, Gao W, et al. New insights of CTLA-4 into its biological function in breast cancer. Curr Cancer Drug Targets. (2010) 10:728–36. doi: 10.2174/156800910793605811

CrossRef Full Text | Google Scholar

61. Yu H, Yang J, Jiao S, Li Y, Zhang W, Wang J. Cytotoxic T lymphocyte antigen 4 expression in human breast cancer: implications for prognosis. Cancer Immunol Immunother. (2015) 64:853–60. doi: 10.1007/s00262-015-1696-2

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Slavik JM, Hutchcroft JE, Bierer BE. CD28/CTLA-4 and CD80/CD86 families: signaling and function. Immunol Res. (1999) 19:1–24. doi: 10.1007/BF02786473

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Scheid MP, Woodgett JR. Unravelling the activation mechanisms of protein kinase B/Akt. FEBS Lett. (2003) 546:108–12. doi: 10.1016/S0014-5793(03)00562-3

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Hanada M, Feng J, Hemmings BA. Structure, regulation and function of PKB/AKT–a major therapeutic target. Biochim Biophys Acta. (2004) 1697:3–16. doi: 10.1016/j.bbapap.2003.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Baroja ML, Vijayakrishnan L, Bettelli E, Darlington PJ, Chau TA, Ling V, et al. Inhibition of CTLA-4 function by the regulatory subunit of serine/threonine phosphatase 2A. J Immunol. (2002) 168:5070–8. doi: 10.4049/jimmunol.168.10.5070

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Gallay N, Dos Santos C, Cuzin L, Bousquet M, Simmonet Gouy V, Chaussade C, et al. The level of AKT phosphorylation on threonine 308 but not on serine 473 is associated with high-risk cytogenetics and predicts poor overall survival in acute myeloid leukaemia. Leukemia. (2009) 23:1029–38. doi: 10.1038/leu.2008.395

CrossRef Full Text | Google Scholar

67. Mo X, Zhang H, Preston S, Martin K, Zhou B, Vadalia N, et al. Interferon-gamma signaling in melanocytes and melanoma cells regulates expression of CTLA-4. Cancer Res. (2018) 78:436–50. doi: 10.1158/0008-5472.CAN-17-1615

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Tian L, Xie XH, Zhu ZH. Calotropin regulates the apoptosis of nonsmall cell cancer by regulating the cytotoxic Tlymphocyte associated antigen 4mediated TGFbeta/ERK signaling pathway. Mol Med Rep. (2018) 17:7683–91. doi: 10.3892/mmr.2018.8853

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Tsatsanis C, Vaporidi K, Zacharioudaki V, Androulidaki A, Sykulev Y, Margioris AN, et al. Tpl2 and ERK transduce antiproliferative T cell receptor signals and inhibit transformation of chronically stimulated T cells. Proc Natl Acad Sci USA. (2008) 105:2987–92. doi: 10.1073/pnas.0708381104

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Wang B, Qin L, Ren M, Sun H. Effects of combination of anti-CTLA-4 and anti-PD-1 on gastric cancer cells proliferation, apoptosis and metastasis. Cell Physiol Biochem. (2018) 49:260–70. doi: 10.1159/000492876

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Nishimura R, Arima N. Is triple negative a prognostic factor in breast cancer? Breast Cancer. (2008) 15:303–8. doi: 10.1007/s12282-008-0042-3

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Rosskopf S, Leitner J, Zlabinger GJ, Steinberger P. CTLA-4 antibody ipilimumab negatively affects CD4(+) T-cell responses in vitro. Cancer Immunol Immunother. (2019) 68:1359–68. doi: 10.1007/s00262-019-02369-x

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Reichert TE, Watkins S, Stanson J, Johnson JT, Whiteside TL. Endogenous IL-2 in cancer cells: a marker of cellular proliferation. J Histochem Cytochem. (1998) 46:603–11. doi: 10.1177/002215549804600506

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Rangel-Corona R, Corona-Ortega T, Soto-Cruz I, Lopez-Labra A, Pablo-Arcos T, Torres-Guarneros CF, et al. Evidence that cervical cancer cells secrete IL-2, which becomes an autocrine growth factor. Cytokine. (2010) 50:273–7. doi: 10.1016/j.cyto.2010.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Mindiola R, Caulejas D, Nunez-Troconis J, Araujo M, Delgado M, Mosquera J. Increased number of IL-2, IL-2 receptor and IL-10 positive cells in premalignant lesions of the cervix. Invest Clin. (2008) 49:533–45.

PubMed Abstract | Google Scholar

76. Kasprzak A, Olejniczak K, Przybyszewska W, Zabel M. Cellular expression of interleukin 2 (IL-2) and its receptor (IL-2R, CD25) in lung tumours. Folia Morphol. (2007) 66:159–66.

PubMed Abstract | Google Scholar

77. Garcia-Tunon I, Ricote M, Ruiz A, Fraile B, Paniagua R, Royuela M. Interleukin-2 and its receptor complex (alpha, beta and gamma chains) in in situ and infiltrative human breast cancer: an immunohistochemical comparative study. Breast Cancer Res. (2004) 6:R1–7. doi: 10.1186/bcr730

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Kane LP, Andres PG, Howland KC, Abbas AK, Weiss A. Akt provides the CD28 costimulatory signal for up-regulation of IL-2 and IFN-gamma but not TH2 cytokines. Nat Immunol. (2001) 2:37–44. doi: 10.1038/83144

CrossRef Full Text | Google Scholar

79. Tischner D, Wiegers GJ, Fiegl H, Drach M, Villunger A. Mutual antagonism of TGF-beta and Interleukin-2 in cell survival and lineage commitment of induced regulatory T cells. Cell Death Differ. (2012) 19:1277–87. doi: 10.1038/cdd.2012.7

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Das L, Levine AD. TGF-beta inhibits IL-2 production and promotes cell cycle arrest in TCR-activated effector/memory T cells in the presence of sustained TCR signal transduction. J Immunol. (2008) 180:1490–8. doi: 10.4049/jimmunol.180.3.1490

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Wang XB, Giscombe R, Yan Z, Heiden T, Xu D, Lefvert AK. Expression of CTLA-4 by human monocytes. Scand J Immunol. (2002) 55:53–60. doi: 10.1046/j.0300-9475.2001.01019.x

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Pech MF, Fong LE, Villalta JE, Chan LJ, Kharbanda S, O'Brien JJ, et al. Systematic identification of cancer cell vulnerabilities to natural killer cell-mediated immune surveillance. Elife. (2019) 8:e47362. doi: 10.7554/eLife.47362

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Lindsten T, Lee KP, Harris ES, Petryniak B, Craighead N, Reynolds PJ, et al. Characterization of CTLA-4 structure and expression on human T cells. J Immunol. (1993) 151:3489–99.

PubMed Abstract | Google Scholar

84. Freeman GJ, Borriello F, Hodes RJ, Reiser H, Hathcock KS, Laszlo G, et al. Uncovering of functional alternative CTLA-4 counter-receptor in B7-deficient mice. Science. (1993) 262:907–9. doi: 10.1126/science.7694362

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Tumeh PC, Harview CL, Yearley JH, Shintaku IP, Taylor EJ, Robert L, et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature. (2014) 515:568–71. doi: 10.1038/nature13954

PubMed Abstract | CrossRef Full Text | Google Scholar

86. Snyder A, Makarov V, Merghoub T, Yuan J, Zaretsky JM, Desrichard A, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. (2014) 371:2189–99. doi: 10.1056/NEJMoa1406498

PubMed Abstract | CrossRef Full Text | Google Scholar

87. Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. (2018) 554:544–8. doi: 10.1038/nature25501

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: triple-negative breast cancer, CTLA-4, immunotherapy, immune landscape, cancer genomics

Citation: Navarrete-Bernal MGC, Cervantes-Badillo MG, Martínez-Herrera JF, Lara-Torres CO, Gerson-Cwilich R, Zentella-Dehesa A, Ibarra-Sánchez MJ, Esparza-López J, Montesinos JJ, Cortés-Morales VA, Osorio-Pérez D, Villegas-Osorno DA, Reyes-Sánchez E, Salazar-Sojo P, Tallabs-Utrilla LF and Romero-Córdoba S and Rocha-Zavaleta L (2020) Biological Landscape of Triple Negative Breast Cancers Expressing CTLA-4. Front. Oncol. 10:1206. doi: 10.3389/fonc.2020.01206

Received: 15 February 2020; Accepted: 15 June 2020;
Published: 05 August 2020.

Edited by:

Laura Cerchia, Italian National Research Council, Italy

Reviewed by:

Martin Hernán Bonamino, National Cancer Institute (INCA), Brazil
Stefania Scala, Istituto Nazionale Tumori IRCCS “Fondazione G. Pascale”, Italy

Copyright © 2020 Navarrete-Bernal, Cervantes-Badillo, Martínez-Herrera, Lara-Torres, Gerson-Cwilich, Zentella-Dehesa, Ibarra-Sánchez, Esparza-López, Montesinos, Cortés-Morales, Osorio-Pérez, Villegas-Osorno, Reyes-Sánchez, Salazar-Sojo, Tallabs-Utrilla, Romero-Córdoba and Rocha-Zavaleta. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sandra Romero-Córdoba, sromero_cordoba@hotmail.co; Leticia Rocha-Zavaleta, lrochaz@iibiomedicas.unam.mx

These authors have contributed equally to this work