Novel microRNAs (miRNAs) encoded by herpesvirus of Turkeys: evidence of miRNA evolution by duplication

Critical roles of non-coding RNAs in lifecycle and biology of Marek's disease herpesvirus

Over the past two decades, numerous non-coding RNAs (ncRNAs) have been identified in different biological systems including virology, especially in large DNA viruses such as herpesviruses. As a representative oncogenic alphaherpesvirus, Marek's disease virus (MDV) causes an important immunosuppressive and rapid-onset neoplastic disease of poultry, namely Marek's disease (MD). Vaccinations can efficiently prevent the onset of MD lymphomas and other clinical disease, often heralded as the first successful example of vaccination-based control of cancer. MDV infection is also an excellent model for research into virally-induced tumorigenesis. Recently, great progress has been made in understanding the functions of ncRNAs in MD biology. Herein, we give a review of the discovery and identification of MDV-encoded viral miRNAs, focusing on the genomics, expression profiles, and emerging critical roles of MDV-1 miRNAs as oncogenic miRNAs (oncomiRs) or tumor suppressor genes involved in the induction of MD lymphomas. We also described the involvements of host cellular miRNAs, lincRNAs, and circRNAs participating in MDV life cycle, pathogenesis, and/or tumorigenesis. The prospects, strategies, and new techniques such as the CRISPR/Cas9-based gene editing applicable for further investigation into the ncRNA-mediated regulatory mechanisms in MDV pathogenesis/oncogenesis were also discussed, together with the possibilities of future studies on antiviral therapy and the development of new efficient MD vaccines.

Efficient Mutagenesis of Marek's Disease Virus-Encoded microRNAs Using a CRISPR/Cas9-Based Gene Editing System

The virus-encoded microRNAs (miRNAs) have been demonstrated to have important regulatory roles in herpesvirus biology, including virus replication, latency, pathogenesis and/or tumorigenesis. As an emerging efficient tool for gene editing, the clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system has been successfully applied in manipulating the genomes of large DNA viruses. Herein, utilizing the CRISPR/Cas9 system with a double-guide RNAs transfection/virus infection strategy, we have established a new platform for mutagenesis of viral miRNAs encoded by the Marek's disease virus serotype 1 (MDV-1), an oncogenic alphaherpesvirus that can induce rapid-onset T-cell lymphomas in chickens. A series of miRNA-knocked out (miR-KO) mutants with deletions of the Meq- or the mid-clustered miRNAs, namely RB-1B∆Meq-miRs, RB-1B∆M9-M2, RB-1B∆M4, RB-1B∆M9 and RB-1B∆M11, were generated from vvMDV strain RB-1B virus. Interestingly, mutagenesis of the targeted miRNAs showed changes in the in vitro virus growth kinetics, which is consistent with that of the in vivo proliferation curves of our previously reported GX0101 mutants produced by the bacterial artificial chromosome (BAC) clone and Rec E/T homologous recombination techniques. Our data demonstrate that the CRISPR/Cas9-based gene editing is a simple, efficient and relatively nondisruptive approach for manipulating the small non-coding genes from the genome of herpesvirus and will undoubtedly contribute significantly to the future progress in herpesvirus biology.

The Roles of MicroRNAs (miRNAs) in Avian Response to Viral Infection and Pathogenesis of Avian Immunosuppressive Diseases

MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in the regulation of various biological processes including cell development and differentiation, apoptosis, tumorigenesis, immunoregulation and viral infections. Avian immunosuppressive diseases refer to those avian diseases caused by pathogens that target and damage the immune organs or cells of the host, increasing susceptibility to other microbial infections and the risk of failure in subsequent vaccination against other diseases. As such, once a disease with an immunosuppressive feature occurs in flocks, it would be difficult for the stakeholders to have an optimal economic income. Infectious bursal disease (IBD), avian leukemia (AL), Marek’s disease (MD), chicken infectious anemia (CIA), reticuloendotheliosis (RE) and avian reovirus infection are on the top list of commonly-seen avian diseases with a feature of immunosuppression, posing an unmeasurable threat to the poultry industry across the globe. Understanding the pathogenesis of avian immunosuppressive disease is the basis for disease prevention and control. miRNAs have been shown to be involved in host response to pathogenic infections in chickens, including regulation of immunity, tumorigenesis, cell proliferation and viral replication. Here we summarize current knowledge on the roles of miRNAs in avian response to viral infection and pathogenesis of avian immunosuppressive diseases, in particular, MD, AL, IBD and RE.

Virulent duck enteritis virus infected DEF cells generate a unique pattern of viral microRNAs and a novel set of host microRNAs

Background

Duck enteritis virus (DEV) belongs to the family Herpesviridae and is an important epornitic agent that causes economic losses in the waterfowl industry. The Chinese virulent (CHv) and attenuate vaccines (VAC) are two different pathogenic DEV strains. MicroRNAs (miRNAs) are a class of non-coding RNAs that regulate gene expression in viral infection. Nonetheless, there is little information on virulent duck enteritis virus (DEV)-encoded miRNAs.

Results

Using high-throughput sequencing, we identified 39 mature viral miRNAs from CHv-infected duck embryo fibroblasts cells. Compared with the reported 33 VAC-encoded miRNAs, only 13 miRNA sequences and 22 “seed sequences” of miRNA were identical, and 8 novel viral miRNAs were detected and confirmed by stem-loop RT-qPCR in this study. Using RNAhybrid and PITA software, 38 CHv-encoded miRNAs were predicted to target 41 viral genes and formed a complex regulatory network. Dual luciferase reporter assay (DLRA) confirmed that viral dev-miR-D8-3p can directly target the 3’-UTR of CHv US1 gene (p < 0.05). Gene Ontology analysis on host target genes of viral miRNAs were mainly involved in biological regulation, cellular and metabolic processes. In addition, 598 novel duck-encoded miRNAs were detected in this study. Thirty-eight host miRNAs showed significant differential expression after CHv infection: 13 miRNAs were up-regulated, and 25 miRNAs were down-regulated, which may affect viral replication in the host cell.

Conclusions

These data suggested that CHv encoded a different set of microRNAs and formed a unique regulatory network compared with VAC. This is the first report of DEF miRNAs expression profile and an analysis of these miRNAs regulatory mechanisms during DEV infection. These data provide a basis for further exploring miRNA regulatory roles in the pathogenesis of DEV infection and contribute to the understanding of the CHv-host interaction at the miRNA level.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1468-2) contains supplementary material, which is available to authorized users.

A phylogenomic analysis of Marek's disease virus reveals independent paths to virulence in Eurasia and North America

Virulence determines the impact a pathogen has on the fitness of its host, yet current understanding of the evolutionary origins and causes of virulence of many pathogens is surprisingly incomplete. Here, we explore the evolution of Marek's disease virus (MDV), a herpesvirus commonly afflicting chickens and rarely other avian species. The history of MDV in the 20th century represents an important case study in the evolution of virulence. The severity of MDV infection in chickens has been rising steadily since the adoption of intensive farming techniques and vaccination programs in the 1950s and 1970s, respectively. It has remained uncertain, however, which of these factors is causally more responsible for the observed increase in virulence of circulating viruses. We conducted a phylogenomic study to understand the evolution of MDV in the context of dramatic changes to poultry farming and disease control. Our analysis reveals evidence of geographical structuring of MDV strains, with reconstructions supporting the emergence of virulent viruses independently in North America and Eurasia. Of note, the emergence of virulent viruses appears to coincide approximately with the introduction of comprehensive vaccination on both continents. The time‐dated phylogeny also indicated that MDV has a mean evolutionary rate of ~1.6 × 10⁻⁵ substitutions per site per year. An examination of gene‐linked mutations did not identify a strong association between mutational variation and virulence phenotypes, indicating that MDV may evolve readily and rapidly under strong selective pressures and that multiple genotypic pathways may underlie virulence adaptation in MDV.

In vivo expression patterns of microRNAs of Gallid herpesvirus 2 (GaHV-2) during the virus life cycle and development of Marek's disease lymphomas

In the past decade, a large number of microRNAs (miRNAs) have been identified in the viral genome of Gallid herpesvirus 2 (GaHV-2), which is historically known as Marek’s disease virus type 1. The biological role of most GaHV-2 miRNAs remains unclear. In the present study, we have performed an overall gene expression profile of GaHV-2 miRNAs during the virus life cycle at each phase of the developing disease, a highly contagious, lymphoproliferative disorder, and neoplastic immunosuppressive disease of poultry known as the Marek’s disease. According to their distinct in vivo expression patterns, the GaHV-2 miRNAs can be divided into three groups: 12 miRNAs in group I, including miR-M4-5p, displayed a typical expression pattern potentially correlated to the latent, late cytolytic, and/or the proliferative phases in the cycle of GaHV-2 pathogenesis; group II consisting of another 12 miRNAs with expression correlated to the early cytolytic and/or latent phases in GaHV-2’s life cycle; while the other two miRNAs in group III showed no identical expression features. Our findings may provide meaningful clues in the search for further potential functions of viral miRNAs in GaHV-2 biology.

Marek's disease virus-encoded analog of microRNA-155 activates the oncogene c-Myc by targeting LTBP1 and suppressing the TGF-β signaling pathway

Marek's disease virus (MDV) is a representative alpha herpes virus able to induce rapid-onset T-cell lymphoma in its natural host and regarded as an ideal model for the study of virus-induced tumorigenesis. Recent studies have shown that the mdv1-miR-M4-5p, a viral analog of cellular miR-155, is critical for MDV׳s oncogenicity. However, the precise mechanism whereby it was involved in MD lymphomagenesis remained unknown. We have presently identified the host mRNA targets of mdv1-miR-M4-5 and identified the latent TGF-β binding protein 1 (LTBP1) as a critical target for it. We found that during MDV infection, down-regulation of LTBP1 expression by mdv1-miR-M4-5p led to a significant decrease of the secretion and activation of TGF-β1, with suppression of TGF-β signaling and a significant activation of expression of c-Myc, a well-known oncogene which is critical for virus-induced tumorigenesis. Our findings reveal a novel and important mechanism of how mdv1-miR-M4-5p potentially contributes to MDV-induced tumorigenesis.

The significance of the individual Meq-clustered miRNAs of Marek's disease virus in oncogenesis

Marek's disease virus (MDV) is an important oncogenic alphaherpesvirus that induces rapid-onset T-cell lymphomas in its natural hosts. The Meq-clustered miRNAs encoded by MDV have been suggested to play potentially critical roles in the induction of lymphomas. Using the technique of bacterial artificial chromosome mutagenesis, we have presently constructed a series of specific miRNA-deleted mutants and demonstrate that these miRNAs are not essential for replication of MDV and have no effects on the early cytolytic or latent phases of the developing disease. However, compared to the parental GX0101, mortality of birds infected with the mutants GXΔmiR-M2, GXΔmiR-M3, GXΔmiR-M5, GXΔmiR-M9 and GXΔmiR-M12 was reduced from 100 % to 18 %, 30 %, 48 %, 24 % and 14 %, coupled with gross tumour incidence reduction from 28 % to 8 %, 4 %, 12 %, 8 % and 0 %, respectively. Our data confirm that except for mdv1-miR-M4, the other Meq-clustered miRNAs also play critical roles in MDV oncogenesis. Further work will be needed to elucidate the miRNA-mediated regulatory mechanisms that trigger the development of MD lymphomas.

Molecular characterization of the complete genome of falconid herpesvirus strain S-18

Falconid herpesvirus type 1 (FaHV-1) is the causative agent of falcon inclusion body disease, an acute, highly contagious disease of raptors. The complete nucleotide sequence of the genome of FaHV-1 has been determined using Illumina MiSeq sequencing. The genome is 204,054 nucleotides in length and has a class E organization. The genome encodes approximately 130 putative protein-coding genes, of which 70 are orthologs of conserved alphaherpesvirus and Mardivirus proteins. Three FaHV-1 genes (UL3.5, UL44.5 and CIRC) were identified that encode protein homologues unique to Mardivirus and Varicellovirus. The genome also encodes homologues to the Mardivirus genes LORF2, LORF3, LORF4, LORF5, SORF3 and SORF4. An opal mutation resulting in premature termination was identified in the FaHV-1 UL43 gene. Phylogenetically, FaHV-1 resides in a monophyletic group with the other Mardiviruses but, along with anatid herpesvirus 1, represents a more distant divergence from the rest of the Mardivirus genus.

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.

Virus-encoded miR-155 ortholog is an important potential regulator but not essential for the development of lymphomas induced by very virulent Marek's disease virus

The microRNA (miRNA) mdv1-miR-M4, a functional miR-155 ortholog encoded by oncogenic Marek's disease virus (MDV), has previously been suggested to be involved in MDV pathogenesis. Using the technique of bacterial artificial chromosome mutagenesis, we have presently evaluated the potential role of mdv1-miR-M4 in the oncogenesis of the very virulent (vv) MDV strain GX0101. Unexpectedly, deletions of the Meq-cluster or mdv1-miR-M4 alone from the viral genome strongly decreased rather than abolished its oncogenicity. Compared to GX0101, mortalities of mutants GXΔmiR-M4 and GXΔMeq-miRs were reduced from 100% to 18% and 4%, coupled with the gross tumor incidence reduction from 28% to 22% and 8%, respectively. Our data suggests that the mdv1-miR-M4 is possibly an important regulator in the development of Marek's disease (MD) lymphomas but is not essential for the oncogenicity of vvMDV. In addition, some of the other Meq-clustered miRNAs may also play potentially critical roles in vvMDV induction of lymphomas.

