|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Division of Microbiology, Institute for Animal Health, Compton RG20 7NN, UK
2 LRF Molecular Haematology Unit, Nuffield Department of Clinical Laboratory Sciences, University of Oxford, Oxford OX3 9DU, UK
3 Bacterial Pathogenesis and Functional Genomics Group, Sir William Dunn School of Pathology, University of Oxford, Oxford OX1 3RE, UK
Correspondence
Venugopal Nair
venu.gopal{at}bbsrc.ac.uk
| ABSTRACT |
|---|
|
|
|---|
Supplementary figures and tables are available with the online version of this paper.
| INTRODUCTION |
|---|
|
|
|---|
Marek's disease virus (MDV) is a highly contagious, oncogenic alphaherpesvirus of the genus Mardivirus; it is associated with Marek's disease (MD), a naturally occurring, rapid-onset T-cell lymphoma of chicken (Calnek, 1986
). The MDV genome encodes several miRNAs that map to the MDV-encoded oncogene Meq and the LAT (latency-associated transcript) regions of the virus (Burnside et al., 2006
; Burnside & Morgan, 2007
). Although the target genes regulated by MDV-encoded miRNAs are yet to be discovered, high levels of their expression in MDV-transformed cell lines and tumours suggest that they play important roles in oncogenesis (Morgan et al., 2008
; Xu et al., 2008
). From a small RNA library generated from MDV-transformed lymphoblastoid cell line MSB-1 (Akiyama & Kato, 1974
), we have previously demonstrated high levels of expression of MDV-encoded miRNAs (Yao et al., 2007
, 2008
). Elevated levels of expression of some of these miRNAs were also confirmed by real-time quantitative PCR in these cells (Xu et al., 2008
). Several host miRNAs that are associated directly with oncogenicity, such as miR-17-92, miR-21 and let-7i, were also present at a high frequency in the MSB-1 library, suggesting a role for these miRNAs in neoplastic transformation of these cells (Yao et al., 2008
). Although analysis of miRNAs by cloning (Yao et al., 2008
) or high-throughput sequencing (Burnside et al., 2008
) is used to identify upregulated miRNAs in tumours or transformed cells, such studies do not provide differential expression profiles of miRNAs in different cell types. Comparisons of miRNA expression profiles between neoplastically transformed and normal cells using miRNA microarrays have enabled the identification of specific miRNA expression signatures in different types of cancer cell, some of which have been shown to be useful indicators of cell type, stage of differentiation and even prognosis of the cancer (Calin & Croce, 2006
; Rosenfeld et al., 2008
). Only a limited number of studies comparing the global miRNA expression profiles of virally transformed tumour cells have been carried out, particularly in tumour cells transformed by oncogenic herpesviruses that themselves encode multiple miRNAs (Sullivan & Grundhoff, 2007
). Here we describe the results of the comparison of the miRNA expression profile of seven different MDV-transformed T-lymphoblastoid cell lines with that of normal chicken splenocyte or CD4+ T-cell populations.
| METHODS |
|---|
|
|
|---|
Chicken splenocytes, CD4+ T cells and magnetic cell sorting.
Single-cell suspensions of lymphocytes were prepared from spleen tissues of uninfected birds by using Histopaque-1083 (Sigma-Aldrich) density-gradient centrifugation. CD4+ T cells were isolated by magnetic cell sorting using mouse anti-chicken CD4 antibodies (Chan et al., 1988
) and goat anti-mouse IgG microbeads (Miltenyi Biotec). After each antibody treatment, cells were washed three times with PBS containing 0.5 % bovine serum albumin. At each wash, the cell suspension was centrifuged at 450 g for 10 min. Positively stained cells were sorted through an AutoMACS Pro Separator (Miltenyi Biotec). Purity of the sorted cells was confirmed to be >99 % by flow cytometry after labelling with monoclonal anti-goat/sheep IgG–fluorescein isothiocyanate (Sigma) antibody (data not shown).
Microarray analysis of miRNA expression.
Preparation of probes and hybridization to the arrays were carried out by using methods described previously (Lawrie et al., 2007
, 2008
). Briefly, 500 ng purified miRNA from lymphoblastoid cell lines, normal splenocytes or CD4+ T-cell populations was labelled with either Cy3 or Cy5 dye using an Array 900microRNA RT kit from Genisphere and hybridized to the µRNA microarray described previously (Lawrie et al., 2008
). The array contains miRNA probe sets (designed from miRBase v. 9.2) spotted in quadruplicate non-adjacently (the sequences of the probes are available at http://www.microRNAworld.com). Where the homologous sequences of miRNAs in different species are identical, miRNA sequences from only one species was spotted on the array. In addition, probes for the mature MDV-1-encoded miRNAs miR-M2-3p, miR-M2-5p, miR-M3-3p, miR-M3-5p, miR-M4-3p, miR-M4-5p, miR-M5-3p, miR-M11-5p and miR-M12-3p (Yao et al., 2008
) were also included on the array. A model design with splenocytes and/or CD4+ T cells as reference was used and the expression values are depicted as log ratios of test and reference samples. Image analysis was done by using BlueFuse software (BlueGnome).
Statistical analysis of microarray data.
Data were normalized within each microarray by subtracting the mean log2 ratio from each measurement. Quantile normalization was then used to standardize the data across arrays, and a linear model was fitted to each miRNA using Limma (Smyth, 2005
). The resultant P-values were adjusted for multiple testing by using the Benjamini–Hochberg correction of the false-discovery rate (Benjamini & Hochberg, 1995
). The P-values and mean expression were calculated for each cell type. Those miRNAs showing consistent differential expression across all cell types were subjected to a second analysis where all samples were treated as virus-infected, regardless of cell type. A similar analysis in Limma was performed, resulting in P-values calculated for expression across all cell types. The data were sorted according to the log2 ratio and heat maps were produced by using R (http://www.R-project.org).
Inducible expression of mature miRNAs.
Constructs for inducible expression of miRNAs were generated in the pRTS-1-SVP-Tom(–) vector, constructed from the pRTS-1 plasmid (Bornkamm et al., 2005
) by replacing the hygromycin B-resistance gene with a puromycin-resistance gene. The inducible bidirectional promoter in this construct expresses monomeric red fluorescent protein td-tomato (Shaner et al., 2004
) and the miRNA of interest simultaneously in a tightly regulated, doxycycline (Dox)-inducible system. The inserted sequence of each miRNA consists of the stem–loop structure and 100–200 bp of upstream and downstream flanking genome sequences, a feature designed to ensure that the expressed miRNAs are processed as naturally as possible. The approximately 500 bp fragment of each miRNA was obtained by RT-PCR using RNA extracted from chicken embryo fibroblasts by using the following primers: gga-miR-155 (5'-AGATCTCTGATGTCTGTACTCTTTATGAC-3' and 5'-CTCGAGCCCAGTGCCCTTAACTTAG-3'); gga-miR-223 (5'-AGATCTGCAACGTCTGTCCTGTCC-3' and 5'-CTCGAGCAGGAACTGTACCAGCAG-3'). As the sequence of the chicken orthologue was not available, we used hsa-miR-150 for generating the expression constructs of miR-150, using RNA prepared from HEK293T cells with primers 5'-AGATCTCTTCCTGCCCTCTTTGATG-3' and 5'-CTCGAGCAATAGAAACAGGTGTACTTTG-3'. More recently, the mature gga-mir-150 sequence was shown to be identical to the hsa-miR-150 sequence except for a single nucleotide substitution (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL6541). Agarose gel-purified PCR fragments were cloned into the BglII and XhoI sites of the pRTS-1-SVP-Tom(–) vector and confirmed by sequence analysis. For the expression of cloned miRNAs, HEK293T cells were transfected with the pRTS-1–miRNA constructs by using Lipofectamine 2000 (Invitrogen) and stable cells were selected with puromycin (1 µg ml–1). Respective miRNAs induced with the addition of Dox (1 µg ml–1) were detected by Northern blotting analysis.
Luciferase reporter assays.
In order to demonstrate that the Dox-induced miRNAs in stably transfected HEK293T cells are functional, we carried out reporter assays on one of the stable cell lines expressing gga-miR-155 by analysing its effect on the previously characterized target Pu.1. For this, we generated a reporter construct in which a 110 bp fragment of the chicken Pu.1 3' untranslated region (UTR) transcript (GenBank accession no. NM_205023
[GenBank]
) that contained the predicted miR-155-response element (MRE) is inserted downstream of Renilla luciferase in the psiCHECK-2 vector (Promega) to generate the Pu.1-3'UTR-wt construct. A mutant construct (Pu.1-3'UTR-mt) with mutations in the MRE sequences was also constructed. The details of the oligonucleotides used for generating the chicken Pu.1-3'UTR reporter constructs are shown elsewhere (Zhao et al., 2009
). HEK293T cells stably expressing gga-miR-223 were used as a control. Transfection into HEK293T cells for luciferase reporter assays was carried out in 96-well plates with Lipofectamine 2000 (Invitrogen). Firefly and Renilla luciferase activities were measured consecutively with the Dual-Luciferase Reporter Assay system (Promega), using a Lucy-1 luminometer (Anthos Labtec). In all cases, the constitutively expressed firefly luciferase activity in the psiCHECK-2 vector served as a normalization control for transfection efficiency.
Preparation of RNA for Northern blotting.
Preparation of RNA for Northern blotting analysis from splenocytes, CD4+ and CD4– populations of T lymphocytes, MDV- and retrovirus-transformed cell lines, and HEK293T cells stably selected for recombinant pRTS-1 vectors expressing miRNAs was carried out by using standard methods as described previously (Yao et al., 2007
, 2008
). Antisense oligonucleotides to the miRNAs gga-miR-155, gga-miR-223 and hsa-miR-150 were used as probes.
| RESULTS |
|---|
|
|
|---|
|
Changes in host miRNA profiles in MD tumour cell lines
Examination of the global miRNA expression profiles of the seven MDV-1-transformed cell lines revealed changes in several host miRNAs (Fig. 1
). Major differences in the host miRNA expression profiles could also be observed between MDV-1-transformed cell lines and the MDV-negative cell line AVOL-1, demonstrating the differences in the molecular oncogenic pathways between these cell lines. Microarray readouts from our studies did demonstrate differences between cell lines with regard to the expression of individual miRNAs. Although such differences could be important for individual cell lines, our main interest was to look for miRNA profiles that are conserved in all MDV cell lines, as this could give insights into the fundamental molecular pathways of miRNA-mediated gene regulation in MDV transformation. Our results showed that several host-encoded miRNAs, such as miR-155, miR-223, miR-150, miR-451, miR-26a and miR-126, were downregulated in all MDV-transformed cell lines relative to the levels in normal splenocytes or CD4+ T cells (Fig. 1
). As miR-223 and miR-150 were also downregulated in retrovirus-transformed AVOL-1 cells, the reduced expression of these two miRNAs is thought to be a broader feature of transformed T cells, irrespective of the viral aetiology. However, this was not the case with miR-155, the levels of which were consistently reduced in all of the seven MDV-1-transformed lymphoblastoid cell lines whilst the expression levels in the AVOL-1 cell line were very high (Fig. 1a
), demonstrating that the downregulation of miR-155 is a feature unique to MDV transformation of T cells.
As all of the MDV-1-transformed cell lines used in this study have a CD4+ T-cell phenotype, we also wondered whether it is appropriate to use whole splenocyte populations as the reference sample in the analysis. In order to rule out the possibility that the altered expression profiles of miRNAs in the MDV-transformed cell lines are not due to the use of splenocytes as the reference, we also repeated the analysis of miRNA expression in four MDV-1-transformed cell lines with purified CD4+ T cells as the reference. Overall, the results were largely consistent with the values obtained by using normal splenocytes as the reference (Fig. 1b
). However, the use of purified CD4+ cells as the reference resulted in the demonstration of increased expression of several more host miRNAs, such as miR-146b, miR-454, miR-7b, miR-34a, miR-18a, miR-133a, miR-29a and gga-1et-7d (Fig. 1b
). As the splenocyte reference data represent the miRNA expression of multiple lymphocyte populations, it was not surprising to see changes in the miRNA expression patterns when purified CD4+ T cells were used as the reference. Data on the log2 fold changes and P-values for each cell line using splenocytes or purified CD4+ T cells as the reference can be seen in Supplementary Tables S1 and S2, respectively.
Analysis of miRNA expression by Northern blotting
For further validation of the microarray data demonstrating the reduced expression of miR-155, miR-223 and miR-150 in MDV-transformed tumour cell lines, we carried out Northern blotting analysis comparing the levels of these miRNAs in four MDV-transformed cell lines with those in normal lymphocyte populations and retrovirus-transformed cell lines AVOL-1 and HP45. Northern blotting analysis confirmed the observations of reduced expression of these three miRNAs in the microarray readouts of MDV-transformed cells. The levels of gga-miR-155 signal detected were very low in the normal splenocyte and CD4+/CD4– T-cell populations (Fig. 2a
). These low levels were reduced further in all MDV-transformed cell lines included in the assay. In contrast, stronger gga-miR-155 signals were evident in avian retrovirus-transformed cell lines HP45 and AVOL-1. The levels of gga-miR-223 were lower in normal CD4+ T cells than in whole spleen-cell or CD4– T-cell populations. However, no hybridization signals for gga-miR-223 expression were evident in any of the cell lines transformed by either MDV-1 or avian retroviruses, suggesting that the downregulation of gga-miR-223 is a broader feature of lymphocyte transformation, irrespective of the viral aetiology. The expression of miR-150 also appeared to be restricted to the untransformed cells, as there was no evidence of its expression in transformed cells. Although miR-150 and miR-223 expression appeared to be similar in this respect, the levels of miR-150 were much higher in CD4+ T cells (Fig. 2a
). We also evaluated the Dox-inducible expression of miR-155, miR-223 and miR-150 from the pRTS-1 vector (Bornkamm et al., 2005
). Northern blotting analysis of HEK293T cells expressing the pRTS-1–miRNA constructs showed high levels of expression of each of the three mature miRNAs, regulated tightly in a Dox-inducible manner (Fig. 2b
).
|
|
| DISCUSSION |
|---|
|
|
|---|
We and others have documented recently that the highly oncogenic MDV-1 encodes several novel miRNAs (Burnside & Morgan, 2007
; Burnside et al., 2008
; Morgan et al., 2008
; Yao et al., 2008
). High levels of expression of these miRNAs in lymphomas and transformed cell lines have been demonstrated by direct cloning, Northern blotting and quantitative RT-PCR (Burnside et al., 2006
; Xu et al., 2008
; Yao et al., 2008
). Although these studies have been valuable in identifying miRNAs that are expressed at high levels in these cells, they do not always provide comprehensive miRNA expression profiles, particularly of those miRNAs that are downregulated in the transformed cells. With a view to examining the global expression of miRNAs in MDV-transformed cells, we carried out miRNA expression profiling of seven independent MDV-transformed tumour cell lines by using miRNA microarray analysis. Each of these cell lines was derived from an independent MD lymphoma. As demonstrated previously with MSB-1 cells (Yao et al., 2007
, 2008
), MDV-1-encoded miRNAs were indeed the most abundant miRNAs in all of the cell lines. High-level expression of virus-encoded miRNAs appears to be a feature common to virus-transformed cell lines, as cells transformed by other oncogenic herpesviruses, such as Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein–Barr virus, also showed high levels of expression of virus-encoded miRNAs (Cai et al., 2005
; Lawrie et al., 2008
; Pratt et al., 2009
). Higher copy numbers of viral genomes and active transcription of miRNA genes may account for the higher expression of virus-encoded miRNAs in the virus-transformed cell lines, although one cannot rule out the possibility of differential processing of virus-encoded miRNAs in these cells. The functions and putative targets of most of the MDV-encoded miRNAs remain unknown. However, we have shown recently that MDV-miR-M4, one of the most abundantly expressed virus-encoded miRNAs in all of the cell lines, is a functional orthologue of miR-155 with the potential to target important lymphocyte-specific transcription factors such as Pu.1 (Zhao et al., 2009
). More efforts in the future to identify the potential targets of other MDV-encoded miRNAs (Morgan et al., 2008
) will unravel more molecular pathways of oncogenesis.
Microarray data analysis of the changes in the global expression profiles of host-encoded miRNAs in MDV-transformed cell lines could be grouped into (i) those that are restricted only to some of the MDV-transformed cell lines and (ii) those that appear to be conserved across all of the cell lines. The changes in miRNA expression in individual cell lines, such as the increased expression of miR-221/miR-222 in MSB-1 cells (Lambeth et al., 2009
), are likely to be important in the regulation of respective target proteins in individual cell lines. However, in this paper, we focus on the global changes in miRNA expression common to all MDV-transformed cell lines.
The miRNA profile of the seven MD tumour cell lines showed changes in the expression of several miRNAs. These included the downregulation of miR-150, miR-223 and miR-155, confirmed by Northern blotting analyses (Figs 1
and 2a
). Of these, miR-150 and miR-223 were also downregulated in AVOL-1 cells, demonstrating it to be a broader feature of lymphocyte transformation, regardless of the viral aetiology. Reduced expression of miR-150 has also been reported in human lymphoid malignancies such as diffuse large B cell lymphoma (Garzon & Croce, 2008
; Landgraf et al., 2007
; Lawrie et al., 2008
), indicating its conserved function across different species. Increasing evidence suggest that miR-150 functions through the transcription factor c-myb (Garcia & Frampton, 2008
; Lin et al., 2008
; Lu et al., 2008
; Xiao et al., 2007
; Zhou et al., 2007
), and the dysregulation of c-myb and its targets could be important in T-cell transformation. In the case of miR-223, although all of the regulatory mechanisms are not yet understood fully, recent studies have indicated a clear role for miR-223 in haematopoiesis, as well as in malignancies (Baek et al., 2008
; Garzon & Croce, 2008
; Johnnidis et al., 2008
; Merkerova et al., 2008
). Whilst identification of the potential targets is important to understand fully the molecular pathways, involvement of miR-223 appears to be logical in MDV-induced lymphocyte transformation.
The downregulation of miR-155 observed by microarray analysis was unique to MDV-transformed cell lines, as it was upregulated in the MDV-negative AVOL-1 cell line (Fig. 1
). Although the levels of miR-155 in the normal lymphocyte populations were low by Northern blotting analysis, it was clear that MDV-transformed cell lines showed a distinct reduction in hybridization signals, especially when compared with the retrovirus-transformed lymphocyte cell lines HP45 and AVOL-1 (Fig. 2a
). Several recent studies have highlighted the potential multiple roles of miR-155 in functions ranging from innate immune responses to oncogenicity (Garzon & Croce, 2008
). The molecular mechanisms that drive down the expression of miR-155 in MDV-transformed cell lines are not known. However, some of its functions on targets such as Pu.1 could be rescued by MDV-miR-M4, a highly expressed MDV-1-encoded functional orthologue of miR-155 (Zhao et al., 2009
). Although the regulatory expression dynamics of miR-155 and MDV-miR-M4 are not understood fully, the existence of autoregulatory mechanisms of miR-155 expression mediated through a common set of targets cannot be ruled out. It is interesting that, in KSHV-infected primary effusion lymphoma cell lines, miR-155 was also found to be downregulated in favour of the KSHV-encoded miR-K12-11 homologue (Skalsky et al., 2007
).
The data from this study have enabled us to characterize the miRNome of MDV-transformed tumours. Although this has provided valuable insights into the expression profiles of miRNAs in these cell lines, the major challenge will be in the identification of the putative targets of the differentially expressed miRNAs in these cells. Although bioinformatic predictions of miRNA targets are valuable, the development of systems for functional characterization of miRNA targets is important to understand the pathways of oncogenesis. The tightly regulated, Dox-inducible miRNA expression system of the differentially expressed miRNAs that we developed in HEK293T cells will be valuable in identifying the putative functional targets of these miRNAs. Demonstration of the expression of mature miRNAs in a Dox-dependent manner clearly showed the proper processing of these miRNAs in this system. For functional validation of the system, we analysed the putative targeting of miR-155 on one of the validated target proteins, Pu.1 (Zhao et al., 2009
). The tightly regulated expression of miR-155 and the specific silencing of the relative luciferase levels with reporter assays with wild-type 3' UTR reporter constructs (Fig. 3
) provide a platform for functional analysis of the putative targets of differentially expressed miRNAs.
In summary, the data presented here demonstrate that miRNA expression profiling using microarrays is a powerful approach for analysing the relative levels of several miRNAs simultaneously. This study, the first of its kind in MDV-transformed cell lines, demonstrates that, in addition to the overexpression of several MDV-encoded miRNAs, downregulation of some of the host-encoded miRNAs is also a hallmark of MDV transformation. Determination of the miRNA profile is a first step towards identification of the regulatory networks of gene expression in these cell types.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Baek, D., Villen, J., Shin, C., Camargo, F. D., Gygi, S. P. & Bartel, D. P. (2008). The impact of microRNAs on protein output. Nature 455, 64–71.[CrossRef][Medline]
Benjamini, Y. & Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol 57, 289–300.
Bornkamm, G. W., Berens, C., Kuklik-Roos, C., Bechet, J. M., Laux, G., Bachl, J., Korndoerfer, M., Schlee, M., Holzel, M. & other authors (2005). Stringent doxycycline-dependent control of gene activities using an episomal one-vector system. Nucleic Acids Res 33, e137
Burnside, J. & Morgan, R. W. (2007). Genomics and Marek's disease virus. Cytogenet Genome Res 117, 376–387.[CrossRef][Medline]
Burnside, J., Bernberg, E., Anderson, A., Lu, C., Meyers, B. C., Green, P. J., Jain, N., Isaacs, G. & Morgan, R. W. (2006). Marek's disease virus encodes microRNAs that map to meq and the latency-associated transcript. J Virol 80, 8778–8786.
Burnside, J., Ouyang, M., Anderson, A., Bernberg, E., Lu, C., Meyers, B. C., Green, P. J., Markis, M., Isaacs, G. K. & other authors (2008). Deep sequencing of chicken microRNAs. BMC Genomics 9, 185[CrossRef][Medline]
Cai, X., Lu, S., Zhang, Z., Gonzalez, C. M., Damania, B. & Cullen, B. R. (2005). Kaposi's sarcoma-associated herpesvirus expresses an array of viral microRNAs in latently infected cells. Proc Natl Acad Sci U S A 102, 5570–5575.
Calin, G. A. & Croce, C. M. (2006). MicroRNA signatures in human cancers. Nat Rev Cancer 6, 857–866.[CrossRef][Medline]
Calin, G. A. & Croce, C. M. (2007). Investigation of microRNA alterations in leukemias and lymphomas. Methods Enzymol 427, 193–213.[CrossRef][Medline]
Calnek, B. W. (1986). Marek's disease: a model for herpesvirus oncology. Crit Rev Microbiol 12, 293–320.[Medline]
Chan, M. M., Chen, C. L., Ager, L. L. & Cooper, M. D. (1988). Identification of the avian homologues of mammalian CD4 and CD8 antigens. J Immunol 140, 2133–2138.[Abstract]
Cosmopoulos, K., Pegtel, M., Hawkins, J., Moffett, H., Novina, C., Middeldorp, J. & Thorley-Lawson, D. A. (2008). Comprehensive profiling of EBV microRNAs in nasopharyngeal carcinoma. J Virol 83, 2357–2367.[CrossRef][Medline]
Cullen, B. R. (2006). Viruses and microRNAs. Nat Genet 38 (Suppl.), S25–S30.[CrossRef][Medline]
Cullen, B. R. (2009). Viral and cellular messenger RNA targets of viral microRNAs. Nature 457, 421–425.[CrossRef][Medline]
Fabbri, M., Garzon, R., Andreeff, M., Kantarjian, H. M., Garcia-Manero, G. & Calin, G. A. (2008). MicroRNAs and noncoding RNAs in hematological malignancies: molecular, clinical and therapeutic implications. Leukemia 22, 1095–1105.[CrossRef][Medline]
Garcia, P. & Frampton, J. (2008). Hematopoietic lineage commitment: miRNAs add specificity to a widely expressed transcription factor. Dev Cell 14, 815–816.[CrossRef][Medline]
Garzon, R. & Croce, C. M. (2008). MicroRNAs in normal and malignant hematopoiesis. Curr Opin Hematol 15, 352–358.[CrossRef][Medline]
Ghosh, Z., Mallick, B. & Chakrabarti, J. (2008). Cellular versus viral microRNAs in host–virus interaction. Nucleic Acids Res 37, 1035–1048.[CrossRef][Medline]
Gottwein, E. & Cullen, B. R. (2008). Viral and cellular microRNAs as determinants of viral pathogenesis and immunity. Cell Host Microbe 3, 375–387.[CrossRef][Medline]
Griffiths-Jones, S., Saini, H. K., van Dongen, S. & Enright, A. J. (2008). miRBase: tools for microRNA genomics. Nucleic Acids Res 36, D154–D158.
Javier, R. T. & Butel, J. S. (2008). The history of tumor virology. Cancer Res 68, 7693–7706.
Johnnidis, J. B., Harris, M. H., Wheeler, R. T., Stehling-Sun, S., Lam, M. H., Kirak, O., Brummelkamp, T. R., Fleming, M. D. & Camargo, F. D. (2008). Regulation of progenitor cell proliferation and granulocyte function by microRNA-223. Nature 451, 1125–1129.[CrossRef][Medline]
Lambeth, L. S., Yao, Y., Smith, L. P., Zhao, Y. & Nair, V. (2009). MicroRNAs 221 and 222 target p27Kip1 in Marek's disease virus-transformed tumour cell line MSB-1. J Gen Virol 90, 1164–1171.
Landgraf, P., Rusu, M., Sheridan, R., Sewer, A., Iovino, N., Aravin, A., Pfeffer, S., Rice, A., Kamphorst, A. O. & other authors (2007). A mammalian microRNA expression atlas based on small RNA library sequencing. Cell 129, 1401–1414.[CrossRef][Medline]
Lawrie, C. H., Soneji, S., Marafioti, T., Cooper, C. D., Palazzo, S., Paterson, J. C., Cattan, H., Enver, T., Mager, R. & other authors (2007). MicroRNA expression distinguishes between germinal center B cell-like and activated B cell-like subtypes of diffuse large B cell lymphoma. Int J Cancer 121, 1156–1161.[CrossRef][Medline]
Lawrie, C. H., Saunders, N. J., Soneji, S., Palazzo, S., Dunlop, H. M., Cooper, C. D., Brown, P. J., Troussard, X., Mossafa, H. & other authors (2008). MicroRNA expression in lymphocyte development and malignancy. Leukemia 22, 1440–1446.[CrossRef][Medline]
Lee, Y. S. & Dutta, A. (2009). MicroRNAs in cancer. Annu Rev Pathol 4, 199–227.[CrossRef][Medline]
Lin, Y. C., Kuo, M. W., Yu, J., Kuo, H. H., Lin, R. J., Lo, W. L. & Yu, A. L. (2008). c-Myb is an evolutionary conserved miR-150 target and miR-150/c-Myb interaction is important for embryonic development. Mol Biol Evol 25, 2189–2198.
Lowery, A. J., Miller, N., McNeill, R. E. & Kerin, M. J. (2008). MicroRNAs as prognostic indicators and therapeutic targets: potential effect on breast cancer management. Clin Cancer Res 14, 360–365.
Lu, J., Guo, S., Ebert, B. L., Zhang, H., Peng, X., Bosco, J., Pretz, J., Schlanger, R., Wang, J. Y. & other authors (2008). MicroRNA-mediated control of cell fate in megakaryocyte–erythrocyte progenitors. Dev Cell 14, 843–853.[CrossRef][Medline]
Martinez, I., Gardiner, A. S., Board, K. F., Monzon, F. A., Edwards, R. P. & Khan, S. A. (2008). Human papillomavirus type 16 reduces the expression of microRNA-218 in cervical carcinoma cells. Oncogene 27, 2575–2582.[CrossRef][Medline]
Merkerova, M., Belickova, M. & Bruchova, H. (2008). Differential expression of microRNAs in hematopoietic cell lineages. Eur J Haematol 81, 304–310.[CrossRef][Medline]
Morgan, R., Anderson, A., Bernberg, E., Kamboj, S., Huang, E., Lagasse, G., Isaacs, G., Parcells, M., Meyers, B. C. & other authors (2008). Sequence conservation and differential expression of Marek's disease virus microRNAs. J Virol 82, 12213–12220.
Nazerian, K. (1987). An updated list of avian cell lines and transplantable tumours. Avian Pathol 16, 527–544.[CrossRef][Medline]
Ozen, M., Creighton, C. J., Ozdemir, M. & Ittmann, M. (2008). Widespread deregulation of microRNA expression in human prostate cancer. Oncogene 27, 1788–1793.[CrossRef][Medline]
Petherbridge, L., Brown, A. C., Baigent, S. J., Howes, K., Sacco, M. A., Osterrieder, N. & Nair, V. K. (2004). Oncogenicity of virulent Marek's disease virus cloned as bacterial artificial chromosomes. J Virol 78, 13376–13380.
Pfeffer, S., Zavolan, M., Grasser, F. A., Chien, M., Russo, J. J., Ju, J., John, B., Enright, A. J., Marks, D. & other authors (2004). Identification of virus-encoded microRNAs. Science 304, 734–736.
Pfeffer, S., Sewer, A., Lagos-Quintana, M., Sheridan, R., Sander, C., Grasser, F. A., van Dyk, L. F., Ho, C. K., Shuman, S. & other authors (2005). Identification of microRNAs of the herpesvirus family. Nat Methods 2, 269–276.[CrossRef][Medline]
Pratt, Z. L., Kuzembayeva, M., Sengupta, S. & Sugden, B. (2009). The microRNAs of Epstein–Barr virus are expressed at dramatically differing levels among cell lines. Virology 386, 387–397.[CrossRef][Medline]
Rosenfeld, N., Aharonov, R., Meiri, E., Rosenwald, S., Spector, Y., Zepeniuk, M., Benjamin, H., Shabes, N., Tabak, S. & other authors (2008). MicroRNAs accurately identify cancer tissue origin. Nat Biotechnol 26, 462–469.[CrossRef][Medline]
Ruike, Y., Ichimura, A., Tsuchiya, S., Shimizu, K., Kunimoto, R., Okuno, Y. & Tsujimoto, G. (2008). Global correlation analysis for micro-RNA and mRNA expression profiles in human cell lines. J Hum Genet 53, 515–523.[CrossRef][Medline]
Schotte, D., Chau, J. C., Sylvester, G., Liu, G., Chen, C., van der Velden, V. H., Broekhuis, M. J., Peters, T. C., Pieters, R. & Boer, M. L. (2008). Identification of new microRNA genes and aberrant microRNA profiles in childhood acute lymphoblastic leukemia. Leukemia 23, 313–322.[Medline]
Shaner, N. C., Campbell, R. E., Steinbach, P. A., Giepmans, B. N., Palmer, A. E. & Tsien, R. Y. (2004). Improved monomeric red, orange and yellow fluorescent proteins derived from Discosoma sp. red fluorescent protein. Nat Biotechnol 22, 1567–1572.[CrossRef][Medline]
Skalsky, R. L., Samols, M. A., Plaisance, K. B., Boss, I. W., Riva, A., Lopez, M. C., Baker, H. V. & Renne, R. (2007). Kaposi's sarcoma-associated herpesvirus encodes an ortholog of miR-155. J Virol 81, 12836–12845.
Smyth, G. K. (2005). limma: linear models for microarray data. In Bioinformatics and Computational Biology Solutions using R and Bioconductor, pp. 397–420. Edited by R. Gentleman, V. J. Carey, W. Huber, R. A. Irizarry & S. Dudoit. New York: Springer.
Sullivan, C. S. & Grundhoff, A. (2007). Identification of viral microRNAs. Methods Enzymol 427, 3–23.[Medline]
Wiemer, E. A. (2007). The role of microRNAs in cancer: no small matter. Eur J Cancer 43, 1529–1544.[CrossRef][Medline]
Xiao, C., Calado, D. P., Galler, G., Thai, T. H., Patterson, H. C., Wang, J., Rajewsky, N., Bender, T. P. & Rajewsky, K. (2007). MiR-150 controls B cell differentiation by targeting the transcription factor c-Myb. Cell 131, 146–159.[CrossRef][Medline]
Xu, H., Yao, Y., Zhao, Y., Smith, L. P., Baigent, S. J. & Nair, V. (2008). Analysis of the expression profiles of Marek's disease virus-encoded microRNAs by real-time quantitative PCR. J Virol Methods 149, 201–208.[CrossRef][Medline]
Yao, Y., Zhao, Y., Xu, H., Smith, L. P., Lawrie, C. H., Sewer, A., Zavolan, M. & Nair, V. (2007). Marek's disease virus type 2 (MDV-2)-encoded microRNAs show no sequence conservation with those encoded by MDV-1. J Virol 81, 7164–7170.
Yao, Y., Zhao, Y., Xu, H., Smith, L. P., Lawrie, C. H., Watson, M. & Nair, V. (2008). MicroRNA profile of Marek's disease virus-transformed T-cell line MSB-1: predominance of virus-encoded microRNAs. J Virol 82, 4007–4015.
Yeung, M. L., Yassunaga, J., Bennasser, Y., Dusetti, N., Harris, D., Ahmad, N., Matsuoka, M. & Jeang, K. T. (2008). Roles for microRNAs, miR-93 and miR-130b, and tumor protein 53-induced nuclear protein 1 tumor suppressor in cell growth dysregulation by human T-cell lymphotrophic virus 1. Cancer Res 68, 8976–8985.
Zhao, Y., Yao, Y., Xu, H., Lambeth, L., Smith, L. P., Kgosana, L., Wang, X. & Nair, V. (2009). A functional MicroRNA-155 ortholog encoded by the oncogenic Marek's disease virus. J Virol 83, 489–492.
Zhou, B., Wang, S., Mayr, C., Bartel, D. P. & Lodish, H. F. (2007). miR-150, a microRNA expressed in mature B and T cells, blocks early B cell development when expressed prematurely. Proc Natl Acad Sci U S A 104, 7080–7085.
Received 22 December 2008;
accepted 15 March 2009.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |