Characterization of pseudorabies virus transcriptome by Illumina sequencing

Identification of herpesvirus transcripts from genomic regions around the replication origins

Long-read sequencing (LRS) techniques enable the identification of full-length RNA molecules in a single run eliminating the need for additional assembly steps. LRS research has exposed unanticipated transcriptomic complexity in various organisms, including viruses. Herpesviruses are known to produce a range of transcripts, either close to or overlapping replication origins (Oris) and neighboring genes related to transcription or replication, which possess confirmed or potential regulatory roles. In our research, we employed both new and previously published LRS and short-read sequencing datasets to uncover additional Ori-proximal transcripts in nine herpesviruses from all three subfamilies (alpha, beta and gamma). We discovered novel long non-coding RNAs, as well as splice and length isoforms of mRNAs. Moreover, our analysis uncovered an intricate network of transcriptional overlaps within the examined genomic regions. We demonstrated that herpesviruses display distinct patterns of transcriptional overlaps in the vicinity of or at the Oris. Our findings suggest the existence of a 'super regulatory center' in the genome of alphaherpesviruses that governs the initiation of both DNA replication and global transcription through multilayered interactions among the molecular machineries.

Hybrid sequencing discloses unique aspects of the transcriptomic architecture in equid alphaherpesvirus 1

This study employed both short-read sequencing (SRS, Illumina) and long-read sequencing (LRS Oxford Nanopore Technologies) platforms to conduct a comprehensive analysis of the equid alphaherpesvirus 1 (EHV-1) transcriptome. The study involved the annotation of canonical mRNAs and their transcript variants, encompassing transcription start site (TSS) and transcription end site (TES) isoforms, in addition to alternative splicing forms. Furthermore, the study revealed the presence of numerous non-coding RNA (ncRNA) molecules, including intergenic and antisense transcripts, produced by EHV-1. An intriguing finding was the abundant production of chimeric transcripts, some of which potentially encode fusion polypeptides. Moreover, EHV-1 exhibited a greater incidence of transcriptional overlaps and splicing compared to related viruses. It is noteworthy that many genes have their unique TESs along with the co-terminal transcription ends, a characteristic scarcely seen in other alphaherpesviruses. The study also identified transcripts that overlap the replication origins of the virus. Moreover, a novel ncRNA, referred to as NOIR, was found to intersect with the 5'-ends of longer transcript isoform specified by the major transactivator genes ORF64 and ORF65, surrounding the OriL. These findings together imply the existence of a key regulatory mechanism that governs both transcription and replication through, among others, a process that involves interference between the DNA and RNA synthesis machineries.

Role of Epitranscriptomic and Epigenetic Modifications during the Lytic and Latent Phases of Herpesvirus Infections

Herpesviruses are double-stranded DNA viruses occurring at a high prevalence in the human population and are responsible for a wide array of clinical manifestations and diseases, from mild to severe. These viruses are classified in three subfamilies (Alpha-, Beta- and Gammaherpesvirinae), with eight members currently known to infect humans. Importantly, all herpesviruses can establish lifelong latent infections with symptomatic or asymptomatic lytic reactivations. Accumulating evidence suggest that chemical modifications of viral RNA and DNA during the lytic and latent phases of the infections caused by these viruses, are likely to play relevant roles in key aspects of the life cycle of these viruses by modulating and regulating their replication, establishment of latency and evasion of the host antiviral response. Here, we review and discuss current evidence regarding epitranscriptomic and epigenetic modifications of herpesviruses and how these can influence their life cycles. While epitranscriptomic modifications such as m⁶A are the most studied to date and relate to positive effects over the replication of herpesviruses, epigenetic modifications of the viral genome are generally associated with defense mechanisms of the host cells to suppress viral gene transcription. However, herpesviruses can modulate these modifications to their own benefit to persist in the host, undergo latency and sporadically reactivate.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

The Activity of Plant-Derived Ren’s Oligopeptides-1 against the Pseudorabies Virus

