Faculty Bibliography

MOTIVATION/BACKGROUND:The chemical modification of ribonucleotides regulates the structure, stability, and interactions of RNAs. Profiling of these modifications using short-read (Illumina) sequencing techniques provides high sensitivity but low-to-medium resolution i.e., modifications cannot be assigned to specific transcript isoforms in regions of sequence overlap. An alternative strategy uses current fluctuations in nanopore-based long read direct RNA sequencing (DRS) to infer the location and identity of nucleotides that differ between two experimental conditions. While highly sensitive, these signal-level analyses require high quality transcriptome annotations and thus are best suited to the study of model organisms. By contrast, the detection of RNA modifications in microbial organisms which typically have no or low-quality annotations requires an alternative strategy. Here, we demonstrate that signal fluctuations directly influence error rates during base calling and thus provides an alternative approach for identifying modified nucleotides. RESULTS:DRUMMER (Detection of Ribonucleic acid Modifications Manifested in Error Rates (i) utilizes a range of statistical tests and background noise correction to identify modified nucleotides with high confidence, (ii) operates with similar sensitivity to signal-level analysis approaches, and (iii) correlates very well with orthogonal approaches. Using well-characterized DRS datasets supported by independent meRIP-Seq and miCLIP-Seq datasets we demonstrate that DRUMMER operates with high sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION/METHODS:DRUMMER is written in Python 3 and is available as open source in the GitHub repository: https://github.com/DepledgeLab/DRUMMER. SUPPLEMENTARY INFORMATION/BACKGROUND:Supplementary data are available at Bioinformatics online.

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.

Primary varicella-zoster virus (VZV) infection leads to varicella and the establishment of lifelong latency in sensory ganglion neurons. Reactivation of latent VZV causes herpes zoster, which is frequently associated with chronic pain. Latent viral gene expression is restricted to the VZV latency-associated transcript (VLT) and VLT-ORF63 (VLT63) fusion transcripts. Since VLT and VLT63 encode proteins that are expressed during lytic infection, we investigated whether pVLT and pVLT-ORF63 are essential for VZV replication by performing VZV genome mutagenesis using CRISPR/Cas9 and BAC technologies. We first established that CRISPR/Cas9 can efficiently mutate VZV genomes in lytically VZV-infected cells through targeting non-essential genes ORF8 and ORF11 and subsequently show recovery of viable mutant viruses. By contrast, the VLT region was markedly resistant to CRISPR/Cas9 editing. Whereas most mutants expressed wild-type or N-terminally altered versions of pVLT and pVLT-ORF63, only a minority of the resulting mutant viruses lacked pVLT and pVLT-ORF63 coding potential. Growth curve analysis showed that pVLT/pVLT-ORF63 negative viruses were viable, but impaired in growth in epithelial cells. We confirmed this phenotype independently using BAC-derived pVLT/pVLT-ORF63 negative and repaired viruses. Collectively, these data demonstrate that pVLT and/or pVLT-ORF63 are dispensable for lytic VZV replication but promote efficient VZV infection in epithelial cells.

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for dissecting the complexity of normal and diseased tissues, enabling characterization of cell diversity and heterogeneous phenotypic states in unprecedented detail. However, this technology has been underutilized for exploring the interactions between the host cell and viral pathogens in latently infected cells. Herein, we use scRNA-seq and single-molecule sensitivity fluorescent in situ hybridization (smFISH) technologies to investigate host single-cell transcriptome changes upon the reactivation of a human neurotropic virus, herpes simplex virus-1 (HSV-1). We identify the stress sensor growth arrest and DNA damage-inducible 45 beta (Gadd45b) as a critical antiviral host factor that regulates HSV-1 reactivation events in a subpopulation of latently infected primary neurons. We show that distinct subcellular localization of Gadd45b correlates with the viral late gene expression program, as well as the expression of the viral transcription factor, ICP4. We propose that a hallmark of a "successful" or "aborted" HSV-1 reactivation state in primary neurons is determined by a unique subcellular localization signature of the stress sensor Gadd45b.

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.

Primary infection with varicella-zoster virus (VZV) causes varicella and the establishment of lifelong latency in sensory ganglion neurons. In one-third of infected individuals VZV reactivates from latency to cause herpes zoster, often complicated by difficult-to-treat chronic pain. Experimental infection of non-human primates with simian varicella virus (SVV) recapitulates most features of human VZV disease, thereby providing the opportunity to study the pathogenesis of varicella and herpes zoster in vivo. However, compared to VZV, the transcriptome and the full coding potential of SVV remains incompletely understood. Here, we performed direct long-read RNA sequencing to annotate the SVV transcriptome in lytically SVV-infected African green monkey (AGM) and rhesus macaque (RM) kidney epithelial cells. We refined structures of canonical SVV transcripts and uncovered numerous RNA isoforms, splicing events, fusion transcripts and non-coding RNAs, mostly unique to SVV. We verified the expression of canonical and newly identified SVV transcripts in vivo, using lung samples from acutely SVV-infected cynomolgus macaques. Expression of selected transcript isoforms, including those located in the unique left-end of the SVV genome, was confirmed by reverse transcription PCR. Finally, we performed detailed characterization of the SVV homologue of the VZV latency-associated transcript (VLT), located antisense to ORF61. Analogous to VZV VLT, SVV VLT is multiply spliced and numerous isoforms are generated using alternative transcription start sites and extensive splicing. Conversely, low level expression of a single spliced SVV VLT isoform defines in vivo latency. Notably, the genomic location of VLT core exons is highly conserved between SVV and VZV. This work thus highlights the complexity of lytic SVV gene expression and provides new insights into the molecular biology underlying lytic and latent SVV infection. The identification of the SVV VLT homolog further underlines the value of the SVV non-human primate model to develop new strategies for prevention of herpes zoster.

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.

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.