Replicating Marek's disease virus (MDV) serotype 2 DNA with inserted MDV serotype 1 DNA sequences in a Marek's disease lymphoblastoid cell line MSB1-41C

Targeted deletion of glycoprotein B gene by CRISPR/Cas9 nuclease inhibits Gallid herpesvirus type 3 in dually-infected Marek's disease virus-transformed lymphoblastoid cell line MSB-1

Marek’s disease virus (MDV) is a member of the genus Mardivirus in the subfamily Alphaherpesvirinae. There are three different serotypes of MDV designated as MDV-1 (Gallid herpesvirus type 2), MDV-2 (Gallid herpesvirus type 3), and MDV-3 (Meleagrid herpesvirus 1, herpesvirus of turkeys, HVT). MDV-1 is the only serotype that induces Marek’s disease (MD), a lymphoproliferative disorder resulting in aggressive T-cell lymphomas and paralytic symptoms. In the lymphomas and lymphoblastoid cell lines (LCL) derived from them, MDV establishes latent infection with limited viral gene expression. The latent viral genome in LCL can be activated by co-cultivation with chicken embryo fibroblast (CEF) monolayers. MSB-1, one of the first MDV-transformed LCL established from the splenic lymphoma, is distinct in harboring both the oncogenic MDV-1 and non-oncogenic MDV-2 viruses. Following the successful application of CRISPR/Cas9 editing approach for precise knockdown of the MDV-1 genes in LCL, we describe here the targeted deletion of MDV-2 glycoprotein B (gB) in MSB-1 cells. Due to the essential nature of gB for infectivity, the production of MDV-2 plaques on CEF was completely abolished in the MDV-2-gB-deleted MSB-1 cells. Our study has demonstrated that the CRISPR/Cas9 system can be used for targeted inactivation of the co-infecting MDV-2 without affecting the MDV-1 in the MSB-1 cell line. Successful inactivation of MDV-2 demonstrated here also points toward the possibility of using targeted gene editing as an antiviral strategy against pathogenic MDV-1 and other viruses infecting chickens.

IMPORTANCE Marek’s disease (MD) is a lymphoproliferative disease of chickens characterized by rapid-onset lymphomas in multiple organs and by infiltration into peripheral nerves, causing paralysis. Lymphoblastoid cell lines (LCL) derived from MD lymphomas have served as valuable resources to improve understanding of distinct aspects of virus-host interactions in transformed cells including transformation, latency, and reactivation. MDV-transformed LCL MSB-1, derived from spleen lymphoma induced by the BC-1 strain of MDV, has a unique feature of harboring an additional non-pathogenic MDV-2 strain HPRS-24. By targeted deletion of essential gene glycoprotein B from the MDV-2 genome within the MSB-1 cells, we demonstrated the total inhibition of MDV-2 virus replication on co-cultivated CEF, with no effect on MDV-1 replication. The identified viral genes critical for reactivation/inhibition of viruses will be useful as targets for development of de novo disease resistance in chickens to avian pathogens.

Epigenetic Silencing of MicroRNA-126 Promotes Cell Growth in Marek's Disease

During latency, herpesvirus infection results in the establishment of a dormant state in which a restricted set of viral genes are expressed. Together with alterations of the viral genome, several host genes undergo epigenetic silencing during latency. These epigenetic dysregulations of cellular genes might be involved in the development of cancer. In this context, Gallid alphaherpesvirus 2 (GaHV-2), causing Marek’s disease (MD) in susceptible chicken, was shown to impair the expression of several cellular microRNAs (miRNAs). We decided to focus on gga-miR-126, a host miRNA considered a tumor suppressor through signaling pathways controlling cell proliferation. Our objectives were to analyze the cause and the impact of miR-126 silencing during GaHV-2 infection. This cellular miRNA was found to be repressed at crucial steps of the viral infection. In order to determine whether miR-126 low expression level was associated with specific epigenetic signatures, DNA methylation patterns were established in the miR-126 gene promoter. Repression was associated with hypermethylation at a CpG island located in the miR-126 host gene epidermal growth factor like-7 (EGFL-7). A strategy was developed to conditionally overexpress miR-126 and control miRNAs in transformed CD4+ T cells propagated from Marek’s disease (MD) lymphoma. This functional assay showed that miR-126 restoration specifically diminishes cell proliferation. We identified CT10 regulator of kinase (CRK), an adaptor protein dysregulated in several human malignancies, as a candidate target gene. Indeed, CRK protein levels were markedly reduced by the miR-126 restoration.

Linc-GALMD1 Regulates Viral Gene Expression in the Chicken

A rapidly increasing number of reports on dysregulated long intergenic non-coding RNA (lincRNA) expression across numerous types of cancers indicates that aberrant lincRNA expression may be a major contributor to tumorigenesis. Marek’s disease (MD) is a T cell lymphoma of chickens induced by Marek’s disease virus (MDV). Although we have investigated the roles of lincRNAs in bursa tissue of MDV-infected chickens in previous studies, the molecular mechanisms of lincRNA functions in T cells remain poorly understood. In the present study, Linc-GALMD1 was identified from CD4+ T cells and MSB1 cells, and its expression was significantly downregulated in MD-resistant line of birds in response to MDV challenge. Furthermore, loss-of-function experiments indicated that linc-GALMD1 significantly affected the expression of 290 genes in trans. Through integrated analysis of differentially expressed genes (DEGs) induced by MDV and linc-GALMD1, we found that IGLL1 gene expression levels had a positive correlation with the degree of MD infection and could potentially serve as an indicator for clinical diagnosis of MD. Moreover, an interaction between MDV and linc-GALMD1 was also observed. Accordingly, chicken embryonic fibroblast cells were inoculated with MDV with and without the linc-GALMD1 knockdown, and the data showed that linc-GALMD1 could repress MDV gene expression during the course of MDV infection. These findings uncovered a role of linc-GALMD1 as a viral gene regulator and suggested a function of linc-GALMD1 contributing to tumor suppression by coordinating expression of MDV genes and tumor-related genes and regulating immune responses to MDV infection.

Alterations in cellular and viral microRNA and cellular gene expression in Marek's disease virus-transformed T-cell lines treated with sodium butyrate

A shared feature of herpesviruses is their ability to enter a latent state following an initially lytic infection. Marek's disease virus serotype 1 (MDV-1) is an oncogenic avian herpesvirus. Small RNA profiling studies have suggested that microRNAs (miRNAs) are involved in viral latency. Sodium butyrate treatment is known to induce herpesvirus reactivation. The present study was undertaken to determine transcriptome and miRNome changes induced by sodium butyrate in 2 MDV-transformed cell lines, RP2 and CU115. In the first 24 h post-treatment, microarray analysis of transcriptional changes in cell lines RP2 and CU115 identified 137 and 114 differentially expressed genes, respectively. Small RNA deep-sequencing analysis identified 17 cellular miRNAs that were differentially expressed. The expression of MDV-encoded miRNAs was also altered upon treatment. Many of the genes and miRNAs that are differentially expressed are involved in regulation of the cell cycle, mitosis, DNA metabolism, and lymphocyte differentiation.

Growth inhibitive effect of betulinic acid combined with tripterine on MSB-1 cells and its mechanism

Marek's disease (MD), a highly infectious lymphoproliferative disease in chickens, is caused by a cell-associated oncogenic herpesvirus, Marek's disease virus (MDV). MSB-1 is a MD-derived lymphoblastoid cell line and can induce tumors when inoculated into susceptible chickens. Betulinic acid, which is present as one of the major effective components in many traditional Chinese medicines, has recently been reported to inhibit growth of cancer cells and employed as a potential anticancer agent. Tripterine, a major active compound extracted from the Chinese herb Tripterygium wilfordii Hook F, has now also shown anti-tumor activities in various cancers. The aim of this study was to investigate the synergistic growth-inhibitive effect of betulinic acid combined with tripterine on MSB-1 cells and its mechanism. Viability of MSB-1 cells was assessed by 3-(4,5-dimethylthiazol-2-y1)-2,5-diphenyltetrazolium bromide (MTT) method. Cell apoptotic analysis was performed by fluorescence detection. NF-κB transcription activity was detected by measuring luciferase activity. Western blotting was used to analyze the expression of p65, IκB and Meq. Our results showed that the proliferation in the combination group was significantly decreased as compared with that of monotherapy using betulinic acid or tripterine, accompanied by an induction of apoptosis, inhibition of NF-κB transcriptional activity and its targeting oncogenic gene Meq. The results suggest that the combination of betulinic acid and tripterine at lower concentration may produce a synergistic inhibitive effect on MSB-1 cells that warrants further investigation for its potential clinical applications.

Genome-wide mapping of DNase I hypersensitive sites and association analysis with gene expression in MSB1 cells

DNase I hypersensitive sites (DHSs) mark diverse classes of cis-regulatory regions, such as promoters and enhancers. MSB-1 derived from chicken Marek's disease (MD) lymphomas is an MDV-transformed CD4+ T-cell line for MD study. Previously, DNase I HS sites were studied mainly in human cell types for mammalian. To capture the regulatory elements specific to MSB1 cells and explore the molecular mechanisms of T-cell transformation caused by MDV in MD, we generated high-quality of DHSs map and gene expression profile for functional analysis in MSB1 cell line. The total of 21,724 significant peaks of DHSs was identified from around 40 million short reads. DHSs distribution varied between chromosomes and they preferred to enrich in the gene-rich chromosomes. More interesting, DHSs enrichments appeared to be scarce on regions abundant in CpG islands. Besides, we integrated DHSs into the gene expression data and found that DHSs tended to enrich on high expressed genes throughout whole gene regions while DHSs did not show significant changes for low and silent expressed genes. Furthermore, the correlation of DHSs with lincRNAs expression was also calculated and it implied that enhancer-associated lincRNAs probably originated from enhancer-like regions of DHSs. Together, our results indicated that DNase I HS sites highly correlate with active genes expression in MSB1 cells, suggesting DHSs can be considered as markers to identify the cis-regulatory elements associated with chicken Marek's disease.

Small RNA cloning and sequencing strategy affects host and viral microRNA expression signatures

The establishment of the microRNA (miRNA) expression signatures is the basic element to investigate the role played by these regulatory molecules in the biology of an organism. Marek's disease virus 1 (MDV-1) is an avian herpesvirus that naturally infects chicken and induces T cells lymphomas. During latency, MDV-1, like other herpesviruses, expresses a limited subset of transcripts. These include three miRNA clusters. Several studies identified the expression of virus and host encoded miRNAs from MDV-1 infected cell cultures and chickens. But a high discrepancy was observed when miRNA cloning frequencies obtained from different cloning and sequencing protocols were compared. Thus, we analyzed the effect of small RNA library preparation and sequencing on the miRNA frequencies obtained from the same RNA samples collected during MDV-1 infection of chicken at different steps of the oncoviral pathogenesis. Qualitative and quantitative variations were found in the data, depending on the strategy used. One of the mature miRNA derived from the latency-associated-transcript (LAT), mdv1-miR-M7-5p, showed the highest variation. Its cloning frequency was 50% of the viral miRNA counts when a small scale sequencing approach was used. Its frequency was 100 times less abundant when determined through the deep sequencing approach. Northern blot analysis showed a better correlation with the miRNA frequencies found by the small scale sequencing approach. By analyzing the cellular miRNA repertoire, we also found a gap between the two sequencing approaches. Collectively, our study indicates that next-generation sequencing data considered alone are limited for assessing the absolute copy number of transcripts. Thus, the quantification of small RNA should be addressed by compiling data obtained by using different techniques such as microarrays, qRT-PCR and NB analysis in support of high throughput sequencing data. These observations should be considered when miRNA variations are studied prior addressing functional studies.

Role of virus-encoded microRNAs in Avian viral diseases

