Faculty Bibliography

Epstein-Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host's innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of world-wide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main sub-populations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV.

During latent infections with herpes simplex virus 1 (HSV-1), viral transcription is restricted and the genomes are mostly maintained in silenced chromatin, whereas in lytically infected cells all viral genes are transcribed and the genomes are dynamically chromatinized. Histones in the viral chromatin bear markers of silenced chromatin at early times in lytic infection or of active transcription at later times. The virion protein VP16 activates transcription of the immediate-early (IE) genes by recruiting transcription activators and chromatin remodelers to their promoters. Two IE proteins, ICP0 and ICP4 which modulate chromatin epigenetics, then activate transcription of early and late genes. Although chromatin is involved in the mechanism of activation of HSV- transcription, its precise role is not entirely understood. In the cellular genome, chromatin dynamics often modulate transcription competence whereas promoter-specific transcription factors determine transcription activity. Here, biophysical fractionation of serially digested HSV-1 chromatin followed by short-read deep sequencing indicates that nuclear HSV-1 DNA has different biophysical properties than protein-free or encapsidated HSV-1 DNA. The entire HSV-1 genomes in infected cells were equally accessible. The accessibility of transcribed or non-transcribed genes under any given condition did not differ, and each gene was entirely sampled in both the most and least accessible chromatin. However, HSV-1 genomes fractionated differently under conditions of generalized or restricted transcription. Approximately 1/3 of the HSV-1 DNA including fully sampled genes resolved to the most accessible chromatin when HSV-1 transcription was active, but such enrichment was reduced to only 3% under conditions of restricted HSV-1 transcription. Short sequences of restricted accessibility separated genes with different transcription levels. Chromatin dynamics thus provide a first level of regulation on HSV-1 transcription, dictating the transcriptional competency of the genomes during lytic infections, whereas the transcription of individual genes is then most likely activated by specific transcription factors. Moreover, genes transcribed to different levels are separated by short sequences with limited accessibility.

The Glasgow s17 syn+ strain of herpes simplex virus 1 (HSV1) is arguably the best characterized strain and has provided the reference sequence for HSV1 genetic studies. Here we show that our original s17 syn+ stock was a mixed population from which we have isolated a minor variant that, unlike other strains in the laboratory, fails to be efficiently released from infected cells and spreads predominantly by direct cell-to-cell transmission. Analysis of other s17-derived viruses that had been isolated elsewhere revealed a number with the same release phenotype. Second-generation sequencing of 8 plaque-purified s17-derived viruses revealed sequences that vary by 50 single-nucleotide polymorphisms (SNPs), including approximately 10 coding SNPs. This compared to interstrain variations of around 800 SNPs in strain Sc16, of which a quarter were coding changes. Amongst the variations found within s17, we identified 13 variants of glycoprotein C within the original stock of virus that were predominantly a consequence of altered homopolymeric runs of C residues. Characterization of seven isolates coding for different forms of gC indicated that all were expressed, despite six of them lacking a transmembrane domain. While the release phenotype did not correlate directly with any of these identified gC variations, further demonstration that nine clinical isolates of HSV1 also fail to spread through extracellular release raises the possibility that propagation in tissue culture had altered the HSV1 s17 transmission phenotype. Hence, the s17 intrastrain variation identified here offers an excellent model for understanding both HSV1 transmission and tissue culture adaptation.

Recent sequencing efforts have led to estimates of human cytomegalovirus (HCMV) genome-wide intrahost diversity that rival those of persistent RNA viruses [Renzette N, Bhattacharjee B, Jensen JD, Gibson L, Kowalik TF (2011) PLoS Pathog 7:e1001344]. Here, we deep sequence HCMV genomes recovered from single and longitudinally collected blood samples from immunocompromised children to show that the observations of high within-host HCMV nucleotide diversity are explained by the frequent occurrence of mixed infections caused by genetically distant strains. To confirm this finding, we reconstructed within-host viral haplotypes from short-read sequence data. We verify that within-host HCMV nucleotide diversity in unmixed infections is no greater than that of other DNA viruses analyzed by the same sequencing and bioinformatic methods and considerably less than that of human immunodeficiency and hepatitis C viruses. By resolving individual viral haplotypes within patients, we reconstruct the timing, likely origins, and natural history of superinfecting strains. We uncover evidence for within-host recombination between genetically distinct HCMV strains, observing the loss of the parental virus containing the nonrecombinant fragment. The data suggest selection for strains containing the recombinant fragment, generating testable hypotheses about HCMV evolution and pathogenesis. These results highlight that high HCMV diversity present in some samples is caused by coinfection with multiple distinct strains and provide reassurance that within the host diversity for single-strain HCMV infections is no greater than for other herpesviruses.

