Transcriptomics: Revisiting the genomes of herpesviruses
Combining integrative genomics and systems biology approaches has revealed new and conserved features in the genome of human herpesvirus 6.
- Julius-Maximilians-Universität Würzburg, Germany
Herpesviruses cause a range of human diseases but many factors complicate the efforts made to precisely map the size and origin of RNA transcripts coded by these pathogens. For example, some mRNAs can code for more than one protein, coding sequences may overlap with each other, and the genes that are expressed may change depending on cell types or stages in the viral cycle. Moreover, the level of expression can greatly vary from gene to gene, which makes it difficult to distinguish between rare viral transcripts and other genetic products that accumulate in infected cells and during viral replication. In fact, in most herpesviruses, the majority of the genome is transcribed to some degree, yet only the most highly expressed or genomically isolated units are readily detectable.
Several new techniques have allowed researchers to bypass these problems to better annotate the genomes of herpesviruses. A tailored RNA sequencing method called cRNA-Seq, which enriches for the 5’ ends of RNA transcripts, has allowed the mapping of transcription start sites; in parallel, ribosome profiling (Ribo-Seq) has helped to highlight translational start sites. Combined, these approaches have revealed dozens to hundreds of new genes in herpesviruses such as the human cytomegalovirus (Stern-Ginossar et al., 2012) and the Kaposi’s sarcoma-associated herpesvirus (Arias et al., 2014). When paired with long-read sequencing platforms (which provide additional information about the 3’ ends of transcripts), the new methods have also led to a better understanding of a number of pathogens in the herpes family. Now, in eLife, Noam Stern-Ginossar and colleagues at the Weizmann Institute of Science and the Hebrew University Hadassah Medical School – including Yaara Finkel as first author – report new insights into human herpesvirus 6A and 6B (Finkel et al., 2020).
The results help to correct and complement previous textbook genome annotations for herpesviruses. Due to the technical limitations of the time, the exact beginnings of many transcripts and coding sequences were assigned a priori, and inclusion into published gene lists relied on rather conservative criteria. For instance, a sequence was classified as an open reading frame (the part of a genetic sequence that can potentially be translated) if it had more than 100 amino acids and started with an AUG codon. Instead, Finkel et al. demonstrate that roughly one-third of open reading frames in human herpesvirus 6A and 6B contain alternative start codons, which are also used by eukaryotes and other herpesviruses (Kearse and Wilusz, 2017; Arias et al., 2014). For instance, strains of human cytomegalovirus can have different start codons for a given gene, which may influence biological properties (Brondke et al., 2007); such questions can now be investigated in herpesvirus 6A and 6B .
Another exciting finding is the identification of hundreds of short, internal or upstream open reading frames (Figure 1). The proteins encoded by many of these sequences are likely to be too small to have direct functions. However, some of these short open reading frames are close to (or overlap with) longer coding sequences, suggesting that they may regulate translation – particularly during the later stage of viral gene expression, when homeostasis in the host cells is most disrupted. Finkel et al. observed that several of these open reading frames are also transcribed in human cytomegalovirus, indicating important conserved roles across the family of viruses that herpesvirus 6A and 6B belong to.
Figure 1
Taking a closer look at the genomes of human herpesviruses 6.
Finkel et al. have used a combination of techniques to reannotate the genomes of human herpesviruses 6A and 6B. They have identified new open reading frames (268 in human herpesvirus 6A and 216 in human herpesvirus 6B) and corrected the annotation of existing frames (10 in human herpesvirus 6A and 11 in human herpesvirus 6B). The figure shows how an open reading frame called U30, which codes for an important protein in both human herpesvirus 6A and 6B, was reannotated. Data from Ribo-Seq (orange) revealed that the start of the open reading frame was downstream of what was expected based on the previous annotation (black) or cRNA-Seq information (blue), leading to a new, more accurate annotation for this sequence (green).
Combining several methods that can pinpoint both translational and transcriptional start sites – as Finkel et al. did – is particularly important because modern sequencing protocols are sensitive enough to identify rare transcription events, but they cannot distinguish between ‘real’ transcriptional units and biological artifacts. Whole-genome conclusions based on one technique or method of analysis are heavily influenced by experimental noise, technical limitations and even the specific algorithm used to interpret the data. For instance, estimates of the exact number of transcriptional start sites in human cytomegalovirus vary by thousands between studies that use different methods (Stern-Ginossar et al., 2012; Parida et al., 2019); in herpes simplex virus, these numbers can vary by over six-fold (Tombácz et al., 2019; Depledge et al., 2019).
While our appreciation of the coding capacity of pathogens increases, efforts must be made to integrate newly identified gene products into already established nomenclatures. The first waves of new annotations using high-throughput techniques will probably be revised as sequencing technology and analysis techniques improve, and the results are validated in the lab. In particular, new algorithms that can better distinguish signal-to-noise values could help to identify hundreds of additional peptides in a second revision of the human cytomegalovirus genome (Erhard et al., 2018). As our ability to sequence deeper develops, multifaceted studies such as the one by Finkel et al. will provide an excellent framework to help distinguish between rare functional events and technical noise when re-examining herpesvirus genome annotations.
References
-
Non-AUG translation: a new start for protein synthesis in eukaryotesGenes & Development 31:1717–1731.https://doi.org/10.1101/gad.305250.117
Article and author information
Author details
Publication history
Copyright
© 2020, Prusty and Whisnant
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 934
- views
-
- 114
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Transcriptomics: Revisiting the genomes of herpesviruses
eLife 9:e54037.
https://doi.org/10.7554/eLife.54037
Further reading
-
- Microbiology and Infectious Disease
Understanding the interplay between the HIV reservoir and the host immune system may yield insights into HIV persistence during antiretroviral therapy (ART) and inform strategies for a cure. Here, we applied machine learning (ML) approaches to cross-sectional high-parameter HIV reservoir and immunology data in order to characterize host–reservoir associations and generate new hypotheses about HIV reservoir biology. High-dimensional immunophenotyping, quantification of HIV-specific T cell responses, and measurement of genetically intact and total HIV proviral DNA frequencies were performed on peripheral blood samples from 115 people with HIV (PWH) on long-term ART. Analysis demonstrated that both intact and total proviral DNA frequencies were positively correlated with T cell activation and exhaustion. Years of ART and select bifunctional HIV-specific CD4 T cell responses were negatively correlated with the percentage of intact proviruses. A leave-one-covariate-out inference approach identified specific HIV reservoir and clinical–demographic parameters, such as age and biological sex, that were particularly important in predicting immunophenotypes. Overall, immune parameters were more strongly associated with total HIV proviral frequencies than intact proviral frequencies. Uniquely, however, expression of the IL-7 receptor alpha chain (CD127) on CD4 T cells was more strongly correlated with the intact reservoir. Unsupervised dimension reduction analysis identified two main clusters of PWH with distinct immune and reservoir characteristics. Using reservoir correlates identified in these initial analyses, decision tree methods were employed to visualize relationships among multiple immune and clinical–demographic parameters and the HIV reservoir. Finally, using random splits of our data as training-test sets, ML algorithms predicted with approximately 70% accuracy whether a given participant had qualitatively high or low levels of total or intact HIV DNA . The techniques described here may be useful for assessing global patterns within the increasingly high-dimensional data used in HIV reservoir and other studies of complex biology.
-
- Microbiology and Infectious Disease
Antibiotic tolerance in Mycobacterium tuberculosis reduces bacterial killing, worsens treatment outcomes, and contributes to resistance. We studied rifampicin tolerance in isolates with or without isoniazid resistance (IR). Using a minimum duration of killing assay, we measured rifampicin survival in isoniazid-susceptible (IS, n=119) and resistant (IR, n=84) isolates, correlating tolerance with bacterial growth, rifampicin minimum inhibitory concentrations (MICs), and isoniazid-resistant mutations. Longitudinal IR isolates were analyzed for changes in rifampicin tolerance and genetic variant emergence. The median time for rifampicin to reduce the bacterial population by 90% (MDK90) increased from 1.23 days (IS) and 1.31 days (IR) to 2.55 days (IS) and 1.98 days (IR) over 15–60 days of incubation, indicating fast and slow-growing tolerant sub-populations. A 6 log10-fold survival fraction classified tolerance as low, medium, or high, showing that IR is linked to increased tolerance and faster growth (OR = 2.68 for low vs. medium, OR = 4.42 for low vs. high, p-trend = 0.0003). High tolerance in IR isolates was associated with rifampicin treatment in patients and genetic microvariants. These findings suggest that IR tuberculosis should be assessed for high rifampicin tolerance to optimize treatment and prevent the development of multi-drug-resistant tuberculosis.
