Epstein-Barr virus (EBV) is a tumor virus that establishes lifelong infection in most of humanity, despite eliciting strong and stable virus-specific immune responses. EBV encodes at least 44 miRNAs, most of them with unknown function. Here, we show that multiple EBV miRNAs modulate immune recognition of recently infected primary B cells, EBV's natural target cells. EBV miRNAs collectively and specifically suppress release of proinflammatory cytokines such as IL-12, repress differentiation of naive CD4+ T cells to Th1 cells, interfere with peptide processing and presentation on HLA class II, and thus reduce activation of cytotoxic EBV-specific CD4+ effector T cells and killing of infected B cells. Our findings identify a previously unknown viral strategy of immune evasion. By rapidly expressing multiple miRNAs, which are themselves nonimmunogenic, EBV counteracts recognition by CD4+ T cells and establishes a program of reduced immunogenicity in recently infected B cells, allowing the virus to express viral proteins required for establishment of life-long infection.
INTRODUCTION
EBV is both ubiquitous and immunogenic. This oncogenic herpesvirus (IARC Working Group on the Evaluation of Carcinogenic Risks to Humans, 2010) has evolved multiple genes to fend off immune responses when its infection is established (Hislop et al., 2002; Rowe et al., 2007; Ressing et al., 2008; Zuo et al., 2009; Qiu et al., 2011; Rancan et al., 2015). Despite these measures, EBV-specific T cells constitute a considerable fraction of the memory T cell repertoire of the latently infected human host (Hislop et al., 2002) and are essential in controlling latent EBV infection (Moosmann et al., 2010). In fact, immunocompromised patients have an increased incidence of EBV-associated malignancies (Gottschalk et al., 2005).
EBV infects nondividing B lymphocytes, activates them, and drives them to proliferate, thus amplifying the load of viral genomes. Once activated, infected B cells acquire properties of antigen-presenting cells. After infection, they rapidly present epitopes of structural proteins from incoming virus particles and transiently express lytic genes that are otherwise characteristic of EBV's productive cycle (Kalla and Hammerschmidt, 2012). This prelatent phase of infection includes expression of two genes coding for viral immunoevasins, BNLF2a and BCRF1 (Jochum et al., 2012), which inhibit the recognition of the infected cells by EBV-specific effector T cells and natural killer cells, respectively. These two viral proteins are insufficient, however, to overcome T cell recognition (Jochum et al., 2012). Within 7–10 d, EBV establishes a latent infection in the infected B cells and expresses only few or no viral genes, which reduces their risk of becoming eliminated by the immune-competent host.
Thus, early infection could be EBV’s Achilles heel, a window when the infected cell expresses and presents many viral antigens to immune cells but is inadequately protected from the host's immune response. We have now established that EBV's miRNAs overcome this vulnerability; they protect newly infected B lymphocytes from immune eradication by CD4+ T cells, supporting EBV’s lifelong success.
EBV encodes at least 44 microRNAs (miRNAs; Barth et al., 2011), which are small RNA regulatory molecules of ∼22 nt in length (Bartel, 2004). miRNAs encoded by herpesviruses are reported to play important roles in cell proliferation, development, immune regulation, and apoptosis in infected cells (Skalsky and Cullen, 2010). The EBV-encoded miRNAs have been found to control expression of several cellular genes with antiapoptotic functions, but they also reportedly down-regulate MICB (Nachmani et al., 2009), CXCL11 (Xia et al., 2008), and NLRP3 (Haneklaus et al., 2012) and thus interfere with innate immune responses and inflammation. Interestingly, MICB, a gene encoding a ligand for the activating receptor NKG2D expressed on T and NK cells, is also targeted by miRNAs of Kaposi sarcoma–associated herpesvirus and human cytomegalovirus (Nachmani et al., 2009; Grundhoff and Sullivan, 2011). These studies imply that certain miRNAs encoded by herpesviruses target pathways involved in innate immune recognition.
EBV’s miRNAs have been studied by several groups with established, EBV-infected cell lines obtained from biopsies of nasopharyngeal carcinoma and Burkitt’s lymphoma, or lymphoblastoid cell lines (LCLs) derived from infecting primary B lymphocytes with EBV in vitro (Dölken et al., 2010; Gottwein et al., 2011; Kuzembayeva et al., 2012; Riley et al., 2012; Erhard et al., 2013). High-throughput target screens using immunoprecipitation of the RNA-induced silencing complex (RISC) and deep sequencing have identified many potential targets of EBV miRNAs, but the catalogs of predicted targets assembled by different groups have a surprisingly small overlap (Klinke et al., 2014). This lack of consensus may be due to the accumulation of profound differences in gene expression between different long-term, cultivated EBV-infected cell lines that do not reflect the impact of EBV’s miRNAs in vivo.
To circumvent these problems, we developed an experimental approach using primary human B lymphocytes, and analyzed them during their initial days of EBV infection (Seto et al., 2010; Vereide et al., 2014). We infected the B lymphocytes with two EBV strains with and without miRNA genes, compared the gene expression in the infected cells, and examined them for their immune recognition.
We found that EBV-encoded miRNAs regulated several immune pathways, which affected CD4+ T cell differentiation and activation. In addition, key molecules important for interactions with CD4+ T cells were down-regulated. EBV miRNAs repressed the secretion of IL-12, which resulted in suppression of type 1 helper T cell (Th1) differentiation. Viral miRNAs controlled gene expression of HLA class II and three lysosomal enzymes important for proteolysis and epitope presentation to CD4+ T cells. Such a wholesale inhibition of adaptive immune responses by multiple miRNAs of a single pathogen is unprecedented. Our findings explain the abundance of miRNAs in complex persisting viruses, and clarify how EBV can escape elimination for the lifetime of its host in spite of intense adaptive immune responses.
RESULTS
EBV miRNAs control immune regulatory pathways
We searched for cellular targets of EBV's miRNAs, using an experimental system that closely mimics human infection in vivo. Two strains of EBV, a laboratory strain (wt/B95-8) that expresses 13 miRNAs, and its derivative (ΔmiR) that expresses none (Seto et al., 2010) were used to infect freshly isolated B lymphocytes from six donors. We used carefully titrated virus stocks and infected the cells with optimal doses of both viruses (Steinbrück et al., 2015). No differences in the percentage of infected cells were seen when comparing cells infected with wt/B95-8 versus ΔmiR EBV. RNAs were isolated on day 5 after infection and sequenced (available from GEO under accession no. GSE75776; see Materials and methods). Genes that were differentially expressed in cells infected with wt/B95-8 versus ΔmiR EBV were identified with those having an absolute z-score >1.6 (Fig. 1 A and Table S1). These genes included the published miRNA targets LY75/DEC205 (Skalsky et al., 2012) and IPO7 (Dölken et al., 2010). Genes that were consistently down-regulated in wt/B95-8 EBV-infected cells were grouped according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway categories (Fig. 1 B). Down-regulated genes were predominant in pathways linked to apoptosis, cell cycle regulation, and p53 signaling, which were previously proposed to be regulated by EBV miRNAs (Seto et al., 2010; Feederle et al., 2011a,b; Vereide et al., 2014). Unexpectedly, EBV's miRNAs also regulated a wide array of genes with functions in immunity, such as cytokine–cytokine receptor interactions, antigen processing, and HLAs and co-stimulatory molecules (Fig. 1, B and C; and Table S1). We immunoprecipitated RISC (RISC-IP) and found that 14.5% (±2.4% SD) of all miRNAs were of viral origin in wt/B95-8 EBV-infected cells, dominated by miRNAs of the BHRF1 gene cluster (Fig. 1 D). No appreciable viral miRNA reads were found in cells infected with ΔmiR EBV (Fig. 1 D), suggesting that the B lymphocytes of six donors were free of EBV field strains. In wt/B95-8 EBV-infected cells, we detected viral miRNAs as early as day 1 after infection, which reached high levels 5 days post infection (dpi; Fig. 1 E). In RISC-IP, detection of mRNAs was variable among infected B cells of the different donors, a phenomenon that was reported earlier using a related model of established infection and PAR-CLIP experiments (Skalsky et al., 2012; GEO accession no. GSE41437). Therefore, we focused our analyses on candidate mRNAs that were uniformly regulated in all samples (Fig. 1 C), and used RISC-IP results to confirm them (Table S1).
EBV miRNAs inhibit secretion of proinflammatory cytokines and antigen presentation
We confirmed that EBV's miRNAs regulate cytokines central to immune function. Supernatants from B cells infected with the two strains of EBV were assayed for the levels of IL-6, IL-10, TNF, IL12B (IL-12p40), IL-12 (p35/p40), and IL-23 (p19/p40). We added CpG DNA, which stimulates TLR9, for the detection of IL-6 secreted from EBV-infected cells (Iskra et al., 2010). B cells infected with wt/B95-8 EBV secreted less IL-6, TNF, and IL-12p40 than B cells infected with ΔmiR EBV. In contrast, release of the anti-inflammatory cytokine IL-10 appeared to be unaffected by viral miRNAs (Fig. 2 A) consistent with our transcriptome analysis (Fig. 1 C). Secretion of IL-12 (p35/p40 or IL-12p70) and IL-23 (p19/p40), both of which contain the IL-12p40 subunit (Szabo et al., 2003), encoded by the IL12B gene, was significantly reduced in wt/B95-8 EBV-infected cells compared with ΔmiR EBV-infected cells (Fig. 2 A). Viral miRNAs also inhibited the secretion of IL-12p40 from PBMCs infected with wt/B95-8 EBV (Fig. 2 B). IL-12p40 secretion from PBMCs infected with ΔmiR EBV was reduced when B cells were removed from the PBMCs, indicating that B cells are the main contributors to release of IL-12p40 in PBMCs. Remarkably, our transcriptome analysis revealed the consistent reduction of IL12B mRNA with EBV's miRNAs, reducing it by 80% in all six donors' B lymphocytes (Fig. 2 C). Quantitative RT-PCR confirmed this finding (Fig. 2 D).
Multiple EBV miRNAs target IL12B and prevent Th1 differentiation of naive CD4+ T cells
We investigated whether IL12B was a direct target of EBV miRNAs. EBV's miR-BART1, miR-BART2, and miR-BHRF1-2 repressed the luciferase activity of the IL12B reporter (Fig. 3 A). The mutation of predicted binding sites of miR-BART1, miR-BART2, or miR-BHRF1-2 abrogated their ability to inhibit the IL12B reporter (Fig. 3 A and Fig. S1), confirming the direct control of IL12B by these miRNAs. We similarly analyzed miR-BART10 and miR-BART22, which are present in field strains of EBV but not in wt/B95-8 EBV. For these miRNAs, mutations of their predicted target sites only partially relieved inhibition (Fig. 3 A and Fig. S1), suggesting the presence of additional binding sites in the IL12B transcript. In summary, these experiments validated IL12B as a direct target of multiple viral miRNAs.
IL-12 is critical for differentiation of Th1 cells (Szabo et al., 2003). Therefore, we co-cultured naive CD4+ T cells with autologous EBV-infected B cells (Fig. 3 B). Relative to ΔmiR EBV, wt/B95-8 EBV-infected B cells repressed Th1 differentiation (Fig. 3, C and D). An antibody that neutralizes the functions of IL12B, but not an isotype control antibody, suppressed Th1 differentiation when T cells were co-cultured with ΔmiR EBV-infected cells (Fig. 3 E), indicating that IL-12, secreted from EBV-infected and activated B cells, was responsible for generation of Th1 cells. Thus, EBV miRNAs suppress the release of IL-12 from infected cells and thereby interfere with formation of Th1 cells, which are important antiviral effectors.
Viral miRNAs directly and indirectly control antigen presentation
Having shown that EBV miRNAs interfere with CD4+ T cell differentiation, we turned our attention to molecules that are involved in recognition of infected cells by specific CD4+ T cells. We quantified levels of surface proteins with a role in HLA class II antigen presentation (Fig. 4). All three subclasses of HLA class II (HLA-DR, HLA-DQ, and HLA-DP) tested were reduced in wt/B95-8 relative to ΔmiR EBV-infected B cells (Fig. 4, A and B), as were many co-stimulatory and adhesion molecules 5 and 15 dpi (Fig. 4 C). MHC class I molecules were also affected, but to a lesser extent.
Among the many co-receptors and adhesion molecules down-regulated in cells infected with wt/B95-8 EBV (Fig. 4 C), we searched for direct targets of EBV miRNAs but found only CD40 and FAS in RISC-IPs (Table S1), which we could not confirm in subsequent luciferase assays using miRNAs encoded by wt/B95-8 EBV. Interestingly, the viral latent membrane protein 1 (LMP1) activates the CD40 pathway, inducing important immune co-receptors (Kieser and Sterz, 2015), but several viral BART miRNAs were reported to control LMP1 expression (Lo et al., 2007; Riley et al., 2012; Verhoeven et al., 2016). We tested these findings in our model of newly infected B cells, and found reduced but highly variable levels of LMP1 transcripts (Fig. 4 D) and protein (Fig. 4 E) in B cells infected with wt/B95-8 EBV compared with ΔmiR EBV 5 dpi. We identified miR-BART3 and miR-BART16 as inhibiting LMP1 in reporter assays (Fig. 4 F). miR-BART3 is encoded in wt/B95-8 EBV, whereas miR-BART16 (Fig. 4 F) is only present in field strains of EBV. We also tested miR-BART1 and miR-BART17, which were reported together with miR-BART16 to target LMP1 3′-UTR (Lo et al., 2007), but failed to confirm that miR-BART1 and miR-BART17 target LMP1 (Fig. 4 F). Collectively, our results showed that viral miRNAs limit LMP1 gene expression and thereby indirectly inhibit surface expression of some immune co-receptors and adhesion molecules.
Viral miRNAs target lysosomal enzymes and inhibit antigen processing
According to our transcriptome analysis, genes encoding lysosomal enzymes actively involved in MHC class II peptide processing (Blum et al., 2013) were inhibited by EBV miRNAs (Fig. 1 C). These included IFI30 (coding for IFN-γ–regulated thiol reductase GILT), LGMN (coding for asparagine endopeptidase AEP alias legumain), and CTSB (coding for the peptidase cathepsin B). Expression of all three genes was reduced by EBV miRNAs (Fig. 1 C), which we verified by quantitative RT-PCR (Fig. 5 A). We found that EBV’s miR-BART1, miR-BART2, and miR-BHRF1-2 could directly regulate IFI30, LGMN, and CTSB gene expression via their 3′-UTRs in luciferase reporter assays (Fig. 5 B and Fig. S1). Importantly, the knock-down of these three genes (Fig. 5 C) resulted in reduced antigen presentation of exogenously loaded protein (Fig. 5 D). Collectively, our results show that EBV miRNAs interfere with processes involved in MHC class II antigen presentation at multiple levels, including lysosomal protein degradation, HLA class II expression, and co-stimulatory molecule expression.
EBV miRNAs inhibit recognition of infected B cells by EBV-specific CD4+ T cells
Next, we asked whether these multiple levels of regulation ultimately resulted in reduced MHC class II–mediated recognition of EBV-infected cells by antiviral CD4+ T cells. CD4+ T cells from EBV-positive individuals were enriched for EBV-specific T cells by repeated stimulation with irradiated wt/B95-8 EBV-infected autologous LCLs. The EBV-specific CD4+ T cells were then co-cultured with autologous B cells that had been infected with the two EBV strains 5 d earlier (Fig. 6 A, top). Release of IFN-γ by EBV-specific CD4+ T cells was substantial when co-cultured with ΔmiR EBV-infected cells as targets, but was consistently reduced when co-cultured with wt/B95-8 EBV-infected B cells at all cell ratios tested (Fig. 6 B). Activation of EBV-specific CD4+ T cells, measured as IFN-γ release, was observed in autologous and partially matched but not in HLA-mismatched conditions (Fig. 6 C and Table S2), indicating that the activation was HLA class II-restricted.
We also tested an antigen-specific CD4+ T cell clone (Fig. 6, A [bottom] and D) directed against the FGQ peptide, an epitope derived from the viral glycoprotein gp350 (Adhikary et al., 2006). We observed dramatically reduced T cell activities with target B cells infected with wt/B95-8 EBV compared with ΔmiR EBVs 5 dpi (Fig. 6 D). T cell activities were much reduced at 15 dpi, but a difference between wt/B95-8 and ΔmiR EBV-infected target cells remained detectable. Weak recognition on day 15 is in accordance with gp350 protein being delivered as a component of the virion (Adhikary et al., 2006) but not synthesized during prelatency or latency (Kalla et al., 2010). Interestingly, expression of cell surface HLA class II levels peaked between 4 to 10 dpi (Fig. 6 E) suggesting the importance of viral miRNAs that counteract CD4+ T cell recognition in the early days of infection.
EBV-specific CD4+ T cells have cytolytic activity (Adhikary et al., 2006). In allogeneic, partially HLA-matched conditions, EBV-specific CD4+ T cells consistently showed stronger cytolysis of target B cells infected with ΔmiR EBV than cells infected wt/B95-8 EBV (Fig. 6 F).
Collectively, we have discovered that EBV miRNAs inhibit the recognition and elimination of infected B cells by HLA class II–restricted CD4+ T cells. Apparently, EBV utilizes multiple miRNAs to interfere with proinflammatory cytokines, antigen processing, and epitope presentation of the infected B lymphocyte to evade EBV-specific and antiviral CD4+ T cells responses early after infection.
DISCUSSION
EBV infects its human hosts for their lifetime, residing in nonproliferating B cells largely invisible to the host's immune response (Thorley-Lawson, 2005). EBV, to be the successful pathogen that it is, however, must both establish a latent infection and produce and disseminate progeny virus, all in the face of robust innate and adaptive immune responses.
Such responses include EBV-specific CD4+ T cells, which have an important role in controlling EBV infection and disease. For example, patients with EBV-associated tumors treated with virus-specific T cell preparations showed better clinical responses if the preparations contained larger fractions of CD4+ T cells (Haque et al., 2007; Icheva et al., 2013). CD4+ T cells target a wide repertoire of EBV antigens from all phases of latent and lytic infection (Adhikary et al., 2007; Long et al., 2011). CD4+ T cells with specificity for structural EBV proteins play a prominent role: they are a universal component of the T cell repertoire, rapidly detecting EBV-infected B cells and killing them directly (Adhikary et al., 2006, 2007). Several EBV proteins expressed during its lytic phase can inhibit recognition of EBV-infected cells by CD4+ T cells (Ressing et al., 2015). The broad-ranging functions of the host shut-off gene product BGLF5 include reduction of HLA class II molecules on the cell surface during EBV’s lytic, productive phase (Rowe et al., 2007). The late glycoprotein gp42, encoded by BZLF2, was shown to associate with HLA class II and to hinder recognition by CD4+ T cells sterically (Ressing et al., 2003). Both mechanisms are unlikely to be operational in newly infected B lymphocytes, because the two viral proteins appear not to be expressed in the prelatent phase (Kalla et al., 2010). Viral IL-10, encoded by BCLF1, is an immunomodulatory protein expressed early in infection, but its effects on B cell elimination by CD4+ T cells were limited (Jochum et al., 2012). Two additional viral gene products reported to affect CD4+ T cell recognition, BDLF3 (Quinn et al., 2015) and BZLF1 (Zuo et al., 2011), may act during EBV primary infection, because BDLF3 is in the virus particle (Johannsen et al., 2004) and BZLF1 has been found to be expressed early during infection (Wen et al., 2007; Kalla et al., 2010). BDLF3 transcripts were present at very low levels, only, whereas BZLF1 transcripts were not mapped in our RNA-Seq analysis. Thus, how EBV infection escapes detection and elimination by EBV-specific T cells during the early phase of infection has remained uncertain.
Here, we present an answer to this question and show that EBV uses its large repertoire of miRNAs to target CD4+ T cell differentiation and recognition of infected cells. It appears that EBV’s immunoevasive strategy uses miRNAs, which are themselves nonimmunogenic (Boss and Renne, 2011), rather than viral proteins which themselves would be antigenic. EBV induces a state of reduced immunogenicity in infected and recently activated B cells with viral miRNAs, which allows the virus to express its latency-associated antigens avoiding the recognition and elimination by CD4+ T cells. Because activated B cells are professional antigen-presenting cells and express multiple immune-activating molecules (Wiesner et al., 2008), the need for EBV to control its host cell is more urgent than for other complex viruses that do not rely on professional antigen-presenting immune cells for their life-cycle.
In our experiments, recognition of early-stage infected B cells by CD4+ T cells was strongly inhibited by multiple mechanisms, pointing to the biological importance of the immunoevasive functions of viral miRNAs. First, IL12B is dramatically repressed in wt/B95-8 EBV-infected B cells compared with ΔmiR EBV-infected B cells, leading to the down-regulation of three IL-12 family cytokines, IL-12B (IL-12p40), IL-12 (p35/p40), and IL-23 (p19/p40). At least five different viral miRNAs control the fate of the IL12B transcript, targeting multiple sites within its 3′-UTR, indicating its critical role in immune regulation by EBV-infected B cells and a redundancy or even cooperativity of viral miRNAs. Interestingly, the miRNAs controlling IL12B originate from different viral transcripts (Barth et al., 2011), suggesting a robust control of IL12B in all phases of EBV’s life cycle and in the many cell types EBV infects, which show different patterns of miRNA expression (Cai et al., 2006; Qiu et al., 2011). Repression of IL12B may not only reduce CD4+ T cell differentiation, as shown here, but also regulate T cell effector functions (Curtsinger and Mescher, 2010).
Second, lysosomal proteolysis is regulated. Transcripts of lysosomal endopeptidases, AEP and CTSB, and a thiol reductase, GILT, which are involved in proteolytic degradation and HLA class II epitope generation (Blum et al., 2013), are direct targets of the miRNAs miR-BART1, miR-BART2, and miR-BHRF1-2. An siRNA-mediated knock-down of these three genes in human B cells reduced their recognition by clonal epitope-specific CD4+ T cells (Fig. 5, C and D; Milosevic et al., 2005) suggesting an important role of the three lysosomal enzymes in antigen processing and presentation via MHC class II molecules in human B cells.
Third, HLA class II surface levels are down-regulated in B cells infected with wt/B95-8 EBV (Fig. 4 B). We also tested if EBV miRNA targeted MHC class II molecules directly. miRNAs encoded in wt/B95-8 EBV failed to inhibit consistently four HLA-DRB1 alleles tested in dual luciferase reporter assays and also lacked functional miRNA-binding sites. Thus, the reduction of HLA class II molecules is likely indirect and may be a consequence of altered lysosomal processing of epitopes. Such phenomena might be particularly important early during infection when the EBV-activated B cells present antigenic peptides of virion components (Fig. 6 D) and the expression of HLA class II molecules peaks at the cell surface (Fig. 6 E). Together, these findings demonstrate that EBV miRNAs redundantly and robustly inhibit specific immune functions in newly infected B cells that may otherwise metabolize viral proteins into HLA class II-presented peptides for recognition by antiviral CD4+ T cells.
A limited number of observations have been made before on immunomodulatory functions of EBV miRNAs. The RNAs encoding the innate immune effector molecules MICB (Nachmani et al., 2009) and NLRP3 (Haneklaus et al., 2012) have been found in cell lines to be inhibited by EBV's miRNAs. They were not down-regulated in our experiments (Fig. 1 C, Table S1, and GSE75776), which examined their levels during early infection of primary B cells. Our studies with primary B cells likely avoided the adaptive changes that arise during long-term culturing in vitro.
It is not immediately apparent why B lymphocytes release proinflammatory cytokines upon infection with EBV. One explanation is based on two viral, noncoding RNAs, termed EBERs, which are contained in virions and are transcribed in all EBV-infected cells. EBERs are known to trigger the endosomal TLR3 receptor or the cytosolic RIG-I sensor signaling pathway and induce type I interferon and IL-6 synthesis in infected B cells (Samanta et al., 2006; Wu et al., 2007; Iwakiri et al., 2009), which might lead to the expression of proinflammatory cytokines found in higher concentrations in the supernatants of B cells infected with ΔmiR than with wt/B95-8 EBV.
Interestingly, LMP1, a viral membrane protein that is predominantly expressed in the latent phase but also early upon B cell infection, activates the CD40 pathway, and induces IL-6, adhesion molecules, and important immune co-receptors (Kieser and Sterz, 2015). Several viral BART miRNAs have been reported to control LMP1 expression (Lo et al., 2007; Verhoeven et al., 2016). We confirmed that miR-BART16, which is not in the wt/B95-8 EBV strain, and miR-BART3, which was not previously known to target LMP1 directly, do target it (Fig. 4 F). In addition, miR-BART3 was recently reported to reduce LMP1 protein levels in HEK 293T cells (Verhoeven et al., 2016). We failed to identify direct targets among the many co-receptors and adhesion molecules down-regulated in cells infected with wt/B95-8 EBV (Fig. 4 C). It is possible that, on average, BART miRNAs repress LMP1 levels in cells infected with wt/B95-8 EBV, and thereby indirectly inhibit surface expression of some immune co-receptors and adhesion molecules, further reducing immune recognition of EBV-infected B lymphocytes.
EBV induces proliferation of the B cells it initially infects, and fosters their survival. We have found that EBV encodes miRNAs that regulate multiple facets of a host's adaptive immune response in newly infected B cells. EBV-infected B cells lacking viral miRNAs are deficient both in regulating these responses and in other miRNA-dependent functions, including an inhibition of apoptosis (Seto et al., 2010). These latter defects have precluded comparisons of B cells newly infected with wt/B95-8 or ΔmiR in humanized mouse models. In infection experiments with these mice, we observed defects in persistence after infection with ΔmiR EBV compared with wt/B95-8 EBV, possibly resulting from a combination of decreased survival and enhanced immune control (unpublished data; C. Münz, personal communication).
Collectively, our studies of a model that closely mimics physiological infection in the early phase show that EBV's miRNAs interfere with CD4+ T cell control through multiple mechanisms. They inhibit the secretion of cytokines, inhibit antigen processing and presentation, inhibit the differentiation of CD4+ T cells, and counteract recognition and elimination of infected B cells by EBV-specific CD4+ effector T cells. The breadth of EBV's use of its miRNAs to inhibit adaptive immune responses is unprecedented and contributes to its efficient establishment of a lifelong infection.
MATERIALS AND METHODS
Patient samples
Surgically removed adenoids and PBMCs were obtained from anonymous patients and anonymous volunteer blood donors, respectively, from Munich, Germany. The use of this human material was approved by the local ethics committee (Ethikkommission bei der LMU München) in writing.
Separation of human primary cells
Human primary B and T cells were prepared from adenoidal mononuclear cells (MNCs) or PBMCs by Ficoll-Hypaque gradient centrifugation with Pancoll (PAN-Biotech). B cells, CD4+ T cells, CD8+ T cells, and naive CD4+ T cells were separated from adenoidal MNCs or PBMCs using MACS separator (Miltenyi Biotec) with CD19 MicroBeads, CD4 MicroBeads, CD8 MicroBeads, and Naive CD4+ T cell Isolation kit II, respectively.
Cell lines and cell culture
Burkitt’s lymphoma cell lines Raji (EBV-positive), DG-75 (EBV-negative), HEK293-based EBV producer cell lines (Seto et al., 2010), infected human primary B cells, and T cells were maintained in RPMI-1640 medium (Thermo Fischer Scientific). HEK293T cells were maintained in DMEM medium. All media were supplemented with 10% FBS (Thermo Fischer Scientific), penicillin (100 U/ml; Thermo Fischer Scientific), and streptomycin (100 mg/ml; Thermo Fischer Scientific). Cells were cultivated at 37°C in a 5% CO2 incubator.
Preparation of infectious EBV stocks and infection of human primary B cells
Infectious EBV stocks were prepared as previously described (Seto et al., 2010). In brief, EBV producer cell lines for ΔmiR (p4027) and wt/B95-8 (p2089) EBV strains were transiently transfected with expression plasmids encoding BZLF1 and BALF4 to induce EBV’s lytic cycle. We collected supernatants 3 d after transfection, and debris was cleared by centrifugation at 3,000 rpm for 15 min. Virus stocks were titrated on Raji cells as previously reported and used at a multiplicity of infection (MOI) of 0.1 Green Raji units (Seto et al., 2010) for infecting primary B lymphocytes with an optimal virus dose (Steinbrück et al., 2015). For virus infection, primary B cells were cultivated with each virus stock for 18 h. After replacement with fresh medium, the infected cells were seeded at an initial density of 5 × 105 cells per ml.
RNA-Seq and RISC-IP
At 5 dpi of human primary B cells, we extracted total RNAs with TRIzol (Thermo Fischer Scientific) and Direct-Zol RNA MiniPrep kit (Zymo Research) from six different donors (Ad1–Ad6; Fig. 1) for RNA-Seq, according to the manufacturers’ protocols. In parallel, we performed RISC immunoprecipitation (RISC-IP) as described previously (Kuzembayeva et al., 2012). In brief, lysed cells were incubated with anti-Ago2 antibody (11A9)-conjugated Dynabeads (Thermo Fischer Scientific), washed, and coprecipitated RNA was extracted. The cDNA libraries were prepared (Vertis Biotechnologie AG). For RNA-Seq, total RNAs were depleted of rRNAs by Ribo-Zero rRNA Removal kit (Illumina), fragmented by ultrasonication, and subjected to first strand synthesis with a randomized primer. For RISC-IP, RNAs were poly (A)-tailed, ligated with an RNA adapter at 5′-phosphates to facilitate Illumina TruSeq sequencing, and subjected to first strand synthesis with an oligo-(dT) primer. The cDNAs were PCR-amplified and sequenced with an Illumina HiSeq2000 instrument at the University of Wisconsin Biotechnology Center DNA Sequencing Facility.
Analysis of deep sequencing
For RNA-Seq, processing of paired-end reads (poly-A tail filtering, N-filtering, and adapter removal) was done using FastQC and R2M (RawReadManipulator). Reads were mapped to the human genome (hg19 ‘core’ chromosome-set) by STAR (Dobin et al., 2013) and feature counts per transcript were determined using featureCounts and GENCODE version 19 annotations, together with EBV’s annotation (available from GenBank under accession no. AJ507799). To screen differentially regulated genes by viral miRNAs, we used a simple but efficient scoring algorithm based on donor/replicate-wise fold changes ranks. For each gene g and replicate k, we calculate the gene-specific rank score as
where n is the number of all replicates, m the number of all genes/transcripts, and rgk is the rank of gene g in sample k. To select highly differentially expressed genes, we transformed the rank score into a z-score and selected all transcripts with an absolute z-score >1.6.
For RISC-IP the mapped reads were normalized using size factors estimated with the R package DEseq2 and filtered for reads mapped to annotated 3′-UTR regions using GENCODE version19. To identify local quantitative differences in the read enrichments on 3′-UTRs between wt EBV compared with ΔmiR EBV-infected B cells, we calculated a donor-wise relative enrichment score. For each genomic position p, the relative expression esp was calculated as
where etp is the enrichment value of sequenced reads at position p in wt/B95-8 EBV-infected cells and ecp the local enrichment value in ΔmiR EBV-infected B cells, respectively. The normalization factor npu = etp/max(eu) was introduced to correct for local maxima in the UTR sequence of interest, where max(eu) is the maximum enrichment value in the UTR sequence u. Finally, we used a Gaussian filter to minimize local noise. To select 3′-UTRs bound by viral miRNAs, we set the threshold as follows: enrichment score >0.6 for a stretch of >20 nt in the 3′-UTRs in two or more donors. To quantify viral miRNAs incorporated into the RISC in infected cells, we mapped reads from the RISC-IP Seq to miRNA entries registered in miRBase 21 and calculated fractions of each viral miRNAs out of total miRNA read counts.
KEGG enrichment pathway
Enrichment of specific pathways was estimated by performing a hypergeometric distribution test via the KEGG API Web Service. All calculations were done using Matlab (Mathworks).
ELISA
To detect cytokine secretion from infected B cells, 106 cells were seeded in 6-well plates at 4 or 11 dpi, cultivated for 4 d with cyclosporine (1 µg/ml; Novartis). Supernatants were harvested and stored at −20°C. ELISAs for IL-6, IL-10, IL12B (IL-12p40), IL-12, IL-23, and TNF were performed following the manufacturer’s protocols (Mabtech). For IL-6 and IL-10, CpG DNA were added as previously described (Iskra et al., 2010) to stimulate infected B cells. ELISA for IFN-γ levels was performed following the manufacturer’s protocol (Mabtech).
To detect IL-12p40 secretion from PBMCs or PBMCs depleted of B cells using the MACS separator and CD19 MicroBeads (Miltenyi Biotec), the cells were infected with either wt/B95-8 or ΔmiR EBV at MOIs of 0.1 Green Raji units (Steinbrück et al., 2015). After 5 d of incubation, supernatants were collected and ELISA for IL-12p40 levels was assessed following the manufacturer’s protocol (Mabtech).
Luciferase reporter assays
The 3′-UTRs of IL12B (Ensembl ENST00000231228), IFI30 (Ensembl ENST00000407280), LGMN (Ensembl ENST00000334869), CTSB (Ensembl ENST00000353047), and LMP1 (available from GenBank under accession no. AJ507799) were cloned downstream of Renilla luciferase (Rluc) in the expression plasmid psiCHECK-2 (Promega). To construct the viral miRNA expression vectors, we cloned TagBFP (Evrogen) under the control of the EF1α promoter into pCDH-EF1-MCS (System Biosciences). Single miRNAs of interest were cloned downstream of the TagBFP-encoding gene. Viral miRNAs were obtained by PCR from the p4080 plasmid (Seto et al., 2010). 50 ng of the psiCHECK-2 reporter and 150 ng of the pCDH-EF1 miRNA expressor plasmid DNAs were cotransfected into 1 × 105 HEK293T cells by Metafectene Pro (Biontex). After 24 h of transfection, we measured luciferase activities with the Dual-Luciferase Assay kit (Promega) and the Orion II Microplate Luminometer (Titertek-Berthold). The activity of Rluc was normalized to the activity of Fluc (Firefly luciferase) encoded in the psiCHECK-2 reporter plasmid. We performed in silico prediction of EBV miRNA-binding sites on 3′-UTRs primarily with TargetScan (Garcia et al., 2011) and used RNAhybrid (Rehmsmeier et al., 2004) to screen for 6mer binding sites (Bartel, 2009). We performed site-directed mutagenesis with overlapping oligo DNAs and Phusion polymerase (NEB).
Quantitative RT-PCR
To quantify mRNA levels RNAs were reverse-transcribed with SuperScript III Reverse transcription (Thermo Fischer Scientific) and quantitative PCR was performed with LightCycler 480 SYBR Green I Mix (Roche) and LightCycler 480 Instrument II (Roche) according to the manufacturers’ instructions. The following primers were used for the detection: HPRT1 5′-TGACCTTGATTTATTTTGCATACC-3′ and 5′-CGAGCAAGACGTTCAGTCCT-3′; HMBS 5′-CTGAAAGGGCCTTCCTGAG-3′ and 5′-CAGACTCCTCCAGTCAGGTACA-3′; IL12B 5′-CCCTGACATTCTGCGTTCA-3′ and 5′-AGGTCTTGTCCGTGAAGACTCTA-3′; IFI30 5′-CTGGGTCACCGTCAATGG-3′ and 5′-GCTTCTTGCCCTGGTACAAC-3′; LGMN 5′-GGAAACTGATGAACACCAATGA-3′ and 5′-GGAGACGATCTTACGCACTGA-3′; CTSB 5′-CTGTGGCAGCATGTGTGG-3′ and 5′-TCTTGTCCAGAAGTTCCAAGC-3′.
To quantify miRNA levels, stem-loop qPCRs were performed with TaqMan MicroRNA Reverse Transcription kit (Thermo Fischer Scientific) and TaqMan Universal Master Mix II (Thermo Fischer Scientific) according to the manufacturer’s protocols. RNU6B was used for normalization. Following TaqMan MicroRNA assays, specific primers (Thermo Fischer Scientific) were used for detection: ebv-miR-BART2: 197238_mat; ebv-miR-BART3: 004578_mat; ebv-miR-BHRF1-3 197221_mat; RNU6B: 197238_mat.
Establishment of EBV-specific effector T cells and T cell clones
EBV-specific CD4+ T cell clones were established from polyclonal T cell lines that were generated by LCLs or mini-LCL stimulation of PBMCs, as previously described (Adhikary et al., 2007).
Flow cytometry and antibodies
After immunostainings with fluorophore-conjugated antibodies, single-cell suspensions were measured with LSRFortessa or FACSCanto (BD) flow cytometers and the FACSDiva software (BD). Acquired data were analyzed with FlowJo software Ver. 9.8 (FlowJo). The following fluorophore-conjugated antibodies reactive to human antigens were used: anti–human IFN-γ APC (4S.B3, IgG1; BioLegend), anti-CD40 PE (5c3, IgG2b; BioLegend), anti–ICOS-L (B7-H2) PE (2D3, IgG2b; BioLegend), anti–PD-L1 (B7-H1) APC (29E.2A3, IgG2b; BioLegend), anti-CD86 (B7-2) PE (37301, IgG1; R&D Systems), anti-CD54 (ICAM-1) APC (HCD54, IgG1; BioLegend), anti–HLA-ABC APC (W6/32, IgG2a; BioLegend), anti-CD80 PE-Cy5 (L307.4; BD), anti-FAS (CD45) PE (Dx2, IgG1; BioLegend), anti–HLA-DR unlabeled (L234, IgG2a; BioLegend), anti–HLA-DQ unlabeled (SPV-L3, IgG2a ; AbD Serotec), anti–HLA-DP unlabeled (B7/21, IgG3; Abcam), anti–mouse F(ab')2 APC (polyclonal, IgG; eBioscience), isotype IgG1 PE (MOPC-21; BioLegend), isotype IgG2b PE (MPC-11; BioLegend), isotype IgG1 APC (MOPC-21; BD), isotype IgG2a APC (MOPC-173; BioLegend), and isotype IgG2b APC (MG2b-57; BioLegend).
Western blotting
We lysed cells with RIPA buffer (50 mM Tris-HCl, pH 8.0, 150 mM NaCl, 0.1% SDS, 1% NP-40, and 0.5% DOC) and boiled the extracts with Laemmli buffer. Proteins were separated on SDS-PAGE gels (Carl Roth) and transferred to nitrocellulose membranes (GE Healthcare) using Mini-PROTEAN Tetra Cell (Bio-Rad Laboratories). Membranes were blocked for 30 min with Roti-Block (Carl Roth), followed by antibody incubation. Secondary antibodies conjugated with horseradish peroxidase were used (Cell Signaling Technology) and exposed to CEA films (Agfa HealthCare). Protein levels were quantified with the software ImageJ. The following primary antibodies reactive to human proteins were used: anti–human Tubulin (B-5-1-2; Santa Cruz Biotechnology, Inc.). The monoclonal antibody (1G6-3) reactive to the EBV protein LMP1 was provided by E. Kremmer (Institute of Molecular Immunology, Helmholtz Zentrum München, München, Germany).
RNAi knock-down and recognition by M1-specific CD4+ T cells
4 × 105 DG-75 cells were incubated in 1 ml Accell Delivery Media (GE Healthcare) and 1 nmol siRNAs directed against GAPDH, IFI30, LGMN, CTSB, or combinations thereof for 48 h. Influenza M1 protein purified as previously described (Nimmerjahn et al., 2003) was added to the medium, and the cells were further incubated for 24 h and co-cultured with M1-specific CD4+ T cell clone E5 for 16 h (Milosevic et al., 2005). IFN-γ levels were detected with ELISA.
T cell differentiation and recognition
Th1 differentiation was assessed by co-culture of sorted naive CD4+ T cells and infected B cells 5 dpi. 1 × 105 naive CD4+ T cells stained with CellTrace Violet (Thermo Fischer Scientific) and 0.5 or 1 × 105 infected B cells were cultured in 96-well plates with Dynabeads Human T-Activator CD3/CD28 (Thermo Fischer Scientific) and cultivated for 7 d. The neutralizing antibody against IL12B (C8.6; BioLegend) or the corresponding isotype control antibody (MOPC-21; BioLegend) were added for certain experiments at 5 µg/ml. Cells were restimulated with PMA and ionomycin (Cell Stimulation Cocktail; eBioscience) for 5 h and treated with Brefeldin A and Monensin (BioLegend) for 2.5 h before fixation. Th1 population was measured by intracellular IFN-γ staining with FIX and PERM Cell Fixation and Cell Permeabilization kit (Thermo Fischer Scientific) and subsequent flow cytometry analysis. The Th1 population was defined as IFN-γ+ T cells in the fraction of proliferating T cells identified via CellTrace Violet staining. EBV-specific effector T cells’ activities were measured with ELISA and Calcein release assays. For IFN-γ detection from T cells, effector and target cells were seeded at 5 × 104 cell per ml (1:1 ratio) each and co-cultured for 16 h in a 96-well plate (V bottom). IFN-γ levels were detected with ELISA. IFN-γ concentrations <16 pg/ml were considered as not detected.
T cell cytotoxicity assays
Primary infected B cells were purified by Ficoll-Hypaque gradient centrifugation, and 5 × 105 target cells were labeled with calcein at 0.5 µg/ml. After three washing steps with PBS, target and effector cells were co-cultured in a 96-well plate (V bottom) with different ratios in RPMI red phenol-free medium to reduce background signals. After 4 h of co-culture, fluorescence intensity of the released calcein was measured by the Infinite F200 PRO fluorometer (Tecan). As controls, spontaneous calcein release of target cells cultivated without effector cells and cells lysed with 0.5% Triton-X100 were used to define the levels of no and fully lysed target cells, respectively.
Statistical analysis
We used Prism 6.0 software (GraphPad) for the statistical analysis. A two-tailed ratio Student’s t test was applied unless otherwise mentioned.
Online supplemental material
Fig. S1 shows the predicted miRNA-binding sites and mutations tested in 3′-UTR reporter assays. Table S1, available as an Excel file, lists the gene transcripts controlled by viral miRNAs. Table S2 lists the HLA allele information of donors used in co-culture experiments.
ACKNOWLEDGMENTS
We thank Christian Münz, Zurich, Elisabeth Kremmer, Anne-Wiebe Mohr, and Liridona Maliqi, Munich, for animal experiments, antibodies, T cell clones, and experimental assistance, respectively. We also thank Dagmar Pich for her experimental expertise and advice.
This work was financially supported by the Deutsche Forschungsgemeinschaft (SFB1054/TP B05 and TP A03, SFB1064/TP A13, and SFB-TR36/TP A04), Deutsche Krebshilfe (107277 and 109661), National Institutes of Health (RO1: CA70723 and PO1: CA022443), and personal grants of Deutscher Akademischer Austauschdienst to T. Tagawa (Studienstipendien für ausländische Graduierte aller wissenschaftlichen Fächer) and European Molecular Biology Organization to M. Bouvet.
The authors declare no competing financial interests.
References
Author notes
T. Tagawa and M. Albanese contributed equally to this paper.