RGS10 deficiency facilitates distant metastasis by inducing epithelial–mesenchymal transition in breast cancer

  1. Yang Liu
  2. Yi Jiang
  3. Peng Qiu
  4. Tie Ma
  5. Yang Bai
  6. Jiawen Bu
  7. Yueting Hu
  8. Ming Jin
  9. Tong Zhu
  10. Xi Gu  Is a corresponding author
  1. Department of Oncology, Shengjing Hospital of China Medical University, China
  2. Department of Anesthesiology, Shengjing Hospital of China Medical University, China
  3. Department of Pathology, Shengjing Hospital of China Medical University, China
  4. Department of Nursing, Shengjing Hospital of China Medical University, China
  5. Breast Surgery of Panjin Central Hospital, China
7 figures, 6 tables and 1 additional file

Figures

The expression and prognostic associations of RGS10 in breast cancer.

(A) RGS10 mRNA levels in 31 normal human tissues. Data were derived from the Genotype-Tissue Expression database. (B) RGS10 mRNA levels in cell lines representing 21 human cancers. Data were derived from the Cancer Cell Line Encyclopedia database. (C) qRT-PCR showing RGS10 mRNA levels in freshly resected breast cancer tissues (n = 20) and matched adjacent normal breast tissues. **p<0.01, Student’s t -test. (D, E) Survival analyses showing disease-free survival (DFS) (D) and overall survival (OS) (E) in patients with breast cancer stratified by high versus low RGS10 mRNA levels. Data were derived from the Kaplan–Meier plotter database. (F) Representative images showing immunohistochemical staining of RGS10 protein expression in breast cancer tissues or normal tissues (n = 133) (magnification: ×200 and ×400). (G, H) Kaplan–Meier analysis showing DFS (G) and OS (H) in patients with breast cancer stratified by presence versus absence of RGS10 protein in breast cancer tissues (n = 133).

RGS10 silencing increases the proliferation and migration of breast cancer cells in vitro.

(A) Western blotting showing RGS10 protein levels in molecularly distinct breast cancer cell lines. The original files of the full raw unedited blots are provided in Figure 2—source data 1. The uncropped gels with the relevant bands labeled are provided in Figure 2—source data 2. The statistical data is provided in Figure 2—source data 3. (B) Western blotting showing RGS10 protein levels in SKBR3 cells transfected with three independent shRNA constructs, shRNA-RGS10-161, shRNA-RGS10-321, and shRNA-RGS10-506, and shRNA-NC. The original files of the full raw unedited blots are provided in Figure 2—source data 4. The uncropped gels with the relevant bands labeled are provided in Figure 2—source data 4. The statistical data is provided in Figure 2—source data 6. (C) CCK-8 assay showing the viability of SKBR3 cells transfected with shRNA-RGS10-161, shRNA-RGS10-506, or shRNA-NC. **p <0.01, one-way ANOVA. (D–F) Colony formation (D) and transwell migration/invasion (E, F) assays in SKBR3 cells transfected with shRNA-RGS10-161, shRNA-RGS10-506, or shRNA-NC. ***p<0.001, one-way ANOVA.

Figure 2—source data 1

Original files for the gels in Figure 2A.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data1-v1.zip
Figure 2—source data 2

Uncropped gels with the relevant bands labeled in Figure 2A.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data2-v1.zip
Figure 2—source data 3

Statistical data for Figure 2A.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data3-v1.pdf
Figure 2—source data 4

Original files for the gels in Figure 2B.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data4-v1.zip
Figure 2—source data 5

Uncropped gels with the relevant bands labeled in Figure 2B.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data5-v1.zip
Figure 2—source data 6

Statistical data for Figure 2B.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data6-v1.pdf
Figure 2—source data 7

Original files for Figure 2C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data7-v1.zip
Figure 2—source data 8

Original files for Figure 2D.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data8-v1.zip
Figure 2—source data 9

Original files for Figure 2E.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data9-v1.zip
Figure 2—source data 10

Original files for Figure 2F.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig2-data10-v1.zip
Protein–protein network interaction network and Gene Ontology (GO) functional- and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway-enrichment analysis of genes co-expressed with RGS10.

