- Proceedings
- Open access
- Published:
c-Myc and viral cofactor Kaposin B co-operate to elicit angiogenesis through modulating miRNome traits of endothelial cells
BMC Systems Biology volume 10, Article number: S1 (2016)
Abstract
Background
MicroRNAs (miRNAs) have emerged as master regulators of angiogenesis and other cancer-related events. Discovering new angiogenesis-regulating microRNAs (angiomiRs) will eventually help in developing new therapeutic strategies for tumor angiogenesis and cardiovascular diseases. Kaposi’s sarcoma (KS), which is induced by the etiological infectious agent KS-associated herpesvirus (KSHV), is a peculiar neoplasm that expresses both blood and lymphatic endothelial markers and possesses extensive neovasculature. Using KSHV and its proteins as baits will be an efficient way to discover new angiomiRs in endothelial cells. Kaposin B is one of the latent viral genes and is expressed in all KSHV tumor cells. Since Kaposin B is a nuclear protein with no DNA-binding domain, it may regulate gene expression by incorporating itself into a transcription complex.
Results
We demonstrated that c-Myc and Kaposin B form a transcription complex and bind to the miR-221/-222 promoter, thereby affecting their expression and anti-angiogenic ability. By small RNA sequencing (smRNA-Seq), we revealed that 72.1 % (173/240) of Kaposin B up-regulated and 46.5 % (113/243) of Kaposin B down-regulated known miRNAs were regulated by c-Myc. We also found that 77 novel miRNA were up-regulated and 28 novel miRNAs were down-regulated in cells expressing both c-Myc and Kaposin B compared with cells expressing Kaposin B only. The result was confirmed by RNA-IP-seq data.
Conclusions
Our study identifies known and novel c-Myc-regulated microRNAs and reveals that a c-Myc-oriented program is coordinated by Kaposin B in KSHV-infected cells.
Background
Kaposi’s sarcoma-associated herpesvirus (KSHV), also known as human herpesvirus 8 (HHV8), is associated with Kaposi’s sarcoma (KS) [1]. KS is the most common cancer in individuals with AIDS and one of the most common cancers in patients who have immunosuppression [2]. KS typically appears as colored lesions or blotches on the skin, although it can spread to internal organs [3]. This disease has been recognized as a peculiar neoplasm that has a diverse cellular makeup and an extensive neovasculature. In host endothelial cells, KSHV infection alters the expression of mRNAs and microRNAs required for many aspects of tumorigenesis, including angiogenesis [4]. Therefore, KS and KSHV are recognized as an excellent model to study angiogenesis and to discover new angiogenic mechanisms, including coding and non-coding RNAs.
Recently, the effect of KSHV infection on cellular microRNA (miRNAs) expression has been reported [5]. It was found that pre-miRNA signatures delineate the stages of KSHV-associated HUVEC transformation [6], [7]. MiRNAs are small RNAs of 18–24 nucleotides in length that have emerged as master regulators of cancer and key modulators of angiogenesis [8]–[10]. An increasing number of studies have shown that some miRNAs, called angiomiRs (microRNAs that regulate angiogenesis), play a crucial role in regulating angiogenesis [11]. For example, miR-221 and miR-222, two closely related miRNAs encoded in cluster from a genomic region on chromosome X, can modulate the angiogenic properties of human umbilical vein endothelial cells (HUVECs) by targeting c-Kit [11] and endothelial nitric oxide synthase (eNOS) [12]. MiR-221 and miR-222 also regulate lymphatic endothelial cell (LEC) motility by targeting ETS2 and ETS1 transcription factors, respectively [4]. Loss of miR-221 expression demarks the transition from merely immortalized to fully tumorigenic HUVECs [6]. KSHV also induces endothelial cell motility by inducing miR-31 expression while suppressing that of miR-221 and miR-222 [4], [13].
Kaposin B, one of the three Kaposin proteins encoded from the KSHV genome K12 locus, is one of the viral genes expressed during latency and is therefore expressed in all KSHV tumor cells [4]. Kaposin B stabilizes various cytokine (such as IL6 and GM-CSF) mRNA containing AU-rich elements (AREs) via the p38/MK2 pathway. In response to LPS, Kaposin B and MK2 were shown to be exported to the cytoplasm, where mRNA stability is regulated [14]. Kaposin B also enhances the PROX1 mRNA stability during lymphatic reprogramming of vascular endothelial cells [15]. Kaposin B can influence cellular gene expression by regulating promoter activities of host genes: both Kaposin B and KSHV viral latency-associated nuclear antigen (LANA) proteins can down-regulate miR-221 and miR-222 levels by repressing the activity of miR-221/-222 cluster promoter [4]. Since there is no predictable DNA-binding domain on Kaposin B, how this nuclear protein can regulate mRNA and miRNA expression remains unclear.
c-Myc achieves its oncogenic effects by regulating transcription of protein-coding genes as well as microRNA genes such as miR-29b-1/miR-29a [16], [17]. c-Myc is also essential for vasculogenesis and angiogenesis during development and tumor progression [18] via inducing the expression of miR-17 ~ 92 angiogenic miRNA cluster [19]. Revealing the angiomiRs regulated by c-Myc and the underlying regulatory mechanisms will help to further understand c-Myc and endothelial cell biology.
Here, we showed that Kaposin B and c-Myc are in the same transcription complex that directly regulates the miR-221/-222 cluster promoter activity. A c-Myc-oriented circuit is therefore formed in the presence of Kaposin B in KSHV-infected endothelial cells. Furthermore, we also provide a global microRNA signature which is regulated by c-Myc and Kaposin B. We hope our roadmap will help the search and development of new therapeutic targets for virus- or cancer-induced angiogenesis, cancer formation and metastasis.
Methods
Cell culture and KSHV infection
Human primary umbilical vein endothelial cells (HUVECs) were purchased from Clonetics Inc. (Walkersville, Md.) and were cultured as described [4]. HMEC1, an immortalized human microvascular endothelial cell line, was cultured in endothelial cell growth medium MV (C-22020; PromoCell, Heidelberg, Germany). A recombinant virus, rKSHV.219, that expresses the red fluorescent protein (RFP) from the KSHV lytic PAN promoter, the green fluorescent protein (GFP) from the EF-1alpha promoter, and with a gene for puromycin resistance as a selectable marker, was constructed using JSC-1 cells as described previously [20].
Plasmid construction
Plasmids expressing KSHV Kaposin B and miR-221, miR-222, or miR-221/-222 were constructed as described previously [4]. c-Myc expression constructs (pcDNA3-HA-c-Myc and pHR-c-Myc) and knock down construct (pLKO.1) was kindly provided by Prof. Kenneth CW Wu [21]. The full-length miR-221/-222 promoter reporter plasmid was constructed as described previously [4], and primers for cloning miR-221/-222 promoter mutants are listed in Additional file 1.
Immunofluorescence assay (IFA)
Cells were fixed with 4 % paraformaldehyde (Sigma-Aldrich), permeabilized with 0.2 % Triton X-100 (Sigma-Aldrich), and then blocked with PBS containing 1 % bovine serum albumin (Sigma-Aldrich). For Kaposin B staining, cells were incubated with the monoclonal antibody anti-FLAG M2 at a 1:500 dilution for 60 to 120 min at 25 °C followed by incubation with FITC-conjugated goat anti-Mouse IgG (1:500, Jackson ImmunoResearch) for 60 to 120 min at 25 °C. Rhodamine-phalloidin (Molecular Probes, Invitrogen) was used to label actin cytoskeleton. Cell nuclei were counterstained with Hoechst 33342 (Sigma-Aldrich), and examined by fluorescence confocal microscopy (Olympus FV1000).
Transwell cell migration and endothelial cell tube formation assays
Cell migration ability was evaluated using Costar Transwell® Polycarbonate Permeable Supports (Corning, NY, USA) as described previously [13]. In brief, 5 × 104 cells in 500 μl of culture medium were applied to the upper chamber of the device, and 750 μl of medium containing 10 ng/ml human VEGF (R&D Systems, Minneapolis, MN, USA) was added to the lower chamber. A polycarbonate membrane with a pore size of 8 μm was placed in between the two chambers. After 6 h of incubation at 37 °C for HMEC1 and HUVEC cells, the membrane was fixed in 4 % paraformaldehyde (Sigma-Aldrich) for 20 min at room temperature and then stained with Hoechst 33342 solution (Sigma-Aldrich) for 30 min. On the upper side of membrane were identified un-migrated cells and removed. The cells under the membrane were counted under a microscope.
In vitro angiogenic activity of vascular endothelial cells was examined by tube formation assays: primary HUVECs were placed onto MatriGel® Basement Membrane Matrix gel (Becton-Dickinson, Franklin Lakes, NJ, USA) for 6 h at 37 °C. The formation of tubule structure was observed every hour under a microscope.
Small RNA sequencing (smRNA-Seq) and data analysis
Total RNA was collected and small RNA fractions were sequenced using Illumina HiSeq2000 sequencer (Illumina, San Diego, CA USA) according to the manufacturer’s instructions. For data analysis, quality Fastq sequences, which were without poly-A, ambiguous nucleotides or a 5′ adapter, yet flanking 6–18 nt of 3′ adapter sequence, had the adapter sequences trimmed and the identical sequences were then collapsed to unique sequences. The resulting unique sequences that did not align to the mRNA database (UCSC genome browsers) but were aligned to known microRNA sequences (miRBase R20; http://www.mirbase.org/) were subjected to further quantification analysis. Sequencing reads were used to obtain a RPM (reads per million mapped reads) [22], [23] value calculated as C/(MN) × 109, where C is “read numbers aligned to a given miRNA chromosomal region”, M is “multiple mapping numbers across all miRNA regions” and N is “total read numbers that map to human genome sequence”.
RNA extraction, real-time quantitative polymerase chain reaction (qPCR) and microRNA expression
RNA extraction and reverse transcription were performed as previously described [13]. The expression of mature human miR-221 and miR-222 was determined by a stem-loop real-time PCR system using the primer pairs previously published [4]. MicroRNA expression data were normalized to that of U6 snRNA. Details of qPCR primers are in Additional file 1.
Coimmunoprecipitation (co-IP), chromatin immunoprecipitation (ChIP) and Immunoblotting
Co-IP and ChIP assays were performed as previously described [24]. Briefly, Kaposin B or c-Myc expression plasmid was co-transfected into 293T cells using TurboFect™ (R0531; Fermentas, Glen Burnie, Maryland, USA). After 48 h, cells were lyzed in NET lysis buffer containing protease and phosphatase inhibitors and cell lysate (500 μl) was immunoprecipitated with anti-FLAG M2 (for detecting FLAG-tagged Kaposin B protein; Sigma-Aldrich), anti-HA or isotype IgG antibody. The immunoprecipitates were probed with anti-FLAG M2 (1:2000) and anti-c-Myc Tag (clone 9E10; MILLIPORE, 05–419, 1:2000) monoclonal antibodies, and visualized with horseradish peroxidase-conjugated secondary antibodies using the ECL reagent (Amersham Pharmacia Biosciences). For chromatin immunoprecipitation assay, chromosomal DNA fragments were prepared as described [25]. Briefly, lysates were incubated with isotype IgG or anti-FLAG antibody. The approximately 100-base pair products from miR-221/-222 proximal promoter (E1 region: −600 (45,607,128-45,607,214); E2/3 region: −1100 (45,607,635-45,607,764); NC region: −2600 (45,609,170-45,609,265)) were detected by real-time PCR.
Luciferase reporter assays
For reporter assays, 293T cells (5 × 104/well) in 24-well plates were transfected with either Kaposin B or c-Myc expression plasmid, along with the firefly luciferase reporter gene construct (500 ng) and 50 ng pRL-TK Renilla luciferase construct (for normalization) using TurboFect™. Cells were harvested 48 h after transfection, and the luciferase activity was measured using the Dual-Luciferase Reporter Assay System (E1910; Promega, Charbonnières, France).
Results
Kaposin B binds to c-Myc for regulating endothelial cell angiogenic activities
Kaposin B is a nuclear protein which can regulate miR-221 and miR-222 expressions without any predictable DNA-binding domain [4], [26]. Therefore, we hypothesized that Kaposin B may be a novel transcription cofactor which regulates promoter activities by binding to other transcription factor(s). To investigate, we overexpressed Kaposin B in HMEC1 and primary HUVEC cells and tested their migration and tube formation abilities. We verified expression of Kaposin B protein in lentivirus-transduced HUVEC using immunoblotting (Additional file 2-A), and immunofluorescence staining showed that localization of Kaposin B was in nucleus (Fig. 1a). Expression of Kaposin B induced an increase in cell migration rate in HMEC1 and primary HUVEC cells (Fig. 1b-c). Kaposin B also induced more microvascular structures formation by HUVEC in the in vitro tube formation assay (Fig. 1d).
To maintain KSHV latency and oncogenesis, KSHV-latency proteins, LANA and vIRF3, stabilized and activated c-Myc proteins [27]–[30] and c-Myc was shown to reduce miR-222 expression in mouse mammary tumors [31]. Therefore, we examined the proximal promoter region of miR-221/-222 and identified three putative c-Myc binding sites (E1, E2 and E3 E-boxes with a motif sequence CANNTG) (Fig. 1e). We hypothesized that Kaposin B may associate with c-Myc in regulating host gene expression. We therefore tested whether Kaposin B and c-Myc are in the same protein complex by co-transfecting FLAG-tagged Kaposin B and HA-tagged c-Myc into 293T cells. c-Myc was detected in immunoprecipitation (IP) assays using anti-FLAG monoclonal antibody (mAb) against Kaposin B. On the other hand, Kaposin B was detected in IP assay using anti-HA mAb against c-Myc (Fig. 1f). The interaction between c-Myc and Kaposin B was independent of promoter DNA since c-Myc was still detected in the anti-FLAG mAb pull-down products even when cell lysates were pre-treated with DNase digestion for removal of genomic DNA fragments (Fig. 1g, Additional file 2-B).
We further checked whether c-Myc is required for Kaposin B functions. Knockdown of endogenous c-Myc in Kaposin B-expressing HUVEC cells abolished Kaposin B-induced cell migration and microvasculature formation. (Fig. 1h-i, Additional file 2-C). We also knocked down c-Myc in Kaposin B-expressing HMEC1 cells and found that cell migration ability was down-regulated. (Additional file 2-D-E).
Small RNA sequencing reveals known miRNAs co-regulated by c-Myc and Kaposin B in primary endothelial cells
Previous study showed that Kaposin B can down-regulate miR-221 and miR-222 levels by repressing the activity of miR-221/-222 cluster promoter [4]. However, there are few studies exploring Kaposin B-induced global miRNAs alteration. Therefore, we overexpressed Kaposin B in HUVEC to examine its regulated miRNA expression. We have demonstrated that c-Myc binds to Kaposin B (Fig. 1f-g) and is required for its functions (Fig. 1h-i), hence we overexpressed c-Myc in Kaposin B(+) HUVEC to examine its ability to augment Kaposin B-regulated miRNA expression.
We analyzed small RNA sequencing (smRNA-Seq) data according to an in-house bioinformatics pipeline [32]. The Illumina next generation sequencing platform generated more than 29 million high-quality sequencing reads for HUVEC transduced with empty lentiviral vector alone, Kaposin B alone and Kaposin B with c-Myc. The initial operations included identifying sequence matches to the coding RNA database for eliminating reads of degraded mRNAs. Non-coding RNA reads that match previously annotated miRNAs deposited in the miRBase database (release 21) were subjected to normalization and quantitative profiling. Firstly, we explored known miRNAs expression. When we overexpressed Kaposin B alone (KapB/vec), we found 240 up-regulated and 243 down-regulated known miRNAs. In contrast, 312 up-regulated and 200 down-regulated miRNAs (≥1.5 fold change; Fig. 2a) were identified upon simultaneous overexpression of Kaposin B and c-Myc (KapB + c-Myc/vec). We further interrogated the contribution of c-Myc in Kaposin B-regulated miRNome changes. Comparison of the c-Myc(+)Kaposin B(+) miRNA signature with that triggered by Kaposin B alone in HUVEC indicated that c-Myc helped to further induced the expression of 72.1 % (173/240) of Kapsoin B up-regulated miRNAs (Fig. 2b, upper panel, fold change ≥1.5). Similarly, the miRNAs down-regulated by Kaposin B and c-Myc were 46.5 % (113/243); Fig. 2b, lower panel, fold change ≥1.5). c-Myc and Kaposin B-regulated miRNAs are enumerated in Tables 1 and 2 (only top 50 of c-Myc(+)Kaposin B(+) miRNA are listed). Gene Set Enrichment Analysis (GSEA) revealed that the expression of up- or down-regulated miRNAs by c-Myc correlated with those induced or repressed by Kaposin B in HUVEC (Fig. 2d). We also performed qPCR to detect miRNAs expression in Kaposin B(+) and Kaposin B(+)c-Myc(+) cells. We identified that miR-193b-5p, miR-197-5p, miR-210, and miR-1246 were up-regulated in Kaposin B(+) and Kaposin B(+)c-Myc(+) cells. In contrast, miR-100, miR-146a-3p and miR-1271-5p were down-regulated (Fig. 2d). RT-qPCR also demonstrated that miR-221 and miR-222 were down-regulated in KaposinB(+) and Kaposin B(+)c-Myc(+) cells (Fig. 2e). The miRNAs detected in RT-qPCR, except miR-1246 and miR-1271-5p, were similar to the ones detected in KSHV-infected HUVECs (Additional file 3). Moreover, we had used up to 3 different kind of tools (TargetScan (http://www.targetscan.org/vert_70/), mirDB (http://mirdb.org/miRDB/index.html) and TargetMiner (http://www.isical.ac.in/~bioinfo_miu/targetminer20.htm)) to analyze the targeted genes of c-Myc and Kaposin B induced/reduced top 5 miRNAs (Table 1). Based on our results, 806 genes were targeted by the up-regulated miRNAs (Additional file 4-A) and 818 genes were targeted by the down-regulated miRNAs (Additional file 4-B), respectively. The biological pathways were further analyzed, genes targeted by the up-regulated miRNAs contribute to axonal guidance, cardiac β-adrenergic signaling and α-adrenergic signaling; while the genes targeted by the down-regulated miRNAs contribute to Wnt/β signaling, reelin signaling and molecular mechanisms of cancer (Fig. 2f). This indicates Kaposin B + c-myc regulated miRNAs may trigger cells to more tumorigenic stage.
Small RNA sequencing identifies novel miRNAs co-regulated by c-Myc and Kaposin B in primary endothelial cells
To identify novel miRNAs, we first removed mRNA contamination and known miRNA loci. Using three independent bioinformatics algorithms (miRDeep2, mireap and miRanalyzer), 20–26 unique genomic loci yielded at the same time in Kaposin B(+)-, Kaposin B(+)c-Myc(+)-expressing cells and HUVEC control cells (Fig. 3a), and we listed all novel miRNA candidates in (Additional file 5). We found 325 and 374 novel miRNA candidates in Kaposin B-expressing cells and Kaposin B(+)c-Myc(+) cells (KapB + c-Myc/vec) that were up-regulated, respectively. In contrast, we identified 246 and 243 novel miRNA candidates in Kaposin B-expressing cells and down-regulated novel miRNA candidates in Kaposin B(+)c-Myc(+) cells (KapB + c-Myc/vec) that were down-regulated, respectively. To characterize c-Myc regulated novel miRNA candidates in Kaposin B cells, we compared the miRNA profiles between Kaposin B(+)c-Myc(+) cells and Kaposin B(+) cells. The miRNA profiles revealed that out of 325 miRNAs, 172 (52.9 %) were up-regulated and out of 246, 169 (68.7 %) were down-regulated novel mature miRNA candidates correlated with c-Myc expression in Kaposin B(+) cells (Fig. 3b). To investigate the RISC-binding properties of the novel miRNA candidates, we acquired nine public Ago1/2-mediated RNA-immunoprecipitation sequencing (Ago1/2-RNA-IP-seq) datasets from the GEO database and re-analyzed using our in-house pipelines [32]. The Ago1/2-RNA-IP-seq datasets verified 127 out of 172 up-regulated and 93 out of 169 down-regulated novel miRNA candidates (Fig. 3c).
To investigate c-Myc-promoting or -reducing novel miRNAs expression in Kaposin B(+) cells, we found that 77 out of 127 were up-regulated and 28 out of 93 were down-regulated novel mature miRNA candidates that were c-Myc-promoting or -reducing novel miRNAs expression in Kaposin B(+) cells (Additional files 6 and 7). We showed 5 up-regulated or down-regulated novel c-Myc-regulating miRNAs candidates (KaposinB + c-Myc > KapsoinB or KaposinB > KaposinB + c-Myc) that expressed higher read numbers in Kaposin B(+) and Kaposin B(+)c-Myc(+) cells (Fig. 3d) and their detailed information, including precursor locus, strand and mature sequence. Based on RNAfold analysis (http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi), we demonstrated that these novel miRNAs could fold into hairpin secondary structure, which is another miRNA property (Fig. 3e).
c-Myc enhances Kaposin B association with and regulation of miRNA promoters
Previous studies showed that miR-221/-222 promoted angiogenesis and metastasis in various cancers, including gliomas, colon cancers and breast cancers [33]–[37] and that c-Myc reduces miR-222 expression in mouse mammary tumors [31]. Furthermore, miR-221/-222 were crucial in KSHV-mediated endothelial cell motility [4]. Therefore, to investigate whether c-Myc can enhance Kaposin B association with and regulation of miR-221/-222 promoters, we knocked down c-Myc expression and we found that miR-221/-222 level were rescued (Fig. 4a).
Next, we performed chromatin immunoprecipitation (ChIP) assays using an anti-FLAG mAb to pull down Kaposin B and its associated genomic DNA fragments. ChIP data identified the presence of c-Myc-interacting miR-221/-222 promoter fragments containing the E1 (−600) or E2/E3 (−1100) E-boxes in HMEC1 cells stably expressing Kaposin B (Fig. 4b). Overexpressing c-Myc further enhanced the association between Kaposin B and the miR-221/-222 promoter (Fig. 4b). A negative control PCR containing a region with no putative E-box motif (the “−2600” group) showed no off-target enrichment from the same ChIP set (Fig. 4b). Direct involvement of E1 and E2/3 motifs in Kaposin B-c-Myc complex binding was verified by ChIP assays using the site-directed mutagenesis of all the E1 and E2/E3 E-boxes plasmids. A wild type −1600 promoter construct, double mutant construct and Kaposin B (empty pLenti4/V5 vector) were co-transfected into 293T cells, and ChIP assays were performed using anti-FLAG mAb. qPCR on IP-purified DNA fragments showed that these E-box motifs were essential for Kaposin B to associate with miR-221/-222 promoter (Fig. 4c). We next examined the critical role of endogenous c-Myc in Kaposin B function. Knockdown of endogenous c-Myc in Kaposin B(+) HMEC1 cells significantly reduced Kaposin B and miR-221/-222 promoter association (Fig. 4d).
Repression of miR-221/-222 promoter activity was observed after transient co-transfection of Kaposin B- or c-Myc-expressing plasmids with the full-length reporter construct into 293T cells (Fig. 4f, the “−1600 bp” group). Ectopic expression of both Kaposin B and c-Myc repressed miR-221/-222 promoter activity significantly more than expression of Kaposin B or c-Myc alone (Fig. 4e), suggesting c-Myc enhances the capability of Kaposin B to regulate the miR-221/-222 promoter. Deleting the region containing E2/3 E-boxes did not dramatically influence Kaposin B or c-Myc function (Fig. 4e, the “−600” group), but additional deletion of the E1 E-box significantly abolished the repression of the miR-221/-222 promoter by Kaposin B, c-Myc or both Kaposin B and c-Myc together (Fig. 4e, the “−200”, “−150” and “−50” group). Site-directed mutagenesis of all the E1 and E2/E3 E-boxes, but not E1 or E2/3 alone, hampered Kaposin B- and/or c-Myc-mediated repression, indicating c-Myc could bind to all the 3 E-boxes to direct Kaposin B to the miR-221/-222 promoter (Fig. 4f). Figure 5 is an illustration of the mechanism whereby Kaposin B and c-Myc in the same transcription complex can directly regulate the miR-221/-222 cluster promoter activity.
Discussion
Kaposin B is known to stabilize various AU-rich elements (ARE)-containing cytokine mRNAs, such as IL6, via the p38/MK2 pathway [14], [38]. AREs are usually 50 to 150 nucleotides in length and located in the 3′ untranslated region (UTR) of short half-life mRNAs. AREs serve as an mRNA-destabilizing determinant by promoting degradation of mRNAs. AREs are now grouped into three classes based on the number and distribution of the core AUUUA pentamers: cluster I ARE genes, such as cell cycle regulators c-Myc and cyclin D1, contain several dispersed copies of the AUUUA motif within the AU-rich region [39]. Cluster II ARE genes (such as cytokines IL1B, IL6, IL8 and granulocyte/macrophage colony-stimulating factor (GM-CSF)) have at least 2 overlapping UUAUUUA(U/A)(U/A) motifs, and cluster III ARE genes (including c-JUN and p53) are less characterized and do not contain the canonical AUUUA motif [39]. Database search revealed that there is no ARE on miR-221-222 pri-miRNA. Whether Kaposin B collaborates with c-Myc to regulate promoter activities of these ARE(−) genes is an interesting direction to follow and can be answered by global genomics approaches, such as ChIP-Seq.
In this study, we use miR-221/-222 as a model to elucidate how Kaposin B regulates promoter activity. We demonstrate that c-Myc together with Kaposin B regulate the activity of miR-221/-222 cluster promoter. c-Myc is known to maintain KSHV latency [40]. KSHV latency proteins, such as LANA and vIRF3, have been demonstrated to regulate c-Myc. LANA stabilizes and activates c-Myc through inhibition of GSK-3 phosphorylation and c-Myc ubiquitination [27], [28]. LANA protein that lacks a known DNA-binding domain, also represses the expression of miR-221/-222 cluster [4]. Nevertheless, pilot studies showed that the interaction of Kaposin B and c-Myc is independent of LANA (not shown). vIRF3 associated with c-Myc suppressor, MM-1α, to stimulate c-Myc function [29]. A c-Myc-oriented angiogenesis program is therefore activated and plays a crucial role in KSHV-infected cells and even in KS. Since KS is a neoplasm with an extensive neovasculature, the induction of endothelial cell migration and microvascular formation by c-Myc and Kaposin B may eventually benefit KS formation and progression.
MYC protein, including c-Myc, N-Myc and L-Myc, have been shown to be involved in pathogenesis of most cancers, such as burkitt’s lymphoma, breast cancer, prostate cancer, gastrointestinal cancer etc. [41]. c-Myc affects most important oncogenic programs through regulating several genes and miRNAs expression [16], [42], [43]. In c-Myc-induced mouse mammary tumors, 50 and 59 miRNAs show increased and decreased expression, respectively. The expression of miR-20a, −20b, and miR-9 were increased, while that of miR-222 was decreased [16], [31]. In angiogenesis, previous report showed that c-Myc activated miR-17 cluster expression to attenuate clusterin expression, thereby promoting angiogenesis [44]. In this report, we provide a comprehensive c-Myc-regulated miRome trait in primary human endothelial cells. Most of the c-Myc and/or Kaposin B-affected miRNAs have not been linked to angiogenesis, cell motility or tumorigenesis. Among miRNAs induced by both c-Myc and Kaposin B (Table 1), miR-210, the most prominent hypoxia-induced microRNA, modulates cellular response to hypoxia through MNT, a MAX-interacting transcriptional repressor that functions as a c-Myc antagonist [45]. miR-483-5p, a microRNA embedded in the intron of insulin-like growth factor 2 (Igf2), inhibits angiogenesis via targeting serum response factor [46]. miR-372 is involved in cell migration, proliferation, apoptosis, and globally affect the regulation circuits centered on MAPK/ERK signaling in colon cancer cells [47]. Knockdown of miR-139 in HUVEC significantly suppressed cell viability [48]. miR-423 promotes cell growth and regulates G(1)/S transition by targeting p21Cip1/Waf1 in hepatocellular carcinoma [49].
Conclusions
In this study, we provide a global picture of KSHV- and/or Kaposin B-induced microRNAs, especially those co-regulated by c-Myc. We have demonstrated that c-Myc is essential for Kaposin B to stimulate endothelial cell activity and to activate miRNA promoters. In addition to providing a global picture, this study also benefits the development of new therapeutic approaches to treat KSHV-induced tumor progression in future.
Additional files
References
Mesri EA, Cesarman E, Boshoff C: Kaposi’s sarcoma and its associated herpesvirus. Nat Rev Cancer. 2010, 10 (10): 707-719. 10.1038/nrc2888.
Hengge UR, Ruzicka T, Tyring SK, Stuschke M, Roggendorf M, Schwartz RA, et al: Update on Kaposi’s sarcoma and other HHV8 associated diseases. Part 2: pathogenesis, Castleman’s disease, and pleural effusion lymphoma. Lancet Infect Dis. 2002, 2 (6): 344-352. 10.1016/S1473-3099(02)00288-8.
Hengge UR, Ruzicka T, Tyring SK, Stuschke M, Roggendorf M, Schwartz RA, et al: Update on Kaposi’s sarcoma and other HHV8 associated diseases. Part 1: epidemiology, environmental predispositions, clinical manifestations, and therapy. Lancet Infect Dis. 2002, 2 (5): 281-292. 10.1016/S1473-3099(02)00263-3.
Wu YH, Hu TF, Chen YC, Tsai YN, Tsai YH, Cheng CC, et al: The manipulation of miRNA-gene regulatory networks by KSHV induces endothelial cell motility. Blood. 2011, 118 (10): 2896-2905. 10.1182/blood-2011-01-330589.
Lagos D, Pollara G, Henderson S, Gratrix F, Fabani M, Milne RS, et al: miR-132 regulates antiviral innate immunity through suppression of the p300 transcriptional co-activator. Nat Cell Biol. 2010, 12 (5): 513-519. 10.1038/ncb2054.
O’Hara AJ, Chugh P, Wang L, Netto EM, Luz E, Harrington WJ, et al: Pre-micro RNA signatures delineate stages of endothelial cell transformation in Kaposi sarcoma. PLoS Pathog. 2009, 5 (4): e1000389-10.1371/journal.ppat.1000389.
O’Hara AJ, Wang L, Dezube BJ, Harrington WJ, Damania B, Dittmer DP: Tumor suppressor microRNAs are underrepresented in primary effusion lymphoma and Kaposi sarcoma. Blood. 2009, 113 (23): 5938-5941. 10.1182/blood-2008-09-179168.
Croce CM, Calin GA: miRNAs, cancer, and stem cell division. Cell. 2005, 122 (1): 6-7. 10.1016/j.cell.2005.06.036.
Schickel R, Boyerinas B, Park SM, Peter ME: MicroRNAs: key players in the immune system, differentiation, tumorigenesis and cell death. Oncogene. 2008, 27 (45): 5959-5974. 10.1038/onc.2008.274.
Zimmerman AL, Wu S: MicroRNAs, cancer and cancer stem cells. Cancer Lett. 2011, 300 (1): 10-19. 10.1016/j.canlet.2010.09.019.
Poliseno L, Tuccoli A, Mariani L, Evangelista M, Citti L, Woods K, et al: MicroRNAs modulate the angiogenic properties of HUVECs. Blood. 2006, 108 (9): 3068-3071. 10.1182/blood-2006-01-012369.
Suarez Y, Fernandez-Hernando C, Pober JS, Sessa WC: Dicer dependent microRNAs regulate gene expression and functions in human endothelial cells. Circ Res. 2007, 100 (8): 1164-1173. 10.1161/01.RES.0000265065.26744.17.
Tsai YH, Wu MF, Wu YH, Chang SJ, Lin SF, Sharp TV, et al: The M type K15 protein of Kaposi’s sarcoma-associated herpesvirus regulates microRNA expression via its SH2-binding motif to induce cell migration and invasion. J Virol. 2009, 83 (2): 622-632. 10.1128/JVI.00869-08.
McCormick C, Ganem D: The kaposin B protein of KSHV activates the p38/MK2 pathway and stabilizes cytokine mRNAs. Science. 2005, 307 (5710): 739-741. 10.1126/science.1105779.
Yoo J, Kang J, Lee HN, Aguilar B, Kafka D, Lee S, et al: Kaposin-B enhances the PROX1 mRNA stability during lymphatic reprogramming of vascular endothelial cells by Kaposi’s sarcoma herpes virus. PLoS Pathog. 2010, 6 (8): e1001046-10.1371/journal.ppat.1001046.
Bui TV, Mendell JT: Myc: Maestro of MicroRNAs. Genes Cancer. 2010, 1 (6): 568-575. 10.1177/1947601910377491.
Mott JL, Kurita S, Cazanave SC, Bronk SF, Werneburg NW, Fernandez-Zapico ME: Transcriptional suppression of mir-29b-1/mir-29a promoter by c-Myc, hedgehog, and NF-kappaB. J Cell Biochem. 2010, 110 (5): 1155-1164. 10.1002/jcb.22630.
Baudino TA, McKay C, Pendeville-Samain H, Nilsson JA, Maclean KH, White EL, et al: c-Myc is essential for vasculogenesis and angiogenesis during development and tumor progression. Genes Dev. 2002, 16 (19): 2530-2543. 10.1101/gad.1024602.
Dews M, Homayouni A, Yu D, Murphy D, Sevignani C, Wentzel E, et al: Augmentation of tumor angiogenesis by a Myc-activated microRNA cluster. Nat Genet. 2006, 38 (9): 1060-1065. 10.1038/ng1855.
Vieira J, O’Hearn PM: Use of the red fluorescent protein as a marker of Kaposi’s sarcoma-associated herpesvirus lytic gene expression. Virology. 2004, 325 (2): 225-240. 10.1016/j.virol.2004.03.049.
Chou YT, Lin HH, Lien YC, Wang YH, Hong CF, Kao YR, et al: EGFR promotes lung tumorigenesis by activating miR-7 through a Ras/ERK/Myc pathway that targets the Ets2 transcriptional repressor ERF. Cancer Res. 2010, 70 (21): 8822-8831. 10.1158/0008-5472.CAN-10-0638.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.
Hsieh JY, Huang TS, Cheng SM, Lin WS, Tsai TN, Lee OK, et al: miR-146a-5p circuitry uncouples cell proliferation and migration, but not differentiation, in human mesenchymal stem cells. Nucleic Acids Res. 2013, 41 (21): 9753-9763. 10.1093/nar/gkt666.
Hwang WL, Yang MH, Tsai ML, Lan HY, Su SH, Chang SC, et al: SNAIL regulates interleukin-8 expression, stem cell-like activity, and tumorigenicity of human colorectal carcinoma cells. Gastroenterology. 2011, 141 (1): 279-291.
Yang MH, Wu MZ, Chiou SH, Chen PM, Chang SY, Liu CJ, et al: Direct regulation of TWIST by HIF-1alpha promotes metastasis. Nat Cell Biol. 2008, 10 (3): 295-305. 10.1038/ncb1691.
Sadler R, Wu L, Forghani B, Renne R, Zhong W, Herndier B, et al: A complex translational program generates multiple novel proteins from the latently expressed kaposin (K12) locus of Kaposi’s sarcoma-associated herpesvirus. J Virol. 1999, 73 (7): 5722-5730.
Bubman D, Guasparri I, Cesarman E: Deregulation of c-Myc in primary effusion lymphoma by Kaposi’s sarcoma herpesvirus latency-associated nuclear antigen. Oncogene. 2007, 26 (34): 4979-4986. 10.1038/sj.onc.1210299.
Liu J, Martin HJ, Liao G, Hayward SD: The Kaposi’s sarcoma-associated herpesvirus LANA protein stabilizes and activates c-Myc. J Virol. 2007, 81 (19): 10451-10459. 10.1128/JVI.00804-07.
Lubyova B, Kellum MJ, Frisancho JA, Pitha PM: Stimulation of c-Myc transcriptional activity by vIRF-3 of Kaposi sarcoma-associated herpesvirus. J Biol Chem. 2007, 282 (44): 31944-31953. 10.1074/jbc.M706430200.
Nador RG, Cesarman E, Chadburn A, Dawson DB, Ansari MQ, Sald J, et al: Primary effusion lymphoma: a distinct clinicopathologic entity associated with the Kaposi’s sarcoma-associated herpes virus. Blood. 1996, 88 (2): 645-656.
Sun Y, Wu J, Wu SH, Thakur A, Bollig A, Huang Y, et al: Expression profile of microRNAs in c-Myc induced mouse mammary tumors. Breast Cancer Res Treat. 2009, 118 (1): 185-196. 10.1007/s10549-008-0171-6.
Cheng WC, Chung IF, Tsai CF, Huang TS, Chen CY, Wang SC, et al: YM500v2: a small RNA sequencing (smRNA-seq) database for human cancer miRNome research. Nucleic Acids Res. 2015, 43 (Database issue): D862-867. 10.1093/nar/gku1156.
Shah MY, Calin GA: MicroRNAs miR-221 and miR-222: a new level of regulation in aggressive breast cancer. Genome Med. 2011, 3 (8): 56-10.1186/gm272.
Ye X, Bai W, Zhu H, Zhang X, Chen Y, Wang L, et al: MiR-221 promotes trastuzumab-resistance and metastasis in HER2-positive breast cancers by targeting PTEN. BMB Rep. 2014, 47 (5): 268-273. 10.5483/BMBRep.2014.47.5.165.
Qin J, Luo M: MicroRNA-221 promotes colorectal cancer cell invasion and metastasis by targeting RECK. FEBS Lett. 2014, 588 (1): 99-104. 10.1016/j.febslet.2013.11.014.
Yang F, Wang W, Zhou C, Xi W, Yuan L, Chen X, et al: MiR-221/-222 promote human glioma cell invasion and angiogenesis by targeting TIMP2. Tumour Biol. 2015, 36 (5): 3763-3773. 10.1007/s13277-014-3017-3.
Hur K, Toiyama Y, Schetter AJ, Okugawa Y, Harris CC, Boland CR, et al: Identification of a metastasis-specific MicroRNA signature in human colorectal cancer. J Natl Cancer Inst. 2015, 107 (3): dju492-10.1093/jnci/dju492.
Wang HW, Boshoff C: Linking Kaposi virus to cancer-associated cytokines. Trends Mol Med. 2005, 11 (7): 309-312. 10.1016/j.molmed.2005.06.002.
Barreau C, Paillard L, Osborne HB: AU-rich elements and associated factors: are there unifying principles?. Nucleic Acids Res. 2005, 33 (22): 7138-7150. 10.1093/nar/gki1012.
Li X, Chen S, Feng J, Deng H, Sun R: Myc is required for the maintenance of Kaposi’s sarcoma-associated herpesvirus latency. J Virol. 2010, 84 (17): 8945-8948. 10.1128/JVI.00244-10.
Nesbit CE, Tersak JM, Prochownik EV: MYC oncogenes and human neoplastic disease. Oncogene. 1999, 18 (19): 3004-3016. 10.1038/sj.onc.1202746.
Dang CV: MYC on the path to cancer. Cell. 2012, 149 (1): 22-35. 10.1016/j.cell.2012.03.003.
Meyer N, Penn LZ: Reflecting on 25 years with MYC. Nat Rev Cancer. 2008, 8 (12): 976-990. 10.1038/nrc2231.
Dews M, Fox JL, Hultine S, Sundaram P, Wang W, Liu YY, et al: The myc-miR-17 ~ 92 axis blunts TGF{beta} signaling and production of multiple TGF{beta}-dependent antiangiogenic factors. Cancer Res. 2010, 70 (20): 8233-8246. 10.1158/0008-5472.CAN-10-2412.
Zhang Z, Sun H, Dai H, Walsh R, Imakura M, Schelter J, et al: MicroRNA miR-210 modulates cellular response to hypoxia through the MYC antagonist MNT. Cell Cycle. 2009, 8 (17): 2756-2768. 10.4161/cc.8.17.9387.
Qiao Y, Ma N, Wang X, Hui Y, Li F, Xiang Y, et al: MiR-483-5p controls angiogenesis in vitro and targets serum response factor. FEBS Lett. 2011, 585 (19): 3095-3100. 10.1016/j.febslet.2011.08.039.
Ragusa M, Statello L, Maugeri M, Majorana A, Barbagallo D, Salito L, et al: Specific alterations of the microRNA transcriptome and global network structure in colorectal cancer after treatment with MAPK/ERK inhibitors. J Mol Med. 2012, 90 (12): 1421-1438. 10.1007/s00109-012-0918-8.
Li YX, Liu DQ, Zheng C, Zheng SQ, Liu M, Li X, et al: miR-200a modulate HUVECs viability and migration. IUBMB Life. 2011, 63 (7): 553-559. 10.1002/iub.486.
Lin J, Huang S, Wu S, Ding J, Zhao Y, Liang L, et al: MicroRNA-423 promotes cell growth and regulates G(1)/S transition by targeting p21Cip1/Waf1 in hepatocellular carcinoma. Carcinogenesis. 2011, 32 (11): 1641-1647. 10.1093/carcin/bgr199.
Acknowledgements
The authors acknowledge the technical services provided by the Microarray & Gene Expression Analysis Core Facility, which is supported by the Ministry of Science and Technology (MOST), of the National Yang-Ming University VGH-YM Genome Research Center. This work is funded by the Ministry of Science and Technology (103-2911-I-010-506, 102-2314-B-010-045, and 101-2320-B-010-059-MY3), Mackay Memorial Hospital (MMH-YM-10304), Veterans General Hospitals University System of Taiwan (VGHUST) Joint Research Program, Tsou’s Foundation (VGHUST102-G7-3-2), National Health Research Institutes (10401-62-014), National Yang-Ming University (Ministry of Education, Aim for the Top University Plan) and Tri-Service General Hospital (TSGH-C104-028). This work is also funded in part by the UST-UCSD International Center for Excellence in Advanced Bioengineering sponsored by the Taiwan MOST I-RiCE Program (102-2911-I-009-101). The authors would like to dedicate this paper to the memory of Dr. Hsei-Wei Wang (Apr. 1963-Sep. 2015), who passed away during the period of this research. Dr. Hsei-Wei Wang focused on bioinformatics and biology studying; without his long-term devotion to KSHV research, this paper could not have been completed.
Declarations
The publication cost for this article was funded by Tri-Service General Hospital, grant number (TSGH-C104-028).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they no competing interests.
Authors’ contributions
HWW and CCC designed the analysis approach; HCC, THH, CFT and YNT performed the research and analyzed the data; HWW provided biologic guidance during the analysis process; HWW and HCC wrote the manuscript; and all authors read and approved the final manuscript.
Electronic supplementary material
12918_2015_242_MOESM2_ESM.pdf
Additional file 2: (A) Primary HUVEC stably transduced with Kaposin B or the vector control by lentivirus. Immunoblotting of Kaposin B proteins with anti-Flag mAb (bottom panel). β-actin was used as an internal control. (B) DNase treatment effects were verified by agarose gel electrophoresis. (C-D) knockdown of endogenous c-Myc levels with shRNA in Kaposin B(+) HUVEC (C) and Kaposin B(+) HMEC1 cells (D) were validated by immunoblotting and qRT-PCR. (E) knockdown of endogenous c-Myc levels with shRNA in Kaposin B(+) HMEC1 cells repressed Kaposin B-induced cellular migration (n = 3). (PDF 311 kb) (PDF 312 KB)
12918_2015_242_MOESM3_ESM.pdf
Additional file 3: Validation of smRNA-Seq data by RT-qPCR in KSHV-infected cells. Mean expression levels of the target miRNAs are compared with the U6 control (n = 3). *: P < .05. (PDF 76 kb) (PDF 76 KB)
12918_2015_242_MOESM4_ESM.pdf
Additional file 4: Genes were targeted by the up-regulated miRNAs (A) and down-regulated miRNAs (B). The genes interactions were analyzed by STRING (http://string-db.org/newstring_cgi/show_input_page.pl?UserId=PZAhu49PCqxs&sessionId=A8hUzxbU17ST). (PDF 8481 kb) (PDF 8 MB)
12918_2015_242_MOESM6_ESM.pdf
Additional file 6: List of Kaposin B and c-Myc down-regulated novel miRNAs which exist in Ago1/2-IP-seq data. (PDF 79 kb) (PDF 80 KB)
12918_2015_242_MOESM7_ESM.pdf
Additional file 7: List of Kaposin B and c-Myc up-regulated novel miRNAs which exist in Ago1/2-IP-seq data. (PDF 83 kb) (PDF 84 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chang, HC., Hsieh, TH., Lee, YW. et al. c-Myc and viral cofactor Kaposin B co-operate to elicit angiogenesis through modulating miRNome traits of endothelial cells. BMC Syst Biol 10 (Suppl 1), S1 (2016). https://doi.org/10.1186/s12918-015-0242-3
Published:
DOI: https://doi.org/10.1186/s12918-015-0242-3