Activation of ribosome biogenesis pathway during the EMT transitioning phase in Tri-PyMT cells.

A, Representative flow cytometry plot displays the percentage of RFP+, Double+ (Doub+, EMT transitioning cells), and GFP+ Tri-PyMT cells in culture.

B, Western blot of EMT markers with flow sorted RFP+ (R), Doub+ (D), and GFP+ (G) Tri-PyMT cells. The Doub+ cells exhibit an intermediate EMT status with expression of epithelial marker (E-cadherin and Epcam) expression and higher mesenchymal marker (Vimentin) expression as compared with RFP+ cells. The numbers indicate normalized intensity of the band according to the β-actin band of the sample.

C, The t-SNE plot of scRNA-seq analysis of Tri-PyMT cells. Two major cell clusters with differential EMT status are shown with the AUC value of 10 epithelial marker genes (in red) and 10 mesenchymal marker genes (in green). The marker genes are shown in Supplementary Fig S4A.

D, E, F, Trajectory analysis using Monocle DDR tree. Cell trajectory analysis was performed with filtered EMT-related genes using Monocle 2 model. Five cell states (D) were identified with differential EMT statuses. From epithelial to mesenchymal phenotype, the Cell States were identified in order of 1, 5, 2, 4, and 3. EMT status (E) of cells was also highlighted by the AUC calculation with epithelial and mesenchymal marker genes. EMT Pseudotime (F) was calculated with Cell State 1 (the most epithelial state) as the root.

G, Classification of cells with EMT states. Cells were classified according to their EMT state as Epi, Trans, and Mes; the cell number in each cluster is indicated in the table.

H, The heatmap of differentially expressed genes in Trans cells compared to Epi and Mes cells. Totally, 313 genes were identified with criteria P <0.01, Fold change > 1.2, and Average expression >=5.

I, The common enriched GO_BP pathways with GSEA when comparing Trans vs. Epi and Trans vs. Mes. There are 49 common pathways, including 5 pathways related to Ribosome Biogenesis (RiBi) or rRNA processing. Pathway names are shown in Supplementary Fig S4B.

J, GSEA plots showing the specific enrichment of Ribosome Biogenesis (RiBi) pathway in EMT transitioning (Trans) cells compared to the Epi or Mes cells.

K, The scatter plot displays the correlation of RiBi activity to EMT pseudotime. The polynomial regression line (order = 3) highlights the elevated RiBi pathway in Trans phase cells, R2 = 0.348, P < 0.001.

Activation of RiBi pathway in MET during lung metastasis outgrowth.

A, Schematic of MET induction. GFP+ Tri-PyMT cells were sorted by flow cytometry (left) and injected into mice through the tail vein. Metastasis-bearing lungs were harvested after 4 weeks and analyzed for the gain of epithelial marker (EpCam) in tumor cells (right). Both Epcam+ and Epcam-cells were sorted and submitted for scRNA-seq analysis.

B, The t-SNE plot of scRNA-seq analysis of GFP+ Tri-PyMT cells sorted from metastatic lungs of 3 individual animals. Two major cell clusters with distinct EMT statuses were identified, as indicated by the AUC values of 10 epithelial genes (in red) and 10 mesenchymal genes (in green).

C, D, E, Monocle 2 model-based cell trajectory analysis identifies nine cell states (C). The MET process is emphasized by AUC values of mesenchymal and epithelial markers (D). MET Pseudotime is calculated with State 1 (the most mesenchymal state) as the root (E).

F, GSEA plot exhibits the enrichment of the GOBP_Ribosome Biogenesis (RiBi) pathway in State 8 cells (the most epithelial phenotype) compared to the other states.

G, Scatter plot illustrates the correlation between RiBi activity and MET pseudotime. A polynomial regression line (order = 3) emphasizes the elevated RiBi pathway during the MET process, R2 = 0.342, P < 0.01.

MAPK and mTOR pathways mediated ribosome biogenesis activation in EMT transitioning phase cells.

