Long-Read Sequencing of Human Cytomegalovirus Transcriptome Reveals RNA Isoforms Carrying Distinct Coding Potentials

Exploring the transcriptomic profile of human monkeypox virus via CAGE and native RNA sequencing approaches

In this study, we employed short- and long-read sequencing technologies to delineate the transcriptional architecture of the human monkeypox virus and to identify key regulatory elements that govern its gene expression. Specifically, we conducted a transcriptomic analysis to annotate the transcription start sites (TSSs) and transcription end sites (TESs) of the virus by utilizing Cap Analysis of gene expression sequencing on the Illumina platform and direct RNA sequencing on the Oxford Nanopore technology device. Our investigations uncovered significant complexity in the use of alternative TSSs and TESs in viral genes. In this research, we also detected the promoter elements and poly(A) signals associated with the viral genes. Additionally, we identified novel genes in both the left and right variable regions of the viral genome.IMPORTANCEGenerally, gaining insight into how the transcription of a virus is regulated offers insights into the key mechanisms that control its life cycle. The recent outbreak of the human monkeypox virus has underscored the necessity of understanding the basic biology of its causative agent. Our results are pivotal for constructing a comprehensive transcriptomic atlas of the human monkeypox virus, providing valuable resources for future studies.

Comprehensive bioinformatics analysis of human cytomegalovirus pathway genes in pan-cancer

BACKGROUND

Human cytomegalovirus (HCMV) is a herpesvirus that can infect various cell types and modulate host gene expression and immune response. It has been associated with the pathogenesis of various cancers, but its molecular mechanisms remain elusive.

METHODS

We comprehensively analyzed the expression of HCMV pathway genes across 26 cancer types using the Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) databases. We also used bioinformatics tools to study immune invasion and tumor microenvironment in pan-cancer. Cox regression and machine learning were used to analyze prognostic genes and their relationship with drug sensitivity.

RESULTS

We found that HCMV pathway genes are widely expressed in various cancers. Immune infiltration and the tumor microenvironment revealed that HCMV is involved in complex immune processes. We obtained prognostic genes for 25 cancers and significantly found 23 key genes in the HCMV pathway, which are significantly enriched in cellular chemotaxis and synaptic function and may be involved in disease progression. Notably, CaM family genes were up-regulated and AC family genes were down-regulated in most tumors. These hub genes correlate with sensitivity or resistance to various drugs, suggesting their potential as therapeutic targets.

CONCLUSIONS

Our study has revealed the role of the HCMV pathway in various cancers and provided insights into its molecular mechanism and therapeutic significance. It is worth noting that the key genes of the HCMV pathway may open up new doors for cancer prevention and treatment.

Viral and host network analysis of the human cytomegalovirus transcriptome in latency

HCMV genes UL135 and UL138 play opposing roles regulating latency and reactivation in CD34+ human progenitor cells (HPCs). Using the THP-1 cell line model for latency and reactivation, we designed an RNA sequencing study to compare the transcriptional profile of HCMV infection in the presence and absence of these genes. The loss of UL138 results in elevated levels of viral gene expression and increased differentiation of cell populations that support HCMV gene expression and genome synthesis. The loss of UL135 results in diminished viral gene expression during an initial burst that occurs as latency is established and no expression of eleven viral genes from the ULb' region even following stimulation for differentiation and reactivation. Transcriptional network analysis revealed host transcription factors with potential to regulate the ULb' genes in coordination with pUL135. These results reveal roles for UL135 and UL138 in regulation of viral gene expression and potentially hematopoietic differentiation.

Hepatitis B virus serum RNA transcript isoform composition and proportion in chronic hepatitis B patients by nanopore long-read sequencing

Introduction

Serum hepatitis B virus (HBV) RNA is a promising new biomarker to manage and predict clinical outcomes of chronic hepatitis B (CHB) infection. However, the HBV serum transcriptome within encapsidated particles, which is the biomarker analyte measured in serum, remains poorly characterized. This study aimed to evaluate serum HBV RNA transcript composition and proportionality by PCR-cDNA nanopore sequencing of samples from CHB patients having varied HBV genotype (gt, A to F) and HBeAg status.

Methods

Longitudinal specimens from 3 individuals during and following pregnancy (approximately 7 months between time points) were also investigated. HBV RNA extracted from 16 serum samples obtained from 13 patients (73.3% female, 84.6% Asian) was sequenced and serum HBV RNA isoform detection and quantification were performed using three bioinformatic workflows; FLAIR, RATTLE, and a GraphMap-based workflow within the Galaxy application. A spike-in RNA variant (SIRV) control mix was used to assess run quality and coverage. The proportionality of transcript isoforms was based on total HBV reads determined by each workflow.

Results

All chosen isoform detection workflows showed high agreement in transcript proportionality and composition for most samples. HBV pregenomic RNA (pgRNA) was the most frequently observed transcript isoform (93.8% of patient samples), while other detected transcripts included pgRNA spliced variants, 3' truncated variants and HBx mRNA, depending on the isoform detection method. Spliced variants of pgRNA were primarily observed in HBV gtB, C, E, or F-infected patients, with the Sp1 spliced variant detected most frequently. Twelve other pgRNA spliced variant transcripts were identified, including 3 previously unidentified transcripts, although spliced isoform identification was very dependent on the workflow used to analyze sequence data. Longitudinal sampling among pregnant and post-partum antiviral-treated individuals showed increasing proportions of 3' truncated pgRNA variants over time.

Conclusions

This study demonstrated long-read sequencing as a promising tool for the characterization of the serum HBV transcriptome. However, further studies are needed to better understand how serum HBV RNA isoform type and proportion are linked to CHB disease progression and antiviral treatment response.

Insights into the Transcriptome of Human Cytomegalovirus: A Comprehensive Review

Human cytomegalovirus (HCMV) is a widespread pathogen that poses significant risks to immunocompromised individuals. Its genome spans over 230 kbp and potentially encodes over 200 open-reading frames. The HCMV transcriptome consists of various types of RNAs, including messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs), with emerging insights into their biological functions. HCMV mRNAs are involved in crucial viral processes, such as viral replication, transcription, and translation regulation, as well as immune modulation and other effects on host cells. Additionally, four lncRNAs (RNA1.2, RNA2.7, RNA4.9, and RNA5.0) have been identified in HCMV, which play important roles in lytic replication like bypassing acute antiviral responses, promoting cell movement and viral spread, and maintaining HCMV latency. CircRNAs have gained attention for their important and diverse biological functions, including association with different diseases, acting as microRNA sponges, regulating parental gene expression, and serving as translation templates. Remarkably, HCMV encodes miRNAs which play critical roles in silencing human genes and other functions. This review gives an overview of human cytomegalovirus and current research on the HCMV transcriptome during lytic and latent infection.

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.

Human cytomegalovirus attenuates AKT activity by destabilizing insulin receptor substrate proteins

The phosphoinositide 3-kinase (PI3K)/AKT pathway plays crucial roles in cell viability and protein synthesis and is frequently co-opted by viruses to support their replication. Although many viruses maintain high levels of AKT activity during infection, other viruses, such as vesicular stomatitis virus and human cytomegalovirus (HCMV), cause AKT to accumulate in an inactive state. To efficiently replicate, HCMV requires FoxO transcription factors to localize to the infected cell nucleus (Zhang et. al. mBio 2022), a process directly antagonized by AKT. Therefore, we sought to investigate how HCMV inactivates AKT to achieve this. Subcellular fractionation and live cell imaging studies indicated that AKT failed to recruit to membranes upon serum-stimulation of infected cells. However, UV-inactivated virions were unable to render AKT non-responsive to serum, indicating a requirement for de novo viral gene expression. Interestingly, we were able to identify that UL38 (pUL38), a viral activator of mTORC1, is required to diminish AKT responsiveness to serum. mTORC1 contributes to insulin resistance by causing proteasomal degradation of insulin receptor substrate (IRS) proteins, such as IRS1, which are necessary for the recruitment of PI3K to growth factor receptors. In cells infected with a recombinant HCMV disrupted for UL38 , AKT responsiveness to serum is retained and IRS1 is not degraded. Furthermore, ectopic expression of UL38 in uninfected cells induces IRS1 degradation, inactivating AKT. These effects of UL38 were reversed by the mTORC1 inhibitor, rapamycin. Collectively, our results demonstrate that HCMV relies upon a cell-intrinsic negative feedback loop to render AKT inactive during productive infection.

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.

Deciphering the Potential Coding of Human Cytomegalovirus: New Predicted Transmembrane Proteome

CMV is a major cause of morbidity and mortality in immunocompromised individuals that will benefit from the availability of a vaccine. Despite the efforts made during the last decade, no CMV vaccine is available. An ideal CMV vaccine should elicit a broad immune response against multiple viral antigens including proteins involved in virus-cell interaction and entry. However, the therapeutic use of neutralizing antibodies targeting glycoproteins involved in viral entry achieved only partial protection against infection. In this scenario, a better understanding of the CMV proteome potentially involved in viral entry may provide novel candidates to include in new potential vaccine design. In this study, we aimed to explore the CMV genome to identify proteins with putative transmembrane domains to identify new potential viral envelope proteins. We have performed in silico analysis using the genome sequences of nine different CMV strains to predict the transmembrane domains of the encoded proteins. We have identified 77 proteins with transmembrane domains, 39 of which were present in all the strains and were highly conserved. Among the core proteins, 17 of them such as UL10, UL139 or US33A have no ascribed function and may be good candidates for further mechanistic studies.

Optimization of a Lambda-RED Recombination Method for Rapid Gene Deletion in Human Cytomegalovirus

Human cytomegalovirus (HCMV) continues to be a major cause of morbidity in transplant patients and newborns. However, the functions of many of the more than 282 genes encoded in the HCMV genome remain unknown. The development of bacterial artificial chromosome (BAC) technology contributes to the genetic manipulation of several organisms including HCMV. The maintenance of the HCMV BAC in E. coli cells permits the rapid generation of recombinant viral genomes that can be used to produce viral progeny in cell cultures for the study of gene function. We optimized the Lambda-Red Recombination system to construct HCMV gene deletion mutants rapidly in the complete set of tested genes. This method constitutes a useful tool that allows for the quick generation of a high number of gene deletion mutants, allowing for the analysis of the whole genome to improve our understanding of HCMV gene function. This may also facilitate the development of novel vaccines and therapeutics.

Time-Course Transcriptome Profiling of a Poxvirus Using Long-Read Full-Length Assay

Viral transcriptomes that are determined using first- and second-generation sequencing techniques are incomplete. Due to the short read length, these methods are inefficient or fail to distinguish between transcript isoforms, polycistronic RNAs, and transcriptional overlaps and readthroughs. Additionally, these approaches are insensitive for the identification of splice and transcriptional start sites (TSSs) and, in most cases, transcriptional end sites (TESs), especially in transcript isoforms with varying transcript ends, and in multi-spliced transcripts. Long-read sequencing is able to read full-length nucleic acids and can therefore be used to assemble complete transcriptome atlases. Although vaccinia virus (VACV) does not produce spliced RNAs, its transcriptome has a high diversity of TSSs and TESs, and a high degree of polycistronism that leads to enormous complexity. We applied single-molecule, real-time, and nanopore-based sequencing methods to investigate the time-lapse transcriptome patterns of VACV gene expression.

Combined nanopore and single-molecule real-time sequencing survey of human betaherpesvirus 5 transcriptome

Long-read sequencing (LRS), a powerful novel approach, is able to read full-length transcripts and confers a major advantage over the earlier gold standard short-read sequencing in the efficiency of identifying for example polycistronic transcripts and transcript isoforms, including transcript length- and splice variants. In this work, we profile the human cytomegalovirus transcriptome using two third-generation LRS platforms: the Sequel from Pacific BioSciences, and MinION from Oxford Nanopore Technologies. We carried out both cDNA and direct RNA sequencing, and applied the LoRTIA software, developed in our laboratory, for the transcript annotations. This study identified a large number of novel transcript variants, including splice isoforms and transcript start and end site isoforms, as well as putative mRNAs with truncated in-frame ORFs (located within the larger ORFs of the canonical mRNAs), which potentially encode N-terminally truncated polypeptides. Our work also disclosed a highly complex meshwork of transcriptional read-throughs and overlaps.

Insight for Immunotherapy of HCMV Infection

Human cytomegalovirus (HCMV), a ubiquitous in humans, has a high prevalence rate. Young people are susceptible to HCMV infection in developing countries, while older individuals are more susceptible in developed countries. Most patients have no obvious symptoms from the primary infection. Studies have indicated that the virus has gradually adapted to the host immune system. Therefore, the control of HCMV infection requires strong immune modulation. With the recent advances in immunotherapy, its application to HCMV infections is receiving increasing attention. Here, we discuss the immune response to HCMV infection, the immune escape mechanism, and the different roles that HCMV plays in various types of immunotherapy, including vaccines, adoptive cell therapy, checkpoint blockade therapy, and targeted antibodies.

Time course profiling of host cell response to herpesvirus infection using nanopore and synthetic long-read transcriptome sequencing

Third-generation sequencing is able to read full-length transcripts and thus to efficiently identify RNA molecules and transcript isoforms, including transcript length and splice isoforms. In this study, we report the time-course profiling of the effect of bovine alphaherpesvirus type 1 on the gene expression of bovine epithelial cells using direct cDNA sequencing carried out on MinION device of Oxford Nanopore Technologies. These investigations revealed a substantial up- and down-regulatory effect of the virus on several gene networks of the host cells, including those that are associated with antiviral response, as well as with viral transcription and translation. Additionally, we report a large number of novel bovine transcript isoforms identified by nanopore and synthetic long-read sequencing. This study demonstrates that viral infection causes differential expression of host transcript isoforms. We could not detect an increased rate of transcriptional readthroughs as described in another alphaherpesvirus. According to our knowledge, this is the first report on the use of LoopSeq for the analysis of eukaryotic transcriptomes. This is also the first report on the application of nanopore sequencing for the kinetic characterization of cellular transcriptomes. This study also demonstrates the utility of nanopore sequencing for the characterization of dynamic transcriptomes in any organisms.

Time-course transcriptome analysis of host cell response to poxvirus infection using a dual long-read sequencing approach

OBJECTIVE

In this study, we applied two long-read sequencing (LRS) approaches, including single-molecule real-time and nanopore-based sequencing methods to investigate the time-lapse transcriptome patterns of host gene expression as a response to Vaccinia virus infection. Transcriptomes determined using short-read sequencing approaches are incomplete because these platforms are inefficient or fail to distinguish between polycistronic RNAs, transcript isoforms, transcriptional start sites, as well as transcriptional readthroughs and overlaps. Long-read sequencing is able to read full-length nucleic acids and can therefore be used to assemble complete transcriptome atlases.

RESULTS

In this work, we identified a number of novel transcripts and transcript isoforms of Chlorocebus sabaeus. Additionally, analysis of the most abundant 768 host transcripts revealed a significant overrepresentation of the class of genes in the "regulation of signaling receptor activity" Gene Ontology annotation as a result of viral infection.

Immune dysregulation and autoreactivity correlate with disease severity in SARS-CoV-2-associated multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV-2 infection. We profiled MIS-C, adult COVID-19, and healthy pediatric and adult individuals using single-cell RNA sequencing, flow cytometry, antigen receptor repertoire analysis, and unbiased serum proteomics, which collectively identified a signature in MIS-C patients that correlated with disease severity. Despite having no evidence of active infection, MIS-C patients had elevated S100A-family alarmins and decreased antigen presentation signatures, indicative of myeloid dysfunction. MIS-C patients showed elevated expression of cytotoxicity genes in NK and CD8+ T cells and expansion of specific IgG-expressing plasmablasts. Clinically severe MIS-C patients displayed skewed memory T cell TCR repertoires and autoimmunity characterized by endothelium-reactive IgG. The alarmin, cytotoxicity, TCR repertoire, and plasmablast signatures we defined have potential for application in the clinic to better diagnose and potentially predict disease severity early in the course of MIS-C.

Genetic Variability of Human Cytomegalovirus Clinical Isolates Correlates With Altered Expression of Natural Killer Cell-Activating Ligands and IFN-γ

Human cytomegalovirus (HCMV) infection often leads to systemic disease in immunodeficient patients and congenitally infected children. Despite its clinical significance, the exact mechanisms contributing to HCMV pathogenesis and clinical outcomes have yet to be determined. One of such mechanisms involves HCMV-mediated NK cell immune response, which favors viral immune evasion by hindering NK cell-mediated cytolysis. This process appears to be dependent on the extent of HCMV genetic variation as high levels of variability in viral genes involved in immune escape have an impact on viral pathogenesis. However, the link between viral genome variations and their functional effects has so far remained elusive. Thus, here we sought to determine whether inter-host genetic variability of HCMV influences its ability to modulate NK cell responses to infection. For this purpose, five HCMV clinical isolates from a previously characterized cohort of pediatric patients with confirmed HCMV congenital infection were evaluated by next-generation sequencing (NGS) for genetic polymorphisms, phylogenetic relationships, and multiple-strain infection. We report variable levels of genetic characteristics among the selected clinical strains, with moderate variations in genome regions associated with modulation of NK cell functions. Remarkably, we show that different HCMV clinical strains differentially modulate the expression of several ligands for the NK cell-activating receptors NKG2D, DNAM-1/CD226, and NKp30. Specifically, the DNAM-1/CD226 ligand PVR/CD155 appears to be predominantly upregulated by fast-replicating (“aggressive”) HCMV isolates. On the other hand, the NGK2D ligands ULBP2/5/6 are downregulated regardless of the strain used, while other NK cell ligands (i.e., MICA, MICB, ULBP3, Nectin-2/CD112, and B7-H6) are not significantly modulated. Furthermore, we show that IFN-γ; production by NK cells co-cultured with HCMV-infected fibroblasts is directly proportional to the aggressiveness of the HCMV clinical isolates employed. Interestingly, loss of NK cell-modulating genes directed against NK cell ligands appears to be a common feature among the “aggressive” HCMV strains, which also share several gene variants across their genomes. Overall, even though further studies based on a higher number of patients would offer a more definitive scenario, our findings provide novel mechanistic insights into the impact of HCMV genetic variability on NK cell-mediated immune responses.

Post-infectious inflammatory disease in MIS-C features elevated cytotoxicity signatures and autoreactivity that correlates with severity

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV2 infection in otherwise healthy children. Here, we define immune abnormalities in MIS-C compared to adult COVID-19 and pediatric/adult healthy controls using single-cell RNA sequencing, antigen receptor repertoire analysis, unbiased serum proteomics, and in vitro assays. Despite no evidence of active infection, we uncover elevated S100A-family alarmins in myeloid cells and marked enrichment of serum proteins that map to myeloid cells and pathways including cytokines, complement/coagulation, and fluid shear stress in MIS-C patients. Moreover, NK and CD8 T cell cytotoxicity genes are elevated, and plasmablasts harboring IgG1 and IgG3 are expanded. Consistently, we detect elevated binding of serum IgG from severe MIS-C patients to activated human cardiac microvascular endothelial cells in culture. Thus, we define immunopathology features of MIS-C with implications for predicting and managing this SARS-CoV2-induced critical illness in children.

Interplay Between CMGC Kinases Targeting SR Proteins and Viral Replication: Splicing and Beyond

Protein phosphorylation constitutes a major post-translational modification that critically regulates the half-life, intra-cellular distribution, and activity of proteins. Among the large number of kinases that compose the human kinome tree, those targeting RNA-binding proteins, in particular serine/arginine-rich (SR) proteins, play a major role in the regulation of gene expression by controlling constitutive and alternative splicing. In humans, these kinases belong to the CMGC [Cyclin-dependent kinases (CDKs), Mitogen-activated protein kinases (MAPKs), Glycogen synthase kinases (GSKs), and Cdc2-like kinases (CLKs)] group and several studies indicate that they also control viral replication via direct or indirect mechanisms. The aim of this review is to describe known and emerging activities of CMGC kinases that share the common property to phosphorylate SR proteins, as well as their interplay with different families of viruses, in order to advance toward a comprehensive knowledge of their pro- or anti-viral phenotype and better assess possible translational opportunities.

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.

Long-read assays shed new light on the transcriptome complexity of a viral pathogen

Characterization of global transcriptomes using conventional short-read sequencing is challenging due to the insensitivity of these platforms to transcripts isoforms, multigenic RNA molecules, and transcriptional overlaps. Long-read sequencing (LRS) can overcome these limitations by reading full-length transcripts. Employment of these technologies has led to the redefinition of transcriptional complexities in reported organisms. In this study, we applied LRS platforms from Pacific Biosciences and Oxford Nanopore Technologies to profile the vaccinia virus (VACV) transcriptome. We performed cDNA and direct RNA sequencing analyses and revealed an extremely complex transcriptional landscape of this virus. In particular, VACV genes produce large numbers of transcript isoforms that vary in their start and termination sites. A significant fraction of VACV transcripts start or end within coding regions of neighbouring genes. This study provides new insights into the transcriptomic profile of this viral pathogen.

Meta-analytic approach for transcriptome profiling of herpes simplex virus type 1

In this meta-analysis, we re-analysed and compared herpes simplex virus type 1 transcriptomic data generated by eight studies using various short- and long-read sequencing techniques and different library preparation methods. We identified a large number of novel mRNAs, non-coding RNAs and transcript isoforms, and validated many previously published transcripts. Here, we present the most complete HSV-1 transcriptome to date. Furthermore, we also demonstrate that various sequencing techniques, including both cDNA and direct RNA sequencing approaches, are error-prone, which can be circumvented by using integrated approaches. This work draws attention to the need for using multiple sequencing approaches and meta-analyses in transcriptome profiling studies to obtain reliable results.

High-resolution annotation of the mouse preimplantation embryo transcriptome using long-read sequencing

The transcriptome of the preimplantation mouse embryo has been previously annotated by short-read sequencing, with limited coverage and accuracy. Here we utilize a low-cell number transcriptome based on the Smart-seq2 method to perform long-read sequencing. Our analysis describes additional novel transcripts and complexity of the preimplantation transcriptome, identifying 2280 potential novel transcripts from previously unannotated loci and 6289 novel splicing isoforms from previously annotated genes. Notably, these novel transcripts and isoforms with transcription start sites are enriched for an active promoter modification, H3K4me3. Moreover, we generate a more complete and precise transcriptome by combining long-read and short-read data during early embryogenesis. Based on this approach, we identify a previously undescribed isoform of Kdm4dl with a modified mRNA reading frame and a novel noncoding gene designated XLOC_004958. Depletion of Kdm4dl or XLOC_004958 led to abnormal blastocyst development. Thus, our data provide a high-resolution and more precise transcriptome during preimplantation mouse embryogenesis.

The Requirement for US28 During Cytomegalovirus Latency Is Independent of US27 and US29 Gene Expression