Newly synthesized Ren’s oligopeptides-1 was found to have an antiviral effect in clinical trials, and the purpose of this study was to further demonstrate the antiviral activity of Ren’s oligopeptides-1 against the PRV 152-GFP strain. We used the real-time cell analysis system (RTCA) to detect the cytotoxicity of different concentrations of Ren’s oligopeptides-1. We then applied high content screening (HCS) to detect the antiviral activity of Ren’s oligopeptides-1 against PRV. Meanwhile, the fluorescence signal of the virus was collected in real time and the expression levels of the related genes in the PK15 cells infected with PRV were detected using real-time PCR. At the mRNA level, we discovered that, at a concentration of 6 mg/mL, Ren’s oligopeptides-1 reduced the expression of pseudorabies virus (PRV) genes such as IE180, UL18, UL54, and UL21 at a concentration of 6 mg/mL. We then determined that Ren’s oligopeptides-1 has an EC50 value of 6 mg/mL, and at this level, no cytotoxicity was observed.

Dual isoform sequencing reveals complex transcriptomic and epitranscriptomic landscapes of a prototype baculovirus

In this study, two long-read sequencing (LRS) techniques, MinION from Oxford Nanopore Technologies and Sequel from the Pacific Biosciences, were used for the transcriptional characterization of a prototype baculovirus, Autographa californica multiple nucleopolyhedrovirus. LRS is able to read full-length RNA molecules, and thereby distinguish between transcript isoforms, mono- and polycistronic RNAs, and overlapping transcripts. Altogether, we detected 875 transcript species, of which 759 were novel and 116 were annotated previously. These RNA molecules include 41 novel putative protein coding transcripts [each containing 5'-truncated in-frame open reading frames (ORFs), 14 monocistronic transcripts, 99 polygenic RNAs, 101 non-coding RNAs, and 504 untranslated region isoforms. This work also identified novel replication origin-associated transcripts, upstream ORFs, cis-regulatory sequences and poly(A) sites. We also detected RNA methylation in 99 viral genes and RNA hyper-editing in the longer 5'-UTR transcript isoform of the canonical ORF 19 transcript.

Integrative profiling of Epstein-Barr virus transcriptome using a multiplatform approach

BACKGROUND

Epstein-Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV.

METHODS

In this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other's data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory.

RESULTS

This study detected novel genes embedded into longer host genes containing 5'-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed.

CONCLUSIONS

An integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.

The role of Meq-vIL8 in regulating Marek's disease virus pathogenesis

Marek's disease virus (MDV) is a highly cell-associated oncogenic alphaherpesvirus that causes T cell lymphoma in chickens. MDV-encoded Meq and vIL8 proteins play important roles in transformation and early cytolytic infection, respectively. Previous studies identified a spliced transcript, meq-vIL8, formed by alternative splicing of meq and vIL8 genes in MDV lymphoblastoid tumour cells. To determine the role of Meq-vIL8 in MDV pathogenesis, we generated a recombinant MDV (MDV-meqΔSD) by mutating the splice donor site in the meq gene to abrogate the expression of Meq-vIL8. As expected, our results show that MDV-meqΔSD virus grows similarly to the parental and revertant viruses in cell culture, suggesting that Meq-vIL8 is dispensable for MDV growth in vitro. We further characterized the pathogenic properties of MDV-meqΔSD virus in chickens. Our results show that lack of Meq-vIL8 did not affect virus replication during the early cytolytic phase, as determined by immunohistochemistry analysis and/or viral genome copy number, but significantly enhanced viral DNA load in the late phase of infection in the spleen and brain of infected chickens. In addition, we observed that abrogation of Meq-vIL8 expression reduced the mean death time and increased the prevalence of persistent neurological disease, common features of highly virulent strains of MDV, in inoculated chickens. In conclusion, our study shows that Meq-vIL8 is an important virulence factor of MDV.

An Integrated Sequencing Approach for Updating the Pseudorabies Virus Transcriptome

In the last couple of years, the implementation of long-read sequencing (LRS) technologies for transcriptome profiling has uncovered an extreme complexity of viral gene expression. In this study, we carried out a systematic analysis on the pseudorabies virus transcriptome by combining our current data obtained by using Pacific Biosciences Sequel and Oxford Nanopore Technologies MinION sequencing with our earlier data generated by other LRS and short-read sequencing techniques. As a result, we identified a number of novel genes, transcripts, and transcript isoforms, including splice and length variants, and also confirmed earlier annotated RNA molecules. One of the major findings of this study is the discovery of a large number of 5′-truncations of larger putative mRNAs being 3′-co-terminal with canonical mRNAs of PRV. A large fraction of these putative RNAs contain in-frame ATGs, which might initiate translation of N-terminally truncated polypeptides. Our analyses indicate that CTO-S, a replication origin-associated RNA molecule is expressed at an extremely high level. This study demonstrates that the PRV transcriptome is much more complex than previously appreciated.

Latent pseudorabies virus infection in medulla oblongata from quarantined pigs

Pseudorabies virus (PRV) is a major pathogen in pig husbandry and is also a risk to human well-being. Pigs with latent PRV infection carry the virus lifelong, and it can be activated under conducive conditions. This poses a very important challenge to the control of the virus and may even prevent its elimination. To investigate latent infection with wild-type (wt) PRV, and also infection due to the use of live attenuated vaccines on farms, 80 pigs from two large-scale swine operations were traced. At 6 months old, the quarantined pigs were slaughtered and brain samples were collected. A PCR assay targeting the gB and gE genes was developed to detect PRV DNA fragments in medulla oblongata. Five of the samples (6.3%) were gB and gE gene fragment double-positive, 60 of the samples (75%) were gB single-positive, and 15 samples (18.7%) showed double-negative. A portion of latency-associated transcripts (LATs), EP0 mRNA, were found to be present in the gB gene fragment positive samples. Furthermore, the five double-positive samples were transmitted blindly, and apparent cytopathic effects were found in three of the five samples in the fourth generation. By means of Western blotting, PCR and sequencing, two of the isolated viruses were found to be related to vaccine strain Bartha-K61. Another was closely related to domestic epidemic strains HN1201 and LA and relatively unrelated to other Asian isolates. These results suggest that the live vaccines are latently present in brains, in a manner similar to wt PRV, and this poses potential safety issues in the pig husbandry industry. Wt PRV and live vaccine viruses were found to co-exist in pigs, demonstrating that the live vaccines were unable to confer complete sterilizing immunity, which may explain outbreaks of pseudorabies on vaccinated farms.

Multiple Long-Read Sequencing Survey of Herpes Simplex Virus Dynamic Transcriptome

Long-read sequencing (LRS) has become increasingly important in RNA research due to its strength in resolving complex transcriptomic architectures. In this regard, currently two LRS platforms have demonstrated adequate performance: the Single Molecule Real-Time Sequencing by Pacific Biosciences (PacBio) and the nanopore sequencing by Oxford Nanopore Technologies (ONT). Even though these techniques produce lower coverage and are more error prone than short-read sequencing, they continue to be more successful in identifying polycistronic RNAs, transcript isoforms including splice and transcript end variants, as well as transcript overlaps. Recent reports have successfully applied LRS for the investigation of the transcriptome of viruses belonging to various families. These studies have substantially increased the number of previously known viral RNA molecules. In this work, we used the Sequel and MinION technique from PacBio and ONT, respectively, to characterize the lytic transcriptome of the herpes simplex virus type 1 (HSV-1). In most samples, we analyzed the poly(A) fraction of the transcriptome, but we also performed random oligonucleotide-based sequencing. Besides cDNA sequencing, we also carried out native RNA sequencing. Our investigations identified more than 2,300 previously undetected transcripts, including coding, and non-coding RNAs, multi-splice transcripts, as well as polycistronic and complex transcripts. Furthermore, we found previously unsubstantiated transcriptional start sites, polyadenylation sites, and splice sites. A large number of novel transcriptional overlaps were also detected. Random-primed sequencing revealed that each convergent gene pair produces non-polyadenylated read-through RNAs overlapping the partner genes. Furthermore, we identified novel replication-associated transcripts overlapping the HSV-1 replication origins, and novel LAT variants with very long 5' regions, which are co-terminal with the LAT-0.7kb transcript. Overall, our results demonstrated that the HSV-1 transcripts form an extremely complex pattern of overlaps, and that entire viral genome is transcriptionally active. In most viral genes, if not in all, both DNA strands are expressed.

Multiplatform next-generation sequencing identifies novel RNA molecules and transcript isoforms of the endogenous retrovirus isolated from cultured cells

In this study, we applied short- and long-read RNA sequencing techniques, as well as PCR analysis to investigate the transcriptome of the porcine endogenous retrovirus (PERV) expressed from cultured porcine kidney cell line PK-15. This analysis has revealed six novel transcripts and eight transcript isoforms, including five length and three splice variants. We were able to establish whether a deletion in a transcript is the result of the splicing of mRNAs or of genomic deletion in one of the PERV clones. Additionally, we re-annotated the formerly identified RNA molecules. Our analysis revealed a higher complexity of PERV transcriptome than it was earlier believed.

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

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.

Long-read sequencing uncovers a complex transcriptome topology in varicella zoster virus

BACKGROUND

Varicella zoster virus (VZV) is a human pathogenic alphaherpesvirus harboring a relatively large DNA molecule. The VZV transcriptome has already been analyzed by microarray and short-read sequencing analyses. However, both approaches have substantial limitations when used for structural characterization of transcript isoforms, even if supplemented with primer extension or other techniques. Among others, they are inefficient in distinguishing between embedded RNA molecules, transcript isoforms, including splice and length variants, as well as between alternative polycistronic transcripts. It has been demonstrated in several studies that long-read sequencing is able to circumvent these problems.

RESULTS

In this work, we report the analysis of the VZV lytic transcriptome using the Oxford Nanopore Technologies sequencing platform. These investigations have led to the identification of 114 novel transcripts, including mRNAs, non-coding RNAs, polycistronic RNAs and complex transcripts, as well as 10 novel spliced transcripts and 25 novel transcription start site isoforms and transcription end site isoforms. A novel class of transcripts, the nroRNAs are described in this study. These transcripts are encoded by the genomic region located in close vicinity to the viral replication origin. We also show that the ORF63 exhibits a complex structural variation encompassing the splice sites of VZV latency transcripts. Additionally, we have detected RNA editing in a novel non-coding RNA molecule.

CONCLUSIONS

Our investigations disclosed a composite transcriptomic architecture of VZV, including the discovery of novel RNA molecules and transcript isoforms, as well as a complex meshwork of transcriptional read-throughs and overlaps. The results represent a substantial advance in the annotation of the VZV transcriptome and in understanding the molecular biology of the herpesviruses in general.

Long-Read Sequencing Revealed an Extensive Transcript Complexity in Herpesviruses

Long-read sequencing (LRS) techniques are very recent advancements, but they have already been used for transcriptome research in all of the three subfamilies of herpesviruses. These techniques have multiplied the number of known transcripts in each of the examined viruses. Meanwhile, they have revealed a so far hidden complexity of the herpesvirus transcriptome with the discovery of a large number of novel RNA molecules, including coding and non-coding RNAs, as well as transcript isoforms, and polycistronic RNAs. Additionally, LRS techniques have uncovered an intricate meshwork of transcriptional overlaps between adjacent and distally located genes. Here, we review the contribution of LRS to herpesvirus transcriptomics and present the complexity revealed by this technology, while also discussing the functional significance of this phenomenon.

Transcriptome-wide survey of pseudorabies virus using next- and third-generation sequencing platforms

Pseudorabies virus (PRV) is an alphaherpesvirus of swine. PRV has a large double-stranded DNA genome and, as the latest investigations have revealed, a very complex transcriptome. Here, we present a large RNA-Seq dataset, derived from both short- and long-read sequencing. The dataset contains 1.3 million 100 bp paired-end reads that were obtained from the Illumina random-primed libraries, as well as 10 million 50 bp single-end reads generated by the Illumina polyA-seq. The Pacific Biosciences RSII non-amplified method yielded 57,021 reads of inserts (ROIs) aligned to the viral genome, the amplified method resulted in 158,396 PRV-specific ROIs, while we obtained 12,555 ROIs using the Sequel platform. The Oxford Nanopore's MinION device generated 44,006 reads using their regular cDNA-sequencing method, whereas 29,832 and 120,394 reads were produced by using the direct RNA-sequencing and the Cap-selection protocols, respectively. The raw reads were aligned to the PRV reference genome (KJ717942.1). Our provided dataset can be used to compare different sequencing approaches, library preparation methods, as well as for validation and testing bioinformatic pipelines.

Third-generation Sequencing Reveals Extensive Polycistronism and Transcriptional Overlapping in a Baculovirus

The Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is an insect-pathogen baculovirus. In this study, we applied the Oxford Nanopore Technologies platform for the analysis of the polyadenylated fraction of the viral transcriptome using both cDNA and direct RNA sequencing methods. We identified and annotated altogether 132 novel transcripts and transcript isoforms, including 4 coding and 4 non-coding RNA molecules, 47 length variants, 5 splice isoforms, as well as 23 polycistronic and 49 complex transcripts. All of the identified novel protein-coding genes were 5'-truncated forms of longer host genes. In this work, we demonstrated that in the case of transcript start site isoforms, the promoters and the initiator sequence of the longer and shorter variants belong to the same kinetic class. Long-read sequencing also revealed a complex meshwork of transcriptional overlaps, the function of which needs to be clarified. Additionally, we developed bioinformatics methods to improve the transcript annotation and to eliminate the non-specific transcription reads generated by template switching and false priming.

Multi-Platform Sequencing Approach Reveals a Novel Transcriptome Profile in Pseudorabies Virus

Third-generation sequencing is an emerging technology that is capable of solving several problems that earlier approaches were not able to, including the identification of transcripts isoforms and overlapping transcripts. In this study, we used long-read sequencing for the analysis of pseudorabies virus (PRV) transcriptome, including Oxford Nanopore Technologies MinION, PacBio RS-II, and Illumina HiScanSQ platforms. We also used data from our previous short-read and long-read sequencing studies for the comparison of the results and in order to confirm the obtained data. Our investigations identified 19 formerly unknown putative protein-coding genes, all of which are 5' truncated forms of earlier annotated longer PRV genes. Additionally, we detected 19 non-coding RNAs, including 5' and 3' truncated transcripts without in-frame ORFs, antisense RNAs, as well as RNA molecules encoded by those parts of the viral genome where no transcription had been detected before. This study has also led to the identification of three complex transcripts and 50 distinct length isoforms, including transcription start and end variants. We also detected 121 novel transcript overlaps, and two transcripts that overlap the replication origins of PRV. Furthermore, in silico analysis revealed 145 upstream ORFs, many of which are located on the longer 5' isoforms of the transcripts.

Evaluation of the impact of ul54 gene-deletion on the global transcription and DNA replication of pseudorabies virus

Pseudorabies virus (PRV) is an animal alphaherpesvirus with a wide host range. PRV has 67 protein-coding genes and several non-coding RNA molecules, which can be classified into three temporal groups, immediate early, early and late classes. The ul54 gene of PRV and its homolog icp27 of herpes simplex virus have a multitude of functions, including the regulation of viral DNA synthesis and the control of the gene expression. Therefore, abrogation of PRV ul54 function was expected to exert a significant effect on the global transcriptome and on DNA replication. Real-time PCR and real-time RT-PCR platforms were used to investigate these presumed effects. Our analyses revealed a drastic impact of the ul54 mutation on the genome-wide expression of PRV genes, especially on the transcription of the true late genes. A more than two hour delay was observed in the onset of DNA replication, and the amount of synthesized DNA molecules was significantly decreased in comparison to the wild-type virus. Furthermore, in this work, we were able to successfully demonstrate the utility of long-read SMRT sequencing for genotyping of mutant viruses.

Oyster RNA-seq Data Support the Development of Malacoherpesviridae Genomics

The family of double-stranded DNA (dsDNA) Malacoherpesviridae includes viruses able to infect marine mollusks and detrimental for worldwide aquaculture production. Due to fast-occurring mortality and a lack of permissive cell lines, the available data on the few known Malacoherpesviridae provide only partial support for the study of molecular virus features, life cycle, and evolutionary history. Following thorough data mining of bivalve and gastropod RNA-seq experiments, we used more than five million Malacoherpesviridae reads to improve the annotation of viral genomes and to characterize viral InDels, nucleotide stretches, and SNPs. Both genome and protein domain analyses confirmed the evolutionary diversification and gene uniqueness of known Malacoherpesviridae. However, the presence of Malacoherpesviridae-like sequences integrated within genomes of phylogenetically distant invertebrates indicates broad diffusion of these viruses and indicates the need for confirmatory investigations. The manifest co-occurrence of OsHV-1 genotype variants in single RNA-seq samples of Crassostrea gigas provide further support for the Malacoherpesviridae diversification. In addition to simple sequence motifs inter-punctuating viral ORFs, recombination-inducing sequences were found to be enriched in the OsHV-1 and AbHV1-AUS genomes. Finally, the highly correlated expression of most viral ORFs in multiple oyster samples is consistent with the burst of viral proteins during the lytic phase.

Long-Read Isoform Sequencing Reveals a Hidden Complexity of the Transcriptional Landscape of Herpes Simplex Virus Type 1

In this study, we used the amplified isoform sequencing technique from Pacific Biosciences to characterize the poly(A)⁺ fraction of the lytic transcriptome of the herpes simplex virus type 1 (HSV-1). Our analysis detected 34 formerly unidentified protein-coding genes, 10 non-coding RNAs, as well as 17 polycistronic and complex transcripts. This work also led us to identify many transcript isoforms, including 13 splice and 68 transcript end variants, as well as several transcript overlaps. Additionally, we determined previously unascertained transcriptional start and polyadenylation sites. We analyzed the transcriptional activity from the complementary DNA strand in five convergent HSV gene pairs with quantitative RT-PCR and detected antisense RNAs in each gene. This part of the study revealed an inverse correlation between the expressions of convergent partners. Our work adds new insights for understanding the complexity of the pervasive transcriptional overlaps by suggesting that there is a crosstalk between adjacent and distal genes through interaction between their transcription apparatuses. We also identified transcripts overlapping the HSV replication origins, which may indicate an interplay between the transcription and replication machineries. The relative abundance of HSV-1 transcripts has also been established by using a novel method based on the calculation of sequencing reads for the analysis.

Deletion of the us7 and us8 genes of pseudorabies virus exerts a differential effect on the expression of early and late viral genes

The pseudorabies virus (PRV; also known as Suid herpesvirus-1) is a neurotropic herpesvirus of swine. The us7 and us8 genes of this virus encode the glycoprotein I and E membrane proteins that form a heterodimer that is known to control cell-to-cell spread in tissue culture and in animals. In this study, we investigated the effect of the deletion of the PRV us7 and us8 genes on the genome-wide transcription and DNA replication using a multi-time-point quantitative reverse transcriptase-based real-time PCR technique. Abrogation of the us7/8 gene function was found to exert a drastic but differential effect on the expression of PRV genes during lytic infection. In the mutant virus, all kinetic classes of viral genes were significantly down-regulated at the first 6 h of infection, while having been upregulated later. The level of upregulation was the highest in the immediate-early (IE) and the early (E) genes; lower in the early-late (E/L) genes; and the lowest in the late (L) genes. The relative contribution of the L transcripts to the global transcriptome became lower, while the rest of the transcripts were expressed at a higher level in the mutant than in the wild-type virus.

Gene Expression Profiling with Cre-Conditional Pseudorabies Virus Reveals a Subset of Midbrain Neurons That Participate in Reward Circuitry

The mesolimbic dopamine pathway receives inputs from numerous regions of the brain as part of a neural system that detects rewarding stimuli and coordinates a behavioral response. The capacity to simultaneously map and molecularly define the components of this complex multisynaptic circuit would thus advance our understanding of the determinants of motivated behavior. To accomplish this, we have constructed pseudorabies virus (PRV) strains in which viral propagation and fluorophore expression are activated only after exposure to Cre recombinase. Once activated in Cre-expressing neurons, the virus serially labels chains of presynaptic neurons. Dual injection of GFP and mCherry tracing viruses simultaneously illuminates nigrostriatal and mesolimbic circuitry and shows no overlap, demonstrating that PRV transmission is confined to synaptically connected neurons. To molecularly profile mesolimbic dopamine neurons and their presynaptic inputs, we injected Cre-conditional GFP virus into the NAc of (anti-GFP) nanobody-L10 transgenic mice and immunoprecipitated translating ribosomes from neurons infected after retrograde tracing. Analysis of purified RNA revealed an enrichment of transcripts expressed in neurons of the dorsal raphe nuclei and lateral hypothalamus that project to the mesolimbic dopamine circuit. These studies identify important inputs to the mesolimbic dopamine pathway and further show that PRV circuit-directed translating ribosome affinity purification can be broadly applied to identify molecularly defined neurons comprising complex, multisynaptic circuits.SIGNIFICANCE STATEMENT The mesolimbic dopamine circuit integrates signals from key brain regions to detect and respond to rewarding stimuli. To further define this complex multisynaptic circuit, we constructed a panel of Cre recombinase-activated pseudorabies viruses (PRVs) that enabled retrograde tracing of neural inputs that terminate on Cre-expressing neurons. Using these viruses and Retro-TRAP (translating ribosome affinity purification), a previously reported molecular profiling method, we developed a novel technique that provides anatomic as well as molecular information about the neural components of polysynaptic circuits. We refer to this new method as PRV-Circuit-TRAP (PRV circuit-directed TRAP). Using it, we have identified major projections to the mesolimbic dopamine circuit from the lateral hypothalamus and dorsal raphe nucleus and defined a discrete subset of transcripts expressed in these projecting neurons, which will allow further characterization of this important pathway. Moreover, the method we report is general and can be applied to the study of other neural circuits.

Characterization of the Dynamic Transcriptome of a Herpesvirus with Long-read Single Molecule Real-Time Sequencing

Herpesvirus gene expression is co-ordinately regulated and sequentially ordered during productive infection. The viral genes can be classified into three distinct kinetic groups: immediate-early, early, and late classes. In this study, a massively parallel sequencing technique that is based on PacBio Single Molecule Real-time sequencing platform, was used for quantifying the poly(A) fraction of the lytic transcriptome of pseudorabies virus (PRV) throughout a 12-hour interval of productive infection on PK-15 cells. Other approaches, including microarray, real-time RT-PCR and Illumina sequencing are capable of detecting only the aggregate transcriptional activity of particular genomic regions, but not individual herpesvirus transcripts. However, SMRT sequencing allows for a distinction between transcript isoforms, including length- and splice variants, as well as between overlapping polycistronic RNA molecules. The non-amplified Isoform Sequencing (Iso-Seq) method was used to analyse the kinetic properties of the lytic PRV transcripts and to then classify them accordingly. Additionally, the present study demonstrates the general utility of long-read sequencing for the time-course analysis of global gene expression in practically any organism.

Molecular epidemiology of pseudorabies virus in Yunnan and the sequence analysis of its gD gene

Outbreaks of pseudorabies (PRs) have occurred in Yunnan, China, which caused significant economic loss. To determine the prevalence and origin of PR in Yunnan, especially among vaccinated pigs, overall 791 samples of blood, tissue, semen, and sera were analyzed by serological methods, PCR, and sequence analysis of gD gene. Detection with viral gI antibody or PCR showed that the yearly positive rates of PR virus (PRV) in Yunnan from 2010 to 2014 were 48.15, 21.26, 2.17, 5.22, and 0.35%, respectively, with an average of 15.43%. In general, the incidence declined through the period of 2010-2014 probably due to the application of PRV eradication strategies. A phylogenetic tree was constructed based on the complete sequence of gD gene, with all strains clustered into two independent clades, i.e., Asian and European-American clades. The virus isolates from Henan, Tianjin, Heilongjiang, Sichuan, Shandong, Fujian, Xinjiang, Hubei, Guangdong, and Yunnan fell into Asian group, which harbored South Korea isolate. Four Yunnan virus isolates together with South Korean Namyangju fell into in the European-American clade. It showed that PR was pandemic as there was not a clear clue about the geographical origin of the PRV isolates in China since 2010.

Full-Length Isoform Sequencing Reveals Novel Transcripts and Substantial Transcriptional Overlaps in a Herpesvirus

Whole transcriptome studies have become essential for understanding the complexity of genetic regulation. However, the conventionally applied short-read sequencing platforms cannot be used to reliably distinguish between many transcript isoforms. The Pacific Biosciences (PacBio) RS II platform is capable of reading long nucleic acid stretches in a single sequencing run. The pseudorabies virus (PRV) is an excellent system to study herpesvirus gene expression and potential interactions between the transcriptional units. In this work, non-amplified and amplified isoform sequencing protocols were used to characterize the poly(A⁺) fraction of the lytic transcriptome of PRV, with the aim of a complete transcriptional annotation of the viral genes. The analyses revealed a previously unrecognized complexity of the PRV transcriptome including the discovery of novel protein-coding and non-coding genes, novel mono- and polycistronic transcription units, as well as extensive transcriptional overlaps between neighboring and distal genes. This study identified non-coding transcripts overlapping all three replication origins of the PRV, which might play a role in the control of DNA synthesis. We additionally established the relative expression levels of gene products. Our investigations revealed that the whole PRV genome is utilized for transcription, including both DNA strands in all coding and intergenic regions. The genome-wide occurrence of transcript overlaps suggests a crosstalk between genes through a network formed by interacting transcriptional machineries with a potential function in the control of gene expression.

Functional Dissection of an Alternatively Spliced Herpesvirus Gene by Splice Site Mutagenesis

Herpesviruses have large and complex DNA genomes. The largest among the herpesviruses, those of the cytomegaloviruses, include over 170 genes. Although most herpesvirus gene products are expressed from unspliced transcripts, a substantial number of viral transcripts are spliced. Some viral transcripts are subject to alternative splicing, which leads to the expression of several proteins from a single gene. Functional analysis of individual proteins derived from an alternatively spliced gene is difficult, as deletion and nonsense mutagenesis, both common methods used in the generation of viral gene knockout mutants, affect several or all gene products at the same time. Here, we show that individual gene products of an alternatively spliced herpesvirus gene can be inactivated selectively by mutagenesis of the splice donor or acceptor site and by intron deletion or substitution mutagenesis. We used this strategy to dissect the essential M112/113 gene of murine cytomegalovirus (MCMV), which encodes the MCMV Early 1 (E1) proteins. The expression of each of the four E1 protein isoforms was inactivated individually, and the requirement for each isoform in MCMV replication was analyzed in fibroblasts, endothelial cells, and macrophages. We show that the E1 p87 isoform, but not the p33, p36, and p38 isoforms, is essential for viral replication in cell culture. Moreover, the presence of one of the two medium-size isoforms (p36 or p38) and the presence of intron 1, but not its specific sequence, are required for viral replication. This study demonstrates the usefulness of splice site mutagenesis for the functional analysis of alternatively spliced herpesvirus genes.

IMPORTANCE

Herpesviruses include up to 170 genes in their DNA genomes. The functions of most viral gene products remain poorly defined. The construction of viral gene knockout mutants has thus been an important tool for functional analysis of viral proteins. However, this strategy is of limited use when viral gene transcripts are alternatively spliced, leading to the expression of several proteins from a single gene. In this study, we showed, as a proof of principle, that each protein product of an alternatively spliced gene can be eliminated individually by splice site mutagenesis. Mutant viruses lacking individual protein products displayed different phenotypes, demonstrating that the products of alternatively spliced genes have nonredundant functions.