(A) Volcano plot showing differentially expressed genes between SKBR3 cells transfected with shRNA-RGS10 or shRNA-NC. (B–E) KEGG pathway analysis and GO enrichment analysis of differentially expressed genes showing the ten most enriched terms. BP: biological processes; MF: molecular function; CC: cellular compartment. (F) LCN2, E-cadherin, and vimentin protein levels in SKBR3 cells transfected with shRNA-RGS10 or shRNA-NC. The original files of the full raw unedited blots are provided in Figure 3—source data 1. The uncropped gels with the relevant bands labeled are provided in Figure 3—source data 2. The statistical data is provided in Figure 3—source data 3. (G) Schematics of predicted MIR539-5p binding sites between wild-type and mutant RGS10 sequences in the 3ʹ-untranslated regions. (H) Relative luciferase activities detected after cotransfection of wild-type or mutant luciferase reporter plasmids and an MIR539-mimic.

MIR539-5p regulates the migration, invasion, proliferation, and epithelial–mesenchymal transition (EMT) of breast cancer cells.

(A) qPCR showing the transfection efficiency of the MIR539-5p mimic. ***p<0.001, one-way ANOVA. (B, C) qRT-PCR and western blotting showing RGS10 mRNA and protein levels in SKBR3 cells transfected with the MIR539-5p mimic, negative control (NC), or wild type (WT). ***p<0.001, one-way ANOVA. The original files of the full raw unedited blots are provided in Figure 4—source data 1. The uncropped gels with the relevant bands labeled are provided in Figure 4—source data 2. The statistical data is provided in Figure 4—source data 3. (D) CCK-8 assay showing the viability of SKBR3 cells transfected with the MIR539-5p mimic, NC, or WT. ***p<0.001, one-way ANOVA. (E–G) Colony formation (E) and transwell migration/invasion (F, G) assays in SKBR3 cells transfected with the MIR539-5p mimic, NC, or WT. *p<0.05, **p<0.01, Student’s t-test. (H) Western blotting showing protein levels of LCN2 and biomarkers of EMT in SKBR3 cells transfected with the MIR539-5p mimic or NC. The original files of the full raw unedited blots are provided in Figure 4—source data 4. The uncropped gels with the relevant bands labeled are provided in Figure 4—source data 5. The statistical data is provided in Figure 4—source data 6. (I) Immunofluorescence staining showing E-cadherin, vimentin, and snail protein expression in SKBR3 cells transfected with the MIR539-5p mimic or NC. Scale bar: 50 µm.

Figure 4—source data 1

Original files for the gels in Figure 4C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data1-v1.zip
Figure 4—source data 2

Uncropped gels with the relevant bands labeled in Figure 4C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data2-v1.zip
Figure 4—source data 3

Statistical data for Figure 4C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data3-v1.pdf
Figure 4—source data 4

Original files for the gels in Figure 4H.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data4-v1.zip
Figure 4—source data 5

Uncropped gels with the relevant bands labeled in Figure 4H.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data5-v1.zip
Figure 4—source data 6

Statistical data for Figure 4H.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data6-v1.pdf
Figure 4—source data 7

Original files for Figure 4A, B and D.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data7-v1.zip
Figure 4—source data 8

Original files for Figure 4E.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data8-v1.zip
Figure 4—source data 9

Original files for Figure 4F.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data9-v1.zip
Figure 4—source data 10

Original files for Figure 4G.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data10-v1.zip
Figure 4—source data 11

Original files for Figure 4I.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig4-data11-v1.zip
MIR539-5p inhibitor suppresses breast cancer cell proliferation and invasion.

(A) qPCR showing transfection efficiency of the MIR539-5p inhibitor after 48 hr. ***p<0.001, one-way ANOVA. (B, C) qRT-PCR and western blotting showing RGS10 mRNA and protein levels in MDA-MB-231 cells transfected with the MIR539-5p inhibitor, negative control (NC), or wild type (WT). ***p<0.001, one-way ANOVA. The original files of the full raw unedited blots are provided in Figure 5—source data 1. The uncropped gels with the relevant bands labeled are provided in Figure 5—source data 2. The statistical data is provided in Figure 5—source data 3. (D) CCK-8 assay showing the viability of MDA-MB-231 cells transfected with the MIR539-5p inhibitor, NC, or WT. ***p<0.001, one-way ANOVA. (E–G) Colony formation (E) and transwell migration/invasion (F, G) assays in MDA-MB-231 cells transfected with the MIR539-5p inhibitor, NC, or WT. *p<0.05, **p<0.01, Student’s t-test. (H) Western blotting showing protein levels of LCN2 and biomarkers of EMT in MDA-MB-231 cells transfected with the MIR539-5p inhibitor or NC. The original files of the full raw unedited blots are provided in Figure 5—source data 4. The uncropped gels with the relevant bands labeled are provided in Figure 5—source data 5. The statistical data is provided in Figure 5—source data 6. (I) Immunofluorescence staining showing E-cadherin, vimentin, and snail protein expression in MDA-MB-231 cells transfected with the MIR539-5p inhibitor or NC. Scale bar: 50 µm.

Figure 5—source data 1

Original files for the gels in Figure 5C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data1-v1.zip
Figure 5—source data 2

Uncropped gels with the relevant bands labeled in Figure 5C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data2-v1.zip
Figure 5—source data 3

Statistical data for Figure 5C.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data3-v1.pdf
Figure 5—source data 4

Original files for the gels in Figure 5H.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data4-v1.zip
Figure 5—source data 5

Uncropped gels with the relevant bands labeled in Figure 5H.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data5-v1.zip
Figure 5—source data 6

Statistical data for Figure 5H.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data6-v1.pdf
Figure 5—source data 7

Original files for Figure 5A, B and D.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data7-v1.zip
Figure 5—source data 8

Original files for Figure 5E.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data8-v1.zip
Figure 5—source data 9

Original files for Figure 5F.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data9-v1.zip
Figure 5—source data 10

Original files for Figure 5G.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data10-v1.zip
Figure 5—source data 11

Original files for Figure 5I.

https://cdn.elifesciences.org/articles/97327/elife-97327-fig5-data11-v1.zip
RGS10 inhibits breast cancer growth by targeting LCN2 in vivo.

(A) Size of tumors derived from RGS10-depleted SKBR3 cells, negative control (NC), and wild type (WT). (B) Volume of tumors derived from RGS10-depleted SKBR3 cells, NC, and WT. ***p<0.001, one-way ANOVA. (C) Hematoxylin and eosin staining of tumors derived from RGS10-depleted SKBR3 cells and NC. (D) Immunohistochemical staining showing LCN2, E-cadherin, snail, and vimentin protein expression in tumors derived from RGS10-depleted SKBR3 cells and NC.

Author response image 1

Tables

Table 1
Clinicopathological characteristics of 20 patients.
Clinicopathological characteristicsPatient 1Patient 2Patient 3Patient 4Patient 5
SexFemaleFemaleFemaleFemaleFemale
Age (years)5840703858
T grade21211
N grade01301
Menopausal statusPostmenopausalPremenopausalPostmenopausalPremenopausalPostmenopausal
SubclassLuminal ALuminal ALuminal ALuminal ALuminal A
MetastasisNoNoNoNoNo
Patient 6Patient 7Patient 8Patient 9Patient 10
SexFemaleFemaleFemaleFemaleFemale
Age (years)5956805537
T grade22321
N grade30223
Menopausal statusPostmenopausalPostmenopausalPostmenopausalPostmenopausalPremenopausal
SubclassLuminal BLuminal BLuminal BLuminal BLuminal B
MetastasisNoNoNoNoNo
Patient 11Patient 12Patient 13Patient 14Patient 15
SexFemaleFemaleFemaleFemaleFemale
Age (years)6257415358
T grade22212
N grade31331
Menopausal statusPostmenopausalPostmenopausalPremenopausalPostmenopausalPostmenopausal
SubclassHER2 positiveHER2 positiveHER2 positiveHER2 positiveHER2 positive
MetastasisNoNoNoNoNo
Patient 16Patient 17Patient 18Patient 19Patient 20
SexFemaleFemaleFemaleFemaleFemale
Age (years)5037355155
T grade22222
N grade00000
Menopausal statusPremenopausalPremenopausalPremenopausalPremenopausalPostmenopausal
SubclassesTriple negativeTriple negativeTriple negativeTriple negativeTriple negative
MetastasisNoNoNoNoNo
  1. HER2, human epidermal growth factor receptor 2.

