Loss of MYC and E-box3 binding contributes to defective MYC-mediated transcriptional suppression of human MC-let-7a-1~let-7d in glioblastoma
Metrics: PDF 765 views | HTML 940 views | ?
Zifeng Wang1, Sheng Lin4, Ji Zhang1, Zhenhua Xu5, Yu Xiang1, Hong Yao6, Lei Ge8, Dan Xie1, Hsiang-fu Kung7, Gang Lu3, Wai Sang Poon3, Quentin Liu1, Marie Chia-mi Lin2
1Sun Yat-Sen University Cancer Center, State Key Laboratory of Oncology in South China, Collaborative Innovation Center of Cancer Medicine, Guangzhou, China; Institute of Cancer Stem Cell, Dalian Medical University, Dalian, China
2Shenzhen Key Laboratory of Translational Medicine of Tumor, School of Medicine, Shenzhen University, Shenzhen, China
3Brain Tumor Centre and Division of Neurosurgery, Department of Surgery, Faculty of Medicine, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, Hong Kong, China
4Laboratory of Medical Genetics, Shenzhen Research Institute of Population and Family Planning, Shenzhen, China
5Biochemistry and Molecular Biology, Medical University of South Carolina, Charleston, SC, USA
6Jiangsu Eng. Laboratory of Cancer Biotherapy, Xuzhou Medical College, Xuzhou, China
7School of Biomedical Science, and State Key Laboratory in Oncology in South China, Chinese University of Hong Kong, Shatin, Hong Kong
8Department of Gastrointestinal Surgery, Tumor Hospital, Xinjiang Medical University, Urumqi, Xinjiang Uyghur Autonomous Region, China
Marie Chia-mi Lin, email: email@example.com
Quentin Liu, email: firstname.lastname@example.org
Keywords: MC-let-7a-1~let-7d, MYC, glioblastoma, transcriptional regulation
Received: October 16, 2015 Accepted: May 20, 2016 Published: July 09, 2016
Previously, we reported that MYC oncoprotein down-regulates the transcription of human MC-let-7a-1~let-7d microRNA cluster in hepatocarcinoma (HCC). Surprisingly, in silico analysis indicated that let-7 miRNA expression levels are not reduced in glioblastoma (GBM). Here we investigated the molecular basis of this differential expression. Using human GBM U87 and U251 cells, we first demonstrated that forced over-expression of MYC indeed could not down-regulate the expression of human MC-let-7a-1~let-7d microRNA cluster in GBM. Furthermore, analysis of MC-let-7a-1~let-7d promoter in GBM indicated that MYC failed to inhibit the promoter activity. Pearson's correlation and Linear Regression analysis using the expression data from GSE55092 (HCC) and GSE4290 (GBM) demonstrated a converse relationship of MC-let-7a-1~let-7d and MYC only in HCC but not in GBM. To understand the underlying mechanisms, we examined whether MYC could bind to the non-canonical E-box 3 located in the promoter of MC-let-7a-1~let-7d. Results from both chromatin immune-precipitation (ChIP) and super-shift assays clearly demonstrated the loss of MYC and E-box 3 binding in GBM, suggesting for the first time that a defective MYC and E-box3 binding in GBM is responsible for the differential MYC mediated transcriptional inhibition of MC-let-7a-1~let-7d and potentially other tumor suppressors. MYC and let-7 are key oncoprotein and tumor suppressor, respectively. Understanding the molecular mechanisms of their regulations will provide new insight and have important implications in the therapeutics of GBM as well as other cancers.
MicroRNAs (miRNAs) are short non-coding RNAs (~22 nt) that regulate protein expression and control diverse aspects of biology, including carcinogenesis . The let-7 family miRNAs are down-regulated in various types of cancer including lung cancer, liver cancer, colon cancer, ovarian cancer, stomach cancer, leiomyoma and melanoma [2–15]. They have become a prototype miRNAs that function as tumor suppressors, as they inhibit the expressions of multiple oncogenes including Ras and MYC [16–19]. The human microRNA cluster MC-let-7a-1~let-7d consists of three members let-7a-1, let-7f-1 and let-7d. These three members together accounted for approximately 24% of the total let-7 precursors . Hence, this cluster is a particularly important component of the let-7 family.
The MYC oncoprotein is a transcription factor which plays a pivotal role in carcinogenesis by transcriptional down-regulation of multiple tumor suppressors. MYC expression is significantly elevated in almost all cancers including HCC and GBM. We have previously demonstrated that MYC inhibits MC-let-7a-1~let-7d expression at transcriptional level through a non-canonical E-box 3 in HCC .
Both HCC and Glioblastoma (GBM, WHO grade IV, the most severe form of glioma) are highly invasive and destructive tumors with no effective treatment [20, 21]. Accumulated evidence has clearly shown the miRNAs play important roles in the carcinogenesis of HCC and GBM [15, 18, 22–26]. Using bioinformatics analysis of human miRNA expression data sets, we showed that the tumor suppressor MC-let-7a-1~let-7d is down-regulated in most types of cancer include HCC. Surprisingly, in silico analysis indicated that let-7 miRNA expression levels are not reduced in GBM [5, 27–35]. As MYC oncoprotein expression is elevated in both HCC and GBM, we investigated why MC-let-7a-1~let-7d expression is down-regulated in HCC but not in GBM.
Contrary to HCC HepG2 and L02 cells, we demonstrated that in human GBM U87 and U251 cells, forced over-expression of MYC could not down-regulate the expression of human MC-let-7a-1~let-7d microRNA cluster in GBM. In addition, we characterized MC-let-7a-1~let-7d promoter in GBM, and showed that MYC failed to inhibit the promoter activity of MC-let-7a-1~let-7d. As MYC down-regulates transcription mainly through the binding with a non-canonical E-box 3, we investigate the binding of MYC and the E-box 3 located in the promoter of MC-let-7a-1~let-7d in GBM. By both chromatin immune-precipitation (ChIP) and super-shift assays, we demonstrated the loss of MYC and E-box 3 binding. Taken together, these results revealed for the first time that GBM exhibited differential MYC mediated transcriptional inhibition on MC-let-7a-1~let-7d due to the defective MYC/E-box3 binding. MYC oncoprotein and MC-let-7a-1~let- 7d both play key roles in carcinogenesis. Understanding the molecular basis of their functions and regulations will provide new insight and important implications in the development of new therapeutics of GBM and other cancers.
The expression levels of let-7 miRNAs are not reduced in GBM
We first analyzed the expression level of 21 mature miRNA of let-7 family members in GBM as compared to normal brain tissues using the microarray data from 11 independent groups [27–35]. Results from ten groups indicated that let-7 miRNA levels are not significantly changed in GBM, and one group showed a modest down-regulation of only 3 miRNAs, let-7d, let-7f and let-7g (Supplementary Table S1).
Next, we compared the expression level of precursor miRNA of let-7 family members in HCC/normal liver tissues and GBM/normal brain tissues. As shown in Table 1 (data from Pablo Landgraf. et al ), in 147 human cells and tissues, both the MC-let-7a-1~let-7d and the let-7 family precursors are significantly down-regulated (more than 10- fold reduction) in HCC. Surprisingly, these precursors are slightly up-regulated (less than 2-fold elevation) in GBM. The miR-21, miR-23a, miR-210 and miR-139, miR- 128 and miR-181 are listed as controls to shown that the deep sequence data (NGS) are consistent with the published results [27, 35–39]. The unbiased data-mining results suggest that MC-let-7a-1~let-7d is not down-regulated in GBM, in sharp contrast to its expression pattern in HCC. Consistent with the published reports, high level of MYC onco-protein expression in GBM was observed by immune-histochemical staining assays using The Human Protein Atlas database  (Figure 1A, left panel). Moreover, we also analyzed the expression of MYC in a large number of clinical glioblastoma specimens (n = 495) in The Cancer Genome Atlas (TCGA) glioblastoma database. Results indicated that MYC expression is significantly elevated in all glioblastoma subtypes relative to normal human cerebrum (Figure 1A, right panel).
Table 1: The ratio of let-7 in total miRNAs in GBM vs normal brain tissues (NBT) and HCC vs normal liver tissues (NLT)
miRNA/ intracellular total miRNAs
Annotation: The ratios are based on 147 human cells and tissues detected by Pablo Landgraf. et al. 2007. miR-21, miR-23a, miR-210 and miR-139, miR-128, miR-181 are listed to show that the deep sequence data (NGS) are consistent with the pubilished results. Log FC, log10-transformed fold change (tumor/normal).
Figure 1: MC-let-7a-1~let-7d is not down-regulated, whereas Myc is up-regulated in GBM. (A) left panel, the expression of c-Myc protein in normal brain tissues and GBM as determined by immunohistochemical staining obtained from Human Protein Atlas database. Right panel, the expression of c-Myc mRNA in normal brain tissues and GBM as determined by high throughput sequencing method from The Cancer Genome Atlas (TCGA) glioblastoma database. Patient specimens were categorized according to TCGA subtypes: classical, mesenchymal, proneural, and neural. (B) the expression levels of MC-let-7a-1~let-7d and MYC in normal brain tissues and GBM according to GEO Profile (GSE4290). FC, fold change.(C) the relative expression levels of primary MC-let-7a-1~let-7d members in normal human brain tissue(NBT), in human glioma cell lines (HG = highgrade; LG = low grade), and in a primitiveneuroectodermal tumor (PNET) cell line. (D) the relative expression levels of c-MycinNBT, in human glioma cell lines, and in PNETcell line. *p < 0.05; **p < 0.01; ***p < 0.001.
To further confirm these findings, we carried out gene-expression analyses of MC-let-7a-1~let-7d and MYC in an independent gene-expression cohort including 81 GBM tissues and 23 normal brain tissues according to GEO Profile (GSE4290). The GPL570 microarray platform provides probes to detect the expression of let-7d host gene (MC-let-7a-1~let-7d). As shown in Figure 1B, while MYC level was significantly up-regulated 2.5-fold (p < 0.0001), the expression of MC-let-7a-1~let-7d was not significantly changed (1.13 fold, p = 0.13).
We further examined the expression of pri-let-7a-1, pri-let-7f-1, pri-let-7d and MYC in three human GBM cell lines from high-grade tumors (U-373 MG, U-87 MG, and U-251 MG), one low-grade glioma (SW-1088), and a malignant neuroectodermal tumor (PFSK-1) by Real-Time PCR assays. Consistently, the pri-let-7a-1, pri-let-7f-1, and pri-let-7d levels were not significantly changed as compared to normal brain tissue (NBT, Figure 1C), while MYC level was significantly elevated in all cell lines (Figure 1D).
Characterization of the promoter of MC-let-7a-1~let-7d cluster in GBM U87 and U251 cells
Previously, we have characterized the MC-let-7a-1~let-7d promoter in HCC . Here, we characterize the promoter and transcription start site(s) of MC-let-7a-1~let-7 in GBM U87 and U251 cells by 5’ Rapid Amplification of cDNA Ends (5’RACE) assays using total RNA as templates, gene-specific primers GSP1-1 for the synthesis of first strand cDNA, and GSP2-1 and GSP2-2 for nested PCR (Figure 2A). A ~700 bp product was amplified (Figure 2B) and cloned. Sequence analysis of twenty 5’RACE products (10 products for each cells) revealed that the TSSs of MC-let-7a-1~let-7d in U87 and U251 are located about 28 bp and 40 bp, respectively, downstream of the TSSs found in HepG2 cells (Figure 2C). These results suggested that MC-let-7a-1~let-7d uses similar promoter regions in HCC and GBM with different TSSs.
Figure 2: MC-let-7a-1~let-7d promoter characterization in GBM as compared to that in HCC. (A) primers used for 5′RACE. Gene-specific primers (GSPs) GSP1-1 was used to synthesize the first strand cDNA. Primary amplification was carried out with abridged anchor primer (AAP) and GSP2-1 primer. Nested PCR was performed with abridged universal amplification primer (AUAP) and GSP2-2.(B) 5′RACE results. The ~700-bp products were amplified in U87 cells and U251 cells by using GSP2-2 and AUAP. (C) sequence analysis of 20 products from each cell lines revealed that variable TSSs were used. The major TSSs of MC-let-7a-1~let-7d in U87 is located about 28 bp, while in U251 is located about 40 bp downstream of the TSSs used in HepG2 cells. (D) schematic diagram of the reporter construct (left panel). Truncation of the downstream region to +1 dramatically improved the promoter activity, whereas truncation of the upstream region to −227 resulted in progressive loss of activity in both U87 cells and L02 cells.
In addition, we compared the promoter activity of the highly conserved region around the TSSs. Our results showed that the promoter regions of MC-let-7a-1~let-7d exhibited similar property in GBM (U87 cells) as compared to that of the liver L02 cells. As shown in Figure 2D, truncation analysis revealed that promoter (−540bp to +1bp) exhibited maximal activity. In addition, the presence of a suppression element from +1317 to +1 and an activator from −540 to −227 were detected in both U87 cells and L02 cells.
MYC failed to inhibit the expression of human MC-let-7a-1~let-7d microRNA cluster in GBM cells.
To investigate whether MYC differentially regulates the transcription of MC-let-7a-1~let-7d in HCC and GBM, we over-expressed MYC in HepG2, L02, U87, and U251 cells. As expected, forced expression of MYC significantly down-regulated MC-let-7a-1~let-7d expression only in L02 and HepG2 cells (Figure3A), but not in U87 and U251 cells (Figure 3B).
Figure 3: MYC oncoprotein down-regulates MC-let-7a-1~let-7d in HCC, but not in GBM. (A) forced expression of MYC significantly down-regulated the MC-let-7a-1~let-7d expression in L02 and HepG2 cells. (B) forced expression of c-Myc failed to down-regulate the MC-let-7a-1~let-7d expression in U87 and U251 cells. (C) the expression correlativity of MYC and MC-let-7a-1~let-7d in HCC patient samples (right panel) and non-tumor control (left panel) according to GSE55092 data set. (D) the expression correlativity of MYC and MC-let-7a-1~let-7d in GBM patient samples (right panel) and non-tumor control (left panel) according to GSE4290 data set. Pearson correlation is calculated to measure the strength of a linear association between the expressions of two genes. R-square is the coefficient of linear regression which measures the linearity between the expressions of two genes. *p < 0.05; **p < 0.01.
We further validated these findings in HCC- and GBM-patient samples. Pearson’s correlation and Linear Regression analysis were employed using data from GSE55092 (HCC, tumor = 49, non-tumor = 91) and GSE4290 (GBM, tumor = 81, non-tumor = 23). The converse relationship of MC-let-7a-1~let-7d and MYC was only observed in HCC samples (Figure 3C, right panel, Pearson = 0.408, p = 0.004), but not in normal liver samples (Figure 3C, left panel, Pearson = −0.066, p = 0.532), normal brain samples (Figure 3D, left panel, Pearson = −060,p = 0.785) or glioblastoma samples (Figure 3D, right panel, Pearson = 0.181, p = 0.106).
MYC failed to down-regulate the promoter activity of MC-let-7a-1~let-7d cluster in GBM
Our previous study showed that MYC down-regulates the transcription of MC-let-7a-1~let-7d in HCC, and identified two functional MYC binding E-box elements (E-box 2, 3) . We showed that over-expression of MYC significantly reduced the expression of MC-let-7a-1~let- 7d primary transcript through the non-canonical E-box 3, in L02 and HepG2 cells. To examine the transcriptional regulation in GBM, we first determined the effect of MYC in MC-let-7a-1~let-7d promoter PPR-10 which only contains canonical E-box 2 (Figure 4A). Our results indicated that over-expression of MYC activated luciferase activity of PPR-10 in all the four cells lines HepG2, L02, U87, and U251, suggesting the similar function of basal activity of canonical E-box 2 (Figure 4B). Furthermore, over-expression of MYC significantly repressed the activity of PPR-3 which contains both canonical E-box 2 and non-canonical E-box 3 in HepG2 and L02 cells, however, it could not change the activity of PPR-3 in GBM U87 and U251 cells (Figure 4C). These results suggested a defective function of the non-canonical E-box 3 in GBM.
Figure 4: MYC failed to bind to the repressive non-canonical E-box 3 in GBM cells. (A) schematic representation of the location of MC-let-7a-1~let-7d promoter PPR-3 and PPR-10 and the E-boxes. (B) overexpression of MYC increased PPR-10 activity in both HCC cells HepG2 and L02 and GBM cells U251 and U87. (C) overexpression of MYC decreased PPR-3 activity in HCC cells HepG2 and L02, whereas failed to decrease its activity in GBM cells U251 and U87. (D) ChIP assays were performed to quantify the MYC binding quantum on E-box 2 and E-box 3. Schematic representation of the location of Ctrl, E-box 2, and E-box 3 (upper panel) is shown. The Real-time PCR results were summarized as relative binding (MYC ChIP relative to input ChIP). The E-box 2 region was significantly enriched, whereas the E-box 3 region was not enriched both in U87 and U251. (E) EMSA and supershift assay analysis of MYC binding to E-box 3 in U87 and U251. Competition assays were performed adding to 100-fold molar unlabeled E-boxes. Shift bands were observed by adding nuclear extract form both U87 and U251, whereas supershift bands was not observed by adding nuclear extract form U87 and U251 with MYC antibody.
Defective MYC and E-box 3 binding in GBM
We determined the binding ability of MYC to E-box 3 in GBM U87 and U251 cells by ChIP assays to quantify the MYC binding quantum on E-box 2 and E-box 3. The results from Real-time PCR experiments were summarized and shown as relative binding (MYC ChIP relative to input ChIP). As shown in Figure 4D, strong MYC binding was observed in E-box 2 but not in E-box 3, suggesting the absence of MYC and E-box 3 binding in GBM.
Furthermore, EMSA assays and super-shift analysis were performed to confirm the loss of MYC binding to E-box 3 in U87 and U251 cells. Competition assays were conducted by adding 100-fold molar unlabeled E-boxes. As shown in Figure 4E, E-box 3 showed shift bands by adding nuclear extract from U87 and U251 cells. However, supershift bands were not observed by adding nuclear extract form U87 and U251 with MYC antibody. These results confirmed that MYC onco-protein could not bind to E-box 3 in GBM.
MYC is an evolutionarily conserved E-box binding transcription factor. It is pathologically activated in many human malignancies . The expression of MYC is important in cell signaling pathways such as cell growth/ proliferation, metabolism and apoptosis . It has been reported that MYC expression is elevated in almost all cancers including both HCC and GBM [18, 41]. As compared to a nearly undetectable level in normal brain tissues, MYC is expressed in 80.6% of primary glioblastoma samples [43, 44], and in 73% of astrocytomas . The wide spread repression of miRNAs including let-7 by MYC and its contribution to tumorigenesis has also been clearly demonstrated .
Previously, we reported that MYC inhibited the transcription of MC-let-7a-1~let-7d in HCC . We proposed that MYC exhibited dual functions in the transcriptional regulation of MC-let-7a-1~let-7d. As shown in Figure 5, in non-cancerous L02 cells with a low level of MYC, MYC dimerizes with Max and bind to the canonical E-box 2 to activate and maintain the basal expression level of let-7. However, in cancerous HepG2 cells with elevated MYC, the extra MYC binds to the non-canonical E-box 3 to suppress the transcription of MC-let-7a-1~let-7d. These mechanisms enable cancer cells to maintain a high level of MYC and a low level of let-7 miRNAs simultaneously in HCC.
Figure 5: Model depicts the mechanism of the differential transcriptional regulation of MC-let-7a-1~let-7d in GBM and HCC. In HCC, a small amount of c-Myc preferentially dimerizes with Max and binds E-box 2 to maintain the basal activity of MC-let-7a-1~let-7d promoter. The extra c-Myc will bind to E-box 3 and repress the promoter activity of MC-let-7a-1~let-7d. In GBM, likewise, a small amount of c-Myc preferentially dimerizes with Max and binds E-box 2 to enhance the transcription of MC-let-7a-1~let-7d. Differently, the extra c-Myc failed to bind to non-canonical E-box 3 to inhibit the expression of MC-let-7a-1~let-7d in GBM cells.
In present study, we reported that MC-let-7a-1~let-7d is not down-regulated in GBM due to the loss of MYC and E-box 3 binding. The detailed molecular mechanism for this defect is not known at present. One explanation is that the binding of MYC to E-box 3 requires additional unknown intermediary trans-acting factors that are defective or absent in GBM.
In the promoter truncation assay, we found that the −540 bp to +1bp region of MC-let-7a-1~let-7d promoter exhibited maximal promoter activity, while longer promoters PPR-1~4 showed decreased promoter activity. These results indicated the presence of suppression elements downstream of the transcription initiation site. Indeed, our study showed that MYC also acts as a transcription suppressor by binding to the non-canonical binding site E-box 3 located downstream of TSSs. However, it is not clear whether and how MYC bind to the non-canonical DNA sequence to suppress the promoter activity. Results from our study inferred that MYC might not bind to E-box3 directly, but is recruited to the promoter by yet to be identified factor, which is present in HCC, but either defective or absent in GBM. Identification and characterization of this MYC-E- box3 binding partner will have important implications as therapeutic target and biomarkers for GBM and cancers.
It has been reported that the expression of let-7 family members are increased during neural differentiation and let-7 members play important roles in neural stem cell proliferation and differentiation [47–50]. Lee et al showed that over-expression of let-7 miRNAs could reduce the proliferation and migration of GBM cells, as well as the tumor size in nude mice xenograft transplanted GBM . Forced expression of let-7 miRNA reduced expression of pan-RAS, N-RAS, and K-RAS , suggesting that let-7 played important tumor suppressor roles in GBM. Our in silico and wet-lab findings suggest that let-7 escapes MYC inhibition in GBM. It is not known whether the co-existence of high levels of onco-MYC and tumor suppressor-let-7 had any impact in the carcinogenesis of brain tumors. As let-7 miRNAs could inhibit the expression of MYC, and how MYC escapes inhibition by let-7 in GBM is also not understood. The existence of recently discovered competing endogenous RNA (ceRNA), and whether they play a role requires further investigations [52, 53].
The present study provided a model to explain the differential transcriptional regulation of tumor suppressor let-7 by oncogene MYC in glioblastoma and hepatocellular carcinoma. These results suggest a potential molecular mechanism for understanding the complexity of molecular regulatory network. The potential implications in cancer treatment warrant further investigations.
MATERIALS AND METHODS
Cell lines and cell culture
Human HCC cell line HepG2 and immortalized human liver cell line L02 were obtained from American Type Culture Collection (ATCC, Manassas, VA). Human high-grade glioma cell lines (U-87 MG, U-251 MG, and U-373 MG), human primitive neuroectodermal PFSK-1 cells, and human low-grade glioma SW-1088 cells were also obtained from ATCC. All these cell lines were cultured in Dulbecco’s modified eagles medium (DMEM, Invitrogen) supplemented with 10% fetal bovine serum (FBS, Gibco) at 37°C in a humidified atmosphere with 5% CO2 as our previous report [15, 54].
Three adjacent normal human brain tissue (NBT) specimens were collected from 3 male GBM patients with 23-, 34- and 35-year-old respectively at Sun Yat-sen University Cancer Center. Patients’ consent and approval from Sun Yatsen University Cancer Center Institute Research Ethics Committee were obtained for the use of these clinical materials. For each sample, total RNA were purified by using AllPrep Mini Kit (Qiagen, Cat number 80004) according to the manufactures’ instructions.
Genomic DNA extraction, RNA extraction, cDNA synthesis and Real-Time PCR
For genomic DNA extraction, total cell genomic DNA was extracted using the DNeasy Blood & Tissue (Qiagen), following the manufacturer’s instructions. For RNA extraction, total RNA was extracted using Trizol (Invitrogen), according to the manufacturer’s instructions. For cDNA synthesis, 2 μg total RNA was reverse transcribed by the SuperScript™ III First-Strand Synthesis System (Invitrogen) after DNase I treatment (Invitrogen). Quantitive Real-Time polymerase chain reaction (qRT-PCR) was subsequently performed in triplicate with a 1:10 dilution of the resultant cDNA using the Brilliant Brilliant® II SYBR® Green QPCR Master Mix (Stratagene) on the Applied Biosystems 7500 Real-Time PCR System Instrument (Applied Biosystems). All mRNA quantification data were normalized to GAPDH. Ct-value for each sample was calculated with the ΔΔCt-method and results were expressed as 2-ΔΔCT .
Protein Extraction and Western blot analysis
Cells were seeded onto a 6-well plate at a density of 106 cells/well. Cells were lysed for 30 min on ice in RIPA lysis buffer (Bioteke) with protease inhibitor cocktail. Protein concentration was determined by Bio-Rad protein assay kit (Bio-Rad). Fifty ug proteins were separated by 12% SDS polyacrylamide gels (SDS-PAGE) and transferred onto PVDF membranes. Membranes were probed with specific antibodies includingMYC (1:2000; Abcam) or β-actin (1:3000; Abcam) respectively. After washing, membranes were incubated with HRP-conjugated secondary antibody (1:5000; Abcam). Membranes were exposed using the enhanced chemiluminescence (ECL) substrate kit (Thermo) and imaged using the ImageQuant™ RT ECL™ system (GE Healthcare). The expression level of β-actin was also determined and served as an internal control.
5′ Rapid Amplification of cDNA Ends (5′ RACE)
The 5′ RACE System for Rapid Amplification of cDNA Ends, Version 2.0 (Invitrogen) was used according to the manufacturer’s instructions. Briefly, cDNA was synthesized using the MC-let-7a-1~let-7d Specific Primers GSP1-1 (Supplementary Table S2). Primary amplification was carried out with abridged anchor primer (AAP, provided in the kit) and GSP2-1 primer (Supplementary Table S2). Nested PCR was performed with abridged universal amplification primer (AUAP, provided in the kit) and GSP2-2 (Supplementary Table S2) using TaKaRa LA Taq® with GC Buffer I (TaKaRa). The 5′RACE PCR products were resolved on 1% agarose gel and stained by EtBr. The bands were excised, cloned into pGL3 Basic Vector (Promega) via restriction endonuclease Mlu I sites and sequenced.
Constructions of promoter reporter, and MYC expression plasmids
A 1.9 kb MC-let-7a-1~let-7d promoter was amplified by PCR using PrimeSTAR™ HS DNA Polymerase (TaKaRa) and cloned into pGL3 Basic Vector (Promega) between Xho I and Hind III sites using primers PPR-1 as listed in Supplementary Table S2. This plasmid was named as PPR-1. Subsequent promoter truncations were generated from this template and cloned in a similar manner using the primers in Supplementary Table S2. To construct MYC expression vector, human MYC open reading frame (ORF) was amplified using MYC primers as indicated in Supplementary Table S2 and inserted into the pcDNA3.1(+) vector (Invitrogen) between Hind III and BamH I sites. The sequence and orientation of all of the inserts were confirmed by sequencing.
Luciferase reporter assays
HepG2 cells, L02 cells, U87 cells or U251 cells cells were plated at a density of 50,000 cells per well in a 24-well plate at 24 hrs before transfection. Cells were co-transfected with 250 ng MC-let-7a-1~let-7d promoter/firefly luciferase reporter plasmids and 5 ng pRL-TK Renilla plasmids (Promega) using FuGENE® HD Transfection Reagent (Promega). For gain and loss-of-function experiments, 750 ng MYC expression plasmids (c-Myc) or its control vectors (pcDNA3.1) were transfected. After 48 h post-transfection, cells were either lysed in Trizol (Invitrogen) for RNA extraction or in passive lysis buffer (Promega) for luciferase assay measured with the Dual-Luciferase Reporter Assay System (Promega) using the TD-20/20 Luminometer (Turner Designs). The relative luciferase activities were determined by calculating the ratio of firefly luciferase activities over Renilla luciferase activities.
Chromatin immunoprecipitation (ChIP)
ChIP assays were carried out by using an EZ-ChIP assay kit (Millipore), in accordance with the manufacturer’s instructions. In brief, U87 and U251 cells were grown to 90% confluence and added 1% formaldehyde at room temperature for 10 min. The cross-link reaction was quenched with 0.125 M glycine for 5 min at room temperature. Cells were then washed, scraped, and resuspended in 1 ml of lysis buffer. DNA was sonicated into around 400-bp pieces at 4°C by using Sonics UibraCellsTM (XINCHEN). Supernatants were recovered by centrifugation and precleared for 1 h at 4°C with 60 μl of protein G-agarose. Then 10 μl (1%) of supernatant was removed as input. Immunoprecipitations were performed overnight with MYC antibody (Abcam, 9E11, ChIP Grade, 6 ug) or IgG antibody (provided in kit, 1 ug). The immune complexes were captured by incubation with 60 ulof protein G-agarose for 2 h at 4°C. The immunoprecipitates were washed sequentially with wash buffers. After that, the immunoprecipitates were eluted from the protein G-agarose by incubating with elution buffer (1% SDS, 100 mMNaHCO3). DNA-protein complexes were reversely cross-linked by a high salt solution at 65°C for 5 h. RNA and protein were eliminated by treating with 10 ug of RNase A at 37°C for 30 min and then with protease K for 2 h at 45°C. Finally, DNA was purified by using the spin column provided in the ChIP kit and eluted with 50 ul of elution buffer. qRT-PCR was performed by using Brilliant Brilliant® II SYBR® Green QPCR Master Mix (Stratagene) on the Applied Biosystems 7500 Real-Time PCR System Instrument (Applied Biosystems). The primer pairs used for PCR analysis are shown in Supplementary Table S2. All data were normalized to input as the previous report .
Electrophoretic mobility shift assay (EMSA)
Nuclear extracts from U87 and U251 cells were prepared in accordance with the manufacturer’s protocol (NE-PER® Nuclear and Cytoplasmic Extraction Reagents, Pierce). Complementary oligonucleotide pairs corresponding to the E-box 2 and E-box 3 embedded in the promoter region of MC-let-7a-1~let-7d and 5′ end-labeled with biotin by Invitrogen. The oligonucleotide sequences were shown in Supplementary Table S2. EMSA assay was performed by using a LightShift® Chemiluminescent EMSA Kit (Pierce). Binding reaction with 10 μg of HepG2 nuclear extracts and 100 fmol 5′ Biotin-labeled oligonucleotide was carried out in accordance with the manufacturer’s instructions. For the competition assay, a 100-fold molar excess of unlabeled oligonucleotide was added to the binding reaction mixture as a specific competitor. For antibody-supershift assay, a nuclear extract was pre-incubated with a 3 μl MYC antibody (Abcam) before adding it to the binding reaction. DNA/protein complexes were separated on a pre-electrophoresed 6% polyacrylamide gel in 0.5 × TBE, transferred to a nylon membrane, and cross-linked at 120 mJ/cm2 for 1 minute and detected by chemiluminescence, in accordance with the manufacturer’s directions. Membranes were exposed by using the enhanced chemiluminescence (ECL) substrate kit (Thermo) and then were photo-documented.
Data were shown as means ± SEM. Statistical analyses for detection of significant differences between the control and experimental groups were carried out using unpaired two-tailed Student’s t test (*p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001).
Hepatocarcinoma, HCC; glioblastoma, GBM; chromatin immunoprecipitation, ChIP; microRNA, miRNAs; normal brain tissues, NBT; normal liver tissues, NLT; 5′ Rapid Amplification of cDNA Ends, 5′RACE; gene-specific primers, GSP; American Type Culture Collection, ATCC; Dulbecco’s modified eagles medium, DMEM; fetal bovine serum, FBS; Electrophoretic Mobility Shift Assay, EMSA; chemiluminescence, ECL; sodium dodecyl sulfate-polyacrylamide gel electrophoresis, SDS-PAGE; polyvinylidene difluoride, PVDF.
We thank other members of Marie Chia-mi Lin lab and Quentin Liu lab for their critical comments and technical support.
CONFLICTS OF INTEREST
All authors claimed no conflicts of interest.
This research work was supported by the National Natural Science Foundation of China (81301734 to Z.W., 81272549 to M. L., 81130040 and 81573025 to Q.L.), Basic Scientific Research Fund of Sun Yat-sen University (14ykpy41 to Z.W.), Key Lab Project of Shenzhen (ZDSY20130329101130496 to M.L.), and RGC (CUHK470911 to M. L.), National Basic Research Program of China (973 Program, No. 2012CB967000 to Q.L.), and Innovative Research Team in University of Ministry of Education of China (No.IRT13049 to Q.L.).
1. Lee YS, Dutta A. MicroRNAs in cancer. Annu Rev Pathol. 2009; 4:199–227.
2. Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, Wilson M, Wang X, Shelton J, Shingara J, Chin L, Brown D, Slack FJ. The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Res. 2007; 67:7713–7722.
3. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, Calin GA, Liu CG, Croce CM, et al. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006; 9:189–198.
4. Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, Mitsudomi T, Takahashi T. Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004; 64:3753–3756.
5. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, Lin C, Socci ND, Hermida L, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007; 129:1401–1414.
6. Akao Y, Nakagawa Y, Naoe T. let-7 microRNA functions as a potential growth suppressor in human colon cancer cells. Biol Pharm Bull. 2006; 29:903–906.
7. Shell S, Park SM, Radjabi AR, Schickel R, Kistner EO, Jewell DA, Feig C, Lengyel E, Peter ME. Let-7 expression defines two differentiation stages of cancer. Proc Natl Acad Sci U S A. 2007; 104:11400–11405.
8. Motoyama K, Inoue H, Nakamura Y, Uetake H, Sugihara K, Mori M. Clinical significance of high mobility group A2 in human gastric cancer and its relationship to let-7 MicroRNA family. Clinical Cancer Research. 2008; 14:2334–2340.
9. Peng Y, Laser J, Shi G, Mittal K, Melamed J, Lee P, Wei JJ. Anti proliferative effects by Let-7 repression of high-mobility group A2 in uterine leiomyoma. Mol Cancer Res. 2008; 6:663–673.
10. Schultz J, Lorenz P, Gross G, Ibrahim S, Kunz M. MicroRNA let-7b targets important cell cycle molecules in malignant melanoma cells and interferes with anchorage-independent growth. Cell Res. 2008; 18:549–557.
11. Yang H, Kong W, He L, Zhao JJ, O’Donnell JD, Wang J, Wenham RM, Coppola D, Kruk PA, Nicosia SV, Cheng JQ. MicroRNA expression profiling in human ovarian cancer: miR-214 induces cell survival and cisplatin resistance by targeting PTEN. Cancer Research. 2008; 68:425–433.
12. Nam EJ, Yoon HJ, Kim SW, Kim HG, Kim YT, Kim JH, Kim JW, Kim SH. MicroRNA expression profiles in serous ovarian carcinoma. Clinical Cancer Research. 2008; 14:2690–2695.
13. Iorio MV, Visone R, Di Leva G, Donati V, Petrocca F, Casalini P, Taccioli C, Volinia S, Liu CG, Alder H, Calin GA, Menard S, Croce CM. MicroRNA signatures in human ovarian cancer. Cancer Research. 2007; 67:8699–8707.
14. Shimizu S, Takehara T, Hikita H, Kodama T, Miyagi T, Hosui A, Tatsumi T, Ishida H, Noda T, Nagano H, Doki Y, Mori M, Hayashi N. The let-7 family of microRNAs inhibits Bcl-xL expression and potentiates sorafenib-induced apoptosis in human hepatocellular carcinoma. J Hepatol. 2010.
15. Wang Z, Lin S, Li JJ, Xu Z, Yao H, Zhu X, Xie D, Shen Z, Sze J, Li K, Lu G, Chan DT, Poon WS, et al. MYC protein inhibits transcription of the microRNA cluster MC-let-7a-1~let-7d via noncanonical E-box. J Biol Chem. 2011; 286:39703–39714.
16. Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ. RAS is regulated by the let-7 microRNA family. Cell. 2005; 120:635–647.
17. Sampson VB, Rong NH, Han J, Yang Q, Aris V, Soteropoulos P, Petrelli NJ, Dunn SP, Krueger LJ. MicroRNA let-7a down-regulates MYC and reverts MYC-induced growth in Burkitt lymphoma cells. Cancer Res. 2007; 67:9762–9770.
18. Lan FF, Wang H, Chen YC, Chan CY, Ng SS, Li K, Xie D, He ML, Lin MC, Kung HF. Hsa-let-7g inhibits proliferation of hepatocellular carcinoma Cells by down-regulation of c-Myc and Up-regulation of p16(INK4A). Int J Cancer. 2010.
19. Wang YC, Chen YL, Yuan RH, Pan HW, Yang WC, Hsu HC, Jeng YM. Lin-28B expression promotes transformation and invasion in human hepatocellular carcinoma. Carcinogenesis. 2010; 31:1516–1522.
20. Sherman M. Epidemiology of hepatocellular carcinoma. Oncology. 2010; 78 Suppl 1:7–10.
21. Omuro A, DeAngelis LM. Glioblastoma and other malignant gliomas: a clinical review. JAMA. 2013; 310:1842–1850.
22. Xia HP, Cheung WKC, Sze J, Lu G, Jiang SS, Yao H, Bian XW, Poon WS, Kung HF, Lin MC. miR-200a Regulates Epithelial-Mesenchymal to Stem-like Transition via ZEB2 and beta-Catenin Signaling. Journal of Biological Chemistry. 2010; 285:36995–37004.
23. Xia H, Ng SS, Jiang S, Cheung WK, Sze J, Bian XW, Kung HF, Lin MC. miR-200a-mediated downregulation of ZEB2 and CTNNB1 differentially inhibits nasopharyngeal carcinoma cell growth, migration and invasion. Biochem Biophys Res Commun. 2010; 391:535–541.
24. Xia H, Qi Y, Ng SS, Chen X, Li D, Chen S, Ge R, Jiang S, Li G, Chen Y, He ML, Kung HF, Lai L, et al. microRNA-146b inhibits glioma cell migration and invasion by targeting MMPs. Brain Res. 2009; 1269:158–165.
25. Xia H, Qi Y, Ng SS, Chen X, Chen S, Fang M, Li D, Zhao Y, Ge R, Li G, Chen Y, He ML, Kung HF, et al. MicroRNA-15b regulates cell cycle progression by targeting cyclins in glioma cells. Biochem Biophys Res Commun. 2009; 380:205–210.
26. Lin S, Cheung WK, Chen S, Lu G, Wang Z, Xie D, Li K, Lin MC, Kung HF. Computational identification and characterization of primate-specific microRNAs in human genome. Comput Biol Chem. 2010; 34:232–241.
27. Ciafre SA, Galardi S, Mangiola A, Ferracin M, Liu CG, Sabatino G, Negrini M, Maira G, Croce CM, Farace MG. Extensive modulation of a set of microRNAs in primary glioblastoma. Biochem Biophys Res Commun. 2005; 334:1351–1358.
28. Lawler S, Godlewski J, Nowicki MO, Bronisz A, Williams S, Otsuki A, Nuovo G, RayChaudhury A, Newton HB, Chiocca EA. Targeting of the Bmi-1 Oncogene/Stem Cell Renewal Factor by MicroRNA-128 Inhibits Glioma Proliferation and Self-Renewal. Cancer Research. 2008; 68:9125–9130.
29. Silber J, Lim DA, Petritsch C, Persson AI, Maunakea AK, Yu M, Vandenberg SR, Ginzinger DG, James CD, Costello JF, Bergers G, Weiss WA, Alvarez-Buylla A, et al. miR-124 and miR-137 inhibit proliferation of glioblastoma multiforme cells and induce differentiation of brain tumor stem cells. BMC Med. 2008; 6:14.
30. Conti A, Aguennouz M, La Torre D, Tomasello C, Cardali S, Angileri FF, Maio F, Cama A, Germano A, Vita G, Tomasello F. miR-21 and 221 upregulation and miR-181b downregulation in human grade II-IV astrocytic tumors. J Neurooncol. 2009; 93:325–332.
31. Malzkorn B, Wolter M, Liesenberg F, Grzendowski M, Stuhler K, Meyer HE, Reifenberger G. Identification and functional characterization of microRNAs involved in the malignant progression of gliomas. Brain Pathol. 2010; 20:539–550.
32. Huse J, Brennan C, Hambardzumyan D, Pena J, Rouhanifard S, Sohn-Lee C, le Sage C, Agami R, Tuschl T, Holland E. The PTEN-Regulating MicroRNA miR-26a is Amplified in High-Grade Glioma and Facilitates Gliomagenesis in vivo. J Neuropath Exp Neur. 2009; 68:558–558.
33. Radlwimmer B, Ernst A, Campos B, Meier J, Devens F, Liesenberg F, Wolter M, Reifenberger G, Herold-Mende C, Lichter P. De-repression of CTGF via the miR-17-92 cluster upon differentiation of human glioblastoma spheroid cultures. Oncogene. 2010; 29:3411–3422.
34. Slaby O, Lakomy R, Fadrus P, Hrstka R, Kren L, Lzicarova E, Smrcka M, Svoboda M, Dolezalova H, Novakova J, Valik D, Vyzula R, Michalek J. MicroRNA-181 family predicts response to concomitant chemoradiotherapy with temozolomide in glioblastoma patients. Neoplasma. 2010; 57:264–269.
35. Zhou X, Ren Y, Moore L, Mei M, You Y, Xu P, Wang B, Wang G, Jia Z, Pu P, Zhang W, Kang C. Downregulation of miR-21 inhibits EGFR pathway and suppresses the growth of human glioblastoma cells independent of PTEN status. Lab Invest. 2010; 90:144–155.
36. Hummel R, Maurer J and Haier J. MicroRNAs in brain tumors: a new diagnostic and therapeutic perspective? Mol Neurobiol. 2011; 44:223–234.
37. Qiu SW, Lin S, Hu D, Feng YM, Tan Y, Peng Y. Interactions of miR-323/miR-326/miR-329 and miR-130a/miR-155/miR-210 as prognostic indicators for clinical outcome of glioblastoma patients. J Transl Med. 2013; 11.
38. Cui JG, Zhao Y, Sethi P, Li YY, Mahta A, Culicchia F, Lukiw WJ. Micro-RNA-128 (miRNA-128) down-regulation in glioblastoma targets ARP5 (ANGPTL6), Bmi- 1 and E2F- 3a, key regulators of brain cell proliferation. J Neurooncol. 2010; 98:297–304.
39. Shi L, Cheng ZH, Zhang JX, Li R, Zhao P, Fu Z, You YP. hsa-mir-181a and hsa-mir-181b function as tumor suppressors in human glioma cells. Brain Research. 2008; 1236:185–193.
40. Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson A, Kampf C, Sjostedt E, Asplund A, Olsson I, Edlund K, Lundberg E, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015; 347:1260419.
41. Shervington A, Cruickshanks N, Wright H, Atkinson-Dell R, Lea R, Roberts G, Shervington L. Glioma: What is the role of c-Myc, hsp90 and telomerase? Mol Cell Biochem. 2006; 283:1–9.
42. Pelengaris S, Khan M. The many faces of c-MYC. Arch Biochem Biophys. 2003; 416:129–136.
43. Kiaris H, Schally AV. Decrease in telomerase activity in U-87MG human glioblastomas after treatment with an antagonist of growth hormone-releasing hormone. P Natl Acad Sci USA. 1999; 96:226–231.
44. Miracco C, De Santi MM, Luzi P, Lalinga AV, Laurini L, De Nisi MC, Angeloni G, Brogi M, Cardone C, Carducci A, Arcuri F, Tosi P, Rubino G, et al. In situ detection of telomeres by fluorescence in situ hybridization and telomerase activity in glioblastoma multiforme: Correlation with p53 status, EGFR, c-myc, MIB1, and Topoisomerase II alpha protein expression. Int J Oncol. 2003; 23:1529–1535.
45. Huang HT, Colella S, Kurrer M, Yonekawa Y, Kleihues P, Ohgaki H. Gene expression profiling of low-grade diffuse astrocytomas by cDNA arrays. Cancer Research. 2000; 60:6868–6874.
46. 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.
47. Zhao C, Sun G, Li S, Lang MF, Yang S, Li W, Shi Y. MicroRNA let-7b regulates neural stem cell proliferation and differentiation by targeting nuclear receptor TLX signaling. Proc Natl Acad Sci U S A. 2010; 107:1876–1881.
48. Kawahara H, Okada Y, Imai T, Iwanami A, Mischel PS, Okano H. Musashi1 cooperates in abnormal cell lineage protein 28 (Lin28)-mediated let-7 family microRNA biogenesis in early neural differentiation. J Biol Chem. 2011; 286:16121–16130.
49. Rybak A, Fuchs H, Smirnova L, Brandt C, Pohl EE, Nitsch R, Wulczyn FG. A feedback loop comprising lin-28 and let- 7 controls pre-let-7 maturation during neural stem-cell commitment. Nat Cell Biol. 2008; 10:987–993.
50. Wulczyn FG, Smirnova L, Rybak A, Brandt C, Kwidzinski E, Ninnemann O, Strehle M, Seiler A, Schumacher S, Nitsch R. Post-transcriptional regulation of the let-7 microRNA during neural cell specification. FASEB J. 2007; 21:415–426.
51. Lee ST, Chu K, Oh HJ, Im WS, Lim JY, Kim SK, Park CK, Jung KH, Lee SK, Kim M, Roh JK. Let-7 microRNA inhibits the proliferation of human glioblastoma cells. J Neurooncol. 2011; 102:19–24.
52. Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, Tramontano A and Bozzoni I. A Long Noncoding RNA Controls Muscle Differentiation by Functioning as a Competing Endogenous RNA. Cell. 2011; 147:358–369.
53. Tay Y, Kats L, Salmena L, Weiss D, Tan SM, Ala U, Karreth F, Poliseno L, Provero P, Di Cunto F, Lieberman J, Rigoutsos I, Pandolfi PP. Coding-Independent Regulation of the Tumor Suppressor PTEN by Competing Endogenous mRNAs. Cell. 2011; 147:344–357.
54. Ng SS, Cheung YT, An XM, Chen YC, Li M, Li GH, Cheung W, Sze J, Lai L, Peng Y, Xia HH, Wong BC, Leung SY, et al. Cell cycle-related kinase: a novel candidate oncogene in human glioblastoma. J Natl Cancer Inst. 2007; 99:936–948.
55. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25:402–408.
56. Chakrabarti SK, James JC, Mirmira RG. Quantitative assessment of gene targeting in vitro and in vivo by the pancreatic transcription factor, Pdx1. Importance of chromatin structure in directing promoter binding. J Biol Chem. 2002; 277:13286–13293.
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.