- Open Access
Comprehensive profiling of Epstein-Barr virus-encoded miRNA species associated with specific latency types in tumor cells
Virology Journalvolume 10, Article number: 314 (2013)
Epstein-Barr virus (EBV) is an etiological cause of many human lymphocytic and epithelial malignancies. EBV expresses different genes that are associated with three latency types. To date, as many as 44 EBV-encoded miRNA species have been found, but their comprehensive profiles in the three types of latent infection that are associated with various types of tumors are not well documented.
In the present study, we utilized poly (A)-tailed quantitative real-time RT-PCR in combination with microarray analysis to measure the relative abundances of viral miRNA species in a subset of representative lymphoid and epithelial tumor cells with various EBV latency types.
Our findings showed that the miR-BHRF1 and miR-BART families were expressed differentially in a tissue- and latency type-dependent manner. Specifically, in nasopharyngeal carcinoma (NPC) tissues and the EBV-positive cell line C666-1, the miR-BART family accounted for more than 10% of all detected miRNAs, suggesting that these miRNAs have important roles in maintaining latent EBV infections and in driving NPC tumorigenesis. In addition, EBV miRNA-based clustering analysis clearly distinguished between the three distinct EBV latency types, and our results suggested that a switch from type I to type III latency might occur in the Daudi BL cell line.
Our data provide a comprehensive profiling of the EBV miRNA transcriptome that is associated with specific tumor cells in the three types of latent EBV infection states. EBV miRNA species represent a cluster of non-encoding latency biomarkers that are differentially expressed in tumor cells and may help to distinguish between the different latency types.
Increasing evidence suggests that EBV is an oncogenic human virus that causes many malignancies [1–3]. The EBV-encoded LMP1 gene is considered a viral oncogene that leads to the development of several tumor types [4, 5]. Viral microRNA species have been found in EBV-infected cells  and recent studies indicate that 44 EBV miRNA species are involved in the establishment of viral latency and tumorigenesis [7–12]. However, the comparable profiling of EBV miRNA transcriptomes in tumor cells in the three types of latency states has not been performed. In this study, we investigated the comprehensive expression patterns of EBV miRNA species in cells with different types of latency and emphasized the viral miRNA patterns in NPC tissues and cells in the EBV latency II state.
MiRNAs play physiological and pathophysiologic roles in cell proliferation, differentiation and apoptosis by regulating protein expression in a post-transcriptional manner. Furthermore, miRNA regulation is an innate defense mechanism against viruses in host cells. Over the course of evolution, a subset of viruses, including the herpesviruses, has obtained the ability to hijack this antiviral mechanism of host cells to help maintain a persistent infection [13–16]. Approximately 20% of human tumors are associated with latent viral infections, and this latency is controlled by various viral genetic and epigenetic regulatory mechanisms [15, 16]. Virus-encoded miRNAs play significant roles in maintaining persistent EBV infections by down-regulating the expression of viral immediate-early genes [7, 17–19] and promoting the tumorigenesis of host cells. In addition, these viral miRNAs help the virus evade host immune defenses [9, 11, 20, 21] by regulating host gene expression at the post-transcriptional level.
The classification of EBV latency is an academic and clinical issue because latent infections of different types are associated with different tumors . The traditional classification system is based on only a few molecular markers of EBV latency, including EBERs, EBNA1-6, LMP1, LMP2 and BARTs [2, 22–25], and typing of EBV latency types has limitations. Since EBV was first found to encode viral miRNAs in 2004 , 44 mature, EBV-encoded miRNA species have been registered in the Sanger miRBase library (Release 16.0). EBV miRNA species are produced from a large intron in BARTs (BamHI A rightward transcripts) prior to splicing [26, 27]. Therefore, these novel, abundant markers of viral latency may help in obtaining further insight into EBV and may be useful as novel, biomarkers for identifying latency types [12, 28].
MiRNAs are a class of difficult-to-detect RNA strands that are 19–23 nt in length. These RNAs can currently be analyzed by high-throughput microarray analysis [29, 30], SOLEXA sequencing , Northern blotting  and quantitative real-time RT-PCR . Given the higher sensitivity and specificity of qPCR , in this study, we mainly used a poly(A)-tailed, qRT-PCR technique in combination with a microarray assay to analyze the expression patterns of the EBV miRNA transcriptomes that are associated with specific latency types in various tumor cells. Additionally, we primarily evaluated a type I to III transition based on clustering analysis of the EBV miRNA species.
Quantitative analysis of viral-coding transcripts in EBV-infected cell lines and NPC tissues
The previous identification of the EBV latency types was mainly based on viral EBERs, nuclear antigens, LMP1, LMP2 and the Bam H1 rightward transcripts [2, 22–25]. EBNA1 is expressed under the alternative promoters Wp, Cp and Qp in a latency type-specific manner [22, 24, 35], and LMP1 exhibits differential abundance when compared to the various latency types , suggesting that the activities of the three latency promoters that control EBNA1 expression and LMP1 abundance can characterize the three latency types. The essential latency characteristics described above are partially illustrated in Figure 1. In the present study, we measured the activities of the Cp/Qp/Wp promoters and the expression level of LMP1 to characterize the different types of latent EBV infections.
We examined the expression levels of the traditional latency biomarkers by quantitative real-time RT-PCR in a subset of typical tumor cell lines and NPC tissues with different EBV latency types. During latency I (Daudi and Akata(+) cells) [36, 37] and latency II (C666-1 and T4-T7 cells) , Qp exhibited a much higher activity than in latency III (B95-8, Raji and Namalwa cells). When compared with high-stringent latency I (Akata(+) cells), latency II (C666-1 and T4-7 cells) and latency III (B95-8 and Raji cells) were associated with slightly higher expression levels of LMP1. During low-stringent latency I (Daudi), the expression of the latency II/III biomarker LMP1 appeared to be elevated. Interestingly, viral latency biomarkers were nearly undetected by both methods in Namalwa cells, implying that the majority of the EBV episomes might be lost during culture (Figure 2).
The comprehensive expression patterns of EBV miRNA species in three types of latency states
EBV miRNA species are expressed as two independent transcripts, BHRF1 and BART. The miR-BHRF1 family is located within the introns of BHRF1 (Bam HI fragment H rightward open reading frame 1) [10, 39, 40], whereas the miR-BART family is located in the intronic regions of the BART and form two separate clusters (cluster I and cluster II) [26, 27]. The prototypic EBV strain that was derived from the B95-8 cell line has a 12-kb deletion that spans the two BART clusters from miR-BART-5 to 14 [41, 42]. A detailed overview of the genomic locations of the EBV miRNA genes is presented in Figure 1.
In our investigation, we mainly used real-time RT-PCR to gain insight into the comparable profiles of the EBV miRNA transcriptomes of the three latency types. To confirm the specificity of the PCR products and verify the reliability of this poly(A)-tailed, real-time qPCR method, we randomly sequenced three qPCR products, submitted them to AT cloning procedures, and the results of DNA sequencing are shown in Figure 3C. Following qPCR analysis, we further evaluated the expression levels of EBV miRNAs using a microarray assay, and the scanning images are shown in Figure 3D. In addition, correlation analysis showed statistically significant correlations between the two methods for the latency I (Akata(+) cells; Pearson’s correlation coefficient = 0.51, p = 0.0004), latency II (C666-1 cells; Pearson’s correlation coefficient = 0.73, p < 0.0001) and latency III profiles (Raji cells; Pearson’s correlation coefficient = 0.47, p = 0.0024) (Figure 3E).
The quantitative analysis of the EBV miRNAs in each latency type yielded a comprehensive expression pattern of the EBV miRNA species (Figure 3A). Generally, the BHRF1 family exhibited a nearly 100-fold higher abundance in latency III (B95-8 and Raji cells) than in latency I (Akata(+) and Daudi cells), whereas the expression of the BHRF1 family was nearly absent in latency II (C666-1 cells), suggesting that the expression of the BHRF1 family is latency-III dependent. In addition, the BART family was differentially expressed in cells with the three latency types (Akata(+), Daudi, C666-1, Raji and B95-8 cells), with the exception of Namalwa cells, suggesting that the BART family plays differential roles in tumors with the three types of latency. Notably, the average expression level of the BART family in the epithelial tumor cell line C666-1, which is a latency II cell line, was the highest of all latency types, indicating that the BART family may contribute more to NPC tumorigenesis. In contrast, the expression level of the BHRF1 family was significantly higher than that of the BART family in the lymphoid tumor cell line Raji, which is a latency III cell line, suggesting that the BART and BHRF1 families are somewhat tissue-specific (Figure 3B). Interestingly, the BART and BHRF1 families were nearly undetectable in Namalwa cells, suggesting that EBV episomes might have been lost by these host cells in vitro.
Although a large difference between the expression levels of the BHRF1 and BART families was found in the same cell line, the difference between BART clusters I and II was relatively small in Akata(+), Daudi, C666-1 and Raji cells. This result suggests that BART clusters I and II are spliced from a common transcript, as was reported previously. Notably, 31 members of the EBV miRNA species that were barely detectable by qPCR are located in the 12-kb deletion in the prototypic B95-8-derived EBV strain that spans nts 139724–151554 of the EBV genome (GenBank accession number AJ507799.2). A detailed schematic of the deletion region in the B95-8-derived strain is shown in Figure 1. A subset of B95-8 cells can produce infectious EBV virions that are released into the medium during culture, which is a sign that the EBV life cycle has switched from latency to lysis. This switch might be related to the miRNAs that are encoded in the deleted region. That is, the miRNAs in the 12-kb deletion region could be involved in maintaining EBV latency. As for Akata(+) cells, which are high-stringent latency I cells, the abundance of the BHRF1 family was 10- to 100-fold lower than that of the BART family and was similar to its abundance in the high-stringent latency II cell line C666-1 (Figure 3A). Recent studies on the EBV methylome have found that hypermethylation of the promoter of the BHRF1 gene leads to its low expression in Akata cells, which is consistent with our present results .
EBV miRNA transcriptome in primary NPC tissues and C666-1 cells with type II latency
NPC tissues and the NPC cell line C666-1, which was the only cell line that persistently carried EBV, exhibited high-stringent, type II EBV latency [2, 25, 38, 44]. In this study, we explored the EBV miRNA expression status of pathologically confirmed, NPC biopsy specimens, and our results indicated that the BHRF1 family was expressed at a low level, whereas the BART family showed much higher abundance in NPC samples and C666-1 cells (Figure 4A). The Pearson and Spearman correlation analyses of BART miRNA expression of these NPC samples and the C666-1 cell line showed great similarity, suggesting that the mechanism by which the EBV BART family regulates latency is the same in different NPC tumors (Figure 4B). The average expression level of the BART family in the NPC tissue samples also showed a strong correlation with its expression level in the NPC cell line C666-1 (R = 0.89, p < 0.0001) (Figure 4B). These results show that the EBV miRNA expression level in C666-1 cells is consistent with its pattern in NPC tissues, indicating that the C666-1 cell line is an ideal model for investigating the relationship between EBV and NPC. In addition, results of previous investigations on the pattern of EBV miRNA expression in NPC tissues using a stem-loop, RT primer-based qPCR assay are generally in agreement with our present results. We also evaluated the ratio of viral miRNAs to human miRNAs using the microarray data, and we found that the ratio was greater than 10% (Figure 4C). This result was also consistent with those of previous studies [31, 45, 46] and further indicates that EBV miRNAs might play important roles in NPC tumorigenesis.
Evaluation of a possible transition from type I to III latency in Daudi cells that is dependent on viral miRNAs
The viral miRNAs that are differentially expressed in different types of EBV-infected cells are expected to help distinguish between the different latency types. In this study, we performed unsupervised, clustering analysis based on the qPCR dataset to examine the two types of indicators that are used to type latency. Interestingly, the human Burkitt’s lymphoma cell line Daudi [36, 37], which was classified into the latency I group based on latency markers (Cp, Wp, Qp and LMP1), was classified into the latency III group based on the EBV miRNA latency markers (Figure 5A). Additionally, clustering analyses based on non-encoding miRNA species and their encoding latent genes categorized Daudi cells into the latency III group (Figure 5A). To evaluate the accuracy of the latter two clustering results, we performed phenotypic and genotypic investigations. Based on phenotypic characteristics (shown in Figure 5B), we observed that, in culture, Daudi cells grew as macroscopic clumps that were similar to those of the typical latency III cell line B95-8, whereas the high-stringent, latency I Akata(+) cells grew as a single-cell suspension, which is a typical characteristic of BL cells in latency I. The same phenomenon regarding the Daudi cell line was also observed in a recent study , suggesting that the EBV latency type might have switched from type I to III during the course of long-term culturing of Daudi cells. Additionally, we examined a latency type-specific change at the viral gene level when comparing high-stringent latency I Akata(+) cells, which have no LMP1 expression and lower Wp activity, to Daudi cells, which began to express the type III-associated marker LMP1 and had elevated Wp activity despite the persistence of higher Qp activity (Figure 5C). In addition to the changes in the encoding latency biomarkers, changes in the expression levels of viral non-coding biomarkers, including the BHRF1 and BART families, occurred in Daudi cells. In the representative type I cell line Akata(+), the average expression level of the BART family was much higher than that of the BHRF1 family, whereas the BHRF1 family had a slightly higher abundance than the BART family in Daudi cells, indicating that a transition from type I to type III latency had possibly occurred in the tested Daudi cells (Figure 5C). It seems that the limited number of traditional viral latency biomarkers is inadequate to clearly distinguish between the differential EBV latency types, whereas a cluster of differentially expressed miRNA species maybe reinforce these demands.
Characterization of latent EBV coding and noncoding transcripts in cells and NPC tissues
The different EBV latency types are closely related to various malignancies, and, consequently, the correct identification of latency types is critical for understanding the virus’s contribution to carcinogenesis. EBV-encoded miRNAs were first found in 2004 , and, currently, as many as 44 EBV miRNA species have been found and registered in the Sanger miRBase. In this study, which is based on our comprehensive and comparable investigation of most EBV latency biomarkers, including coding and noncoding transcripts of the three latency types, we provide additional information on the in-depth characterization of EBV latent infections, as summarized in Table 1. In light of our analysis of viral latency biomarkers, these viral, non-encoding miRNA species should be considered when classifying the latency types of EBV-infected tumor cells.
In light of the roles of human miRNAs under physiological and pathological conditions, especially in the process of carcinogenesis , the discovery of EBV-encoded miRNAs is an important event and a milestone in molecular virology that provided us with new insight into the relationship between a virus and its host cell. Earlier investigations showed that, in addition to EBV, other viruses such as HSV1/2 [19, 48], HCMV , SV40 , KSHV  and even HIV  encode their own miRNAs, suggesting that, during long-term evolution, some viruses have hijacked this innate antivirus mechanism to serve themselves or to manipulate human cellular miRNAs to target viral genes . Theoretically, the ability of a subset of viruses to encode their own miRNAs is an adaptive response to environmental stress because of the low immunogenic risk to the virus that is posed by those noncoding miRNAs [16, 51].
Epstein-Barr virus infects over 90 percent of adults worldwide, and EBV infections are strongly associated with multiple human diseases, especially some cancers. In tumor cells, EBV genomes exist as episomes in a state of latency that does not produce lytic viral particles that can be detected by the host immune system. The EBV genome encodes approximately 90 genes but only expresses a few genes, such as EBERs, EBNA1-6, LMP1, LMP2 and BARTs , in latent stages, and these genes serve as markers for classical latency typing [23–25]. Based on the differential expression patterns of these viral genes, latent infections are divided into three types. In addition, tumors have different latency types, and the latency type can change during in vivo and in vitro culture [12, 52]. Since the identification of EBV-encoded miRNAs in 2004 , researchers have increasingly focused their attention on investigating the relationships between viral miRNAs and the establishment of EBV latency and host carcinogenesis.
In this study, we analyzed the comprehensive expression patterns of EBV miRNA species using quantitative real-time PCR and microarray assays in representative tumor cell lines with the three latency types. The cell lines included the type I latency cell lines Akata(+) and Daudi, the type II latency cell line C666-1 and the type III latency cell lines Namalwa, Raji and B95-8, and we analyzed several NPC tissues. In addition, we explored the profile of the EBV miRNA transcriptome in-depth. Generally, the expression of the BHRF1 miRNA family is dependent on the type of viral latency and the cell’s histological origin, whereas the BART miRNA family is expressed differentially in cells during all forms of latency. Moreover, qPCR analysis indicated a difference of nearly 100-fold between all tested cell lines, and the expression of individual BART miRNA species within one cell line could differ by 50-fold or more, although those miRNAs were transcribed together as a single, primary transcript . In contrast, the BART family was expressed constantly at a high level without a significant difference between its expression in EBV-positive NPC tissues and the C666-1 cell line (latency II), suggesting that the BART family is involved in maintaining EBV latency and promoting NPC tumorigenesis. In fact, recent research has demonstrated that the BART family plays important roles in the establishment of EBV latency and NPC initiation and development. Noticeably, ebv-miR-bart6-5p can target the human dicer gene to control the latency switch from type I/II to type III . Briefly, our results and emerging evidence indicate that EBV-encoded miRNA species have a functional role in target gene regulation during carcinogenesis in host cells.
NPC, which is prevalent in southern China, is a typical model of type II EBV latency, in which only three protein-coding genes are expressed: EBNA1, LMP1 and LMP2A. Our results indicate that the noncoding BART family was expressed highly and accounted for more than 10% of all known human miRNAs, suggesting that, like the viral oncogene LMP1, these abundant noncoding miRNA species are likely involved in NPC pathogenesis. As in other viruses, the use of miRNAs is one of the strategies used by the virus to evade host-immune defenses because noncoding miRNAs, which have no immunogenicity, cannot trigger the body’s immune response. Therefore, the production of miRNAs is a relatively low-risk and efficient regulatory mechanism for maintaining a persistent infection and stability of its genome . In addition, we found that the EBV miRNA transcriptome was somewhat stable in the NPC tissues and the type II latency cell line C666-1, which further suggested that the EBV BART family plays key roles in virus latency and host tumorigenesis . The results also imply that the NPC cell line C666-1 is an appropriate cellular model to examine the EBV infection status in vivo and is an ideal cell line for exploring the relationship between EBV and NPC .
It is generally agreed that Burkitt’s lymphoma exhibits EBV type I latency in vivo and that NPC exhibits type II latency . Latency I and II are low-risk types because fewer immunogenic viral proteins are expressed during those times than during high-risk latency III . EBV-positive cells cultured in vitro are not under surveillance by the immune system of the host; therefore, these cells experience little or no immune pressure, which may result in the loss of the virus or a change in its latency type. This could explain why EBV is missing from nearly all NPC cell lines, with the exception of the C666-1 cell line. Similarly, lymphoma cells in long-term culture occasionally lose EBV episomes, as was observed in the Namalwa BL cell line in this study, where the BHRF1 and BART families were expressed at very low levels (Figure 3A and D). Moreover, EBV can spontaneously switch latency types in some lymphoma cells in vitro, and the most common switch is a type I-to-type III transition . However, certain subtle latency transitions are difficult to detect and identify due to the absence of more specific and sensitive latency biomarkers. Because EBV miRNA species are a cluster of novel latency markers consisting of 44 members, including the BHRF1 family and the BART family, that are differentially expressed in host cells, we profiled the patterns of the EBV miRNA transcriptome in tumor cells with the three different latency types, and we expected that clustering analysis would distinguish between the three latency types. We primarily observed a latency switch from type I to type III in the Daudi BL cell line, which was previously regarded to be a cell line in the latency I state, during the clustering analysis, and we further investigated this switch by phenotypic and genotypic analyses. It seems that these numerous, noncoding viral miRNA species might be informative and valuable for determining latency type.
Viral noncoding miRNAs seem to be preferential regulatory tools that are extensively utilized by EBV due to their low immunogenicity. Thus, EBV is a good model for studying the role of viral miRNA regulation in the maintenance of latency. Our results suggest that these differentially expressed viral miRNAs are molecular biomarkers of latent infection and may be useful in identifying latency types. In particular, the BART family contributes greatly to epithelial tumors and to some lymphomas, whereas the BHRF1 family is primarily responsible for a subset of lymphocyte tumors. Our data provide fine profiles of the EBV miRNA transcriptomes that are associated with the three EBV latency types and suggest that different viral miRNA species are involved in various host-cell carcinogenesis processes.
Materials and methods
Cell culture and tissues
The cell lines and tissues used in this study are listed in Table 2. The Daudi, Raji and B95-8 cell lines were obtained from the American Type Culture Collection (http://www.atcc.org). The Namalwa, Akata, Ramos and CNE2 cell lines were kindly provided by Prof. Musheng Zeng (Cancer Center, Sun Yat-Sen University), and the C666-1 cell line was a gift from The Chinese University of Hong Kong (CUHK, Hong Kong SAR, China) . All of the cells were maintained in RPMI media with 10% fetal bovine serum (Gibco, USA) and 1% penicillin-streptomycin and were incubated in 5% CO2 at 37°C in a humidified atmosphere. The primary NPC tissues used in this research were biopsies obtained from the Cancer Center at Sun Yat-Sen University (Guangzhou, China) prior to any treatments and were pathologically confirmed as nondifferentiated non-keratinizing carcinomas (UNKC, WHO type III). This project was performed in agreement with the clinical ethics guidelines, and all clinical samples were collected and analyzed with prior written, informed consent from the patients. This project was approved by the Research Ethics Committee of the Sun Yat-Sen University Cancer Center.
Small RNA-enriched total RNA isolation
Total RNA containing low-molecular-weight RNAs was prepared using the TRIzol method with minor modifications. Briefly, to extract more small RNAs, following the first precipitation, the resulting upper aqueous phase was mixed with 0.5 ml pre-cooled, isopropyl alcohol per 1 ml of TRIzol Reagent, and this mixture was then placed at −20°C for 45 min and then spun at 12,000 × g at 4°C for 45 min. The pellets were washed in 75% ethanol and then allowed to air dry for 5 min. The RNA pellets were resuspended in nuclease-free water, quantified by UV spectrophotometry (Nanodrop, Thermo, USA), aliquoted and stored at −80°C.
Quantitative RT-PCR analysis of viral-coding gene expression
The three EBNA1 promoter activity patterns (known as Wp, Cp and Qp) and LMP1 expression were analyzed by quantitative RT-PCR assays using an ABI7900HT Fast Real-Time PCR System (Applied Biosystems, USA). The housekeeping gene GAPDH was used as an internal control to verify the RNA quality and loading accuracy. Reverse transcription of 2 μg of total RNA was performed using M-MLV reverse transcriptase (Promega, USA) according to the manufacturer’s guidelines. The EBV-negative tumor cell lines Akata(−), CNE2 and Ramos were used as negative controls for type I, II and III latency, respectively. In the qPCR assay, each reaction was run in triplicate tubes in the presence of SYBR green I, and the raw Ct values were normalized to that of the housekeeping gene GAPDH. The PCR primer sequences for all analyzed genes are listed in Additional file 1: Table S1. The DNA oligonucleotides that were used as PCR primers and the chip probes in the study were synthesized and HPLC-purified by Invitrogen Inc. (Guangzhou, China).
Poly (A)-tailed, quantitative real-time RT-PCR analysis of the EBV miRNA species
The All-in-One miRNA qRT-PCR Detection Kit (GeneCopoeia, USA) was used to determine the miRNA levels in the samples. Briefly, poly(A) polymerase was used to add poly(A) tails to the 3′ ends of the miRNAs, and M-MLV RTase and a unique oligo-dT adaptor primer were used to transcribe the poly(A)-tailed miRNAs to generate a cDNA library. The All-in-One qPCR Mix containing SYBR green I along with the miRNA-specific primers and the general primer were used to produce the respective cDNAs for all EBV-encoded miRNA species using an ABI7900HT Fast Real-Time PCR System (Applied Biosystems, USA). Thermal cycling was performed as follows: initial denaturation at 95°C for 10 min; 40 cycles of denaturation at 95°C for 10 s, annealing at 60°C for 20 s and extension at 72°C for 10 s; and a final extension at 72°C for 10 min (GeneCopoeia, USA). The specificity of the qPCR products was evaluated and then verified by melting curve analysis, gel electrophoresis and AT cloning-based DNA sequencing (pMD18-T, TAKARA, Japan). The genomic locations of the EBV miRNA genes are illustrated in Figure 1, and all of the miRNA-specific forward primers that were used are listed in Additional file 1: Table S2. The EBV-negative NPC cell line CNE2 was used as a negative control, and the housekeeping gene U6 was used as an internal control for normalization. Each qPCR reaction was run in triplicate, and the relative expression levels were calculated using the comparative method (2-ΔΔCt) .
Microarray evaluation of human and viral miRNA expression
The oligonucleotide microarray for the analysis of microRNAs was fabricated in-house, as described previously with minor modifications . Briefly, pre-cleaned, Gold Seal slides (Thermo Fisher Scientific Inc.) were washed with detergents and baked at 140°C for 4 hours, and then the probes (40 nM final concentration) were mixed with printing buffer at a 1:4 ratio and were printed onto the slides at a constant temperature and humidity using a SmartArrayer 136 printer (CapitalBio Inc., Beijing, China). The microRNA probe set included human and EBV miRNAs that were registered in Sanger miRBase Release 12.0 (http://www.mirbase.org), and several positive and negative controls were designed and labeled, in part, according to other published documents. The EBV miRNA probe sequences are listed in Additional file 1: Table S2. The probe labeling and hybridization were performed, in part, as described previously with minor modifications [30, 61]. Briefly, 2.5 μg of total RNA was labeled with 2 nmol of pCp-DY647 or pCp-DY547 (Dharmacon, USA) using 15 units of T4 RNA ligase (GE Healthcare, USA) in a total reaction volume of 20 μL at 16°C overnight. The resulting 2.5 μg of labeled, small RNA-enriched total RNA was purified using Bio-Spin 6 desalting columns (BIO-RAD, USA) and then condensed in vacuum. These miRNAs were hybridized to pretreated chips in 1× hybridization solution (5× Denhardt’s solution, 0.5% SDS, 5× SSC) in a hybridization chamber at 46°C for 24 hours, which was followed by slide washing and scanning with a LuxScan™ 10 K Microarray Scanner (CapitalBio Inc., China). Subsequently, the scanned images were analyzed using GenPix Pro 6.0 (Axon Instruments, U.S.A.), all data were MIAME compliant, and the raw data have been deposited into a MIAME-compliant database (GSE33225).
Given the higher sensitivity, specificity and reliability of qPCR analysis compared to microarray evaluation, we performed unsupervised clustering analysis using the normalized qPCR data [28, 62]. The tumor cell lines and tissues with different EBV latency types were clustered using TMEV (TIGR multiarray experiment viewer, version 4.2.01, http://www.tm4.org/) based on the EBV miRNA expression levels that were quantified by poly(A)-tailed, quantitative real-time RT-PCR, and the data were log2 transformed. The hierarchical clustering of the EBV miRNA species and the classical typing markers, including Cp, Wp and Qp activity and LMP1 expression, was performed using average linkage clustering, and the Euclidian distance matrix was calculated.
Baumforth KR, Young LS, Flavell KJ, Constandinou C, Murray PG: The Epstein-Barr virus and its association with human cancers. Mol Pathol 1999, 52: 307-322. 10.1136/mp.52.6.307
Niedobitek G: Epstein-Barr virus infection in the pathogenesis of nasopharyngeal carcinoma. Mol Pathol 2000, 53: 248-254. 10.1136/mp.53.5.248
Pathmanathan R, Prasad U, Sadler R, Flynn K, Raab-Traub N: Clonal proliferations of cells infected with Epstein-Barr virus in preinvasive lesions related to nasopharyngeal carcinoma. N Engl J Med 1995, 333: 693-698. 10.1056/NEJM199509143331103
Lo AK, Lo KW, Tsao SW, Wong HL, Hui JW, To KF, Hayward DS, Chui YL, Lau YL, Takada K, Huang DP: Epstein-Barr virus infection alters cellular signal cascades in human nasopharyngeal epithelial cells. Neoplasia 2006, 8: 173-180. 10.1593/neo.05625
Lo AK, To KF, Lo KW, Lung RW, Hui JW, Liao G, Hayward SD: Modulation of LMP1 protein expression by EBV-encoded microRNAs. Proc Natl Acad Sci USA 2007, 104: 16164-16169. 10.1073/pnas.0702896104
Pfeffer S, Zavolan M, Grasser FA, Chien M, Russo JJ, Ju J, John B, Enright AJ, Marks D, Sander C, Tuschl T: Identification of virus-encoded microRNAs. Science 2004, 304: 734-736. 10.1126/science.1096781
Barth S, Pfuhl T, Mamiani A, Ehses C, Roemer K, Kremmer E, Jaker C, Hock J, Meister G, Grasser FA: Epstein-Barr virus-encoded microRNA miR-BART2 down-regulates the viral DNA polymerase BALF5. Nucleic Acids Res 2008, 36: 666-675.
Choy EY, Siu KL, Kok KH, Lung RW, Tsang CM, To KF, Kwong DL, Tsao SW, Jin DY: An Epstein-Barr virus-encoded microRNA targets PUMA to promote host cell survival. J Exp Med 2008, 205: 2551-2560. 10.1084/jem.20072581
Nachmani D, Stern-Ginossar N, Sarid R, Mandelboim O: Diverse herpesvirus microRNAs target the stress-induced immune ligand MICB to escape recognition by natural killer cells. Cell Host Microbe 2009, 5: 376-385. 10.1016/j.chom.2009.03.003
Seto E, Moosmann A, Gromminger S, Walz N, Grundhoff A, Hammerschmidt W: Micro RNAs of Epstein-Barr virus promote cell cycle progression and prevent apoptosis of primary human B cells. PLoS Pathog 2010, 6: e1001063. 10.1371/journal.ppat.1001063
Stern-Ginossar N, Elefant N, Zimmermann A, Wolf DG, Saleh N, Biton M, Horwitz E, Prokocimer Z, Prichard M, Hahn G, et al.: Host immune system gene targeting by a viral miRNA. Science 2007, 317: 376-381. 10.1126/science.1140956
Xia T, O’Hara A, Araujo I, Barreto J, Carvalho E, Sapucaia JB, Ramos JC, Luz E, Pedroso C, Manrique M, et al.: EBV microRNAs in primary lymphomas and targeting of CXCL-11 by ebv-mir-BHRF1-3. Cancer Res 2008, 68: 1436-1442. 10.1158/0008-5472.CAN-07-5126
Hansen A, Henderson S, Lagos D, Nikitenko L, Coulter E, Roberts S, Gratrix F, Plaisance K, Renne R, Bower M, et al.: KSHV-encoded miRNAs target MAF to induce endothelial cell reprogramming. Genes Dev 2010, 24: 195-205. 10.1101/gad.553410
Mahajan VS, Drake A, Chen J: Virus-specific host miRNAs: antiviral defenses or promoters of persistent infection? Trends Immunol 2009, 30: 1-7. 10.1016/j.it.2008.08.009
Speck SH, Ganem D: Viral latency and its regulation: lessons from the gamma-herpesviruses. Cell Host Microbe 2010, 8: 100-115. 10.1016/j.chom.2010.06.014
Swaminathan S: Noncoding RNAs produced by oncogenic human herpesviruses. J Cell Physiol 2008, 216: 321-326. 10.1002/jcp.21480
Bellare P, Ganem D: Regulation of KSHV lytic switch protein expression by a virus-encoded microRNA: an evolutionary adaptation that fine-tunes lytic reactivation. Cell Host Microbe 2009, 6: 570-575. 10.1016/j.chom.2009.11.008
Murphy E, Vanicek J, Robins H, Shenk T, Levine AJ: Suppression of immediate-early viral gene expression by herpesvirus-coded microRNAs: implications for latency. Proc Natl Acad Sci U S A 2008, 105: 5453-5458. 10.1073/pnas.0711910105
Umbach JL, Kramer MF, Jurak I, Karnowski HW, Coen DM, Cullen BR: MicroRNAs expressed by herpes simplex virus 1 during latent infection regulate viral mRNAs. Nature 2008, 454: 780-783.
Cullen BR: Viral and cellular messenger RNA targets of viral microRNAs. Nature 2009, 457: 421-425. 10.1038/nature07757
Sullivan CS, Grundhoff AT, Tevethia S, Pipas JM, Ganem D: SV40-encoded microRNAs regulate viral gene expression and reduce susceptibility to cytotoxic T cells. Nature 2005, 435: 682-686. 10.1038/nature03576
Brooks LA, Lear AL, Young LS, Rickinson AB: Transcripts from the Epstein-Barr virus BamHI a fragment are detectable in all three forms of virus latency. J Virol 1993, 67: 3182-3190.
Crawford DH: Biology and disease associations of Epstein-Barr virus. Philos Trans R Soc Lond B Biol Sci 2001, 356: 461-473. 10.1098/rstb.2000.0783
Szeles A, Falk KI, Imreh S, Klein G: Visualization of alternative Epstein-Barr virus expression programs by fluorescent in situ hybridization at the cell level. J Virol 1999, 73: 5064-5069.
Young LS, Dawson CW, Eliopoulos AG: The expression and function of Epstein-Barr virus encoded latent genes. Mol Pathol 2000, 53: 238-247. 10.1136/mp.53.5.238
Chen H, Huang J, Wu FY, Liao G, Hutt-Fletcher L, Hayward SD: Regulation of expression of the Epstein-Barr virus BamHI-a rightward transcripts. J Virol 2005, 79: 1724-1733. 10.1128/JVI.79.3.1724-1733.2005
Edwards RH, Marquitz AR, Raab-Traub N: Epstein-Barr virus BART microRNAs are produced from a large intron prior to splicing. J Virol 2008, 82: 9094-9106. 10.1128/JVI.00785-08
Cameron JE, Fewell C, Yin Q, McBride J, Wang X, Lin Z, Flemington EK: Epstein-Barr virus growth/latency III program alters cellular microRNA expression. Virology 2008, 382: 257-266. 10.1016/j.virol.2008.09.018
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature 2005, 433: 769-773. 10.1038/nature03315
Thomson JM, Parker J, Perou CM, Hammond SM: A custom microarray platform for analysis of microRNA gene expression. Nat Methods 2004, 1: 47-53. 10.1038/nmeth704
Chen SJ, Chen GH, Chen YH, Liu CY, Chang KP, Chang YS, Chen HC: Characterization of Epstein-Barr virus miRNAome in nasopharyngeal carcinoma by deep sequencing. PLoS One 2010,5(9):e12745. 10.1371/journal.pone.0012745
Cai X, Schafer A, Lu S, Bilello JP, Desrosiers RC, Edwards R, Raab-Traub N, Cullen BR: Epstein-Barr virus microRNAs are evolutionarily conserved and differentially expressed. PLoS Pathog 2006, 2: e23. 10.1371/journal.ppat.0020023
Reichenstein I, Aizenberg N, Goshen M, Bentwich Z, Avni YS: A novel qPCR assay for viral encoded microRNAs. J Virol Methods 2010, 163: 323-328. 10.1016/j.jviromet.2009.10.018
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc 2008, 3: 1101-1108. 10.1038/nprot.2008.73
Kelly GL, Long HM, Stylianou J, Thomas WA, Leese A, Bell AI, Bornkamm GW, Mautner J, Rickinson AB, Rowe M: An Epstein-Barr virus anti-apoptotic protein constitutively expressed in transformed cells and implicated in burkitt lymphomagenesis: the Wp/BHRF1 link. PLoS Pathog 2009, 5: e1000341. 10.1371/journal.ppat.1000341
Kitagawa N, Goto M, Kurozumi K, Maruo S, Fukayama M, Naoe T, Yasukawa M, Hino K, Suzuki T, Todo S, Takada K: Epstein-Barr virus-encoded poly(A)(−) RNA supports Burkitt’s lymphoma growth through interleukin-10 induction. EMBO J 2000, 19: 6742-6750. 10.1093/emboj/19.24.6742
Pratt ZL, Kuzembayeva M, Sengupta S, Sugden B: The microRNAs of Epstein-Barr Virus are expressed at dramatically differing levels among cell lines. Virology 2009, 386: 387-397. 10.1016/j.virol.2009.01.006
Cheung ST, Huang DP, Hui AB, Lo KW, Ko CW, Tsang YS, Wong N, Whitney BM, Lee JC: Nasopharyngeal carcinoma cell line (C666-1) consistently harbouring Epstein-Barr virus. Int J Cancer 1999, 83: 121-126. 10.1002/(SICI)1097-0215(19990924)83:1<121::AID-IJC21>3.0.CO;2-F
Feederle R, Linnstaedt SD, Bannert H, Lips H, Bencun M, Cullen BR, Delecluse HJ: A viral microRNA cluster strongly potentiates the transforming properties of a human herpesvirus. PLoS Pathog 2011, 7: e1001294. 10.1371/journal.ppat.1001294
Xing L, Kieff E: Epstein-Barr virus BHRF1 micro- and stable RNAs during latency III and after induction of replication. J Virol 2007, 81: 9967-9975. 10.1128/JVI.02244-06
Bankier AT, Deininger PL, Farrell PJ, Barrell BG: Sequence analysis of the 17,166 base-pair EcoRI fragment C of B95-8 Epstein-Barr virus. Mol Biol Med 1983, 1: 21-45.
Raab-Traub N, Dambaugh T, Kieff E: DNA of Epstein-Barr virus VIII: B95-8, the previous prototype, is an unusual deletion derivative. Cell 1980, 22: 257-267. 10.1016/0092-8674(80)90173-7
Fernandez AF, Rosales C, Lopez-Nieva P, Grana O, Ballestar E, Ropero S, Espada J, Melo SA, Lujambio A, Fraga MF, et al.: The dynamic DNA methylomes of double-stranded DNA viruses associated with human cancer. Genome Res 2009, 19: 438-451.
Sengupta S, den Boon JA, Chen IH, Newton MA, Dahl DB, Chen M, Cheng YJ, Westra WH, Chen CJ, Hildesheim A, et al.: Genome-wide expression profiling reveals EBV-associated inhibition of MHC class I expression in nasopharyngeal carcinoma. Cancer Res 2006, 66: 7999-8006. 10.1158/0008-5472.CAN-05-4399
Cosmopoulos K, Pegtel M, Hawkins J, Moffett H, Novina C, Middeldorp J, Thorley-Lawson DA: Comprehensive profiling of Epstein-Barr virus microRNAs in nasopharyngeal carcinoma. J Virol 2009, 83: 2357-2367. 10.1128/JVI.02104-08
Zhu JY, Pfuhl T, Motsch N, Barth S, Nicholls J, Grasser F, Meister G: Identification of novel Epstein-Barr virus microRNA genes from nasopharyngeal carcinomas. J Virol 2009, 83: 3333-3341. 10.1128/JVI.01689-08
Sotiropoulou G, Pampalakis G, Lianidou E, Mourelatos Z: Emerging roles of microRNAs as molecular switches in the integrated circuit of the cancer cell. RNA 2009, 15: 1443-1461. 10.1261/rna.1534709
Gupta A, Gartner JJ, Sethupathy P, Hatzigeorgiou AG, Fraser NW: Anti-apoptotic function of a microRNA encoded by the HSV-1 latency-associated transcript. Nature 2006, 442: 82-85.
Bennasser Y, Le SY, Yeung ML, Jeang KT: HIV-1 encoded candidate micro-RNAs and their cellular targets. Retrovirology 2004, 1: 43. 10.1186/1742-4690-1-43
Huang J, Wang F, Argyris E, Chen K, Liang Z, Tian H, Huang W, Squires K, Verlinghieri G, Zhang H: Cellular microRNAs contribute to HIV-1 latency in resting primary CD4+ T lymphocytes. Nat Med 2007, 13: 1241-1247. 10.1038/nm1639
Gitlin L, Andino R: Nucleic acid-based immune system: the antiviral potential of mammalian RNA silencing. J Virol 2003, 77: 7159-7165. 10.1128/JVI.77.13.7159-7165.2003
Iizasa H, Wulff BE, Alla NR, Maragkakis M, Megraw M, Hatzigeorgiou A, Iwakiri D, Takada K, Wiedmer A, Showe L, et al.: Editing of Epstein-Barr virus-encoded BART6 microRNAs controls their dicer targeting and consequently affects viral latency. J Biol Chem 2010, 285: 33358-33370. 10.1074/jbc.M110.138362
Cullen BR: Viruses and microRNAs. Nat Genet 2006,38(Suppl):S25-S30.
Gottwein E, Cullen BR: Viral and cellular microRNAs as determinants of viral pathogenesis and immunity. Cell Host Microbe 2008, 3: 375-387. 10.1016/j.chom.2008.05.002
Fresen KO, Hausen H: Establishment of EBNA-expressing cell lines by infection of Epstein-Barr virus (EBV)-genome-negative human lymphoma cells with different EBV strains. Int J Cancer 1976, 17: 161-6. 10.1002/ijc.2910170203
Takada K, Horinouchi K, Ono Y, Aya T, Osato T, Takahashi M, Hayasaka S: An Epstein-Barr virus-producer line Akata: establishment of the cell line and analysis of viral DNA. Virus Genes 1991, 5: 147-56. 10.1007/BF00571929
Pulvertaft JV: A study of malignant tumours in nigeria by short-term tissue culture. J Clin Pathol 1965, 18: 261-73. 10.1136/jcp.18.3.261
Klein G, Dombos L, Gothoskar B: Sensitivity of Epstein-Barr virus (EBV)producer and non-producer human lymphoblastoid cell lines to superinfection with EB-virus. Int J Cancer 1972, 10: 44-57. 10.1002/ijc.2910100108
Miller G, Robinson J, Heston L, Lipman M: Differences between laboratory strains of Epstein-Barr virus based on immortalization, abortive infection and interference. IARC Sci Publ 1975,11(1):395-408.
Wang HY, Luo M, Tereshchenko IV, Frikker DM, Cui X, Li JY, Hu G, Chu Y, Azaro MA, Lin Y, et al.: A genotyping system capable of simultaneously analyzing >1000 single nucleotide polymorphisms in a haploid genome. Genome Res 2005, 15: 276-283. 10.1101/gr.2885205
Wang H, Ach RA, Curry B: Direct and sensitive miRNA profiling from low-input total RNA. RNA 2007, 13: 151-159.
Chen HC, Chen GH, Chen YH, Liao WL, Liu CY, Chang KP, Chang YS, Chen SJ: MicroRNA deregulation and pathway alterations in nasopharyngeal carcinoma. Br J Cancer 2009, 100: 1002-1011. 10.1038/sj.bjc.6604948
We thank Prof. Musheng Zeng (Cancer Center, Sun Yat-Sen University) for his technical assistance and the gift of the cell lines. Financial supports for this study were provided by the Chinese National Natural Science Foundation Project (30770108, 31170151, 81071822), the Guangdong Provincial Key Program of Science and Technology (2009A030331005), the Natural Science Foundation of Guangdong (06021203), the Guangzhou Key Program of Science and Technology (2012Y2-00026) and National Key Basic Research Program of China (2012CB519003). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
HBJ and QCN contributed to the conception of the idea, design, data collection, drafting and writing of the manuscript. YHJ, HTJ and YCF contributed to conception of the idea, design and writing of the manuscript. YCF and YGD contributed to data collection and culture of the cells. CQQ and HJL contributed to the laboratory work and drafting of the manuscript. LRY and PLX contributed to data analysis and writing of the manuscript. HHB and ZZY contributed to conception of the idea, data collection and writing of the manuscript. LN contributed to conception and writing of the manuscript. All read and approved the manuscript.
Hong-Jie Yang, Tie-Jun Huang, Chang-Fu Yang contributed equally to this work.