With total dependence on the host cell, several viruses have adopted strategies to modulate the host cellular environment, including the modulation of microRNA (miRNA) pathway through virus-encoded miRNAs. Several avian viruses, mostly herpesviruses, have been shown to encode a number of novel miRNAs. These include the highly oncogenic Marek’s disease virus-1 (26 miRNAs), avirulent Marek’s disease virus-2 (36 miRNAs), herpesvirus of turkeys (28 miRNAs), infectious laryngotracheitis virus (10 miRNAs), duck enteritis virus (33 miRNAs) and avian leukosis virus (2 miRNAs). Despite the closer antigenic and phylogenetic relationship among some of the herpesviruses, miRNAs encoded by different viruses showed no sequence conservation, although locations of some of the miRNAs were conserved within the repeat regions of the genomes. However, some of the virus-encoded miRNAs showed significant sequence homology with host miRNAs demonstrating their ability to serve as functional orthologs. For example, mdv1-miR-M4-5p, a functional ortholog of gga-miR-155, is critical for the oncogenicity of Marek’s disease virus. Additionally, we also describe the potential association of the recently described avian leukosis virus subgroup J encoded E (XSR) miRNA in the induction of myeloid tumors in certain genetically-distinct chicken lines. In this review, we describe the advances in our understanding on the role of virus-encoded miRNAs in avian diseases.

Analysis of the mRNA targetome of microRNAs expressed by Marek's disease virus

Marek’s disease virus 1 (MDV-1), an oncogenic α-herpesvirus that induces T-cell lymphomas in chickens, serves as model system to study transformation by lymphotropic herpesviruses. Like the oncogenic human γ-herpesviruses Kaposi’s sarcoma-associated herpesvirus (KSHV) and Epstein-Barr virus (EBV), MDV-1 encodes several viral microRNAs (miRNAs). One MDV-1 miRNA, miR-M4, shares the same “seed” targeting sequence with both a KSHV miRNA, miR-K11, and cellular miR-155. Importantly, miR-M4 plays a critical role in T-cell transformation by MDV-1, while miR-K11 and cellular miR-155 are thought to play key roles in B-cell transformation by KSHV and EBV, respectively. Here, we present an analysis of the mRNAs targeted by viral miRNAs expressed in the chicken T-cell line MSB1, which is naturally coinfected with MDV-1 and the related nonpathogenic virus MDV-2. Our analysis identified >1,000 endogenous mRNAs targeted by miRNAs encoded by each virus, many of which are targeted by both MDV-1 and MDV-2 miRNAs. We present a functional analysis of an MDV-1 gene, RLORF8, targeted by four MDV-1 miRNAs and a cellular gene, encoding interleukin-18 (IL-18) and targeted by both MDV-1 and MDV-2 miRNAs, and show that ectopic expression of either protein in a form resistant to miRNA inhibition results in inhibition of cell proliferation. Finally, we present a restricted list of 9 genes targeted by not only MDV-1 miR-M4 but also KSHV miR-K11 and human miR-155. Given the critical role played by miR-155 seed family members in lymphomagenesis in humans and chickens, these mRNA targets may contain genes whose inhibition plays a conserved role in herpesvirus transformation.

Herpesviruses cause lymphomas in both humans and chickens, and in both cases, evidence indicates that virally encoded miRNAs, or virally subverted cellular miRNAs, belonging to the miR-155 seed family, play a critical role in this process. However, because each miRNA regulates numerous cellular mRNAs species, it has been difficult to elucidate which miRNA targets are important. Given the evolutionary distance between chickens and humans and the observation that miR-155 is nevertheless highly conserved in both species, we reasoned that the identification of shared miR-155 targets might shed light on this process. Here, we present an analysis of the mRNAs targeted by miRNAs encoded by the oncogenic avian herpesvirus MDV-1 in transformed chicken T cells, including a short list of mRNAs that are also targeted by miR-155 seed family miRNAs in EBV- or KSHV-transformed human B cells, and present an initial functional analysis of some of these miRNA targets.

The oncogenic microRNA OncomiR-21 overexpressed during Marek's disease lymphomagenesis is transactivated by the viral oncoprotein Meq