A, Western blots show phospho-ERK, total ERK, phospho-mTORC1, total mTORC1 and p-rps6 in sorted RFP+, Doub+, and GFP+ Tri-PyMT cells following stimulation with 10% FBS at 0, 10, 30 and 60 mins. Activation of p-ERK, p-mTORC1 and p-Rps6 are quantified by the band intensity after normalizing to loading controls (right).

B, Cell size analysis. Tri-PyMT cells (p5) were analyzed by flow cytometry. The relative cell size was indicated by FSC-A, for individual RFP+, GFP+, and Doub+ cells. Data from three biological replicates.

C, Nascent protein synthesis assay. Histogram of OPP incorporation, as measured by flow cytometry, demonstrates enhanced nascent protein synthesis in Doub+ cells. The rate of new protein synthesis is assessed by adding fluorescently labeled O-propargyl-puromycin (OPP). Three biological replicates, One-way ANOVA, *P < 0.0001, Doub+ vs RFP+, and Doub+ vs GFP+.

D, Fluorescent images reveal increased RiBi activity in Doub+ Tri-PyMT cells, as evidenced by Fibrillarin staining. RFP+, Doub+, and GFP+ Tri-PyMT cells were sorted and stained for nucleoli using an anti-Fibrillarin antibody. The quantification of nucleoli per cell is shown on the right. n = 522 (RFP+), 141 (Doub+), 289 (GFP+). One-way ANOVA, *P < 0.0001.

RiBi inhibition by Pol I inhibitor reduces the EMT/MET capability of tumor cells.

A, Flow cytometry plots display the percentage of GFP+ cells following BMH21 treatment. RFP+/EpCam+ Tri-PyMT cells were sorted and cultured in growth medium with or without BMH21 for 5 days. Percentage of GFP+ cells was analyzed by flow cytometry. Biological repeat, n=2, One-way ANOVA, *P=0.0039, 100nM vs Cont; *P=0.0029, 200nM vs Cont.

B, Western blots reveal the expression of the epithelial marker (E-cad) and mesenchymal markers (Vim and Snail) in Tri-PyMT cells following a 5-day treatment with BMH21 (100nM) and CX5461 (20nM). The numbers indicate normalized intensity of the band according to the β-actin band of the sample.

C, Flow cytometry plots show the percentage of Epcam+/GFP+ Tri-PyMT cells treated with BMH21. GFP+/EpCam-Tri-PyMT cells were sorted and culture in 3D to induce MET for 14 days. The gain of Epcam expression was analyzed by flow cytometry. n=4, One way ANOVA, *P<0.0001, 100nM vs Cont; *P<0.0001, 200nM vs Cont.

RNA Pol I inhibitor synergizes with chemo drug in vitro.

A, Flow cytometry plots display the fluorescence switch of Tri-PyMT cells. Tri-PyMT cells (p5) were treat with cyclophosphamide (CTX, 2μM), Pol I inhibitor (BMH21, 0.1μM) or vehicle control for 7 days. Flow cytometry was performed to analyze the percentage of GFP+ cells. n=3 wells/treatment, One way ANOVA, *P=0.0050, Cont vs. CTX; *P=0.0002, Cont vs. BMH21.

B, Cytotoxic assay of BMH21 and CTX treatments. Tri-PyMT cells were treated with serial concentrations of BMH21 (0, 0.1, 0.2, 0.4, 0.8, and 1.6 μM) in combination with serial concentrations of CTX (0, 0.25, 0.5, 1.0, 2.0, 4.0, 8.0, 16.0 μM). Bars represent the mean value of cell viabilities from duplicated treatments, n=2 wells/treatment.

C, Synergy plots of BMH21 and Cyclophosphamide (CTX). Synergic scores (δ) for each combination were calculated using the ZIP reference model in SynergyFinder 3.0. Deviations between observed and expected responses indicate synergy (red) for positive values and antagonism (green) for negative values.

RNA Pol I inhibition synergizes with chemo drug in treating lung metastatic tumors.

A-C, The same number of RFP+ and GFP+ Tri-PyMT cells were injected in animals via tail vein. Animals were treated CTX (100mg/kg, once a week, i.p.), BMH21 (20mg/kg, 5 times a week, i.p.) or both in combination.