The interplay between MDV and HVT affects viral miRNa expression

It is well established that herpesviruses encode numerous microRNAs (miRNAs) and that these virally encoded small RNAs play multiple roles in infection. The present study was undertaken to determine how co-infection of a pathogenic MDV serotype one (MDV1) strain (MD5) and a vaccine strain (herpesvirus of turkeys [HVT]) alters viral miRNA expression in vivo. We first used small RNA deep sequencing to identify MDV1-encoded miRNAs that are expressed in tumorigenic spleens of MDV1-infected birds. The expression patterns of these miRNAs were then further assessed at an early time point (7 days postinfection [dpi]) and a late time point (42 dpi) in birds with and without HVT vaccination using real-time PCR (RT-PCR). Additionally, the effect of MDV1 co-infection on HVT-encoded miRNAs was determined using RT-PCR. A diverse population of miRNAs was expressed in MDV-induced tumorigenic spleens at 42 dpi, with 18 of the 26 known mature miRNAs represented. Of these, both mdv1-miR-M4-5p and mdv1-miR-M2-3p were the most highly expressed miRNAs. RT-PCR analysis further revealed that nine MDV miRNAs were differentially expressed between 7 dpi and 42 dpi infected spleens. At 7 dpi, three miRNAs were differentially expressed between the spleens of birds co-infected with HVT and MD5 compared with birds singly infected with MD5, whereas at 42 dpi, nine miRNAs were differentially expressed. At 7 dpi, the expression of seven HVT-encoded miRNAs was affected in the spleens of co-infected birds compared with birds only receiving the HVT vaccine. At 42 dpi, six HVT-encoded miRNAs were differentially expressed between the two groups. Target prediction analysis suggests that these differentially expressed viral miRNAs are involved in regulating several cellular processes, including cell proliferation and the adaptive immune response.

Novel microRNAs encoded by duck enteritis virus

Duck enteritis virus (DEV) is an important herpesvirus pathogen associated with acute, highly contagious lethal disease in waterfowls. Using a deep sequencing approach on RNA from infected chicken embryo fibroblast cultures, we identified several novel DEV-encoded micro (mi)RNAs. Unlike most mardivirus-encoded miRNAs, DEV-encoded miRNAs mapped mostly to the unique long region of the genome. The precursors of DEV miR-D18 and miR-D19 overlapped with each other, suggesting similarities to miRNA-offset RNAs, although only the DEV-miR-D18-3p was functional in reporter assays. Identification of these novel miRNAs will add to the growing list of virus-encoded miRNAs enabling the exploration of their roles in pathogenesis.

Roles of avian herpesvirus microRNAs in infection, latency, and oncogenesis

MicroRNAs have been reported for the avian herpesviruses Marek's disease virus 1 (MDV1; oncogenic), Marek's disease virus 2 (MDV2; non-oncogenic), herpesvirus of turkeys (HVT), and infectious laryngotracheitis virus (ILTV). No obvious phylogenetic relationships exist among the avian herpesvirus microRNAs, but the general genomic locations of microRNA clusters are conserved, with microRNAs being located in the repeat regions of the genomes. In some cases, microRNAs are antisense to open reading frames. Among MDV1 field isolates with different virulence properties, microRNAs are highly conserved, and variations that have been observed lie in putative promoter regions. One cluster of MDV1 microRNAs lies upstream of the meq gene, and this cluster is more highly expressed in tumors caused by an extremely virulent MDV1 isolate compared to tumors caused by a less virulent isolate. Several of the avian herpesvirus microRNAs are orthologs of microRNAs in other species. For example, mdv1-miR-M4 shares a seed sequence with gga-miR-155 (also shared with Kaposi sarcoma herpesvirus (KSHV) kshv-miR-K12), mdv2-miR-M21 shares a seed with miR-29b, and hvt-miR-H14 shares a seed sequence with miR-221. Functional analyses of avian herpesvirus microRNAs include a variety of in vitro assays to demonstrate potential function as well as the use of mutants that can exploit the ability to assess phenotypes experimentally in the natural host. This article is part of a Special Issue entitled:MicroRNA's in viral gene regulation.

