Host-encoded CTCF regulates human cytomegalovirus latency via chromatin looping.

Human cytomegalovirus (HCMV) is a prevalent pathogen that establishes life-long latent infection in hematopoietic cells. While this infection is usually asymptomatic, immune dysregulation leads to viral reactivation, which can cause significant morbidity and mortality. However, the mechanisms underpinning reactivation remain incompletely understood. The HCMV major immediate early promoter (MIEP)/enhancer is a key factor in this process, as its transactivation from a repressed to active state helps drive viral gene transcription necessary for reactivation from latency. Numerous host transcription factors bind the MIE locus and recruit repressive chromatin modifiers, thus impeding virus reactivation. One such factor is CCCTC-binding protein (CTCF), a highly conserved host zinc finger protein that mediates chromatin conformation and nuclear architecture. However, the mechanisms by which CTCF contributes to HCMV latency were previously unexplored. Here, we confirm that CTCF binds two convergent sites within the MIE locus during latency in primary CD14+ monocytes, and following cellular differentiation, CTCF association is lost as the virus reactivates. While mutation of the MIE enhancer CTCF binding site does not impact viral lytic growth in fibroblasts, this mutant virus fails to maintain latency in myeloid cells. Furthermore, we show the two convergent CTCF binding sites allow looping to occur across the MIEP, supporting transcriptional repression during latency. Indeed, looping between the two sites diminishes during virus reactivation, concurrent with activation of MIE transcription. Taken together, our data reveal that three-dimensional chromatin looping aids in the regulation of HCMV latency and provides insight into promoter/enhancer regulation that may prove broadly applicable across biological systems.

Chromatin control of human cytomegalovirus infection.

Human cytomegalovirus (HCMV) is a betaherpesvirus that establishes lifelong infection in its host and can cause severe comorbidities in individuals with suppressed or compromised immune systems. The lifecycle of HCMV consists of lytic and latent phases, largely dependent upon the cell type infected and whether transcription from the major immediate early locus can ensue. Control of this locus, which acts as a critical "switch" region from where the lytic gene expression cascade originates, as well as regulation of the additional ~235 kilobases of virus genome, occurs through chromatinization with cellular histone proteins after infection. Upon infection of a host cell, an initial intrinsic antiviral response represses gene expression from the incoming genome, which is relieved in permissive cells by viral and host factors in concert. Latency is established in a subset of hematopoietic cells, during which viral transcription is largely repressed while the genome is maintained. As these latently infected cells differentiate, the cellular milieu and epigenetic modifications change, giving rise to the initial stages of virus reactivation from latency. Thus, throughout the cycle of infection, chromatinization, chromatin modifiers, and the recruitment of specific transcription factors influence the expression of genes from the HCMV genome. In this review, we discuss epigenetic regulation of the HCMV genome during the different phases of infection, with an emphasis on recent reports that add to our current perspective.

Cytomegalovirus US28 regulates cellular EphA2 to maintain viral latency.

Cytomegalovirus (CMV) reactivation from latency following immune dysregulation remains a serious risk for patients, often causing substantial morbidity and mortality. Here, we demonstrate the CMV-encoded G protein-coupled receptor, US28, in coordination with cellular Ephrin receptor A2, attenuates mitogen-activated protein kinase signaling, thereby limiting viral replication in latently infected primary monocytes. Furthermore, treatment of latently infected primary monocytes with dasatinib, a Food and Drug Association-approved kinase inhibitor used to treat a subset of leukemias, results in CMV reactivation. These ex vivo data correlate with our retrospective analyses of the Explorys electronic health record database, where we find dasatinib treatment is associated with a significant risk of CMV-associated disease (odds ratio 1.58, P = 0.0004). Collectively, our findings elucidate a signaling pathway that plays a central role in the balance between CMV latency and reactivation and identifies a common therapeutic cancer treatment that elevates the risk of CMV-associated disease.

Short- and long-range cis interactions between integrated HPV genomes and cellular chromatin dysregulate host gene expression in early cervical carcinogenesis.

Development of cervical cancer is directly associated with integration of human papillomavirus (HPV) genomes into host chromosomes and subsequent modulation of HPV oncogene expression, which correlates with multi-layered epigenetic changes at the integrated HPV genomes. However, the process of integration itself and dysregulation of host gene expression at sites of integration in our model of HPV16 integrant clone natural selection has remained enigmatic. We now show, using a state-of-the-art 'HPV integrated site capture' (HISC) technique, that integration likely occurs through microhomology-mediated repair (MHMR) mechanisms via either a direct process, resulting in host sequence deletion (in our case, partially homozygously) or via a 'looping' mechanism by which flanking host regions become amplified. Furthermore, using our 'HPV16-specific Region Capture Hi-C' technique, we have determined that chromatin interactions between the integrated virus genome and host chromosomes, both at short- (<500 kbp) and long-range (>500 kbp), appear to drive local host gene dysregulation through the disruption of host:host interactions within (but not exceeding) host structures known as topologically associating domains (TADs). This mechanism of HPV-induced host gene expression modulation indicates that integration of virus genomes near to or within a 'cancer-causing gene' is not essential to influence their expression and that these modifications to genome interactions could have a major role in selection of HPV integrants at the early stage of cervical neoplastic progression.

A BMPR2/YY1 Signaling Axis Is Required for Human Cytomegalovirus Latency in Undifferentiated Myeloid Cells.

Human cytomegalovirus (HCMV) presents a major health burden in the immunocompromised and in stem cell transplant medicine. A lack of understanding about the mechanisms of HCMV latency in undifferentiated CD34+ stem cells, and how latency is broken for the virus to enter the lytic phase of its infective cycle, has hampered the development of essential therapeutics. Using a human induced pluripotent stem cell (iPSC) model of HCMV latency and patient-derived myeloid cell progenitors, we demonstrate that bone morphogenetic protein receptor type 2 (BMPR2) is necessary for HCMV latency. In addition, we define a crucial role for the transcription factor Yin Yang 1 (YY1) in HCMV latency; high levels of YY1 are maintained in latently infected cells as a result of BMPR2 signaling through the SMAD4/SMAD6 axis. Activation of SMAD4/6, through BMPR2, inhibits TGFbeta receptor signaling, which leads to the degradation of YY1 via induction of a cellular microRNA (miRNA), hsa-miR-29a. Pharmacological targeting of BMPR2 in progenitor cells results in the degradation of YY1 and an inability to maintain latency and renders cells susceptible to T cell killing. These data argue that BMPR2 plays a role in HCMV latency and is a new potential therapeutic target for maintaining or disrupting HCMV latency in myeloid progenitors. IMPORTANCE Understanding the mechanisms which regulate HCMV latency could allow therapeutic targeting of the latent virus reservoir from where virus reactivation can cause severe disease. We show that the BMPR2/TGFbeta receptor/YY1 signaling axis is crucial to maintain HCMV latency in undifferentiated cells and that pharmacological reduction of BMPR2 in latently infected cells leads to reactivation of the viral lytic transcription program, which renders the infected cell open to immune detection and clearance in infected individuals. Therefore, this work identifies key host-virus interactions which regulate HCMV latent infection. It also demonstrates a potential new therapeutic approach to reduce HCMV reactivation-mediated disease by the treatment of donor stem cells/organs prior to transplantation, which could have a major impact in the transplant disease setting.

Bromodomain proteins regulate human cytomegalovirus latency and reactivation allowing epigenetic therapeutic intervention.

Reactivation of human cytomegalovirus (HCMV) from latency is a major health consideration for recipients of stem-cell and solid organ transplantations. With over 200,000 transplants taking place globally per annum, virus reactivation can occur in more than 50% of cases leading to loss of grafts as well as serious morbidity and even mortality. Here, we present the most extensive screening to date of epigenetic inhibitors on HCMV latently infected cells and find that histone deacetylase inhibitors (HDACis) and bromodomain inhibitors are broadly effective at inducing virus immediate early gene expression. However, while HDACis, such as myeloid-selective CHR-4487, lead to production of infectious virions, inhibitors of bromodomain (BRD) and extraterminal proteins (I-BETs), including GSK726, restrict full reactivation. Mechanistically, we show that BET proteins (BRDs) are pivotally connected to regulation of HCMV latency and reactivation. Through BRD4 interaction, the transcriptional activator complex P-TEFb (CDK9/CycT1) is sequestered by repressive complexes during HCMV latency. Consequently, I-BETs allow release of P-TEFb and subsequent recruitment to promoters via the superelongation complex (SEC), inducing transcription of HCMV lytic genes encoding immunogenic antigens from otherwise latently infected cells. Surprisingly, this occurs without inducing many viral immunoevasins and, importantly, while also restricting viral DNA replication and full HCMV reactivation. Therefore, this pattern of HCMV transcriptional dysregulation allows effective cytotoxic immune targeting and killing of latently infected cells, thus reducing the latent virus genome load. This approach could be safely used to pre-emptively purge the virus latent reservoir prior to transplantation, thereby reducing HCMV reactivation-related morbidity and mortality.