A, Lung metastasis growth curves as measured by bioluminescent imaging (BLI). n=5, two-way ANOVA with Tukey’s test, *P=0.0154, Cont vs. CTX+BMH; *P=0.0151, BMH21 vs. CTX+BMH21; *P=0.0372, CTX vs. CTX+BMH21 at Day 21 post-inoculation.

B, Fluorescent images display RFP+ and GFP+ metastatic nodules in the lungs treated with CTX and BMH21. Cell nuclei were stained with DAPI. Scale bar: 200μm.

C, Fluorescent images exhibit EMT statuses of RFP+ and GFP+ metastases. Sections were stained with anti-EpCam antibody; cell nuclei were stained with DAPI. Scale bar: 50μm.

D, Lung metastasis growth curves for LM2 model. Lung metastatic breast tumor cells (LM2 cells) were injected into animals via tail vein. Animals were treated CTX (100mg/kg, once a week, i.p.), BMH21 (20mg/kg, 5 times a week, i.p.) or both in combination. n=5, two-way ANOVA with Tukey’s test, *P=0.0404, Cont vs. CTX+BMH21; *P=0.0161, BMH21 vs. CTX+BMH21; *P=0.0187, CTX vs. CTX+BMH at Day 32 post-inoculation. Representative BLI images of Day 32 were presented on the right.

Characterization of fluorescent marker switch in Tri-PyMT cells.

A, RFP+ Tri-PyMT cells were sorted via flow cytometry and cultured in DMEM with 10% FBS. Plot shows the percentage of RFP+, GFP+ and Doub+ cells at different generations of cell culture. The experiment was repeated twice with similar results.

B, Cell cycle analysis of Tri-PyMT cells. Tri-PyMT cells (p7) composed by both RFP+ and GFP+ cells were cultured in growth medium. Cell cycle analysis was performed by using EdU incorporation and Click-chemistry. Cells in G0/G1, S and G2/M phases were quantified by flow cytometry. n =4, two-way Anova, P=0.0242, for S phase RFP+ vs GFP+; P < 0.0001, for S phases Doub+ vs RFP+ and Doub+ vs GFP+ cells.

RNA-sequencing analysis of Tri-PyMT cells.

RFP+ Doub+ and GFP+ Tri-PyMT cells were purified using flow cytometry sorting, total RNA was extracted and submitted for RNA-sequencing analysis.

A, Heatmap of selected EMT marker genes in RFP+, Doub+ and GFP+ cells. Of note, the intermediate expression of both epithelial and mesenchymal markers in Doub+ cells.

B, Heatmap of differentially expressed genes in Doub+ cells compared to RFP+ and GFP+ cells. In total, 213 genes are presented in 4 clusters, P<0.01, Doub+ vs RFP+ and GFP+.

C, Gene set over-representative assay with upregulated genes in Doub+ cells showing the significant enrichment of Ribosome pathway (326 KEGG pathways analyzed in total).

The scRNA-sequencing analysis of Tri-PyMT cells.

RFP+, Doub+, and GFP+ Tri-PyMT cells were sorted via flow cytometry and submitted for sc-RNA sequencing analysis.

A. Quality controls of the scRNA-seq. Plots show the single-cell QA/QC analyses including total counts of UMI, detected gene features, and the percentage of mitochondria genes. Totally, 6167 cells out of 7214 cells were filtered according to the criteria for downstream analyses.

B. C. EMT trajectory analysis plots. Cell trajectory analysis was performed with filtered EMT-related genes using Monocle 2 model. Pairs of epithelial and mesenchymal marker genes, S100a4 vs Epcam (B) and Vim vs Krt18 (C), are used to highlight the EMT status of individual cells.

D. Scatter plot shows the expression of Ribosome biogenesis genes (AUC values of RiBi pathway genes) according to EMT cell states. A trend of highest RiBi activation was detected in State 5 cells (p=0.072 for 5 vs 1, and p<0.01 for 5 vs 2, 3, 4, One-way ANOVA).

RT-PCR analysis of RiBi genes in Tri-PyMT cells.

