Faculty Bibliography

UNLABELLED:Nanopore direct RNA sequencing (DRS) enables the capture and full-length sequencing of native RNAs, without recoding or amplification bias. Resulting data sets may be interrogated to define the identity and location of chemically modified ribonucleotides, as well as the length of poly(A) tails, on individual RNA molecules. The success of these analyses is highly dependent on the provision of high-resolution transcriptome annotations in combination with workflows that minimize misalignments and other analysis artifacts. Existing software solutions for generating high-resolution transcriptome annotations are poorly suited to small gene-dense genomes of viruses due to the challenge of identifying distinct transcript isoforms where alternative splicing and overlapping RNAs are prevalent. To resolve this, we identified key characteristics of DRS data sets that inform resulting read alignments and developed the nanopore guided annotation of transcriptome architectures (NAGATA) software package (https://github.com/DepledgeLab/NAGATA). We demonstrate, using a combination of synthetic and original DRS data sets derived from adenoviruses, herpesviruses, coronaviruses, and human cells, that NAGATA outperforms existing transcriptome annotation software and yields a consistently high level of precision and recall when reconstructing both gene sparse and gene-dense transcriptomes. Finally, we apply NAGATA to generate the first high-resolution transcriptome annotation of the neglected pathogen human adenovirus type F41 (HAdV-41) for which we identify 77 distinct transcripts encoding at least 23 different proteins. IMPORTANCE/OBJECTIVE:The transcriptome of an organism denotes the full repertoire of encoded RNAs that may be expressed. This is critical to understanding the biology of an organism and for accurate transcriptomic and epitranscriptomic-based analyses. Annotating transcriptomes remains a complex task, particularly in small gene-dense organisms such as viruses which maximize their coding capacity through overlapping RNAs. To resolve this, we have developed a new software nanopore guided annotation of transcriptome architectures (NAGATA) which utilizes nanopore direct RNA sequencing (DRS) datasets to rapidly produce high-resolution transcriptome annotations for diverse viruses and other organisms.

Most individuals are latently infected with herpes simplex virus type 1 (HSV-1), and it is well-established that HSV-1 establishes latency in sensory neurons of peripheral ganglia. However, it was recently proposed that latent HSV-1 is also present in immune cells recovered from the ganglia of experimentally infected mice. Here, we reanalyzed the single-cell RNA sequencing (scRNA-Seq) data that formed the basis for that conclusion. Unexpectedly, off-target priming in 3' scRNA-Seq experiments enabled the detection of non-polyadenylated HSV-1 latency-associated transcript (LAT) intronic RNAs. However, LAT reads were near-exclusively detected in mixed populations of cells undergoing cell death. Specific loss of HSV-1 LAT and neuronal transcripts during quality control filtering indicated widespread destruction of neurons, supporting the presence of contaminating cell-free RNA in other cells following tissue processing. In conclusion, the reported detection of latent HSV-1 in non-neuronal cells is best explained using compromised scRNA-Seq datasets.IMPORTANCEMost people are infected with herpes simplex virus type 1 (HSV-1) during their life. Once infected, the virus generally remains in a latent (silent) state, hiding within the neurons of peripheral ganglia. Periodic reactivation (reawakening) of the virus may cause fresh diseases such as cold sores. A recent study using single-cell RNA sequencing (scRNA-Seq) proposed that HSV-1 can also establish latency in the immune cells of mice, challenging existing dogma. We reanalyzed the data from that study and identified several flaws in the methodologies and analyses performed that invalidate the published conclusions. Specifically, we showed that the methodologies used resulted in widespread destruction of neurons which resulted in the presence of contaminants that confound the data analysis. We thus conclude that there remains little to no evidence for HSV-1 latency in immune cells.

This chapter first details the structure, organization and coding content of the VZV genome to provide a foundation on which the molecular evolution of the virus can be projected. We subsequently describe the evolution of molecular profiling approaches from restriction fragment length polymorphisms to single nucleotide polymorphism profiling to modern day high-throughput sequencing approaches. We describe how the application of these methodologies led to our current model of VZV phylogeograpy including the number and structure of geographic clades and the role of recombination in reshaping these.

Adenovirus is a common human pathogen that relies on host cell processes for transcription and processing of viral RNA and protein production. Although adenoviral promoters, splice junctions, and polyadenylation sites have been characterized using low-throughput biochemical techniques or short read cDNA-based sequencing, these technologies do not fully capture the complexity of the adenoviral transcriptome. By combining Illumina short-read and nanopore long-read direct RNA sequencing approaches, we mapped transcription start sites and RNA cleavage and polyadenylation sites across the adenovirus genome. In addition to confirming the known canonical viral early and late RNA cassettes, our analysis of splice junctions within long RNA reads revealed an additional 35 novel viral transcripts that meet stringent criteria for expression. These RNAs include fourteen new splice junctions which lead to expression of canonical open reading frames (ORFs), six novel ORF-containing transcripts, and 15 transcripts encoding for messages that could alter protein functions through truncation or fusion of canonical ORFs. In addition, we detect RNAs that bypass canonical cleavage sites and generate potential chimeric proteins by linking distinct gene transcription units. Among these chimeric proteins we detected an evolutionarily conserved protein containing the N-terminus of E4orf6 fused to the downstream DBP/E2A ORF. Loss of this novel protein, E4orf6/DBP, was associated with aberrant viral replication center morphology and poor viral spread. Our work highlights how long-read sequencing technologies combined with mass spectrometry can reveal further complexity within viral transcriptomes and resulting proteomes.

Attenuation of the live varicella Oka vaccine (vOka) has been attributed to mutations in the genome acquired during cell culture passage of pOka (parent strain); however, the precise mechanisms of attenuation remain unknown. Comparative sequence analyses of several vaccine batches showed that over 100 single-nucleotide polymorphisms (SNPs) are conserved across all vaccine batches; 6 SNPs are nearly fixed, suggesting that these SNPs are responsible for attenuation. By contrast, prior analysis of chimeric vOka and pOka recombinants indicates that loci other than these six SNPs contribute to attenuation. Here, we report that pOka consists of a heterogenous population of virus sequences with two nearly equally represented bases, guanine (G) or adenine (A), at nucleotide 2096 of the ORF31 coding sequence, which encodes glycoprotein B (gB) resulting in arginine (R) or glutamine (Q), respectively, at amino acid 699 of gB. By contrast, 2096A/699Q is dominant in vOka (>99.98%). gB699Q/gH/gL showed significantly less fusion activity than gB699R/gH/gL in a cell-based fusion assay. Recombinant pOka with gB669Q (rpOka_gB699Q) had a similar growth phenotype as vOka during lytic infection in cell culture including human primary skin cells; however, rpOka_gB699R showed a growth phenotype similar to pOka. rpOka_gB699R entered neurons from axonal terminals more efficiently than rpOka_gB699Q in the presence of cell membrane-derived vesicles containing gB. Strikingly, when a mixture of pOka with both alleles equally represented was used to infect human neurons from axon terminals, pOka with gB699R was dominant for virus entry. These results identify a variant allele in gB that contributes to attenuation of vOka. IMPORTANCE The live-attenuated varicella vaccine has reduced the burden of chickenpox. Despite its development in 1974, the molecular basis for its attenuation is still not well understood. Since the live-attenuated varicella vaccine is the only licensed human herpesvirus vaccine that prevents primary disease, it is important to understand the mechanism for its attenuation. Here we identify that a variant allele in glycoprotein B (gB) selected during generation of the varicella vaccine contributes to its attenuation. This variant is impaired for fusion, virus entry into neurons from nerve terminals, and replication in human skin cells. Identification of a variant allele in gB, one of the essential herpesvirus core genes, that contributes to its attenuation may provide insights that assist in the development of other herpesvirus vaccines.

Understanding the molecular mechanisms of herpes simplex virus 1 (HSV-1) latent infection and reactivation in neurons requires the use of in vitro model systems. Establishing a quiescent infection in cultured neurons is problematic, as any infectious virus released can superinfect the cultures. Previous studies have used the viral DNA replication inhibitor acyclovir to prevent superinfection and promote latency establishment. Data from these previous models have shown that reactivation is biphasic, with an initial phase I expression of all classes of lytic genes, which occurs independently of histone demethylase activity and viral DNA replication but is dependent on the cell stress protein DLK. Here, we describe a new model system using HSV-1 Stayput-GFP, a reporter virus that is defective for cell-to-cell spread and establishes latent infections without the need for acyclovir. The establishment of a latent state requires a longer time frame than previous models using DNA replication inhibitors. This results in a decreased ability of the virus to reactivate using established inducers, and as such, a combination of reactivation triggers is required. Using this system, we demonstrate that biphasic reactivation occurs even when latency is established in the absence of acyclovir. Importantly, phase I lytic gene expression still occurs in a histone demethylase and viral DNA replication-independent manner and requires DLK activity. These data demonstrate that the two waves of viral gene expression following HSV-1 reactivation are independent of secondary infection and not unique to systems that require acyclovir to promote latency establishment. IMPORTANCE Herpes simplex virus-1 (HSV-1) enters a latent infection in neurons and periodically reactivates. Reactivation manifests as a variety of clinical symptoms. Studying latency and reactivation in vitro is invaluable, allowing the molecular mechanisms behind both processes to be targeted by therapeutics that reduce the clinical consequences. Here, we describe a novel in vitro model system using a cell-to-cell spread-defective HSV-1, known as Stayput-GFP, which allows for the study of latency and reactivation at the single neuron level. We anticipate this new model system will be an incredibly valuable tool for studying the establishment and reactivation of HSV-1 latent infection in vitro. Using this model, we find that initial reactivation events are dependent on cellular stress kinase DLK but independent of histone demethylase activity and viral DNA replication. Our data therefore further validate the essential role of DLK in mediating a wave of lytic gene expression unique to reactivation.

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.

The chemical modification of ribonucleotides plays an integral role in the biology of diverse viruses and their eukaryotic host cells. Mapping the precise identity, location, and abundance of modified ribonucleotides remains a key goal of many studies aimed at characterizing the function and importance of a given modification. While mapping of specific RNA modifications through short-read sequencing approaches has powered a wealth of new discoveries in the past decade, this approach is limited by inherent biases and an absence of linkage information. Moreover, in viral contexts, the challenge is increased due to the compact nature of viral genomes giving rise to many overlapping transcript isoforms that cannot be adequately resolved using short-read sequencing approaches. The recent emergence of nanopore sequencing, specifically the ability to directly sequence native RNAs from virus-infected host cells, provides not just a new methodology for mapping modified ribonucleotides but also a new conceptual framework for what can be derived from the resulting sequencing data. In this minireview, we provide a detailed overview of how nanopore direct RNA sequencing works, the computational approaches applied to identify modified ribonucleotides, and the core concepts underlying both. We further highlight recent studies that have applied this approach to interrogating viral biology and finish by discussing key experimental considerations and how we predict that these methodologies will continue to evolve.

Regulated loading of eIF3-bound 40S ribosomes on capped mRNA is generally dependent upon the translation initiation factor eIF4E; however, mRNA translation often proceeds during physiological stress, such as virus infection, when eIF4E availability and activity are limiting. It remains poorly understood how translation of virus and host mRNAs are regulated during infection stress. While initially sensitive to mTOR inhibition, which limits eIF4E-dependent translation, we show that protein synthesis in human cytomegalovirus (HCMV)-infected cells unexpectedly becomes progressively reliant upon eIF3d. Targeting eIF3d selectively inhibits HCMV replication, reduces polyribosome abundance, and interferes with expression of essential virus genes and a host gene expression signature indicative of chronic ER stress that fosters HCMV reproduction. This reveals a strategy whereby cellular eIF3d-dependent protein production is hijacked to exploit virus-induced ER stress. Moreover, it establishes how switching between eIF4E and eIF3d-responsive cap-dependent translation can differentially tune virus and host gene expression in infected cells.