Gallid herpesvirus 2 (GaHV-2) is an oncogenic herpesvirus that causes T lymphoma in chicken. GaHV-2 encodes a basic leucine zipper (bZIP) protein of the AP-1 family, Meq. Upon formation of homo- or heterodimers with c-Jun, Meq may modulate the expression of viral and cellular genes involved in lymphomagenesis. GaHV-2 also encodes viral microRNAs (miRNAs) involved in latency and apoptosis escape. However, little is known about cellular miRNA deregulation during the development of GaHV-2-associated lymphoma. We determined the cellular miRNA expression profiles of chickens infected with a very virulent strain (RB-1B) or a vaccine strain (CVI988) or noninfected. Among the most deregulated cellular miRNAs, we focused our efforts on gga-miR-21, which is upregulated during GaHV-2 infection. We mapped the gga-miR-21 promoter to the 10th intron of the TMEM49 gene and found it to be driven by AP-1- and Ets-responsive elements. We show here that the viral oncoprotein Meq binds to this promoter, thereby transactivating gga-miR-21 expression. We confirmed that this miRNA targets chicken programmed death cell 4 (PDCD4) and promotes tumor cell growth and apoptosis escape. Finally, gga-miR-21 was overexpressed only during infection with a very virulent strain (RB-1B) and not during infection with a nononcogenic strain (CVI988), providing further evidence for its role in GaHV-2 lymphomagenesis. Our data therefore suggest an additional role for Meq in GaHV-2-mediated lymphomagenesis through the induction of miR-21 expression.

A p53-dependent promoter associated with polymorphic tandem repeats controls the expression of a viral transcript encoding clustered microRNAs

The tumor suppressor protein p53 plays a role in cellular responses to cancer-initiating events by regulating progress through the cell cycle. Several recent studies have shown that p53 transactivates expression of the members of the proapoptotic microRNA-34 family, which are underexpressed in several cancers. We demonstrate here that the latency-associated cluster of microRNAs (miRNA) encoded by an oncogenic herpesvirus, gallid herpesvirus 2 (GaHV-2), is a direct target of p53. Robust transcriptional activity was induced in three avian cell lines by a sequence mapping 600 base pairs (bp) upstream of the cluster of miRNAs. We found transcription start sites for the pri-miRNA transcript at the 3' end of this transcription-inducing sequence. The promoter has no consensus core promoter element, but is organized into a variable number of tandem repeats of 60-bp harboring p53-responsive elements (RE). The minimal functional construct consists of two tandem repeats. Mutagenesis to change the sequence of the p53 RE abolished transcriptional activity, whereas p53 induction enhanced mature miRNA expression. The identification of a viral miRNA promoter regulated by p53 is biologically significant, because all avirulent GaHV-2 strains described to date lack the corresponding regulatory sequence, whereas all virulent, very virulent, and hypervirulent strains possess at least two tandem repeats harboring the p53 RE.

MicroRNAs of Gallid and Meleagrid herpesviruses show generally conserved genomic locations and are virus-specific

Many herpesviruses, including Marek's disease viruses (MDV1 and MDV2), encode microRNAs. In this study, we report microRNAs of two related herpesviruses, infectious laryngotracheitis virus (ILTV) and herpesvirus of turkeys (HVT), as well as additional MDV2 microRNAs. The genome locations, but not microRNA sequences, are conserved among all four of these avian herpesviruses. Most are clustered in the repeats flanking the unique long region (I/TR(L)), except in ILTV which lacks these repeats. Two abundant ILTV microRNAs are antisense to the immediate early gene ICP4. A homologue of host microRNA, gga-miR-221, was found among the HVT microRNAs. Additionally, a cluster of HVT microRNAs was found in a region containing two locally duplicated segments, resulting in paralogous HVT microRNAs with 96-100% identity. The prevalence of microRNAs in the genomic repeat regions as well as in local repeats suggests the importance of genetic plasticity in herpesviruses for microRNA evolution and preservation of function.

MicroRNAs 221 and 222 target p27Kip1 in Marek's disease virus-transformed tumour cell line MSB-1