The ability to establish a latent infection with periodic reactivation events ensures herpesviruses, like human cytomegalovirus (HCMV), lifelong infection, and serial passage. The host-pathogen relationship throughout HCMV latency is complex, though both cellular and viral factors influence the equilibrium between latent and lytic infection. We and others have shown one of the viral-encoded G protein-coupled receptors, US28, is required for HCMV latency. US28 potentiates signals both constitutively and in response to ligand binding, and we previously showed deletion of the ligand binding domain or mutation of the G protein-coupling domain results in the failure to maintain latency similar to deletion of the entire US28 open reading frame (ORF). Interestingly, a recent publication detailed an altered phenotype from that previously reported, showing US28 is required for viral reactivation rather than latency, suggesting the US28 ORF deletion impacts transcription of the surrounding genes. Here, we show an independently generated US28-stop mutant, like the US28 ORF deletion mutant, fails to maintain latency in hematopoietic cells. Further, we found US27 and US29 transcription in each of these mutants was comparable to their expression during wild type infection, suggesting neither US28 mutant alters mRNA levels of the surrounding genes. Finally, infection with a US28 ORF deletion virus expressed US27 protein comparable to its expression following wild type infection. In sum, our new data strongly support previous findings from our lab and others, detailing a requirement for US28 during HCMV latent infection.

Bright and Early: Inhibiting Human Cytomegalovirus by Targeting Major Immediate-Early Gene Expression or Protein Function

The human cytomegalovirus (HCMV), one of eight human herpesviruses, establishes lifelong latent infections in most people worldwide. Primary or reactivated HCMV infections cause severe disease in immunosuppressed patients and congenital defects in children. There is no vaccine for HCMV, and the currently approved antivirals come with major limitations. Most approved HCMV antivirals target late molecular processes in the viral replication cycle including DNA replication and packaging. “Bright and early” events in HCMV infection have not been exploited for systemic prevention or treatment of disease. Initiation of HCMV replication depends on transcription from the viral major immediate-early (IE) gene. Alternative transcripts produced from this gene give rise to the IE1 and IE2 families of viral proteins, which localize to the host cell nucleus. The IE1 and IE2 proteins are believed to control all subsequent early and late events in HCMV replication, including reactivation from latency, in part by antagonizing intrinsic and innate immune responses. Here we provide an update on the regulation of major IE gene expression and the functions of IE1 and IE2 proteins. We will relate this insight to experimental approaches that target IE gene expression or protein function via molecular gene silencing and editing or small chemical inhibitors.

Comprehensive annotations of human herpesvirus 6A and 6B genomes reveal novel and conserved genomic features

Human herpesvirus-6 (HHV-6) A and B are ubiquitous betaherpesviruses, infecting the majority of the human population. They encompass large genomes and our understanding of their protein coding potential is far from complete. Here, we employ ribosome-profiling and systematic transcript-analysis to experimentally define HHV-6 translation products. We identify hundreds of new open reading frames (ORFs), including upstream ORFs (uORFs) and internal ORFs (iORFs), generating a complete unbiased atlas of HHV-6 proteome. By integrating systematic data from the prototypic betaherpesvirus, human cytomegalovirus, we uncover numerous uORFs and iORFs conserved across betaherpesviruses and we show uORFs are enriched in late viral genes. We identified three highly abundant HHV-6 encoded long non-coding RNAs, one of which generates a non-polyadenylated stable intron appearing to be a conserved feature of betaherpesviruses. Overall, our work reveals the complexity of HHV-6 genomes and highlights novel features conserved between betaherpesviruses, providing a rich resource for future functional studies.

Functional Profile of Human Cytomegalovirus Genes and Their Associated Diseases: A Review

The human cytomegalovirus (HCMV), whose genome is 235 ± 1.9 kbp long, is a common herpesvirus. However, the functions of many of its genes are still unknown. HCMV is closely associated with various human diseases and infects 60-90% of the global population. It can infect various human cells, including fibroblasts, epithelial cells, endothelial cells, smooth muscle cells, and monocytes. Although HCMV infection is generally asymptomatic and causes subtle clinical symptoms, it can generate a robust immune response and establish a latent infection in immunocompromised individuals, including those with AIDS, transplant recipients, and developing fetuses. Currently available antivirals approved for the treatment of HCMV-associated diseases are limited by dose-limiting toxicity and the emergence of resistance; however, vaccines and immunoglobulins are unavailable. In this review, we have summarized the recent literature on 43 newly identified HCMV genes. We have described their novel functions on the viral replication cycle, latency, and host immune evasion. Further, we have discussed HCMV-associated diseases and current therapeutic targets. Our review may provide a foundational basis for studies aiming to prevent and develop targeted therapies for HCMV-associated diseases.

Human cytomegalovirus interactome analysis identifies degradation hubs, domain associations and viral protein functions

Human cytomegalovirus (HCMV) extensively modulates host cells, downregulating >900 human proteins during viral replication and degrading ≥133 proteins shortly after infection. The mechanism of degradation of most host proteins remains unresolved, and the functions of many viral proteins are incompletely characterised. We performed a mass spectrometry-based interactome analysis of 169 tagged, stably-expressed canonical strain Merlin HCMV proteins, and two non-canonical HCMV proteins, in infected cells. This identified a network of >3400 virus-host and >150 virus-virus protein interactions, providing insights into functions for multiple viral genes. Domain analysis predicted binding of the viral UL25 protein to SH3 domains of NCK Adaptor Protein-1. Viral interacting proteins were identified for 31/133 degraded host targets. Finally, the uncharacterised, non-canonical ORFL147C protein was found to interact with elements of the mRNA splicing machinery, and a mutational study suggested its importance in viral replication. The interactome data will be important for future studies of herpesvirus infection.

Template-switching artifacts resemble alternative polyadenylation

BACKGROUND

Alternative polyadenylation is commonly examined using cDNA sequencing, which is known to be affected by template-switching artifacts. However, the effects of such template-switching artifacts on alternative polyadenylation are generally disregarded, while alternative polyadenylation artifacts are attributed to internal priming.

RESULTS

Here, we analyzed both long-read cDNA sequencing and direct RNA sequencing data of two organisms, generated by different sequencing platforms. We developed a filtering algorithm which takes into consideration that template-switching can be a source of artifactual polyadenylation when filtering out spurious polyadenylation sites. The algorithm outperformed the conventional internal priming filters based on comparison to direct RNA sequencing data. We also showed that the polyadenylation artifacts arise in cDNA sequencing at consecutive stretches of as few as three adenines. There was no substantial difference between the lengths of poly(A) tails at the artifactual and the true transcriptional end sites even though it is expected that internal priming artifacts have shorter poly(A) tails than genuine polyadenylated reads.

