|Home | About | Journals | Submit | Contact Us | Français|
High levels of SOX2 protein are correlated with increased dissemination of breast cancer. However, the underlying molecular mechanisms are not fully understood.
In this study we investigate the role of SOX2 in breast cancer metastasis using multiple in vitro and in vivo assays including cell culture, shRNA-mediated knockdown, wound healing, colony formation, transwell chamber, xenograft and tail vein injection. Moreover, western blot, immunostaining, microarray and real-time PCR were used to determine the change of protein and miRNA levels. Luciferase assays were also used to evaluate activity which TUSC3 is a target of miR-181a-5p and miR-30e-5p, and the clinical survival relevance was analyzed by Kaplan-Meier analysis.
We identified a novel pathway involving SOX2 regulation of microRNAs to control the proliferation and migration of breast cancer cells. shRNA-mediated knockdown of SOX2 inhibits breast cancer cell expansion and migration. More importantly, we found that these changes are accompanied by significant reduction in the levels of two microRNAs, miR-181a-5p and miR-30e-5p. Overexpression of these two microRNAs leads to reduced protein levels of Tumor Suppressor Candidate 3 (TUSC3) in breast cancer cells; mutations of the potential binding sites in the 3’-UTR of TUSC3 abrogate the inhibitory effects of the microRNAs. We further found that upregulation of TUSC3 expression leads to reduced proliferation and migration of breast cancer cells. In human breast cancer samples the levels of TUSC3 protein are inversely correlated with those of SOX2 protein.
Taken together, our work reveals a novel SOX2-mediated regulatory axis that plays critical roles in the proliferation, migration and invasiveness of breast cancer cells. Targeting this axis may provide beneficial effect in the treatment of breast cancer.
The online version of this article (doi:10.1186/s12943-017-0632-9) contains supplementary material, which is available to authorized users.
Breast cancer is a leading cause of cancer-related death in women worldwide. Despite the significant progress that has been made in cancer treatment and early detection, approximately 30 to 70% of breast cancer patients die from breast cancer recurrence or metastasis . Therefore it is important to investigate the mechanisms underlying cancer cell proliferation and metastasis-related cellular activities such as cell migration and invasion. Genetic mutations or epigenetic modifications that provide growth advantage are critical for tumour initiation and progression. In breast cancer, activation mutations of genes such as EGFR, HER2 have been found to transmit signals into the nucleus to up-regulate the transcription of downstream genes, promoting cancer cell growth [2–4]. However, little is known about the mediators in the nucleus that are involved in gene regulation. Interestingly, the transcription factor SOX2 has recently been associated with breast cancer formation and metastasis . SOX2 belongs to the sex determining region Y (SRY)-like box (SOX) gene family, members of which are pivotal for normal development and maintenance of stemness [6–11], cell proliferation and differentiation . Emerging studies demonstrated that the SOX2 protein is also involved in tumourigenesis [13–15]. High levels of SOX2 closely correlate with multiple processes during tumour development, including initiation , maintenance [17, 18], invasion and metastasis [19–21]. Accordingly, aberrant SOX2 expression was detected in multiple types of cancers at different stages. We and others have shown high levels of SOX2 in lung cancer, esophageal squamous cell carcinomas (ESCC) and ovarian cancer [13, 22–26]. Interestingly, high levels of SOX2 proteins have also been found in 19% of breast cancer patients . Moreover, SOX2 expression has been linked to tamoxifen resistance and relapse in breast cancer treatment . Nevertheless, the mechanisms by which high levels of SOX2 regulate the progression and metastasis of breast cancer remain largely unexplored.
MicroRNAs (miRNAs) are small, noncoding RNAs that are important for the regulation of gene expressions at the post-transcriptional level. miRNAs are involved in almost all key cellular activities such as proliferation, differentiation and migration . Dysregulation of miRNAs is also involved in the initiation and progression of human cancers including breast cancer [29–31]. For example, miR-124 is down-regulated in breast cancer cell lines. Ectopic expression of miR-124 reduces the protein levels of Myc and phospho-Rb while up-regulating p27 in the breast cancer cell line MDA-MB-231 , suggesting that miRNAs function as a tumour suppressor. In this study we examined the molecular mechanisms used by SOX2 in the proliferation and migration of breast cancer cell lines and identified that SOX2 regulates two miRNAs (miR-181a-5p and miR30e-5p) which in turn influence the expression of a common downstream target TUSC3. Importantly, inversed expression of SOX2 and TUSC3 is associated with poor prognosis of a subpopulation of breast cancer patients. Therefore we demonstrate here a clinically relevant axis which involves SOX2/miRNAs/TUSC3 in breast cancer development.
ZR7530 (Cat NO: CRL-1504™) and MDA-MB-231 (Cat NO: HTB-26™) breast cancer cell lines were purchased from ATCC (Manassas, VA, USA), Cell line authentication was conducted by short tandem repeat analysis, and mycoplasma contamination in cell lines was tested using multiplex PCR. They were infected with lentivirus carrying pLK0.1SOX2 (SOX2 knockdown plasmid), pCDH-SOX2-IRES-GFP (SOX2 overexpression plasmid), pLK0.1Scramble (control plasmid) , and pCDH-TUSC3-IRES-GFP (TUSC3 overexpression plasmid, primers used for construction were shown in Additional file 1: Table S1). Puromycin (Santa Cruz Biotechnology, Cat No: sc-108071, Dallas, TX, USA) selection was used to establish cell lines stably expressing these constructs. Cell lines were cultured in Dulbecco’s Modified Eagle Medium (DMEM, HyClone, Cat No: SH30022.01B, Beijing, China) with 10% foetal bovine serum (FBS, Gibco, Cat No:10270-106, Grand Island, NY, USA). Nude mice at 6-weeks old were purchased from SLRC Laboratory Company (Shanghai, China).
Western blot was performed as previously described . PVDF membrane was used for protein transfer and probed with antibodies against human SOX2 (polyclonal rabbit anti-SOX2, SEVEN HILLS, Cat No: WRAB-1236, Cincinnati, USA), human β-actin (mouse anti-β-actin, Beyotime, Cat No: AA128, Nantong, China), CDK4 (BBI, Cat No: AB20396b, Shanghai, China), CDK6 (BBI, Cat No: AB20398a, Shanghai, China), CCND1 (BBI, Cat No: AB60236b, Shanghai, China), TUSC3 (Novus Biologicals, Cat No: NBP1-55630, Littleton, USA). HRP-conjugated goat anti-rabbit IgG (Abcam, Cat No: ab136817, Cambridge, MA, USA) and HRP-conjugated goat anti-mouse IgG (ZSGB-Bio, Cat No: ZB-2305, Beijing, China) were used as secondary antibodies to detect the proteins. The density of protein band was quantified with Image-Pro Plus 6.0 software, and the ratio of target protein to housekeeping protein, which reflects the change of expression level, was calculated.
ZR7530 or MDA-MB-231 cells with genetic manipulation (knockdown or overexpression) were seeded in 96-well plates. The OD value of cells was measured under 450 nm and 630 nm at 24, 48, 72 and 96 h using the Cell Counting Kit 8 (CCK-8) (Beyotime) assay, and the absorbance difference between 450 nm and 630 nm represents cell proliferation rates. For soft agar colony formation assays, ZR7530 or MDA-MB-231 cells were seeded in 6-well plates and maintained for 2 weeks and colonies were fixed with 4% paraformaldehyde (DingGuo, Cat No: AR-0211, Beijing, China) and stained with 0.5% crystal violet. The number of total colonies and colonies with a diameter greater than 0.5 mm were calculated with Image J (Fiji-win32 software, University of Wisconsin-Madison), and the data were analysed as previously described [35, 36].
1×106 ZR7530 or MDA-MB-231 cells (SOX2 knockdown, scramble and TUSC3 overexpression) were seeded in 6-well plates and grown to 90% confluence. Cell monolayers were scraped with a sterile pipette tip. Floating cells were removed and the cultures were maintained in DMEM supplemented with 5% FBS. The wound area was recorded 0 h and 48 h after the scrape. The healing index was calculated and analysed using the formula: (S0−Sn)/S0×100%, which represents the ability of cellular migration. S0 represents the blank area 0 h after scraping, and Sn represents the blank area n h after scraping .
Invasion assays were performed in chambers with a 6.5-mm insert in 24-well plates (Corning company, Cat No: 3422, Corning, NY, USA). An 8-μm polycarbonate membrane was coated with 12 μl of ice-cold Matrigel™ Basement Membrane Matrix (BD company, Cat No: 356234, Bedford, MA, USA). Cells (5x104 per well) were added to the upper chamber in 200 μl of the medium containing 5% FBS. The lower chamber was filled with 500 μl of medium containing 20% FBS. The culture was maintained for 24 h. Cells were then fixed in methanol and stained with 1% crystal violet. Cells at the lower side of the filter were counted under a light microscope and dissolved with 50% acetic acid.
1x105 ZR7530 cells were re-suspended in 200 μl of DMEM and subcutaneously injected into the flanks of 6-week old BALB/c nude mice, which were maintained in the SPF environment. Five female nude mice per cohort (SOX2 knockdown and control) were randomly divided to receive cancer cells. Tumour size was measured every 3 days using a Vernier calliper after injection, and tumour volumes were calculated with the formula 0.52×width2×length [38, 39]. Tumours were collected 7 weeks after cell inoculation.
For lung metastatic experiments, 2×106 tumour cells were re-suspended in 100 μl of PBS and injected into BALB/c nude mice via the tail vein. Five female nude mice per cohort, which were also maintained in the SPF environment, received the tumour cells through tail vein. These nude mice were sacrificed 8 weeks after tail vein injection. Lungs were fixed with 4% PFA and sectioned for haematoxylin and eosin (HE) staining to determine the metastatic capabilities of cancer cells. The metastatic capabilities were assessed as previously described [40, 41].
Total RNAs of ZR7530 cells (SOX2 knockdown and control) were extracted and purified using Trizol (Invitrogen, Carlsbad, CA, USA) and miRNeasy kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The human miRNA microarray miRCURY™ LNA Expression Array (v.18.0, Exiqon) was used to detect the level changes after SOX2 knockdown. The low intensity miRNAs were filtered, and original signal intensities were normalised using the Median Normalization Method. The miRNA data were analysed with Significance Analysis of Microarrays (SAM) version 2.1 to generate a list of miRNAs that up-regulate or down-regulate with a fold change ≥1.5 and a p-value ≤0.05. Quantitative real-time PCR analysis was used to validate the level changes of miRNA. Additional file 1: Table S2 and Table S3 describe the stem-loop RT primers and specific PCR primers used for quantitative PCR.
Differentially expressed miRNAs were clustered using the Cluster software, and we performed miRNA target gene predictions using the microRNA target prediction tools miRBase, TargetScan, and miRanda, and all common target genes which are regulated by two or multiple miRNAs were analyzed with cytoscape program. We also performed Gene Ontology (GO) analysis (http://www.geneontology.org) of the target genes that were predicted using the three databases, and the KEGG database (http://www.genome.jp/kegg/) was selected for pathway analysis based on the target genes of these differentially expressed miRNAs. The significance of the GO terms and Pathway was selected based on a p-value cut-off of 0.05.
Luciferase reporter activities were determined using the Dual-Luciferase Reporter Assay System (Promega, Cat No. E1910, Madison, USA). The primers used for obtaining mutated 3’UTR of TUSC3 were described in Additional file 1: Table S4. pGL3 with wild-type or mutated 3’UTR of TUSC3, miRNA mimics and pRL-TK were co-transfected into cells for all luciferase assays, and Renilla luciferase activities were measured to assess the transfection efficiency, sequences of miRNA mimics, inhibitors and relevant controls were shown in Additional file 1: Table S5.
The data represent at least three independent experiments using cells or extracts from a minimum of three separate isolations. Differences between groups were compared using analysis of variance for repeated measures. All statistical analyses were performed using GraphPad PRISM. v5.0 software (San Diego, CA, USA). Data were presented as mean±SD (standard deviation). Statistical significance between two groups was calculated by unpaired student’s t-test, and P value<0.05 is considered significant.
The breast cancer cells ZR-7530 and MDA-MB-231 express high levels of SOX2 protein. To test whether SOX2 is required for the proliferation of these cell lines, we performed lentivirus-delivered SOX2 shRNA knockdown in ZR-7530 as previously described . Knockdown leads to up to 85–92% reduction in the levels of SOX2 protein (p<0.01 for ZR-7530 and p<0.001 for MDA-MB-231) (Fig. 1a). More importantly, decreased protein levels are associated with decreased proliferation of cancer cells, especially 96 h after knockdown (p<0.001 for ZR-7530 and p<0.001 for MDA-MB-231) (Fig. 1b). The decrease in cell proliferation is correlated with reduced levels of cell cycle protein CCND1, CDK4 and CDK6 (Fig. 1d). Similarly, SOX2 knockdown inhibits the proliferation of MDA-MB-231 cells (Fig. 1b right). We also performed colony forming assay to test whether SOX2 knockdown affects colony formation from single cells. Reduced levels of SOX2 affect the number and size of colonies. The number of colonies with a diameter greater than 0.5 mm formed from ZR7530 knockdown cells was significantly reduced in the SOX2-knockdwon (SOX2-Kd) group compared to scramble controls (P<0.05) (Fig. 1c). Similarly, SOX2 knockdown in MDA-MB-231 cells also significantly reduced the number of colonies with diameters greater than 0.5 mm (P<0.05). Together these data confirm that SOX2 is required for the proliferation of breast cancer cells.
High levels of SOX2 have been found in metastasized breast cancer nodules . We next asked whether SOX2 is involved in cellular activities (e.g. cell migration) relevant to metastasis. SOX2 knockdown results in a decrease in the migration capability of ZR7530 and MDA-MB-231 cells (Fig. 2a). Wound-healing assay showed that the healing index was reduced from 24.75% (controls) to 3.33% for ZR7530 cells (p<0.001), and from 57.39 to 43.8% for MDA-MB-231 cells (p<0.05) (Fig. 2b). In addition, SOX2 knockdown inhibits the invasion of breast cancer cells in transwell assay. The number of cancer cells passing through the filter of the transwell are significantly reduced after SOX2 knockdown (Fig. 2c, ,d,d, p<0.01 for ZR7530 and p<0.05 for MDA-MB-231). The changed tendency of the OD value of eluted crystal violet was consistent with the number of migrated cells in these groups (p<0.01 for ZR7530 and p<0.05 for MDA-MB-231) (Fig. 2e). Together, these in vitro results demonstrate that high levels of SOX2 protein are critical for migration and invasion of breast cancer cells.
We next used a xenograft model to test whether SOX2 knockdown reduces tumour growth in vivo. SOX2 knockdown ZR7530 and control cells were transplanted into immunocompromised nude mice. Significantly, SOX2 knockdown leads to 4.4 and 4.7 fold decreases in the volume and weight of tumour, respectively (Fig. 3a–c) (n=5). In addition, SOX2 overexpression enables higher tumourigenesis capability of breast cancer cells. We transplanted SOX2-overexpressing ZR7530 cells into nude mice, and SOX2 overexpression results in a 2.3 fold increase in tumour volume and a 3.6 fold increase in tumour weight (Fig. 3a–c) (n=5). These results suggest that the levels of SOX2 are correlated with increased growth capability of breast cancer cells in vivo.
Our in vitro results show that SOX2 knockdown reduces the migration of breast cancer cell lines. To test whether decreased levels of SOX2 also block metastasis capability in vivo, we performed tail vein injection of ZR7530 and studied tumour seeding and development in the lung. Upon SOX2 knockdown, the number of tumours seeded in the lung are significantly reduced (33±2 Vs 7±3) (p<0.001) (n=5) (Fig. 3d, ,e).e). In addition, SOX2 knockdown also reduces the ratio of tumour/total lung area (10.4% Vs 1.8%) (p<0.001) (Fig. 3f). By contrast, ectopic SOX2 expression enables high metastastic capability of ZR7530 cells, leading to an increased number of seeding tumours (33±2 Vs 51±4) (p<0.01) and enlarged metastatic foci in the lung (10.4% vs 27.5%) (p<0.001) (Fig. 3d–f). Together, these in vivo findings support the important roles of SOX2 in promoting metastasis of breast cancer cells.
miRNAs have been implicated in the regulation of multiple steps in breast cancer progression [43–45]. To test whether SOX2 regulates proliferation and metastasis through miRNAs, we performed microRNA microarray analysis (miRCURY™ LNA Expression Array v.18.0, Exiqon). Total RNAs were isolated from scramble control and SOX2 knockdown ZR7530 cells. SOX2 knockdown significantly changed the levels of multiple miRNAs, among which 11 were up-regulated and 10 were significantly down-regulated (Fig. 4a and Additional file 1: Table S6, p<0.05). We validated the expression levels with real-time PCR (Fig. 4b). Although the levels of miR-338-3p are significantly down-regulated in the array analysis, the levels of miR-30e-5p are most dramatically down-regulated upon SOX2 knockdown, which was validated by real-time PCR (Fig. 4b). Notably, the levels of miR-181a-5p which has previously been associated with breast cancer metastasis are also down-regulated (p<0.001) (Fig. 4a, ,b)b) . Down-regulation of miR-30e-5p and miR-181a-5p was also observed in MDA-MB-231 cells after SOX2 knockdown (data not shown).
Several miRNAs down-regulated upon SOX2 knockdown, have been implicated in the regulation of genes involved in metabolism, signalling pathway activation (Fig. 4c, ,d).d). Intriguingly, miR-181a-5p and miR-30e-5p share a common downstream target, TUSC3 as identified by the cytoscape program. TUSC3 is a putative tumor suppressor gene, and hypermethylation or loss of this gene has been found in cancers such as ovarian and breast cancer [47–49]. Intriguingly, two potential binding sites of miR-181a-5p and miR-30e-5p were identified at the 3’-UTR of TUSC3 (Fig. 5a). Consistently, transfection of miR-181a-5p or miR-30e-5p mimics reduces the protein levels of TUSC3 in ZR7530 cells (Fig. 5b). Conversely, inhibition of these two miRNAs with specific inhibitors leads to increased protein levels of TUSC3 (Fig. 5c). We then cloned the 3’-UTR of TUSC3 containing the potential miRNA binding sites into the luciferase reporter vectors and tested whether these two miRNAs regulate TUSC3 through the potential binding sites (Fig. 5d, Additional file 2: Figure S1). Interestingly, co-transfection of miR-181a-5p or miR-30e-5p mimics reduces the luciferase activities to 47 and 55.4% of control, respectively (Fig. 5e). More importantly, mutation of these potential binding sites leads to attenuated luciferase activities, suggesting that miR-181a-5p or miR-30e-5p controls TUSC3 expression through the binding sites. As expected, SOX2 knockdown leads to increased protein levels of TUSC3 protein (Fig. 5f). Moreover, in the xenografts of SOX2 knockdown group the transcript levels of SOX2 remain low and the level of miR-181a-5p or miR-30e-5p are also significantly lower than the control group (p<0.05 for hsa-miR-30e-5p and p<0.01 for hsa-miR-181a-5p), suggesting that SOX2 is required for maintaining the expression of two microRNAs in vivo. In addition, we did not see significant change of TUSC3 transcript levels, which is consistent with the finding that the two microRNA regulate translation but not transcription of TUSC3 (Fig. 5g). Taken together, these data indicate a regulatory axis where SOX2 modulates the levels of the tumor suppressor protein TUSC3 through miR-181a-5p and miR-30e-5p.
SOX2 knockdown leads to decreased levels of miR-181a-5p and miR-30e-5p and increased levels of TUSC3. To test whether TUSC3 acts as a negative regulator of proliferation and migration of breast cancer cells, we performed TUSC3 gain-of-function study (Fig. 6a). TUSC3 overexpression leads to decreased proliferation of ZR7530 cells (p<0.01) (Fig. 6b). In addition, TUSC3 overexpression significantly inhibits the colony formation efficiency (274±6 Vs 119±4 colonies/well) in 6-well plates (p<0.01) and reduces the size of individual colonies (p<0.05) (Fig. 6c), suggesting that TUSC3 acts as a negative regulator of proliferation of breast cancer cells.
We next asked whether increased levels of TUSC3 are also able to reduce the migration of breast cancer cells. We first used a wound-healing assay to address the issue and observed that TUSC3 overexpression significantly reduced the healing index of ZR7530 cells as compared to control (24.3% Vs 12.4%) (p<0.01) (Fig. 6d). We then used the transwell assay to further test the effect of TUSC3 overexpression on cell invasion. Similar to SOX2 knockdown, TUSC3 overexpression dramatically decreased the invasion capability of ZR7530 (p<0.01) (Fig. 6e). Collectively, these results demonstrated that high levels of TUSC3 protein reduce the invasive potential of breast cancer cells.
Previous studies revealed that approximately 16% of breast cancer patients exhibit SOX2-positive expression in clinical samples, and SOX2 is expressed in the early stage of breast tumours . However, other studies revealed that SOX2 expression was associated with subtypes and tumour grades of breast cancer [5, 50]. To analyse whether Sox2 expression level (by probe 214178_s_at) correlates with distant metastasis free survival in breast cancer patients, we used the Kaplan-Meier analysis from IBM SPSS software to analyse the relationship between SOX2 protein levels and 5-year survival rates in 471 patients, and the patients data used for analysis derives from previous study . As shown in Fig. 7a, SOX2 expression level significantly converts poor prognosis (Chisquare p-value=0.014, Logrank p-value=0.02268), and the distant metastasis-free survival of breast cancer patients expressing high, middle and low levels of SOX2 protein were 50.4, 70.5 and 78.6%, suggesting that high levels of SOX2 are correlated with poor prognosis in breast cancer patients. Further analysis revealed that patients with high and middle levels of SOX2 have a significant worse survival rate than those with low levels (Fig. 7a, ,b)b) (p<0.05). However, the survival difference between the middle and low level groups is not significant (p>0.05). Moreover, we found that SOX2 levels are inversely correlated with those of TUSC3. We examined 12 pairs of normal and breast cancer biopsies, and found that 9 breast cancer biopsies have increased expression of SOX2 proteins concomitant with low levels of TUSC3 (Fig. 7c). In keeping with this line, high levels of TUSC3 protein were detected in adjacent normal tissue samples but were very low in cancer biopsies (Fig. 7d).
High levels of SOX2 proteins have been linked to breast cancer metastasis. However, the underlying mechanisms remain largely unknown. Here, we use in vitro and in vivo models to demonstrate that SOX2 is required for multiple steps in metastatic processes. We showed that SOX2 regulates the proliferation, migration and seeding of breast cancer cells. More importantly, we have identified a miRNA-mediated regulatory axis downstream of SOX2 to control these processes. Specifically, we found that SOX2 regulates the levels of miR-181a-5p and miR-30e-5p which in turn modulate the levels of the tumor suppressor TUSC3. Consistently, overexpression of TUSC3 reduces the proliferation and migration of SOX2-high breast cancer cells.
SOX2 is a key transcription factor important for the self-renewal of stem cells in multiple tissues including the trachea and brain . We and others have shown that SOX2 overexpression promotes the proliferation of stem cells in the esophagus and lung, leading to malignancy [15, 33], suggesting that SOX2 acts as an oncogene. In this study, we found that SOX2 is also required for the proliferation of breast cancer cells. SOX2 knockdown leads to decreased levels of cell cycle proteins including CCND1, CDK4 and CDK6. Interestingly, SOX2 knockdown also leads to reduced levels of multiple microRNAs. We found that the levels of miR-30e-5p (previously named as miR-30e) are most dramatically down-regulated upon SOX2 knockdown, which was validated by real-time PCR. Interestingly, we also found that miR-181a-5p (previously named as miR-181a) which has previously been linked to breast cancer development is also dramatically downregulated. More importantly, further analysis using the cytoscape program we found that miR-30e-5p together with miR-181a-5p co-regulate a common downstream target, TUSC3. Both miR-30e-5p and miR-181a-5p have been implicated in the regulation of cancer cell growth and migration, invasion and metastasis [52–54]. miR-181a-5p is a member of the miR-181 s family which includes three other highly conserved miRNAs, miR-181b, c, and d . Previous studies have shown increased levels of miR-181a-5p in multiple types of cancer including ovarian cancer and gastric cancer [56, 57]. The levels of miR-181a-5p in breast cancer were recently assessed. While Li et al. demonstrated that miR-181a-5p is down-regulated in aggressive human breast and colon cancers , a recent study showed that the levels of miR-181a-5p are high in breast cancer samples . By contrast, we have just begun to elicit the role of miR-30e-5p in tumour development. Tissue microarray analysis suggests that high levels of miR-30e-5p in primary tumours are associated with a favourable prognosis in breast cancer . In our study, we observed increased levels of miR-181a-5p and miR-30e-5p in breast cancer biopsies as compared to normal mammary gland tissue. We found that SOX2 knockdown leads to decreased levels of miR-181a-5p and miR-30e-5p accompanied by decreased cell proliferation, suggesting that these two microRNAs are required for the proliferation of breast cancer cells.
Consistently, we identified that miR-181a-5p and miR-30e-5p target the tumor suppressor TUSC3. Overexpression of these two miRNAs reduces the protein levels of TUSC3. Moreover, the luciferase assay showed that miR-181a-5p and miR-30e-5p regulate the 3’-UTR of TUSC3 transcripts. TUSC3 is a putative tumour suppressor often lost in epithelial cancers e.g. ovarian cancer and head and neck squamous cell carcinomas [47, 61]. Promoter methylation leading to silencing of TUSC3 gene was considered as a poor prognostic factor in ovarian cancer . Here, we show decreased but not complete loss of TUSC3 protein in breast cancers, suggesting that promoter methylation is unlikely a significant contributor to disease progression. We found that the protein levels of TUSC3 are decreased in breast tumour tissues as compared to surrounding normal tissues. In addition, overexpression of TUSC3 decreases the proliferation and migration capabilities of breast cancer cells. These findings are consistent with the tumor suppressor function of TUSC3 . TUSC3 is an important component of the oligosaccharyltransferase (OST) complex that catalyzes N-linked glycosylation of proteins in the endoplasmic reticulum (ER) [63, 64]. Loss of TUSC3 in mammalian cells leads to accumulation of unfolded or misfolded proteins in the ER, resulting in ER stress. In ovarian cancer cells loss of TUSC3 induces ER stress and morphological alterations in the ER. Interestingly, ovarian cancer cells lacking TUSC3 seem able to alleviate the massive ER stress response through the PERK-mediated adaptive pathway , which is consistent with the finding that accumulated unfolded proteins in cancer cells induce the unfolded protein response (UPR), facilitating cancer cells to adapt to ER stress . miR-181a-5p overexpression has been shown to promote ER stress and leads to myogenic differentiation in C2C12 cells . Along these lines our study indicates that high levels of TUSC3 expression may induce a cellular response of breast cancer cells, resulting in reduced cell proliferation and migration.
Previous studies revealed that SOX2 protein levels are correlated with the overall survival of patients with multiple types of cancer including esophageal  and lung cancer . Moreover, SOX2 protein has been assessed as a diagnostic marker and therapeutic target for multiple cancer treatments [68–70]. Our results suggest that high levels of SOX2 protein are associated with poor prognosis of breast cancer.
In summary, we showed in this study that SOX2 protein is critical for the proliferation, migration and metastasis of breast cancer. We identified that SOX2 exerts these functions through miR-181a-5p and miR-30e-5p which share a common downstream target, the putative tumour suppressor TUSC3 (Fig. 8). Knockdown SOX2 leads to decreased levels of miR-181a-5p and miR-30e-5p accompanied by increased levels of TUSC3. Consistently, overexpression of TUSC3 suppresses the proliferation and migration of breast cancer cells. Together these findings support a regulatory axis orchestrated by high levels of SOX2, and this axis is critical for breast cancer development. In future studies it will be interesting to explore potential avenues to target this SOX2/microRNA/TUSC3 axis for therapeutic gains.
The authors are grateful for inputs from colleagues in Liu and Que laboratory.
This work was supported by the National Natural Science Foundation of China (No.81302068 to K.L), the National High Technology Research and Development Program of China (863 Program, No. 2014AA020541 to K. L), and the Program for the Top Young Innovative Talents of Fujian Province (to K. L), the National Institutes of Health (R01-DK100342 to J. Q), the New York State Stem Cell Research Program (C029555 to J. Q), and the National Key Research and Development Program (2016YFA0502500 to L.Z).
The datasets supporting the conclusions of this article are included within the article, and the materials will be available once request.
KCL did the most experiments and contributed the writing of manuscript. FAX, ADG, RZ, QH, BSL conducted cell culture, western blot and qPCR experiments. KCL, JWQ and XPL analysed the data and designed the study, prepared the figures and performed the statistical analysis. LZ, ZWX, WFH, QJH, JZ and HKW planned the project and participated coordination. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
All authors read and approved the final manuscript.
We confirmed that all the protocols on human tissue examination and animal experiments were approved by the Medical Ethics Committee of Fuzhou General Hospital.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.(84K, doc) Primers used for generating TUSC3 overexpression construct. Table S2. Stem-loop reverse transcription primers used for detecting miRNAs expression. Table S3. qPCR Primers used for detecting miRNAs expression. Table S4. Primers used for generating Luciferase reporters. Table S5. Sequences of miRNA mimics and inhibitors. Table S6. microRNAs that were significantly up-regulated in ZR7530 following SOX2 knockdown as detected by microarray analysis. (DOC 83 kb) Construction of Luciferase reporters. (DOC 158 kb)
Primers used for generating TUSC3 overexpression construct. Table S2. Stem-loop reverse transcription primers used for detecting miRNAs expression. Table S3. qPCR Primers used for detecting miRNAs expression. Table S4. Primers used for generating Luciferase reporters. Table S5. Sequences of miRNA mimics and inhibitors. Table S6. microRNAs that were significantly up-regulated in ZR7530 following SOX2 knockdown as detected by microarray analysis. (DOC 83 kb)Additional file 2: Figure S1.(159K, doc)
Construction of Luciferase reporters. (DOC 158 kb)
Kuancan Liu, Phone: +86-591-22859490, Email: moc.361@nacnaukuiL.
Jianwen Que, Phone: +1-212-305-5961, Email: ude.aibmuloc.cmuc@0422qj.
Xiaopeng Lan, Phone: +86-591-22859490, Email: moc.anis@pxnaL.