Table 2
Correlations between RGS10 expression and clinicopathological characteristics.
VariableRGS10 expression (n, %)RGS10 expression (n, %)p-Value
No. of patients34 (25.6)99 (74.4)
Age (years)0.590
≤458 (23.5)28 (28.3)
>4526 (76.5)71 (71.7)
Other disease0.746
Yes28 (82.4)79 (79.8)
No6 (17.6)20 (20.2)
Histological grade0.744
I12 (35.3)30 (30.3)
II19 (55.9)56 (56.6)
III3 (8.8)13 (13.1)
Tumor size (cm)0.101
Median (range)2.191 (1.0–8.0)1.861 (0.6–4.0)
No. of positive axillary lymph nodes0.316
Median (range)1.06 (0–19)1.93 (0–27)
No. of pregnancies0.706
Median (range)2.50 (1–6)2.40 (0–7)
No. of births0.497
Median (range)1.62 (1–3)1.75 (0–7)
Age at first pregnancy (years)0.489
<211 (2.9)6 (6.1)
21–2515 (44.1)42 (42.4)
26–3018 (52.9)46 (46.5)
>300 (0)5 (5.0)
Menopausal status0.216
Premenopausal21 (61.8)49 (49.5)
Postmenopausal13 (38.2)50 (50.5)
Estrogen receptor status0.397
Positive25 (73.5)65 (65.7)
Negative9 (26.5)34 (34.3)
Progesterone receptor status0.527
Positive25 (73.5)78 (78.8)
Negative9 (26.5)21 (21.2)
HER2 status0.496
Positive6 (17.6)23 (23.2)
Negative28 (82.4)76 (76.8)
Ki67 status0.773
>2014 (41.2)38 (38.4)
≤2020 (58.8)61 (61.6)
Molecular type0.043
Luminal A6 (17.6)30 (30.3)
Luminal B19 (55.9)36 (36.3)
HER2+2 (5.9)16 (16.2)
Triple-negative breast cancer7 (20.6)17 (17.2)
Distant metastasis0.008
Yes7 (20.6)46 (46.5)
No27 (79.4)53 (53.5)
Disease-free survival (months)<0.001
Median (range)114.71 (53–140)89.66 (9–132)
Death0.024
Yes5 (14.7)35 (35.4)
No29 (85.3)64 (64.6)
Overall survival (months)0.030
Median (range)120.41 (80–140)107.64 (24–140)
  1. HER2, human epidermal growth factor receptor 2.

Table 3
Univariate and multivariate Cox regression analyses of clinicopathological risk factors for disease-free survival (DFS).
VariableDFS
Univariate analysisMultivariate analysis
HR (95% CI)p-ValueHR (95% CI)p-Value
Age (years)0.720 (0.400–1.296)0.273NA
Other disease0.885 (0.445–1.762)0.728NA
Histological grade
I0.0030.004
II2.636 (1.263–5.503)0.0102.704 (1.294–5.649)0.008
III4.638 (1.877–11.462)0.0014.566 (1.846–11.294)0.001
Tumor size (cm)1.065 (0.840–1.350)0.605NA
No. of positive axillary lymph nodes1.012 (0.953–1.074)0.698NA
No. of pregnancies0.835 (0.662–1.054)0.129NA
No. of births0.694 (0.481–1.001)0.051NA
Age at first pregnancy (years)NA
<210.501
21–253.336 (0.449–24.810)0.239
26–303.886 (0.530–28.499)0.182
>301.838 (0.115–29.429)0.667
Menopausal status0.844 (0.492–1.451)0.540NA
Estrogen receptor status1.534 (0.820–2.869)0.180NA
Progesterone receptor status1.098 (0.577–2.089)0.776NA
HER2 status0.841 (0.433–1.635)0.610NA
Ki67 status0.816 (0.465–1.431)0.478NA
Molecular typeNA
Luminal A0.307
Luminal B1.940 (0.819–4.593)0.132
HER2+1.147 (0.479–2.747)0.758
Triple-negative breast cancer1.453 (0.527–4.007)0.470
RGS10 expression0.323 (0.146–0.716)0.0050.321 (0.144–0.713)0.005
  1. HR, hazard ratio; HER2, human epidermal growth factor receptor 2; NA, not analyzed.

Table 4
Univariate and multivariate Cox regression analyses of clinicopathological risk factors for overall survival.
VariableOverall survival
Univariate analysisMultivariate analysis
HR (95% CI)p-ValueHR (95% CI)p-Value
Age (years)0.423 (0.226–0.792)0.0070.305 (0.154–0.603)0.001
Other disease0.984 (0.453–2.136)0.967NA
Histological grade
I0.0670.008
II1.519 (0.703–3.284)0.2881.834 (0.840–4.007)0.128
III3.062 (1.180–7.943)0.0215.090 (1.812–14.300)0.002
Tumor size (cm)1.045 (0.788–1.385)0.760NA
No. of positive axillary lymph nodes1.011 (0.943–1.083)0.759NA
No. of pregnancies0.941 (0.731–1.213)0.640NA
No. of births0.896 (0.625–1.285)0.550NA
Age at first pregnancy (years)NA
<210.580
21–252.818 (0.378–21.001)0.312
26–302.021 (0.269–15.190)0.494
>303.255 (0.295–35.934)0.335
Menopausal status0.769 (0.405–1.458)0.420NA
Estrogen receptor status0.769 (0.405–1.458)0.420NA
Progesterone receptor status1.091 (0.519–2.293)0.818NA
HER2 status0.587 (0.246–1.398)0.229NA
Ki67 status0.718 (0.370–1.392)0.327NA
Molecular typeNA
Luminal A0.589
Luminal B0.639 (0.295–1.383)0.255
HER2+0.756 (0.272–2.100)0.592
Triple-negative breast cancer1.079 (0.467–2.493)0.860
RGS10 expression0.353 (0.138–0.900)0.0290.364 (0.154–0.603)0.001
  1. HR, hazard ratio; HER2, human epidermal growth factor receptor 2; NA, not analyzed.