CONCLUSIONS

Our findings suggest that template switching plays an important role in the generation of spurious polyadenylation and support the need for more rigorous filtering of artifactual polyadenylation sites in cDNA data, or that alternative polyadenylation should be annotated using native RNA sequencing.

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.

Alternative promoters drive human cytomegalovirus reactivation from latency

Reactivation from latency requires reinitiation of viral gene expression and culminates in the production of infectious progeny. The major immediate early promoter (MIEP) of human cytomegalovirus (HCMV) drives the expression of crucial lytic cycle transactivators but is silenced during latency in hematopoietic progenitor cells (HPCs). Because the MIEP has poor activity in HPCs, it is unclear how viral transactivators are expressed during reactivation. It has been presumed that viral gene expression is reinitiated via de-repression of the MIEP. We demonstrate that immediate early transcripts arising from reactivation originate predominantly from alternative promoters within the canonical major immediate early locus. Disruption of these intronic promoters results in striking defects in re-expression of viral genes and viral genome replication in the THP-1 latency model. Furthermore, we show that these promoters are necessary for efficient reactivation in primary CD34+ HPCs. Our findings shift the paradigm for HCMV reactivation by demonstrating that promoter switching governs reactivation from viral latency in a context-specific manner.

Identification and characterization of a novel group of natural anti-sense transcripts from RNA1.2 gene locus of human cytomegalovirus

Background:

Natural anti-sense transcripts (NATs), which are transcribed from the complementary DNA strand of annotated genes, exert regulatory function of gene expression. Increasing studies recognized anti-sense transcription widespread throughout human cytomegalovirus (HCMV) genome, whereas the anti-sense transcription of RNA1.2 gene locus has never been investigated. In this study, the transcription of the RNA1.2 anti-sense strand was investigated in clinically isolated HCMV strain.

Methods:

Strand-specific high-through RNA-sequencing (RNA-seq) was performed to find possible anti-sense transcripts (ASTs). For analyzing and visualization of RNA-seq data sets, Integrative Genomics Viewer software was applied. To confirm these possibilities, Northern blotting and rapid amplification of cDNA ends (RACE) were used.

Results:

Transcription of the opposite strand of RNA1.2 gene locus was detected by RNA-sequencing using RNAs extracted from human embryonic lung fibroblasts infected with HCMV clinical isolate HAN. At least three HCMV NATs, named RNA1.2 AST 1, RNA1.2 AST2, and RNA1.2 AST3, were characterized by Northern blotting and RACE analyses. These RNA1.2 ASTs orientated from the complementary strand of RNA1.2 locus during the late phase of HCMV infection. The 5′- and 3′-termini of these transcripts were located within the opposite sequence of the predicted RNA1.2 gene.

Conclusion:

A cluster of novel NATs was transcribed from the opposite sequence of the HCMV RNA1.2 gene region.

Hidden regulation of herpes simplex virus 1 pre-mRNA splicing and polyadenylation by virally encoded immediate early gene ICP27

In contrast to human cells, very few HSV-1 genes are known to be spliced, although the same pre-mRNA processing machinery is shared. Here, through global analysis of splice junctions in cells infected with HSV-1 and an HSV-1 mutant virus with deletion of infectious cell culture protein 27 (ICP27), one of two viral immediate early (IE) genes essential for viral replication, we identify hundreds of novel alternative splice junctions mapping to both previously known HSV-1 spliced genes and previously unknown spliced genes, the majority of which alter the coding potential of viral genes. Quantitative and qualitative splicing efficiency analysis of these novel alternatively spliced genes based on RNA-Seq and RT-PCR reveals that splicing at these novel splice sites is efficient only when ICP27 is absent; while in wildtype HSV-1 infected cells, the splicing of these novel splice junctions is largely silenced in a gene/sequence specific manner, suggesting that ICP27 not only promotes accumulation of ICP27 targeted transcripts but also ensures correctness of the functional coding sequences through inhibition of alternative splicing. Furthermore, ICP27 toggles expression of ICP34.5, the major viral neurovirulence factor, through inhibition of splicing and activation of a proximal polyadenylation signal (PAS) in the newly identified intron, revealing a novel regulatory mechanism for expression of a viral gene. Thus, through the viral IE protein ICP27, HSV-1 co-opts both splicing and polyadenylation machinery to achieve optimal viral gene expression during lytic infection. On the other hand, during latent infection when ICP27 is absent, HSV-1 likely takes advantages of host splicing machinery to restrict expression of randomly activated antigenic viral genes to achieve immune evasion.

Little is known regarding to how HSV, a large DNA virus and known to contain very few spliced genes, escapes host pre-mRNA splicing machinery. Here, by establishing a high throughput splice junction identification platform and quantitative analysis method to assess splicing efficiency based on high throughput data, we find that HSV-1 encodes hundreds of previously unknown alternative splice junctions; however, splicing of these novel spliced genes is largely silenced in wild-type HSV-1 infected cells, explaining why only very few spliced genes have been previously identified in HSV-1. Moreover, ICP27 is required for splicing inhibition and 3’ end formation of ICP34.5, the major viral neurovirulence factor and also the major target of latently expressed viral miRNAs. These findings not only fundamentally change the view of HSV gene structure, but also reveal a mechanism by which HSV employs host splicing and polyadenylation machineries to achieve optimal gene expression during acute infection and may also contribute to immune evasion during latency when ICP27 is not expressed.

Cytomegalovirus: Shape-Shifting the Immune System

Systems-based based approaches have begun to shed light on extrinsic factors that contribute to immune system variation. Among these, CMV (HHV-5, a β-herpesvirus) imposes a surprisingly profound impact. Most of the world's population is CMV⁺, and the virus goes through three distinct infection phases en route to establishing lifelong détente with its host. Immune control of CMV in each phase recruits unique arms of host defense, and in turn the virus employs multiple immune-modulatory strategies that help facilitate the establishment of lifelong persistence. In this review, we explain how CMV shapes immunity and discuss the impact it may have on overall health.

