Faculty Bibliography

N⁶-methyladenosine (m⁶A) is an abundant internal RNA modification, influencing transcript fate and function in uninfected and virus-infected cells. Installation of m⁶A by the nuclear RNA methyltransferase METTL3 occurs cotranscriptionally; however, the genomes of some cytoplasmic RNA viruses are also m⁶A-modified. How the cellular m⁶A modification machinery impacts coronavirus replication, which occurs exclusively in the cytoplasm, is unknown. Here we show that replication of SARS-CoV-2, the agent responsible for the COVID-19 pandemic, and a seasonal human β-coronavirus HCoV-OC43, can be suppressed by depletion of METTL3 or cytoplasmic m⁶A reader proteins YTHDF1 and YTHDF3 and by a highly specific small molecule METTL3 inhibitor. Reduction of infectious titer correlates with decreased synthesis of viral RNAs and the essential nucleocapsid (N) protein. Sites of m⁶A modification on genomic and subgenomic RNAs of both viruses were mapped by methylated RNA immunoprecipitation sequencing (meRIP-seq). Levels of host factors involved in m⁶A installation, removal, and recognition were unchanged by HCoV-OC43 infection; however, nuclear localization of METTL3 and cytoplasmic m⁶A readers YTHDF1 and YTHDF2 increased. This establishes that coronavirus RNAs are m⁶A-modified and host m⁶A pathway components control β-coronavirus replication. Moreover, it illustrates the therapeutic potential of targeting the m⁶A pathway to restrict coronavirus reproduction.

The contributions of the viral component of the microbiome-the virome-to the development of innate and adaptive immunity are largely unknown. Here, we systematically defined the host response in mice to a panel of eukaryotic enteric viruses representing six different families. Infections with most of these viruses were asymptomatic in the mice, the magnitude and duration of which was dependent on the microbiota. Flow cytometric and transcriptional profiling of mice mono-associated with these viruses unveiled general adaptations by the host, such as lymphocyte differentiation and IL-22 signatures in the intestine, as well as numerous viral-strain-specific responses that persisted. Comparison with a dataset derived from analogous bacterial mono-association in mice identified bacterial species that evoke an immune response comparable with the viruses we examined. These results expand an understanding of the immune space occupied by the enteric virome and underscore the importance of viral exposure events.

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm³ and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

Members of the family Herpesviridae have enveloped, spherical virions with characteristic complex structures consisting of symmetrical and non-symmetrical components. The linear, double-stranded DNA genomes of 125-241 kbp contain 70-170 genes, of which 43 have been inherited from an ancestral herpesvirus. In general, herpesviruses have coevolved with and are highly adapted to their hosts, which comprise many mammalian, avian and reptilian species. Following primary infection, they are able to establish lifelong latent infection, during which there is limited viral gene expression. Severe disease is usually observed only in the foetus, the very young, the immunocompromised or following infection of an alternative host. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Herpesviridae, which is available at ictv.global/report/herpesviridae.

Varicella-zoster virus (VZV) establishes lifelong neuronal latency in most humans world-wide, reactivating in one-third to cause herpes zoster and occasionally chronic pain. How VZV establishes, maintains and reactivates from latency is largely unknown. VZV transcription during latency is restricted to the latency-associated transcript (VLT) and RNA 63 (encoding ORF63) in naturally VZV-infected human trigeminal ganglia (TG). While significantly more abundant, VLT levels positively correlated with RNA 63 suggesting co-regulated transcription during latency. Here, we identify VLT-ORF63 fusion transcripts and confirm VLT-ORF63, but not RNA 63, expression in human TG neurons. During in vitro latency, VLT is transcribed, whereas VLT-ORF63 expression is induced by reactivation stimuli. One isoform of VLT-ORF63, encoding a fusion protein combining VLT and ORF63 proteins, induces broad viral gene transcription. Collectively, our findings show that VZV expresses a unique set of VLT-ORF63 transcripts, potentially involved in the transition from latency to lytic VZV infection.

Adenovirus is a nuclear replicating DNA virus reliant on host RNA processing machinery. Processing and metabolism of cellular RNAs can be regulated by METTL3, which catalyzes the addition of N6-methyladenosine (m⁶A) to mRNAs. While m⁶A-modified adenoviral RNAs have been previously detected, the location and function of this mark within the infectious cycle is unknown. Since the complex adenovirus transcriptome includes overlapping spliced units that would impede accurate m⁶A mapping using short-read sequencing, here we profile m⁶A within the adenovirus transcriptome using a combination of meRIP-seq and direct RNA long-read sequencing to yield both nucleotide and transcript-resolved m⁶A detection. Although both early and late viral transcripts contain m⁶A, depletion of m⁶A writer METTL3 specifically impacts viral late transcripts by reducing their splicing efficiency. These data showcase a new technique for m⁶A discovery within individual transcripts at nucleotide resolution, and highlight the role of m⁶A in regulating splicing of a viral pathogen.

Varicella-zoster virus (VZV), a double-stranded DNA virus, causes varicella, establishes lifelong latency in ganglionic neurons, and reactivates later in life to cause herpes zoster, commonly associated with chronic pain. The VZV genome is densely packed and produces multitudes of overlapping transcripts deriving from both strands. While 71 distinct open reading frames (ORFs) have thus far been experimentally defined, the full coding potential of VZV remains unknown. Here, we integrated multiple short-read RNA sequencing approaches with long-read direct RNA sequencing on RNA isolated from VZV-infected cells to provide a comprehensive reannotation of the lytic VZV transcriptome architecture. Through precise mapping of transcription start sites, splice junctions, and polyadenylation sites, we identified 136 distinct polyadenylated VZV RNAs that encode canonical ORFs, noncanonical ORFs, and ORF fusions, as well as putative noncoding RNAs (ncRNAs). Furthermore, we determined the kinetic class of all VZV transcripts and observed, unexpectedly, that transcripts encoding the ORF62 protein, previously designated Immediate-Early, were expressed with Late kinetics. Our work showcases the complexity of the VZV transcriptome and provides a comprehensive resource that will facilitate future functional studies of coding RNAs, ncRNAs, and the biological mechanisms underlying the regulation of viral transcription and translation during lytic VZV infection.IMPORTANCE Transcription from herpesviral genomes, executed by the host RNA polymerase II and regulated by viral proteins, results in coordinated viral gene expression to efficiently produce infectious progeny. However, the complete coding potential and regulation of viral gene expression remain ill-defined for the human alphaherpesvirus varicella-zoster virus (VZV), causative agent of both varicella and herpes zoster. Here, we present a comprehensive overview of the VZV transcriptome and the kinetic class of all identified viral transcripts, using two virus strains and two biologically relevant cell types. Additionally, our data provide an overview of how VZV diversifies its transcription from one of the smallest herpesviral genomes. Unexpectedly, the transcript encoding the major viral transactivator protein (pORF62) was expressed with Late kinetics, whereas orthologous transcripts in other alphaherpesviruses are typically expressed during the immediate early phase. Therefore, our work both establishes the architecture of the VZV transcriptome and provides insight into regulation of alphaherpesvirus gene expression.

The varicella-zoster virus (VZV) genome, comprises both unique and repeated regions. The genome also includes reiteration regions, designated R1 to R5, which are tandemly repeating sequences termed elements. These regions represent an understudied feature of the VZV genome. The R4 region is duplicated, with one copy in the internal repeat short (IRs) which we designated R4A and a second copy in the terminal repeat short (TRs) termed R4B. We developed primers to amplify and Sanger sequence these regions, including independent amplification of both R4 regions. Reiteration regions from >80 cases of PCR-confirmed shingles were sequenced and analyzed. Complete genome sequences for the remaining portions of these viruses were determined using Illumina MiSeq. We identified 28 elements not previously reported, including at least one element for each R region. Length heterogeneity was substantial in R3, R4A and R4B. Length heterogeneity between the two copies of R4 was common.

The genomes of DNA viruses encode deceptively complex transcriptomes evolved to maximize coding potential within the confines of a relatively small genome. Defining the full range of viral RNAs produced during an infection is key to understanding the viral replication cycle and its interactions with the host cell. Traditional short-read (Illumina) sequencing approaches are problematic in this setting due to the difficulty of assigning short reads to individual RNAs in regions of transcript overlap and to the biases introduced by the required recoding and amplification steps. Additionally, different methodologies may be required to analyze the 5' and 3' ends of RNAs, which increases both cost and effort. The advent of long-read nanopore sequencing simplifies this approach by providing a single assay that captures and sequences full length RNAs, either in cDNA or native RNA form. The latter is particularly appealing as it reduces known recoding biases whilst allowing more advanced analyses such as estimation of poly(A) tail length and the detection of RNA modifications including N⁶ -methyladenosine. Using herpes simplex virus (HSV)-infected primary fibroblasts as a template, we provide a step-by-step guide to the production of direct RNA sequencing libraries suitable for sequencing using Oxford Nanopore Technologies platforms and provide a simple computational approach to deriving a high-quality annotation of the HSV transcriptome from the resulting sequencing data. © 2020 by John Wiley & Sons, Inc. Basic Protocol 1: Productive infection of primary fibroblasts with herpes simplex virus Support Protocol: Cell passage and plating of primary fibroblasts Basic Protocol 2: Preparation and sequencing of dRNA-seq libraries from virus-infected cells Basic Protocol 3: Processing, alignment, and analysis of dRNA-seq datasets.