Identification of Novel Epstein-Barr Virus MicroRNA Genes from Nasopharyngeal Carcinomas

doi:10.1128/JVI.01689-08

J Virol. 2009 April; 83(7): 3333–3341.

Published online 2009 January 14. doi: 10.1128/JVI.01689-08

PMCID: PMC2655542

Identification of Novel Epstein-Barr Virus MicroRNA Genes from Nasopharyngeal Carcinomas ^†

Jia Yun Zhu,¹ Thorsten Pfuhl,² Natalie Motsch,² Stephanie Barth,² John Nicholls,³ Friedrich Grässer,²^* and Gunter Meister¹^*

Center for Integrated Protein Sciences Munich (CIPSM), Laboratory of RNA Biology, Max Planck Institute of Biochemistry, Am Klopferspitz 18, 82152 Martinsried, Germany,¹ Institute of Microbiology and Hygiene, Department of Virology, University of Homburg/Saar, House 47, 66421 Homburg/Saar, Germany,² Department of Pathology, The University of Hong Kong, Pok Fu Lam, Hong Kong³

^*Corresponding author. Mailing address for Gunter Meister: Center for Integrated Protein Sciences Munich (CIPSM), Laboratory of RNA Biology, Max Planck Institute of Biochemistry, Am Klopferspitz 18, 82152 Martinsried, Germany. Phone: 49 89 8578 3420. Fax: 49 8578 3430. E-mail: meister/at/biochem.mpg.de. Mailing address for Friedrich Grässer: Institute of Microbiology and Hygiene, Department of Virology, University of Homburg/Saar, House 47, 66421 Homburg/Saar, Germany. Phone: 4968411623983. Fax: 4968411623980. E-mail: graesser/at/uks.eu

Received August 8, 2008; Accepted January 1, 2009.

This article has been cited by other articles in PMC.

Abstract

MicroRNAs (miRNAs) represent a conserved class of small noncoding RNAs that are found in all higher eukaryotes as well as some DNA viruses. miRNAs are 20 to 25 nucleotides in length and have important regulatory functions in biological processes such as embryonic development, cell differentiation, hormone secretion, and metabolism. Furthermore, miRNAs have been implicated in the pathology of various diseases, including cancer. miRNA expression profiles not only classify different types of cancer but also may even help to characterize distinct tumor stages, therefore constituting a valuable tool for prognosis. Here we report the miRNA profile of Epstein-Barr virus (EBV)-positive nasopharyngeal carcinoma (NPC) tissue samples characterized by cloning and sequencing. We found that all EBV miRNAs from the BART region are expressed in NPC tissues, whereas EBV miRNAs from the BHRF1 region are not found. Moreover, we identified two novel EBV miRNA genes originating from the BART region that have not been found in other tissues or cell lines before. We also identified three new human miRNAs which might be specific for nasopharyngeal tissues. We further show that a number of different cellular miRNAs, including miR-15a and miR-16, are up- or downregulated in NPC tissues compared to control tissues. We found that the tumor suppressor BRCA-1 is a target of miR-15a as well as miR-16, suggesting a miRNA role in NPC pathogenesis.

MicroRNAs (miRNAs) constitute a conserved class of endogenously expressed small noncoding RNAs of 20 to 25 nucleotides (nt) in size with regulatory functions in various cellular processes (1, 4, 7, 14). miRNA genes are transcribed by RNA polymerase II or III as primary transcripts that are processed to stem-loop structured precursors (pre-miRNAs) by the nuclear microprocessor complex containing the RNase III Drosha and its cofactor DGCR8 (33, 51). Pre-miRNAs are transported to the cytoplasm by the export receptor exportin-5, where the RNase III Dicer cleaves off the loop of the hairpin, thereby creating a short double-stranded RNA (1, 4, 39). Such intermediates are subsequently unwound, and one strand is incorporated as mature miRNA into a miRNA-protein complex often referred to as miRNP, while the other strand, referred to as miRNA*, is rapidly degraded (22, 39). The miRNA guides miRNPs to partially complementary sequences located in the 3′-untranslated region (UTR) of target mRNAs and regulates their expression by translational inhibition and/or mRNA destabilization (20, 23, 38, 40, 43).

Soon after the discovery of mammalian miRNAs, viruses were analyzed for miRNA expression as well (17). Indeed, DNA viruses of the herpesvirus family, including Epstein-Barr virus (EBV) (42), Kaposi's sarcoma-associated herpesvirus (9-11, 41), and human cytomegalovirus (41), encode and express miRNAs that are easily detectable in infected cells. The individual functions of the majority of the viral miRNAs are still unknown. However, some viral miRNAs have been characterized functionally in more detail. For example, human cytomegalovirus miR-UL112 targets genes of the host immune system, leading to reduced killing of natural killer cells (44). Moreover, simian virus 40 expresses a miRNA that regulates viral gene expression to reduce susceptibility to cytotoxic T cells (45). It has also been reported that Kaposi's sarcoma-associated herpesvirus miR-K12-11 functions as an ortholog of cellular miR-155 and may therefore exploit preexisting pathways in B cells (24). Very recently, it was shown that herpes simplex virus type 1 expresses miRNAs that regulate and maintain latent infection of neurons of the sensory ganglia (47).

miRNAs have been implicated in a variety of diseases, including cancer. Initially, it was reported that miR-15a and miR-16 expression is frequently downregulated or even deleted in B-cell chronic lymphocytic leukemia (13). To date, a variety of different types of cancer have been analyzed for miRNA expression (12, 19). In many cases, distinct miRNA profiles have been reported. Moreover, oncogenes such as RAS and BCL2 are controlled by the miRNA pathway, and alterations in miRNA expression result in upregulation of these oncogenes (15, 31). miRNAs can therefore function as tumor suppressors. Reciprocally, miRNAs have also been reported to function as oncogenes (12, 19, 48).