Table 5
Gene sets enriched in phenotype ‘high’.
MSigDB collectionGene set nameNESNOMp-valueFDR q-value
h.all.v6.2.symbols.gmt [Hallmarks]HALLMARK_ALLOGRAFT_REJECTION–2.5500
HALLMARK_APOPTOSIS–2.4100
HALLMARK_IL6_JAK_STAT3_SIGNALING–2.4000
HALLMARK_IL2_STAT5_SIGNALING–2.3900
HALLMARK_INFLAMMATORY_RESPONSE–2.3700
HALLMARK_INTERFERON_GAMMA_RESPONSE–2.3300
HALLMARK_COMPLEMENT–2.2700
HALLMARK_KRAS_SIGNALING_UP–2.2500
HALLMARK_INTERFERON_ALPHA_RESPONSE–2.2200
HALLMARK_APICAL_JUNCTION–2.1300
  1. Gene sets with NOM p-value <0.05 and FDR q-value <0.25 were considered as significant.

  2. MSigDB, Molecular Signatures Database; NES, normalized enrichment score; NOM, nominal; FDR, false discovery rate.

Key resources table
Reagent type (species) or resourceDesignationSourceor referenceIdentifiersAdditional information
Cell line (SKBR3)HER2+ breast cancer cell lineATCCCell line cultured in McCoy’s 5A
Culture medium supplemented with 10% FBS
Cell line (MCF-7)ER+/PR+ breast cancer cell lineATCCCell line cultured in DMEM high glucose
Culture medium supplemented with 10% FBS
Cell line (MDA-MB-
231)
Triple-negative breast cancer cell lineATCCCell line cultured in L-15
Culture medium supplemented with 10% FBS
Transfected construct (human)RGS10 shRNA#161,321,506Syngentech BiotechnologyLentiviral construct to transfect and express the shRNA
Transfected construct (human)MIRNA-539-5p mimics and inhibitorGenePharma
Biotechnology
AntibodyAnti-
RGS10 (rabbit monoclonal)
AbcamCat# ab154172WB (1:250)
IHC (1:50)
AntibodyAnti- lipocalin-2 (rabbit monoclonal)ABclonalCat#
A2092
WB (1:1000)
IHC (1:500)
AntibodyAnti-
E-cadherin
(rabbit monoclonal)
ProteintechCat#
20874-1-AP
WB (1:1000)
IHC (1:500)
IF (1:100)
AntibodyAnti-vimentin
(mouse monoclonal)
HUABIOCat#
EM0401
WB (1:1000)
IHC (1:200)
IF (1:100)
AntibodyAnti-snail
(rabbit monoclonal)
ABclonalCat#
A5243
WB (1:1000)
IHC (1:00)
IF (1:100)
AntibodyAnti-GAPDH
(rabbit monoclonal)
CSTCat#
5174
WB (1:1000)
AntibodyAnti-
β-tubulin (rabbit monoclonal)
ABclonalCat#
AC008
WB (1:5000)
Sequence-based reagentRGS10_FThis paperPCR primersCACGACAGCGATGGCAGTTCC
Sequence-based reagentRGS10_RThis paperPCR primersCTTTTCACGCCTTCTGGGTCTTCC
Sequence-based reagentMIR539-5p_FThis paperPCR primersCGCTGCATGGAGAAATTATCCT
Sequence-based reagentMIR539-5p_RThis paperPCR primersCAGTGCAGGGTCCGAGGT
Software, algorithmSPSSSPSSSPSS, version 25.0
Armonk, NY

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Yang Liu
  2. Yi Jiang
  3. Peng Qiu
  4. Tie Ma
  5. Yang Bai
  6. Jiawen Bu
  7. Yueting Hu
  8. Ming Jin
  9. Tong Zhu
  10. Xi Gu
(2024)
RGS10 deficiency facilitates distant metastasis by inducing epithelial–mesenchymal transition in breast cancer
eLife 13:RP97327.
https://doi.org/10.7554/eLife.97327.3