Human cytomegalovirus major immediate early transcripts arise predominantly from the canonical major immediate early promoter in reactivating progenitor-derived dendritic cells.

Human cytomegalovirus latency and reactivation is a major source of morbidity in immune-suppressed patient populations. Lifelong latent infections are established in CD34+progenitor cells in the bone marrow, which are hallmarked by a lack of major lytic gene expression, genome replication and virus production. A number of studies have shown that inhibition of the major immediate early promoter (MIEP) - the promoter that regulates immediate early (IE) gene expression - is important for the establishment of latency and that, by extension, reactivation requires reversal of this repression of the MIEP. The identification of novel promoters (termed ip1 and ip2) downstream of the MIEP that can drive IE gene expression has led to speculation over the precise role of the MIEP in reactivation. In this study we show that IE transcripts arise from both the MIEP and ip2 promoter in the THP1 cell macrophage cell line and also CD14+monocytes stimulated with phorbol ester. In contrast, we show that in in vitro generated dendritic cells or macrophages that support HCMV reactivation IE transcripts arise predominantly from the MIEP and not the intronic promoters. Furthermore, inhibition of histone modifying enzyme activity confirms the view that the MIEP is predominantly regulated by the activity of cellular chromatin. Finally, we observe that ip2-derived IE transcription is cycloheximide-sensitive in reactivating DCs, behaviour consistent with an early gene designation. Taken together, these data argue that MIEP activity is still important for HCMV reactivation but ip2 activity could play cell-type-specific roles in reactivation.

Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription.

The complex life cycle of oncogenic human papillomavirus (HPV) initiates in undifferentiated basal epithelial keratinocytes where expression of the E6 and E7 oncogenes is restricted. Upon epithelial differentiation, E6/E7 transcription is increased through unknown mechanisms to drive cellular proliferation required to support virus replication. We report that the chromatin-organising CCCTC-binding factor (CTCF) promotes the formation of a chromatin loop in the HPV genome that epigenetically represses viral enhancer activity controlling E6/E7 expression. CTCF-dependent looping is dependent on the expression of the CTCF-associated Yin Yang 1 (YY1) transcription factor and polycomb repressor complex (PRC) recruitment, resulting in trimethylation of histone H3 at lysine 27. We show that viral oncogene up-regulation during cellular differentiation results from YY1 down-regulation, disruption of viral genome looping, and a loss of epigenetic repression of viral enhancer activity. Our data therefore reveal a key role for CTCF-YY1-dependent looping in the HPV life cycle and identify a regulatory mechanism that could be disrupted in HPV carcinogenesis.

Human papillomavirus genome integration in squamous carcinogenesis: what have next-generation sequencing studies taught us?

Human papillomavirus (HPV) infection is associated with ∼5% of all human cancers, including a range of squamous cell carcinomas. Persistent infection by high-risk HPVs (HRHPVs) is associated with the integration of virus genomes (which are usually stably maintained as extrachromosomal episomes) into host chromosomes. Although HRHPV integration rates differ across human sites of infection, this process appears to be an important event in HPV-associated neoplastic progression, leading to deregulation of virus oncogene expression, host gene expression modulation, and further genomic instability. However, the mechanisms by which HRHPV integration occur and by which the subsequent gene expression changes take place are incompletely understood. The advent of next-generation sequencing (NGS) of both RNA and DNA has allowed powerful interrogation of the association of HRHPVs with human disease, including precise determination of the sites of integration and the genomic rearrangements at integration loci. In turn, these data have indicated that integration occurs through two main mechanisms: looping integration and direct insertion. Improved understanding of integration sites is allowing further investigation of the factors that provide a competitive advantage to some integrants during disease progression. Furthermore, advanced approaches to the generation of genome-wide samples have given novel insights into the three-dimensional interactions within the nucleus, which could act as another layer of epigenetic control of both virus and host transcription. It is hoped that further advances in NGS techniques and analysis will not only allow the examination of further unanswered questions regarding HPV infection, but also direct new approaches to treating HPV-associated human disease. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Epstein-Barr virus nuclear antigen EBNA-LP is essential for transforming naïve B cells, and facilitates recruitment of transcription factors to the viral genome.