Nasopharyngeal carcinoma (NPC) is a human malignancy derived from the epithelium of the nasopharyngeal recess which shows a particularly high incidence in the area of Southeast Asia and is constantly associated with EBV infection. The full-length viral genome is contained in the nuclei of all malignant NPC cells, and its encoded viral RNAs and proteins probably contribute to the malignant phenotype (8, 50). So far, up to 23 EBV miRNAs from two distinct clusters (BHRF1 and BART) have been identified in EBV-positive cell lines and show different expression patterns depending on various infectious samples and stages (11, 26, 42).

Here we report a small RNA cloning and deep sequencing approach with different NPC samples as well as healthy tissues. In addition to three new human miRNA genes, we identify two new EBV miRNA genes from the BART region that have not been found in other tissues or cell lines before. Moreover, we demonstrate that miR-15a and miR-16 are overexpressed in NPC samples and inhibit expression of the tumor suppressor BRCA-1.

MATERIALS AND METHODS

Small RNA cloning.

The small RNA fraction from tissue samples was isolated using a mirVana-microRNA isolation kit (Ambion) and separated in a denaturing 12.5% polyacrylamide (PAA) gel stained with SYBR green II. After passive elution, RNAs with a length of 15 to 30 bases were concentrated by ethanol precipitation and dissolved in water. Next, RNAs were poly(A) tailed using poly(A) polymerase, and an adapter with the following sequence was ligated to the 5′ phosphate of the miRNAs: 5′-end adapter (43 nt), 5′-GCCTCCCTCGCGCCATCAGCTNNNNGACCTTGGCTGTCACTCA-3′. NNNN represents a “bar code” sequence for the individual samples, as follows: CAAT for NPC1, ATCG for Control-1, ACTA for NPC-2, and AGGT for Control-2. Next, first-strand cDNA synthesis was performed using an oligo(dT) linker primer [3′-end oligo(dT) linker primer (61 nucleotides), 5′-GCCTTGCCAGCCCGCTCAGACGAGACATCGCCCCGC(T)₂₅-3′] and Moloney murine leukemia virus RNase H⁻ reverse transcriptase. The resulting cDNAs were PCR amplified in 22 cycles, using the high-fidelity Phusion polymerase (Finnzymes). Amplification products (120 to 135 bp) were confirmed by polyacrylamide gel electrophoresis (PAGE) analysis. All cDNA pools were mixed in equal amounts and subjected to gel fractionation. The 120- to 135-bp fraction was electroeluted from 6% PAA gels. After isolation with Nucleospin extract II (Macherey and Nagel), cDNA pools were dissolved in 5 mM Tris-HCl, pH 8.5, at a concentration of 10 ng/μl and used in single-molecule sequencing. Massively parallel sequencing was performed by 454 Life Sciences, using a Genome Sequencer 20 system.

Oligonucleotides and plasmids.

The dual-luciferase reporter plasmid pMIR-RL was modified from the firefly luciferase reporter plasmid pMIR-REPORT (Ambion) (6, 28). The 3′ UTR of BRCA-1 mRNA was cloned via PCR amplification from cDNA libraries and ligated into the corresponding SacI and NaeI sites of pMIR-RL.

2′-O-Methylated miRNA (2′OMe-miRNA) inhibitors were designed as antisense oligonucleotides to the mature miRNAs according to sequences in the miRNA registry (microrna.sanger.ac.uk/sequences/index.shtml).

Transfections and luciferase assays.

Plasmid transfections for luciferase assays were performed with 100 ng pMIR-RL and 20 pmol 2′OMe-miRNA inhibitor per 1 × 10⁴ cells in a 96-well plate, using Lipofectamine 2000 reagent (Invitrogen) as described by the manufacturer.

Luciferase activity was measured 48 h after transfection, using a dual-luciferase reporter system as described by the manufacturer (Promega).

Sequence analysis.