Characterizing complex viral transcriptomes by conventional RNA sequencing approaches is complicated by high gene density, overlapping reading frames, and complex splicing patterns. Direct RNA sequencing (direct RNA-seq) using nanopore arrays offers an exciting alternative whereby individual polyadenylated RNAs are sequenced directly, without the recoding and amplification biases inherent to other sequencing methodologies. Here we use direct RNA-seq to profile the herpes simplex virus type 1 (HSV-1) transcriptome during productive infection of primary cells. We show how direct RNA-seq data can be used to define transcription initiation and RNA cleavage sites associated with all polyadenylated viral RNAs and demonstrate that low level read-through transcription produces a novel class of chimeric HSV-1 transcripts, including a functional mRNA encoding a fusion of the viral E3 ubiquitin ligase ICP0 and viral membrane glycoprotein L. Thus, direct RNA-seq offers a powerful method to characterize the changing transcriptional landscape of viruses with complex genomes.

Transcriptome profiling has become routine in studies of many biological processes. However, 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 viral transcripts made during infection.

BACKGROUND:Human herpesvirus-6A and -6B (HHV-6) are betaherpesviruses that reach > 90% seroprevalence in the adult population. Unique among human herpesviruses, HHV-6 can integrate into the subtelomeric regions of human chromosomes; when this occurs in germ line cells it causes a condition called inherited chromosomally integrated HHV-6 (iciHHV-6). Only two complete genomes are available for replicating HHV-6B, leading to numerous conflicting annotations and little known about the global genomic diversity of this ubiquitous virus. RESULTS:Using a custom capture panel for HHV-6B, we report complete genomes from 61 isolates of HHV-6B from active infections (20 from Japan, 35 from New York state, and 6 from Uganda), and 64 strains of iciHHV-6B (mostly from North America). HHV-6B sequence clustered by geography and illustrated extensive recombination. Multiple iciHHV-6B sequences from unrelated individuals across the United States were found to be completely identical, consistent with a founder effect. Several iciHHV-6B strains clustered with strains from recent active pediatric infection. Combining our genomic analysis with the first RNA-Seq and shotgun proteomics studies of HHV-6B, we completely reannotated the HHV-6B genome, altering annotations for more than 10% of existing genes, with multiple instances of novel splicing and genes that hitherto had gone unannotated. CONCLUSION:Our results are consistent with a model of intermittent de novo integration of HHV-6B into host germline cells during active infection with a large contribution of founder effect in iciHHV-6B. Our data provide a significant advance in the genomic annotation of HHV-6B, which will contribute to the detection, diversity, and control of this virus.

Meeting Report on the 8th Annual Symposium of the Colorado Alphaherpesvirus Latency Society (CALS), held on May 16-19, 2018, in Vail, Colorado.

Varicella-zoster virus (VZV), an alphaherpesvirus, establishes lifelong latent infection in the neurons of >90% humans worldwide, reactivating in one-third to cause shingles, debilitating pain and stroke. How VZV maintains latency remains unclear. Here, using ultra-deep virus-enriched RNA sequencing of latently infected human trigeminal ganglia (TG), we demonstrate the consistent expression of a spliced VZV mRNA, antisense to VZV open reading frame 61 (ORF61). The spliced VZV latency-associated transcript (VLT) is expressed in human TG neurons and encodes a protein with late kinetics in productively infected cells in vitro and in shingles skin lesions. Whereas multiple alternatively spliced VLT isoforms (VLT_ly) are expressed during lytic infection, a single unique VLT isoform, which specifically suppresses ORF61 gene expression in co-transfected cells, predominates in latently VZV-infected human TG. The discovery of VLT links VZV with the other better characterized human and animal neurotropic alphaherpesviruses and provides insights into VZV latency.

Background/UNASSIGNED:Adenoviruses are significant pathogens for the immunocompromised, arising from primary infection or reinfection. Serotyping is insufficient to support nosocomial transmission investigations. We investigate whether whole-genome sequencing (WGS) provides clinically relevant information on transmission among patients in a pediatric tertiary hospital. Methods/UNASSIGNED:We developed a target-enriched adenovirus WGS technique for clinical samples and retrospectively sequenced 107 adenovirus-positive residual diagnostic samples, including viremias (>5 × 104 copies/mL), from 37 patients collected January 2011-March 2016. Whole-genome sequencing was used to determine genotype and for phylogenetic analysis. Results/UNASSIGNED:Adenovirus sequences were recovered from 105 of 107 samples. Full genome sequences were recovered from all 20 nonspecies C samples and from 36 of 85 species C viruses, with partial genome sequences recovered from the rest. Whole-genome phylogenetic analysis suggested linkage of 3 genotype A31 cases and uncovered an unsuspected epidemiological link to an A31 infection first detected on the same ward 4 years earlier. In 9 samples from 1 patient who died, we identified a mixed genotype adenovirus infection. Conclusions/UNASSIGNED:Adenovirus WGS from clinical samples is possible and useful for genotyping and molecular epidemiology. Whole-genome sequencing identified likely nosocomial transmission with greater resolution than conventional genotyping and distinguished between adenovirus disease due to single or multiple genotypes.