Research | Open | Published:
BRD7 expression and c-Myc activation forms a double-negative feedback loop that controls the cell proliferation and tumor growth of nasopharyngeal carcinoma by targeting oncogenic miR-141
Journal of Experimental & Clinical Cancer Researchvolume 37, Article number: 64 (2018)
miR-141 is up-regulated and plays crucial roles in nasopharyngeal carcinoma (NPC). However, the molecular mechanism underlying the dysregulation of miR-141 is still obscure.
Thus, the ChIP-PCR was performed to identify the c-Myc-binding sites in miR-141 and BRD7. qRT-PCR, western blot and immunohistochemistry assays were used to detect the expression of miR-141 and its up/down stream molecules. The rescue experiments on the c-Myc/miR-141 axis were performed in vitro and in vivo.
Our results showed that the levels of mature miR-141, pre-miR-141 and pri-miR-141 were downregulated in c-Myc knockdown NPC cells. Meanwhile, c-Myc transactivates the expression of miR-141 by binding its promoter region. Moreover, BRD7 was identified as a co-factor of c-Myc to negatively regulate the activation of c-Myc/miR-141 axis, as well as a direct target of c-Myc. Moreover, restoration of miR-141 in c-Myc knockdown NPC cells notably rescued the effect of c-Myc on cell proliferation and tumor growth, as well as the blocking of PTEN/AKT pathway. Additionally, the expression of c-Myc was positively correlated with that of miR-141 and the clinical stages of NPC patients and negatively associated with the expression of BRD7. Our findings demonstrated that BRD7 expression and c-Myc activation forms a negative feedback loop to control the cell proliferation and tumor growth by targeting miR-141.
These observations provide new mechanistic insights into the dysregulation of miR-141 expression and a promising therapeutic option for NPC.
Nasopharyngeal carcinoma (NPC) is a leading form of head and neck cancer, especially in South China. It has high metastatic potential that contributes to a high rate of local invasion and early metastasis . Thus, there is a clear need to identify some sensitive biomarkers and novel therapeutic targets for diagnosis and treatment options of NPC in a premalignant stage.
Recently, non-coding RNAs have been demonstrated to represent most of the human transcriptome and are frequently dysregulated in disease pathogenesis, including cancer . MicroRNAs (miRNAs) are a class of functional, short non-coding RNAs that are associated with all hallmarks in cancer initiation, progression and metastasis [3, 4]. Therefore, considering their tissue-specific profiles and stability, miRNA-based strategies are considered as promising options for early detection, accurate diagnosis and the prediction of responses to treatment in malignancies. miR-141 is a member of the miR-200 family, and it has been shown to be dysregulated in a wide variety of cancers [5, 6]. However, the mechanisms underlying the dysregulation of miR-141 are only beginning to be understood in cancer.
c-Myc is a transcription factor that generally heterodimerizes with other basic-helix–loop–helix (bHLH) proteins at thousands of genomic loci , which affects more than 15% of the human transcriptome . Hence, c-Myc plays vital roles in almost every aspect of basic cellular processes. Aberrant activation of c-Myc is known to be a crucial hallmark of numerous cancers  and often leads to widespread dysregulation of miRNAs . In our previous study, we found that knockdown of c-Myc significantly inhibited cell proliferation and tumor growth of NPC  and downregulated the level of miR-141 in NPC cells . BRD7, a member of the bromodomain-containing protein family, was identified as a critical tumor suppressor in multiple types of cancers, including NPC and breast cancer [13,14,15] and also involved in many physiological processes [16,17,18]. Growing evidence has revealed that BRD7, as a ubiquitously expressed nucleic protein [17, 19], is generally involved in the transcription of specific target genes through a protein-protein interaction (PPI) with other transcription factors [20,21,22]. Moreover, BRD7 could negatively regulate the transcription of miR-141 in an indirect manner in NPC . However, the detail mechanism underlying BRD7 and c-Myc involvement in the regulation of miR-141 is still obscure and remains to be further investigated.
In this study, we estimated the expression of BRD7 and miR-141 in c-Myc knockdown NPC cells and explored the mechanism of BRD7 and c-Myc in the transcription of miR-141. Furthermore, we performed rescue experiments to identify the functional role and molecular mechanism of the c-Myc/miR-141 axis in NPC in vitro and in vivo. Additionally, we investigated the clinical significance of c-Myc and the association of c-Myc expression with BRD7 and miR-141 in NPC patients. In summary, our findings demonstrated that BRD7 expression and c-Myc activation forms a negative feedback loop that regulates miR-141 transcription and controls cell growth and proliferation in NPC. These observations provide new mechanistic insights into the dysregulation of miR-141 expression and a promising therapeutic option for NPC in early stages.
Cell lines and tissue samples
The human NPC cell line 5-8F was kindly provided by the Cancer Center of Sun Yet-Sen University (Guangzhou, China). The HNE1 cell line was purchased from the cell center of the Xiangya School of Medicine (Central South University, Changsha, China). Both cell lines were preserved in our laboratory and cultured in RPMI-1640 (Hyclone, Logan, UT, USA) supplemented with 10% fetal bovine serum (FBS) (Gibco-BRL, Invitrogen, Paisley, UK) in a humidified incubator with 5% CO2 at 37 °C. The embryonic kidney cell line HEK293 were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA) and maintained in Dulbecco’s modified Eagle’s medium (DMEM) (Hyclone) with 10% FBS (Gibco) at 37 °C in a humidified atmosphere with 5% CO2.
NPC samples (n = 41) and non-cancerous nasopharyngeal tissues (n = 27) from healthy donors were collected at the Second Xiangya Hospital of Central South University (Changsha, China). The non-cancerous nasopharyngeal tissues were collected from independent patients with chronic inflammation of nasopharyngeal mucosa or polyps. The profile of clinicopathological characteristics of the NPC patients is shown in Table 1. All tissue samples were immediately snap-frozen in liquid nitrogen and stored in a freezer at − 80 °C. Tissues for immunohistochemistry were fixed in 4% paraformaldehyde and then paraffin embedded. Clinicopathological data were reviewed, and TNM staging classification was based on the criteria of the American Joint Committee on Cancer (AJCC, 6th edition).
Plasmids, oligonucleotides and transfections
The interfering plasmid (pRNAT-U6.1/shc-Myc) expressing short hairpin RNA (shRNA) for c-Myc knockdown, the expression plasmid encoding the full-length open reading frame of human c-Myc with 2Flag tags (pIRES2-EGFP/2Flag-c-Myc), and the plasmid encoding the full-length open reading frame of human BRD7 with 3Flag tags (pIRESneo3/3Flag-BRD7) or GFP tag (pEGFP-C2/BRD7) were constructed by our laboratory. The hsa-miR-141 mimic (miR-141) and negative control (miR-NC) were purchased from RiBoBio (Guangdong, China).
Cells were transfected with plasmids and/or oligonucleotides using Lipofectamine3000 Reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. Briefly, cells in logarithmic growth phase were trypsinized, counted and seeded in a 6-well plate to ensure 60–80% cell confluence on the next day for transfection. Transfection of cells with oligonucleotides was performed at a final concentration of 50 nM.
The c-Myc knockdown stable cell lines (5-8F/shc-Myc and HNE1/shc-Myc) and the control cell lines (5-8F/shCtrl and HNE1/shCtrl) were transfected with pRNAT-U6.1/shc-Myc and pRNAT-U6.1/shCtrl, respectively, and then screened under G418 (Invitrogen). Finally, stable pool clones were obtained.
Total RNA was isolated from tissue samples and cell lines using TRIzol reagent (Invitrogen) according to the manufacturer’s protocol. The levels of mature miR-141 (miR-141) and precursor miR-141 (pre-miR-141) were evaluated using the miDETECT A Track miRNA qRT-PCR Kit (RiBoBio). The U6 small nuclear RNA (RNU6B) (RiBoBio) was used for normalization. The expression of primary miR-141 transcript (pri-miR-141), c-Myc and BRD7 was measured by qRT-PCR according to the instructions of the SYBR Premix Ex Taq (TaKaRa, Dalian, China). The GAPDH mRNA level was used for normalization. The relative expression ratio was calculated using the 2−ΔΔCT method. The primers for miR-141, pre-miR-141, U6, pri-miR-141, c-Myc, BRD7, and GAPDH were described previously . PCRs of each sample were conducted at least in triplicate.
Western blot analysis
Cells and tissues were lysed in RIPA buffer in the presence of Protease Inhibitor Cocktail and PhoSTOP (Roche, Basel, Switzerland). Protein was quantified using a BCA Protein Assay Kit (Pierce Biotechnology, Rockford, IL, USA). Protein (30–80 μg) was separated by 8–12% sodium dodecyl sulfate polyacrylamide gel electrophoresis, and transferred onto polyvinylidene fluoride (PVDF) membranes (Millipore, Billerica, MA, USA). The membranes were blocked with 5% non-fat milk in Tris-buffered saline and then incubated with primary antibodies at 4 °C overnight. The primary antibodies used were anti-c-Myc (dilution 1:1000; CST, Danvers, MA, USA), anti-Dicer (dilution 1:1000; CST), anti-Drosha (dilution 1:1000; CST), anti-BRD7 (dilution 1:500; ProteinTech, Wuhan, China), anti-Flag (dilution 1:2000; Sigma, St. Louis, MO, USA), anti-GFP (dilution 1:2000; ProteinTech), anti-β-actin (dilution 1:1000; Santa Cruz, CA, USA), anti-PARP (dilution 1:1000; CST), anti-cleaved-PRAP (C-PARP) (dilution 1:1000; CST), anti-CCND1 (dilution 1:500; Santa Cruz), anti-CDK4 (dilution 1:500; Santa Cruz), anti-PTEN (dilution 1:500; Bioworld Technology, Atlanta, GA, USA), anti-pS473AKT (dilution 1:500; Bioworld Technology), anti-AKT (dilution 1:500; Bioworld Technology), anti-p27 (dilution 1:1000; CST), anti-Caspase9 (dilution 1:500; Santa Cruz) and anti-cleaved-Caspase9 (dilution 1:500; Santa Cruz). Membranes were then washed three times in TBST solution for 10 min each time and incubated with secondary antibodies. Signals were detected with an enhanced chemiluminescence detection system (Bio-Rad, Hercules, CA, USA).
The luciferase reporter assays
The recombinant reporter vectors of the miR-141 promoter (pGL3/miR-141P) and BRD7 promoter (pGL3/BRD7P) were constructed by our laboratory as described previously [23, 24]. The cells were seeded in 24-well plates. After 24 h, the cells were transfected with recombinant reporter vector or pGL3-Enhancer empty vector (Promega, Madison, WI, USA), together with the pRL-TK vector (Promega) containing the Renilla luciferase gene. Transfection was performed using Lipofectamine 3000 (Invitrogen). Cells were harvested at 36 h post-transfection. Firefly and Renilla luciferase activities were measured using a Dual-luciferase reporter kit (Promega) according to the manufacturer’s protocol. Firefly luciferase activity was normalized to Renilla luciferase activity.
Chromatin immunoprecipitation (ChIP) assays
HEK293 cells transfected with pIRES2-EGFP/2Flag-c-Myc plasmid were cross-linked in 1% formaldehyde for 10 min at 37 °C. DNA from fixed chromatin cells was then subjected to immunoprecipitation using a ChIP assay kit (Millipore) and antibodies against Flag tags or mouse IgG according to the manufacturer’s protocol. The purified DNA was used for the following PCRs. The primers for confirmation of the c-Myc-binding sites in the miR-141 and BRD7 promoter in PCR reactions are described in Additional file 1: Table S1.
Immunofluorescence (IF) and co-immunoprecipitation (co-IP)
Subcellular localization of endogenous c-Myc and exogenous BRD7 was analyzed by IF. Briefly, 5-8F and HNE1 cells were seeded on glass coverslips, transfected with pEGFP-C2/BRD7, fixed with 4% paraformaldehyde for 15 min at room temperature and washed with TBS. Cells were then permeabilized with 0.3% Triton X-100 for 5 min at room temperature. Cells were washed 3 times with TBS and incubated for 15 min in 5% goat serum and then incubated at 4 °C overnight with the primary antibody anti-c-Myc (dilution 1:800; CST). After three 5-min washes with TBS, cells were incubated for 1 h at 37 °C in the dark with the secondary antibody Alexa Fluor 568 Donkey anti-Rabbit IgG (H + L) (dilution 1:1000; Invitrogen). Coverslips were washed with TBS and visualized using an epifluorescence microscope.
For Co-IP, HEK293 cells transfected with pIRES2-EGFP/2Flag-c-Myc plasmid and pEGFP-BRD7 plasmid together were lysed with cell lysis buffer for western blotting and IP (Beyotime, Beijing, China). They were then pre-mixed with protease inhibitor cocktail (Roche) on ice for 20 min and centrifuged at 12,000 g for 15 min. The supernatants were incubated at 4 °C overnight with 30 μl of protein A/G magnetic beads (Biomake, Houston, TX, USA) pre-coated with anti-Flag antibodies (Sigma) or anti-mouse IgG (Santa Cruz). The immunocomplexes were isolated and purified according to the manufacturer’s protocol, and then were subjected to western blot analysis.
Cell proliferation (MTS) and colony-forming assays
The cell proliferation capacity was estimated with the MTS kit (Biomake) following the manufacturer’s protocol. Briefly, cells were seeded at a density of 1 × 103 cells/well in a 96-well plate and cultured for 48 or 96 h. Ten microliters of MTS solution was added to each well, and the samples were incubated at 37 °C for 3 h before the absorbance was measured at 450 nm.
For colony-forming assays, approximately 150 cells per well were added to a 6-well plate, with three wells per sample. After a 14-day incubation, the cells were washed twice with PBS and stained with Giemsa solution (Beyotime, Beijing, China). The plate colony-forming efficiency was calculated as (number of colonies/number of cells inoculated) × 100%. All experiments were performed at least in triplicate.
Flow cytometry for cell cycle and apoptosis analysis
For cell cycle analysis, cells were harvested and fixed in 70% ethanol for 24 h at − 20 °C. The cells were then treated with RNase A and stained with 25 μg/ml propidium iodide (PI). For cell apoptosis analysis, cells were treated with serum starvation for 24 h, and then the ratio of apoptotic cells was determined using an Annexin V-PE/7-AAD double staining kit (BD Biosciences, MD, USA). Samples were all analyzed using a MoFlo™ XDP High-Performance Cell Sorter (Beckman Coulter, CA, USA), and the data were analyzed using Summit v.5.2 software according to the manufacturer’s protocol. At least three independent experiments were performed.
Immunohistochemical (IHC) staining, in situ hybridization (ISH) and scores
For IHC staining, the sections were incubated with anti-c-Myc (dilution 1:100; Millipore), anti-PTEN (dilution 1:100; Bioworld), anti-pS473AKT (dilution 1:100; Bioworld), anti-p27 (dilution 1:100; CST), and anti-CCND1 (dilution 1:100; Santa Cruz) antibodies. The immune complex was visualized by the MaxVision HRP-polymer IHC Kit Detection System, Peroxidase/DAB, Rabbit/Mouse (MaxVision, Fuzhou, China) according to the manufacturer’s protocol. The nuclei were counterstained with hematoxylin (Beyotime). For ISH staining, the mature Hsa-miR-141-specific probe and negative control (Scramble) were purchased from Sangon (Sangon Biotech, Shanghai, China). An in situ hybridization kit (Boster, Wuhan, China) was used for hybridization of the probe according to the manufacturer’s instructions.
Scores for IHC and ISH analysis. IHC and ISH staining were independently evaluated at 200× magnification using light microscopy by two pathologists who were blinded to the clinicopathological data. A semiquantitative evaluation of c-Myc, BRD7 protein and mature miR-141 was performed using a method described in our previous work .
Tumor growth xenograft model
A total volume of 200 μl of 5-8F/shc-Myc and 5-8F/shCtrl cells (1× 107 cells) transfected with miR-141 mimic (miR-141) or negative control (miR-NC) was inoculated subcutaneously into the left flanks of 6-week-old female nude mice. Mice were checked every 4 days. Tumor volume was evaluated using the following formula: volume = (width + length)/2 × width × length × 0.5236. All three groups were killed after 28 days. All tumor grafts were excised, weighed, harvested, fixed and embedded. The anti-Ki67 antibody (dilution 1:100, Bioworld) was used to detect the proliferation marker Ki67 using IHC procedures, and a TUNEL assay was used to detect apoptosis in situ with paraffin-embedded xenograft sections using the DeadEnd Colorimetric TUNEL Detection Kit (Promega) according to the manufacturer’s protocols. Samples were observed using an Olympus microscope (Olympus, Tokyo, Japan). The proliferative and apoptosis index scores were measured as the mean percentage of nuclei that stained positive for Ki67 and TUNEL cells in 10 different 200× fields.
The relationships between the expression of c-Myc, BRD7, miR-141 and clinicopathological characteristics in NPC were tested using a chi-squared test. The Spearman’s rank correlation coefficient was used to assess the significance of the association among c-Myc, BRD7 and miR-141 expression in NPC patients. The differences between the groups were analyzed using a Student’s t-test when there were only two groups or using a one-way ANOVA when there were more than two groups. All statistical analyses were performed using SPSS software (SPSS, Chicago, IL, USA). A two-tailed value of P < 0.05 was considered statistically significant.
C-Myc transactivates miR-141 expression
We first evaluated the level of miR-141 in c-Myc knockdown stable NPC cells (Fig. 1a). The results showed that the level of miR-141 was significantly downregulated in both c-Myc knockdown 5-8F and HNE1 cells (Fig. 1b). Considering the general principles of miRNA biosynthesis, we further examined the abundance of pri-miR-141 and pre-miR-141 and the expression of Dicer and Drosha, two key enzymes in the miRNA maturation. As a result, the levels of both pri-miR-141 and pre-miR-141 were correlated with the miR-141 level in c-Myc knockdown cells (Fig. 1c and Additional file 2: Figure S1a). Nevertheless, knocking down c-Myc had no obvious effect on the expression of Dicer and Drosha protein (Additional file 2: Figure S1b). Conversely, in HEK293 cells with endogenous low-level expression of c-Myc, overexpression of c-Myc significantly elevated the abundance of miR-141, pre-miR-141 and pri-miR-141 (Fig. 1d-e and Additional file 2: Figure S1c) but did not obviously influence the Dicer and Drosha protein expression (Additional file 2: Figure S1d).
Considering that c-Myc is a transcription factor and that its aberrant activation often leads to widespread dysregulation of miRNAs, we further detected the effect of c-Myc on the promoter activity of miR-141. The results showed that the miR-141 promoter activity was significantly reduced in both c-Myc knockdown cells (Fig. 1f). Meanwhile, the miR-141 promoter activity was significantly increased in c-Myc-overexpressing HEK293 cells (Fig. 1g). Furthermore, we mapped the c-Myc-binding profile using the JASPAR Database . Within the miR-141 promoter region, 10-putative c-Myc-binding sites were present (Additional file 3: Table S2). Considering the location of the miR-141-encoding gene on chromosome 12 forward strand, we conducted ChIP-PCR assays and tested the 5-putative c-Myc-binding sites (named BS1, BS2, BS3, BS4 and BS5) (Fig. 1h). The assays validated that c-Myc binds directly to the region covering the BS3 and BS4 sites in the miR-141 promoter (Fig. 1i). Altogether, these data demonstrated that c-Myc activates miR-141 expression by directly binding to its promoter region.
A negative feedback loop formed by BRD7 and c-Myc is involved in the regulation of miR-141 expression
c-Myc plays a dual role in the transcription, activation or repression of its target genes, and this transition practically depends on its partner protein . Likewise, BRD7 often regulates specific target genes through interaction with other transcription factors [20,21,22]. Additionally, we have identified that BRD7, as a transcriptional regulator, could repress the transcription of miR-141 in an indirect manner . Therefore, we want to know whether BRD7 could serve as a cofactor of c-Myc to regulate the expression of miR-141. Firstly, we predicted that BRD7 might be one of c-Myc interacting proteins by interrogating the annotated PPI predictions data from the Fpclass database (Additional file 4: Table S3). Then, the IF analysis showed that BRD7 was colocalized with c-Myc in a punctate distribution pattern (Fig. 2a), and the Co-IP experiments further revealed that BRD7 was co-immunoprecipitated with c-Myc (Fig. 2b). All of these data supported that BRD7 can indeed serve as an interactor of c-Myc. Furthermore, we detected the effect of BRD7 on the transcription activity and expression of miR-141 regulated by c-Myc. As a result, ectopic expression of BRD7 significantly repressed the c-Myc-induced activation of the miR-141 promoter (Fig. 2c and Additional file 2: Figure S2a) as well as the expression of miR-141 (Fig. 2d and Additional file 2: Figure S 2b). Additionally, we examined that ectopic expression of BRD7 did not affect the expression of c-Myc protein and mRNA (Additional file 2: Figure S2c and d). Importantly, we have previously performed an integrative analysis of genome-wide chromatin occupancy of BRD7 by ChIP-sequencing  and noted that there were 8 binding peaks of BRD7 on chromosome 12 but no more than one in the miR-141 promoter region (Fig. 2e). These observations demonstrated that BRD7 could serve as a cofactor of c-Myc to be negatively involved in the transcription of miR-141 regulated by c-Myc.
Since BRD7, as co-factor of c-Myc, is negatively involved in the transcription of miR-141, we still want to know whether c-Myc feedbackly affect the expression of BRD7. As a result, knockdown of c-Myc observably enhanced the abundance of BRD7 protein and mRNA as well as the BRD7 promoter activity in both 5-8F and HNE1 cells (Fig. 2f-h and Additional file 2: Figure S2e-g). On the contrary, overexpression of c-Myc could remarkably downregulate the expression of BRD7 protein and mRNA as well as the BRD7 promoter activity in HEK293 cells (Fig. 2i, j and k). However, overexpression of c-Myc could not change the expression and the promoter activity of BRD7 in 5-8F and HNE1 cells (Additional file 2: Figure S2 h, i and j), suggesting that the high level of endogenous c-Myc is sufficient to inactivate the transcription of BRD7 in NPC cells. Beyond that, the ChIP-PCR experiments further revealed that c-Myc could directly bind to the BRD7 promoter region (Fig. 2l). These results suggest that BRD7 expression and c-Myc activation form a negative feedback loop responsible for the regulation of miR-141 transcription.
Restoring the expression of miR-141 significantly rescues the tumor suppressive effects of c-Myc knockdown in NPC cells
To further investigate the role of c-Myc transactivating miR-141 expression, we performed rescue experiments in c-Myc knockdown NPC cells. The level of miR-141 restoration was confirmed by qRT-PCR assays (Additional file 2: Figure S3). The MTS assays and colony-forming experiments showed that restoring the expression of miR-141 significantly reversed cell proliferation inhibition by c-Myc knockdown (Fig. 3a and b). Meanwhile, the cell cycle analysis further validated that the restoration of miR-141 levels elicited a notable reduction of the G0/G1 population compared to the c-Myc knockdown controls (Fig. 3c and d). Furthermore, restoring the expression of miR-141 significantly recovered the percentage of apoptotic cells (Fig. 3e) and the expression of cleaved-PARP (C-PARP) (Fig. 3f) compared to the c-Myc knockdown controls. These results together demonstrated that the c-Myc/miR-141 axis has a critical role in cell proliferation and apoptotic responses, suggesting that c-Myc exerted its oncogenic through transactivating miR-141 expression in NPC progression.
miR-141 can reverse c-Myc knockdown-mediated tumor growth suppression in vivo
To corroborate our in vitro observations, we conducted in vivo experiments with xenograft tumor models in nude mice. As a result, the growth rate and weight of the xenograft tumors with miR-141 restoration were significantly increased compared with those of the c-Myc knockdown group (Fig. 4a and b). In addition, the expression of c-Myc and miR-141 in xenografts was confirmed by western blotting and qRT-PCR assays, respectively (Fig. 4c and d). To further validate the effect of miR-141 restoration on tumor growth in vivo, we determined the percentage of cells expressing the proliferation marker Ki67 and the percentage of TUNEL-positive cells. Consistent with our in vitro data, the restoration of miR-141 caused a significant increase in the number of Ki67+ cells and a reduction in TUNEL+ cells in c-Myc knockdown xenograft tumors when compared with the c-Myc knockdown controls (Fig. 4e and f). All of these in vivo findings corresponded with our in vitro results and supported that c-Myc, as an oncogenic transcription factor, promoted tumor growth through transactivating miR-141 expression in NPC.
C-Myc promotes NPC progression through the miR-141/PTEN/AKT pathway
We have previously demonstrated that miR-141 could directly target the 3′ untranslated region (3’UTR) of PTEN , and the restoration of miR-141 could significantly reverse the tumor suppressive effect of BRD7 on tumor growth through the PTEN/AKT pathway in NPC . Thus, we reasoned that c-Myc can promote NPC growth through the PTEN/AKT pathway by activating miR-141 expression. First, we confirmed the restoration of miR-141 levels in c-Myc knockdown NPC cells (Fig. 5a). Then, we detected the expression of PTEN, AKT and pAKT. The immunoblot showed that knockdown of c-Myc in either 5-8F or HNE1 cells led to an increase in the expression of PTEN protein and induced a notable dephosphorylation of pAKT compared with the controls (Fig. 5b). Conversely, restoring miR-141 levels reduced the protein expression of PTEN and stimulated AKT phosphorylation in c-Myc knockdown NPC cells (Fig. 5b), suggesting that the activation of the AKT pathway by c-Myc is at least partially dependent on the repression of PTEN by miR-141 in NPC cells. Notably, IHC staining of xenografts showed that knockdown of c-Myc induced the expression of PTEN and decreased the phosphorylation level of pAKT compared with the control group (Fig. 5c). Collectively, these data illustrated that c-Myc activates AKT signals by downregulating the expression of PTEN protein at least partially via transcriptionally modulating miR-141 expression in vitro and in vivo.
To identify the effectors responsible for the phenotypes observed in vitro and in vivo, we examined alterations of the key components downstream in the PTEN/AKT pathway. The immunoblot results revealed that the restoration of miR-141 in c-Myc knockdown cells resulted in a reduction of p27 protein and an increase of CCND1, while CDK4 expression remained unaltered compared to the controls (Fig. 5b). Consistent with the in vitro findings, IHC staining of xenografts further confirmed that the expression of p27 was decreased and that CCND1 was increased in miR-141-restoring tumors when compared to the c-Myc knockdown controls (Fig. 5c). Additionally, Caspase9 cleavage was blocked when miR-141 was restored in c-Myc knockdown tumors (Fig. 5c). Collectively, these data further validated that the c-Myc/miR-141 axis controls NPC growth at least partially through modulating the PTEN/AKT pathway.
Clinical significance of the feedback loop in miR-141 transcription in human NPC
To assess the significance of the feedback loop between BRD7 and c-Myc in NPC progression, we detected their expression in NPC tissues and non-cancerous NP samples (Table 1, Fig. 6a). As a result, most patients exhibited significantly enhanced expression of c-Myc (P < 0.01) (Fig. 6a and b). Importantly, the high-level expression of c-Myc was positively correlated with the clinical stages of NPC patients (P = 0.014) (Table 1), indicating that c-Myc, as an oncogenic transcription factor, plays crucial roles in NPC progression.
Consistent with our previous data, the expression of BRD7 protein was significantly decreased (P < 0.001), and on the contrary, the level of miR-141 was overexpressed in NPC tissues compared with the normal controls (P < 0.001) (Fig. 6a and b). Considering the mechanism of BRD7 and c-Myc in miR-141 transcription, we further explored the association of c-Myc, BRD7 and miR-141 expression in NPC patients. The observations validated that the expression of c-Myc was significantly positively correlated with the level of miR-141 (P = 0.008) and negatively associated with the expression of BRD7 (P = 0.001) in NPC patients (Table 2). More importantly, we have previously demonstrated that BRD7 or miR-141 could serve as an independent prognostic factor of NPC . In total, these results suggested that the feedback loop formed by BRD7 and c-Myc in miR-141 transcription is of great importance in NPC progression (Fig. 6c) and could serve as potential biomarkers for early detection and prognosis of NPC.
MiR-141 is located on chromosome 12p13.31, and as a member of the miR-200 family, it is dysregulated in a wide variety of cancers [28, 29]. The genetic basis for aberrant miRNA expression is complicated in cancers. In FANTOM5 project, an integrated expression atlas of miRNAs and their promoters has been recently created and revealed that transcription factors are broadly involved in the regulation of specific miRNAs . For example, TP53 regulated the miR-15/16 cluster , c-Myc modulated the miR-17-92 cluster , and HIF-1alpha targeted miR-210 . In our previous study, we documented that BRD7 could negatively regulate miR-141 expression in an indirect manner at the transcription level . Herein, we found that c-Myc knockdown decreased the levels of pri-miR-141 and pre-miR-141 in NPC cells, and downregulation was correlated with the expression of miR-141, which is consistent with our previous results derived from the differentially expressed miRNA profile in c-Myc knockdown NPC cells . Moreover, we demonstrated that c-Myc bound directly to the miR-141 promoter and activated its transcription. In addition, c-Myc knockdown in NPC cells did not influence the expression of Dicer and Drosha proteins, two key enzymes that are indispensable to miRNA biogenesis and maturation. All of these observations together suggest that miR-141 is tightly regulated by transcription factors through a proximal promoter and that aberrant activation of c-Myc is at least partially responsible for the dysregulation of miR-141 in NPC.
c-Myc is a member of the proto-oncogenic bHLH transcription factor family. Numerous studies have identified that no monomeric c-Myc protein is present in vivo, and instead, c-Myc usually interacts with a partner protein, such as Max, Mad1, Mxi1 and Mnt, which is essential for the dimer to directly bind to the general E-box [26, 34, 35]. Indeed, c-Myc plays a dual role in the transcription activation or repression of its target genes, and this transition practically depends on the interaction partner of c-Myc . For example, p19ARF inhibits cell proliferation in a p53 stabilization-dependent manner and also functions in a p53-independent way by binding to c-Myc [36, 37]. BRD7, as a ubiquitous nuclear factor, is frequently involved in transcriptional regulation through protein-protein interactions with other transcription factors, such as TP53, BRCA1 and Smads [20, 22, 38]. In this study, we report that BRD7 colocalized with c-Myc in a punctate distribution pattern and interacted with c-Myc. Meanwhile, overexpression of BRD7 did not affect the mRNA level and protein stability of c-Myc. Importantly, we have previously performed an integrative analysis of genome-wide chromatin occupancy of BRD7 by ChIP-sequencing and digital gene expression profiling by RNA-sequencing and identified that there was no more than one binding-site of BRD7 within the miR-141 promoter region . Interestingly, BRD7, as a direct target of c-Myc, could negatively regulate the transcription of miR-141 in an indirect manner in NPC [23, 24]. More importantly, we now document BRD7 being able to impair the transcriptional activation of miR-141 by c-Myc in NPC cells. Consequently, in light of the results observed above, it is apparent that the activity of miR-141 transcription was enhanced via BRD7 inactivation by c-Myc, and BRD7 could also serve as a cofactor of c-Myc to negatively regulate the activity of c-Myc/miR-141 transcription axis, suggesting that BRD7 and c-Myc might form a negative feedback loop responsible for the regulation of miR-141 expression.
Aberrantly elevated expression of c-Myc was highly prevalent in a variety of cancers and has been widely demonstrated as contributing to almost every aspect of tumor cell biology . Several studies have shown that in c-Myc-driven malignancies, blocking continuous activation of c-Myc leads to growth inhibition and apoptosis . Here, we confirm that knockdown of c-Myc causes NPC cell proliferation repression and tumor growth arrest. Additionally, the restoration of miR-141 can significantly rescue the tumor suppressive effect of c-Myc knockdown in NPC cells in vitro and in vivo. PTEN is a well-known lipid phosphatase and bona fide tumor suppressor, and numerous studies have identified a variety of cancers that arise following loss of PTEN function, which finally leads to a potent derepression of the PI3K/AKT pathway that stimulates cell proliferation and survival . We have described previously that BRD7 has potent anti-proliferative and tumor suppressor activity at least partially through transcriptionally repressing miR-141, which directly targets the 3’UTR of PTEN mRNA and subsequently contributes to AKT pathway activation in NPC [12, 23]. Here, we document that knockdown of c-Myc activates PTEN/AKT signaling via transcriptional downregulation of miR-141 in NPC cells. Numerous observations indicated that following the downregulation of c-Myc, p27 is activated and cyclin D/CDK activity is repressed, resulting in cell proliferation inhibition by arresting G0/G1 phase progression [42, 43]. In accordance with previous identification, we validated that c-Myc knockdown induces G0/G1 phase arrest in NPC cells and that miR-141 can significantly reverse the tumor suppressive effects of c-Myc knockdown partially through restoring the expression of cyclin D1 and p27 in vitro and in vivo via the PTEN/AKT pathway. All of these results are consistent with the previous realization that minute variations of PTEN protein possess potent consequences for tumor initiation and cancer susceptibility  and suggest that the c-Myc/miR-141 axis regulates NPC growth at least partially through modulating the PTEN/AKT pathway.
We have previously reported that miR-141 was significantly upregulated while BRD7 was downregulated in NPC tissues, and both of them could serve as an independent prognostic factor of NPC patients . In this work, we report that the expression of c-Myc is significantly upregulated in NPC patients, and the high-level expression of c-Myc is significantly correlated with advanced clinical stages of NPC. More importantly, the expression of c-Myc is positively associated with the level of miR-141 and negatively associated with the expression of BRD7 in NPC patients. Therefore, it is tempting to speculate that abnormalities in the feedback loop between BRD7 and c-Myc might be particularly required for NPC initiation and progression.
Our findings demonstrated that c-Myc transactivates miR-141 expression and that BRD7, a direct target of c-Myc, serves as a cofactor of c-Myc and forms a negative feedback loop with c-Myc in miR-141 transcription in NPC. Moreover, the c-Myc/miR-141 axis promotes NPC growth through modulating the PTEN/AKT pathway. Furthermore, the feedback loop in miR-141 transcription is of great significance for early detection and prognosis of NPC. All of these observations provide new mechanistic insights into the dysregulation of miR-141 expression and a promising therapeutic option to interrupt the aberrant activation of miR-141 contributing to NPC progression.
3′ untranslated region
Functional Annotation of Mammalian Genome
In situ hybridization
Primary miR-141 transcript
Quantitative Real-Time polymerase chain reaction
Short hairpin RNA
Chua MLK, Wee JTS, Hui EP, Chan ATC. Nasopharyngeal carcinoma. Lancet. 2016;387:1012–24.
Beermann J, Piccoli MT, Viereck J, Thum T. Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches. Physiol Rev. 2016;96:1297–325.
Hammond SM. An overview of microRNAs. Adv Drug Deliv Rev. 2015;87:3–14.
Yang R, Liu M, Liang H, Guo S, Guo X, Yuan M, Lian H, Yan X, Zhang S, Chen X, Fang F, Guo H, Zhang C. miR-138-5p contributes to cell proliferation and invasion by targeting Survivin in bladder cancer cells. Mol Cancer. 2016;15:82.
Yang D, Sun Y, Hu L, Zheng H, Ji P, Pecot CV, Zhao Y, Reynolds S, Cheng H, Rupaimoole R, Cogdell D, Nykter M, Broaddus R, Rodriguez-Aguayo C, Lopez-Berestein G, Liu J, Shmulevich I, Sood AK, Chen K, Zhang W. Integrated analyses identify a master microRNA regulatory network for the mesenchymal subtype in serous ovarian cancer. Cancer Cell. 2013;23:186–99.
Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, Vadas MA, Khew-Goodall Y, Goodall GJ. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol. 2008;10:593–601.
Kress TR, Sabo A, Amati B. MYC: connecting selective transcriptional control to global RNA production. Nat Rev Cancer. 2015;15:593–607.
Eilers M, Eisenman RN. Myc's broad reach. Genes Dev. 2008;22:2755–66.
Gonzalez V, Hurley LH. The c-MYC NHE III(1): function and regulation. Annu Rev Pharmacol Toxicol. 2010;50:111–29.
Chang TC, Yu D, Lee YS, Wentzel EA, Arking DE, West KM, Dang CV, Thomas-Tikhonenko A, Mendell JT. Widespread microRNA repression by Myc contributes to tumorigenesis. Nat Genet. 2008;40:43–50.
Niu Z, Liu H, Zhou M, Wang H, Liu Y, Li X, Xiong W, Ma J, Li X, Li G. Knockdown of c-Myc inhibits cell proliferation by negatively regulating the Cdk/Rb/E2F pathway in nasopharyngeal carcinoma cells. Acta Biochim Biophys Sin. 2015;47:183–91.
Zhang L, Deng T, Li X, Liu H, Zhou H, Ma J, Wu M, Zhou M, Shen S, Li X, Niu Z, Zhang W, Shi L, Xiang B, Lu J, Wang L, Li D, Tang H, Li G. microRNA-141 is involved in a nasopharyngeal carcinoma-related genes network. Carcinogenesis. 2010;31:559–66.
Zhou M, Liu H, Xu X, Zhou H, Li X, Peng C, Shen S, Xiong W, Ma J, Zeng Z, Fang S, Nie X, Yang Y, Zhou J, Xiang J, Cao L, Peng S, Li S, Li G. Identification of nuclear localization signal that governs nuclear import of BRD7 and its essential roles in inhibiting cell cycle progression. J Cell Biochem. 2006;98:920–30.
Zhou J, Ma J, Zhang BC, Li XL, Shen SR, Zhu SG, Xiong W, Liu HY, Huang H, Zhou M, Li GY. BRD7, a novel bromodomain gene, inhibits G1-S progression by transcriptionally regulating some important molecules involved in ras/MEK/ERK and Rb/E2F pathways. J Cell Physiol. 2004;200:89–98.
Drost J, Mantovani F, Tocco F, Elkon R, Comel A, Holstege H, Kerkhoven R, Jonkers J, Voorhoeve PM, Agami R, Del Sal G. BRD7 is a candidate tumour suppressor gene required for p53 function. Nat Cell Biol. 2010;12:380–9.
Zhao R, Liu Y, Wang H, Yang J, Niu W, Fan S, Xiong W, Ma J, Li X, Phillips JB, Tan M, Qiu Y, Li G, Zhou M. BRD7 plays an anti-inflammatory role during early acute inflammation by inhibiting activation of the NF-kB signaling pathway. Cell Mol Immunol. 2017;14:830–41.
Wang H, Zhao R, Guo C, Jiang S, Yang J, Xu Y, Liu Y, Fan L, Xiong W, Ma J, Peng S, Zeng Z, Zhou Y, Li X, Li Z, Li X, Schmitt DC, Tan M, Li G, Zhou M. Knockout of BRD7 results in impaired spermatogenesis and male infertility. Sci Rep. 2016;6:21776.
Park SW, Herrema H, Salazar M, Cakir I, Cabi S, Basibuyuk Sahin F, Chiu YH, Cantley LC, Ozcan U. BRD7 regulates XBP1s' activity and glucose homeostasis through its interaction with the regulatory subunits of PI3K. Cell Metab. 2014;20:73–84.
Liu H, Li X, Niu Z, Zhang L, Zhou M, Huang H, He J, Zhang W, Xiao L, Tang Y, Wang L, Li G. Preparation of polyclonal antibody specific for BRD7 and detection of its expression pattern in the human fetus. J Histochem Cytochem. 2008;56:531–8. official journal of the Histochemistry Society
Harte MT, O'Brien GJ, Ryan NM, Gorski JJ, Savage KI, Crawford NT, Mullan PB, Harkin DP. BRD7, a subunit of SWI/SNF complexes, binds directly to BRCA1 and regulates BRCA1-dependent transcription. Cancer Res. 2010;70:2538–47.
Tae S, Karkhanis V, Velasco K, Yaneva M, Erdjument-Bromage H, Tempst P, Sif S. Bromodomain protein 7 interacts with PRMT5 and PRC2, and is involved in transcriptional repression of their target genes. Nucleic Acids Res. 2011;39:5424–38.
Liu T, Zhao M, Liu J, He Z, Zhang Y, You H, Huang J, Lin X, Feng XH. Tumor suppressor bromodomain-containing protein 7 cooperates with Smads to promote transforming growth factor-beta responses. Oncogene. 2017;36:362–72.
Liu Y, Zhao R, Wang H, Luo Y, Wang X, Niu W, Zhou Y, Wen Q, Fan S, Li X, Xiong W, Ma J, Li X, Tan M, Li G, Zhou M. miR-141 is involved in BRD7-mediated cell proliferation and tumor formation through suppression of the PTEN/AKT pathway in nasopharyngeal carcinoma. Cell Death Dis. 2016;7:e2156.
Liu H, Zhou M, Luo X, Zhang L, Niu Z, Peng C, Ma J, Peng S, Zhou H, Xiang B, Li X, Li S, He J, Li X, Li G. Transcriptional regulation of BRD7 expression by Sp1 and c-Myc. BMC Mol Biol. 2008;9:111.
Mathelier A, Fornes O, Arenillas DJ, Chen CY, Denay G, Lee J, Shi W, Shyr C, Tan G, Worsley-Hunt R, Zhang AW, Parcy F, Lenhard B, Sandelin A, Wasserman WW. JASPAR 2016: a major expansion and update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 2016;44:D110–5.
Adhikary S, Eilers M. Transcriptional regulation and transformation by Myc proteins. Nat Rev Mol Cell Biol. 2005;6:635–45.
Xu K, Xiong W, Zhou M, Wang H, Yang J, Li X, Chen P, Liao Q, Deng H, Li X, Li G, Zeng Z. Integrating ChIP-sequencing and digital gene expression profiling to identify BRD7 downstream genes and construct their regulating network. Mol Cell Biochem. 2016;411:57–71.
Liu C, Liu R, Zhang D, Deng Q, Liu B, Chao HP, Rycaj K, Takata Y, Lin K, Lu Y, Zhong Y, Krolewski J, Shen J, Tang DG. MicroRNA-141 suppresses prostate cancer stem cells and metastasis by targeting a cohort of pro-metastasis genes. Nat Commun. 2017;8:14270.
Hur K, Toiyama Y, Takahashi M, Balaguer F, Nagasaka T, Koike J, Hemmi H, Koi M, Boland CR, Goel A. MicroRNA-200c modulates epithelial-to-mesenchymal transition (EMT) in human colorectal cancer metastasis. Gut. 2013;62:1315–26.
de Rie D, Abugessaisa I, Alam T, Arner E, Arner P, Ashoor H, Astrom G, Babina M, Bertin N, Burroughs AM, Carlisle AJ, Daub CO, Detmar M, Deviatiiarov R, Fort A, Gebhard C, Goldowitz D, Guhl S, Ha TJ, Harshbarger J, Hasegawa A, Hashimoto K, Herlyn M, Heutink P, Hitchens KJ, Hon CC, Huang E, Ishizu Y, Kai C, Kasukawa T, Klinken P, Lassmann T, Lecellier CH, Lee W, Lizio M, Makeev V, Mathelier A, Medvedeva YA, Mejhert N, Mungall CJ, Noma S, Ohshima M, Okada-Hatakeyama M, Persson H, Rizzu P, Roudnicky F, Saetrom P, Sato H, Severin J, Shin JW, Swoboda RK, Tarui H, Toyoda H, Vitting-Seerup K, Winteringham L, Yamaguchi Y, Yasuzawa K, Yoneda M, Yumoto N, Zabierowski S, Zhang PG, Wells CA, Summers KM, Kawaji H, Sandelin A, Rehli M, Consortium F, Hayashizaki Y, Carninci P, Forrest ARR, de Hoon MJL. An integrated expression atlas of miRNAs and their promoters in human and mouse. Nat Biotechnol. 2017;35:872–8.
Fabbri M, Bottoni A, Shimizu M, Spizzo R, Nicoloso MS, Rossi S, Barbarotto E, Cimmino A, Adair B, Wojcik SE, Valeri N, Calore F, Sampath D, Fanini F, Vannini I, Musuraca G, Dell'Aquila M, Alder H, Davuluri RV, Rassenti LZ, Negrini M, Nakamura T, Amadori D, Kay NE, Rai KR, Keating MJ, Kipps TJ, Calin GA, Croce CM. Association of a microRNA/TP53 feedback circuitry with pathogenesis and outcome of B-cell chronic lymphocytic leukemia. JAMA. 2011;305:59–67.
Psathas JN, Thomas-Tikhonenko A. MYC and the art of microRNA maintenance. Cold Spring Harbor perspectives in medicine. 2014;4
Kulshreshtha R, Davuluri RV, Calin GA, Ivan M. A microRNA component of the hypoxic response. Cell Death Differ. 2008;15:667–71.
Blackwood EM, Eisenman RN. Max: a helix-loop-helix zipper protein that forms a sequence-specific DNA-binding complex with Myc. Science. 1991;251:1211–7.
Blackwood EM, Luscher B, Eisenman RN. Myc and max associate in vivo. Genes Dev. 1992;6:71–80.
Qi Y, Gregory MA, Li Z, Brousal JP, West K, Hann SR. p19ARF directly and differentially controls the functions of c-Myc independently of p53. Nature. 2004;431:712–7.
Datta A, Nag A, Pan W, Hay N, Gartel AL, Colamonici O, Mori Y, Raychaudhuri P. Myc-ARF (alternate reading frame) interaction inhibits the functions of Myc. J Biol Chem. 2004;279:36698–707.
Burrows AE, Smogorzewska A, Elledge SJ. Polybromo-associated BRG1-associated factor components BRD7 and BAF180 are critical regulators of p53 required for induction of replicative senescence. Proc Natl Acad Sci U S A. 2010;107:14280–5.
Dang CV. MYC on the path to cancer. Cell. 2012;149:22–35.
Kirschner AN, Wang J, van der Meer R, Anderson PD, Franco-Coronel OE, Kushner MH, Everett JH, Hameed O, Keeton EK, Ahdesmaki M, Grosskurth SE, Huszar D, Abdulkadir SA. PIM kinase inhibitor AZD1208 for treatment of MYC-driven prostate cancer. J Natl Cancer Inst. 2015;107
Song MS, Salmena L, Pandolfi PP. The functions and regulation of the PTEN tumour suppressor. Nat Rev Mol Cell Biol. 2012;13:283–96.
Amendola D, De Salvo M, Marchese R, Verga Falzacappa C, Stigliano A, Carico E, Brunetti E, Moscarini M, Bucci B. Myc down-regulation affects cyclin D1/cdk4 activity and induces apoptosis via Smac/Diablo pathway in an astrocytoma cell line. Cell Prolif. 2009;42:94–109.
Vlach J, Hennecke S, Alevizopoulos K, Conti D, Amati B. Growth arrest by the cyclin-dependent kinase inhibitor p27Kip1 is abrogated by c-Myc. EMBO J. 1996;15:6595–604.
Alimonti A, Carracedo A, Clohessy JG, Trotman LC, Nardella C, Egia A, Salmena L, Sampieri K, Haveman WJ, Brogi E, Richardson AL, Zhang J, Pandolfi PP. Subtle variations in Pten dose determine cancer susceptibility. Nat Genet. 2010;42:454–8.
This work was supported by the National Natural Science Foundation of China (81772990, 81572748); the Natural Science Foundation of Hunan Province (2015JJ2148); the National Key R&D Program of China (2016YFC1201803); the “111” project (111–2-12); the Fundamental Research Funds for the Central Universities of Central South University (2016zzts114) and the Open-End Fund for the Valuable and Precision Instruments of Central South University (CSUZC201743).
Availability of data and materials
In our study, ethical approval was obtained from the Ethics Review Committees/Institutional Review Boards of Central South University. NPC samples and non-cancerous nasopharyngeal tissues from healthy donors were collected at the Second Xiangya Hospital of Central South University (Changsha, China). Written informed consent was obtained from all patients and donors. All animal procedures were conducted in accordance with protocols approved by the Institutional Animal Care and Use Committee (IACUC) of Central South University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The primers for confirmation of the c-Myc-binding sites in the miR-141 and BRD7 promoter. (DOC 34 kb)
Figure S1. c-Myc did not influence the processing of pri-miR-141 and pre-miR-141 in miR-141 biosynthesis. Figure S2. BRD7 was inactivated by c-Myc and did not affect the stability of c-Myc protein. Figure S3. The levels of miR-141 restoration in c-Myc knockdown NPC cells. (DOC 2382 kb)
Table S2. Ten putative c-Myc-binding sites were predicted in the miR-141 promoter region from the JASPAR Database. (DOC 40 kb)
Table S3. Protein-Protein Interaction predictions for BRD7. (XLS 49 kb)