Next Article in Journal
Processing Phage Therapy Requests in a Brussels Military Hospital: Lessons Identified
Next Article in Special Issue
Gram-Positive Bacteria-Like DNA Binding Machineries Involved in Replication Initiation and Termination Mechanisms of Mimivirus
Previous Article in Journal
Early Multiphasic HBV Infection Initiation Kinetics Is Not Clone-Specific and Is Not Affected by Hepatitis D Virus (HDV) Infection
Previous Article in Special Issue
A Needle in A Haystack: Tracing Bivalve-Associated Viruses in High-Throughput Transcriptomic Data
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 11
Issue 3
10.3390/v11030264
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Transcriptional Landscape of Marek’s Disease Virus in Primary Chicken B Cells Reveals Novel Splice Variants and Genes

by
Luca D. Bertzbach
1,†,
Florian Pfaff
2,*,†,
Viktoria I. Pauker
3,
Ahmed M. Kheimar
1,4,
Dirk Höper
2,
Sonja Härtle
5,
Axel Karger
3,* and
Benedikt B. Kaufer
1,*
1
Institute of Virology, Freie Universität Berlin, Robert von Ostertag-Straße 7-13, 14163 Berlin, Germany
2
Institute of Diagnostic Virology, Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Südufer 10, 17493 Greifswald-Insel Riems, Germany
3
Institute of Molecular Virology and Cell Biology, Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Südufer 10, 17493 Greifswald-Insel Riems, Germany
4
Department of Poultry Diseases, Faculty of Veterinary Medicine, Sohag University, Sohag 82424, Egypt
5
Department of Veterinary Sciences, Institute for Animal Physiology, Ludwig-Maximilians-Universität München, 80539 Munich, Germany
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2019, 11(3), 264; https://doi.org/10.3390/v11030264
Submission received: 18 February 2019 / Revised: 12 March 2019 / Accepted: 13 March 2019 / Published: 16 March 2019
(This article belongs to the Special Issue Virus Bioinformatics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Marek’s disease virus (MDV) is an oncogenic alphaherpesvirus that infects chickens and poses a serious threat to poultry health. In infected animals, MDV efficiently replicates in B cells in various lymphoid organs. Despite many years of research, the viral transcriptome in primary target cells of MDV remained unknown. In this study, we uncovered the transcriptional landscape of the very virulent RB1B strain and the attenuated CVI988/Rispens vaccine strain in primary chicken B cells using high-throughput RNA-sequencing. Our data confirmed the expression of known genes, but also identified a novel spliced MDV gene in the unique short region of the genome. Furthermore, de novo transcriptome assembly revealed extensive splicing of viral genes resulting in coding and non-coding RNA transcripts. A novel splicing isoform of MDV UL15 could also be confirmed by mass spectrometry and RT-PCR. In addition, we could demonstrate that the associated transcriptional motifs are highly conserved and closely resembled those of the host transcriptional machinery. Taken together, our data allow a comprehensive re-annotation of the MDV genome with novel genes and splice variants that could be targeted in further research on MDV replication and tumorigenesis.

1. Introduction

Marek’s disease virus (MDV), also known as Gallid alphaherpesvirus 2, causes a deadly lymphoproliferative disease in chickens. Typical clinical symptoms include immunosuppression, paralysis and polyneuritis, acute brain edema, and lymphoma that develop as early as 3 weeks post infection [1,2]. MDV has a big economic impact on the poultry industry worldwide due to animal losses, reduced growth, decreased egg production, and cost of vaccination [3]. Vaccines are crucial for the protection against MDV, as very virulent strains can cause mortalities of up to 100% in susceptible unvaccinated chickens [4]. Live attenuated MDV vaccines such as the gold standard Rispens strain (CVI988) are highly effective in preventing tumor formation [3,5], but do not provide sterilizing immunity.
MDV infection is initiated by the inhalation of the virus from a contaminated environment. Macrophages and other phagocytic cells are thought to take up the virus and transfer it to lymphatic tissues, where B and T cells are infected [6]. B cells are efficiently infected during the initial lytic virus replication, whereas T cell subsets play key roles in MDV pathogenesis [7,8,9,10,11]. MDV establishes latency predominantly in CD4+ T cells, which can also transport the virus to the feather follicle epithelium. These cells efficiently produce the infectious virus and shed it into the environment [3]. Latently infected T cells can also be transformed, resulting in deadly lymphomas.
MDV has a 180 kilo base pairs class E genome consisting of a unique long (UL) and a unique short (US) sequence that are flanked by terminal (TRL and TRS) and internal (IRL and IRS) inverted repeat regions [12,13]. MDV encodes about 100 proteins that orchestrate the virus life cycle and/or contribute to pathogenesis [12,14]. Until now, analyses of the viral transcriptome has been limited to chicken fibroblasts that are not infected in chickens, ex vivo samples [15,16] and tumor cells [17]; however, the mRNA expression in the primary target cells of lytic replication in vivo remained elusive. This is mainly due to the short lifespan of B and T cells in culture and the low quantity of infected cells in lymphoid organs of chickens [7,18]. To overcome this obstacle, we recently developed an in vitro infection system for primary B and T cells that allows for a prolonged survival and efficient infection of these cells [19]. We used this system to analyze the MDV transcriptome in the most frequent lytically infected cell type in vivo, the B cells.
In this study, we performed next generation RNA-sequencing and protein profiling in primary B cells infected with the very virulent RB1B strain or the vaccine strain CVI988. Our data reveal that the coding capacity of the MDV genome is larger than expected. We identified novel MDV genes and splice variants, and confirmed them either on the protein level or by RT-PCR. This comprehensive approach provides novel insights into the transcriptome of MDV in the natural target cells and supply a basis for future research on MDV pathogenesis.

2. Materials and Methods

2.1. Ethics Statement

Valo specific-pathogen-free (SPF) chickens (VALO BioMedia GmbH, Osterholz-Scharmbeck, Germany) were housed for 6 to 11 weeks of age and humanely euthanized prior to the isolation of the bursa of Fabricius. The animal work was approved by the governmental agency, the Landesamt für Gesundheit und Soziales (LAGeSo) in Berlin, Germany (approval number T0245/14, approval date 23 October 2014).

2.2. Cells

Embryonated SPF Valo chicken eggs (VALO BioMedia GmbH,) were used for the preparation of chicken embryo cells (CEC). CEC were maintained in minimal essential medium (MEM, PAN Biotech; Aidenbach, Germany) supplemented with 1-10% fetal bovine serum (FBS) and penicillin/streptomycin as previously described [20]. B cells were obtained from the bursa of Fabricius by dissociation of the organ and subsequent isolation of the cells by density gradient centrifugation as previously described [21]. Briefly, the bursa of Fabricius was homogenized through a 40 μm cell filter to obtain a uniform single cell suspension. Suspension cells were carefully applied on a Biocoll separating solution (Biochrom; Berlin, Germany), centrifuged for 12 min at 650 × g with slow acceleration, and deactivated deceleration. Lymphocytes at the interphase were carefully transferred to a new tube, washed with PBS, and maintained in RPMI 1640 (PAN Biotech) supplemented with 10% FBS and penicillin [100 U/mL]/streptomycin [100 µg/mL] at 41 °C under a 5% CO2 atmosphere. B cells were activated using recombinant soluble chicken CD40 ligand (chCD40L) [22], which was expressed in HEK293 cells and purified using a Vivacell 250 ultrafiltration concentrator (Sartorius; Göttingen, Germany).

2.3. Viruses

All viruses were reconstituted by calcium phosphate transfection of CEC with purified bacterial artificial chromosome (BAC) DNA as previously described [23]. The very virulent RB1B strain and the vaccine strain CVI988 both express a green fluorescent protein (GFP) under the control of the early thymidine kinase promotor. The viruses were propagated on CEC for up to six passages, and infected cells were stored in liquid nitrogen until further use. All virus stocks were titrated on fresh CEC.

2.4. Infection of Primary Chicken B Cells

Primary chicken B cells were infected by co-cultivation with infected CEC due to the strictly cell-associated nature of MDV. One million CEC were infected with 30,000 plaque-forming units (PFU) of CVI988, RB1B, or mock infected. After 4 days, one million B cells were seeded on the infected 6-well-plates in the presence of CD40L for 16 h at 41 °C. All B cells were then carefully removed from the CEC monolayer, washed with phosphate buffered saline (PBS), and prepared for fluorescence-activated cell sorting.

2.5. Flow Cytometry

Viable bursal B cells were detected using the eFluor780 fixable viability dye at a dilution of 1:1000 (Affymetrix eBioscience; San Diego, CA, USA) as previously described [19]. Cells were sorted using a FACS Aria III sorter and the FACSDiva software (Becton Dickinson; Franklin Lakes, NJ, USA). For each sample, approximately 105 to 106 infected B cells were sorted at 4 °C and stored at −80 °C until further analysis. The purity of GFP+ sorted fractions was determined by FACS reanalysis and yielded 99.73% (±0.46 SD) for mock-infected cells, 95.33% (±1.29 SD) for RB1B infected cells and 97.47% (±1.33 SD) for CVI988 infected cells.

2.6. High-Throughput RNA-Sequencing

RNA was isolated from three independent experiments of CVI988 or RB1B infected chicken B cell cultures and sequenced as described [24]. Briefly, total RNA was extracted using TRIzol reagent (Life Technologies; Carlsbad, CA, USA) in combination with the RNeasy Mini Kit (Qiagen; Hilden, Germany) following the manufacturer’s instructions. Additionally, RNA was treated with DNase using the RNase-Free DNase Set (Qiagen). Subsequently, ERCC ExFold RNA Spike-In mix 1 (Invitrogen; Carlsbad, CA, USA) was added to the total RNA as an internal control and the polyadenylated (poly(A)) RNA fraction was extracted using the Dynabeads mRNA DIRECT Micro kit (Invitrogen). Whole transcriptome libraries were prepared using the Ion Total RNA-Seq Kit v2 (Life Technologies) following the manufacturer’s instructions. Quality and quantity of the nucleic acids was controlled at each step using the NanoDrop 1000 spectrophotometer (Peqlab) or Agilent 2100 Bioanalyzer (Agilent Technologies; Böblingen, Germany) in combination with appropriate chips, respectively. The resulting libraries were finally quantified using the KAPA Library Quantification Kit for Ion Torrent (Kapa Biosystems; Wilmington, MA, USA) on a CFX96 Real-Time PCR Detection System (BioRad Laboratories) and sequenced on an Ion S5XL system (Life Technologies) using the Ion 540 OT2 and Chip kit (Life Technologies).

2.7. Sequence and Bioinformatic Analyses

Reads from separate sequencing runs (technical replicates) of the same library (biological replicates) were combined and quality-trimmed using the 454 Sequencing System Software (v. 3.0; Roche; Mannheim, Germany) along with appropriate Ion Torrent specific adapter sequences. Each quality-trimmed data set was then mapped to a non-redundant version of the MDV reference NC_002229.3 [12] (only segments UL-IRL-IRS-US) using STAR (version 2.6.1a; [25]), running in basic two-pass mode. In this manuscript, the MDV genes were designated according to the current gene nomenclature used for the prototype alphaherpesvirus herpes simplex virus 1 (HSV-1) [26,27]. As the mRNA libraries were amplified for several rounds, the sequence duplicates were removed prior to de-novo assembly and coverage analysis. Therefore, the duplicated reads were marked and removed in each mapped dataset using Picard (version 2.18.20; http://broadinstitute.github.io/picard). Subsequently, the unique aligned reads were directionally sorted using samtools (version 1.9; [28]) and the sequence depth was deduced from each dataset using bedtools (version 2.15.0; [29]) and samtools. The directionally sorted reads were then used for de-novo assembly using the 454 Sequencing System Software (v. 3.0; Roche) running in “-cdna” mode. Assembly was done for each biological replicate separately and with the combined read data from all replicates. Deduplicated and directionally sorted reads were also used as basis for the coverage plots. In order to receive high quality splice junctions, all assembled “isotigs” (transcript variants) were then re-mapped to the non-redundant version of MDV reference NC_002229.3 [12] using STARlong (version 2.6.1a; [25]) in basic two-pass mode. Options were set to allow a single mapping isotig to yield a splice junction with a maximum intron length of 10,000 bp. Only positions corresponding to the high-quality splice junctions were then selected from the splice junctions of the initial mapping approach for further analysis. The frequency of spliced reads was calculated by dividing the number of reads with splice junction by the total number of reads at the respective donor site. The overall splice frequency at a single donor site for CVI988 and RB1B was then averaged from the individual splice frequencies in the replicates. Based on the deduced splice junctions, the positions up- and downstream of these were extracted and visualized using the R (version 3.4.1; [30]) package “ggseqlogo” (version 0.1; [31]) in combination with RStudio (version 1.0.153) in order to receive information on donor and acceptor motif sequences. Polyadenylation cleavage clusters were determined with ContextMap (version 2.7.9; [32]) in combination with bowtie2 (version 2.2.9; [33]) using the complete trimmed dataset as input and the “--polyA” parameter. The resulting poly(A) cleavage sites were then combined into clusters, as the exact position of mRNA cleavage downstream of a cleavage signal can be heterogeneous [34]. Starting with the first poly(A) cleavage site, all other sites on the same strand within a window of the next 30 nt were combined into a single cluster. The window was then moved to the next cleavage site that was not within the last cluster. This was repeated for all cleavage sites. To scan for enriched regulatory motifs within the three prime untranslated region (3′-UTR) of MDV transcripts, the sequences 50 nt up- and 20 nt downstream of the start position of each identified poly(A) cleavage cluster were extracted and analyzed using DREME as part of the MEME suite (version 4.9.0; [35]) using default settings. Subsequently, the identified enriched motif (AWTAAA) was searched in the non-redundant version of MDV reference NC_002229.3 [12] using FIMO (version 4.9.0; [36]) with default settings. Poly(A) cleavage clusters and regulatory motifs were then grouped into relevant pairs based on the FIMO p-value and the their distance, allowing a maximum distance of 50 nt. Differential gene expression between RB1B and the vaccine strain CVI988 was conducted using Salmon (version 0.12.0; [37]) in combination with DESeq2 (version 1.18.1; [38]) as described earlier [24]. All relevant MDV CDS sequences were used as transcript reference and genes with an adjusted p-value > 0.01 were considered significant. Potential phosphorylation sites in novel protein SORF6 were predicted using the NetPhos 3.1 Server [39] and its DNA-protein binding probability using the DNABIND server [40].

2.8. LC-MALDI TOF/TOF Mass Spectrometry

Infected FACS sorted primary chicken B cells and mock infected primary chicken B cells were lysed in batches of 1.5 × 106 cells using 150 µl of a lysis buffer containing 0.1 M DL-Dithothreitol (DTT) and 2% SDS in 0.1 M Tris-HCl (pH 8.0) at 99 °C for 5 min. Protein contents were determined by densitometry of Coomassie stained SDS gels [41,42]. After cell lysis, 20 µg aliquots were digested using the FASP protocol as described [43]. Samples were differentially labeled by dimethylation [44] using unlabeled and 13C-labeled formaldehyde, respectively, and subjected to nano-LC MALDI-TOF/TOF mass spectrometry as described previously [45]. Briefly, peptides were separated by nano reversed-phase liquid chromatography (EASY-nLC II, Bruker; Bremen, Germany), spotted to a MALDI target (Proteineer fcII, Bruker), and analyzed with an UltrafleXtreme MALDI-TOF/TOF mass spectrometer (Bruker) as described previously [46]. Peptide spectra were acquired in the m/z range 700 to 3.500 Da with a minimum signal-to-noise (S/N) ratio of 7. Proteins were identified with a Mascot server (version 2.4.1; Matrix Science Ltd; London, UK) and analyzed using ProteinScape software (version 3; Bruker). Oxidation of methionine, acetylation of protein N-termini, and dimethylation of lysine and peptide N-termini with either isotopomer were set as variable modifications, whereas the carbamydomethylation of cysteine residues was set as a fixed modification. Two independent experiments were performed with inverted labeling. As database for the protein identification with the MASCOT search engine (Matrix Science Ltd), the Gallus gallus proteome was downloaded from the ENSEMBL website [47] and the viral sequences were added to the FASTA file. Viral protein content was calculated in mol% using the exponentially modified protein abundance index (emPAI) [48]. To identify peptides covering the potential new splicing sites discovered by RNA-sequencing, a database with sequence fragments covering a 67 amino acid region centered on the splicing site was constructed and used for the database search with MASCOT.

2.9. Reverse Transcription and PCR over Splice Junctions

RNA was isolated as described above. cDNA was synthesized after DNAse treatment (Promega; Fitchburg, WI, USA) with the Applied Biosystems High-Capacity cDNA Reverse Transcription Kit (Thermo Fischer; Waltham, MA, USA). Conventional Taq-PCR was performed with primers specific to the respective viral gene (Table S2). Amplification of BAC DNA was used as a positive control. Mock-infected cells and samples without reverse transcriptase to exclude a contamination with genomic DNA were included as negative controls.

2.10. Data Availability

The RNA-seq raw data were deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-7772. A supplementary GFF file for the reference sequence NC_002229.3 containing annotations for all detected introns, poly(A) cleavage sites and associated motifs, as well as the novel CDS for SORF6 can be found in the Supplementary File S1.

3. Results and Discussion

3.1. The MDV Transcriptional Landscape

To assess the transcriptional landscape of MDV in the primary target cells of lytic replication in vivo, we used a previously established in vitro infection system for primary chicken B cells [19]. B cells were infected for 16 h with the very virulent MDV strain RB1B and the vaccine strain CVI988 and analyzed by high-throughput RNA sequencing.
The overall RNA-seq dataset consisted of 82.6 million reads from three biological replicates of CVI988 (48.2 million reads) and two biological replicates of RB1B (34.4 million reads). A third replicate of RB1B did not yield sufficient amount of reads and was therefore excluded from our analysis. An average of 10.5% and 11.2% of the CVI and RB1B datasets respectively could be mapped to the MDV reference sequence (see Table S3).
The position and direction of mapped reads fitted very well to the previously annotated MDV genes (Figure 1A). Highly abundant genes like the immediate early gene SORF1 (ICP4) or the UL49 tegument protein (VP22) correlate well with previously published data [49]. Surprisingly, only minor differences were detected between the transcriptome of CVI988 and RB1B. Comparing the 94 detected MDV genes, two variants of MDV075 encoding the 14-kDa polypeptides (pp14), were significantly higher expressed in RB1B infected primary B cells (Figure 1B). These phosphorylated cytoplasmic proteins arise from splice variants of the same gene and are thought to be involved in transcriptional regulation and increased neurovirulence [50,51,52]. Furthermore, the hypothetical gene MDV082 [53] that is located on the same transcript as the ICP4 gene, was significantly higher expressed in CVI988 infected B cells. However, the transcriptome of RB1B and CVI988 only shows subtle differences in primary chicken B cells, suggesting that the differences in their pathogenesis might be due to sequence changes on the protein level and/or functional differences of virulence factors. Similarly, we only detected marginal quantitative differences between the expression levels of viral proteins in a proteome analysis using LC-MALDI TOF/TOF mass spectrometry (Figure S1 and Table S1). Taken together, our data indicate that there are only minor differences in the mRNA and protein expression levels after B cell infections with the very virulent RB1B or the CVI988 vaccine strain.

3.2. Splicing of Polycistronic MDV Transcripts

In addition to the transcriptional profile, we could readily identify 71 introns that were represented by at least one de novo transcript (Figure 2A). Some of the introns and associated spliced genes have been previously described such as the viral lipase (vLIP) [54], LORF2 (MDV012) [55], UL15 [12], UL44 (glycoprotein C) [56], vIL8 [57] and pp14 [51]. However, analysis of the MDV transcriptome revealed a number of novel splice forms (Table S4). These results are in line with previous RNA-seq analysis for other alphaherpesviruses [58] that also revealed a plethora of novel splice products. The detected splice variants could contribute to viral proteomic diversity and could prevent viral mRNA degradation through the virion host shutoff UL41 endoribonuclease [59]. For HSV-1 it has been shown that UL41 not only targets many cellular but also viral mRNAs. Spliced mRNAs are protected from UL41-mediated degradation by bound exon junction complexes (EJCs) [60].
The identified splice site sequences mostly represent canonical splicing motifs, containing the GT at the donor and AG at the acceptor sites (Figure 2B,C). The intron length varied between 70 and 8651 nt (Figure 2E). Intriguingly, the intron frequencies differed between RB1B and CVI988 in infected primary chicken B cells (Figure 3, Figure 4 and Figure 5). By matching the intron positions with our MS data, we identified a peptide that spans the exon-exon junction of UL15 (Figure S2).
The analysis of poly(A) cleavage signals within the RNA-seq data revealed abundant bicistronic and polycistronic MDV transcripts (Table S5). These transcripts encode for two or more proteins and were characterized as regions of high coverage that were not separated by a poly(A) cleavage site. Here we found that the canonical AATAAA motif is the most frequent and functional polyadenylation signal in MDV, followed by ATTAAA (Figure 2D). Interestingly, we also found evidence for alternative non-canonical polyadenylation signals in MDV mRNA 3′ UTRs (Figure 2 and Table S5). The distances between the detected AWTAAA polyadenylation signal motif and the actual poly(A) cleavage site was 13.8 nt (±4.4 SD) and confirmed that not only the polyadenylation signal sequence, but also its distance from the poly(A) is highly conserved (Figure 2D,F) [61,62].

3.3. The Transcriptional Makeup of the MDV Unique Long Region (UL).

The MDV unique regions mainly harbor genes that are conserved among alphaherpesviruses and are involved in DNA replication and production of progeny virus [26]. The UL spans over roughly 113,000 base pairs and harbors the majority of the MDV-encoded genes [26]. Within the UL, we could detect high transcription rates of nearly all annotated genes. Splicing was identified in multiple genes including RLORF14 (pp24), vLIP [63], and LORF2 (MDV012) [55] and in a transcript antisense to UL5 (MDV017) (Figure 3). Only minor differences were observed in the intron frequencies between RB1B and CVI988. To confirm the splice events and frequencies detected by RNA-seq, we performed RT-PCR analyses on several randomly selected genes (Figure 3C, Figure 4B, Figure 5C, Figures S2 and S3). All analyzed genes showed a comparable splice pattern in both RNA-seq and RT-PCR. In addition, we confirmed a novel splice site of UL15 by MS and RT-PCR (Figure S2). UL15 encodes the tripartite terminase subunit that is involved in DNA packaging into the viral capsid. Splicing of UL15 mRNA has already been observed in herpes simplex virus type 1 (HSV-1) [64] and duck enteritis virus (DEV) [65]. However, the observed UL15 isoforms in MDV are to our knowledge unknown and expand the number of potential proteins encoded by UL15 to at least five.
We could also confirm known splice sites in UL44 (gC) in our analysis (Figure 4) [56]. These splice variants lead to a gC protein that lacks the transmembrane domain and is secreted into the supernatant [56]. Beyond that, we also confirmed novel splice sites like in RLORF14a (pp38) (Figure S3). Several capsid and tegument protein-encoding genes encoded in the UL and US regions, like UL18 (triplex capsid protein 2), UL19 (VP5), UL49 (VP22), UL49.5 (gN), US7 (gI) and US8 (gE) also undergo different splicing events (Table S4), contradicting the long-standing paradigm that splicing is a rare phenomenon throughout the alphaherpesvirus family [13].
Some identified splice variants would result in proteins with an altered membrane topology (TMHMM Server, v. 2.0). This is for example the case for splice variants of pp24 and of pp38, which results in changes of the previously assessed hydrophobic anchor domains of both proteins [66]. While pp38 seems to exist as splice variants with and without a membrane anchor, splicing of pp28 could retain its membrane association while altering its function (Figure S5).

3.4. The Transcriptional Makeup of the MDV Unique Short Region (US).

The MDV US region contains many genes that play important roles in the viral life cycle. Intriguingly, we detected splice variants of several envelope glycoproteins as described above. In addition, we identified a hitherto uncharacterized spliced transcript of a gene located downstream of SORF2A, termed SORF6 (Figure 5). This novel gene possesses an upstream TATA box in the transcriptional regulatory region, an intron and exon with respective donor and acceptor site and a downstream polyadenylation signal with the poly(A) cleavage cluster (as described in Figure 2). Furthermore, the resulting protein is predicted to be 85 amino acids in size, shows several predicted phosphorylation sites (Figure S6) and may act as a DNA-binding protein (DNABIND server [40]). The region containing this novel gene (Figure 5) has previously been associated with the virulence of the virus [67]. However, more work needs to be done to understand the contribution of this novel gene and the region in MDV virulence.

3.5. The Transcriptional Makeup of the MDV Repeat Regions

The repeat regions mostly contain MDV-specific genes encoding for proteins or RNA that play a role in the cell tropism, MDV pathogenesis, latency, and transformation [14]. Here, we observed excessive splicing antisense to ICP4. These transcripts are part of the latency associated transcript (LAT) region, have a complex splice pattern, and their functions remain largely unknown [68,69]. Some of these RNAs function as MDV-encoded micro RNAs and are described elsewhere [70,71].
Only moderate RLORF7 (Meq) and vIL8 splicing was detected in infected primary B cells 16 hpi although more extensive splicing activity has been observed in this region of the MDV genome in vitro and in vivo [72,73,74]. These splice variants are likely higher expressed in latently infected and transformed cells. The splice variants of the neurovirulence factor pp14 encoded by MDV075 were efficiently detected as published previously [15,44].
In the RNA-seq data, we did not detect any reads complementary to vTR. However, this region is annotated as a hypothetical MDV gene termed RLORF1 (an arginine-rich protein/ICP0-like protein) [75]. RLORF1 is discussed as a potential positional orthologue of alphaherpesviral ICP0 proteins; however, it does not contain typical ICP0 features such as a C3HC4 zinc RING finger at the N-terminus or a nuclear localization signal (NLS). To assess if ICP0 protein is expressed and if it plays a role in replication, we generated recombinant MDV mutants harboring an HA-tagged ICP0 (RB1B_ICP0-HA) or an ICP0 knockout (RB1B_ΔMetICP0). The knockout did not affect MDV replication and cell to cell spread in vitro (Method S1 and Figure S4A) and no ICP0 was detected by western blotting (Method S1 and Figure S4B), suggesting that ICP0 is not expressed and therefore does not play a role in the virus life cycle.
Furthermore, we detected several poly(A) cleavage sites in combination with appropriate motifs within the MDV repeat regions, that are to our knowledge undescribed. The presence of these transcriptional signals in combination with sufficient read coverage suggest the existence of hitherto hypothetical protein coding regions, including RLORF11 and MDV082.

3.6. MDV Noncoding RNAs

Although we enriched for poly(A) mRNA, some newly identified introns do not lie in annotated MDV coding sequences and some splice donor and acceptor sites do not give rise to conclusive protein-encoding mRNAs. Such sequences could be easily regarded as ‘nonsense’ transcripts that are rapidly degraded. However, the importance of noncoding RNAs (ncRNAs) in MDV infections is expanding. Several viral ncRNAs were found to be expressed by MDV [75,76,77,78] but the multitude of functions played by viral ncRNAs, and especially by long ncRNAs (lncRNAs) and stable intronic sequence RNAs (sisRNAs), have not been thoroughly investigated yet and the unexpected transcriptomic complexity may have been overlooked in MDV research so far. Of note is, that antisense transcription was also observed in related alphaherpesviruses like HSV-1 or pseudorabies virus (PRV), and in human herpesvirus 6 (HHV-6) RNA-seq data [79,80,81].
Taken together, our MDV RNA-seq data provide novel insights into the transcriptional profile of the RB1B and CVI988 strains. Despite stark differences in their pathogenicity, the two viruses show a similar transcriptomic profile in primary chicken B cells.

4. Conclusions

B cells are a major target for lytic MDV replication in vivo [8,10]; however, it remained impossible to assess the MDV transcriptome in primary B cells, due to the short-lived nature of these cell. The aim of this study was to evaluate the gene expression profiles of the very virulent RB1B strain and the commercial MDV vaccine CVI988 in primary chicken B cells by RNA-seq using our recently established in vitro infection system [19]. We developed a bioinformatics pipeline that can be easily transferred to other herpesviruses or large DNA viruses to identify unknown transcript isoforms and associated motifs.
The RNA-seq revealed the expression of 94 MDV transcripts and the presence of 71 introns that lead to mostly novel splice forms and antisense transcripts. In addition, we could identify a novel gene in the US region of the MDV genome that we will characterize in future studies. While some of the detected splice sites were previously published, we identified several novel splice variants and confirmed some of them by RT-PCR and/or MS. However, more work is certainly required to dissect their relevance in the MDV life cycle.
We found that MDV produces bicistronic and polycistronic transcripts as a mechanism to maximize its coding capacities. Poly(A) cleavage after the upstream AATAAA motif seems to be the most frequent and functional polyadenylation signal in MDV. The identification of possible alternative transcript termination (ATT) needs further experimental evidence (Figure 2 and Table S5). ATT is a strong regulatory factor in eukaryotes [82], but there is only limited data for ATT in herpesvirus transcription.
The comparison of the transcriptome between the very virulent RB1B strain and the CVI988 vaccine revealed differences in only a few transcripts (Figure 1). However, more work needs to be done to unravel significant differences that could possibly point towards a mechanism of attenuation or provide valuable information for the development of diagnostic tools.
In summary, our data demonstrate that the MDV genome is more complex than previously assumed. It provides a source of reference for MDV transcripts expressed in primary chicken B cells and lays the foundation for future research on MDV-encoded gene products and splice variants.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/3/264/s1: Method S1: Generation and in vitro characterization of RB1BΔMetICP0, Figure S1: Protein concentration (mol%) scatterplot comparing levels of detected viral proteins in RB1B and CVI988 infected primary chicken B cells, Figure S2: Detected splice variants of UL15, Figure S3: RT-PCR confirmation of splicing events in pp38, Figure S4: Plaques size assay and western blot analysis for the ICP0 mutant RB1BΔMetICP0, Figure S5: Transmembrane helix prediction for spliced and unspliced pp24 and pp38, Figure S6: Prediction of serine, threonine or tyrosine phosphorylation sites in the hypothetical MDV protein SORF6 encoded on the US segment, Table S1: Marek’s disease virus (MDV) proteins detected by mass spectrometry, Table S2: Primers used in this study, Table S3: Summary of RNA-seq read mapping, Table S4: Introns identified from RNA-seq data, Table S5: Poly(A) cleavage sites, polyadenylation signals and polycistronic transcripts identified in MDV transcriptomes, Supplementary File S1: Annotations for NC_002229.3 in GFF format.

Author Contributions

Conceptualization: L.D.B., F.P., A.M.K., A.K. and B.B.K.; methodology: L.D.B., F.P., V.I.P., A.K. and B.B.K.; formal analysis: L.D.B., F.P., A.K. and B.B.K..; investigation: L.D.B., F.P., V.I.P. and A.M.K.; resources: S.H.; writing—original draft preparation: L.D.B. and F.P.; writing—review and editing: L.D.B., F.P., A.K. and B.B.K.; visualization: L.D.B., F.P., A.K. and B.B.K.; supervision: D.H., A.K. and B.B.K.; project administration: A.K. and B.B.K.; funding acquisition: A.K. and B.B.K.; all authors read and approved the final manuscript.

Funding

This research was funded by the Deutsche Forschungsgemeinschaft (DFG), grant numbers KA 3492/3-1 and KA 2900/3-1 awarded to B.B.K. and A.K. respectively.

Acknowledgments

We thank Ann Reum and Yu You (Institut of Virology, Freie Universität Berlin), Marina Kohn (Institute for Animal Physiology, Ludwig-Maximilians-Universität München) and Patrick Zitzow (Institute of Diagnostic Virology, Friedrich-Loeffler-Institut) for outstanding technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nair, V. Evolution of Marek’s disease—A paradigm for incessant race between the pathogen and the host. Vet. J. 2005, 170, 175–183. [Google Scholar] [CrossRef]
  2. Davison, T.F.; Nair, V. Marek’s Disease: An Evolving Problem; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  3. Davison, F.; Nair, V. Use of Marek’s disease vaccines: Could they be driving the virus to increasing virulence? Expert Rev. Vaccines 2005, 4, 77–88. [Google Scholar] [CrossRef]
  4. Witter, R.L. Increased virulence of Marek’s disease virus field isolates. Avian Dis. 1997, 41, 149–163. [Google Scholar] [CrossRef]
  5. Baigent, S.J.; Nair, V.K.; Le Galludec, H. Real-time PCR for differential quantification of CVI988 vaccine virus and virulent strains of Marek’s disease virus. J. Virol. Methods 2016, 233, 23–36. [Google Scholar] [CrossRef] [Green Version]
  6. Jarosinski, K.W.; Tischer, B.K.; Trapp, S.; Osterrieder, N. Marek’s disease virus: Lytic replication, oncogenesis and control. Expert Rev. Vaccines 2006, 5, 761–772. [Google Scholar] [CrossRef]
  7. Baigent, S.J.; Ross, L.J.; Davison, T.F. Differential susceptibility to Marek’s disease is associated with differences in number, but not phenotype or location, of pp38+ lymphocytes. J. Gen. Virol. 1998, 79 Pt 11, 2795–2802. [Google Scholar] [CrossRef]
  8. Baaten, B.J.; Staines, K.A.; Smith, L.P.; Skinner, H.; Davison, T.F.; Butter, C. Early replication in pulmonary B cells after infection with Marek’s disease herpesvirus by the respiratory route. Viral Immunol. 2009, 22, 431–444. [Google Scholar] [CrossRef]
  9. Shek, W.; Calnek, B.; Schat, K.; Chen, C. Characterization of Marek’s disease virus-infected lymphocytes: Discrimination between cytolytically and latently infected cells. J. Natl. Cancer Inst. 1983, 70, 485–491. [Google Scholar]
  10. Bertzbach, L.D.; Laparidou, M.; Hartle, S.; Etches, R.J.; Kaspers, B.; Schusser, B.; Kaufer, B.B. Unraveling the role of B cells in the pathogenesis of an oncogenic avian herpesvirus. Proc. Natl. Acad. Sci. USA 2018, 115, 11603–11607. [Google Scholar] [CrossRef]
  11. Calnek, B.W.; Schat, K.A.; Ross, L.J.; Shek, W.R.; Chen, C.L. Further characterization of Marek’s disease virus-infected lymphocytes. I. In vivo infection. Int. J. Cancer 1984, 33, 389–398. [Google Scholar] [CrossRef]
  12. Tulman, E.R.; Afonso, C.L.; Lu, Z.; Zsak, L.; Rock, D.L.; Kutish, G.F. The genome of a very virulent Marek’s disease virus. J. Virol. 2000, 74, 7980–7988. [Google Scholar] [CrossRef]
  13. Davison, A.J. Comparative analysis of the genomes. In Human Herpesviruses: Biology, Therapy, and Immunoprophylaxis; Arvin, A., Campadelli-Fiume, G., Mocarski, E., Moore, P.S., Roizman, B., Whitley, R., Yamanishi, K., Eds.; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  14. Bertzbach, L.D.; Kheimar, A.; Ali, F.A.Z.; Kaufer, B.B. Viral Factors Involved in Marek’s Disease Virus (MDV) Pathogenesis. Curr. Clin. Microbiol. Rep. 2018, 5, 238–244. [Google Scholar] [CrossRef]
  15. Heidari, M.; Huebner, M.; Kireev, D.; Silva, R.F. Transcriptional profiling of Marek’s disease virus genes during cytolytic and latent infection. Virus Genes 2008, 36, 383–392. [Google Scholar] [CrossRef]
  16. Neerukonda, S.N.; Tavlarides-Hontz, P.; McCarthy, F.; Pendarvis, K.; Parcells, M.S. Comparison of the Transcriptomes and Proteomes of Serum Exosomes from Marek’s Disease Virus-Vaccinated and Protected and Lymphoma-Bearing Chickens. Genes 2019, 10, 116. [Google Scholar] [CrossRef]
  17. Mwangi, W.N.; Vasoya, D.; Kgosana, L.B.; Watson, M.; Nair, V. Differentially expressed genes during spontaneous lytic switch of Marek’s disease virus in lymphoblastoid cell lines determined by global gene expression profiling. J. Gen. Virol. 2017, 98, 779–790. [Google Scholar] [CrossRef]
  18. Baigent, S.J.; Ross, L.J.; Davison, T.F. A flow cytometric method for identifying Marek’s disease virus pp38 expression in lymphocyte subpopulations. Avian Pathol. 1996, 25, 255–267. [Google Scholar] [CrossRef]
  19. Schermuly, J.; Greco, A.; Hartle, S.; Osterrieder, N.; Kaufer, B.B.; Kaspers, B. In vitro model for lytic replication, latency, and transformation of an oncogenic alphaherpesvirus. Proc. Natl. Acad. Sci. USA 2015, 112, 7279–7284. [Google Scholar] [CrossRef] [Green Version]
  20. Schat, K.; Purchase, H. Cell-Culture Methods. A Laboratory Manual for the Isolation and Identification of Avian Pathogens; American Association of Avian Pathologists: Kennett Square, PA, USA, 1998. [Google Scholar]
  21. Martin, A.; Lillehoj, H.S.; Kaspers, B.; Bacon, L.D. Antigen-specific T cell proliferation following coccidia infection. Poult Sci. 1993, 72, 2084–2094. [Google Scholar] [CrossRef]
  22. Kothlow, S.; Morgenroth, I.; Tregaskes, C.A.; Kaspers, B.; Young, J.R. CD40 ligand supports the long-term maintenance and differentiation of chicken B cells in culture. Dev. Comp. Immunol. 2008, 32, 1015–1026. [Google Scholar] [CrossRef]
  23. Engel, A.T.; Selvaraj, R.K.; Kamil, J.P.; Osterrieder, N.; Kaufer, B.B. Marek’s disease viral interleukin-8 promotes lymphoma formation through targeted recruitment of B cells and CD4+ CD25+ T cells. J. Virol. 2012, 86, 8536–8545. [Google Scholar] [CrossRef]
  24. Pfaff, F.; Hägglund, S.; Zoli, M.; Blaise-Boisseau, S.; Laloy, E.; Koethe, S.; Zühlke, D.; Riedel, K.; Zientara, S.; Bakkali-Kassimi, L.; et al. Proteogenomics Uncovers Critical Elements of Host Response in Bovine Soft Palate Epithelial Cells Following In Vitro Infection with Foot-And-Mouth Disease Virus. Viruses 2019, 11, 53. [Google Scholar] [CrossRef]
  25. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  26. Lee, L.F.; Wu, P.; Sui, D.; Ren, D.; Kamil, J.; Kung, H.J.; Witter, R.L. The complete unique long sequence and the overall genomic organization of the GA strain of Marek’s disease virus. Proc. Natl. Acad. Sci. USA 2000, 97, 6091–6096. [Google Scholar] [CrossRef]
  27. Davison, A.J. Herpesvirus systematics. Vet. Microbiol. 2010, 143, 52–69. [Google Scholar] [CrossRef] [Green Version]
  28. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  29. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef]
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2017. [Google Scholar]
  31. Wagih, O. ggseqlogo: A versatile R package for drawing sequence logos. Bioinformatics 2017, 33, 3645–3647. [Google Scholar] [CrossRef]
  32. Bonfert, T.; Kirner, E.; Csaba, G.; Zimmer, R.; Friedel, C.C. ContextMap 2: Fast and accurate context-based RNA-seq mapping. BMC Bioinform. 2015, 16, 122. [Google Scholar] [CrossRef]
  33. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  34. Pauws, E.; van Kampen, A.H.; van de Graaf, S.A.; de Vijlder, J.J.; Ris-Stalpers, C. Heterogeneity in polyadenylation cleavage sites in mammalian mRNA sequences: Implications for SAGE analysis. Nucl. Acids Res. 2001, 29, 1690–1694. [Google Scholar] [CrossRef]
  35. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucl. Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  36. Grant, C.E.; Bailey, T.L.; Noble, W.S. FIMO: Scanning for occurrences of a given motif. Bioinformatics 2011, 27, 1017–1018. [Google Scholar] [CrossRef]
  37. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417–419. [Google Scholar] [CrossRef]
  38. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  39. Blom, N.; Gammeltoft, S.; Brunak, S. Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J. Mol. Biol. 1999, 294, 1351–1362. [Google Scholar] [CrossRef]
  40. Szilagyi, A.; Skolnick, J. Efficient prediction of nucleic acid binding function from low-resolution protein structures. J. Mol. Biol. 2006, 358, 922–933. [Google Scholar] [CrossRef]
  41. Neuhoff, V.; Arold, N.; Taube, D.; Ehrhardt, W. Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie Brilliant Blue G-250 and R-250. Electrophoresis 1988, 9, 255–262. [Google Scholar] [CrossRef]
  42. Laemmli, U.K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970, 227, 680–685. [Google Scholar] [CrossRef]
  43. Wisniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359–362. [Google Scholar] [CrossRef]
  44. Boersema, P.J.; Aye, T.T.; van Veen, T.A.; Heck, A.J.; Mohammed, S. Triplex protein quantification based on stable isotope labeling by peptide dimethylation applied to cell and tissue lysates. Proteomics 2008, 8, 4624–4632. [Google Scholar] [CrossRef]
  45. Pauker, V.I.; Bertzbach, L.D.; Hohmann, A.; Kheimar, A.; Teifke, J.P.; Mettenleiter, T.C.; Karger, A.; Kaufer, B.B. Imaging Mass Spectrometry and Proteome Analysis of Marek’s Disease Virus-Induced Tumors. mSphere 2019, 4. [Google Scholar] [CrossRef]
  46. Henning, A.K.; Groschup, M.H.; Mettenleiter, T.C.; Karger, A. Analysis of the bovine plasma proteome by matrix-assisted laser desorption/ionisation time-of-flight tandem mass spectrometry. Vet. J. 2014, 199, 175–180. [Google Scholar] [CrossRef]
  47. Aken, B.L.; Achuthan, P.; Akanni, W.; Amode, M.R.; Bernsdorff, F.; Bhai, J.; Billis, K.; Carvalho-Silva, D.; Cummins, C.; Clapham, P.; et al. Ensembl 2017. Nucl. Acids Res. 2017, 45, D635–D642. [Google Scholar] [CrossRef]
  48. Ishihama, Y.; Oda, Y.; Tabata, T.; Sato, T.; Nagasu, T.; Rappsilber, J.; Mann, M. Exponentially modified protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein. Mol. Cell. Proteom. 2005, 4, 1265–1272. [Google Scholar] [CrossRef]
  49. Chakraborty, P.; Kuo, R.; Vervelde, L.; Dutia, B.M.; Kaiser, P.; Smith, J. Macrophages from Susceptible and Resistant Chicken Lines have Different Transcriptomes following Marek’s Disease Virus Infection. Genes 2019, 10, 74. [Google Scholar] [CrossRef]
  50. Tahiri-Alaoui, A.; Smith, L.P.; Kgosana, L.; Petherbridge, L.J.; Nair, V. Identification of a neurovirulence factor from Marek’s disease virus. Avian Dis. 2013, 57, 387–394. [Google Scholar] [CrossRef] [PubMed]
  51. Hong, Y.; Coussens, P.M. Identification of an immediate-early gene in the Marek’s disease virus long internal repeat region which encodes a unique 14-kilodalton polypeptide. J. Virol. 1994, 68, 3593–3603. [Google Scholar] [PubMed]
  52. Tahiri-Alaoui, A.; Smith, L.P.; Baigent, S.; Kgosana, L.; Petherbridge, L.J.; Lambeth, L.S.; James, W.; Nair, V. Identification of an intercistronic internal ribosome entry site in a Marek’s disease virus immediate-early gene. J. Virol. 2009, 83, 5846–5853. [Google Scholar] [CrossRef]
  53. Hearn, C.; Preeyanon, L.; Hunt, H.D.; York, I.A. An MHC class I immune evasion gene of Mareks disease virus. Virology 2015, 475, 88–95. [Google Scholar] [CrossRef]
  54. Becker, Y.; Asher, Y.; Tabor, E.; Davidson, I.; Malkinson, M. Open reading frames in a 4556 nucleotide sequence within MDV-1 BamHI-D DNA fragment: Evidence for splicing of mRNA from a new viral glycoprotein gene. Virus Genes 1994, 8, 55–69. [Google Scholar] [CrossRef]
  55. Schippers, T.; Jarosinski, K.; Osterrieder, N. The ORF012 gene of Marek’s disease virus type 1 produces a spliced transcript and encodes a novel nuclear phosphoprotein essential for virus growth. J. Virol. 2015, 89, 1348–1363. [Google Scholar] [CrossRef] [PubMed]
  56. Jarosinski, K.W.; Osterrieder, N. Marek’s disease virus expresses multiple UL44 (gC) variants through mRNA splicing that are all required for efficient horizontal transmission. J. Virol. 2012, 86, 7896–7906. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, J.L.; Lin, S.F.; Xia, L.; Brunovskis, P.; Li, D.; Davidson, I.; Lee, L.F.; Kung, H.J. MEQ and V-IL8: Cellular genes in disguise? Acta Virol. 1999, 43, 94–101. [Google Scholar] [PubMed]
  58. Tombacz, D.; Csabai, Z.; Szucs, A.; Balazs, Z.; Moldovan, N.; Sharon, D.; Snyder, M.; Boldogkoi, Z. Long-Read Isoform Sequencing Reveals a Hidden Complexity of the Transcriptional Landscape of Herpes Simplex Virus Type 1. Front. Microbiol. 2017, 8, 1079. [Google Scholar] [CrossRef] [PubMed]
  59. Gimeno, I.; Silva, R.F. Deletion of the Marek’s disease virus UL41 gene (vhs) has no measurable effect on latency or pathogenesis. Virus Genes 2008, 36, 499–507. [Google Scholar] [CrossRef] [PubMed]
  60. Sadek, J.; Read, G.S. The Splicing History of an mRNA Affects Its Level of Translation and Sensitivity to Cleavage by the Virion Host Shutoff Endonuclease during Herpes Simplex Virus Infections. J. Virol. 2016, 90, 10844–10856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Mandel, C.R.; Bai, Y.; Tong, L. Protein factors in pre-mRNA 3′-end processing. Cell. Mol. Life Sci. 2008, 65, 1099–1122. [Google Scholar] [CrossRef]
  62. Proudfoot, N.J. Ending the message: Poly(A) signals then and now. Genes Dev. 2011, 25, 1770–1782. [Google Scholar] [CrossRef]
  63. Kamil, J.P.; Tischer, B.K.; Trapp, S.; Nair, V.K.; Osterrieder, N.; Kung, H.J. vLIP, a viral lipase homologue, is a virulence factor of Marek’s disease virus. J. Virol. 2005, 79, 6984–6996. [Google Scholar] [CrossRef]
  64. Hardy, W.R.; Sandri-Goldin, R.M. Herpes-Simplex Virus Inhibits Host-Cell Splicing, and Regulatory Protein Icp27 Is Required for This Effect. J. Virol. 1994, 68, 7790–7799. [Google Scholar] [PubMed]
  65. Zhu, H.; Li, H.; Han, Z.; Shao, Y.; Wang, Y.; Kong, X. Identification of a spliced gene from duck enteritis virus encoding a protein homologous to UL15 of herpes simplex virus 1. Virol. J. 2011, 8, 156. [Google Scholar] [CrossRef] [PubMed]
  66. Zhu, G.S.; Iwata, A.; Gong, M.; Ueda, S.; Hirai, K. Marek’s disease virus type 1-specific phosphorylated proteins pp38 and pp24 with common amino acid termini are encoded from the opposite junction regions between the long unique and inverted repeat sequences of viral genome. Virology 1994, 200, 816–820. [Google Scholar] [CrossRef]
  67. Lv, H.; Zhang, Y.; Sun, G.; Bao, K.; Gao, Y.; Qi, X.; Cui, H.; Wang, Y.; Li, K.; Gao, L.; et al. Genetic evolution of Gallid herpesvirus 2 isolated in China. Infect. Genet. Evol. 2017, 51, 263–274. [Google Scholar] [CrossRef] [PubMed]
  68. Morgan, R.W.; Xie, Q.; Cantello, J.L.; Miles, A.M.; Bernberg, E.L.; Kent, J.; Anderson, A. Marek’s disease virus latency. Curr. Top. Microbiol. Immunol. 2001, 255, 223–243. [Google Scholar]
  69. Cantello, J.L.; Parcells, M.S.; Anderson, A.S.; Morgan, R.W. Marek’s disease virus latency-associated transcripts belong to a family of spliced RNAs that are antisense to the ICP4 homolog gene. J. Virol. 1997, 71, 1353–1361. [Google Scholar] [PubMed]
  70. Burnside, J.; Bernberg, E.; Anderson, A.; Lu, C.; Meyers, B.C.; Green, P.J.; Jain, N.; Isaacs, G.; Morgan, R.W. Marek’s disease virus encodes MicroRNAs that map to meq and the latency-associated transcript. J. Virol. 2006, 80, 8778–8786. [Google Scholar] [CrossRef] [PubMed]
  71. Yao, Y.; Zhao, Y.; Xu, H.; Smith, L.P.; Lawrie, C.H.; Watson, M.; Nair, V. MicroRNA profile of Marek’s disease virus-transformed T-cell line MSB-1: Predominance of virus-encoded microRNAs. J. Virol. 2008, 82, 4007–4015. [Google Scholar] [CrossRef]
  72. Jarosinski, K.W.; Schat, K.A. Multiple alternative splicing to exons II and III of viral interleukin-8 (vIL-8) in the Marek’s disease virus genome: The importance of vIL-8 exon I. Virus Genes 2007, 34, 9–22. [Google Scholar] [CrossRef]
  73. Anobile, J.M.; Arumugaswami, V.; Downs, D.; Czymmek, K.; Parcells, M.; Schmidt, C.J. Nuclear localization and dynamic properties of the Marek’s disease virus oncogene products Meq and Meq/vIL8. J. Virol. 2006, 80, 1160–1166. [Google Scholar] [CrossRef]
  74. Okada, T.; Takagi, M.; Murata, S.; Onuma, M.; Ohashi, K. Identification and characterization of a novel spliced form of the meq transcript in lymphoblastoid cell lines derived from Marek’s disease tumours. J. Gen. Virol. 2007, 88, 2111–2120. [Google Scholar] [CrossRef]
  75. Trapp, S.; Parcells, M.S.; Kamil, J.P.; Schumacher, D.; Tischer, B.K.; Kumar, P.M.; Nair, V.K.; Osterrieder, N. A virus-encoded telomerase RNA promotes malignant T cell lymphomagenesis. J. Exp. Med. 2006, 203, 1307–1317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Tycowski, K.T.; Guo, Y.E.; Lee, N.; Moss, W.N.; Vallery, T.K.; Xie, M.; Steitz, J.A. Viral noncoding RNAs: More surprises. Genes Dev. 2015, 29, 567–584. [Google Scholar] [CrossRef] [PubMed]
  77. Zhuang, G.; Sun, A.; Teng, M.; Luo, J. A Tiny RNA that Packs a Big Punch: The Critical Role of a Viral miR-155 Ortholog in Lymphomagenesis in Marek’s Disease. Front. Microbiol. 2017, 8, 1169. [Google Scholar] [CrossRef] [PubMed]
  78. Figueroa, T.; Boumart, I.; Coupeau, D.; Rasschaert, D. Hyperediting by ADAR1 of a new herpesvirus lncRNA during the lytic phase of the oncogenic Marek’s disease virus. J. Gen. Virol. 2016, 97, 2973–2988. [Google Scholar] [CrossRef]
  79. Wyler, E.; Menegatti, J.; Franke, V.; Kocks, C.; Boltengagen, A.; Hennig, T.; Theil, K.; Rutkowski, A.; Ferrai, C.; Baer, L.; et al. Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection. Genome Biol. 2017, 18, 209. [Google Scholar] [CrossRef] [PubMed]
  80. Greninger, A.L.; Knudsen, G.M.; Roychoudhury, P.; Hanson, D.J.; Sedlak, R.H.; Xie, H.; Guan, J.; Nguyen, T.; Peddu, V.; Boeckh, M.; et al. Comparative genomic, transcriptomic, and proteomic reannotation of human herpesvirus 6. BMC Genom. 2018, 19, 204. [Google Scholar] [CrossRef]
  81. Olah, P.; Tombacz, D.; Poka, N.; Csabai, Z.; Prazsak, I.; Boldogkoi, Z. Characterization of pseudorabies virus transcriptome by Illumina sequencing. BMC Microbiol. 2015, 15, 130. [Google Scholar] [CrossRef] [PubMed]
  82. Shabalina, S.A.; Ogurtsov, A.Y.; Spiridonov, N.A.; Koonin, E.V. Evolution at protein ends: Major contribution of alternative transcription initiation and termination to the transcriptome and proteome diversity in mammals. Nucl. Acids Res. 2014, 42, 7132–7144. [Google Scholar] [CrossRef] [PubMed]
Viruses 11 00264 g001 550
Figure 1. The Marek’s disease virus (MDV) transcriptome and proteome in in vitro infected primary chicken B cells. (A) Visualization of the deduplicated and strand-specific RNA-seq read coverage of plus (+) and minus (-) strand-encoded genes across the MDV genome. Orange curves indicate CVI988 reads and grey curves indicate reads for RB1B (with the respective annotated genes as grey arrow bars). Red bars depict proteins identified by MS. The two unique regions, unique long (UL) and short (US) are flanked by terminal (TRL and TRS) and internal (IRL and IRS) inverted repeat regions. Nucleotide position numbers are derived from [12]. (B) Gene expression scatterplot comparing normalized expression levels (rlog) in RB1B and CVI988 infected primary chicken B cells. Red dots indicate significantly differentially expressed genes.
Figure 1. The Marek’s disease virus (MDV) transcriptome and proteome in in vitro infected primary chicken B cells. (A) Visualization of the deduplicated and strand-specific RNA-seq read coverage of plus (+) and minus (-) strand-encoded genes across the MDV genome. Orange curves indicate CVI988 reads and grey curves indicate reads for RB1B (with the respective annotated genes as grey arrow bars). Red bars depict proteins identified by MS. The two unique regions, unique long (UL) and short (US) are flanked by terminal (TRL and TRS) and internal (IRL and IRS) inverted repeat regions. Nucleotide position numbers are derived from [12]. (B) Gene expression scatterplot comparing normalized expression levels (rlog) in RB1B and CVI988 infected primary chicken B cells. Red dots indicate significantly differentially expressed genes.
Viruses 11 00264 g001
Viruses 11 00264 g002 550
Figure 2. Overview of Marek’s disease virus (MDV) mRNA splicing and poly(A) cleavage. (A) Visualization of the cumulative CVI988 and RB1B RNA-seq coverage of plus (+) and minus (-) strand-encoded genes across the MDV genome. Black arrow bars indicate introns and dashed lines indicate poly(A) cleavage sites. More detailed information is shown in Table S4 and Table S5. (B) Nucleotide frequency maps (sequence logo) of splice donor sites in MDV-encoded mRNAs. The relative heights of letters correspond to frequencies of bases at each position. (C) Sequence logo of splice acceptor sites in MDV-encoded mRNAs. (D) Sequence logo of polyadenylation signals in MDV-encoded three prime untranslated regions (3′-UTR). (E) Histogram depicting MDV intron length distributions. (F) Histogram depicting the distance from AWTAAA-like motifs to the poly(A) cleavage site.
Figure 2. Overview of Marek’s disease virus (MDV) mRNA splicing and poly(A) cleavage. (A) Visualization of the cumulative CVI988 and RB1B RNA-seq coverage of plus (+) and minus (-) strand-encoded genes across the MDV genome. Black arrow bars indicate introns and dashed lines indicate poly(A) cleavage sites. More detailed information is shown in Table S4 and Table S5. (B) Nucleotide frequency maps (sequence logo) of splice donor sites in MDV-encoded mRNAs. The relative heights of letters correspond to frequencies of bases at each position. (C) Sequence logo of splice acceptor sites in MDV-encoded mRNAs. (D) Sequence logo of polyadenylation signals in MDV-encoded three prime untranslated regions (3′-UTR). (E) Histogram depicting MDV intron length distributions. (F) Histogram depicting the distance from AWTAAA-like motifs to the poly(A) cleavage site.
Viruses 11 00264 g002
Viruses 11 00264 g003 550
Figure 3. The Marek’s disease virus (MDV) unique long (UL) region. (A) Visualization of RNA-Seq coverage across parts of the MDV UL region with respective introns in black. Green and red arrows indicate the polyadenylation signal and the poly(A) cleavage site respectively. Underlying black arrows suggest unspliced (mono-, bi-, or polycistronic) mRNAs. (B) Comparison of intron frequencies in RB1B and CVI988 infected primary chicken B cells. I1: pp24, I2: vLIP, I3: LORF2 (MDV012), I4: transcript antisense to UL5 (MDV017). (C) RT-PCR was performed to validate the splicing event. PCR products were derived using forward/reverse primers to amplify the respective intron-flanking regions. The representative gel images illustrate the results of RT-PCR analysis. The black arrows indicate the spliced form of the respective gene.
Figure 3. The Marek’s disease virus (MDV) unique long (UL) region. (A) Visualization of RNA-Seq coverage across parts of the MDV UL region with respective introns in black. Green and red arrows indicate the polyadenylation signal and the poly(A) cleavage site respectively. Underlying black arrows suggest unspliced (mono-, bi-, or polycistronic) mRNAs. (B) Comparison of intron frequencies in RB1B and CVI988 infected primary chicken B cells. I1: pp24, I2: vLIP, I3: LORF2 (MDV012), I4: transcript antisense to UL5 (MDV017). (C) RT-PCR was performed to validate the splicing event. PCR products were derived using forward/reverse primers to amplify the respective intron-flanking regions. The representative gel images illustrate the results of RT-PCR analysis. The black arrows indicate the spliced form of the respective gene.
Viruses 11 00264 g003
Viruses 11 00264 g004 550
Figure 4. Glycoprotein C (gC) splicing. (A) Visualization of RNA-seq coverage across the Marek’s disease virus (MDV) gC gene UL44 with respective introns in black. The green and red arrow indicate the canonical ATTAAA polyadenylation signal and the poly(A) cleavage site respectively. (B) RT-PCR was performed to validate gC splicing. PCR products were derived using forward/reverse primers to amplify the respective intron-flanking regions. (C) Comparison of the gC intron frequencies in RB1B and CVI988 infected primary chicken B cells.
Figure 4. Glycoprotein C (gC) splicing. (A) Visualization of RNA-seq coverage across the Marek’s disease virus (MDV) gC gene UL44 with respective introns in black. The green and red arrow indicate the canonical ATTAAA polyadenylation signal and the poly(A) cleavage site respectively. (B) RT-PCR was performed to validate gC splicing. PCR products were derived using forward/reverse primers to amplify the respective intron-flanking regions. (C) Comparison of the gC intron frequencies in RB1B and CVI988 infected primary chicken B cells.
Viruses 11 00264 g004
Viruses 11 00264 g005 550
Figure 5. The Marek’s disease virus (MDV) unique short (US) region. (A) Visualization of RNA-seq coverage across the MDV US region with respective introns in black. (B) Zoom into the far-right region of the MDV US with depiction of the novel US gene SORF6. Green and red arrows indicate the polyadenylation signals and the poly(A) cleavage sites, respectively. (C) RT-PCR was performed to validate the splicing event in the novel gene SORF6. PCR products were derived using forward/reverse primers to amplify the respective intron-flanking regions as full-length (upper band) and spliced (lower band). (D) Comparison of the novel US gene (SORF6) intron frequencies in RB1B and CVI988 infected primary chicken B cells.
Figure 5. The Marek’s disease virus (MDV) unique short (US) region. (A) Visualization of RNA-seq coverage across the MDV US region with respective introns in black. (B) Zoom into the far-right region of the MDV US with depiction of the novel US gene SORF6. Green and red arrows indicate the polyadenylation signals and the poly(A) cleavage sites, respectively. (C) RT-PCR was performed to validate the splicing event in the novel gene SORF6. PCR products were derived using forward/reverse primers to amplify the respective intron-flanking regions as full-length (upper band) and spliced (lower band). (D) Comparison of the novel US gene (SORF6) intron frequencies in RB1B and CVI988 infected primary chicken B cells.
Viruses 11 00264 g005

Share and Cite

MDPI and ACS Style

Bertzbach, L.D.; Pfaff, F.; Pauker, V.I.; Kheimar, A.M.; Höper, D.; Härtle, S.; Karger, A.; Kaufer, B.B. The Transcriptional Landscape of Marek’s Disease Virus in Primary Chicken B Cells Reveals Novel Splice Variants and Genes. Viruses 2019, 11, 264. https://doi.org/10.3390/v11030264

AMA Style

Bertzbach LD, Pfaff F, Pauker VI, Kheimar AM, Höper D, Härtle S, Karger A, Kaufer BB. The Transcriptional Landscape of Marek’s Disease Virus in Primary Chicken B Cells Reveals Novel Splice Variants and Genes. Viruses. 2019; 11(3):264. https://doi.org/10.3390/v11030264

Chicago/Turabian Style

Bertzbach, Luca D., Florian Pfaff, Viktoria I. Pauker, Ahmed M. Kheimar, Dirk Höper, Sonja Härtle, Axel Karger, and Benedikt B. Kaufer. 2019. "The Transcriptional Landscape of Marek’s Disease Virus in Primary Chicken B Cells Reveals Novel Splice Variants and Genes" Viruses 11, no. 3: 264. https://doi.org/10.3390/v11030264

APA Style

Bertzbach, L. D., Pfaff, F., Pauker, V. I., Kheimar, A. M., Höper, D., Härtle, S., Karger, A., & Kaufer, B. B. (2019). The Transcriptional Landscape of Marek’s Disease Virus in Primary Chicken B Cells Reveals Novel Splice Variants and Genes. Viruses, 11(3), 264. https://doi.org/10.3390/v11030264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop