Next Article in Journal
Comparison of the Transcriptomic Signatures in Pediatric and Adult CML
Next Article in Special Issue
Targeted Deletion of Kindlin-2 in Mouse Mammary Glands Inhibits Tumor Growth, Invasion, and Metastasis Downstream of a TGF-β/EGF Oncogenic Signaling Pathway
Previous Article in Journal
The Estrogen Receptor α Signaling Pathway Controls Alternative Splicing in the Absence of Ligands in Breast Cancer Cells
Previous Article in Special Issue
Anthocyanidins Inhibit Growth and Chemosensitize Triple-Negative Breast Cancer via the NF-κB Signaling Pathway
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 24
10.3390/cancers13246262
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

YB1 Is a Major Contributor to Health Disparities in Triple Negative Breast Cancer

by
Priyanka Shailendra Rana
1,2,†,
Wei Wang
1,2,†,
Akram Alkrekshi
1,2,
Vesna Markovic
2,
Amer Khiyami
2,
Ricky Chan
3,
Adam Perzynski
1,2,
Natalie Joseph
1,2 and
Khalid Sossey-Alaoui
1,2,3,*
1
Department of Medicine, MetroHealth Medical Center, Cleveland, OH 44109, USA
2
Case Western Reserve University School of Medicine, Case Western Reserve University, Cleveland, OH 44016, USA
3
Case Comprehensive Cancer Center, Case Western Reserve University, Cleveland, OH 44106, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2021, 13(24), 6262; https://doi.org/10.3390/cancers13246262
Submission received: 11 November 2021 / Revised: 6 December 2021 / Accepted: 9 December 2021 / Published: 14 December 2021
(This article belongs to the Special Issue Breast Development and Cancer)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Triple negative breast cancer (TNBC) is a devastating disease that affects many women, due to the lack of FDA-approved targeted therapy. In the absence of cell surface receptors ER, PR, and Her2 that can be targeted with hormonal and antibody treatments, cytotoxic chemotherapy remains the major course of treatment, with a dismal response and rapid recurrence due to the acquisition of resistance. TNBC is also twice as more prevalent in African American (AA) when compared to Caucasian American (CA) women. This study investigated the role of the YB1 gene in the disparities in TNBC between AA and CA women. We found that YB1 is highly expressed in TNBC tumors of AA origin when compared to CAs. Increased expression levels and activity of YB1 correlates with poor disease outcomes, resistance to chemotherapy, and the activation of the cancer stem cell (CSC) phenotype, with higher levels in AA than in CA TNBC tumors. More importantly, we found that the targeted inhibition of the expression and activity of YB1 significantly inhibited the oncogenic behavior of AA tumors through sensitization to chemotherapy and inhibition of CSCs. Our study is the first to show that YB1 activity may be a major biological contributor to the health disparities in TNBC, and that development of therapies that specifically target YB1 could reduce these disparities.

Abstract

Triple negative breast cancer (TNBC) is the most aggressive amongst all breast cancer (BC) subtypes. While TNBC tumors represent less than 20% of all BC subtypes, they are responsible for the most BC-related deaths. More significantly, when considering TNBC incidence across all racial/ethnic groups, TNBC accounts for less than 20% of all BCs. However, in non-Hispanic black women, the incidence rate of TNBC is more than 40%, which may be a contributing factor to the higher BC-related death rate in this population. These disparities remain strong even after accounting for differences in socioeconomic status, healthcare access, and lifestyle factors. Increased evidence now points to biological mechanisms that are intrinsic to the tumor that contribute to disparate TNBC disease burdens. Here, we show that YB1, a multifunction gene, plays a major role in the TNBC disparities between African American (AA) and Caucasian American (CA) women. We show in three independent TNBC tumors cohorts, that YB1 is significantly highly expressed in AA TNBC tumors when compared to CAs, and that increased levels of YB1 correlate with poor survival of AA patients with TNBC. We used a combination of genetic manipulation of YB1 and chemotherapy treatment, both in vitro and in animal models of TNBC to show that YB1 oncogenic activity is more enhanced in TNBC cell lines of AA origin, by increasing their tumorigenic and aggressive behaviors, trough the activation of cancer stem cell phenotype and resistance to chemotherapeutic treatments.

1. Introduction

In women, ~90% of breast cancer-related deaths are caused by metastasis [1]. In 2021, breast cancer (BC) surpassed all other cancers affecting women, and ranks number two as the cause of cancer-related deaths in women [2]. BC is a heterogeneous disease that can be divided into at least five genetically distinct subtypes, where more than 10 subtypes can be identified based on their molecular signatures [2,3,4,5,6,7]. Triple-negative BC (TNBC), a BC subtype that lacks the expression of hormone receptors (ER-a and PR) and ErbB2/HER2 [8,9,10,11,12], is notoriously known for its aggressiveness and lethality because of its high metastatic and recurrence potentials, coupled with the lack of FDA-approved targeted therapies that could be effective against this BC subtype. Additionally, when TNBC tumors recur, they are generally associated with the acquisition of resistance to standard-of-care chemotherapies through mechanisms that remain to be resolved. Interestingly, while across all racial/ethnic groups, TNBC accounts for less than 20%, its incidence rate more than doubles in non-Hispanic black women, which may contribute to the higher breast cancer death rate in this population. There is now increasing evidence that African American (AA) women with TNBC have worse clinical outcomes when compared with Caucasian American (CA) women who have the same disease. At MetroHealth Medical Center, the site of this study and the largest SafetyNet hospital in north-east Ohio that is also serving the most disadvantaged populations in the Cleveland area, 39% of our patients treated for BC are AAs when compared to the national average of 11%. In fact, of the 113 women with TNBC that were treated at MetroHealth Medical Center between 2014 and 2018, 58 (51%) were AAs. These statistics, that are characteristic of our sociodemographic area, place us in a unique position to study cancer health disparities in our TNBC patient population. More importantly, the disparities in TNBC still remain poorly understood, in part due to the complex interactions between genetics, social environment, and lifestyle that may contribute to the observed differences. Moreover, there is growing evidence that the observed disparities in TNBC outcomes between AAs and CAs are not exclusively the result of socioeconomic or healthcare access factors, but rather reflect biological mechanisms intrinsic to the tumor that readily contribute to unequal TNBC disease burdens [13,14,15].
The Y Box binding protein 1 (YB1) is a multifunctional protein, found both in the cytoplasm and inside the nucleus, that belongs to the highly conserved Cold Shock Domain protein family [16]. Among its multitude of functions, YB1 is more known to be involved in the regulation of transcription, mRNA stability, and splicing [17,18,19,20,21,22,23,24,25,26,27,28]. More recently, YB1 was found to play key roles in the regulation of tumor development, progression and metastasis, as well as drug resistance in several cancers, including the ones generated in the breast [17,24,26,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45]. In this study, we show that dysregulated YB1 expression to be associated primarily with breast cancer cell lines in primary tumors of the TNBC subtype as compared to their luminal or HER2+ counterparts. We also show that YB1 abundance correlates with the metastatic potential of TNBC cell lines, and robust YB1 expression correlates with reduced overall survival in BC patients. We also show that YB1 expression levels to be significantly (p < 0.01) higher in TNBC cell lines and tumors of AA origin when compared to their CA counterparts. Importantly, this association remained significant even after controlling for known prognostic factors, such as tumor stage, age, and node status, suggesting that YB1 activation is likely to be a population-specific prognostic factor in TNBC. Mechanistically, we show that aggressiveness of TNBC tumors of AA origin is driven by the YB1-mediated activation of resistance to chemotherapeutic drugs, by activating genes that are associated with the cancer stem cell phenotype and chemoresistance. Our data identify YB1 as a major biological factor of cancer health disparities in AA women with TNBC tumors, and open opportunities for the development of targeted therapies against TNBC tumors that are disparately affecting AA women.

2. Materials and Methods

2.1. Ethics Statement

All animal studies were performed under protocols approved by the Institutional Animal Care and Use Committee and conducted in accordance with the guidelines and regulations set and approved by the MetroHealth Medical Center, Case Western Reserve University, and NIH. For this study, we used six to eight week-old female NSG or Balb/C mice (Jackson Laboratory, Farmington, CT, USA).

2.2. Cell Lines and Reagents

TNBC cell lines 4T1, MDA-MB-231, MDA-MB-468, HCC1806, HCC38, HCC70, and BT549 were obtained from American Type Culture Collection (ATCC) and maintained according to the manufacturer’s protocols. Cell lines were also routinely authenticated by STR DNA fingerprinting analysis. The Yb1-KO cells were generated by electroporation as described below.

2.3. Antibodies and Reagents

We used the following antibodies: rabbit antibodies against pAKT-S473 (1:1000), pAKT-T308 (1:1000), AKT (1:1000), p-Yb1 (Cell Signaling Technology, Danvers, MA, USA) (1:1000); Yb-1 (Abcam, Cambridge UK) (1:1000); goat horseradish-peroxidase-conjugated anti-mouse IgG and goat horseradish-peroxidase-conjugated anti-rabbit IgG (Biorad, Hercules, CA, USA) (1:2000) and mouse monoclonal anti-actin (Sigma, St. Louis, MO, USA) (1:5000). ECL reagent was from Thermo Scientific (Hanover Park, IL, USA). The antibodies were dissolved to a working concentration either in 5% BSA (Primary antibodies) or 5% non-fat dry milk (secondary antibodies). Cisplatin was obtained from Sigma and doxorubicin was obtained from European Pharmacopoeia, and used at 30 μM and 200 nM, respectively, to treat cell cultures for 48 h.

2.4. Immunoblotting

