Abstract
Up to 30% of patients with metastatic breast cancer eventually develop brain metastasis, yet the pathologic mechanism behind this development remains poorly understood. Here, we profiled long noncoding RNAs in brain metastatic tumors from patients with breast cancer and found that the X-inactive–specific transcript (XIST) was significantly downregulated in these tissues. XIST expression levels inversely correlated with brain metastasis, but not with bone metastasis in patients. Silencing of XIST preferentially promoted brain metastatic growth of XISThigh cells in our xenograft models. Moreover, knockout of XIST in mice mammary glands accelerated primary tumor growth as well as metastases in the brain. Decreased expression of XIST stimulated epithelial–mesenchymal transition and activated c-Met via MSN-mediated protein stabilization, which resulted in the promotion of stemness in the tumor cells. Loss of XIST also augmented secretion of exosomal miRNA-503, which triggered M1–M2 polarization of microglia. This M1–M2 conversion upregulated immune suppressive cytokines in microglia that suppressed T-cell proliferation. Furthermore, we screened an FDA-approved drug library and identified fludarabine as a synthetic lethal drug for XISTlow breast tumor cells and found that fludarabine blocked brain metastasis in our animal model. Our results indicate that XIST plays a critical role in brain metastasis in breast cancer by affecting both tumor cells and the tumor microenvironment and that the XIST-mediated pathway may serve as an effective target for treating brain metastasis.
Significance: These findings describe mechanisms of how loss of the lncRNA XIST promotes brain metastasis in breast cancer and identify fludarabine as a potential therapeutic agent that specifically eliminates XISTlow tumor cells in the brain. Cancer Res; 78(15); 4316–30. ©2018 AACR.
Introduction
Around 30% of patients with metastatic breast cancer will eventually develop brain metastasis, which profoundly affects the cognitive and sensory function as well as morbidity of the patients (1). The development of brain metastasis is a complex, multistep process, including invasion of tumor cells through the blood–brain barrier (BBB) to reach the brain parenchyma, adaptation to the brain microenvironment to acquire growth factors, and signaling that are critical for their survival and sustained growth. Adaptation is mainly achieved by reprograming of the cells within the brain metastatic niche that consists of astrocytes, microglia, and various immune cells (2). Therefore, only a small population of disseminated tumor cells, that have the appropriate genetic profile, will eventually outgrow in the brain (3, 4). Although many genes such as COX2, Cx43, and CTSS have been shown to play roles in brain metastasis by affecting some of these steps, it is not clear whether there exists a single master regulatory oncogene or tumor suppressor whose expression controls the entire process of brain metastasis (3, 5). Our group has previously shown the importance of noncoding RNAs, especially miRNAs in breast cancer brain metastasis, through their targeting of multiple oncogenes that regulate the invasive ability of cancer cells, as well as cancer stem cell–like populations (6, 7). Another group of noncoding RNAs called long noncoding RNAs (lncRNA) have recently drawn much attention due to their potential roles in tumor progression (8). LncRNAs are more complex in structure, and their mechanisms of gene modulation include chromatin remolding, as well as transcriptional and posttranscriptional regulation (9). Therefore, it is likely that this group of genes also play critical roles in the key steps of brain metastasis, including stem cell growth and reprograming tumor microenvironment.
Microglia are the major innate immune cells found in the brain that become activated under many pathologic conditions including infection, injury, and cancer of the central nervous system (10). Activated microglia are abundant source of inflammatory molecules that can affect the development of neurodegenerative diseases as well as tumor progression (11). Therefore, microglia are key components of the tumor microenvironment in brain metastasis. Similar to macrophages, activated microglia have both tumor-suppressive (M1) and tumor-promoting (M2) roles depending on the activation of specific signal pathways (12). However, it is not yet clear how metastatic tumor cells evade the cytotoxic effect of M1 microglia, while at the same time inducing the M1–M2 phenotypic change that supports their growth. In addition to microglia, tumor-infiltrating lymphocytes (TIL), which are known to be able to efficiently eliminate tumor cells by triggering a series of antitumor responses, are frequently found at tumor sites in central nervous system cancer and metastatic brain tumors, even though a normal brain is considered to be an immune privileged site (13). However, how tumor cells evade the antitumor effect of TILs is yet to be clarified. In this study, we show that the loss of lncRNA X-inactive specific transcript (XIST) promotes breast cancer brain metastasis by enhancing both stemness and aggressiveness of tumor cells through induction of epithelial-mesenchymal transition (EMT)– and MSN-mediated upregulation of c-Met. We also show that loss of XIST in tumor cells causes local immune suppression by converting the microglia to the M2 phenotype through the transport of exosomal miR-503 from the tumor cells.
Materials and Methods
Cell culture and reagents
Human breast carcinoma cell lines, MCF7, ZR75-1, SKBR3, and MDA-MB-231 (MDA-231), were purchased from ATCC. SIM-A9 was purchased from Kumi Nagamoto-Combs, through Kerafast.com. MDA-MB-231BrM2a (231BrM) was a kind gift from Dr. Massague (Memorial Sloan-Kettering Cancer Center, New York, NY; ref. 4). SKBrM3 and MDA-MB-231BrM2a are brain metastatic cell lines derived from parental SKBR3 and MDA-MB-231 cells through three rounds of in vivo selections (14). SKBR3, SKBrM3, MCF7, MDA231, and 231BrM were cultured in DMEM supplemented with 10% FBS, streptomycin (100 mg/mL) and penicillin (100 U/mL). ZR75-1 was cultured in RPMI medium with 10%FBS. SIM-A9 was cultured in DMEM/F12 medium with 5% FBS. T cells were isolated from mouse spleen by using the pan T-cell isolation kit (Miltenyi Biotec) and cultured in RPMI with 10% FBS. All cells were grown at 37°C in a 5% CO2 atmosphere. All cell lines were obtained between 2010 and 2016, and they were authenticated by qRT-PCR analysis for the expression of 20 signature genes.
Plasmids and reagents
Lentiviral vectors expressing shRNA for XIST and MSN were obtained from abm Inc, and Dharmacon, respectively. siRNA targeting MPP1, MSN, and MID1 were purchased from Dharmacon. The plasmids expressing dCAS9 (lenti dCAS-VP64_Blast) and MS2 (lenti MS2-P65-HSF1_Hygro), and gRNA backbone plasmid (lenti sgRNA (MS2)_zeo backbone) were obtained from Addgene. The gRNAs targeting XIST promoter were cloned into the backbone plasmid as described previously (15). gRNA-11 and gRNA-8 target 5′TGTCCGGCTTTCAATCTTCT3′ and 5′GCAGCGCTTTAAGAACTGAA3′, respectively. DNA oligomers were synthesized by Invitrogen. The plasmid-expressing miR503 was purchased from System Biology Inc.
Quantitative real-time PCR
Total RNAs were isolated from cells and reverse transcribed as described previously (16). The cDNA was then amplified with a pair of forward and reverse primers for the following genes:
XIST-F 5′ GTGATCAGCACCCCAGCTAT3′
XIST-R 5′ TCTTTCTTTTCCTCCCAGCA3′
c-Met-F 5′ GGTCCTTTGGCGTCGTCCTC3′
c-Met-R 5′ CTCATCATCAGCGTTATCTTC3′
MSN-F 5′ ACCGGGAAGCAGCTATTTGA3′
MSN-R 5′ GAACTTGGCACGGAACTTAA3′
MID1-F 5′ CCTGTCAACATGTTGAAGTC3′
MID1-R 5′ GCAATCTGCTGAGCCAGTTT3′
MPP1-F 5′ ATTGAATACTGTGACCGAGG3′
MPP1-R 5′ TTCTAGGATCTCATCCCCCA3′
E-Cad-F 5′ GACCAGGACTATGACTACTTGAAACG3′
E-Cad-R 5′ ATCTGCAAGGTGCTGGGTGAACCTT3′
Vim-F 5′ ATTCCACTTTGCGTTCAAGG3′
Vim-R 5′ CTTCAGAGAGAGGAAGCCGA3′
mT-antigen-F 5′ GGAAGCAAGTACTTCACAAGGG3′
mT-antigen-R 5′ GGAAAGTCACTAGGAGCAGGG3′
Mouse-CD86-F 5′TTACGGAAGCACCCACGATG3′
Mouse-CD86-R 5′ CCTGTTACATTCTGAGCCAGT3′
Mouse-CD204-F 5′TTCAATGACAGCATCCCTTCC3′
Mouse-CD204-R 5′ GCTTTCGATTCTCTCCTCCAT3′
Mouse-Arg1-F 5′ TCCTTAGAGATTATCGGAGCG3′
Mouse-Arg1-R 5′ GTCTTTGGCAGATATGCAGG3′
Mouse-CD74-F 5′ AGAGCCAGAAAGGTGCAGC3′
Mouse-CD74-R 5′ GATGCATCACATGGTCCTGG3′
Mouse-CD80-F 5′ TACACCACTCCTCAAGTTTCC3′
Mouse-CD80-R 5′ CAGGTAATCCTTTTAGTGTCTG3′
Mouse-Arg2-F 5′ CAAATTCCTTGCGTCCTGACG3′
Mouse-Arg2-R 5′ GGTACCTATTGCCAGGCTGT3′
Mouse-IL10-F 5′ AGCAGAGTGAAGACTTTCTTTC3′
Mouse-IL10-R 5′ CCTTGCTCTTGTTTTCACAGG3′
PCR reactions were performed using CFX Connect (Bio-Rad) and iTaq Universal SYBR Green Supermix (Bio-Rad). The thermal cycling conditions composed of an initial denaturation step at 95°C for 5 minutes, followed by 40 cycles of PCR using the following profile: 94°C, 30 seconds; 60°C, 30 seconds; and 72°C, 30 seconds. TaqMan real-time PCR reagent and primer were purchased from Invitrogen. TaqMan PCR was performed as described by the manufacturer.
Human tissue samples
Human primary breast cancer specimens were obtained from surgical pathology archives of the, Wake Forest Baptist Comprehensive Cancer Center (WFBCCC), and Cooperative Human Tissue Network (CHTN). Human breast cancer brain metastasis samples were obtained from CHTN and Tumor Tissue and Pathology Shared Resource at WFBCCC. All tissue sections were obtained by surgical resection.
Immunocytochemistry and immunofluorescence
Cells or frozen section tissues were fixed with methanol, washed with PBS, and blocked with 2% BSA for 1 hour. The sections were then washed with PBS and incubated with anti-c-Met (1/500: R&D Systems), anti-E-Cad (1/500, Cell Signaling Technology), anti-vimentin (1/500; Cell Signaling Technology), anti-mT-antigen-FITC (1/200; Novus Biologicals), anti-CK14 (1:200, Abcam), anti-CK8/18 (1:200, Abcam), anti-IBA-1 (1:200, Abcam), and anti-CD163 (1/500; Cell Signaling Technology) antibodies for 12 hours at 4°C. Samples were then incubated with anti-rabbit IgG Alexa Fluor 555 molecular probe (Cell Signaling Technology) for 1 hour at room temperature if primary antibody is not conjugated with fluorescence. Fluorescence images were taken by a fluorescent microscope (Olympus IX71).
Gene expression microarray profiling
RNAs were extracted, labeled, and hybridized to Human Gene 1.0 ST array (Affymetrix) using the manufacturer's protocol. Normalization of the data was performed using the RMA algorithm. Existing breast cancer cohort data were analyzed as described previously. For the cancer cohort data analysis, we compiled a microarray dataset of 710 patients from GEO (accession numbers: GSE12276, GSE2034, GSE2603, GSE5327, and GSE14020) as we published previously (14). These datasets were all normalized using MAS5.0, and each microarray was centered to the median of all probes. For each patient, brain metastasis–free survival was defined as the time interval between the surgery and the diagnosis of metastasis.
Gene set enrichment analysis
The Gene MatriX file (.gmx) was generated as described previously (14) by combining the top 200 downregulated genes when MCF7 cells were treated with pathway-specific inhibitors (GSE31912) and also by adopting oncogenic signatures that were deposited in Molecular Signatures Database (MsigDB). Gene Cluster Text file (.gct) was generated from The Cancer Genome Atlas (TCGA) mRNA data of patients with invasive breast cancer. Top 300 patients with highest XIST expression were defined as XIST-high patients, and top 300 patients with lowest XIST expression were defined as XISTlow patients. The number of permutations was set to 1,000.
Animal experiments
All animal experiments were done in accordance with a protocol approved by the Wake Forest Institutional Animal Care and Use Committee. For experimental metastasis assay, nude mice (7–8 weeks) were injected with luciferase-labeled tumor cells in PBS into left cardiac ventricle in a total volume of 100 μL. To confirm a successful injection, the photon flux from whole body of the mice was immediately measured using IVIS Xenogen bioimager (Caliper). The brain metastasis progression was monitored and the luminescence was quantified at the indicated time points. At the endpoint of the study, whole brain was removed, incubated in RPMI1640 medium with 0.6 mg/mL luciferin for 5 minutes, and photon flux was measured. Intracranial injection was performed as previously described. In brief, luciferase-labeled cancer cells in PBS were injected into mouse brain. Photon flux of mice was measured using IVIS Xenogeny bioimager weekly. XISTlox/lox mouse strain is a kind gift of Dr. Rudolf Jaenisch (Whitehead Institute, Cambridge, MA; ref. 17) and was obtained from the MMRRC repository. Tg(MMTV-cre)4Mam/J (MMTV-Cre/Line D) mouse strain was purchased from Jackson laboratory. B6.FVB-Tg(MMTV-PyVT)634Mul/LellJ (MMTV-PyMT) mouse strain was purchased from Jackson laboratory. The MMTV-PyMT and XISTlox/lox were back crossed to B6129F1 for five generations before cross breeding. XISTlox/lox and MMTV-Cre were crossbred to generate MMTV-CRE/XISTlox/lox. MMTV-CRE/XISTlox/lox was further crossbred with MMTV-PyMT to generate XISTlox/lox/MMTV-PyMT or MMTV-CRE/XISTlox/lox/MMTV-PyMT. Siblings were used as a control for each experiment. Mammary gland whole mount slides were generated as described previously (18). The early stage (4 weeks) hyperplastic nodules were counted under microscope. Tumor areas at age of 8 weeks were measured by ImageJ. Tumor length, height, and width were measured by caliper, and tumor size was calculated using the formula, 1/2 × L × W × H. The degree of lung metastases was measured by counting the number of metastatic lesions on the lungs at the endpoint (24 weeks). Micrometastases in the brain were defined as mT-antigen–positive clusters, and the number of lesions was counted under fluorescent microscope.
Exosome isolation
Exosomes were isolated by sequential ultracentrifugation (19). Briefly, CM harvested from cell culture was centrifuged at 300 × g for 10 minutes to remove cells. Exosomes were collected by centrifugation at 2,000 × g for 20 minutes. The supernatant was again centrifuged at 16,500 × g for 20 minutes to remove microvesicles. The supernatant was then passed through a 0.2-μm filter (Sarstedt) to remove particles larger than 200 nm. Exosomes were then collected by ultracentrifugation at 120,000 × g for 70 minutes. Pellets were washed in PBS and centrifuged 70 minutes at 120,000 × g for 70 minutes and finally suspended in PBS.
Plasma exosomes were isolated as described previously (19). In brief, 1 mL of serum was loaded onto the qEV size exclusion columns (Izon), and exosomes were elution with PBS. Isolated exosomes were concentrated by centrifugation using the Amicon Ultra-4 10 kDa column.
The size distribution of isolated exosomes was analyzed by nanoparticle tracking analysis using Nanosight NS50 and by electron microscopy.
Exosome uptake analysis
Exosomes were isolated from cells infected with lentivirus expressing PalmGFP (a kind gift from Dr. Breakefield; ref. 20). Human microglial cells were stained with PKH26 according to product manual and seeded onto cell culture slides. Cells were treated with 50 μg/mL exosomes for 24 hours and they were fixed and stained with DAPI (Invitrogen). The exosome uptake was observed by Olympus FV1200 spectral laser scanning confocal microscope.
Carboxyfluorescein succinimidyl ester assay
T cells were isolated from the spleen of BALB/c mouse by using Pan T Cell Isolation Kit II (Miltenyi Biotec), followed by incubation with Dynabeads Mouse T-Activator (Thermo Fisher) and 2 U/mL IL2 for 2 days. A total of 106 activated T cells were labeled by CellTrace carboxyfluorescein succinimidyl ester (CFSE) dye (Thermo Fisher) and cocultured with same amount of microglia for 5 days, followed by FACS analysis.
Analysis of H3K4me3 chromatin immunoprecipitation sequencing data
We acquired the publicly available chromatin immunoprecipitation sequencing (ChIP-seq) data from Sequence Read Archive (SRA) of NCBI (https://www.ncbi.nlm.nih.gov/sra). The H3K4me3 ChIP-seq data were derived from the following projects: PRJNA266312 (21), PRJNA142885 (22), and PRJNA363069. The raw fastq files were derived from SRA, by the fastq-dump tool. We checked the data quality of the fastq files by using Fastqc software (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and identified all of the fastq files passed the quality control. Then we aligned the raw fastq files to hg19 genome by bowtie2. The peak calling was enabled by MACS2 software using default parameters. The visualization is enabled by IGV software from Broad Institute (http://software.broadinstitute.org/software/igv/).
Statistical analysis
All analysis was calculated by GraphPad Prism. Data are presented as mean ±SD. The P value was calculated by an unpaired Student t test. Gene correlation analysis was calculated by liner regression analysis. The Kaplan–Meier survival analysis was calculated by log-rank (Mantel–Cox) test. Significance between each groups were represented as *, P < 0.05; **, P < 0.01.
Results
XIST expression is downregulated in breast cancer brain metastasis
To identify key lncRNAs in brain metastasis of breast cancer, we isolated RNA from microdissected primary tumors and brain metastatic tumors (Fig. 1A). The expression profiles of lncRNAs were examined in pooled samples using a quantitative real-time PCR (qRT-PCR)-based lncRNA array that we recently developed. The cut-off threshold was set to over 10-fold difference between the two groups, followed by validation of selected lncRNAs in individual sample. We found that XIST was the most highly differentially expressed lncRNA between the two groups (Fig. 1B; Supplementary Fig. S1A and S1B). XIST gene is located on the X chromosome and is known to be involved in X chromosome inactivation by suppressing the genes on this chromosome in a cis-acting manner (23). We also examined the expression of XIST in two paired brain tropic tumor cell lines and found that XIST expression was also significantly lower in these cells compared with their parental cells (Fig. 1C). The results of Kaplan–Meier analyses of 710 patients with breast cancer using combined cohort data (14) indicated that low levels of XIST was associated with a poor brain metastasis-free survival, but not bone metastasis–free survival (Fig. 1D; Supplementary Fig. S1C). Furthermore, XIST expression was found significantly downregulated in brain metastatic lesions compared with other metastatic tumors based on an organ-specific cohort analysis (Fig. 1E). We also examined XIST expression in various subtypes of breast cancer using the TCGA database and found that XIST expression levels were significantly decreased in basal-type breast tumors (Fig. 1F), which are the most common types of breast cancer that known to metastasize to the brain. Interestingly, our analysis of TCGA ChIP-seq database for various breast cancer cell lines revealed that XIST expression was strongly correlated to the histone modification of H3K4me3 at the promoter region of XIST gene (Supplementary Fig. S1D). These results strongly indicate that the loss of XIST expression plays a key role in brain metastasis.
XIST preferentially suppresses brain metastasis in vivo
To address the question how downregulation of XIST is involved in brain metastasis, we first established cell lines expressing short hairpin RNAs (shRNA) that targeted four different locations in the XIST gene. The knockdown effects of shXIST in MCF7 cells, which have high endogenous levels of XIST were confirmed by qRT-PCR (Fig. 2A). To examine the effects of XIST on breast cancer metastasis in vivo, we inoculated luciferase-labeled MCF7 cells with or without knocking down of XIST into nude mice by intracardiac injection, followed by monitoring the tumor growth by measuring the amount of bioluminescence (BLI). We found that the knockdown of XIST in both cell lines significantly promoted their brain metastatic abilities (Fig. 2B–D). This metastasis-suppressive function of XIST was also verified with another two XIST high breast cancer cell lines, SKBR3 and ZR75-1 in vivo (Fig. 2E–H; Supplementary Fig. S2A–S2C). To further examine whether gain of function of XIST suppresses metastasis, we used the dCas9 system (15) to activate and restored XIST expression in 231BrM cells using gRNAs that target key transcriptional activation regions. We have designed 20 gRNAs based on the prediction by the Benchling website and successfully restored the expression of XIST by gRNA dXIST#11 (Fig. 2I). We found that the overexpression of XIST in 231BrM cell significantly attenuated its metastatic ability to the brain but not to the bone (Fig. 2J–L; Supplementary Fig. S2D and S2E). These in vivo results suggest that knockdown and over expression of XIST preferentially promotes and suppresses brain metastasis, respectively (Supplementary Fig. S2F). Notably, knockdown of XIST in MCF7 and ZR75-1 significantly promoted their growth in vivo (Fig. 2M). We also observed a high spontaneous brain metastasis incidence in shXIST groups, suggesting that knockdown of XIST enhances tumor invasion and extravasation at the primary site (Fig. 2N; Supplementary Fig. S2G).To further test the ability of XIST-knockdown cells to colonize the brain, we transplanted MCF7-shCTRL or MCF7-shXIST cells into nude mice by intracranial injection, followed by monitoring tumor growth by bioluminescence imaging. Knockdown of XIST in MCF7 cells significantly promoted tumor growth in the brain (Fig. 2O and P). These results suggest that XIST contributes to multiple metastatic events that are crucial for colonization and outgrowth of tumor cells preferentially in the brain.
Mammary-specific knockout of XIST enhances brain metastasis
To gain a further insight into the role of XIST in tumorigenesis and metastasis in immunocompetent animals, we generated a mammary-specific XIST knockout mouse by crossbreeding XISTlox/lox and MMTV-Cre mouse. However, the knockout of XIST was insufficient to induce tumorigenesis even at 12 months of age. Therefore, we further crossbred the hybrid mouse with the MMTV-PyMT mouse, which is known to spontaneously develop lung metastasis but not brain metastasis (24). The resultant PyMT-XISTΔ/Δ mouse developed mammary tumors at a significantly accelerated rate compared with the PyMT-XISTWT mouse. As early as 3 weeks after birth when the mammary gland starts to develop, hyperplastic nodules appeared in the mice of both genotypes. However, the XIST knockout mice had a significantly higher number of nodules compared with the control mice, suggesting an early onset of tumor development and a higher frequency of tumor initiation in the XIST knockout mouse (Fig. 3A and B). At 8 weeks, by which time the mammary glands have developed and the end buds have started to branch, XIST knockout mice developed significantly larger tumors compared with control mice in their respective mammary glands (Supplementary Fig. S3A and S3B). At 20 weeks, the volumes of palpable tumors in XIST knockout mice were significantly larger than that of control mice (Fig. 3C). The MMTV-PyMT mouse model represents luminal breast cancer, and we have shown that decreased XIST expression was more prominent in patients with triple-negative breast cancer with a basal-like phenotype. Therefore, we speculated that the loss of XIST promotes the phenotypic switch from luminal-like to basal-like tumor in MMTV-PyMT model. Interestingly, we found a mixture of luminal- and basal-type tumor cell populations in the XIST knockout tumors stained by anti-CK18 (luminal) and anti-CK-14 (basal) antibodies, while the majority of tumors in control mice were CK18+ luminal type (Fig. 3D). Importantly, 69% of the XIST knockout mice developed micrometastasis in the brain at 24 weeks, which could be detected by anti-mT-antigen antibody and majority of them were CK14 positive (Fig. 3E–G). This observation was further verified by qRT-PCR analysis of mouse brains using mT antigen–specific primers (Fig. 3H). We also established two syngenic cell lines isolated from XIST WT (Py216) and XISTΔ/Δ (Py222) tumors and verified the XIST expression by qRT-PCR (Supplementary Fig. S3C). We then tested the metastatic ability of these cells in wild-type mice by intracardiac injection and found that the knockout of XIST preferentially promoted brain metastasis but not bone metastasis (Fig. 3I and J; Supplementary Fig. S3D). Hematoxylin and eosin analysis also demonstrated a clear perivascular growth pattern of tumor cells in the brain metastatic lesions, which was consistent with other established brain metastatic models (Fig. 3K). Our data strongly suggest that the knockout of XIST drastically promoted brain metastasis in both MMTV-PyMT transgenic and syngenic models that were consistent with the xenograft data. These results further support our notion that the loss of XIST plays a critical role in brain metastasis.
Knockdown of XIST promotes EMT and stemness
Our in vivo data demonstrated that loss of XIST promotes malignancy of tumor cells, which enable them metastasize to the brain. However, the molecular mechanism underneath is not clear. During the establishment of XIST knockdown cell lines from MCF7 and ZR75-1, we observed a morphologic change of those cells with the acquisition of a mesenchymal phenotype (Fig. 4A), suggesting that XIST may play a role in the EMT. Tumor cells are known to gain invasiveness, migratory ability, and stemness following the EMT (25). We therefore examined the expression of EMT markers on these cells and found that E-cadherin and vimentin were significantly suppressed and upregulated, respectively, when XIST was knocked down (Fig. 4B–D; Supplementary Fig. S4A and S4B). We also observed the same results in Py216 and Py222 cell lines (Supplementary Fig. S4C and S4D). Because the EMT is closely associated with the properties of cancer stem-like cells that are known to be involved in tumor progression, we performed mammosphere assays to test the self-renewal ability of cells with or without knockdown of XIST. As shown in Fig. 4E; Supplementary Fig. S4E, knockdown of XIST significantly increased their sphere-forming ability. In addition, an in vivo limiting dilution assay indicated that knockdown of XIST significantly enhanced the tumor-initiating ability of these cells (Fig. 4F). Furthermore, because cancer stem cells have characteristic of intrinsic therapy resistance, we performed clonogenic assays in cells with or without knockdown of XIST, followed by radiation, and found that knockdown of XIST significantly increased their survival after radiation (Supplementary Fig. S4F). Knockdown of XIST also promoted cell invasion and migration, as shown using a transmigration assay that mimics the BBB (Fig. 4G).These results strongly indicate that a loss of XIST enhances the EMT and tumor stemness, and hence, promotes metastatic progression.
Knockdown of XIST activates the c-Met pathway by upregulating MSN
To identify key downstream effectors of XIST in patients with breast cancer that affect brain metastasis, we performed a pathway screening based on gene set enrichment analysis (GSEA) by using the TCGA database of 937 patients and 119 breast-specific oncogenic signatures established previously (14). This approach identifies up and downregulated pathways based on the enrichment of signature genes in patients with two phenotypes, for example, patients with high or low expression of XIST. Therefore, the result of the screening is not dependent on the expression of a single gene and rather represents overall activation of a pathway, which is a consequence of both transcriptional and posttranscriptional regulation. Out of those pathways tested, four of them, c-Met, EIF4E, Myc, and CCND1 pathways, were significantly enriched in patients with low XIST levels (Fig. 5A; Supplementary Fig. S5A). Among these four pathways, the c-Met pathway showed the highest enrichment score (NES = 2.1), and it was the only pathway significantly correlated with brain metastasis–free survival in patients with breast cancer (P = 0.004, Supplementary Fig. S5A). Western blot analysis indicated that there was a strong upregulation of total c-Met in three pairs of XIST knockdown cell lines (Fig. 5B). Furthermore, knockdown of XIST strongly increased c-Met phosphorylation even without HGF stimulation, indicating that c-Met signaling is constitutively activated in XISTlow cells (Supplementary Fig. S5B). We also found a negative correlation between the RNA levels of XIST and the protein expression of c-Met in brain metastasis samples from patients (Fig. 5C). These data suggest that loss of XIST upregulates c-Met and its downstream signaling pathway. To understand how XIST regulates c-Met, we performed a gene expression array analysis on MCF7 cells with or without knockdown of XIST and found that 168 genes were significantly upregulated in MCF7-shXIST cells [>2-fold change, P < 0.05, false discovery rate (FDR) < 0.05; Fig. 5D]. Interestingly, several EMT drivers and mesenchymal markers were among these genes (Supplementary Fig. S5C), which is consistent with our finding in Fig. 4. Notably, out of the 168 genes, 8 genes were found to be significantly correlated with poor brain metastasis–free survival (Fig. 5D). Because XIST functions only in a cis-acting manner on the X chromosome, from among the 8 genes, we chose to further examine MID1, MPP1, and MSN, which are located on the X- chromosome (Supplementary Fig. S5D). We then confirmed upregulation of these genes after knocking down of XIST in MCF7 (Fig. 5E). Among those genes, MSN has been previously reported to indirectly bind c-Met through CD44, and mediates c-Met internalization (26). Upon internalization, c-Met could sort into early endosomes for degradation or traffic to the nondegradative perinuclear compartment (27). Therefore, we explored the possibility that XIST-mediated MSN expression contributes to c-Met upregulation by mediating its accumulation in the perinuclear compartment. Notably, c-Met expression was drastically decreased by knockdown of MSN but not by MPP1 or MID1in MCF7-shXIST cells (Fig. 5F; Supplementary Fig. S5E). In addition, gene correlation analyses in breast cancer cell lines demonstrated that MSN is the only gene whose expression is inversely correlated with XIST in triple-negative tumor cell lines, suggesting a distinct role of MSN in triple-negative tumors, which are known to express high amounts of c-Met (Supplementary Fig. S5F). Our immunocytochemical analysis indicated that c-Met was localized to the perinuclear region in MCF7-shXIST cells (Fig. 5G). In addition, knockdown of MSN also suppressed c-Met expression and reduced its accumulation in the perinuclear compartment in 231BrM cells (Fig. 5H and I). These results indicate that down regulation of XIST enhances c-Met expression by upregulating MSN, which then facilitates perinuclear localization and stabilization of the c-Met protein. We also generated SKBrM3 cell subline that overexpressed XIST using the dCas9 approach and found that the expression of both MSN and c-Met were strongly suppressed in the XIST-overexpressing cells (Fig. 5J; Supplementary Fig. S5G). Moreover, the restoration of XIST expression suppressed c-Met expression by diminishing the perinuclear accumulation of the c-Met protein (Fig. 5K). To further validate the notion that c-Met is a crucial downstream of XIST to enhance brain metastasis, we suppressed c-Met in MCF7-shXIST, 231BrM, and SKBrM3 cells by doxycycline-inducible c-Met shRNA. Loss of c-Met in these cells significantly decreases their transmigration abilities (Fig. 5L and M; Supplementary Fig. S5H and S5I). We also found that knockdown of c-Met in MCF7-shXIST cells significantly decreased their metastatic abilities on brain but not on bone in vivo (Fig. 5N-P). IHC analysis also showed strong c-Met expression in tumors from XIST-knockout mice (Supplementary Fig. S5J).These results further support our notion that induction of c-Met in XISTlow cells is mediated by MSN.
XISTlow cells secrete exosomal miR-503 to reprogram microglia in the brain
The effect of the loss of XIST on metastasis is more profound in the brain compared with other distant organs in the clinical settings and this selectivity was also evident in our xenograft model. These observations suggest that the tumor microenvironment may affect the colonization ability of cells with low levels of XIST expression in the brain. Because microglia are one of the major components of the brain microenvironment, and they have been previously shown to suppress brain metastasis in lung cancer, it is plausible that XISTlow cells reprogram microglia toward a prometastatic phenotype through cell–cell communication (28). Interestingly, when we examined the relationship between XIST expression and 95 immune-related signaling pathways in MSigDB by GSEA analysis, the innate cell immune signature was highly enriched in XISTlow patients (Fig. 6A). IHC analysis of brain metastatic lesions in patients with breast cancer revealed a strong localization of M2 microglia (both IBA-1 and CD163 positive) around the tumor invasion front (Fig. 6B). Metastatic cells are known to educate the tumor microenvironment by various factors including cytokines and exosomes (29). We found that conditioned medium of MCF7-shXIST cells significantly promotes the M1–M2 conversion of microglia as shown by the expression of M1 gene CD86 and M2 gene ARG1. However, depletion of exosomes by ultracentrifuge significantly attenuated such effect, indicating that exosomes are major contributor to the M1–M2 conversion of microglia (Fig. 6C). Although it is not clear how M2 microglia is generated in the tumor areas by exosomes, recent studies showed that miRNAs carried by exosomes play critical roles in the generation of metastatic niche even before the presence of tumor cells (30). Therefore, we hypothesized that exosomes secreted from tumor cells mediate the generation of M2 microglia in the tumor lesions. To test this hypothesis, we first prepared exosomes from MCF7 and MCF7-shXIST cells, followed by size validation by the nanoparticle tracking analysis and electron microscopy (Supplementary Fig. S6A and S6B). RNAs were extracted from the exosomes and they were subjected to miRNA expression profiling (Fig. 6D). Our results revealed that miR-503 was among the most highly expressed miRNAs in the exosomes from MCF7-shXIST cells, and it ranked as the most highly expressed miRNA encoded by the X chromosome. We further validated the induction of miR-503 in cells with knockdown of XIST (Fig. 6E). As shown in Fig. 6F, exosomes were robustly uptake by microglials at a physiologic relevant concentration (10 mg/mL). We also examined the relationship between miR-503 expression and relapse-free survival in patients with breast cancer using the GSE22220 database and found that a higher expression of miR-503 was negatively correlated with the relapse-free survival (Fig. 6G), suggesting that miR-503 plays a role in metastatic progression. Furthermore, we isolated exosomes from the serum of patients with breast cancer with or without brain metastasis, followed by examine the miR-503 levels. We detected a significantly higher level of miR-503 in the serum-derived exosomes from patients with brain metastasis than in those from other patients, suggesting that exosomal miR-503 is a potential biomarker for brain metastasis (Fig. 6H). To test the effect of exosomal miR-503 on microglia, exosomes were prepared from (i) MCF7 cells, (ii) MCF7-shXIST cells, or (iii) MCF7-shXIST cells expressing anti-miR-503, and then incubated with SIM-A9, a mouse microglia cell line. We examined the expression of the M1 and M2 markers in these treated microglia, and found that the treatment of microglia with exosomes from MCF7-shXIST cells significantly up- and downregulated the M2 and M1 markers, respectively. This effect of the exosomes was rescued by suppressing miR-503 expression in the MCF7-shXIST cells (Fig. 6I). In addition, CM derived from brain tropic tumor cells (231BrM and SKBRM3) significantly promoted the microglial M1–M2 transition (Supplementary Fig. S6C). When miR-503 was ectopically expressed in SIM-A9 cells, M2 markers were significantly upregulated while the expression of M1 markers was suppressed (Fig. 6J). In addition, the ectopic expression of miR-503 suppressed ROS in SIM-A9 cells, which is also indicative of the M1–M2 conversion. (Supplementary Fig. S6D). These results strongly suggest that exosomal miR-503 is capable of converting microglia from the M1 to the M2 phenotype. STAT3 and NFκB are two major signaling pathways that known to play critical roles in the M1–M2 conversion of macrophage (31). We found that miR-503 strongly increased phosphorylation of STAT3 and decreased phosphorylation of the p65 subunit of NFκB, suggesting a systematic reprogramming of M1 and M2 pathways by miR-503 (Fig. 6K). Treatment of microglia with CM from brain tropic breast cancer cells (231BrM and SKBrM3) further validated our hypothesis (Supplementary Fig. S6E). Furthermore, this miR-503–mediated M1–M2 conversion was significantly suppressed by the STAT3 inhibitor, STATTIC (Supplementary Fig. S6F and S6G). These results indicate that exosomal miR-503 promotes the M1–M2 conversion through manipulating STAT3 and NFκB pathways. Similar to M2 macrophages, M2 microglia have been reported to promote progression of brain tumors by releasing cytokines and growth factors (28, 32). Tumor-associated M2 macrophages are also known to suppress T-cell immunity by expressing PD-L1 (33, 34). Indeed, ectopic expression of miR-503 in both human and mouse microglial cells increased the expression of PD-L1 (Fig. 6L). IHC analysis of brain metastatic lesions from patients with breast cancer also revealed that PD-L1 expression was colocalized with microglial cells adjacent to the tumor, but not with astrocytes (Fig. 6M; Supplementary Fig. S6H). Importantly, when mouse T cells were cocultured with SIM-A9–expressing miR-503, the proliferation rate of T cells was significantly suppressed (Fig. 6N). These results strongly suggest that exosomal miR-503 derived from tumor cells promotes M1–M2 conversion of microglia, followed by enhancing their PD-L1 expression to suppress local immunity and thereby enhance tumor growth.
Fludarabine selectively inhibits the growth of XISTlow tumor cells
We have shown that cells with low levels of XIST are highly aggressive and metastasize to the brain. Therefore, identifying a drug that specifically targets these cells would be a promising therapeutic approach to treat brain metastasis. To accomplish this goal, we used a FDA-approved drug library (n = 1250) and a natural compound library (n = 140) to screen compounds that selectively suppress the growth of XIST low cells, followed by MTS assay (Fig. 7A). Seven drugs were picked up by the first round of screening (Supplementary Fig. S7A). Among those seven drugs, we found that fludarabine was the most effective drug that kills MCF7-shXIST cells compared with MCF7-shCTRL (Fig. 7B). In addition, fludarabine was found to be highly selective for cells with low levels of XIST expression as it showed selective inhibition of 231BrM and SKBrM3 cells, compared with SKBR3, MCF7, and ZR75-1 cells, which have high levels of XIST (Fig. 7C). This inhibitory effect was rescued by restoring the expression of XIST by dCas9 in 231BrM cells (Fig. 7D). This selective cytotoxic effect on XISTlow cells was not observed for cisplatin and paclitaxel, two chemotherapy drugs commonly used for treating patients with breast cancer (Supplementary Fig. S7B). It should be noted that the IC50 for fludarabine on XISTlow cells was around 20 nmol/L, which is more than 10-fold lower than the effective dose for leukemic cells, for which fludarabine is currently in clinical use (35). The in vitro toxicity of fludarabine on primary neuronal cells was minimal (Supplementary Fig. S7C). Because fludarabine is known to be BBB permeable (36), we tested the efficacy of this drug in vivo on tumors with low levels of XIST expression by transplanting MCF7-shXIST cells into nude mice by intracardiac injection, followed by administration of fludarabine intraperitoneally (10 mg/kg) every 2 days for 28 days. We found that fludarabine significantly delayed the onset of brain metastasis, and suppressed the growth of tumor cells in the brain (Fig. 7E–G) without notable toxicity (Supplementary Fig. S7D). Although the specific mechanism of action of this drug on cells with low levels of XIST expression is not yet clear, our results suggest a promising utility of this drug for the patients with XISTlow breast cancer cells that spread to the brain.
Discussion
Brain metastasis is a highly complex and selective process due to the unique microenvironment of this specialized organ. In this study, we identified XIST as a key lncRNA, which functions as a metastasis suppressor gene in the brain. Downregulation of XIST activates three distinct pathways, EMT, MSN/c-Met, and release of exosomal miR-503 (Fig. 7H). Together, they promote the expression of metastatic traits in tumor stem cells and trigger the reprogramming of the microenvironment through exosomal-mediated communication. The result of our analysis of the TCGA database indicates that XIST expression is decreased in 48% of basal-like breast cancer, while it is downregulated in 28% of Her2+ and 19% of luminal A subtype cancers. It should be noted that most of the brain metastases of breast cancer belong to basal-like and Her2+ subtypes (37). Importantly, XIST was found to be downregulated in 78% of tumor tissues from patients with brain metastases, compared with 32% in patients with bone metastases and 41% with lung metastasis. These observations are consistent with the results of our study and further support our notion that XIST exerts a metastasis-suppressive function preferentially in the brain.
XIST silences X chromosome genes in a cis-acting manner, mainly through PRC2 complex–mediated H3K27me3 methylation (38). How this lncRNA functions only on the X chromosome in a cis-acting manner is not well understood. However, when XIST is artificially introduced to autosomes, a similar cis-acting effect of chromosome inactivation was observed, indicating that autosomal chromatin is fully competent to be inactivated by XIST (39). Aberrant expression of XIST has been implicated to inactivate miRNAs by a “sponge” effect and also to dysregulate the AKT signal (40). In this study, we showed that both MSN and miR-503 are upregulated following knockdown of XIST. Importantly, both genes are encoded by the X chromosome, and hence, they are under the direct epigenetic control of XIST. Therefore, the XIST axis is considered to be a gender-specific oncogenic pathway.
We previously demonstrated that metastatic growth of breast cancer in the brain is strongly associated with the self-renewal abilities of cancer stem cell (41, 42). Moreover, the link between the EMT and the properties of CSCs has been well established in breast cancer (43). Therefore, the increased stemness of cells with low XIST expression can in part be attributed to the XIST-mediated EMT. Upon knockdown of XIST, we observed a systematic upregulation of multiple transcription factors that target the E-box (ZEB1, ZEB2, and TWIST2), as well as components of the Notch and Wnt pathways (HEY1, TCF4, and WNT5A) that are known to be potent EMT promoters (44). However, the upregulation of EMT genes is not directly mediated either by MSN or miR-503 on the X chromosome (data not shown). Therefore, we speculate that the EMT phenotype is caused by a net effect of XIST downregulation of the X chromosome, which is consistent with the observation that XIST is decreased in many of the basal-like breast cancer samples.
We have shown that the loss of XIST upregulates c-Met expression through activation of MSN on the X chromosome. It is noteworthy that high levels of c-Met signaling in pancreatic cancer and glioblastoma have also been reported to be associated with an increased expression of stem cell–related genes (45). Moreover, we previously showed that c-Met signaling specifically promotes brain metastasis by inducing angiogenesis through CXCL1 and IL8 secretion, suggesting that activation of this pathway also modifies the tumor microenvironment (14). Therefore, c-Met signaling is likely to contribute to brain metastasis at multiple levels that affect the malignancy of tumor cells as well as in the local microenvironment. MSN encodes moesin, which is one of the three components of the ERM complex (ezrin/radixin/moesin). In addition to c-Met internalization, ERM can also serve as a signal transducer that binds to Grb2 and SOS to further amplify c-Met signaling (46). Moreover, perinuclear c-Met has been shown to be able to activate another key oncogenic signaling, STAT3 pathway, by direct phosphorylation and facilitates STAT3 to enter into the nucleus (47).
Successful colonization of brain metastatic cells requires reciprocal communications with microenvironmental cells and local immune cells. We have previously shown that a group of inflammatory cytokines that activate astrocytes and brain endothelial cells are highly secreted by brain metastatic cells (14). In addition to soluble factors, tumor cell–secreted miRNAs in exosomes are considered to be important mediators between the tumor and the host microenvironment (30). Here, we have shown that the loss of XIST enhances the secretion of exosomal miR-503, which is also encoded by the X chromosome, and that exosomal miR-503 is capable of reprogramming microglia from M1 to M2 type by modulating the STAT3 and NFκB pathways. This M1–M2 conversion is reported to occur through activation of signaling pathways by either direct interaction with tumor cells or through soluble factors (48). However, when exosomes were depleted in our M1–M2 conversion assay system, this phenotypic change was drastically diminished, indicating that exosomal miR-503, rather than a soluble factor, plays a major role in the microglial M1–M2 transition in XIST-mediated brain metastasis.
Normal brain is considered to be “immune-privileged,” and normally cytotoxic T cells are unable to penetrate the BBB. However, recent mounting evidence has indicated that both T and B cells are frequently found in brain tumors, including glioma and metastatic tumors (13). The amount of T cells and the level of expression of immune checkpoint genes in the brain tumor areas have been associated with prognosis in patients with brain metastases (13). We have shown here that immune suppression is one of the mechanisms by which M2 microglia promote tumor progression, and this effect is achieved by the upregulation of PD-L1 in the microglia. Although several studies have shown that many growth factors that directly promote the survival and migration of tumor cells are also secreted from M2 microglia, it is likely that both mechanisms contribute to tumor growth in the brain and that these two mechanisms are not mutually exclusive. Nevertheless, our data provide a promising clue of how to exploit immunotherapy for the treatment of brain metastasis. A recent study has shown that there is a synergetic effect between immune checkpoint inhibitors and an oncolytic virus, which increases intratumoral M1-like macrophages in a mouse GBM model (49), suggesting a potential immune therapy to treat brain metastasis could be developed by boosting T-cell activities by pharmacologically manipulating the population of these microglial subtypes in brain metastatic sites.
Several groups, including ours, have shown that there is a strong correlation between the amount of a particular exosomal miRNA in circulation and the prognosis of patients with cancer (50). In this study, we have demonstrated that exosomal miR-503 levels are positively correlated with brain metastasis status. Exosomes secreted from cancer cells can be taken up by nontumor cells in metastatic niches even before the arrival of tumor cells (51). Although the effect of exosomes is likely to be maximum at the local site, it is possible that exosomal miR-503, shed by the primary tumor, can also be taken up by microglia to generate a premetastatic niche at an early stage of tumor progression. Furthermore, exosomal miR-503 present in patient serum may serve as a promising biomarker, with the level of miR-503 in circulation helping to predict the risk of future brain metastasis.
We have shown that fludarabine acts as a synthetic lethal drug in tumor cells with low XIST levels. Importantly, the IC50 of fludarabine against tumor cells having low XIST levels is at least 10 times lower than the clinical dose used in treating leukemia, and furthermore fludarabine showed significant efficacy in our animal model of brain metastasis without any obvious adverse effect. Therefore, it is possible that the synthetic lethality of fludarabine to XISTlow cells, which are mainly triple-negative, is in part due to a synergistic effect with ROS production. Fludarabine is a BBB-permeable compound whose neurotoxicity is tolerable even at a relatively high dose (52) and it is already approved by the FDA. Therefore, fludarabine could be used as a potential therapeutic agent for the treatment of patients with breast cancer with brain metastases that have low levels of XIST expression. In summary, we have shown that the lncRNA, XIST, plays a critical role in brain metastasis in breast cancer in a gender-specific manner by affecting both tumor cells and the tumor microenvironment, and that the XIST-mediated pathway may serve as an effective target for the treatment of brain metastasis.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: F. Xing, Y. Liu, L.J. Metheny-Barlow, K. Watabe
Development of methodology: F. Xing, Y. Liu, K. Wu, K. Watabe
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): F. Xing, S.-Y Wu, Y. Liu, K. Wu, S. Sharma, J. Feng, S. Sanders, R. Singh, P.-A. Vidi, A. Tyagi, H.-W. Lo, K. Watabe
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): F. Xing, Y. Liu, K. Wu, G. Jin, R. Singh, J. Ruiz, R.B. D'Agostino Jr, K. Watabe
Writing, review, and/or revision of the manuscript: F. Xing, Y. Liu, G. Jin, M.D. Chan, J. Ruiz, W. Debinski, B.C. Pasche, L.J. Metheny-Barlow, K. Watabe
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): F. Xing, Y. Liu, Y.-Y. Mo, M.D. Chan, K. Watabe
Study supervision: F. Xing, K. Watabe
Acknowledgments
We thank Dr. Joan Massagué for providing the brain metastatic cell lines. The XISTlox/lox mouse is kind gift from Dr. Rudolf Jaenisch (Whitehead Institute). This work was supported by NIH grant R01CA173499, R01CA185650, and R01CA205067 (to K. Watabe), a pilot grant from Comprehensive Cancer Center of Wake Forest University (to F. Xing). Y. Liu is supported by F31CA200286. The Tumor Tissue and Pathology Shared Resources, and Biostatistics/Bioinformatics Shared Resource are supported by the Comprehensive Cancer Center of Wake Forest University NCI, NIH grant (P30CA012197).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.