Cellular Cullin RING Ubiquitin Ligases: Druggable Host Dependency Factors of Cytomegaloviruses

Human cytomegalovirus (HCMV) is a ubiquitous betaherpesvirus that frequently causes morbidity and mortality in individuals with insufficient immunity, such as transplant recipients, AIDS patients, and congenitally infected newborns. Several antiviral drugs are approved to treat HCMV infections. However, resistant HCMV mutants can arise in patients receiving long-term therapy. Additionally, side effects and the risk to cause birth defects limit the use of currently approved antivirals against HCMV. Therefore, the identification of new drug targets is of clinical relevance. Recent work identified DNA-damage binding protein 1 (DDB1) and the family of the cellular cullin (Cul) RING ubiquitin (Ub) ligases (CRLs) as host-derived factors that are relevant for the replication of human and mouse cytomegaloviruses. The first-in-class CRL inhibitory compound Pevonedistat (also called MLN4924) is currently under investigation as an anti-tumor drug in several clinical trials. Cytomegaloviruses exploit CRLs to regulate the abundance of viral proteins, and to induce the proteasomal degradation of host restriction factors involved in innate and intrinsic immunity. Accordingly, pharmacological blockade of CRL activity diminishes viral replication in cell culture. In this review, we summarize the current knowledge concerning the relevance of DDB1 and CRLs during cytomegalovirus replication and discuss chances and drawbacks of CRL inhibitory drugs as potential antiviral treatment against HCMV.

Human cytomegalovirus genomics and transcriptomics through the lens of next-generation sequencing: revision and future challenges

The human cytomegalovirus (HCMV) genome was sequenced by hierarchical shotgun almost 30 years ago. Over these years, low and high passaged strains have been sequenced, improving, albeit still far from complete, the understanding of the coding potential, expression dynamics and diversity of wild-type HCMV strains. Next-generation sequencing (NGS) platforms have enabled a huge advancement, facilitating the comparison of differentially passaged strains, challenging diagnostics and research based on a single or reduced gene set genotyping. In addition, it allowed to link genetic features to different viral phenotypes as for example, correlating large genomic re-arrangements to viral attenuation or different mutations to antiviral resistance and cell tropism. NGS platforms provided the first high-resolution experiments to HCMV dynamics, allowing the study of intra-host viral population structures and the description of rare transcriptional events. Long-read sequencing has recently become available, helping to identify new genomic re-arrangements, partially accounting for the genetic variability displayed in clinical isolates, as well as, in changing the understanding of the HCMV transcriptome. Better knowledge of the transcriptome resulted in a vast number of new splicing events and alternative transcripts, although most of them still need additional validation. This review summarizes the sequencing efforts reached so far, discussing its approaches and providing a revision and new nuances on HCMV sequence variability in the sequencing field.

Interactions between the transcription and replication machineries regulate the RNA and DNA synthesis in the herpesviruses

The temporal coordination of viral gene expression is imperative for the regulation of the herpesvirus replication cycle. While the main factors of this transcriptional coordination are known, the subtler control mechanisms of gene expression remain elusive. Recent long read sequencing-based approached have revealed an intricate meshwork of overlaps between the herpesvirus transcripts and the overlap of the replication origins with noncoding RNAs. It has been shown that the transcriptional apparatuses can physically interfere with one another while transcribing overlapping regions. We hypothesize that transcriptional interference regulates the global gene expression across the herpesvirus genome. Additionally, an overall decrease in transcriptional activity in individual viral genes has been observed following the onset of DNA replication. An overlap of the replication origins with specific transcripts has also been described in several herpesviruses. The genome-wide interactions between the transcriptional apparatuses and between the replication and transcriptional machineries suggest the existence of novel layers of genetic regulation.

Going the Distance: Optimizing RNA-Seq Strategies for Transcriptomic Analysis of Complex Viral Genomes

Transcriptome profiling has become routine in studies of many biological processes. However, the favored approaches such as short-read Illumina RNA sequencing are giving way to long-read sequencing platforms better suited to interrogating the complex transcriptomes typical of many RNA and DNA viruses. Here, we provide a guide-tailored to molecular virologists-to the ins and outs of viral transcriptome sequencing and discuss the strengths and weaknesses of the major RNA sequencing technologies as tools to analyze the abundance and diversity of the viral transcripts made during infection.

Transcriptome-wide analysis of a baculovirus using nanopore sequencing

Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is a prototypic baculovirus infecting specific insects. AcMNPV contains a large double-stranded DNA genome encoding a complex transcriptome. This virus has a widespread application as a vector for the expression of heterologous proteins. Here, we present a dataset, derived from Oxford Nanopore Technologies (ONT) long-read sequencing platform. We used both cDNA and direct RNA sequencing techniques. The dataset contains 520,310 AcMNPV and 1,309,481 host cell reads using the regular cDNA-sequencing method of ONT technique, whereas altogether 6,456 reads were produced by using direct RNA-sequencing. We also used a Cap-selection protocol for certain ONT samples, and obtained 2,568,669 reads by using this method. The raw reads were aligned to the AcMNPV reference genome (KM667940.1). Here, we openly released the ‘static’ and the dynamic transcript catalogue of AcMNPV. This dataset can be used for deep analyses of the transcriptomic and epitranscriptomic patterns of the AcMNPV and the host cell. The data can be also useful for the validation of different bioinformatics software packages and analysis tools.

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.

Dynamic transcriptome profiling dataset of vaccinia virus obtained from long-read sequencing techniques

Background

Poxviruses are large DNA viruses that infect humans and animals. Vaccinia virus (VACV) has been applied as a live vaccine for immunization against smallpox, which was eradicated by 1980 as a result of worldwide vaccination. VACV is the prototype of poxviruses in the investigation of the molecular pathogenesis of the virus. Short-read sequencing methods have revolutionized transcriptomics; however, they are not efficient in distinguishing between the RNA isoforms and transcript overlaps. Long-read sequencing (LRS) is much better suited to solve these problems and also allow direct RNA sequencing. Despite the scientific relevance of VACV, no LRS data have been generated for the viral transcriptome to date.