Immunoblotting analyses were performed as described previously [46,47]. Cells cultured on 6-well plates were washed with ice-cold PBS twice before lysis. Ice-cold RIPA lysis buffer with protease and phosphatase inhibitors was used to lyse the adherent cells and the lysate was transferred to prechilled 1.5 mL microcentrifuge tubes. The samples were then centrifuged at 15,000 RPM for 20 min at 4 °C. Cell lysates with an equal amount of total protein (25 mg) were first digested in SDS sample buffer, then resolved on a 10% SDS-polyacrylamide gel, followed by transferring the proteins to PVDF membrane. The membrane was incubated in a blocking solution (5% BSA) to prevent non-specific binding at room temperature for 60 min. After blocking, the membrane was then incubated with primary antibody overnight at 4 °C followed by washing and incubating again for an hour with an appropriate secondary antibody at room temperature. The membrane was washed, and signals were developed with a Pierce ECL Western Blot chemiluminescence detection kit and further visualized and imaged using ChemiDoc MP Imaging System. Signals were quantified using ImageJ software, according to the parameters described in ImageJ user guide (http://rsbweb.nih.gov/ij/docs/guige/146.html, accessed on 2 July 2021. Average values from 3 different blots are presented.

2.5. RNA Extraction and Quantitative Real Time Reverse Transcription PCR (qt-RT-PCR)

RNA was extracted from cell lines as described previously, according to the manufacturer’s instructions [48]. Cells were lysed in a 6 well plate using TRIzol reagent (Invitrogen), and total RNA resuspended in RNase-free water was quantified using a Nanodrop 2000 spectrophotometer (ThermoFisher Scientific, Waltham, MA, USA). To generate cDNA for qtRT -PCR, 1 mg RNA was used with the SuperScript III First-Strand Synthesis System RT-PCR kit (Invitrogen, Waltham, MA, USA). qtRT -PCR was done as described previously, using the C1000 Touch Thermal Cycler CFX96 Real-Time System (Bio-Rad, Hercules, CA, USA). Oligonucleotide primers used for qtRT-PCR were described in [49].

2.6. Colony Formation Assay

Clonogenic growth was assessed by seeding 5000 cells (for MDA-MB-231, MDA-MB-157, HCC70, HCC38) or 10,000 cells (for MDA-MB-468 and BT549) onto a 6 well plate and incubating them for 5 days until macroscopically visible colonies were formed. Growth media were replaced once every other day. Colonies were visualized by fixing cells in 4% PFA for 20 min, were washed with PBS, and then stained with 0.05% crystal violet [50]. Plates were washed with water, and colonies were counted under a dissecting scope.

2.7. Generation of Cisplatin- and Doxorubicin-Resistant TNBC Cells

MDA-MB-231 or MDA-MB-468 cells were cultured with increasing concentrations of either cisplatin or doxorubicin over a 60-day period, starting with 1 µM for 7 days, 10 µM for 14 days, 21 µM, 30 µM (MDA-MB-231), or 40 µM (MDA-MB-468) for 30 days, for cisplatin. For doxorubicin, cells were cultured with 50 nM for 7 days, 100 nM for 14 days, 150 nM for 21 days, followed by 300 nM (MDA-MB-231) or 450 nM (MDA-MB-468) for 30 days. After each cycle, the IC50 was determined by MTT to determine the concentration for the next cycle. After the 60 day period, resistance to cisplatin or doxorubicin was established and the derived cells were used for subsequent analyses.

2.8. sgRNA Preparation and Electroporation for Generating eYb1KO Cell Lines

The sgRNA pool of 3 guide RNAs for mouse Yb1 (A*C*G*GGCAGCGGCGCGGGUAG, G*G*C*GGGGACAAGAAGGUCAU, and G*G*C*CCGAGCCACGGACUACG) and 2 guide RNAs for human Yb1 (U*U*U*UCCAGCAACGAAGGUUU and U*U*C*AUCAACAGGUGAGCUGC), with their respective Synthego modified EZ scaffolds targeting Yb1, were obtained from Synthego. The sgRNA pellets were rehydrated in 1X TE buffer (provided by Synthego, Menlo Park, CA, USA) to make a stock of 100 µM. A working solution of 30 µM sgRNA was made (in nuclease-free water) fresh before electroporating the cells. For every reaction, the RNP (Ribonucleoprotein) complex was assembled by adding 3 µl of 30 µM sgRNA to 0.5 µl of 20 µM Cas9 (provided by Synthego) at a ratio of 9:1 in 3.5 µl resuspension buffer R, provided by Neon Transfection System; Invitrogen Basel, Switzerland, and was incubated for 10 min at room temperature.

2.9. Generation of Yb-1 Knockout Cell Lines Using Electroporation

Electroporation of all 7 TNBC cell lines (4T1, MDA-MB-231, MDA-MB-468, HCC1806, HCC38, HCC70 and BT549) was achieved by an implemented electroporation device system, according to manufacturer’s instructions (Neon Transfection System; Invitrogen, Basel, Switzerland). The Neon Transfection System 10 µL kit was used for the transfection of human and mouse TNBC cells. Cells were subcultured 48 h before electroporation and harvested at nearly 80% confluency. Cells at a density of 2 × 105 were washed with PBS and resuspended in 5 µl Resuspension buffer R (Neon Transfection System; Invitrogen Basel, Switzerland). Within 15 min of resuspension, the cells were added to the tube containing RNP and the cell-RNP complex was electroporated with the Neon Transfection System. Per electroporation, 2 × 105 cells were taken up in a 10 µL Neon tip using the Neon Transfection System pipette (Invitrogen). The electroporation was performed by applying 4 pulses at 1400 Volts for 10 ms to all the TNBC cell lines. The respective control cells were incubated with the resuspension buffer without the sgRNA and electroporated at the same settings. After electroporation, the cells were seeded in 6 well plate by adding 2.5 mL DMEM (Cytiva) or RPMI (ATCC) with 10% FBS without antibiotic supplements. The cells were cultured for 48–72 h and subsequently proceeded for further analysis. We used Western Blot to verify the loss of YB1 expression.

2.10. Animal Experiments

For MDA-MB-231 and MDA-MB-468 human TNBC cell lines, parental (control) or Yb1-KO cells (1 × 106 cells/injection) were injected into mammary fat pads in both sites (n = 10 tumors) of NSG mice (n = 5 mice). For the 4T1 murine TNBC cell lines, 100 K cells per injection were used to inject syngeneic BalbC mice. Mice were anaesthetized using isofluorane prior to injecting tumor cells. Tumor growth was followed by twice weekly monitoring of tumor volume with digital Vernier calipers. At the start of the experiment, mice were euthanized, and tumors were excised, weighed, and snap-frozen in liquid N2 and stored in a −80 °C freezer for subsequent analyses. Lung metastasis nodules were counted, and the results were plotted as the average number of foci per lobe. For the animal treatment with cisplatin, 21 days post injection of parental or cis-resistant MDA-MB-468 cells, when the first tumors grew to ~100 mm3, mice were then treated with cisplatin (1 mg/kg, IP) 2 times a week for 4 weeks. When tumor volume had reached 1000 mm3, all animals were euthanized.

2.11. Human TNBC Specimens

A total of 45 TNBC tumors (25 AA and 20 CA, self-reported ethnicity) were identified from the MetroHealth breast cancer tumor biobank of FFPE tumors. All FFPE TNBC tumor blocks were confirmed for tumor content and TNBC diagnosis by a subspecialty breast pathologist. Total RNA was isolated from the tumor sections using the Qiagen AllPrep FFPE DNA/RNA kit in accordance with manufacturer’s protocols. After RNA extraction, RNA purity was assessed using the 260/280 ratio, as well as the RNA integrity number (RIN). Only those RNA samples with an RIN > 7 were retained for subsequent qt-RT-PCR as described above.

2.12. Analyses of RNA-Seq Data from the TCGA Datasets

Data generated from the illumina sequencing platform were qualitatively assessed and trimmed to remove the adapter sequences, with the help of TrimGalore! v0.4.2 (Babraham Bioinformatics, Cambrigde, UK), a wrapper script for FastQC, and cutadapt. Sequence reads that passed quality control were aligned with the human reference genome (GRCh38) using the STAR aligner v2.5.1 [51]. The sequence alignment in this case was guided with the help of GENCODE annotation for hg38. The aligned reads were then analyzed for differential expression using an RNASeq analysis package—Cufflinks v2.2.1 [52], which reports the fragments per kb of exon per million fragments mapped (FPKM) for each gene. Cuffdiff was used to generate the differential analysis report. A significance cutoff of q-value < 0.05 (after Benjamini-Hochberg correction for multiple testing) was used to identify differential genes, which were then subjected to gene set enrichment and pathway analysis using iPathwayGuide (AdvaitaBio, Ann Arbor, MI, USA) to further determine any relevant biological processes that may be differentially represented between the control and the experimental groups.

2.13. Bioinformatic Methods for Datamining of TCGA PanCancer Dataset

BC patient data from TCGA PanCancer Atlas [53] was accessed using cBioPortal [54,55]. The PanCancer Atlas includes 1082 patients with mRNA expression data and 151 individuals with data available for survival analysis. Race information for the samples was obtained through cBioPortal and 182 patients identified as African American and 749 identified as Caucasian. The distribution of BC subtypes by race is provided in Table 1. Gene expression differences between normal and tumor pairs were tested for significance by the Wilcoxon Signed Rank Sum test. Differences between subtypes were tested for significance using the Wilcoxon Rank Sum test. The mRNA co-expressions were reported and graphed for the genes of interest based on RSEM values which were batch-normalized from Illumina HiSeq RNASeq data. The protein expression levels were reported and graphed for the genes of interest based on protein expression (RPPA). Additionally, the Spearman correlation between genes of interest were reported, along with their p-value and multiple testing corrected q-value, obtained using cBioPortal. Overall survival graphs for thegenes of interest were generated using R (‘survival v3.2-3′ package from CRAN) using the survival data obtained from cBioPortal.

2.14. Other Statistical Analyses

Experiments were done in triplicate and analyzed using the Student’s t-test. In calculating the two-tailed significance levels for equality of means, equal variances were assumed for the two populations. Results were considered significant at p < 0.05.

3. Results

3.1. YB1 Expression Levels Are Increased in TNBC Cell Lines and Tumors and Are Associated with Poor Survival and Worst Patients’ Outcomes

We previously reported that expression levels of YB1 were significantly elevated in aggressive BC cancer cell lines when compared to their less aggressive counterparts [31]. More importantly, we showed that YB1 expression levels are significantly higher in TNBC cell lines when compared to all other BC subtypes [31], which indicated an association between YB1 expression and the aggressiveness of TNBC cell lines. We also showed that YB1 expression correlated with the metastatic potential of three TNBC metastasis progression models, namely the human MCA10A series [56], the murine 4T1 series [57], and the NME series [58,59]. To validate these findings in human patients’ specimens, we interrogated the datasets of two large, independent BC cohorts. The first cohort, the TCGA PanCancer BC cohort, which contains clinical information on 1082 BC patients, including 181 TNBCs (Table 1), showed a significant (p < 0.0001) positive association between YB1 expression levels and the basal (TNBC) BC subtype (Figure 1A). Elevated expression levels of YB1 in the TNBC subtype also correlated (p = 0.03) with poor disease outcome and reduced survival probability in this cohort (Figure 1B). Next, we interrogated the BC KM-Plotter (https://kmplot.com/analysis, accessed on 6 January 2021) cohort, which contains clinical information on ~5000 BC patients, including 405 TNBCs, and this showed a very significant (p = 1.9 × 108), positive correlation between YB1 mRNA expression and reduced survival probability (Figure 1C, upper panel). In this cohort, BC patients with high levels of YB1 in their tumors had worst clinical outcomes when compared to patients with low levels of YB1, and patients with high tumor YB1 have an average reduced survival of 31 months when compared to patients with low tumor YB1 (Figure 1C, lower panel). This correlation between YB1 expression levels and poor prognosis was also confirmed at the protein expression level (Figure 1D). Finally, we found YB1 mRNA expression levels in the TNBC cohort to also correlate with reduced survival (Figure 1E). Thus, the data support the function of YB1 as a promoter of BC aggressiveness in patients with TNBC tumors.

3.2. YB1 Is Disparately Highly Expressed in AA vs. CA TNBC Cell Lines and Tumors and Correlates with Poor Outcome in AA TNBC

The literature provides overwhelming evidence that incidence of TNBC and TNBC-related death are disparately elevated in AA women when compared to their CA counterparts. The socioeconomic status and lack of access to health care by the AA population cannot alone explain these disparities, which suggests a difference in the biology of TNBC tumors between AA and CA might also contribute to the observed disparities. We therefore set out to investigate whether the oncogenic activity of YB1 in TNBC tumors contributes to these disparities. First, we assessed the expression levels of YB1 by RT-PCR in six established TNBC cell lines; three from AA and three from CA patients, and found YB1 to be significantly (p < 0.001) higher in AA cell lines MDA-MB-468, HCC70, and HCC180 when, compared to CA TNBC cell lines MDA-MB-231, BT549, and HCC38 (Figure 2A,B). We confirmed the results that were obtained in our cell lines in the human specimens from the TCGA PanCancer cohort, and found YB1 expression levels to be significantly (p < 0.0001) higher in AA patients (n = 182) when compared to CA patients (n = 749) (Figure 2C). Since mRNA levels may not always correlate with protein levels, we interrogated the relationship between YB1 mRNA and protein levels in the TCGA PanCancer cohort and found a significant positive correlation (p = 1.77 × 10−15, Spearman) between mRNA and protein levels (Figure S1). We also found YB1 protein expression levels to be significantly (p = 2.03 × 10−3) higher in AA BC tumors when compared to Cas (Figure S2). We further confirmed these findings in a separate TNBC cohort that we acquired from the MetroHealth tumor biobank which contains clinical and survival information on 45 TNBC patients (25 AA and 20 CA). RT-PCR analyses showed that YB1 expression levels were significantly (p = 0.002) higher in AA TNBC tumors when compared to their CA counterparts (Figure 2D). Thus, we show that YB1 expression levels are associated with the AA race.

3.3. Loss of YB1 Inhibits the Oncogenic Behavior of AA TNBC Cell Lines In Vitro

To better understand how the oncogenic activity of YB1 contributes to TNBC disparities, we sought to determine the effect of inhibition of YB1 expression on the behavior of TNBC cell lines of both AA and CA origins. We applied the Synthego CRISPR-based knockout strategy that is based on a ribonucleoprotein (RNP) mix which contains both sgRNAs oligos and the translated Cas9 protein, and is used to deliver the sgRNAs to the nucleus of target cells via electroporation without the need for virus infection and drug selection. We used a pool of three different sgRNAs to target YB1 in three AA and three CA TNBC cell lines. Western blot analyses showed a KO efficiency of more than 90% in all six cell lines (Figure 3A and Figure S5). The KO efficiency was sustained after more than 15 passages in all six lines supporting the effectiveness the KO strategy. In some cases, after serial passaging of the KO cells, we observed progressive re-expression of YB1, and a simple second electroporation with the sgRNA RNP mix was sufficient to eliminate the residual YB1 expression.
The AKT pathway is widely documented to be affected by the loss of YB1 expression [60], which we confirmed in our YB1-KO cells where phosphorylation levels of AKT were decreased in all six cell lines as a result of the loss of expression of YB1 (Figure 3A). Next, we used the colony formation assay as an in vitro indicator assay for the oncogenicity of cancer cells, and found the loss of YB1 expression to significantly inhibit the colony formation potential of all six cell lines (Figure 3B,C). While colony formation was reduced in all six TNBC cell lines, the fold reduction in the AA cell lines ranged from 1.8 to 3 and averaged a 2.4-fold inhibition when compared to the CA cell lines, where the fold reduction ranged from 2.5 to 3.25 with and an average of 2.95-fold inhibition (Figure 3C). These data indicate that the CA TNBC cells are more sensitive to the loss of YB1 activity when compared to their AA counterparts, which show more addiction to YB1.

3.4. Loss of YB1 Inhibits Tumor Growth and Metastasis of AA TNBC Cell Lines In Vivo

To further assess the effects of loss of YB1 expression on tumor growth and metastasis of TNBC cells in vivo, mammary fat pads of female Balb/C mice were inoculated with parental syngeneic 4T1 TNBC cells or their YB1-KO derivatives, and tumor growth was assessed over 5 weeks. Loss of YB1 significantly inhibited the growth of primary tumors (Figure 4A,B). While every mouse implanted with the parental 4T1 or the YB1-KO cells developed tumors (100% incidence), tumor latency was more than 1 week longer in the mice implanted with the YB1-KO 4T1 cells (Figure 4A). Tumor burden, as assessed by tumor volume (Figure 4B), was significantly lower (p  <  0.0001) in mice implanted with the YB1-deficient cells compared to their parental counterparts. Therefore, we show that loss of YB1 expression inhibits tumor growth, in vivo. Proliferation assays in 2D cultures showed no significant difference in cell growth between parental and YB1-KO (Figure S3). Therefore, the differences in tumor growth were rather the result of inhibition of YB1 oncogenic activity. We also used Western Blot analyses of tumor lysates to confirm the sustained loss of YB1 expression in the tumors derived from the YB1-KO cells as these tumors grew in a slower rate as compared to the tumors derived from the parental cells (Figure S4).
Next, we sought to investigate whether inhibition of YB1 expression influences BC metastasis to the lungs. In the spontaneous metastasis assay, after the mice were sacrificed and the lungs were removed, and metastasis foci were counted. The number of lung metastasis foci was more than five-fold (p  <  0.001) lower in mice injected with the YB1-deficient 4T1 cells when compared to their parental counterparts (Figure 4C). These results show that YB1 plays a critical role in the growth, invasion, and metastasis of TNBC tumors, as demonstrated here in the tumor growth and the spontaneous metastasis assays.
Next, we assessed whether loss of YB1 would differentially affect tumor growth and metastasis in human AA vs CA TNBC cell lines. We used one CA TNBC cell line (MDA-MB-231) and one AA TNBC cell line (MDA-MB-468) and the immuno-compromised mouse BC model NOD-scid-IL2Rgamma knockout (NSG) mice. While mice implanted by either the parental MDA-MB-231 or MDA-MB-468 developed tumors (Figure 4D), tumor burden was significantly higher in the mice implanted with the AA MDA-MB-468 TNBC cells (Figure 4D, black graph) when compared to the AA MDA-MB-231 cells (Figure 4D, grey graph), which suggests the more aggressive nature of TNBC tumors of AA origin. Additionally, the loss of YB1 expression inhibited the tumor growth of both MDA-MB-231-YB1-KO cells (Figure 4D, yellow graph) and MDA-MB-468-YB1-KO cells (Figure 4D, red graph). Tumor growth inhibition was, however, more pronounced in the CA MDA-MB-231 cells that are more sensitive to loss of YB1 when compared to AA MDA-MB-468 cells that are more resistant to loss of YB1 (Figure 4D, compare red graph to yellow graph). Thus, our data suggest a difference in sensitivity to loss of YB1 in tumor growth between AA and CA TNBC cells in vivo. These findings were also confirmed with respect to spontaneous lung metastasis (Figure 4E,F), where lung metastasis burden was ~2-fold higher with the AA MDA-MB-468 cells when compared to the CA MDA-MB-231 cells (Compare black bar of Figure 4E to grey bar of Figure 4F). Finally, while loss of YB1 significantly inhibited the metastasis of both AA MDA-MB-468 cells (Figure 4E, red bar) and CA MDA-MB-231 cells (Figure 4F, yellow bar), lung metastasis was reduced by ~3-fold in the AA MDA-MB-468 cells (Figure 4E) when compared to more than 8-fold in the sensitive CA MDA-MB-231 cells (Figure 4F). Together, these data, strongly suggest a critical role of YB1 in the aggressiveness of AA TNBC in vitro and in vivo.

3.5. YB1 Is Associated with the Cancer Stem Cell Phenotype in AA TNBC

In our previously published studies [31,49] we showed that YB1 is required for the regulation of cancer stem cell phenotypes by modulating the expression levels of cancer stem cell (CSC) transcription machinery, including Oct4, Sox2, and Nanog. To determine whether these CSC transcription factors are differentially expressed between AA and CA TNBC tumors and cell lines, we interrogated the TCGA PanCancer cohort and found Oct4 (Figure 5A), Nanog (Figure 5B), and Sox2 (Figure 5C) to be significantly (p < 0.001) highly expressed in tumors of the basal (TNBC) breast cancer subtype when compared to tumors of Her2, Luminal-A, or Luminal-B subtypes. Oct4 (Figure 5D), Nanog (Figure 5E), and Sox2 (Figure 5F) were also significantly disparately highly expressed in AA tumors when compared to their CA counterparts. We have also shown that the WAVE3-YB1 signaling axis is required for the activation of the CSC phenotype in TNBC [31,49]. Here, WAVE3 is also significantly (p < 0.001) highly expressed in the tumors of the basal subtype (Figure 5G) and is significantly (p < 0.001) highly expressed in AA than in CA tumors (Figure 5H). Thus, we show that not only YB1 is highly expressed in AA versus CA TNBC tumors, but YB1-associated CSC genes are also disparately associated with TNBC tumors of AA origin, which may in part explain the aggressivity of these tumors in the AA population. We confirmed these findings in our MetroHealth TNBC cohort, where YB1 (Figure 2E), Oct4 (Figure 5I), Nanog (Figure 5J), Sox2 (Figure 5K), and WAVE3 (Figure 5L) are significantly highly expressed in AA tumors when compared to their CA counterparts.

3.6. YB1 Is Associated with Chemoresistance in AA TNBC Both In Vitro and In Vivo

YB1 was shown to be associated with chemotherapy resistance in several types of cancer, including BC (reviewed in [16]). We interrogated the Cancer Cell Line Encyclopedia (CCLE) for the IC50 of three of the most commonly used chemotherapeutic agents for the treatment of patients with TNBC (Cisplatin, Docetaxel and Doxorubicin) in cell lines of the TNBC subtype; these included three AA and three CA TNBC cell lines. We found the IC50 values, which are indicative of level of resistance to a given drug, to be consistently higher in AA cell lines when compared to their CA counterparts for cisplatin (Figure 6A,B), doxorubicin (Figure 6C,D) and docetaxel (Figure 6E,F). Since YB1 expression levels were also significantly elevated in AA cell lines when compared to their CA counterparts (Figure 2E), these results show a clear positive correlation between YB1 expression and resistance to chemotherapies. Gene set enrichment analysis (GSEA) of the TCGA breast cancer cohort confirmed our findings, where we found a significant (p = 0.003) enrichment of genes that are associated with docetaxel resistance in BC tumors that have increased levels of YB1 expression (Figure 6G) [61]. To investigate the role of YB1 in the chemoresistance in TNBC, we treated MDA-MB-468 with doxorubicin for 48 h. and found that this transient treatment resulted in the inhibition of YB1 expression and its phosphorylation levels (Figure 6H and Figure S6). However, in MDA-MB-468 that were generated to be permanently resistant to doxorubicin (Dox-Res), we found this chronic exposure to doxorubicin to instead activate both the expression and phosphorylated levels of YB1 (Figure 6H and Figure S6). These finding were confirmed in MDA-MB-468 that were generated to be resistant to cisplatin (Cis-Res) and found both the expression and phosphorylation of YB1 to be increased in the Cis-Res cells when compared to their control counterparts (Figure 6I and Figure S6). Thus, we show that chronic exposure to chemotherapy induces YB1 expression and phosphorylation. Next, we injected either the parental or the cisplatin-resistant (Cis-Res) MDA-MB-468 cells into the mammary fat pads of female NSG mice. At 21 days post implantation, mice from both the groups received two doses of cisplatin (1 mg/kg, IP) per week for 3 weeks, and tumor growth was followed for several weeks (Figure 6J). We found that treatment with cisplatin resulted in inhibition of tumor growth in both groups for the first three weeks following the treatment; after this three-week period, the mice that were injected with the Cis-res cells showed a tumor growth that was much faster than the mice injected with the parental control cells. After 70 days of follow-up, the average volume of the tumors derived from the parental cells did not exceed 150 mm3, while that of the tumors derived from the cisplatin-resistant cells rapidly reached the maximal tumor burden of more than 1000 mm3 (Figure 6J). Therefore, we speculate that the increased YB1 expression in the Cis-resistant cells, which is induced by chronic exposure to chemotherapy, confers these cells with an enhanced tumor growth potential and makes them refractory to chemotherapy treatment. Together, our data suggest a novel positive feedback loop between YB1 expression and chemoresistance, therefore creating a vicious oncogenic environment which activates a sustained tumor growth and metastasis, which is more active in AA TNBC tumors, and ultimately contributes to the disparities in health outcomes in AA women with TNBC disease.

4. Discussion

Breast cancer is the second leading cause of cancer-related death in both Caucasian American (CA) and African American (AA) women. The most devastating subtype of BC is TNBC. TNBC accounts for about 15–20% of all invasive breast cancer cases in the United States [62]. By definition, TNBCs lack the expression of estrogen and progesterone hormone receptors and presents with low copies of Her2 (EGFR2). Whereas targeted therapies are available for other breast cancer subtypes, including hormone receptor-positive and HER-2 amplified breast cancers, there aren’t any FDA-approved targeted therapies for TNBC, and treatment options are limited to chemotherapy that is usually ineffective especially in later stages. There have not been any meaningful additions to our TNBC armamentarium for a while, and the recent excitement of a breakthrough, with the accelerated approval of the immune checkpoint inhibitor atezolizumab, was short-lived when the manufacturer voluntarily withdrew the approval. The cancer mortality disparity between AA and CA women is alarming, with death rates being 40% higher in AAs (28.4 vs. 20.3 deaths per 100,000), and in AAs aged less than 50 years, it is tragic with the death rate that is double that of CA and is the highest among all race groups [63]. Reasons for these disparities include AA women presenting at a younger age, where 33% are diagnosed aged less than 50 years when compared to 21.9% of CA women [64]. AA women are also more likely to delay breast cancer screening mammography, as well as to present with an advanced stage, face an average waiting time to the initiation of treatment that is double that of CA women face, and less likely to have surgery [65].
While environmental factors and socioeconomic status play undeniable roles in driving poor outcomes in AAs, a key gap in knowledge remains regarding the biological and molecular mechanisms that interact with these extrinsic factors to impact treatment response and clinical outcomes. Studying disparities in cancer-related morbidity and mortality between different racial and ethnic groups is complex and multifaceted. Examining the potential relationships between social and biological factors is of paramount importance. From the time of the inception of the Surveillance, Epidemiology, and End Results (SEER) program in 1973 until 1983, breast cancer mortality in CA and AA women was almost the same; however, a divergence in mortality started following 1983, which coincided with the introduction of tamoxifen for hormone receptor-positive cancers [66]. A combination of these factors contributed to the fact that AA TNBC patients end up with significantly worse outcomes when compared to their CA counterparts. While environmental factors and socioeconomic status play undeniable roles in driving poor outcomes in AAs, a key gap in knowledge remains regarding the biological and molecular mechanisms that interact with these extrinsic factors to impact treatment response and clinical outcomes. In the present study, we used bioinformatics and biostatistics analyses to interrogate both public and local BC cohorts to show that YB1 may be a driver of BC health disparities in AA women. We found that not only YB1, but its associated cancer stem cell and invasion-metastasis gene signatures, including WAVE3, are also associated with worst disease outcomes in AA patients with TNBC tumors.
We supported our findings with a combination of in vitro genetic manipulations of YB1 (CRISPR/Cas-mediated knockout) and treatment with chemotherapeutic agents, as well as in vivo mouse models of TNBC to show that the loss of YB1 expression inhibits its oncogenic activity both in vitro (colony formation) and in vivo (inhibition of tumor growth and metastasis) in TNBC cell lines of both CA and AA origin. However, inhibition of YB1 oncogenic activity was more pronounced in cell lines of CA origin when compared to their AA counterparts. This suggested to us that YB1 might be more active in AA TNBC cell lines than in their CA counterparts, since the residual YB1 expression in the AA cell lines after CRISPR-KO may be sufficient to continue driving YB1 oncogenic activity in AA cell lines. This suggestion is supported by the fact that YB1 is also involved in several oncogenic activities, including in the regulation of the CSC phenotype [31,45,49,61,67]. In fact, our data show that expression levels of CSC markers Nanog, Sox2, and Oct4, as well as their upstream regulator, WAVE3, are significantly higher in AA TNBC tumors when compared to CAs (Figure 5). Our published studies [31] also showed that loss of either YB1 or WAVE3 expression inhibits the expression of these CSC markers, albeit not completely. Therefore, we posit that the residual levels of YB1 and its associated CSC markers is responsible for the enhanced aggressivity of the AA TNBC cell lines both in vitro (Figure 3) and in vivo (Figure 4). YB1 is also widely documented as a major player in the regulation of chemoresistance in several types of cancer, including BC (reviewed in [16]), which reenforces the association between YB1 and chemoresistance in AA TNBC tumors, given that AA TNBC are more prone to develop resistance to standards of care chemotherapies than their CA counterparts (reviewed in [16]). Our data show a clear and significant positive correlation between YB1 expression and the IC50 values of the chemotherapeutic drugs that are commonly used to treat TNBC tumors, in TNBC cell lines of AA origin (Figure 6). We further identify a positive aggressive feedback loop between YB1 expression and resistance to chemotherapy to drive the tumor growth of AA tumors. While short-term treatment with chemotherapy seems to inhibit both YB1 expression and its phosphorylation levels, chronic exposure to these therapeutic drugs results in the reactivation of YB1 expression and its oncogenic activities: Cisplatin-and doxorubicin-resistant MDA-MB-468 cells express more YB1 and phopho-YB1 when compared to their parental cells, and tumors derived from the cisplatin-resistant AA MDA-MB-468 TNBC cells grow faster and bigger than those derived from their parental counterparts, even after chemotherapy treatment (Figure 6). These data strongly support the interplay between YB1 and chemoresistance that drives the aggressivity of AA TNBC tumors.

5. Conclusions

Taken together, our study supports the role of YB1 as a major molecular contributor to the health outcome disparities observed in AA women with TNBC tumors. A simple PubMed search with the keywords YB1 and cancer resulted in more than 600 publications, which demonstrate the significant interest in investigating YB1 in the context of cancer biology. However, none of these studies attempted to investigate the potential role of YB1 in health disparities. Our study is the first to demonstrate the potential contribution of YB1 to TNBC health disparities in AA women. In fact, the interrogation of public databases of transcriptomic TCGA datasets revealed that YB1 is associated with health disparities in AA women not only in TNBC tumors, but in several other cancers of epithelial origin, such as colon and ovarian cancers. This clearly demonstrates a global association of YB1 with AA health disparities and justifies the urgent need for the development of targeted therapies against YB1 that would benefit not only the general population but the AA patient population, which is disparately affected by the cancer burden.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13246262/s1, Figure S1: Correlation between YB1 mRNA and protein levels in the PanCancer TCGA data set. mRNA levels are plotted as log2 RSEM values batch normalized from Illumina HiSeq RNASeq data, while protein levels are plotted as RPPA, Figure S2: Comparison of YB1 protein phosphorylation levels, based on RPPA, between AA and CA BC patients from the TCGA BRCA PanCancer Atlas cohort. YB1 protein phosphorylation levels are significantly higher (Wilcoxon) in AA compared to CA TNBC tumors, Figure S3: Cell proliferation over 5 days of the indicated TNBC cell lines and their YB1-KO derivatives. Viable cell were Identified with PI staining, Figure S4: Western Blot analysis of protein lysates from tumors derived from the indicated parental (Ctrl) and their YB1-deficient (YB1KO) derivatives TNBC cell lines, that were probed with YB1 antibodies. β-Actin was used a loading control, Figure S5: Uncropped Western blots for panel A of Figure 3, Figure S6: Uncropped WBs for panel E and panel F of Figure 6.

Author Contributions

Conceptualization, N.J. and K.S.-A.; Data curation, R.C.; Formal analysis, K.S.-A.; Funding acquisition, K.S.-A.; Investigation, W.W., A.P. and K.S.-A.; Methodology, P.S.R., W.W., A.A. and V.M.; Project administration, K.S.-A.; Resources, A.P., A.K. and N.J.; Software, R.C.; Supervision, A.P. and N.J.; Validation, P.S.R., W.W. and A.A.; Visualization, R.C.; Writing—original draft, P.S.R. and K.S.-A.; Writing—review & editing, W.W. and V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by NIH grant R01CA226921 and R01CA226921-S1 to K. Sossey-Alaoui.

Institutional Review Board Statement

All animals were maintained according to the Institutional Animal Care and Use Committee (IACUC) guidelines (IACUC # 2019-0001; approval date: April 17 2019) of Case Western Reserve University and MetroHealth System, Cleveland, OH, USA.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank members of Sossey-Alaoui Lab for their valuable suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nguyen, D.X.; Bos, P.D.; Massague, J. Metastasis: From dissemination to organ-specific colonization. Nat. Rev. Cancer 2009, 9, 274–284. [Google Scholar] [CrossRef] [PubMed]
  2. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer Statistics, 2021. CA Cancer J. Clin. 2021, 71, 7–33. [Google Scholar] [CrossRef] [PubMed]
  3. Perou, C.M.; Sorlie, T.; Eisen, M.B.; van de Rijn, M.; Jeffrey, S.S.; Rees, C.A.; Pollack, J.R.; Ross, D.T.; Johnsen, H.; Akslen, L.A.; et al. Molecular portraits of human breast tumours. Nature 2000, 406, 747–752. [Google Scholar] [CrossRef] [PubMed]
  4. Sorlie, T.; Perou, C.M.; Tibshirani, R.; Aas, T.; Geisler, S.; Johnsen, H.; Hastie, T.; Eisen, M.B.; van de Rijn, M.; Jeffrey, S.S.; et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc. Natl. Acad. Sci. USA 2001, 98, 10869–10874. [Google Scholar] [CrossRef] [Green Version]
  5. Sorlie, T.; Tibshirani, R.; Parker, J.; Hastie, T.; Marron, J.S.; Nobel, A.; Deng, S.; Johnsen, H.; Pesich, R.; Geisler, S.; et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc. Natl. Acad. Sci. USA 2003, 100, 8418–8423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lehmann, B.D.; Bauer, J.A.; Chen, X.; Sanders, M.E.; Chakravarthy, A.B.; Shyr, Y.; Pietenpol, J.A. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J. Clin. Invest. 2011, 121, 2750–2767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Chen, X.; Li, J.; Gray, W.H.; Lehmann, B.D.; Bauer, J.A.; Shyr, Y.; Pietenpol, J.A. TNBCtype: A Subtyping Tool for Triple-Negative Breast Cancer. Cancer Inform. 2012, 11, 147–156. [Google Scholar] [CrossRef] [PubMed]
  8. Finnegan, T.J.; Carey, L.A. Gene-expression analysis and the basal-like breast cancer subtype. Future Oncol. 2007, 3, 55–63. [Google Scholar] [CrossRef] [PubMed]
  9. Schneider, B.P.; Winer, E.P.; Foulkes, W.D.; Garber, J.; Perou, C.M.; Richardson, A.; Sledge, G.W.; Carey, L.A. Triple-negative breast cancer: Risk factors to potential targets. Clin. Cancer Res. 2008, 14, 8010–8018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Anders, C.K.; Carey, L.A. Biology, metastatic patterns, and treatment of patients with triple-negative breast cancer. Clin. Breast Cancer 2009, 9, S73–S81. [Google Scholar] [CrossRef] [PubMed]
  11. Carey, L.; Winer, E.; Viale, G.; Cameron, D.; Gianni, L. Triple-negative breast cancer: Disease entity or title of convenience? Nat. Rev. Clin. Oncol. 2010, 7, 683–692. [Google Scholar] [CrossRef] [PubMed]
  12. Foulkes, W.D.; Smith, I.E.; Reis-Filho, J.S. Triple-negative breast cancer. N. Engl. J. Med. 2010, 363, 1938–1948. [Google Scholar] [CrossRef] [Green Version]
  13. Dawood, S.; Broglio, K.; Kau, S.W.; Green, M.C.; Giordano, S.H.; Meric-Bernstam, F.; Buchholz, T.A.; Albarracin, C.; Yang, W.T.; Hennessy, B.T.; et al. Triple receptor-negative breast cancer: The effect of race on response to primary systemic treatment and survival outcomes. J. Clin. Oncol. 2009, 27, 220–226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Albain, K.S.; Unger, J.M.; Crowley, J.J.; Coltman, C.A., Jr.; Hershman, D.L. Racial disparities in cancer survival among randomized clinical trials patients of the Southwest Oncology Group. J. Natl. Cancer Inst. 2009, 101, 984–992. [Google Scholar] [CrossRef] [Green Version]
  15. Dietze, E.C.; Sistrunk, C.; Miranda-Carboni, G.; O’Regan, R.; Seewaldt, V.L. Triple-negative breast cancer in African-American women: Disparities versus biology. Nat. Rev. Cancer 2015, 15, 248–254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Alkrekshi, A.; Wang, W.; Rana, P.S.; Markovic, V.; Sossey-Alaoui, K. A comprehensive review of the functions of YB-1 in cancer stemness, metastasis and drug resistance. Cell Signal. 2021, 85, 110073. [Google Scholar] [CrossRef] [PubMed]
  17. Tiwari, A.; Iida, M.; Kosnopfel, C.; Abbariki, M.; Menegakis, A.; Fehrenbacher, B.; Maier, J.; Schaller, M.; Brucker, S.Y.; Wheeler, D.L.; et al. Blocking Y-Box Binding Protein-1 through Simultaneous Targeting of PI3K and MAPK in Triple Negative Breast Cancers. Cancers 2020, 12, 2795. [Google Scholar] [CrossRef]
  18. Wang, X.; Li, L.; Zhao, K.; Lin, Q.; Li, H.; Xue, X.; Ge, W.; He, H.; Liu, D.; Xie, H.; et al. A novel LncRNA HITT forms a regulatory loop with HIF-1alpha to modulate angiogenesis and tumor growth. Cell Death Differ. 2020, 27, 1431–1446. [Google Scholar] [CrossRef]
  19. Anastasiadou, E.; Jacob, L.S.; Slack, F.J. Non-coding RNA networks in cancer. Nat. Rev. Cancer 2018, 18, 5–18. [Google Scholar] [CrossRef] [PubMed]
  20. Kosnopfel, C.; Sinnberg, T.; Sauer, B.; Niessner, H.; Schmitt, A.; Makino, E.; Forschner, A.; Hailfinger, S.; Garbe, C.; Schittek, B. Human melanoma cells resistant to MAPK inhibitors can be effectively targeted by inhibition of the p90 ribosomal S6 kinase. Oncotarget 2017, 8, 35761–35775. [Google Scholar] [CrossRef] [Green Version]
  21. Wu, S.L.; Fu, X.; Huang, J.; Jia, T.T.; Zong, F.Y.; Mu, S.R.; Zhu, H.; Yan, Y.; Qiu, S.; Wu, Q.; et al. Genome-wide analysis of YB-1-RNA interactions reveals a novel role of YB-1 in miRNA processing in glioblastoma multiforme. Nucleic Acids Res. 2015, 43, 8516–8528. [Google Scholar] [CrossRef] [PubMed]
  22. Luo, C.; Tetteh, P.W.; Merz, P.R.; Dickes, E.; Abukiwan, A.; Hotz-Wagenblatt, A.; Holland-Cunz, S.; Sinnberg, T.; Schittek, B.; Schadendorf, D.; et al. miR-137 inhibits the invasion of melanoma cells through downregulation of multiple oncogenic target genes. J. Invest. Derm. 2013, 133, 768–775. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Evdokimova, V.; Tognon, C.; Ng, T.; Ruzanov, P.; Melnyk, N.; Fink, D.; Sorokin, A.; Ovchinnikov, L.P.; Davicioni, E.; Triche, T.J.; et al. Translational activation of snail1 and other developmentally regulated transcription factors by YB-1 promotes an epithelial-mesenchymal transition. Cancer Cell 2009, 15, 402–415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Stratford, A.L.; Fry, C.J.; Desilets, C.; Davies, A.H.; Cho, Y.Y.; Li, Y.; Dong, Z.; Berquin, I.M.; Roux, P.P.; Dunn, S.E. Y-box binding protein-1 serine 102 is a downstream target of p90 ribosomal S6 kinase in basal-like breast cancer cells. Breast Cancer Res. 2008, 10, R99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Evdokimova, V.; Ruzanov, P.; Anglesio, M.S.; Sorokin, A.V.; Ovchinnikov, L.P.; Buckley, J.; Triche, T.J.; Sonenberg, N.; Sorensen, P.H. Akt-mediated YB-1 phosphorylation activates translation of silent mRNA species. Mol. Cell. Biol. 2006, 26, 277–292. [Google Scholar] [CrossRef] [Green Version]
  26. Sutherland, B.W.; Kucab, J.; Wu, J.; Lee, C.; Cheang, M.C.; Yorida, E.; Turbin, D.; Dedhar, S.; Nelson, C.; Pollak, M.; et al. Akt phosphorylates the Y-box binding protein 1 at Ser102 located in the cold shock domain and affects the anchorage-independent growth of breast cancer cells. Oncogene 2005, 24, 4281–4292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Evdokimova, V.; Ruzanov, P.; Imataka, H.; Raught, B.; Svitkin, Y.; Ovchinnikov, L.P.; Sonenberg, N. The major mRNA-associated protein YB-1 is a potent 5’ cap-dependent mRNA stabilizer. EMBO J. 2001, 20, 5491–5502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Ceman, S.; Nelson, R.; Warren, S.T. Identification of mouse YB1/p50 as a component of the FMRP-associated mRNP particle. Biochem. Biophys. Res. Commun. 2000, 279, 904–908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Bargou, R.C.; Jurchott, K.; Wagener, C.; Bergmann, S.; Metzner, S.; Bommert, K.; Mapara, M.Y.; Winzer, K.J.; Dietel, M.; Dorken, B.; et al. Nuclear localization and increased levels of transcription factor YB-1 in primary human breast cancers are associated with intrinsic MDR1 gene expression. Nat. Med. 1997, 3, 447–450. [Google Scholar] [CrossRef] [PubMed]
  30. Bergmann, S.; Royer-Pokora, B.; Fietze, E.; Jurchott, K.; Hildebrandt, B.; Trost, D.; Leenders, F.; Claude, J.C.; Theuring, F.; Bargou, R.; et al. YB-1 provokes breast cancer through the induction of chromosomal instability that emerges from mitotic failure and centrosome amplification. Cancer Res. 2005, 65, 4078–4087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Bledzka, K.; Schiemann, B.; Schiemann, W.P.; Fox, P.; Plow, E.F.; Sossey-Alaoui, K. The WAVE3-YB1 interaction regulates cancer stem cells activity in breast cancer. Oncotarget 2017, 8, 104072–104089. [Google Scholar] [CrossRef]
  32. Darb-Esfahani, S.; Loibl, S.; Muller, B.M.; Roller, M.; Denkert, C.; Komor, M.; Schluns, K.; Blohmer, J.U.; Budczies, J.; Gerber, B.; et al. Identification of biology-based breast cancer types with distinct predictive and prognostic features: Role of steroid hormone and HER2 receptor expression in patients treated with neoadjuvant anthracycline/taxane-based chemotherapy. Breast Cancer Res. 2009, 11, R69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dhillon, J.; Astanehe, A.; Lee, C.; Fotovati, A.; Hu, K.; Dunn, S.E. The expression of activated Y-box binding protein-1 serine 102 mediates trastuzumab resistance in breast cancer cells by increasing CD44+ cells. Oncogene 2010, 29, 6294–6300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Fujii, T.; Seki, N.; Namoto-Matsubayashi, R.; Takahashi, H.; Inoue, Y.; Toh, U.; Kage, M.; Shirouzu, K. YB-1 prevents apoptosis via the mTOR/STAT3 pathway in HER-2-overexpressing breast cancer cells. Future Oncol. 2009, 5, 153–156. [Google Scholar] [CrossRef] [PubMed]
  35. Fujita, T.; Ito, K.; Izumi, H.; Kimura, M.; Sano, M.; Nakagomi, H.; Maeno, K.; Hama, Y.; Shingu, K.; Tsuchiya, S.; et al. Increased nuclear localization of transcription factor Y-box binding protein 1 accompanied by up-regulation of P-glycoprotein in breast cancer pretreated with paclitaxel. Clin. Cancer Res. 2005, 11, 8837–8844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Guay, D.; Evoy, A.A.; Paquet, E.; Garand, C.; Bachvarova, M.; Bachvarov, D.; Lebel, M. The strand separation and nuclease activities associated with YB-1 are dispensable for cisplatin resistance but overexpression of YB-1 in MCF7 and MDA-MB-231 breast tumor cells generates several chemoresistance signatures. Int. J. Biochem. Cell Biol. 2008, 40, 2492–2507. [Google Scholar] [CrossRef] [PubMed]
  37. Habibi, G.; Leung, S.; Law, J.H.; Gelmon, K.; Masoudi, H.; Turbin, D.; Pollak, M.; Nielsen, T.O.; Huntsman, D.; Dunn, S.E. Redefining prognostic factors for breast cancer: YB-1 is a stronger predictor of relapse and disease-specific survival than estrogen receptor or HER-2 across all tumor subtypes. Breast Cancer Res. 2008, 10, R86. [Google Scholar] [CrossRef] [Green Version]
  38. Ito, T.; Kamijo, S.; Izumi, H.; Kohno, K.; Amano, J.; Ito, K. Alteration of Y-box binding protein-1 expression modifies the response to endocrine therapy in estrogen receptor-positive breast cancer. Breast Cancer Res. Treat. 2012, 133, 145–159. [Google Scholar] [CrossRef] [PubMed]
  39. Janz, M.; Harbeck, N.; Dettmar, P.; Berger, U.; Schmidt, A.; Jurchott, K.; Schmitt, M.; Royer, H.D. Y-box factor YB-1 predicts drug resistance and patient outcome in breast cancer independent of clinically relevant tumor biologic factors HER2, uPA and PAI-1. Int. J. Cancer 2002, 97, 278–282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Lee, C.; Dhillon, J.; Wang, M.Y.; Gao, Y.; Hu, K.; Park, E.; Astanehe, A.; Hung, M.C.; Eirew, P.; Eaves, C.J.; et al. Targeting YB-1 in HER-2 overexpressing breast cancer cells induces apoptosis via the mTOR/STAT3 pathway and suppresses tumor growth in mice. Cancer Res. 2008, 68, 8661–8666. [Google Scholar] [CrossRef] [Green Version]
  41. Lettau, K.; Zips, D.; Toulany, M. Simultaneous Targeting of RSK and AKT Efficiently Inhibits YB-1-Mediated Repair of Ionizing Radiation-Induced DNA Double-Strand Breaks in Breast Cancer Cells. Int. J. Radiat. Oncol. Biol. Phys. 2021, 109, 567–580. [Google Scholar] [CrossRef] [PubMed]
  42. Lim, J.P.; Nair, S.; Shyamasundar, S.; Chua, P.J.; Muniasamy, U.; Matsumoto, K.; Gunaratne, J.; Bay, B.H. Silencing Y-box binding protein-1 inhibits triple-negative breast cancer cell invasiveness via regulation of MMP1 and beta-catenin expression. Cancer Lett. 2019, 452, 119–131. [Google Scholar] [CrossRef] [PubMed]
  43. Lim, J.P.; Shyamasundar, S.; Gunaratne, J.; Scully, O.J.; Matsumoto, K.; Bay, B.H. YBX1 gene silencing inhibits migratory and invasive potential via CORO1C in breast cancer in vitro. BMC Cancer 2017, 17, 201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Stope, M.B.; Weiss, M.; Koensgen, D.; Popp, S.L.; Joffroy, C.; Mustea, A.; Buck, M.B.; Knabbe, C. Y-box Binding Protein-1 Enhances Oncogenic Transforming Growth Factor beta Signaling in Breast Cancer Cells via Triggering Phospho-Activation of Smad2. Anticancer Res. 2017, 37, 6745–6748. [Google Scholar] [CrossRef] [PubMed]
  45. Yang, F.; Chen, S.; He, S.; Huo, Q.; Hu, Y.; Xie, N. YB-1 interplays with ERalpha to regulate the stemness and differentiation of ER-positive breast cancer stem cells. Theranostics 2020, 10, 3816–3832. [Google Scholar] [CrossRef] [PubMed]
  46. Davuluri, G.; Augoff, K.; Schiemann, W.P.; Plow, E.F.; Sossey-Alaoui, K. WAVE3-NFkappaB interplay is essential for the survival and invasion of cancer cells. PLoS ONE 2014, 9, e110627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Davuluri, G.; Schiemann, W.P.; Plow, E.F.; Sossey-Alaoui, K. Loss of WAVE3 sensitizes triple-negative breast cancers to chemotherapeutics by inhibiting the STAT-HIF-1alpha-mediated angiogenesis. JAKSTAT 2014, 3, e1009276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Sossey-Alaoui, K.; Safina, A.; Li, X.; Vaughan, M.M.; Hicks, D.G.; Bakin, A.V.; Cowell, J.K. Down-regulation of WAVE3, a metastasis promoter gene, inhibits invasion and metastasis of breast cancer cells. Am. J. Pathol. 2007, 170, 2112–2121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Kansakar, U.; Wang, W.; Markovic, V.; Sossey-Alaoui, K. Phosphorylation of the proline-rich domain of WAVE3 drives its oncogenic activity in breast cancer. Sci. Rep. 2021, 11, 3868. [Google Scholar] [CrossRef] [PubMed]
  50. Wang, W.; Li, Y.; Zhu, J.Y.; Fang, D.; Ding, H.F.; Dong, Z.; Jing, Q.; Su, S.B.; Huang, S. Triple negative breast cancer development can be selectively suppressed by sustaining an elevated level of cellular cyclic AMP through simultaneously blocking its efflux and decomposition. Oncotarget 2016, 7, 87232–87245. [Google Scholar] [CrossRef] [Green Version]
  51. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  52. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Cancer Genome Atlas Research, N.; Weinstein, J.N.; Collisson, E.A.; Mills, G.B.; Shaw, K.R.; Ozenberger, B.A.; Ellrott, K.; Shmulevich, I.; Sander, C.; Stuart, J.M. The Cancer Genome Atlas Pan-Cancer analysis project. Nat. Genet. 2013, 45, 1113–1120. [Google Scholar] [CrossRef]
  54. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Santner, S.J.; Dawson, P.J.; Tait, L.; Soule, H.D.; Eliason, J.; Mohamed, A.N.; Wolman, S.R.; Heppner, G.H.; Miller, F.R. Malignant MCF10CA1 cell lines derived from premalignant human breast epithelial MCF10AT cells. Breast Cancer Res. Treat. 2001, 65, 101–110. [Google Scholar] [CrossRef]
  57. Aslakson, C.J.; Miller, F.R. Selective events in the metastatic process defined by analysis of the sequential dissemination of subpopulations of a mouse mammary tumor. Cancer Res. 1992, 52, 1399–1405. [Google Scholar]
  58. Wendt, M.K.; Smith, J.A.; Schiemann, W.P. Transforming growth factor-beta-induced epithelial-mesenchymal transition facilitates epidermal growth factor-dependent breast cancer progression. Oncogene 2010, 29, 6485–6498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Wendt, M.K.; Taylor, M.A.; Schiemann, B.J.; Schiemann, W.P. Down-regulation of epithelial cadherin is required to initiate metastatic outgrowth of breast cancer. Mol. Biol. Cell 2011, 22, 2423–2435. [Google Scholar] [CrossRef] [PubMed]
  60. Bader, A.G.; Felts, K.A.; Jiang, N.; Chang, H.W.; Vogt, P.K. Y box-binding protein 1 induces resistance to oncogenic transformation by the phosphatidylinositol 3-kinase pathway. Proc. Natl. Acad. Sci. USA 2003, 100, 12384–12389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Honma, K.; Iwao-Koizumi, K.; Takeshita, F.; Yamamoto, Y.; Yoshida, T.; Nishio, K.; Nagahara, S.; Kato, K.; Ochiya, T. RPN2 gene confers docetaxel resistance in breast cancer. Nat. Med. 2008, 14, 939–948. [Google Scholar] [CrossRef] [PubMed]
  62. Prakash, O.; Hossain, F.; Danos, D.; Lassak, A.; Scribner, R.; Miele, L. Racial Disparities in Triple Negative Breast Cancer: A Review of the Role of Biologic and Non-biologic Factors. Front. Public Health 2020, 8, 576964. [Google Scholar] [CrossRef] [PubMed]
  63. DeSantis, C.E.; Ma, J.; Gaudet, M.M.; Newman, L.A.; Miller, K.D.; Goding Sauer, A.; Jemal, A.; Siegel, R.L. Breast cancer statistics, 2019. CA Cancer J. Clin. 2019, 69, 438–451. [Google Scholar] [CrossRef]
  64. Daly, B.; Olopade, O.I. A perfect storm: How tumor biology, genomics, and health care delivery patterns collide to create a racial survival disparity in breast cancer and proposed interventions for change. CA Cancer J. Clin. 2015, 65, 221–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Foy, K.C.; Fisher, J.L.; Lustberg, M.B.; Gray, D.M.; DeGraffinreid, C.R.; Paskett, E.D. Disparities in breast cancer tumor characteristics, treatment, time to treatment, and survival probability among African American and white women. NPJ Breast Cancer 2018, 4, 7. [Google Scholar] [CrossRef] [PubMed]
  66. Newman, L.A. Parsing the Etiology of Breast Cancer Disparities. J. Clin. Oncol. 2016, 34, 1013–1014. [Google Scholar] [CrossRef] [PubMed]
  67. Davies, A.H.; Reipas, K.; Hu, K.; Berns, R.; Firmino, N.; Stratford, A.L.; Dunn, S.E. Inhibition of RSK with the novel small-molecule inhibitor LJI308 overcomes chemoresistance by eliminating cancer stem cells. Oncotarget 2015, 6, 20570–20577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 13 06262 g001 550
Figure 1. (A) Quantification of Yb1 mRNA expression levels based on log2 RSEM values batch normalized from Illumina HiSeq RNASeq data by breast cancer subtype in the breast cancer BRCA patient data from TCGA PanCancer Atlas. Expression levels of YB1 are significantly higher (*** p < 0.0001, Wilcoxon) in the basal (TNBC) subtype when compared to other BC subtypes. (B) Survival analysis based on YB1 expression levels (median values) in the basal subtype of the TCGA BRCA PanCancer Atlas cohort. YB1 activation in TNBCs is associated with poorer overall survival (p = 0.049). (C) KM plot correlating survival of 4929 BC patients with YB1 mRNA expression levels. High YB1 expression levels correlate with poor survival probability in BC patients (p < 1.9 × 10−8). (D) KM plot correlating the survival of BC patients with YB1 protein expression levels. High YB1 expression levels correlate with poor survival probability in BC patients (p = 0.047). (E) KM plot correlating the survival of TNBC patients with YB1 mRNA expression levels. High YB1 expression levels correlate with poor survival probability in TNBC patients (p = 0.03).
Figure 1. (A) Quantification of Yb1 mRNA expression levels based on log2 RSEM values batch normalized from Illumina HiSeq RNASeq data by breast cancer subtype in the breast cancer BRCA patient data from TCGA PanCancer Atlas. Expression levels of YB1 are significantly higher (*** p < 0.0001, Wilcoxon) in the basal (TNBC) subtype when compared to other BC subtypes. (B) Survival analysis based on YB1 expression levels (median values) in the basal subtype of the TCGA BRCA PanCancer Atlas cohort. YB1 activation in TNBCs is associated with poorer overall survival (p = 0.049). (C) KM plot correlating survival of 4929 BC patients with YB1 mRNA expression levels. High YB1 expression levels correlate with poor survival probability in BC patients (p < 1.9 × 10−8). (D) KM plot correlating the survival of BC patients with YB1 protein expression levels. High YB1 expression levels correlate with poor survival probability in BC patients (p = 0.047). (E) KM plot correlating the survival of TNBC patients with YB1 mRNA expression levels. High YB1 expression levels correlate with poor survival probability in TNBC patients (p = 0.03).
Cancers 13 06262 g001
Cancers 13 06262 g002 550
Figure 2. (A) YB1 mRNA transcript quantification in AA and CA TNBC cell lines was determined by qt-RT-PCR and data was plotted as the fold change in HCC1806. (B) Average YB1 mRNA transcripts in AA and CA TNBC cell lines was plotted as the fold change in AA cell lines. YB1 transcript expression levels are significantly higher (* p < 0.001, Student ’t test) in AA TNBC cell lines when compared to their CA counterparts. (C) Comparison of YB1 expression levels, levels based on log2 RSEM values batch normalized from Illumina HiSeq RNASeq data, between AA and CA BC patients from the TCGA BRCA PanCancer Atlas cohort. YB1 expression levels are significantly higher (* p < 0.01, Wilcoxon) in AA when compared to CA BC tumors. (D) Comparison of YB1 expression levels between CA and AA TNBC from the MetroHealth TNBC cohort, as determined by qt-RT-PCR. Data was plotted as the fold change to AA tumors. YB1 expression levels are significantly higher (* p < 0.01, Student ’t test) in AA when compared to CA BC tumors.
Figure 2. (A) YB1 mRNA transcript quantification in AA and CA TNBC cell lines was determined by qt-RT-PCR and data was plotted as the fold change in HCC1806. (B) Average YB1 mRNA transcripts in AA and CA TNBC cell lines was plotted as the fold change in AA cell lines. YB1 transcript expression levels are significantly higher (* p < 0.001, Student ’t test) in AA TNBC cell lines when compared to their CA counterparts. (C) Comparison of YB1 expression levels, levels based on log2 RSEM values batch normalized from Illumina HiSeq RNASeq data, between AA and CA BC patients from the TCGA BRCA PanCancer Atlas cohort. YB1 expression levels are significantly higher (* p < 0.01, Wilcoxon) in AA when compared to CA BC tumors. (D) Comparison of YB1 expression levels between CA and AA TNBC from the MetroHealth TNBC cohort, as determined by qt-RT-PCR. Data was plotted as the fold change to AA tumors. YB1 expression levels are significantly higher (* p < 0.01, Student ’t test) in AA when compared to CA BC tumors.
Cancers 13 06262 g002
Cancers 13 06262 g003 550
Figure 3. (A) Western Blot analysis of protein lysates from parental (Ctrl) AA (black font) and CA (red font) TNBC cell lines and their YB1-KO derivatives, that were probed with the indicated antibodies. β-Actin was used a loading control. (B,C) Colony formation assay. The numbers under the bands represent the fold change of the Western Blot signal normalized to that of the parental cell line, as determined by densitometry analyses from one representative of 3 replicate blots. (B) Representative images of colony formation of parental and YB1-KO AA (black font) and CA (red font) TNBC cell lines. (C) Quantification of the number of colonies in each plate. Data is plotted as the percentage change from the control cells. Data shown are representative of 3 replicates (***, p < 0.0001, Student ’t test).
Figure 3. (A) Western Blot analysis of protein lysates from parental (Ctrl) AA (black font) and CA (red font) TNBC cell lines and their YB1-KO derivatives, that were probed with the indicated antibodies. β-Actin was used a loading control. (B,C) Colony formation assay. The numbers under the bands represent the fold change of the Western Blot signal normalized to that of the parental cell line, as determined by densitometry analyses from one representative of 3 replicate blots. (B) Representative images of colony formation of parental and YB1-KO AA (black font) and CA (red font) TNBC cell lines. (C) Quantification of the number of colonies in each plate. Data is plotted as the percentage change from the control cells. Data shown are representative of 3 replicates (***, p < 0.0001, Student ’t test).
Cancers 13 06262 g003
Cancers 13 06262 g004 550
Figure 4. (A) Incidence (top) and latency (bottom) of tumor growth of prenatal 4T1 or its YB1-KO derivative cells that were inoculated in the mammary fat pads of female Balb/C mice (n = 5 mice per group). (B) Quantification of tumor volume of tumors generated from the inoculation of control 4T1 cells or their YB1-KO derivatives into the mammary fat pads of female Balb/C mice. (C, top panel) Quantification of lung metastasis nodules from the animal experiment described in (B). The bottom panel shows representative lungs from each mouse group. (D) Quantification of the tumor volume of tumors generated from the inoculation of control MDA-MB-231 (grey graph), MDA-MB-468 (black graph) cells or their respective YB1-KO derivatives (yellow and red graphs, respectively) into the mammary fat pads of female NSG mice (n = 5 mice per group). (E,F) Quantification of lung metastasis nodules from the MDA-MB-468 (E) and the MDA-MB-231 animal experiment described in (D) (**, p < 0.001, Student ’t-test).
Figure 4. (A) Incidence (top) and latency (bottom) of tumor growth of prenatal 4T1 or its YB1-KO derivative cells that were inoculated in the mammary fat pads of female Balb/C mice (n = 5 mice per group). (B) Quantification of tumor volume of tumors generated from the inoculation of control 4T1 cells or their YB1-KO derivatives into the mammary fat pads of female Balb/C mice. (C, top panel) Quantification of lung metastasis nodules from the animal experiment described in (B). The bottom panel shows representative lungs from each mouse group. (D) Quantification of the tumor volume of tumors generated from the inoculation of control MDA-MB-231 (grey graph), MDA-MB-468 (black graph) cells or their respective YB1-KO derivatives (yellow and red graphs, respectively) into the mammary fat pads of female NSG mice (n = 5 mice per group). (E,F) Quantification of lung metastasis nodules from the MDA-MB-468 (E) and the MDA-MB-231 animal experiment described in (D) (**, p < 0.001, Student ’t-test).
Cancers 13 06262 g004
Cancers 13 06262 g005 550
Figure 5. (A–D) Comparison of the expression levels of Oct4 (A), Nanog (B), Sox2 (C), and WAVE3 (D) between the BC subtypes in the TCGA BRCA PanCancer Atlas cohort. (**, p < 0.001; ***, p < 0.0001, Wilcoxon). (E–F) Comparison of the expression levels of Oct4 (E), Nanog (F), Sox2 (G), and WAVE3 (H) between AA and CA BC patients in the TCGA BRCA PanCancer Atlas cohort. (*, p < 0.01; **, p < 0.001; ***, p < 0.0001, Wilcoxon). In all cases, mRNA expression levels were based on log2 RSEM values batch normalized from Illumina HiSeq RNASeq data. (I–L) Comparison of the expression levels of Oct4 (I), Nanog (J), Sox2 (K), and WAVE3 (L) between AA and CA BC patients in the MetroHealth cohort, as determined by qt-RT-PCR. Data were plotted as fold change to AA tumors. (* p < 0.01, **, p < 0.001; ***, p < 0.0001, Wilcoxon).
Figure 5. (A–D) Comparison of the expression levels of Oct4 (A), Nanog (B), Sox2 (C), and WAVE3 (D) between the BC subtypes in the TCGA BRCA PanCancer Atlas cohort. (**, p < 0.001; ***, p < 0.0001, Wilcoxon). (E–F) Comparison of the expression levels of Oct4 (E), Nanog (F), Sox2 (G), and WAVE3 (H) between AA and CA BC patients in the TCGA BRCA PanCancer Atlas cohort. (*, p < 0.01; **, p < 0.001; ***, p < 0.0001, Wilcoxon). In all cases, mRNA expression levels were based on log2 RSEM values batch normalized from Illumina HiSeq RNASeq data. (I–L) Comparison of the expression levels of Oct4 (I), Nanog (J), Sox2 (K), and WAVE3 (L) between AA and CA BC patients in the MetroHealth cohort, as determined by qt-RT-PCR. Data were plotted as fold change to AA tumors. (* p < 0.01, **, p < 0.001; ***, p < 0.0001, Wilcoxon).
Cancers 13 06262 g005
Cancers 13 06262 g006 550
Figure 6. (A–E) Comparisons of the IC50 values of cisplatin (A,B), doxorubicin (C,D) and docetaxel (E,F) between AA (black bars) and CA (red bars) TNBC cell lines. The IC50 values were derived from the CCLE database. Average IC50 in the AA and CA TNBC cell lines was plotted as fold change to AA cell lines for cisplatin (B), doxorubicin (D) and docetaxel (F), (*, p < 0.01, Fisher’s Exact). (G) Gene set enrichment analysis of the TCGA BRCA PanCancer Atlas dataset shows that the gene signature associated with docetaxel resistance is enriched in the BC tumors that express high levels of YB1. (H) Western Blots with the indicated antibodies of protein lysates from parental MDA-MB-MB-468 cells or those that were treated with doxorubicin for 48 h (transient treatment), or MDA-MB-468 cells that were generated to be permanently resistant (chronic) to doxorubicin (468-Dox-res). The numbers under the bands represent the fold change of the Western Blot signal normalized to that of the parental cell line. (I) Western Blots with the indicated antibodies of protein lysates from parental MDA-MB-468, their cisplatin resistant (Cis-R) or doxorubicin (Dox-R) derivatives. β-Actin was used as a loading control. The numbers under the bands represent the fold change of the Western Blot signal normalized to that of the parental cell line. (J) Quantification of the tumor volume of tumors generated from the inoculation of parental MDA-MB-468 cells (black graph), or their cisplatin resistant (Cis-R) derivatives into the mammary fat pads of female NSG mice (n = 5 mice per group). Treatment with cisplatin started 21 days post inoculation. (***, p < 0.0001, Student ’t test).
Figure 6. (A–E) Comparisons of the IC50 values of cisplatin (A,B), doxorubicin (C,D) and docetaxel (E,F) between AA (black bars) and CA (red bars) TNBC cell lines. The IC50 values were derived from the CCLE database. Average IC50 in the AA and CA TNBC cell lines was plotted as fold change to AA cell lines for cisplatin (B), doxorubicin (D) and docetaxel (F), (*, p < 0.01, Fisher’s Exact). (G) Gene set enrichment analysis of the TCGA BRCA PanCancer Atlas dataset shows that the gene signature associated with docetaxel resistance is enriched in the BC tumors that express high levels of YB1. (H) Western Blots with the indicated antibodies of protein lysates from parental MDA-MB-MB-468 cells or those that were treated with doxorubicin for 48 h (transient treatment), or MDA-MB-468 cells that were generated to be permanently resistant (chronic) to doxorubicin (468-Dox-res). The numbers under the bands represent the fold change of the Western Blot signal normalized to that of the parental cell line. (I) Western Blots with the indicated antibodies of protein lysates from parental MDA-MB-468, their cisplatin resistant (Cis-R) or doxorubicin (Dox-R) derivatives. β-Actin was used as a loading control. The numbers under the bands represent the fold change of the Western Blot signal normalized to that of the parental cell line. (J) Quantification of the tumor volume of tumors generated from the inoculation of parental MDA-MB-468 cells (black graph), or their cisplatin resistant (Cis-R) derivatives into the mammary fat pads of female NSG mice (n = 5 mice per group). Treatment with cisplatin started 21 days post inoculation. (***, p < 0.0001, Student ’t test).
Cancers 13 06262 g006
Table
Table 1. Breast cancer patients’ tumors breakdown by race and BC subtype from the TCGA PanCancer Atlas database.
Table 1. Breast cancer patients’ tumors breakdown by race and BC subtype from the TCGA PanCancer Atlas database.
RaceBasalHer2LumALumBNANormal LikeGrand Total
AfAm53155626239182
Asian71520161160
NA51039296190
Native American 1 1
White106373841267125749
Grand Total17178499197101361082
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rana, P.S.; Wang, W.; Alkrekshi, A.; Markovic, V.; Khiyami, A.; Chan, R.; Perzynski, A.; Joseph, N.; Sossey-Alaoui, K. YB1 Is a Major Contributor to Health Disparities in Triple Negative Breast Cancer. Cancers 2021, 13, 6262. https://doi.org/10.3390/cancers13246262

AMA Style

Rana PS, Wang W, Alkrekshi A, Markovic V, Khiyami A, Chan R, Perzynski A, Joseph N, Sossey-Alaoui K. YB1 Is a Major Contributor to Health Disparities in Triple Negative Breast Cancer. Cancers. 2021; 13(24):6262. https://doi.org/10.3390/cancers13246262

Chicago/Turabian Style

Rana, Priyanka Shailendra, Wei Wang, Akram Alkrekshi, Vesna Markovic, Amer Khiyami, Ricky Chan, Adam Perzynski, Natalie Joseph, and Khalid Sossey-Alaoui. 2021. "YB1 Is a Major Contributor to Health Disparities in Triple Negative Breast Cancer" Cancers 13, no. 24: 6262. https://doi.org/10.3390/cancers13246262

APA Style

Rana, P. S., Wang, W., Alkrekshi, A., Markovic, V., Khiyami, A., Chan, R., Perzynski, A., Joseph, N., & Sossey-Alaoui, K. (2021). YB1 Is a Major Contributor to Health Disparities in Triple Negative Breast Cancer. Cancers, 13(24), 6262. https://doi.org/10.3390/cancers13246262

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop