Previous Article | Next Article ![]()
Journal of Virology, December 2007, p. 12836-12845, Vol. 81, No. 23
0022-538X/07/$08.00+0 doi:10.1128/JVI.01804-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Department of Molecular Genetics and Microbiology, and Shands Cancer Center, University of Florida, Gainesville, Florida
Received 16 August 2007/ Accepted 11 September 2007
|
|
|---|
|
|
|---|
Investigating miRNA targets is complicated by the fact that miRNAs require only limited complementarity for 3'UTR binding (5, 6, 39). Consequently, a single miRNA can modulate the expression of multiple genes. Analysis of the entire human transcriptome in silico in combination with experimental verification has revealed that the binding of an miRNA to a specific 3'UTR is critically dependent on 5' nucleotides (nt) 2 to 7 of the mature miRNA. This sequence, termed the "seed," classifies miRNAs into families (38).
We hypothesized that virally encoded miRNAs, like virokines, may have been captured from the host genome and/or coevolved within herpesvirus genomes. Since the seed represents the most important determinant for biological activity, we aligned seed sequences of all human herpesvirus-encoded miRNAs with the human miRNA database. Here, we report on a KSHV-encoded miRNA, which has an identical seed sequence to that of hsa-miR-155, and provide evidence that both miRNAs can regulate a common set of cellular genes.
|
|
|---|
3'UTR sequences were obtained from the Ensembl database (www.ensembl.org) or the NCBI and were PCR amplified (TripleMaster; Eppendorf) from either BCBL-1 or 293 genomic DNA (DNAzol; Molecular Research Center, Inc., OH). Primers were designed using Vector NTI (Invitrogen) and are indicated in Table S2 in the supplemental material. PCR products were cloned into pCRII-TOPO (Invitrogen), excised, and inserted into the 3'UTR of pGL3-Promoter. Constructs were confirmed by restriction enzyme digestion and/or sequencing.
For the construction of pmiR-K12-11 and pmiR-155,
180 nt encompassing the stem-loop pre-miRNA were PCR amplified from BCBL-1 genomic DNA using primers indicated in Table S2 in the supplemental material. PCR products were TOPO cloned, excised with HindIII and XhoI, and inserted into pcDNA3.1V5/HisA (Invitrogen) at corresponding sites. KSHV miRNA expression vectors were described previously by Samols et al. (57).
2' O-methylated (2'OMe) RNA oligonucleotides were synthesized by Dharmacon, Inc., and are antisense to the mature miRNA sequence. All bases were modified at the 2' position.
QuikChange site-directed mutagenesis was performed using the primers indicated in Table S3 in the supplemental material according to the manufacturer's protocols (Stratagene). Primer design was done using PrimerX (www.bioinformatics.com). Each 6-mer seed match site was changed to an XhoI restriction enzyme site, and mutants were analyzed by restriction enzyme digestion.
Cell lines and transfections. BCP-1 (a gift from Denise Whitbey at the NCI), BC-1, VG-1, JSC-1 (gifts from Dirk Dittmer at the University of North Carolina), RAJI (a gift from Sankar Swaminathan at the University of Florida), BCBL-1, and BJAB cell lines were grown in RPMI 1640 supplemented with 10% fetal bovine serum and 5% penicillin-streptomycin (Gibco). 293 cells were maintained in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum and 5% penicillin-streptomycin. For stable cell lines, cells were transfected with 10 µg plasmid DNA using Lipofectamine 2000 (Invitrogen) and selected with 500 µg/ml G418 (Mediatech, Inc.) for 4 weeks. Transfections for luciferase assays were performed in 6-well or 24-well plates with Lipofectamine 2000 according to the manufacturer's protocols. Cells were cultured 48 to 72 h prior to harvest. Transfections for BJAB cells were done by nucleofection using cell line solution T, program O-17, according to the manufacturer's instructions (Amaxa, Inc.).
Immunoblotting. Whole-cell lysates were harvested 72 h posttransfection. Lysates corresponding to 1 x 105 cells were separated on 10% sodium dodecyl sulfate-polyacrylamide gels and transferred onto Immobilon-P membranes (Millipore). Blots were probed for BACH-1 or actin (Santa-Cruz Biotechnology, Inc.) and developed with peroxidase-conjugated antibodies and chemiluminescent substrate (Pierce).
Affymetrix array-based gene expression profiling. RNA isolation was performed using the RNeasy kit according to the manufacturer's instructions (QIAGEN). RNA labeling was done using the GeneChip eukaryotic one-cycle target labeling assay as directed by Affymetrix and used to probe Affymetrix U133 2.0 Plus Human GeneChips. Hybridization was performed at 45°C for 16 h, and the chips were then washed and stained using Affymetrix protocol EukGE-WS2v5_450. Array scanning and analysis were performed as described previously (57).
Luciferase assays. Luciferase activity was quantified using the Luciferase assay system (Promega) according to the manufacturer's protocols. Briefly, transfected 293 cells were lysed in cell culture lysis reagent (Promega), and 20% of each cell lysate was assayed for firefly luciferase activity. Light units were normalized to total protein, determined using the BCA protein assay kit (Pierce) according to the manufacturer's instructions, or Renilla luciferase, using a dual luciferase reporter kit (Promega).
RNA extraction for RT-PCR. Total RNA was prepared using RNA-Bee (Tel-Test, Inc., TX) according to the manufacturer's protocols. For detection of BIC expression, 1 µg of DNase-treated RNA was reverse transcribed using SuperScript III reverse transcriptase (Invitrogen) and random hexamers in the presence of RNase-OUT recombinant RNase inhibitor (Invitrogen) to create a cDNA pool. Ten percent of each reverse transcription (RT) reaction was used for PCR amplification. BIC and ß-actin primers are listed in Table S2 in the supplemental material.
Northern blot analysis.
Thirty micrograms of total RNA was loaded onto 15% 8 M urea polyacrylamide gel and transferred onto Genescreen Plus (Perkin-Elmer) following electrophoresis. Probe labeling was performed using T4 polynucleotide kinase (New England Biolabs) in the presence of [
-32P-
]ATP.
Microarray data accession number. The Gene Expression Omnibus accession number for the microarray data is GSE9264.
|
|
|---|
![]() View larger version (22K): [in a new window] |
FIG. 1. Herpesvirus miRNAs share seed sequence homology with human miRNAs. (A) miRNA sequences from KSHV, EBV, and human cytomegalovirus were aligned with sequences in the human miRNA database. Shown are those viral miRNAs exhibiting the greatest sequence homology to human miRNAs. (B) miR-155 is expressed from exon 3 of the non-protein-coding BIC. Primers (forward [fwd] and reverse [rev]) for RT-PCR analysis of BIC expression are indicated. (C) RT-PCR analysis of BIC expression in KSHV-infected PEL cell lines. The RAJI cell line, an EBV-infected BL cell line, was used as a positive control. H2O indicates the no-template control. (D) Northern blot analysis of KSHV-miR-K12-11 and hsa-miR-155 expression in PEL and BL cells. Twenty-five micrograms of total RNA was loaded per lane and hybridized to probes for either miR-K12-11 or miR-155.
|
miR-155 is not expressed in primary effusion lymphomas, which express high levels of miR-K12-11. Based on the identical seed sequences of miR-155 and miR-K12-11 and the association of KSHV with primary effusion lymphoma (PEL) and multicentric Castleman's disease (11), we hypothesized that miR-K12-11 and miR-155 might regulate a common set of cellular genes, thereby contributing to pathogenesis. First, we analyzed BIC expression in latently KSHV-infected PEL-derived cell lines by RT-PCR. EBV-positive Burkitt's lymphoma (BL) cells (RAJI) expressed readily detectable levels of BIC (Fig. 1C). In contrast, all KSHV-positive PEL lines tested did not express BIC (Fig. 1C) and did not express miR-155, as determined by Northern blot analysis (Fig. 1D). PEL cells do, however, highly express miR-K12-11, as shown in Fig. 1D and as previously reported (8, 51, 56). Hence, we hypothesized that miR-K12-11 could mimic miR-155-dependent regulation in PELs.
Generation of miR-155 and miR-K11 expression and sensor vectors. To test this idea, we generated pmiR-155 and pmiR-K12-11 expression vectors (Fig. 2A). To monitor miRNA expression from these constructs, miRNA sensor vectors (pGL3-155 and pGL3-K12-11) containing two antisense complementary binding sites within the 3'UTR of a luciferase reporter (Fig. 2B) were generated. Cotransfection of each miRNA expression vector with the respective sensor showed a 70% to 80% knockdown of luciferase expression (Fig. 2C and D, middle bars). Note that this inhibition reflects a small interfering RNA-like transcript-mediated effect since the miRNA and target 3'UTR are 100% complementary.
![]() View larger version (14K): [in a new window] |
FIG. 2. Importance of the seed sequence for miRNA targeting. (A) Schematic of miRNA expression vector. A region encompassing the pre-miRNA for either miR-155 or miR-K12-11 was PCR amplified and cloned downstream of the cytomegalovirus promoter in pcDNA3.1. (B) Schematic of miRNA sensor vector. Two antisense complementary binding sites for each miRNA were inserted into the 3'UTR of pGL3-Promoter (Promega) downstream of the luciferase gene. (C and D) Luciferase expression from miRNA sensors (pGL3-155 or pGL3-K12-11) is down-regulated in response to ectopic miRNA expression. 293 cells were transfected with 40 ng of the indicated sensor and 800 ng miRNA expression vector using Lipofectamine 2000 (Invitrogen). Lysates were analyzed at 72 h posttransfection (Promega). Relative light units (RLU) are normalized to total protein determined by BCA (Pierce). * indicates a P value of <0.01, and ** indicates a P value of <0.01 by Student's t test. (E) Luciferase assays to control for ectopic miRNA effects. miRNA expression and sensor vectors for KSHV-miR-K12-3-5p were cotransfected with plasmids as indicated in C and D.
|
miR-155 and miR-K12-11 can target the BACH-1 3'UTR. For miR-155, target genes associated with hematopoietic stem cell differentiation (19) and the angiotensin II type I receptor have been reported (42). No targets have been identified for KSHV miR-K12-11. Using three target prediction algorithms (miRanda, PicTar, and TargetScan) (17, 30, 35, 38), we identified BACH-1 (Btb and CNC homolog 1) as a candidate target for both miR-155 and miR-K12-11, as it contains four perfect seed match sites within its 3,315-nt-long 3'UTR (Fig. 3A). BACH-1 is a transcriptional regulator that binds NF-E2 sites and coordinates transcription with small Maf (musculoaponeurotic fibrosarcoma oncogene homolog) proteins (26a).
![]() View larger version (42K): [in a new window] |
FIG. 3. BACH-1 is targeted by miR-K12-11 and miR-155. (A) The BACH-1 3'UTR contains four seed match sites (highlighted) based on three target prediction programs (miRanda, PicTar, and TargetScan). (B and C) A region of the BACH-1 3'UTR (nt 658 to 2495) encompassing the predicted miRNA binding sites was cloned downstream of luciferase (pGL3-BACH1). 293 cells were cotransfected with 40 ng pGL3-BACH-1 and increasing amounts (100, 200, 400, and 800 ng) of pmiR-155 or pmiR-K12-11. Ratios on the x axis indicate the amount of miRNA vector to reporter. pcDNA3.1 was used as filler. Relative light units (RLU) are normalized to total protein determined by BCA (Pierce). (D) Additional KSHV miRNAs do not target the BACH-1 3'UTR. 293 cells were cotransfected with 40 ng pGL3-BACH-1, 10 ng pEF-RL (Renilla), and 800 ng of either pcDNA3.1 or KSHV miRNA expression vector. Light units are normalized to Renilla luciferase. For B to D, average data from three independent transfections are shown.
|
We next asked whether the additional KSHV miRNAs could target the BACH-1 3'UTR. Cotransfection of pGL3-BACH-1 with individual KSHV miRNA expression vectors (described in reference 57) revealed that of the nine different KSHV miRNAs tested, only miR-K12-11 targeted the BACH-1 3'UTR (Fig. 3D). To further confirm miRNA specificity, derepression assays utilizing miRNA antagomirs were performed. 2'OMe oligoribonucleotides antisense to miRNAs are potent inhibitors of RNA-induced silencing complex-associated miRNA activity (25, 45). Cotransfection of miR-155 or miR-K12-11 expression vectors with the BACH-1 reporter, in addition to an miRNA-specific antagomir, released the inhibitory miRNA effect on the reporter (data not shown). Together, these data validate the observed expression differences from the BACH-1 reporter as being miR155/miR-K12-11 dependent.
Seed match site 2 is essential for combinatorial regulation of BACH-1 by miR-155 and miR-K12-11. To confirm that the observed BACH-1 inhibition depends on the presence of miR-155 and miR-K12-11 seed match sites, a mutational analysis of the BACH-1 3'UTR was performed. First, 3'UTR truncation mutants were generated by deleting one or more of the seed match sites (Fig. 4A). Deletion of all four sites or sites 2, 3, and 4 abrogated the miRNA-dependent inhibition, while deletion of site 1 only had no effect on luciferase expression (Fig. 4B), demonstrating that seed match sites 2, 3, and 4 within the BACH-1 3'UTR are likely targeted.
![]() View larger version (33K): [in a new window] |
FIG. 4. Seed match site 2 is essential for combinatorial regulation of BACH-1 by miR-155 and miR-K12-11. (A) Schematic of the BACH-1 3'UTR with potential miRNA seed match binding sites (sites 1 to 4) and restriction enzyme sites. "A" to "D" are truncation mutants generated by restriction enzyme digestion. (B) Mapping of miRNA binding sites. BACH-1 fragments A to D were inserted into pGL3-Promoter, and 40 ng of each reporter was cotransfected with 10 ng pEF-RL and 800 ng either pcDNA3.1, pmiR-155, or pmiR-K12-11 into 293 cells. (C) Site-directed mutagenesis to validate mapping data. Individual or double seed match sites within pGL3-BACH-1 were mutated to an XhoI site. Forty nanograms of each reporter was cotransfected with 10 ng pEF-RL and 800 ng miRNA expression vectors as described above (B). For B and C, lysates were harvested at 72 h, and relative light units (RLU) were normalized to Renilla luciferase. Shown are average data for five to eight independent transfections. (D) Potential miR-155 and miR-K12-11 binding to site 2 within the BACH-1 3'UTR (RNA-hybrid) (53). The binding position within the 3'UTR and the minimum free hybridization energy (mfe) for each miRNA are shown.
|
Site 2 represents an 8-mer site for both miRNAs containing the seed match (nt 2 to 7) plus a match at position 8 and an A at position 1 (Fig. 4D) (38). Very recently, Grimson et al. reported a set of criteria that further contribute to the affinity by which an miRNA binds to target mRNAs. We note that the binding of miR-155/miR-K12-11 to site 2 within the BACH-1 3'UTR fulfills several of these criteria, which include an AU-rich environment surrounding the miRNA binding site (>70%), a seed sequence preceded and followed by an A (Fig. 4D), and the overall positioning within the periphery of the 3'UTR (Fig. 3A) (20).
Analysis of BACH-1 expression in PEL and BL cell lines. Following the identification of BACH-1 as a potential target using bioinformatic tools and the confirmation of its regulation by luciferase reporter assays, we asked whether both miRNAs can regulate endogenous BACH-1 protein levels in BJAB cells, an EBV- and KSHV-negative BL-cell line which does not express either miRNA (Fig. 1D). Western blot analysis revealed that the transient transfection of either miR-155 or miR-K12-11 expression vectors leads to 70% and 80% decreases in BACH-1 protein levels, respectively, compared to a control vector, further confirming BACH-1 as a target (Fig. 5A and B).
![]() View larger version (12K): [in a new window] |
FIG. 5. BACH-1 is down-regulated in miR-155- and miR-K12-11-expressing cells. (A) Transient expression of miR-155 or miR-K12-11 affects BACH-1 protein levels. BJAB cells were transfected with 2 µg miRNA expression vector or pcDNA3.1 control (nucleofection) (Amaxa, Inc.). Lysates were harvested at 72 h and probed for endogenous BACH-1 or actin. (B) Graph of BACH-1 levels normalized to actin (analyzed by NIH ImageJ) for transfected BJAB cells shown in A. (C) PEL cells expressing miR-K12-11 and BL cells expressing miR-155 express low amounts of BACH-1 protein. Lysates from 1 x 105 BJAB, RAJI, BCBL-1, BC-1, or VG-1 cells were probed for endogenous BACH-1 and actin.
|
Gene expression profiling reveals a common set of down-regulated genes in response to miRNA expression. To investigate the effect of miR-155/miR-K12-11 expression on the cellular transcriptome, we generated 293 cells stably expressing miR-155, miR-K12-11, or the vector control (pcDNA3.1) and performed Affymetrix-based gene expression profiling. miRNA expression in 293/pmiR-155 and 293/pmiR-K12-11 cells was confirmed both by luciferase derepression assays using specific miRNA sensor vectors and antagomirs (Fig. 6A) and Northern blot analysis (Fig. 6B and C). To avoid the risk of analyzing integration events, we examined cell pools rather than individual cell clones. RNA was harvested from two independent cultures of each cell line at two time points. A total of 12 samples were used to synthesize cRNA and hybridized to Affymetrix U133 2.0 Plus Human GeneChips containing over 47,000 human probe sets. RNA extraction, quantification, cRNA synthesis, hybridization, and washing steps were performed as recommended by the manufacturer and as previ ously described (57). Figure 7A shows a heat map representing the expression profiles of genes significantly altered in the presence of miR-K12-11 or miR-155. In both 293/pmiR-155 and 293/pmiR-K12-11 cells, a total of 195 annotated genes, represented by 208 probe sets, were altered compared to the vector control. Of those genes altered, we found 137 genes (70%) to be down-regulated. A cross-validation analysis showed that of these 137 genes, 78 exhibited a cross-validation of 75% or greater. This analysis, based on t values of signal intensities, gives statistically robust data while including probe sets with relatively small changes (n-fold). A common set of 66 genes was significantly down-regulated between –1.11- and –11.05-fold in response to both miR-155 and miR-K12-11 (see Table S1 in the supplemental material).
![]() View larger version (35K): [in a new window] |
FIG. 6. miRNA expression in stable cell lines. (A) Luciferase derepression assays were performed using 293/pmiR-K12-11, 293/pmiR-155, or 293/pcDNA3.1 cells using 40 ng of the indicated sensor vector and increasing amounts (20 to 80 pmol) of a 2'OMe antagomir specific to the miRNA of interest. 2'OMe K12-10 and 100 ng of pcDNA3.1 were used as filler. Light units are normalized to total protein by BCA (Pierce). (B and C) Northern blot analysis of miRNA expression in stable cell lines. Thirty micrograms of total RNA was loaded per lane and hybridized to a probe for either miR-K12-11 (B) or miR-155 (C). RNA from BCBL-1 and RAJI cells was used as a positive control.
|
![]() View larger version (34K): [in a new window] |
FIG. 7. Gene expression profiling reveals a common set of target genes. (A) Genes are down-regulated in response to miRNA expression. Colors represent changes in variance-normalized gene expression differences for individual genes represented by the probe sets as indicated on the color scale. Four samples were used for each cell line. One hundred ninety-five genes were changed in response to miRNA expression: 137 genes were down-regulated, and 48 were up-regulated. Sixty-four genes showed >75% cross-validation. (B) The 3'UTRs of down-regulated genes are enriched for seed match sites. 3'UTRs of altered genes were scanned for seed match sites. Shown is the percentage of 3'UTRs with seed match sites for either the total number of altered genes (total), up-regulated (up) or down-regulated (down) genes, or the 66 commonly down-regulated genes (66 common). Thirty-three out of 137 down-regulated genes contained 6-mer seed match sites within their 3'UTRs, while 20 out of the 66 commonly down-regulated genes contained seed match sites.
|
Following the identification of these new potential miRNA targets, we examined pathways affected by miRNA expression using Biocarta. The majority of pathways included cell signaling, cell division, and T-cell activation. Additionally, we identified several genes with known roles in apoptosis, such as LDOC1, PHF17, BCL6B, BCLAF1, and BCL2L11 (27, 46). The significance and potential roles of these genes in KSHV biology will be discussed below.
Luciferase reporter assays confirm miRNA targeting of LDOC1, MATR3, and TM6SF1. Based on the bioinformatic data, we chose five genes for further examination and cloned their 3'UTRs into pGL3-Promoter downstream of the luciferase gene. Reporter constructs were cotransfected into 293 cells with pmiR-155, pmiR-K12-11, or pcDNA3.1 as a control. We observed a significant inhibition of luciferase from reporters containing the 3'UTRs of LDOC1, MATR3, and TM6SF1 in the presence of miR-155 or miR-K12-11 (Fig. 8). In contrast, neither miRNA had any effect on the PHF17 or NFAT2CIP 3'UTRs (Fig. 8). Thus, these assays confirmed LDOC1, MATR3, and TM6SF1 as being targets of both miR-155 and miR-K12-11.
![]() View larger version (19K): [in a new window] |
FIG. 8. Validation of candidate target genes. The 3'UTRs of LDOC1, MATR3, TM6SF1, PHF17, and NFAT2CIP contained at least one seed match site and were cloned into pGL3-Promoter downstream of luciferase. 3'UTR reporters were cotransfected into 293 cells with 10 ng pEF-RL and either 800 ng pcDNA3.1 or the indicated amounts of miRNA expression vector to determine miRNA targeting. Four hundred nanograms of pcDNA3.1 was used as filler for the transfection of 400 ng miRNA expression vector. Lysates were harvested at 72 h, and relative light units (RLU) were normalized to Renilla luciferase. For pGL3-LDOC1, pGL3-MATR3, and pGL3-TM6SF1, a 1.5- to 2.5-fold decrease in luciferase was observed with miR-155 and miR-K12-11 expression. Shown are average data from three to six independent transfections.
|
|
|
|---|
Interestingly, this is not an anomaly, since alignments revealed additional herpesvirus-encoded miRNAs that share seed sequence homology with human miRNAs (Fig. 1A). KSHV miR-K12-6-5p has homology to miR-15a and miR-16, both associated with tumor suppressor activity (13), while EBV BART5 has homology to miR-18, a member of the oncogenic miR-17-92 cluster (24). This suggests that a subset of herpesvirus miRNAs represent cellular homologues that have arisen either by capturing host miRNA genes or through coevolution with their hosts. The fact that the pre-miRNAs of miR-155 and miR-K12-11 are vastly different would suggest coevolution similar to precursors of some let-7 miRNA family members, who share a common seed sequence and can target similar 3'UTRs (1, 31) (www.microrna.sanger.ac.uk).
By utilizing reporter assays and detailed mutational analyses, we found that both miR-155 and miR-K12-11 target seed match sites within the BACH-1 3'UTR (Fig. 3 and 4). In contrast to these reporter assays, the RNA levels of BACH-1 were not significantly affected in miRNA-expressing 293 cells as determined by expression profiling. Importantly, BACH-1 protein levels were down-regulated in BJAB cells transiently transfected with miRNA expression vectors (Fig. 5A). Furthermore, we observed low levels of BACH-1 in PEL or BL lines expressing both miRNAs compared to BJAB cells (Fig. 5C). Additionally, BACH-1 mRNA levels are down-regulated in latently KSHV-infected endothelial cells (4). BACH-1 is a broadly expressed transcriptional repressor, which regulates genes involved in the hypoxia response, such as heme-oxygenase 1 (HMOX1) (26a). Increased levels of HMOX1, which enhance cell survival and proliferation, have been reported following de novo KSHV infection of endothelial cells (44). While the increase in HMOX1 has been attributed in part to the KSHV G-protein-coupled receptor (41), the observed down-regulation of BACH-1 may be partially due to miR-K12-11 expression during latency.
The seed sequence homology of these two naturally occurring miRNA orthologs also provides an opportunity to analyze the contribution of sequences outside the seed match to mRNA targeting and target site selection. A detailed analysis of the BACH-1 3'UTR demonstrated that both miR-155 and miR-K12-11 preferentially utilized seed match site 2 for targeting; however, sites 3 and 4 also contributed (Fig. 4B and C). Further analysis of the 66 down-regulated genes using miRanda showed that >30% of 3'UTRs are predicted to contain miR-155 or miR-K12-11 binding sites, a majority of which (>50%) contain exact seed match sites (not shown) (17). Very recently, the Bartel laboratory published a new model outlining additional determinants of miRNA target site selection for more accurate target prediction (20). Analysis of the four seed match sites within the BACH-1 3'UTR revealed that target site 2, which conveyed miRNA-dependent inhibition (Fig. 4C), fulfilled several of those new criteria. Hence, a detailed analysis of additional targets will also aid a better understanding of how sequences outside the seed sequence contribute to miRNA targeting.
The biological consequences of these newly identified miRNA targets are currently under investigation. Using Biocarta, we found that the majority of pathways affected by genes that changed in response to miRNA expression included cell signaling, cell division, apoptosis, and T-cell activation. Additionally, we identified several miR-155- and miR-K12-11-responsive genes known to be aberrantly expressed in human malignancies. These include LDOC1, a regulator of apoptosis (27, 46), as well as several Bcl-2 and Bcl-6 family members (BCL2L11, BCL6B, and BCLAF1).
BCL2L11, known as Bim, inactivates the Bcl-2-like proteins, thereby inducing apoptosis (2, 14). BCLAF1 (Bcl-2 associated transcription factor 1 [Btf]) has been shown to interact with the antiapoptotic Bcl-2 and Bcl-xL family members and, when overexpressed, to induce apoptosis (32). BCL6B, a transcriptional repressor also known as BAZF, is an NF-
B regulator and homolog of the Bcl-6 oncogene; however, unlike Bcl-6, which is ubiquitously expressed, BAZF is expressed predominantly in the heart, lung, and, importantly, stimulated lymphocytes (49). Both Bim and BAZF have roles in cellular immunity. Bim is required for the apoptosis of activated T cells following viral infection, which is critical for limiting antiviral immune responses (50). BAZF is required for secondary responses of memory CD8+ T cells (40) and was recently linked to the homeostasis of hematopoietic progenitors (7). Thus, the observed down-regulation of these genes by miR-155 or miR-K12-11 suggests important roles in the regulation of immune responses to infection as well as in targeting apoptosis to regulate cell survival during lymphocyte differentiation. It is interesting that, to date, of the only
50 experimentally determined human miRNA targets, a significant proportion are regulators of apoptosis, such as targets of miR-15a and miR-16 (13). Furthermore, KSHV miRNAs have recently been reported to target apoptosis regulators (22, 57).
While there are only a few experimentally confirmed targets for miR-155, it was the first human miRNA suggested to be oncogenic (61) and was recently shown to be important for lymphocyte differentiation and immunity. miR-155 is expressed in activated macrophages following treatment with Toll-like receptor ligands, suggesting that miR-155 has a role in innate immunity (48). Recently, it was shown that BIC/miR-155 plays a significant role in B- and T-lymphocyte maturation (54, 63). BIC is absent in resting and progenitor B cells but induced upon activation; miR-155 is expressed shortly thereafter and regulates the germinal center reaction during B-cell maturation (63). The late stages of B-cell development, particularly the affinity maturation in conjunction with helper T-cell signaling, are critically dependent upon a timely miR-155 expression whereby both the induction and down-regulation of miR-155 are important (63).
Intriguingly, KSHV-associated PELs are of B-cell lineage and have a distinct developmental phenotype. PELs have rearranged immunoglobulin genes as well as somatic point mutations within rearranged immunoglobulin variable genes (18, 43), indicating that B-cell activation during the germinal center reaction has occurred. Accordingly, the PEL developmental phenotype was classified as plasmablastic (33), and expression profiling studies suggest that PEL represent B cells blocked late in their differentiation pathway towards an antibody-secreting plasma cell (29). We showed that PEL cells do not express miR-155 but do express high levels of miR-K12-11 (Fig. 1D). Exactly some of these late steps seem to be dependent upon a short miR-155 expression burst during normal B-cell development (63).
Based on our data demonstrating common targets for both miRNAs, we propose a model in which sustained miR-K12-11 expression continues to regulate a set of miR-155 targets at a stage when miR-155 is not expressed and thereby directly contributes to dysregulated nonmature B-cell proliferation and, potentially, lymphomagenesis (Fig. 9). Clearly, this hypothesis will have to be tested in appropriate in vivo models, such as NOD/SCID mice, that allow the study of human B-cell maturation.
![]() View larger version (8K): [in a new window] |
FIG. 9. Model of a potential role for miR-K12-11 in lymphomagenesis. Germinal-center-dependent B-cell maturation is dependent on a precise miR-155 expression burst whereby both induction and shutoff are equally important (54, 63). In contrast, after naïve B cells are infected with KSHV, sustained miR-K12-11 expression could contribute to proliferation in the absence of terminal differentiation into a mature antibody-producing B cell, a phenotype which is congruent with that characterized for PEL cells (29).
|
Finally, we found a total of four viral miRNAs (Fig. 1A) harboring seed match sequence identity to cellular miRNAs. Most importantly, three of these are potential orthologs of human miRNAs known to have either oncogenic or tumor suppressor activity. Hence, these observations strongly suggest that DNA tumor viruses have utilized molecular mimicry of miRNA genes to regulate host cellular environments in the same fashion as the long list of homologous cytokines, growth factors, and immunomodulatory genes commonly found in DNA tumor virus genomes (47, 55).
Published ahead of print on 19 September 2007. ![]()
Supplemental material for this article may be found at http://jvi.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»