The Epstein-Barr virus (EBV) nuclear antigen leader protein (EBNA-LP) is the first viral latency-associated protein produced after EBV infection of resting B cells. Its role in B cell transformation is poorly defined, but it has been reported to enhance gene activation by the EBV protein EBNA2 in vitro. We generated EBNA-LP knockout (LPKO) EBVs containing a STOP codon within each repeat unit of internal repeat 1 (IR1). EBNA-LP-mutant EBVs established lymphoblastoid cell lines (LCLs) from adult B cells at reduced efficiency, but not from umbilical cord B cells, which died approximately two weeks after infection. Adult B cells only established EBNA-LP-null LCLs with a memory (CD27+) phenotype. Quantitative PCR analysis of virus gene expression after infection identified both an altered ratio of the EBNA genes, and a dramatic reduction in transcript levels of both EBNA2-regulated virus genes (LMP1 and LMP2) and the EBNA2-independent EBER genes in the first 2 weeks. By 30 days post infection, LPKO transcription was the same as wild-type EBV. In contrast, EBNA2-regulated cellular genes were induced efficiently by LPKO viruses. Chromatin immunoprecipitation revealed that EBNA2 and the host transcription factors EBF1 and RBPJ were delayed in their recruitment to all viral latency promoters tested, whereas these same factors were recruited efficiently to several host genes, which exhibited increased EBNA2 recruitment. We conclude that EBNA-LP does not simply co-operate with EBNA2 in activating gene transcription, but rather facilitates the recruitment of several transcription factors to the viral genome, to enable transcription of virus latency genes. Additionally, our findings suggest that EBNA-LP is essential for the survival of EBV-infected naïve B cells.

Identification of host transcriptional networks showing concentration-dependent regulation by HPV16 E6 and E7 proteins in basal cervical squamous epithelial cells.

Development of cervical squamous cell carcinoma requires increased expression of the major high-risk human-papillomavirus (HPV) oncogenes E6 and E7 in basal cervical epithelial cells. We used a systems biology approach to identify host transcriptional networks in such cells and study the concentration-dependent changes produced by HPV16-E6 and -E7 oncoproteins. We investigated sample sets derived from the W12 model of cervical neoplastic progression, for which high quality phenotype/genotype data were available. We defined a gene co-expression matrix containing a small number of highly-connected hub nodes that controlled large numbers of downstream genes (regulons), indicating the scale-free nature of host gene co-expression in W12. We identified a small number of 'master regulators' for which downstream effector genes were significantly associated with protein levels of HPV16 E6 (n = 7) or HPV16 E7 (n = 5). We validated our data by depleting E6/E7 in relevant cells and by functional analysis of selected genes in vitro. We conclude that the network of transcriptional interactions in HPV16-infected basal-type cervical epithelium is regulated in a concentration-dependent manner by E6/E7, via a limited number of central master-regulators. These effects are likely to be significant in cervical carcinogenesis, where there is competitive selection of cells with elevated expression of virus oncoproteins.

Pathogenesis of human papillomavirus-associated mucosal disease.

Human papillomaviruses (HPVs) are a necessary cause of carcinoma of the cervix and other mucosal epithelia. Key events in high-risk HPV (HRHPV)-associated neoplastic progression include persistent infection, deregulated expression of virus early genes in basal epithelial cells and genomic instability causing secondary host genomic imbalances. There are multiple mechanisms by which deregulated virus early gene expression may be achieved. Integration of virus DNA into host chromosomes is observed in the majority of cervical squamous cell carcinomas (SCCs), although in ∼15% of cases the virus remains extrachromosomal (episomal). Interestingly, not all integration events provide a growth advantage to basal cervical epithelial cells or lead to increased levels of the virus oncogenes E6 and E7, when compared with episome-containing basal cells. The factors that provide a competitive advantage to some integrants, but not others, are complex and include virus and host contributions. Gene expression from integrated and episomal HRHPV is regulated through host epigenetic mechanisms affecting the virus long control region (LCR), which appear to be of functional importance. New approaches to treating HRHPV-associated mucosal neoplasia include knockout of integrated HRHPV DNA, depletion of virus transcripts and inhibition of virus early gene transcription through targeting or use of epigenetic modifiers. Copyright © 2014 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Virus transcript levels and cell growth rates after naturally occurring HPV16 integration events in basal cervical keratinocytes.