MicroRNAs (miRNAs) are a class of short RNAs that function as post-transcriptional suppressors of protein expression and are involved in a variety of biological processes, including oncogenesis. Several recent studies have implicated the involvement of miR-221 and miR-222 in tumorigenesis as these miRNAs are upregulated in a number of cancers and affect the expression of cell cycle regulatory proteins such as the cyclin-dependent kinase (cdk) inhibitor p27(Kip1). Marek's disease virus (MDV) is a highly oncogenic herpesvirus that affects poultry, causing acute neoplastic disease with lymphomatous lesions in several organs. MDV-encoded oncogenes such as Meq are directly implicated in the neoplastic transformation of T cells and have been well studied. More recently, however, the involvement of both host and virus-encoded miRNAs in the induction of MD lymphomas is being increasingly recognized. We analysed the miRNA expression profiles in the MDV-transformed lymphoblastoid cell line MSB-1 and found that endogenous miRNAs miR-221 and miR-222 were significantly upregulated. Demonstration of the conserved binding sites for these miRNAs in the chicken p27(Kip1) 3'-untranslated region sequence and the repression of luciferase activity of reporter constructs indicated that miR-221 and miR-222 target p27(Kip1) in these cells. We also found that overexpression of miR-221 and miR-222 decreased p27(Kip1) levels and that treatment with retrovirally expressed antagomiRs partially alleviated this suppression. These data show that an oncogenic herpesvirus, as in the case of many cancers, can exploit the miRNA machinery for suppressing cell cycle regulatory molecules such as p27(Kip1) in the induction and progression of T-cell lymphomas.

MicroRNA profile of Marek's disease virus-transformed T-cell line MSB-1: predominance of virus-encoded microRNAs

Research over the last few years has demonstrated the increasing role of microRNAs (miRNAs) as major regulators of gene expression in diverse cellular processes and diseases. Several viruses, particularly herpesviruses, also use the miRNA pathway of gene regulation by encoding their own miRNAs. Marek's disease (MD) is a widespread lymphomatous neoplastic disease of poultry caused by the highly contagious Marek's disease virus type 1 (MDV-1). Recent studies using virus-infected chicken embryo fibroblasts have identified at least eight miRNAs that map to the R(L)/R(S) region of the MDV genome. Since MDV is a lymphotropic virus that induces T-cell lymphomas, analysis of the miRNA profile in T-cell lymphoma would be more relevant for examining their role in oncogenesis. We determined the viral and host miRNAs expressed in MSB-1, a lymphoblastoid cell line established from an MDV-induced lymphoma of the spleen. In this paper, we report the identification of 13 MDV-1-encoded miRNAs (12 by direct cloning and 1 by Northern blotting) from MSB-1 cells. These miRNAs, five of which are novel MDV-1 miRNAs, map to the Meq and latency-associated transcript regions of the MDV genome. Furthermore, we show that miRNAs encoded by MDV-1 and the coinfected MDV-2 accounted for >60% of the 5,099 sequences of the MSB-1 "miRNAome." Several chicken miRNAs, some of which are known to be associated with cancer, were also cloned from MSB-1 cells. High levels of expression of MDV-1-encoded miRNAs and potentially oncogenic host miRNAs suggest that miRNAs may have major roles in MDV pathogenesis and neoplastic transformation.

Marek's disease virus type 2 (MDV-2)-encoded microRNAs show no sequence conservation with those encoded by MDV-1

MicroRNAs (miRNAs) are increasingly being recognized as major regulators of gene expression in many organisms, including viruses. Among viruses, members of the family Herpesviridae account for the majority of the currently known virus-encoded miRNAs. The highly oncogenic Marek's disease virus type 1 (MDV-1), an avian herpesvirus, has recently been shown to encode eight miRNAs clustered in the MEQ and LAT regions of the viral genome. The genus Mardivirus, to which MDV-1 belongs, also includes the nononcogenic but antigenically related MDV-2. As MDV-1 and MDV-2 are evolutionarily very close, we sought to determine if MDV-2 also encodes miRNAs. For this, we cloned, sequenced, and analyzed a library of small RNAs from the lymphoblastoid cell line MSB-1, previously shown to be coinfected with both MDV-1 and MDV-2. Among the 5,099 small RNA sequences determined from the library, we identified 17 novel MDV-2-specific miRNAs. Out of these, 16 were clustered in a 4.2-kb long repeat region that encodes R-LORF2 to R-LORF5. The single miRNA outside the cluster was located in the short repeat region, within the C-terminal region of the ICP4 homolog. The expression of these miRNAs in MSB-1 cells and infected chicken embryo fibroblasts was further confirmed by Northern blotting analysis. The identification of miRNA clusters within the repeat regions of MDV-2 demonstrates conservation of the relative genomic positions of miRNA clusters in MDV-1 and MDV-2, despite the lack of sequence homology among the miRNAs of the two viruses. The identification of these novel miRNAs adds to the growing list of virus-encoded miRNAs.