Noncoding RNPs of viral origin

Like their host cells, many viruses produce noncoding (nc)RNAs. These show diversity with respect to time of expression during viral infection, length and structure, protein-binding partners and relative abundance compared with their host-cell counterparts. Viruses, with their limited genomic capacity, presumably evolve or acquire ncRNAs only if they selectively enhance the viral life cycle or assist the virus in combating the host's response to infection. Despite much effort, identifying the functions of viral ncRNAs has been extremely challenging. Recent technical advances and enhanced understanding of host-cell ncRNAs promise accelerated insights into the RNA warfare mounted by this fascinating class of RNPs.

Expression profiles of microRNAs encoded by the oncogenic Marek's disease virus reveal two distinct expression patterns in vivo during different phases of disease

Marek's disease virus (MDV) is a long-recognized oncogenic herpesvirus, which induces lymphoma in its natural host that can be prevented by vaccination. MDV infection provides an excellent biological model for investigating the biology, genetics and immunology of viral oncogenesis. Recently discovered microRNAs (miRNAs) in the MDV genome have been suggested to have regulatory roles during MDV oncogenesis. We have examined the expression profiles of all 22 previously reported miRNAs encoded by MDV-1 in chickens artificially challenged with MDV-GX0101. We found that a subset of the miRNAs was differentially expressed during different phases of the developing disease. These miRNAs show early or late expression during disease progression, accompanied by obvious tissue-specific and differential expression patterns. This temporal and differential tissue distribution suggest that these miRNAs may perform different regulatory roles in switching from latency to lytic replication, immunosuppression, neoplastic transformation or other aspects of lymphoma formation. These reported in vivo expression profiles indicate the potentially functional MDV-1-encoded miRNAs that should be selected for further investigation of their functions in MDV oncogenesis.

miRNAs in the pathogenesis of oncogenic human viruses

Tumor viruses are a class of pathogens with well established roles in the development of malignant diseases. Numerous bodies of work have highlighted miRNAs (microRNAs) as critical regulators of tumor pathways and it is clear that the dysregulation of cellular miRNA expression can promote tumor formation. Tumor viruses encode their own miRNAs and/or manipulate the expression of cellular miRNAs to modulate their host cell environment, thereby facilitating their respective infection cycles. The modulation of these miRNA responsive pathways, however, often influences certain signal transduction cascades in ways that favor tumorigenesis. In this review, we discuss the roles of virally-encoded and virally-regulated cellular miRNAs in the respective viral life cycles and in virus associated pathogenesis.

MicroRNA-26a-mediated regulation of interleukin-2 expression in transformed avian lymphocyte lines

Background

Micro(mi)RNAs are a class of small non-coding RNAs that play critical roles in the induction of various cancers, including lymphomas induced by oncogenic viruses. While some of the miRNAs are oncogenic, miRNAs such as miR-26a are consistently downregulated in a number of cancers, demonstrating their potential tumor suppressor functions. Global miRNA expression profiles of a number of virus-transformed avian lymphoma cell lines have shown downregulation of gga-miR-26a expression, irrespective of molecular mechanisms of transformation or the viral aetiology. The neoplastic transformation of lymphocytes by many viruses accompanies high levels of proliferative responses, mostly mediated through cytokines such as IL-2. Chicken IL-2 can modulate T-cell proliferation and cytotoxicity in vitro and in vivo and dysregulation of IL-2 expression is observed in diseases such as leukaemia.

Results

The expression levels of gga-miR-26a in chicken lymphoma cells transformed by 3 distinct avian oncogenic viruses, viz Marek's disease virus (MDV), avian leukosis virus (ALV) and Reticuloendotheliosis virus (REV) were consistently downregulated compared to the levels in the normal lymphocytes. This downregulation of miR-26a regardless of the viral etiology and molecular mechanisms of transformation was consistent with the tumor suppressor role of this miRNA. Notwithstanding this well-established role in cancer, we demonstrate the additional role of this miRNA in directly targeting chicken IL-2 through reporter and biochemical assays. The downregulation of miR-26a can relieve the suppressive effect of this miRNA on IL-2 expression.

Conclusions

We show that miR-26a is globally downregulated in a number of avian lymphoma cells irrespective of the mechanisms of transformation, reiterating the highly conserved tumor suppressor function of this miRNA. However, with the potential for directly targeting chicken IL-2, the downregulation of miR-26a in these tumor cells could relieve the inhibitory effect on IL-2 expression assisting in the proliferative features of the transformed lymphocyte lines.