Previous Article | Next Article ![]()
Journal of Virology, July 2007, p. 7164-7170, Vol. 81, No. 13
0022-538X/07/$08.00+0 doi:10.1128/JVI.00112-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Avian Oncogenic Virus Group, Institute for Animal Health, Compton, Berkshire RG20 7NN, United Kingdom,1 LRF Molecular Haematology Unit, Nuffield Department of Clinical Laboratory Sciences, University of Oxford, Oxford OX3 9DU, United Kingdom,2 Division of Bioinformatics, Biozentrum, University of Basel, Klingelbergstr 50-70, CH-4056 Basel, Switzerland3
Received 17 January 2007/ Accepted 18 April 2007
|
|
|---|
|
|
|---|
Genome sequences of the MDV-1 strains Md5 and GA (19, 36) and the MDV-2 strain HPRS24 (17) show a similar organization composed of two unique sequences, one long (UL) and one short (US), flanked by inverted long and short internal repeats (IRL and IRS, respectively) and by long and short terminal repeats (TRL and TRS, respectively). The similarity of genome structures and sequences is also reflected in the close phylogenetic relationships of these viruses (22). While the Mardivirus lineage is thought to have diverged from the mammalian herpesviruses belonging to the genera Simplexvirus and Varicellovirus about 131 million years ago, MDV-1 and MDV-2 are estimated to have separated only about 26 million years ago (D. J. McGeoch, personal communication). The vast majority of the genes within the UL and US regions of these two viruses are homologous with a high degree of conservation (25). However, many open reading frames (ORFs), particularly those in the IRL/TRL region, are unique to MDV-1. As some of the genes within this region such as the meq and vIL-8 genes are known to be associated with oncogenicity (23, 26), the absence of these ORFs may account for the nononcogenic features of MDV-2.
In addition to the better-known gene regulatory pathways mediated through encoded proteins, noncoding RNAs, particularly microRNAs (miRNAs), are increasingly being recognized as important regulators of gene expression. Because of their tiny structures and nonimmunogenic features, several viruses have exploited this new regulatory mechanism for survival within cells (14, 35). Among the virus-encoded miRNAs, the vast majority have been identified in herpesviruses. The numbers of herpesvirus-encoded miRNAs listed in the latest Release 9.1 of the miRBase database (http://microrna.sanger.ac.uk) include 23 in Epstein-Barr virus (EBV), 16 in rhesus lymphocryptovirus (rLCV), 11 in human cytomegalovirus, 13 in Kaposi's sarcoma herpesvirus, 9 in murine herpesvirus 68, 1 in herpes simplex virus (HSV), and 8 in MDV-1 (13). It has been suggested that miRNA-mediated regulation of gene expression is particularly suited for the distinct biological features of herpesviruses, including nuclear replication and latency (30).
The eight MDV-1-specific miRNAs mapped to the MEQ and the latency-associated transcript (LAT) regions (8), and their expression in primary lymphomas as well as in transformed lymphoblastoid cell lines paralleled the pattern of the MEQ gene expression, suggesting a major role for these miRNAs in MDV latency and transformation. Since MDV-2 is nononcogenic but is antigenically and evolutionarily close to MDV-1, we wanted to find out whether MDV-2 encoded any miRNAs and, if so, whether there was any conservation of miRNA sequences between MDV-1 and MDV-2. For this purpose, we constructed a library using small RNA fractionated from an MSB-1 lymphoblastoid cell line known to be coinfected with the BC-1 strain of MDV-1 and the HPRS24 strain of MDV-2 (1, 15). In this study, we report the identification of novel MDV-2-specific miRNAs from the MSB-1 cell line. We also demonstrate that although there was conservation of the relative genomic positions of MDV-1 and MDV-2 miRNA clusters, there was no sequence conservation despite the close phylogenetic and antigenic relationship between the two viruses.
|
|
|---|
Construction of a small RNA library and sequencing of miRNAs.
The cloning of small RNAs from MSB-1 cells was conducted by following the protocols described previously (29). Briefly, total RNA was extracted from MSB-1 cells using TRIzol reagent (Invitrogen, Paisley, United Kingdom), and 500 µg of total RNA was spiked with 0.5 nM of radiolabeled ([
-32P]ATP [Amersham]) 19- and 24-nucleotide oligoribonucleotides containing the PmeI restriction site (5'-CGUACGCGGGUUUAAACGA-3' and 5'-CGUACGCGGAAUAGUUUAAACUGU-3') and then separated on a 15% denaturing polyacrylamide gel. A gel slice containing RNAs of 19- to 24-nucleotide size was excised and eluted overnight in 0.3 M NaCl at 4°C. RNA was recovered after ethanol precipitation and ligated sequentially to a 5'-end-adenylated 3' adapter oligonucleotide (5'-rAppTTTAACCGCGAATTCCAG/3ddC/-3', where r is ribonucleotide [Integrated DNA Technologies]) and a 5' adapter chimeric DNA/RNA oligonucleotide (5'-ACGGAATTCCTCACTrArArA-3'). Reverse transcription (RT)-PCR was performed with a 3' primer (5'-GACTAGCTGGAATTCGCGGTTAAA-3') and a 5' primer (5'-CAGCCAACGGAATTCCTCACTAAA-3'). The purified PCR products were digested with PmeI to eliminate size marker sequences. In order to introduce BanI restriction sites, a second PCR was performed using the primer pair 5'-CAGCCAACAGGCACCGAATTCCTCACTAAA-3' and 5'-GACTAGCTTGGTGCCGAATTCGCGGTTAAA-3', followed by concatamerization after BanI digestion and T4 DNA ligation. Subsequently, the 3' ends of the concatemers were filled in by incubation for 15 min at 72°C with Taq polymerase in a standard PCR mixture and used directly for ligation into pGEM-T Easy vector (Promega, Southampton, United Kingdom). Selected clones were sequenced and analyzed.
Northern blotting analysis of miRNAs.
Total RNA was extracted from different cells using TRIzol reagent (Invitrogen). Samples of 20 µg total RNA were resolved using a 15% polyacrylamide-1x Tris-borate-EDTA-8 M urea gel and blotted to a GeneScreen Plus membrane (Perkin-Elmer). DNA oligonucleotides with sequences exactly complementary to candidate miRNAs were end labeled with [
-32P]ATP (Amersham, Bucks, United Kingdom) and T4 polynucleotide kinase (New England Biolabs, Herts, United Kingdom) to generate high-specific-activity probes. Hybridization, washing, and autoradiography were carried out as previously described (29).
Immunofluorescence and Western blotting. CEF were cocultivated with MSB-1 cells for 6 to 7 days, and MDV-specific plaques were stained with MDV-1-specific monoclonal antibody (MAb) BD1 (20) or MDV-2-specific MAb Y5 (18). The specific cell populations detected by the two antibodies were visualized with anti-mouse immunoglobulin G (IgG) conjugates of Alexa Fluor 488 or Alexa Fluor 568 (Molecular Probes), using a DM IRB model microscope (Leica Microsystems, Milton Keynes, United Kingdom). For Western blotting, 106 MSB-1 cells, uninfected CEF, or CEF infected with MSB-1 MDV isolates or SB-1 were lysed in protein gel sample buffer (8 M urea, 2% sodium dodecyl sulfate, 10 mM Tris-HCl [pH 6.8], 0.05% bromophenol blue) and separated on a NuPAGE 10% bis-Tris gel (Invitrogen) and transferred onto nitrocellulose membranes using an iBlot gel transfer system (Invitrogen). Western blotting (WB) was performed with MAb Y5, followed by an anti-mouse IgG peroxidase conjugate (Sigma-Aldrich, Dorset, United Kingdom). The membrane was finally developed with an ECL Western blotting analysis system (Amersham).
miRNA sequence accession numbers. MDV-2 miRNAs miRNA-1 to miRNA-17 described here are annotated as mdv2-miR-M14 to mdv2-miR-M30, respectively, in the miRBase Sequence Database (http://microrna.sanger.ac.uk/sequences).
|
|
|---|
Identification of MDV-2-specific miRNAs from the MSB-1 cell line.
For the construction of the small RNA library from MSB-1 cells, a size-fractionated 19- to 24-nucleotide small RNA population was used. After adapter ligation, amplification, and restriction digestion with BanI, the amplicons were cloned as cDNA concatemers into pGEM-T Easy vector. Sequencing of
1,200 clones using vector-specific primers identified a total of 5,099 high-quality reads containing small RNA sequences with both the 5' and 3' adapters. BLAST (2) homology searches of the sequences were performed against the published HPRS24 strain sequence (GenBank accession number AB049735). A total of 518 sequences (about 10% of the total), representing 17 candidate miRNAs, perfectly matched the HPRS24 sequence. In addition, the MSB-1 library also revealed all the recently described (8) MDV-1 miRNAs (not shown).
During the biogenesis of miRNAs, only one of the strands of the
20-bp miRNA duplex is incorporated into the RNA-induced silencing complex (RISC) and guides gene regulation. Even though the non-miRNA passenger strand of the duplex is rapidly degraded (6), these strands could still be captured during cloning. Moreover, for some miRNAs, both strands can be incorporated into the RISC with comparable frequency. In the HPRS24 library, two mature forms, representing both strands of the duplex, were found for 10 out of the 17 candidate miRNAs, increasing the total number of cloned mature MDV-2-encoded miRNAs to 27. In most of these cases, the relative frequency of the passenger strand was much lower than that of the miRNA strand (Table 1). However, in some cases, both mature forms were incorporated with similar frequencies, suggesting the possibility that both strands may be functional and are incorporated into the RISC. Designations of the miRNAs were made sequentially according to their locations in the HPRS24 genome, and the suffix 5p or 3p was added to indicate the 5' or 3' arm, respectively, of the stem-loop precursor from which it was derived (Table 1). Some of the mature miRNAs showed variations in length (Table 1, variations are shown in parentheses). For further validation of HPRS24-encoded candidate miRNAs, each of the sequences including 60 to 80 bases of the surrounding sequence in the viral genome was subjected to MFOLD calculation (39), and secondary structures were drawn according to RNADRAW software (21). As shown in Fig. 1, the miRNA precursors showed an average length of
70 nucleotides and were able to form a hairpin structure. Moreover, 8 of the 17 candidates (Table 1) reside in pre-miRNAs predicted from this virus by using a computational method that was previously used for identification of miRNAs in other herpesviruses (28). These findings support the conclusion that the cloned sequences are MDV-2-encoded miRNAs.
|
View this table: [in a new window] |
TABLE 1. Sequences and genomic positions of MDV-2 miRNAs
|
![]() View larger version (17K): [in a new window] |
FIG. 1. Predicted secondary structures of MDV-2 pre-miRNAs These structure predictions were derived using the MFOLD algorithm (39). The mature miRNAs are indicated in red.
|
![]() View larger version (20K): [in a new window] |
FIG. 2. Genomic location of MDV-1 and MDV-2 miRNAs. (A) Schematic diagram showing the loci in the MDV-2 genome to which the miRNAs (miR) identified in this report map. The TRL and IRL flanking the unique long region as well as the IRS and TRS flanking the unique short region are shown. Genomic positions and the orientation of the ORFs contained in the miRNA loci are indicated. Small arrowheads indicate the locations of the MDV-2 miRNAs. (B) Schematic diagram showing the relative position of the miRNA-containing loci of MDV-1 genome reported recently (8).
|
Demonstration of miRNA expression by Northern blotting. For further confirmation of expression of MDV-2 miRNAs in MSB-1 and MDV-infected CEF, Northern blotting hybridization with individual miRNA probes was carried out. For those miRNAs where both strands were cloned, probes only from the 5p strand was used in Northern blotting. As a negative control for the dual-virus-infected MSB-1 cells, a lymphoblastoid cell line 769T infected only with MDV-1 (RB-1B strain) was used. Uninfected CEF as well as CEF lytically infected with MDV-2 strain HPRS24 were also included in the analysis. Northern blotting detected 15 out of the 17 MDV-2 miRNAs in MSB-1 cells as well as in HPRS24-infected CEF. No miRNAs were detected with RNA extracted from the uninfected CEF or the MDV-1-infected 769T cell line (Fig. 3). There were differences in the levels of expression (based on the intensity of the signals) of miRNAs between the two cell types, with MSB-1 cell line expressing most of the miRNAs at higher levels. Both pre-miRNA and mature miRNAs could be detected in both cell types with Northern blotting. In some cases, such as miRNA-5, miRNA-12-5p, and miRNA-14-5p, the amount of pre-miRNAs was equal to or more abundant than the mature miRNAs, suggesting inefficient processing. This was also evident with miRNA-7, where such differences in the efficiency of processing were evident between the two cell types. In the case of miRNA-15-5p, two bands possibly representing two mature miRNA species were observed with both MSB-1- and HPRS24-infected CEF. In order to confirm that the MDV-2 miRNAs reported here are not unique to just the HPRS24 strain of MDV-2, we also carried out Northern blotting analysis of CEF infected with another MDV-2 strain, SB-1 (33). These studies confirmed that the SB-1 strain also expressed these miRNAs (not shown).
![]() View larger version (82K): [in a new window] |
FIG. 3. Northern blotting analysis for determining the expression of HPRS24 miRNAs. Twenty micrograms of total RNA from MSB-1 (lane 1), pRB-1B virus-transformed 769T cell line (lane 2), HPRS24-infected CEF (lane 3), and uninfected CEF (lane 4) was separated on a 15% denaturing polyacrylamide gel, blotted, and probed with an end-labeled antisense oligonucleotide to the indicated miRNA. Size markers to indicate the positions of the pre-miRNA and the mature miRNA are shown. The cellular U6 snRNA served as the loading control (a representative blot of this set is shown).
|
|
|
|---|
The MDV-transformed lymphoblastoid cell line MSB-1 (1) has been previously reported to be infected with both MDV-1 strain BC-1 and MDV-2 strain HPRS24 (15), and we have confirmed this using IF and WB with specific MAbs. Analysis of the small RNA library from the MSB-1 cell line showed that nearly 10% (518 of 5,099) of the cloned sequences were attributed to 17 candidate HPRS24-encoded miRNAs (Table 1). These candidate sequences conformed to most of the miRNA validation criteria (4) such as (i) their cloning from a cDNA library, (ii) the detection of either strand of the putative double-stranded miRNA precursor, and (iii) their identification by Northern blotting. Nearly all of the predicted miRNAs could also be folded into hairpin structures (Fig. 1). However, some of the miRNAs, such as MDV-2 miRNA-6, which were cloned at least nine times in the library, showed slightly uncharacteristic central bulging and a small terminal bulge using an MFOLD structure prediction. The identification of these novel MDV-2-encoded miRNAs adds to the growing list of herpesvirus-encoded miRNAs.
MDV-2 has a typical alphaherpesvirus genome structure with a colinear arrangement of genes, as in MDV-1 (25). Although both viruses have numerous genes showing conservation in primary sequence and in genomic positions across the genome, several genes located in the IRL/TRL region are unique to each virus. For example, the genes such as meq and vIL-8 that are located in this region and are associated with the oncogenicity of MDV-1 are absent in MDV-2. On the other hand, MDV-2 has several ORFs within this region that are absent in MDV-1. Previous studies have indicated that the IRL/TRL region in MDV-1 is transcriptionally active in latently infected tumors and transformed cells (31, 34). Furthermore, five out of the eight miRNAs identified recently in MDV-1 within this region are also expressed at very high levels in transformed cell lines and tumor cells (8), confirming that this region is transcriptionally active even during stages of latency and transformation. It is remarkable that all but one of the 17 MDV-2-encoded miRNAs reported here are clustered in a 4.2-kb fragment of the IRL/TRL region spanning the putative ORFs R-LORF2, R-LORF3, R-LORF4, and R-LORF5 (Fig. 2). Although there have not been any experimental data to confirm the existence, the transcriptional profile, or the functions of these ORFs, high levels of expression of MDV-2 miRNAs from this region, as demonstrated by direct cloning and Northern blotting, would indicate that this region is transcriptionally active during latency. Compared to the location of the majority of the herpesvirus-encoded miRNAs in the noncoding region of different genes, the identification of 8 out of 17 MDV-2 miRNAs in the coding regions of the ORFs was quite striking. R-LORF2 and R-LORF3 have opposite orientations to R-LORF4 and R-LORF5. The miRNA cluster shares the same orientation as R-LORF5 and R LORF4, with MDV-2-miR-1 located upstream of the R-LORF5 promoter region.
The single other miRNA identified outside the cluster, MDV-2 miRNA-17, is located in the IRS/TRS region. Interestingly, MDV-1 also has three miRNAs located in an equivalent position (8). These three MDV-1 miRNAs map to a large intron in the 5' end of the latency-associated MSR transcript, antisense to the ICP4 gene (10). The MDV-2 miRNA-17 is also in an antisense orientation to ICP4 but overlaps with the coding region, in a region corresponding to amino acids 1,629 to 1,645 that are predicted to form a loop structure between two alpha helices not conserved between MDV-1 and MDV-2. No data for either the existence or for the transcriptional regulation of LAT/ICP4 in MDV-2 have been examined, and until such data are available, the role of this single MDV-2 miRNA in the regulation of MDV-2 latency remains speculative.
Despite the demonstration of relative conservation of miRNA loci in these two closely related viruses, there was no evidence of conservation of primary miRNA sequences of MDV-1 and MDV-2, suggesting their independent evolution. As the majority of the herpesvirus-encoded miRNAs did not show sequence conservation beyond very closely related viruses (9, 24, 28), the lack of sequence homology between MDV-1 and MDV-2 miRNAs was not entirely surprising. Even between closely related EBV and rLCV, the only example where such conservation of miRNA sequences has been demonstrated, only 7 out of 16 miRNAs showed sequence conservation, indicative of a selection probably imposed by the target genes (9). Although MDV-1 and MDV-2 are closely related and separated by only approximately 26 million years, the absence of sequence conservation in any of the miRNA sequences may reflect rapid evolution toward new targets, possibly related to the differences in pathogenicity between the two viruses. Specifically, while MDV-1 is oncogenic, MDV-2 shows other distinctive features such as the induction of cytolytic lesions in immunosuppressed birds (32), protective synergism with serotype 3 vaccines (37), enhancement of lymphoid leukosis (5), and reticuloendotheliosis (3), as well as persistence in transformed B cells (11, 12). Future studies will elucidate whether the virus-encoded miRNAs participate in these processes. Prediction of the putative target genes of virus-encoded miRNAs is particularly difficult, and despite the availability of numerous computational programs for target prediction, targets of only a very few have been validated. Nevertheless, as MDV-2 miRNA-10, miRNA-13, miRNA-14, and miRNA-15 are transcribed antisense to R-LORF2 and R-LORF3 transcripts, these miRNAs could potentially target these transcripts for cleavage as small interfering RNA because of the perfect complementarity of the mRNA sequences. However, as the functions of these novel ORFs are largely unknown at this point, it is difficult to evaluate the impact that the miRNAs could have on the biology of these viruses. Our recent success with the construction of the complete genome of the SB-1 strain of MDV-2 as an infectious bacterial artificial chromosome clone (unpublished data) will allow us to carry out functional analysis of the role of these miRNAs in the biology of MDV-2 in natural infection models.
This work was funded by the Biotechnology and Biological Sciences Research Council (BBSRC), United Kingdom.
Published ahead of print on 25 April 2007. ![]()
|
|
|---|
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»