RFP+, Doub+ and GFP+ cells were purified from Tri-PyMT (p7) culture using flow cytometry sorting. Total RNA was extracted and analyzed by RT-PCR. The heatmap shows z-score of RiBi related genes, epithelial and mesenchymal marker genes expression with Gapdh as internal control, n=3.

The scRNA-sequencing analysis of GFP+ Tri-PyMT cells.

GFP+ Tri-PyMT cells, including Epcam+ and Epcam-cells were sorted from metastasis-bearing lungs and submitted for single cells RNA sequencing analysis.

A. Quality controls of the scRNA-seq. Plots show the single-cell QA/QC analyses of the total counts of UMI, detected gene features and the percentage of mitochondria genes. Totally, 4961 cells out of 6934 cells were included in the downstream analyses.

B. Plots of MET trajectory analysis. Cell trajectory analysis was performed with filtered EMTome genes using Monocle 2 model. Pairs of epithelial and mesenchymal marker genes, S100a4 vs Epcam (left) and Vim vs Krt18 (right) are used to highlight the MET status of individual cells.

C. Elevated ribosome biogenesis pathway during MET. Cells were classified by their MET state. The activation of RiBi pathway was indicated by AUC value of genes in the ribosome biogenesis pathway. Of note, the highest RiBi expression in MET State 8 cells, which exhibit the most epithelial phenotypes.

D. Scatter plot illustrates the correlation between RiBi activity and MET pseudotime. Proliferating cells (in S phase) are highlighted in red color.

Elevated RNA transcription activity in the Doub+ Tri-PyMT cells.

Tri-PyMT cells (p7) including both RFP+ and GFP+ cells were cultured in growth medium. EU incorporation assay was performed and detected with click chemistry using Alexa A647-azide. Labeled cells were analyzed by flow cytometry (Left). The median fluorescence intensity (MFI) of EU staining of RFP+, GFP+ and Doub+ cells were quantified, n =3, One-way ANOVA, *P < 0.0001, Doub+ vs RFP+; *P < 0.0001, Doub+ vs GFP+.

RiBi inhibition impairs RNA transcription activity and reduces EMT capability of tumor cells.

A, Fluorescence switch in Tri-PyMT cells with CX5461 treatment. Flow cytometry plots exhibit the percentage of GFP+ Tri-PyMT cells after treatment with the Pol I inhibitor CX5461. Biological repeat, n=2, One-way ANOVA, *P=0.0015, 20nM vs Cont; *P=0.0011, 100nM vs Cont.

B, EU incorporation assay. Tri-PyMT cells (p10) were treated with RiBi inhibitors, BMH21 at 100nM, or CX5461 at 20nM, for 2 days. EU incorporation assay was performed and detected with click chemistry using Alexa A647-azide. Labeled cells were analyzed by flow cytometry (Left). The median fluorescence intensity (MFI) of EU staining of cells were quantified, n =3, One-way ANOVA, P = 0.0001, +BMH21 vs Cont; P = 0.0032, +CX5461 vs Cont.

C, Cell migration assay. Tri-PyMT cells (p7) were seeded in the inserts of migration plate and treated BMH21 (100nM) or CX5461 (20nM) for 5 days. Cell migration was induced by FBS gradient for 24 hours. Migrating cells were counted from the bottom of the insert. n=6, One-way ANOVA, * P < 0.001 BMH21 vs Cont and CX5461 vs Cont.

RiBi modulating by knocking down ribosome proteins reduces EMT/MET capability of tumor cells.

A, Western blots show the knockdown expression of Rps24 and Rps28 in targeted Tri-PyMT cells. The numbers indicate relative intensity of the band normalizing to the β-actin band of the sample.

B, Impact on nucleoli numbers by Rps knockdown in RFP+ Tri-PyMT cells. Fluorescent images show the reduced number of nucleoli in Rps24 and Rps28 knockdown cells. n = 623 cells (shScr), 490 cells (shRps24), 553 cells (shRps28). One-way ANOVA with Dunnett multiple comparisons, *P < 0.0001 for shScr vs shRps24 and shScr vs shRps28.

C, Flow cytometry plots show the percentage of GFP+ cells in Rps24 and Rps28 knockdown Tri-PyMT cells. 4 biological repeats, One-way ANOVA, *P=0.0085, shScr vs. shRps24, *P=0.0117, shScr vs. shRps28.

D, Flow cytometry plots show the Epcam+/GFP+ cells in Rps24 or Rps28 knockdown Tri-PyMT cells. 3 biological repeats, one-way ANOVA, *P=0.0055, shScr vs. shRps24, *P=0.0124, shScr vs. shRps28.

Combination therapy of RNA Pol I inhibitor and chemo drugs.

Left, Bar graphs display Tri-PyMT cell viability following combination treatment with BMH21 and chemotherapeutic drugs at various concentrations. Bars represent the mean value of duplicated treatments, n=2 wells/treatment. Right, Synergy plots of BMH21 and chemotherapeutic drugs. Synergistic scores (δ) for each combination were calculated using the ZIP model in SynergyFinder 3.0. The highest δ numbers among combinations is shown.

Differential impact of Pol I inhibitor and chemo drug on epithelial and mesenchymal Tri-PyMT cells in vivo.

The same number of RFP+ and GFP+ Tri-PyMT cells were injected in animals via tail vein. Animals were treated CTX (100mg/kg, once a week, i.p.), BMH21 (20mg/kg, 5 times a week, i.p.) or both in combination. Number of RFP+ and GFP+ cells were quantified by flow cytometry at day 21 after inoculation.

A, Plots show the number of RFP+ and GFP+ Tri-PyMT cells in the lung treated with CTX and/or BMH21. n=3 mice/treatment, One way ANOVA, *P<0.0001 for Cont vs CTX, Cont vs BMH21 and Cont vs CTX+BMH21.

B, The ratio of GFP+ vs RFP+ cells in the lung treated with CTX and/or BMH21. n=3 mice/treatment, One way ANOVA, *P=0.0001 for Cont vs CTX; *P=0.0661 for Cont vs BMH21 and *P<0.0001 for Cont vs CTX+BMH21.

E-cadherin expression by LM2 cells in lung metastases.

qRT-PCR analysis of Ecadherin expression in LM2 cells isolated from lung metastases. LM2 cells in culture served as control. GAPDH served as internal control for PCR. n=3 experiments, unpaired t test, * P=0.0025.

Correlation of ribosome biogenesis pathway and EMT status in human breast cancer cells.

Single-nucleus cell RNA sequencing analysis of primary tumor cells from breast cancer patients (BC-A and BC-B, GEO database, GSE 198745).

A,D, The t-SNE plots show the heterogeneity in EMT statuses of tumor cells with AUC values of epithelial genes (in red) and mesenchymal genes (in green).

B,E, Cell trajectory analysis was performed with filtered EMT-related genes (EMTome genes) using Monocle 2 model. Cell states were identified with the most epithelial cells as root for calculation of EMT pseudotime (PsT).

C,F, The scatter plot displays the correlation of Ribi activity to EMT pseudotime. The polynomial regression line (order = 3) highlights the higher RiBi activity in cells toward epithelial phenotypes.

Survival curves of breast cancer patients with different RiBi activities.

Gene expression data of METABRIC and TCGA-PanCancer was obtained from CBioPortal. Patients were grouped according to the average z-score of RiBi gene expression, RiBi High (>1) and RiBi Low (RiBi score <0.5). Upper, METABRIC dataset, Gehan-Breslow-Wilcoxon test P-value = 0.0019. Lower, TCGA data, Gehan-Breslow-Wilcoxon test P-value = 0.0225.

Working model of ribosome biogenesis in EMP plasticity.

Tumor cells exhibit plasticity by undergoing EMT (epithelial-mesenchymal transition) and MET (mesenchymal-epithelial transition), contributing to the development of resistance to chemotherapies. Elevated ribosome biogenesis is essential for maintaining this EMT plasticity. Inhibition of ribosome biogenesis by RNA Pol I inhibitors synergizes with chemotherapeutic drugs to overcome resistance development.