Cervical carcinogenesis is characterized by a clonal selection process in which the high-risk human papillomavirus (HRHPV) genome usually changes from the extra-chromosomal (episomal) state seen in productive infections to DNA that is integrated into host chromosomes. However, it is not clear whether all HRHPV integration events provide cells with a selective growth advantage compared with the episome-containing cells from which they originate. It is also unclear whether selection of cells containing a particular integrant from a mixed population simply reflects the highest levels of virus oncogene expression or has additional determinants. These early events in cervical carcinogenesis cannot readily be addressed by cross-sectional studies of clinical samples. We used the W12 model system to generate a panel of cervical squamous cell clones that were derived from an identical background under non-competitive conditions and differed only by the genomic site of HPV16 integration. Compared with the 'baseline' episome-containing cells from which they were isolated, only 9/17 clones (53%) showed significantly greater growth rates and only 7/17 (41%) showed significantly greater expression of the major virus oncogenes E7/E6. There were significant variations in levels of HPV16 transcription per DNA template, changes that were associated with histone modifications in the integrated virus chromatin. Cell growth rates showed only weak and non-significant associations with protein and mRNA levels for E7, E6, and the mean E7/E6 values. We conclude that HPV16 integration in basal cervical cells does not necessarily lead to increased levels of virus oncogenes, or to a competitive growth advantage, when compared with the initiating episome-containing cells.

Depletion of HPV16 early genes induces autophagy and senescence in a cervical carcinogenesis model, regardless of viral physical state.

In cervical carcinomas, high-risk human papillomavirus (HR-HPV) may be integrated into host chromosomes or remain extra-chromosomal (episomal). We used the W12 cervical keratinocyte model to investigate the effects of HPV16 early gene depletion on in vitro cervical carcinogenesis pathways, particularly effects shared by cells with episomal versus integrated HPV16 DNA. Importantly, we were able to study the specific cellular consequences of viral gene depletion by using short interfering RNAs known not to cause phenotypic or transcriptional off-target effects in keratinocytes. We found that while cervical neoplastic progression in vitro was characterized by dynamic changes in HPV16 transcript levels, viral early gene expression was required for cell survival at all stages of carcinogenesis, regardless of viral physical state, levels of early gene expression or histology in organotypic tissue culture. Moreover, HPV16 early gene depletion induced changes in host gene expression that were common to both episome-containing and integrant-containing cells. In particular, we observed up-regulation of autophagy genes, associated with enrichment of senescence and innate immune-response pathways, including the senescence-associated secretory phenotype (SASP). In keeping with these observations, HPV16 early gene depletion induced autophagy in both episome-containing and integrant-containing W12 cells, as evidenced by the appearance of autophagosomes, punctate expression of the autophagy marker LC3, conversion of LC3B-I to LC3B-II, and reduced levels of the autophagy substrate p62. Consistent with the reported association between autophagy and senescence pathways, HPV16 early gene depletion induced expression of the senescence marker beta-galactosidase and increased secretion of the SASP-related protein IGFBP3. Together, these data indicate that depleting HR-HPV early genes would be of potential therapeutic benefit in all cervical carcinogenesis pathways, regardless of viral physical state. In addition, the senescence/SASP response associated with autophagy induction may promote beneficial immune effects in bystander cells.

LIN28 Expression in malignant germ cell tumors downregulates let-7 and increases oncogene levels.