Findings

For the deep characterization of the VACV RNA profile, various LRS platforms and library preparation approaches were applied. The raw reads were mapped to the VACV reference genome and also to the host (Chlorocebus sabaeus) genome. In this study, we applied the Pacific Biosciences RSII and Sequel platforms, which altogether resulted in 937,531 mapped reads of inserts (1.42 Gb), while we obtained 2,160,348 aligned reads (1.75 Gb) from the different library preparation methods using the MinION device from Oxford Nanopore Technologies.

Conclusions

By applying cutting-edge technologies, we were able to generate a large dataset that can serve as a valuable resource for the investigation of the dynamic VACV transcriptome, the virus-host interactions, and RNA base modifications. These data can provide useful information for novel gene annotations in the VACV genome. Our dataset can also be used to analyze the currently available LRS platforms, library preparation methods, and bioinformatics pipelines.

New therapies for human cytomegalovirus infections

The recent approval of letermovir marks a new era of therapy for human cytomegalovirus (HCMV) infections, particularly for the prevention of HCMV disease in hematopoietic stem cell transplant recipients. For almost 30 years ganciclovir has been the therapy of choice for these infections and by today's standards this drug exhibits only modest antiviral activity that is often insufficient to completely suppress viral replication, and drives the selection of drug-resistant variants that continue to replicate and contribute to disease. While ganciclovir remains the therapy of choice, additional drugs that inhibit novel molecular targets, such as letermovir, will be required as highly effective combination therapies are developed not only for the treatment of immunocompromised hosts, but also for congenitally infected infants. Sustained efforts, largely in the biotech industry and academia, have identified additional highly active lead compounds that have progressed into clinical studies with varying levels of success and at least two have the potential to be approved in the near future. Some of the new drugs in the pipeline inhibit new molecular targets, remain effective against isolates that have developed resistance to existing therapies, and promise to augment existing therapeutic regimens. Here, we will describe some of the unique features of HCMV biology and discuss their effect on therapeutic needs. Existing drugs will also be discussed and some of the more promising candidates will be reviewed with an emphasis on those progressing through clinical studies. The in vitro and in vivo antiviral activity, spectrum of antiviral activity, and mechanism of action of new compounds will be reviewed to provide an update on potential new therapies for HCMV infections that have progressed significantly in recent years.

Expression Levels of Glycoprotein O (gO) Vary between Strains of Human Cytomegalovirus, Influencing the Assembly of gH/gL Complexes and Virion Infectivity

The tropism of human cytomegalovirus (HCMV) is influenced by the envelope glycoprotein complexes gH/gL/gO and gH/gL/UL128-131. During virion assembly, gO and the UL128-131 proteins compete for binding to gH/gL in the endoplasmic reticulum (ER). This assembly process clearly differs among strains, since Merlin (ME) virions contain abundant gH/gL/UL128-131 and little gH/gL/gO, whereas TR contains much higher levels of total gH/gL, mostly in the form of gH/gL/gO, but much lower levels of gH/gL/UL128-131 than ME. Remaining questions include (i) what are the mechanisms behind these assembly differences, and (ii) do differences reflect in vitro culture adaptations or natural genetic variations? Since the UL74(gO) open reading frame (ORF) differs in 25% of amino acids between TR and ME, we analyzed recombinant viruses in which the UL74(gO) ORF was swapped. TR virions were >40-fold more infectious than ME. Transcriptional repression of UL128-131 enhanced the infectivity of ME to the level of TR, despite still far lower levels of gH/gL/gO. Swapping the UL74(gO) ORF had no effect on either TR or ME. A quantitative immunoprecipitation approach revealed that gH/gL expression levels were within 4-fold between TR and ME, but the gO expression level was 20-fold lower for ME, which suggested differences in mRNA transcription, translation, or rapid ER-associated degradation of gO. trans-Complementation of gO expression during ME replication gave a 6-fold enhancement of infectivity beyond the 40-fold effect of UL128-131 repression alone. Overall, strain variations in the assembly of gH/gL complexes result from differences in the expression of gO and UL128-131, and selective advantages for reduced UL128-131 expression during fibroblast propagation are much stronger than those for higher gO expression.IMPORTANCE Specific genetic differences between independently isolated HCMV strains may result from purifying selection on de novo mutations arising during propagation in culture or random sampling among the diversity of genotypes present in clinical specimens. Results presented indicate that while reduced UL128-131 expression may confer a powerful selective advantage during cell-free propagation of HCMV in fibroblast cultures, selective pressures for increased gO expression are much weaker. Thus, variation in gO expression among independent strains may represent natural genotype variability present in vivo This may have important implications for virus-host interactions, such as immune recognition, and underscores the value of studying molecular mechanisms of replication using multiple HCMV strains.

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.

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.

Long-read sequencing of the human cytomegalovirus transcriptome with the Pacific Biosciences RSII platform

Long-read RNA sequencing allows for the precise characterization of full-length transcripts, which makes it an indispensable tool in transcriptomics. The human cytomegalovirus (HCMV) genome has been first sequenced in 1989 and although short-read sequencing studies have uncovered much of the complexity of its transcriptome, only few of its transcripts have been fully annotated. We hereby present a long-read RNA sequencing dataset of HCMV infected human lung fibroblast cells sequenced by the Pacific Biosciences RSII platform. Seven SMRT cells were sequenced using oligo(dT) primers to reverse transcribe poly(A)-selected RNA molecules and one library was prepared using random primers for the reverse transcription of the rRNA-depleted sample. Our dataset contains 122,636 human and 33,086 viral (HMCV strain Towne) reads. The described data include raw and processed sequencing files, and combined with other datasets, they can be used to validate transcriptome analysis tools, to compare library preparation methods, to test base calling algorithms or to identify genetic variants.