Base calling and quality trimming of sequence chromatograms were done by the publicly available Phred software (21). The adapter sequences and poly(A) tails were removed from the sequences. The remaining sequences of 15 bases or longer were subjected to later analysis. Due to the disturbing signal from poly(A) tails in 454 sequencing, the sequences were analyzed semimanually, distinguishing the false “A” signals in the sequences. The set of sequences was first compared with the mature sequences of known miRNAs from miRBase v11.0 (25). The sequences that were not perfectly matched may have been artifacts of the cloning procedure or a result of nontemplated modification of mature miRNAs (3) and were annotated according to the best BLAST hit to the database. The rest of the sequences were again compared with the hairpin sequences of known miRNAs from miRBase v11.0, and the star sequences of known miRNAs that have not been annotated in the database were identified. The sequences were then subjected to a BLAST search against noncoding RNA sequences, including rRNAs, tRNAs, snoRNAs, etc., retrieved from GenBank (http://www.ncbi.nlm.nih.gov), and the annotated sequences were discarded from the sequence set. The noncoding RNA sequences from EBV were identified by a BLAST search against the EBV genome retrieved from GenBank.

The rest of the sequences, including the unannotated sequences from the EBV and human genomes, were then subjected to a search for putative novel miRNA sequences. Genomic regions containing inserts with 100-nt flanks were retrieved from Ensembl to calculate RNA secondary structures by mfold (53). Only regions that folded into hairpins and contained an insert in one of the hairpin arms were considered putative novel miRNA sequences.

Cell culture.

BL41, BL41/B95.8, Jijoye, and EREB2.5 cells were grown in RPMI 1640 medium with l-glutamine supplemented with 10% heat-inactivated fetal bovine serum, 100 units/ml penicillin, and 100 μg/ml streptomycin at 37°C in an atmosphere containing 5% CO₂. C666.1 cells were grown in the above-mentioned RPMI medium supplemented with 25 mM HEPES (Sigma) in flask collagen I (Biocoast BD) and digested by Acutase (PAA). HEK 293 cells were grown in Dulbecco's modified Eagle's medium with l-glutamine. MCF7 cells were grown in Dulbecco's modified Eagle's medium with l-glutamine supplemented with 18 μg/ml bovine insulin (Sigma).

Northern blotting.

Total RNAs were isolated from BL41, BL41/B95.8, Jijoye, EREB2.5, and C666.1 cells by using Trizol reagent (Invitrogen) according to the manufacturer's manual and were separated by 15% denaturing RNA PAGE. Northern blotting against miRNAs was performed as described before (35). The probe sequences against new miRNA candidates were as follows: ebv-mir-BART-21-5P, 5′-GTTAGTTGCCTTCACTAGTGA; ebv-mir-BART-21-3P, 5′-AAACACCAGTGGGCACAACTAG; and ebv-mir-BART-22, 5′-ACTACTAGACCATGACTTTGTA.

Ago complex immunoprecipitation and RNA extraction.

Ago complex immunoprecipitation was performed as described previously (6). In short, HEK 293 cells were grown in 10-cm plates and transfected with 1 nmol 2′OMe-miRNA inhibitor, using RNAiMAX reagent (Invitrogen) as described by the manufacturer. Forty-eight hours later, cells were lysed and centrifuged. For immunoprecipitation, 2 ml of hybridoma supernatant from monoclonal anti-Ago1-4B8 and anti-Ago2-5D4 antibodies (6) was mixed and coupled to approximately 50 μl protein G-Sepharose (GE Healthcare). Beads were subsequently incubated with HEK 293 cell lysate for 2.5 h with constant rotation at 4°C. After incubation, the beads were washed five times with RIPA-300 buffer. Finally, the beads were washed once with phosphate-buffered saline. Coprecipitated RNA was extracted using 1 volume of phenol and subsequently precipitated from the aqueous phase using 3 volumes of ethanol. The RNA pellet was treated with DNase I (Fermentas) and used for reverse transcription.

RT and qRT-PCR.

DNase I-treated RNA was reverse transcribed using a first-strand cDNA synthesis kit (Fermentas) with (dT)₁₈ oligonucleotides. Quantitative real-time reverse transcription-PCR (qRT-PCR) was performed using MESA Green qPCR Mastermix Plus for SYBR green assay with fluorescein (Eurogentec) and carried out on a MyiQ single-color real-time PCR detection system (Bio-Rad) with supplied software. The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mRNA level was used for normalization. The primers used were as follows: BRCA-1 for, 5′TAAGCCAGAATCCAGAAGGC; BRCA-1 rev, 5′GGGATGACCTTTCCACTCCT; GAPDH for, 5′AATGGAAATCCCATCACCATCT; and GAPDH rev, 5′CACCCCACTTGATTTTGG.

Western blotting.

MCF7 cells were grown in 6-cm plates and transfected with 0.5 nmol 2′OMe-miRNA inhibitor, using RNAiMAX reagent (Invitrogen). Forty-eight hours later, cells were lysed with RIPA-150 buffer. Extracted proteins were separated by 10% sodium dodecyl sulfate-PAGE and transferred to a nitrocellulose membrane (Amersham Biosciences) by semidry electroblotting. The primary antibodies used were anti-BRCA-1 (Calbiochem) and anti-β-actin (Abcam), and the secondary antibodies used were peroxidase-conjugated anti-mouse antibodies (Sigma).

Immunohistochemistry.

Immunohistochemistry with BRCA-1 antibody was performed as follows. Eight cases of nonkeratinizing, undifferentiated NPC (according to WHO criteria) were retrieved from the files of the Department of Anatomical Pathology, Queen Mary Hospital. Five-micrometer sections were cut and then microwaved in 10 mM citrate buffer, pH 6.0, for 15 min. They were blocked with 3% H₂O₂ in Tris-buffered saline for 10 min, washed, blocked with an avidin-biotin blocking system (Vector Laboratories) for 15 min each with reagent A and reagent B, washed, blocked with 10% normal donkey serum for 10 min at room temperature, and then incubated with primary antibody (anti-BRCA-1; Zytomed) at 4°C overnight. The following day, they were incubated with secondary antibody at 1/100 for 30 min at room temperature, incubated with strep-ABC complex (Vector Laboratories) at 1/100 for 30 min at room temperature, and then developed with an AEC substrate kit (Vector Laboratories) at room temperature for 30 min, followed by counterstaining with Mayer's hematoxylin and mounting with Sigma Crystal Mount reagent.

RESULTS

Small RNA cloning and sequencing from NPC and control samples.

In order to investigate a possible contribution of cellular or EBV-derived miRNAs to NPC pathogenesis, we cloned and sequenced small RNAs from EBV-positive NPC samples as well as control tissues (GEO number GSE14738). Biopsies were taken from the nasopharynx in an area where there was clinical evidence of tumor, and control tissues were taken from clinically normal mucosa from the opposite side of the nasopharynx in the same patient. We obtained 25,223 sequence reads for NPC sample 1 (NPC-1) and 22,636 reads for NPC sample 2 (NPC-2) that were further grouped according to their origins (Table (Table1).1). For control tissues, we analyzed 16,249 (Control-1) and 19,052 (Control-2) sequences. Interestingly, only for the NPC tissues were we able to map >50% of the sequences to the human or the EBV genome. The libraries obtained from control tissues contained only a minor portion of human- or EBV-derived sequences. The majority of the sequences were derived from bacterial genomes, probably due to the bacterial flora in the throat of the patients or to contamination of the samples. About 20,000 sequences of NPC-1 and about 6,000 sequences of NPC-2 were identified as known cellular or viral miRNAs. Notably, the libraries from NPC-2 and Control-2 contained significant numbers of other noncoding RNAs, which may indicate sample-specific lower RNA quality (Table (Table11).

TABLE 1.

Small RNA read numbers obtained from NPC libraries and control samples

Identification of ebv-miR-BART21 and ebv-miR-BART22.

Next, we analyzed our data sets for novel miRNA candidates. First, we investigated sequences that match the EBV genome but have not yet been annotated functionally. We found three different small RNAs derived from two individual precursors that fold into hairpins and originate from the EBV BART region (Table (Table2;2; Fig. Fig.1A).1A). We refer to these miRNA candidates as ebv-miR-BART21 and ebv-miR-BART22. Interestingly, miR-BART21 is conserved in rhesus lymphocryptovirus, a common model system for EBV infection (Table (Table2).2). We further validated the expression levels of the novel EBV miRNA candidates by Northern blotting with the EBV-positive cell lines Jijoye, EREB2.5, C666.1, and Bl41/B95.8 as well as the EBV-negative B-cell line BL-41 (Fig. 1B and C). Probes specific to the two arms of miR-BART21 as well as miR-BART22 were readily detectable in Jijoye, EREB2.5, and C666.1 cells, whereas no signal was observed in BL41/B95.8 and BL41 cells. Notably, the BL41/B95.8 cell line carries a deletion within the BART region where the two novel miRNAs are located (Fig. (Fig.1A).1A). The C666.1 cell line originated from an NPC, and its miRNA pattern closely resembles the EBV miRNA pattern in NPC tissues. It is therefore a model cell line for NPC pathogenesis (16). Furthermore, we detected star sequences for miR-BART4 and miR-BART5 in NPC samples, suggesting that some viral miRNA* molecules might be functional as well (Table (Table3).3). We used Target Scan to predict putative cellular mRNA targets of the new EBV miRNAs (36). Target Scan predicts several strong mRNA target candidates with multiple target sites (see Table S1 in the supplemental material).

TABLE 2.

Novel miRNA genes identified from NPC samples

FIG. 1.

(A) Schematic locations of putative new miRNAs in the BART miRNA cluster of the EBV genome. Each box with an arrow indicates an EBV miRNA precursor and its transcription direction. The putative new miRNA precursors are depicted in red. The numbers in (more ...)

TABLE 3.

miRNAs from known novel and precursors identified from NPC samples

To further confirm the existence of the new EBV miRNAs in NPC samples, we measured the new miRNA levels in six undifferentiated NPC samples and one control tissue. HEK 293 cells served as a negative control (see Fig. S1 in the supplemental material). Since the amount of extracted RNA from the tissues was very low, we polyadenylated and reverse transcribed the RNA as described previously (29), generating cDNA from miRNA suitable for qRT-PCR validation. GAPDH mRNA was taken for normalization. The background signal from the control tissue sample was set to 1 for calculation. On that basis, all six NPC samples gave strong signals for one known EBV miRNA (ebv-BART1-5p) and the three new EBV miRNAs (ebv-BART21-5p, ebv-BART21-3p, and ebv-BART22) tested, indicating that the new EBV miRNAs can be detected in a number of NPC samples.

In order to find novel human miRNAs, we analyzed the pool of unknown sequences which match the human genome. Indeed, we identified three new miRNA genes. Two of them were recently annotated as miR-1301 and miR-1307, while this work was in progress. Using miRBase v11.0 (25), we named the third one miR-2110 (Table (Table2).2). miR-1301 is located in an intergenic region, whereas miR-1307 and miR-2110 reside in the 5′ UTRs of protein coding genes (Table (Table2).2). Additionally, so far we have found six unidentified miRNA* sequences that derive from known miRNA genes and have not been cloned before (Table (Table3).3). In summary, we have identified two new EBV miRNA genes that give rise to mature miR-BART21-5p, miR-BART21-3p, and miR-BART22. In addition, we have found three novel human miRNA genes, referred to as miR-1301, miR-1307, and miR-2110.

miRNA expression analysis suggests an NPC-specific miRNA signature.

miRNA expression analyses revealed distinct miRNA signatures for different types of tumors (12). In order to establish a miRNA expression profile characteristic of NPCs, we examined the most abundant miRNAs in the individual libraries obtained from NPC or control tissues (Fig. (Fig.2;2; see Table S2 in the supplemental material). About 5% of the NPC-1 and 19% of the NPC-2 miRNAs were identified as EBV miRNAs, indicating that viral miRNAs are highly expressed in NPCs (Fig. (Fig.2A).2A). Notably, small numbers of EBV miRNAs were present in Control-1 as well, which might be due to a minor contamination of the control tissue with tumor tissue (Table (Table1).1). Among the EBV miRNAs, we found high levels of ebv-miR-BART1, ebv-miR-BART4, ebv-miR-BART6, ebv-miR-BART7, ebv-miR-BART11, ebv-miR-BART12, and ebv-miR-BART19 as well as the new miRNAs ebv-miR-BART21 and ebv-miR-BART22. Interestingly, no EBV miRNA from the BHRF1 region was identified in our libraries, suggesting that miRNAs from this region might not be involved in NPC pathogenesis (Fig. (Fig.2B).2B). The most abundant human miRNAs are presented in Fig. Fig.2C.2C. The 20 most frequently cloned miRNAs were similar between NPC and control tissues (see Table S2 in the supplemental material). However, miR-320, miR-17-5p, and miR-652 were less frequent in the libraries from the NPCs than in those from control tissues, suggesting that these miRNAs might be downregulated in NPC tissue. Furthermore, the read numbers for miR-23a/b (both miRNAs were indistinguishable in our sequencing data sets), miR-200c, and miR-27a/b were significantly increased in the NPC tumor samples, suggesting specific upregulation of these miRNAs. Consistent with the cloning data, miR-23a/b and miR-27a/b were processed from the same primary transcript.

FIG. 2.

miRNA expression in NPC and control tissue. (A) Schematic representation of cellular and EBV-derived miRNA fractions in NPC tissue samples. (B) Relative abundances of EBV miRNAs in the NPC-1 and NPC-2 libraries. Individual read numbers are presented as (more ...)

Misregulated miRNAs miR-15a and miR-16 inhibit expression of the tumor suppressor BRCA-1.

Next, we analyzed potential targets of miRNAs that are misexpressed in the two NPC libraries analyzed above. In both NPC samples, expression of miR-15a and miR-16, which are derived from one primary transcript, was upregulated, and indeed, bioinformatic predictions identified target sites for miR-15a and miR-16 in the 3′ UTR of the tumor suppressor BRCA-1 (Fig. 3A and B). Therefore, we hypothesized that high levels of miR-15a and miR-16 may lead to low levels of BRCA-1, resulting in low tumor suppressor activity and tumor growth. In order to analyze effects of miR-15a and miR-16 on BRCA-1 expression, we fused the BRCA-1 3′ UTR to a firefly luciferase reporter. The reporter plasmid was transfected into HEK 293 cells together with 2′OMe inhibitors against endogenous miR-15a or miR-16 (Fig. (Fig.3C).3C). Indeed, luciferase expression was elevated when endogenous miR-15a or miR-16 was inhibited, suggesting that BRCA-1 is regulated by miR-15a or -16.

FIG. 3.

miR-15a and miR-16 target the BRCA-1 mRNA. (A) Elevated miR-15a and miR-16 expression in two NPC samples compared to that in healthy tissues from the same individuals. miRNA read numbers are shown as percentages of the total miRNA reads in the libraries. (more ...)

To further validate a functional interaction of the miRNA pathway with the BRCA-1 mRNA, we investigated binding of Ago proteins to the BRCA-1 mRNA (Fig. (Fig.3D).3D). We and others have shown before that miRNA target mRNAs can be immunoprecipitated and identified using antibodies against Ago proteins (6, 18, 27, 32). Endogenous Ago1 and Ago2 complexes were immunoprecipitated from a HEK 293 cell lysate by using Ago1- and Ago2-specific antibodies, and coimmunoprecipitated mRNAs were analyzed by qRT-PCR (Fig. (Fig.3D).3D). Indeed, BRCA-1 mRNA was coimmunoprecipitated with anti-Ago1 and anti-Ago2 antibodies. Moreover, the BRCA-1 mRNA was released from Ago protein complexes when 2′OMe inhibitors complementary to miR-15a or miR-16, which interrupt miRNA-mRNA interactions, were transfected, indicating that BRCA-1 is targeted by the miRNA pathway.

Finally, we measured BRCA-1 mRNA as well as protein levels in the presence or absence of miR-15a or miR-16 inhibitor (Fig. 3E and F). It has been shown before that miRNAs not only function in translation but also have effects on mRNA stability. Endogenous BRCA-1 mRNA was robustly increased when miR-15a or miR-16 was inhibited. Notably, inhibition of let-7a, a miRNA that does not target BRCA-1, led to a decrease of BRCA-1 mRNA levels, which was most likely due to unrelated, secondary let-7a effects. Consistent with the mRNA expression data, we also observed an increase in BRCA-1 protein levels when miR-15a and miR-16 were inhibited (Fig. (Fig.3F3F).

In summary, our data suggest that miR-15a and miR-16 are upregulated in NPC samples and target the tumor suppressor BRCA-1.

Reduced BRCA-1 signals in NPC tissue sections.

To further support our hypothesis that BRCA-1 expression is reduced in NPC tissues, we performed immunohistochemistry experiments using paraffin-embedded tissue sections (Fig. (Fig.4).4). We analyzed archived tissue biopsy samples from eight NPCs and tested them for the presence or absence of BRCA-1, using a monoclonal antibody which identifies the C-terminal region of the BRCA-1 protein, which according to the manufacturer's instructions is cytoplasmic in location. In eight undifferentiated NPCs, there was strong cytoplasmic staining in two cases, medium to weak staining in four cases, and no staining in two cases (Fig. (Fig.4A).4A). The adjacent stroma epithelium showed medium positivity, which was cytoplasmic (Fig. (Fig.4B).4B). Thus, in six of eight NPCs, there was a degree of reduction or absence of staining with BRCA-1. Taken together, these data further support our finding that BRCA-1 expression is decreased in a number of NPC samples.

FIG. 4.

Analysis of BRCA-1 protein levels in various NPC tissue sections by immunohistochemistry. Paraffin-embedded sections were stained with antibodies against BRCA-1. (A) Summary of BRCA-1 protein detection in eight different NPC samples. (B) Example of weak (more ...)

DISCUSSION

Using small RNA cloning and sequencing, we have characterized the miRNA expression profile of NPC tissues. EBV is reported to be present in all NPCs and can transform cells, which subsequently induces cell proliferation and tumor growth. Indeed, we find that EBV expresses all miRNAs from the BART cluster in NPC tissues, suggesting that these miRNAs may have an important function in maintaining the virus in NPC tissues. Moreover, we identified two novel and highly abundant EBV miRNA genes, namely, miR-BART21 and miR-BART22. Notably, no miRNA originating from the BHRF1 region of the EBV genome was found in our data sets, suggesting that these miRNAs might not be required for NPC pathogenesis. More provocatively, the absence of BHRF1 miRNA expression might be important for NPC pathogenesis, and only viruses with impaired BHRF1 miRNA expression might be capable of NPC tumor formation. Interestingly, BHRF1 miRNA expression was not found in EBV-positive gastric cancer samples as well (34). The detailed characterization of the small RNA profiles of other EBV-positive tumor types, such as T-cell lymphomas, B-cell lymphomas, and gastric cancer, will help to elucidate the contribution of EBV miRNAs to cancer pathogenesis.

Like other tumors, NPC samples are characterized by distinct miRNA expression profiles. In particular, we demonstrate that miR-23a/b, miR-200c, and miR-27a/b are significantly upregulated, whereas miR-320, miR-17-5p, and miR-652 are downregulated in NPC samples. Bioinformatic predictions suggest that miR-23a and -b, which are the most abundant miRNAs in NPC samples, may regulate the tumor suppressors ITGB1 (integrin-beta1), LIG4 (ATP-dependent DNA ligase IV), and SAFB (scaffold attachment factor B). It is tempting to speculate that other upregulated miRNAs might target tumor suppressors as well. Indeed, we find that the expression of miR-15a and miR-16 is elevated in NPC tissues and show that both miRNAs can target BRCA-1, a tumor suppressor gene important for breast cancer (49). Interestingly, miR-15a and miR-16 have been implicated in cancer before. In B-cell lymphocytic leukemia, miR-15a and miR-16 expression is downregulated, and it has been suggested that both miRNAs target the bcl-2 oncogene (12). Together with our data, these findings suggest that miR-15a and miR-16 can target oncogenes as well as tumor suppressors, depending on individual tissues.

BRCA-1 was shown to interact with the E6 and E7 oncoproteins of human papillomaviruses, and the association of BRCA-1 with E6 and E7 reverses the inhibition of expression of estrogen receptor alpha (52). It is thought that the growth arrest normally brought about by BRCA-1 is overcome through interaction with hypophosphorylated pRB (2). So far, no involvement of BRCA-1 in EBV-mediated induction of NPC has been reported. However, it was shown that the BRCA-1-associated ZRBK1/KAP-1 corepressor forms a tight complex with the BBLF2/3 protein of EBV (37). BBLF2/3 in turn is the bridging factor for the tripartite EBV helicase-primase complex, consisting of the helicase BBLF4, the primase BSLF1, and BBLF2/3 as the linking molecule. The ectopic expression of ZRKB1, KAP-1, or BRCA-1 stimulated oriLyt-driven replication of EBV DNA, establishing a functional link between the presence of BRCA-1 and lytic viral replication. Sera of NPC patients exhibit peculiar signs of lytic replication, while the tumor cells are conspicuously devoid of lytic viral proteins. The downregulation of BRCA-1 by miR-15a and miR-16 might inhibit the induction of lytic cycle replication in the tumor cells. This notion is in line with our previous demonstration that the EBV-encoded miRNA ebv-miR-BART2 downregulates the viral DNA polymerase BALF5, resulting in reduced virus production (5). A function for miRNAs encoded by the LAT latency transcript of herpes simplex virus in the reduction of virus production was recently demonstrated (46, 47).

The identification of cellular miRNA targets as well as the analysis of more NPC samples will help to elucidate the molecular basis of NPC pathogenesis. Our data might be a starting point for identifying and classifying NPCs, probably even at very early stages. Moreover, interfering with overrepresented miRNAs or artificial expression of underrepresented miRNAs might be useful for NPC therapy in the future.

Supplementary Material

[Supplemental material]

Click here to view.

Acknowledgments

We thank Adam Wong Chun Nin, Sabine Rottmüller, and Bernd Haas for technical assistance; Lasse Weinmann, Sabine Rüdel, Julia Höck, and Anne Frohn for comments on the manuscript; and Vertis Biotechnology (Weihenstephan, Germany) for cloning and help with data analysis.

This work was supported by the Max-Planck Society, in part by a grant from the Deutsche Krebshilfe (grant 107166 to G.M. and F.G.), and by Deutsche Forschungsgemeinschaft grant GR950/11-2 (to F.G.).

Footnotes

Published ahead of print on 14 January 2009.

^†Supplemental material for this article may be found at http://jvi.asm.org/.

REFERENCES

1. Ambros, V. 2004. The functions of animal microRNAs. Nature 431350-355. [PubMed]

2. Aprelikova, O. N., B. S. Fang, E. G. Meissner, S. Cotter, M. Campbell, A. Kuthiala, M. Bessho, R. A. Jensen, and E. T. Liu. 1999. BRCA1-associated growth arrest is RB-dependent. Proc. Natl. Acad. Sci. USA 9611866-11871. [PubMed]

3. Aravin, A., and T. Tuschl. 2005. Identification and characterization of small RNAs involved in RNA silencing. FEBS Lett. 5795830-5840. [PubMed]

4. Bartel, D. P. 2004. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116281-297. [PubMed]

5. Barth, S., T. Pfuhl, A. Mamiani, C. Ehses, K. Roemer, E. Kremmer, C. Jaker, J. Hock, G. Meister, and F. A. Grasser. 2008. Epstein-Barr virus-encoded microRNA miR-BART2 down-regulates the viral DNA polymerase BALF5. Nucleic Acids Res. 36666-675. [PMC free article] [PubMed]

6. Beitzinger, M., L. Peters, J. Y. Zhu, E. Kremmer, and G. Meister. 2007. Identification of human microRNA targets from isolated Argonaute protein complexes. RNA Biol. 476-84. [PubMed]

7. Bushati, N., and S. M. Cohen. 2007. microRNA functions. Annu. Rev. Cell Dev. Biol. 23175-205. [PubMed]

8. Busson, P., C. Keryer, T. Ooka, and M. Corbex. 2004. EBV-associated nasopharyngeal carcinomas: from epidemiology to virus-targeting strategies. Trends Microbiol. 12356-360. [PubMed]

9. Cai, X., and B. R. Cullen. 2006. Transcriptional origin of Kaposi's sarcoma-associated herpesvirus microRNAs. J. Virol. 802234-2242. [PMC free article] [PubMed]

10. Cai, X., S. Lu, Z. Zhang, C. M. Gonzalez, B. Damania, and B. R. Cullen. 2005. Kaposi's sarcoma-associated herpesvirus expresses an array of viral microRNAs in latently infected cells. Proc. Natl. Acad. Sci. USA 1025570-5575. [PubMed]

11. Cai, X., A. Schafer, S. Lu, J. P. Bilello, R. C. Desrosiers, R. Edwards, N. Raab-Traub, and B. R. Cullen. 2006. Epstein-Barr virus microRNAs are evolutionarily conserved and differentially expressed. PLoS Pathog. 2e23. [PMC free article] [PubMed]

12. Calin, G. A., and C. M. Croce. 2006. MicroRNA signatures in human cancers. Nat. Rev. Cancer 6857-866. [PubMed]

13. Calin, G. A., C. D. Dumitru, M. Shimizu, R. Bichi, S. Zupo, E. Noch, H. Aldler, S. Rattan, M. Keating, K. Rai, L. Rassenti, T. Kipps, M. Negrini, F. Bullrich, and C. M. Croce. 2002. Frequent deletions and down-regulation of microRNA genes miR-15 and miR-16 at 13q14 in chronic lymphocytic leukemia. Proc. Natl. Acad. Sci. USA 9915524-15529. [PubMed]

14. Chen, P. Y., and G. Meister. 2005. microRNA-guided posttranscriptional gene regulation. Biol. Chem. 3861205-1218. [PubMed]

15. Cimmino, A., G. A. Calin, M. Fabbri, M. V. Iorio, M. Ferracin, M. Shimizu, S. E. Wojcik, R. I. Aqeilan, S. Zupo, M. Dono, L. Rassenti, H. Alder, S. Volinia, C. G. Liu, T. J. Kipps, M. Negrini, and C. M. Croce. 2005. miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc. Natl. Acad. Sci. USA 10213944-13949. [PubMed]

16. Cosmopoulos, K., M. Pegtel, J. Hawkins, H. Moffett, C. Novina, J. Middeldorp, and D. A. Thorley-Lawson. 2009. Comprehensive profiling of Epstein-Barr virus microRNAs in nasopharyngeal carcinoma. J. Virol. 832357-2367. [PMC free article] [PubMed]

17. Cullen, B. R. 2006. Viruses and microRNAs. Nat. Genet. 38(Suppl.)S25-S30. [PubMed]

18. Easow, G., A. A. Teleman, and S. M. Cohen. 2007. Isolation of microRNA targets by miRNP immunopurification. RNA 131198-1204. [PubMed]

19. Esquela-Kerscher, A., and F. J. Slack. 2006. Oncomirs—microRNAs with a role in cancer. Nat. Rev. Cancer 6259-269. [PubMed]

20. Eulalio, A., E. Huntzinger, and E. Izaurralde. 2008. Getting to the root of miRNA-mediated gene silencing. Cell 1329-14. [PubMed]

21. Ewing, B., L. Hillier, M. C. Wendl, and P. Green. 1998. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 8175-185. [PubMed]

22. Farazi, T. A., S. A. Juranek, and T. Tuschl. 2008. The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members. Development 1351201-1214. [PubMed]

23. Filipowicz, W., L. Jaskiewicz, F. A. Kolb, and R. S. Pillai. 2005. Post-transcriptional gene silencing by siRNAs and miRNAs. Curr. Opin. Struct. Biol. 15331-341. [PubMed]

24. Gottwein, E., N. Mukherjee, C. Sachse, C. Frenzel, W. H. Majoros, J. T. Chi, R. Braich, M. Manoharan, J. Soutschek, U. Ohler, and B. R. Cullen. 2007. A viral microRNA functions as an orthologue of cellular miR-155. Nature 4501096-1099. [PMC free article] [PubMed]

25. Griffiths-Jones, S. 2004. The microRNA Registry. Nucleic Acids Res. 32D109-D111. [PMC free article] [PubMed]

26. Grundhoff, A., C. S. Sullivan, and D. Ganem. 2006. A combined computational and microarray-based approach identifies novel microRNAs encoded by human gamma-herpesviruses. RNA 12733-750. [PubMed]

27. Hendrickson, D. G., D. J. Hogan, D. Herschlag, J. E. Ferrell, and P. O. Brown. 2008. Systematic identification of mRNAs recruited to argonaute 2 by specific microRNAs and corresponding changes in transcript abundance. PLoS ONE 3e2126. [PMC free article] [PubMed]

28. Hock, J., L. Weinmann, C. Ender, S. Rudel, E. Kremmer, M. Raabe, H. Urlaub, and G. Meister. 2007. Proteomic and functional analysis of Argonaute-containing mRNA-protein complexes in human cells. EMBO Rep. 81052-1060. [PubMed]

29. Hurteau, G. J., S. D. Spivack, and G. J. Brock. 2006. Potential mRNA degradation targets of hsa-miR-200c, identified using informatics and qRT-PCR. Cell Cycle 51951-1956. [PubMed]

30. John, B., A. J. Enright, A. Aravin, T. Tuschl, C. Sander, and D. S. Marks. 2004. Human microRNA targets. PLoS Biol. 2e363. [PMC free article] [PubMed]

31. Johnson, S. M., H. Grosshans, J. Shingara, M. Byrom, R. Jarvis, A. Cheng, E. Labourier, K. L. Reinert, D. Brown, and F. J. Slack. 2005. RAS is regulated by the let-7 microRNA family. Cell 120635-647. [PubMed]

32. Karginov, F. V., C. Conaco, Z. Xuan, B. H. Schmidt, J. S. Parker, G. Mandel, and G. J. Hannon. 2007. A biochemical approach to identifying microRNA targets. Proc. Natl. Acad. Sci. USA 10419291-19296. [PubMed]

33. Kim, V. N. 2005. MicroRNA biogenesis: coordinated cropping and dicing. Nat. Rev. Mol. Cell. Biol. 6376-385. [PubMed]

34. Kim Do, N., H. S. Chae, S. T. Oh, J. H. Kang, C. H. Park, W. S. Park, K. Takada, J. M. Lee, W. K. Lee, and S. K. Lee. 2007. Expression of viral microRNAs in Epstein-Barr virus-associated gastric carcinoma. J. Virol. 811033-1036. [PMC free article] [PubMed]

35. Lagos-Quintana, M., R. Rauhut, A. Yalcin, J. Meyer, W. Lendeckel, and T. Tuschl. 2002. Identification of tissue-specific microRNAs from mouse. Curr. Biol. 12735-739. [PubMed]

36. Lewis, B. P., I. Shih, M. W. Jones-Rhoades, D. P. Bartel, and C. B. Burge. 2003. Prediction of mammalian microRNA targets. Cell 115787-798. [PubMed]

37. Liao, G., J. Huang, E. D. Fixman, and S. D. Hayward. 2005. The Epstein-Barr virus replication protein BBLF2/3 provides an origin-tethering function through interaction with the zinc finger DNA binding protein ZBRK1 and the KAP-1 corepressor. J. Virol. 79245-256. [PMC free article] [PubMed]

38. Meister, G. 2007. miRNAs get an early start on translational silencing. Cell 13125-28. [PubMed]

39. Meister, G., and T. Tuschl. 2004. Mechanisms of gene silencing by double-stranded RNA. Nature 431343-349. [PubMed]

40. Peters, L., and G. Meister. 2007. Argonaute proteins: mediators of RNA silencing. Mol. Cell 26611-623. [PubMed]

41. Pfeffer, S., A. Sewer, M. Lagos-Quintana, R. Sheridan, C. Sander, F. A. Grasser, L. F. van Dyk, C. K. Ho, S. Shuman, M. Chien, J. J. Russo, J. Ju, G. Randall, B. D. Lindenbach, C. M. Rice, V. Simon, D. D. Ho, M. Zavolan, and T. Tuschl. 2005. Identification of microRNAs of the herpesvirus family. Nat. Methods 2269-276. [PubMed]

42. Pfeffer, S., M. Zavolan, F. A. Grasser, M. Chien, J. J. Russo, J. Ju, B. John, A. J. Enright, D. Marks, C. Sander, and T. Tuschl. 2004. Identification of virus-encoded microRNAs. Science 304734-736. [PubMed]

43. Pillai, R. S., S. N. Bhattacharyya, and W. Filipowicz. 2007. Repression of protein synthesis by miRNAs: how many mechanisms? Trends Cell Biol. 17118-126. [PubMed]

44. Stern-Ginossar, N., N. Elefant, A. Zimmermann, D. G. Wolf, N. Saleh, M. Biton, E. Horwitz, Z. Prokocimer, M. Prichard, G. Hahn, D. Goldman-Wohl, C. Greenfield, S. Yagel, H. Hengel, Y. Altuvia, H. Margalit, and O. Mandelboim. 2007. Host immune system gene targeting by a viral miRNA. Science 317376-381. [PubMed]

45. Sullivan, C. S., A. T. Grundhoff, S. Tevethia, J. M. Pipas, and D. Ganem. 2005. SV40-encoded microRNAs regulate viral gene expression and reduce susceptibility to cytotoxic T cells. Nature 435682-686. [PubMed]

46. Tang, S., A. Patel, and P. R. Krause. 2009. Novel less-abundant viral microRNAs encoded by herpes simplex virus 2 latency-associated transcript and their roles in regulating ICP34.5 and ICP0 mRNAs. J. Virol. 831433-1442. [PMC free article] [PubMed]

47. Umbach, J. L., M. F. Kramer, I. Jurak, H. W. Karnowski, D. M. Coen, and B. R. Cullen. 2008. MicroRNAs expressed by herpes simplex virus 1 during latent infection regulate viral mRNAs. Nature 454780-783. [PMC free article] [PubMed]

48. Voorhoeve, P. M., C. le Sage, M. Schrier, A. J. Gillis, H. Stoop, R. Nagel, Y. P. Liu, J. van Duijse, J. Drost, A. Griekspoor, E. Zlotorynski, N. Yabuta, G. De Vita, H. Nojima, L. H. Looijenga, and R. Agami. 2006. A genetic screen implicates miRNA-372 and miRNA-373 as oncogenes in testicular germ cell tumors. Cell 1241169-1181. [PubMed]

49. Wang, W. 2007. Emergence of a DNA-damage response network consisting of Fanconi anaemia and BRCA proteins. Nat. Rev. Genet. 8735-748. [PubMed]

50. Wei, W. I., and J. S. Sham. 2005. Nasopharyngeal carcinoma. Lancet 3652041-2054. [PubMed]

51. Zamore, P. D., and B. Haley. 2005. Ribo-gnome: the big world of small RNAs. Science 3091519-1524. [PubMed]

52. Zhang, Y., S. Fan, Q. Meng, Y. Ma, P. Katiyar, R. Schlegel, and E. M. Rosen. 2005. BRCA1 interaction with human papillomavirus oncoproteins. J. Biol. Chem. 28033165-33177. [PubMed]

53. Zuker, M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 313406-3415. [PMC free article] [PubMed]

Articles from Journal of Virology are provided here courtesy of
American Society for Microbiology (ASM)