Despite their clinicopathologic heterogeneity, malignant germ cell tumors (GCT) share molecular abnormalities that are likely to be functionally important. In this study, we investigated the potential significance of downregulation of the let-7 family of tumor suppressor microRNAs in malignant GCTs. Microarray results from pediatric and adult samples (n = 45) showed that LIN28, the negative regulator of let-7 biogenesis, was abundant in malignant GCTs, regardless of patient age, tumor site, or histologic subtype. Indeed, a strong negative correlation existed between LIN28 and let-7 levels in specimens with matched datasets. Low let-7 levels were biologically significant, as the sequence complementary to the 2 to 7 nt common let-7 seed "GAGGUA" was enriched in the 3' untranslated regions of mRNAs upregulated in pediatric and adult malignant GCTs, compared with normal gonads (a mixture of germ cells and somatic cells). We identified 27 mRNA targets of let-7 that were upregulated in malignant GCT cells, confirming significant negative correlations with let-7 levels. Among 16 mRNAs examined in a largely independent set of specimens by quantitative reverse transcription PCR, we defined negative-associations with let-7e levels for six oncogenes, including MYCN, AURKB, CCNF, RRM2, MKI67, and C12orf5 (when including normal control tissues). Importantly, LIN28 depletion in malignant GCT cells restored let-7 levels and repressed all of these oncogenic let-7 mRNA targets, with LIN28 levels correlating with cell proliferation and MYCN levels. Conversely, ectopic expression of let-7e was sufficient to reduce proliferation and downregulate MYCN, AURKB, and LIN28, the latter via a double-negative feedback loop. We conclude that the LIN28/let-7 pathway has a critical pathobiologic role in malignant GCTs and therefore offers a promising target for therapeutic intervention.

Depletion of polycistronic transcripts using short interfering RNAs: cDNA synthesis method affects levels of non-targeted genes determined by quantitative PCR.

BACKGROUND: Short interfering RNAs (siRNAs) are often used to deplete viral polycistronic transcripts, such as those encoded by human papillomavirus (HPV). There are conflicting data in the literature concerning how siRNAs targeting one HPV gene can affect levels of other genes in the polycistronic transcripts. We hypothesised that the conflict might be partly explained by the method of cDNA synthesis used prior to transcript quantification. FINDINGS: We treated HPV16-positive cervical keratinocytes with siRNAs targeting the HPV16 E7 gene and used quantitative PCR to compare transcript levels of E7 with those of E6 and E2, viral genes located upstream and downstream of the target site respectively. We compared our findings from cDNA generated using oligo-dT primers alone with those from cDNA generated using a combination of random hexamer and oligo-dT primers. Our data show that when polycistronic transcripts are targeted by siRNAs, there is a period when untranslatable cleaved mRNA upstream of the siRNA binding site remains detectable by PCR, if cDNA is generated using random hexamer primers. Such false indications of mRNA abundance are avoided using oligo-dT primers. The period corresponds to the time taken for siRNA activity and degradation of the cleaved transcripts. Genes downstream of the siRNA binding site are detectable during this interval, regardless of how the cDNA is generated. CONCLUSIONS: These data emphasise the importance of the cDNA synthesis method used when measuring transcript abundance following siRNA depletion of polycistronic transcripts. They provide a partial explanation for erroneous reports suggesting that siRNAs targeting HPV E7 can have gene-specific effects.

Extensive co-operation between the Epstein-Barr virus EBNA3 proteins in the manipulation of host gene expression and epigenetic chromatin modification.

Epstein-Barr virus (EBV) is able to drive the transformation of B-cells, resulting in the generation of lymphoblastoid cell lines (LCLs) in vitro. EBV nuclear proteins EBNA3A and EBNA3C are necessary for efficient transformation, while EBNA3B is dispensable. We describe a transcriptome analysis of BL31 cells infected with a series of EBNA3-knockout EBVs, including one deleted for all three EBNA3 genes. Using Affymetrix Exon 1.0 ST microarrays analysed with the MMBGX algorithm, we have identified over 1000 genes whose regulation by EBV requires one of the EBNA3s. Remarkably, a third of the genes identified require more than one EBNA3 for their regulation, predominantly EBNA3C co-operating with either EBNA3B, EBNA3A or both. The microarray was validated by real-time PCR, while ChIP analysis of a selection of co-operatively repressed promoters indicates a role for polycomb group complexes. Targets include genes involved in apoptosis, cell migration and B-cell differentiation, and show a highly significant but subtle alteration in genes involved in mitosis. In order to assess the relevance of the BL31 system to LCLs, we analysed the transcriptome of a set of EBNA3B knockout (3BKO) LCLs. Around a third of the genes whose expression level in LCLs was altered in the absence of EBNA3B were also altered in 3BKO-BL31 cell lines.Among these are TERT and TCL1A, implying that EBV-induced changes in the expression of these genes are not required for B-cell transformation. We also identify 26 genes that require both EBNA3A and EBNA3B for their regulation in LCLs. Together, this shows the complexity of the interaction between EBV and its host, whereby multiple EBNA3 proteins co-operate to modulate the behaviour of the host cell.