Comparative genomic, transcriptomic, and proteomic reannotation of human herpesvirus 6

Liver-Humanized NSG-PiZ Mice Support the Study of Chronic Hepatitis B Virus Infection and Antiviral Therapies

Hepatitis B virus (HBV) is a pathogen of major public health importance that is largely incurable once a chronic infection is established. Only humans and great apes are fully permissive to HBV infection, and this species restriction has impacted HBV research by limiting the utility of small animal models. To combat HBV species restrictions and enable more in vivo studies, liver-humanized mouse models have been developed that are permissive to HBV infection and replication. Unfortunately, these models can be difficult to establish and are expensive commercially, which has limited their academic use. As an alternative mouse model to study HBV, we evaluated liver-humanized NSG-PiZ mice and showed that they are fully permissive to HBV. HBV selectively replicates in human hepatocytes within chimeric livers, and HBV-positive (HBV+) mice secrete infectious virions and hepatitis B surface antigen (HBsAg) into blood while also harboring covalently closed circular DNA (cccDNA). HBV+ mice develop chronic infections lasting at least 169 days, which should enable the study of new curative therapies targeting chronic HBV, and respond to entecavir therapy. Furthermore, HBV+ human hepatocytes in NSG-PiZ mice can be transduced by AAV3b and AAV.LK03 vectors, which should enable the study of gene therapies that target HBV. In summary, our data demonstrate that liver-humanized NSG-PiZ mice can be used as a robust and cost-effective alternative to existing chronic hepatitis B (CHB) models and may enable more academic research labs to study HBV disease pathogenesis and antiviral therapy. IMPORTANCE Liver-humanized mouse models have become the gold standard for the in vivo study of hepatitis B virus (HBV), yet their complexity and cost have prohibited widespread use of existing models in research. Here, we show that the NSG-PiZ liver-humanized mouse model, which is relatively inexpensive and simple to establish, can support chronic HBV infection. Infected mice are fully permissive to hepatitis B, supporting both active replication and spread, and can be used to study novel antiviral therapies. This model is a viable and cost-effective alternative to other liver-humanized mouse models that are used to study HBV.

Designing a multi-epitopes subunit vaccine against human herpes virus 6A based on molecular dynamics and immune stimulation

Human Herpesvirus 6A (HHV-6A) is a prevalent virus associated with various clinical manifestations, including neurological disorders, autoimmune diseases, and promotes tumor cell growth. HHV-6A is an enveloped, double-stranded DNA virus with a genome of approximately 160-170 kb containing a hundred open-reading frames. An immunoinformatics approach was applied to predict high immunogenic and non-allergenic CTL, HTL, and B cell epitopes and design a multi-epitope subunit vaccine based on HHV-6A glycoprotein B (gB), glycoprotein H (gH), and glycoprotein Q (gQ). The stability and correct folding of the modeled vaccines were confirmed through molecular dynamics simulation. Molecular docking found that the designed vaccines have a strong binding network with human TLR3, with K_d values of 1.5E-11 mol/L, 2.6E-12 mol/L, 6.5E-13 mol/L, and 7.1E-11 mol/L for gB-TLR3, gH-TLR3, gQ-TLR3, and the combined vaccine-TLR3, respectively. The codon adaptation index values of the vaccines were above 0.8, and their GC content was around 67 % (normal range 30-70 %), indicating their potential for high expression. Immune simulation analysis demonstrated robust immune responses against the vaccine, with approximately 650,000/ml combined IgG and IgM antibody titer. This study lays a strong foundation for developing a safe and effective vaccine against HHV-6A, with significant implications for treating associated conditions.

The Telomeric Repeats of HHV-6A Do Not Determine the Chromosome into Which the Virus Is Integrated

Human herpes virus 6A (HHV-6A) is able to integrate into the telomeric and subtelomeric regions of human chromosomes representing chromosomally integrated HHV-6A (ciHHV-6A). The integration starts from the right direct repeat (DRR) region. It has been shown experimentally that perfect telomeric repeats (pTMR) in the DRR region are required for the integration, while the absence of the imperfect telomeric repeats (impTMR) only slightly reduces the frequency of HHV-6 integration cases. The aim of this study was to determine whether telomeric repeats within DRR may define the chromosome into which the HHV-6A integrates. We analysed 66 HHV-6A genomes obtained from public databases. Insertion and deletion patterns of DRR regions were examined. We also compared TMR within the herpes virus DRR and human chromosome sequences retrieved from the Telomere-to-Telomere consortium. Our results show that telomeric repeats in DRR in circulating and ciHHV-6A have an affinity for all human chromosomes studied and thus do not define a chromosome for integration.

Comprehensive Investigations Relationship Between Viral Infections and Multiple Sclerosis Pathogenesis

Multiple sclerosis (MS) is a chronic autoimmune disease that affects the central nervous system (CNS). Compared to other types of self-limiting myelin disorders, MS compartmentalizes and maintains chronic inflammation in the CNS. Even though the exact cause of MS is unclear, it is assumed that genetic and environmental factors play an important role in susceptibility to this disease. The progression of MS is triggered by certain environmental factors, such as viral infections. The most important viruses that affect MS are Epstein-Barr virus (EBV), human herpes virus 6 (HHV-6), human endogenous retrovirus (HERV), cytomegalovirus (CMV), and varicella zoster virus (VZV). These viruses all have latent stages that allow them to escape immune detection and reactivate after exposure to various stimuli. Furthermore, their tropism for CNS and immune system cells explains their possible deleterious function in neuroinflammation. In this study, the effect of viral infections on MS disease focuses on the details of viruses that can change the risk of the disease. Paying attention to the most recent articles on the role of SARS-CoV-2 in MS disease, laboratory indicators show the interaction of the immune system with the virus. Also, strategies to prevent viruses that play a role in triggering MS are discussed, such as EBV, which is one of the most important.

Transmission dynamics of human herpesvirus 6A, 6B and 7 from whole genome sequences of families

While hundreds of thousands of human whole genome sequences (WGS) have been collected in the effort to better understand genetic determinants of disease, these whole genome sequences have less frequently been used to study another major determinant of human health: the human virome. Using the unmapped reads from WGS of over 1000 families, we present insights into the human blood DNA virome, focusing particularly on human herpesvirus (HHV) 6A, 6B, and 7. In addition to extensively cataloguing the viruses detected in WGS of human whole blood and lymphoblastoid cell lines, we use the family structure of our dataset to show that household drives transmission of several viruses, and identify the Mendelian inheritance patterns characteristic of inherited chromsomally integrated human herpesvirus 6 (iciHHV-6). Consistent with prior studies, we find that 0.6% of our dataset's population has iciHHV, and we locate candidate integration sequences for these cases. We document genetic diversity within exogenous and integrated HHV species and within integration sites of HHV-6. Finally, in the first observation of its kind, we present evidence that suggests widespread de novo HHV-6B integration and HHV-7 integration and reactivation in lymphoblastoid cell lines. These findings show that the unmapped read space of WGS is a promising source of data for virology research.

Assessing the potential correlation of polymorphisms in the IL6R with relative IL6 elevation in severely ill COVID-19 patients'

Background

Elevated Interleukin-6 (IL-6) may play an important role in the pathophysiology of COVID-19 yet attenuated response is not seen across all severe patients. We aimed to determine the effect of IL-6 baseline levels and other clinical variables on mortality and outcomes in hospitalized COVID-19 patients as well as to explore genetic variants associated with attenuated IL-6 response.

Methods

Baseline IL-6 cytokine levels were measured in hospitalized patients participating in ongoing ODYSSEY phase 3 randomized study of tradipitant and placebo in hospitalized patients with severe COVID-19 who are receiving supplemental oxygen support. Furthermore blood samples for whole genome sequencing analysis were collected from 150 participants.

Results

We report significantly elevated IL-6 in COVID-19 infected hospitalized patients, n = 100 (p-value < 0.0001) when compared to controls n = 324. We also report a significantly increased level of IL-6 (p-value < 0.01) between the severe and mild COVID-19 patients with severity defined on a WHO scale. Excessive IL-6 plasma levels correlate with higher mortality (p-value 0.001). Additionally, based on our classification analysis, combination of IL-6 elevation and high levels of serum glucose can identify highest risk-group of COVID19 patients.

Furthermore, we explore the role of genetic regulatory variants affecting baseline IL-6 levels specifically in COVID-19 patients. We have directly tested the association between variants in the IL6 and IL6R gene region and IL6 plasma levels. We provide results for a common IL-6 variant previously associated with pneumonia, rs1800795, and rs2228145 that was previously shown to affect IL-6 plasma levels, as well as report on novel variants associated with IL-6 plasma levels detected in our study patients.

Conclusions

While it is unlikely that “cytokine storm” is the norm in severe COVID19, baseline elevations above 150 pg/ml may be associated with worst outcomes and as such may warrant treatment considerations. So far no clinical studies used IL-6 baseline assessment to stratify the patient population participating in clinical studies. We believe that careful examination and interpretation of the IL-6 levels and genetic variants can help to determine a patient population with a potentially very robust clinical response to IL-6 inhibition.

Trial Registration: Clinicaltrials.gov: NCT04326426.

Evolutionary History of Endogenous Human Herpesvirus 6 Reflects Human Migration out of Africa

Human herpesvirus 6A and 6B (HHV-6) can integrate into the germline, and as a result, ∼70 million people harbor the genome of one of these viruses in every cell of their body. Until now, it has been largely unknown if 1) these integrations are ancient, 2) if they still occur, and 3) whether circulating virus strains differ from integrated ones. Here, we used next-generation sequencing and mining of public human genome data sets to generate the largest and most diverse collection of circulating and integrated HHV-6 genomes studied to date. In genomes of geographically dispersed, only distantly related people, we identified clades of integrated viruses that originated from a single ancestral event, confirming this with fluorescent in situ hybridization to directly observe the integration locus. In contrast to HHV-6B, circulating and integrated HHV-6A sequences form distinct clades, arguing against ongoing integration of circulating HHV-6A or “reactivation” of integrated HHV-6A. Taken together, our study provides the first comprehensive picture of the evolution of HHV-6, and reveals that integration of heritable HHV-6 has occurred since the time of, if not before, human migrations out of Africa.

The Role of Viral Proteins in the Regulation of Exosomes Biogenesis

Exosomes are membrane-bound vesicles of endocytic origin, secreted into the extracellular milieu, in which various biological components such as proteins, nucleic acids, and lipids reside. A variety of external stimuli can regulate the formation and secretion of exosomes, including viruses. Viruses have evolved clever strategies to establish effective infections by employing exosomes to cloak their viral genomes and gain entry into uninfected cells. While most recent exosomal studies have focused on clarifying the effect of these bioactive vesicles on viral infection, the mechanisms by which the virus regulates exosomes are still unclear and deserve further attention. This article is devoted to studying how viral components regulate exosomes biogenesis, composition, and secretion.

Mapping the Human Herpesvirus 6B transcriptome

The "omics" revolution of recent years has simplified the study of RNA transcripts produced during viral infection and under specific defined conditions. In the quest to find new and differentially expressed transcripts during the course of human Herpesvirus 6B (HHV-6B) infection, we made use of large-scale RNA sequencing to analyze the HHV-6B transcriptome during productive infection of human Molt-3 T-cells. Analyses were performed at different time points following infection and specific inhibitors were used to classify the kinetic class of each open reading frame (ORF) reported in the annotated genome of HHV-6B Z29 strain. The initial search focussed on HHV-6B-specific reads matching new HHV-6B transcripts. Differential expression of new HHV-6B transcripts were observed in all samples analyzed. The presence of many of these new HHV-6B transcripts were confirmed by RT-PCR and Sanger sequencing. Many of these transcripts represented new splice variants of previously reported ORFs, including some transcripts that have yet to be defined. Overall, our work demonstrates the diversity and the complexity of the HHV-6B transcriptome.IMPORTANCERNA sequencing (RNA-seq) is an important tool for studying RNA transcripts, particularly during active viral infection. We made use of RNA-seq to study human Herpesvirus 6B (HHV-6B) infection. Using six different time points, we were able to identify the presence of differentially spliced genes at 6, 9, 12, 24, 48 and 72 hours post-infection. Determination of the RNA profiles in the presence of cycloheximide (CHX) or phosphonoacetic acid (PAA) also permitted identification of the kinetic class of each ORF described in the annotated GenBank file. We also identified new spliced transcripts for certain genes and evaluated their relative expression over time. These data and next-generation sequencing (NGS) of the viral DNA have led us to propose a new version of the HHV-6B Z29 GenBank annotated file, without changing ORF names in order to facilitate trace back and correlate our work with previous studies on HHV-6B.

The role of micropeptides in biology

Micropeptides are small polypeptides coded by small open-reading frames. Progress in computational biology and the analyses of large-scale transcriptomes and proteomes have revealed that mammalian genomes produce a large number of transcripts encoding micropeptides. Many of these have been previously annotated as long noncoding RNAs. The role of micropeptides in cellular homeostasis maintenance has been demonstrated. This review discusses different types of micropeptides as well as methods to identify them, such as computational approaches, ribosome profiling, and mass spectrometry.

The tetraspanin protein CD9 modulates infection with human herpesvirus 6A and 6B in a CD46-dependent manner

Tetraspanins are four-span transmembrane proteins that organize the membrane by forming tetraspanin-enriched microdomains. These have been shown to be important for virus entry. The human herpesvirus (HHV)-6A receptor CD46 is known to form complexes with the tetraspanin CD9 and β1-integrins, however the significance of this for HHV-6A infection remains unexplored. Using a genetic approach, we demonstrate that knock out of CD46 abolishes binding to and infection of SupT1 cells by both HHV-6A and HHV-6B, establishing CD46 as a necessary receptor for productive infection of these cells. Knock out of CD9 in SupT1 cells had no effect on binding of either virus to the cell surface, but it reduced expression of immediate early transcripts to between 25-60% compared with the wild type cells. Although HHV-6B required CD46 for infection of SupT1, infection of Molt3 cells was independent of CD46 expression. Conversely, the absence of CD9 expression promoted infection of Molt3 cells with HHV-6B, indicating a negative role of CD9 for CD46-independent infection. Taken together, these data demonstrate that CD9 modulates infection with HHV-6A/B by promoting CD46-dependent infection and impairing CD46-independent infection. This also suggests that HHV-6A is strictly dependent on CD46 for entry, although other proteins, like CD9, may enhance the infection, whereas HHV-6B is more promiscuous and may use CD134, as demonstrated by others, CD46 in SupT1, and a novel yet unidentified receptor in Molt3 cells.Importance The mechanisms of entry of human herpesvirus (HHV)-6A and HHV-6B into host cells are of significance in order to develop novel drugs that may inhibit infection. To elucidate the contribution of the membrane proteins CD9 and CD46, we employed a genetic approach that eliminated these molecules from the host cell. This demonstrated that CD46 is critical for infection by HHV-6A, whereas infection by HHV-6B appeared to be more promiscuous. The infection of a T-cell line in the absence of CD46 and CD134 strongly suggest that an additional receptor for HHV-6B entry exists. Moreover, elimination of CD9 and subsequent reconstitution experiments demonstrated that CD9 promoted infection with HHV-6A and HHV-6B mediated by CD46, but inhibited infection with HHV-6B that occurred independent of CD46. Together, this demonstrated a CD46-dependent role of CD9 during infection with HHV-6A and HHV-6B and emphasized that HHV-6B may employ different entry mechanisms in various cells.

Variation in human herpesvirus 6B telomeric integration, excision, and transmission between tissues and individuals

Human herpesviruses 6A and 6B (HHV-6A/6B) are ubiquitous pathogens that persist lifelong in latent form and can cause severe conditions upon reactivation. They are spread by community-acquired infection of free virus (acqHHV6A/6B) and by germline transmission of inherited chromosomally integrated HHV-6A/6B (iciHHV-6A/6B) in telomeres. We exploited a hypervariable region of the HHV-6B genome to investigate the relationship between acquired and inherited virus and revealed predominantly maternal transmission of acqHHV-6B in families. Remarkably, we demonstrate that some copies of acqHHV-6B in saliva from healthy adults gained a telomere, indicative of integration and latency, and that the frequency of viral genome excision from telomeres in iciHHV-6B carriers is surprisingly high and varies between tissues. In addition, newly formed short telomeres generated by partial viral genome release are frequently lengthened, particularly in telomerase-expressing pluripotent cells. Consequently, iciHHV-6B carriers are mosaic for different iciHHV-6B structures, including circular extra-chromosomal forms that have the potential to reactivate. Finally, we show transmission of an HHV-6B strain from an iciHHV-6B mother to her non-iciHHV-6B son. Altogether, we demonstrate that iciHHV-6B can readily transition between telomere-integrated and free virus forms.

The U94 Gene of Human Herpesvirus 6: A Narrative Review of Its Role and Potential Functions

Human herpesvirus 6 (HHV-6) is a β-herpesvirus that is highly prevalent in the human population. HHV-6 comprises two recognized species (HHV-6A and HHV-6B). Despite different cell tropism and disease association, HHV-6A/B show high genome homology and harbor the conserved U94 gene, which is limited to HHV-6 and absent in all the other human herpesviruses. U94 has key functions in the virus life cycle and associated diseases, having demonstrated or putative roles in virus replication, integration, and reactivation. During natural infection, U94 elicits an immune response, and the prevalence and extent of the anti-U94 response are associated with specific diseases. Notably, U94 can entirely reproduce some virus effects at the cell level, including inhibition of cell migration, induction of cytokines and HLA-G expression, and angiogenesis inhibition, supporting a direct U94 role in the development of HHV-6-associated diseases. Moreover, specific U94 properties, such as the ability to modulate angiogenesis pathways, have been exploited to counteract cancer development. Here, we review the information available on this key HHV-6 gene, highlighting its potential uses.

Donor-to-recipient transmission and reactivation in a kidney transplant recipient of an inherited chromosomally integrated HHV-6A: Evidence and outcomes

Human herpesvirus (HHV)-6A can be inherited and chromosomally integrated (iciHHV-6A), and donor-to-recipient transmission has been reported in solid organ transplant. However, when HHV-6A reactivation happens after transplant, the source of HHV-6A is often not evident and its pathogenicity remains unclear. Here, we present an exhaustive case of donor-to-recipient transmission and reactivation of iciHHV-6A through kidney transplant. The absence of HHV-6A genome from the nails of the recipient excluded a recipient-related iciHHV-6A. Viral loads > 7 log₁₀ copies/10⁶ cells in donor blood samples and similarities of U38, U39, U69, and U100 viral genes between donor, recipient, and previously published iciHHV-6A strains are proof of donor-related transmission. Detection of noncoding HHV-6 snc-RNA14 using fluorescence in situ hybridization analysis and immunofluorescence staining of HHV-6A gp82/gp105 late proteins on kidney biopsies showed evidence of reactivation in the transplanted kidney. Because HHV-6A reactivation can be life threatening in immunocompromised patients, we provide several tools to help during the complete screening and diagnosis.

rs73185306 C/T Is Not a Predisposing Risk Factor for Inherited Chromosomally Integrated Human Herpesvirus 6A/B

Approximately 1% of people worldwide carry a copy of the human herpesvirus 6A or 6B (HHV-6A/B) in every cell of their body. This condition is referred to as inherited chromosomally integrated HHV-6A/B (iciHHV-6A/B). The mechanisms leading to iciHHV-6A/B chromosomal integration are yet to be identified. A recent report suggested that the rs73185306 C/T single-nucleotide polymorphism (SNP) represents a favorable predisposing factor leading to HHV-6A/B integration. After genotype analysis of an independent cohort (N = 11 967), we report no association between the rs73185306 C/T SNP and HHV-6A/B chromosomal integration (odds ratio, 0.90 [95% confidence interval, .54-1.51]; P = .69).

Unbiased optical mapping of telomere-integrated endogenous human herpesvirus 6

Next-generation sequencing technologies allowed sequencing of thousands of genomes. However, there are genomic regions that remain difficult to characterize, including telomeres, centromeres, and other low-complexity regions, as well as transposable elements and endogenous viruses. Human herpesvirus 6A and 6B (HHV-6A and HHV-6B) are closely related viruses that infect most humans and can integrate their genomes into the telomeres of infected cells. Integration also occurs in germ cells, meaning that the virus can be inherited and result in individuals harboring the virus in every cell of their body. The integrated virus can reactivate and cause disease in humans. While it is well established that the virus resides in the telomere region, the integration locus is poorly defined due to the low sequence complexity (TTAGGG)n of telomeres that cannot be easily resolved through sequencing. We therefore employed genome imaging of the integrated HHV-6A and HHV-6B genomes using whole-genome optical site mapping technology. Using this technology, we identified which chromosome arm harbors the virus genome and obtained a high-resolution map of the integration loci of multiple patients. Surprisingly, this revealed long telomere sequences at the virus-subtelomere junction that were previously missed using PCR-based approaches. Contrary to what was previously thought, our technique revealed that the telomere lengths of chromosomes harboring the integrated virus genome were comparable to the other chromosomes. Taken together, our data shed light on the genetic structure of the HHV-6A and HHV-6B integration locus, demonstrating the utility of optical mapping for the analysis of genomic regions that are difficult to sequence.

Cryptic transmission of SARS-CoV-2 in Washington state

After its emergence in Wuhan, China, in late November or early December 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus rapidly spread globally. Genome sequencing of SARS-CoV-2 allows the reconstruction of its transmission history, although this is contingent on sampling. We analyzed 453 SARS-CoV-2 genomes collected between 20 February and 15 March 2020 from infected patients in Washington state in the United States. We find that most SARS-CoV-2 infections sampled during this time derive from a single introduction in late January or early February 2020, which subsequently spread locally before active community surveillance was implemented.

The role of herpesvirus 6A and 6B in multiple sclerosis and epilepsy

Human herpesvirus 6A (HHV‐6A) and 6B (HHV‐6B) are two closely related viruses that can infect cells of the central nervous system (CNS). The similarities between these viruses have made it difficult to separate them on serological level. The broad term HHV‐6 remains when referring to studies where the two species were not distinguished, and as such, the seroprevalence is over 90% in the adult population. HHV‐6B has been detected in up to 100% of infants with the primary infection roseola infantum, but less is known about the primary infection of HHV‐6A. Both viruses are neurotropic and have capacity to establish lifelong latency in cells of the central nervous system, with potential to reactivate and cause complications later in life. HHV‐6A infection has been associated with an increased risk of multiple sclerosis (MS), whereas HHV‐6B is indicated to be involved in pathogenesis of epilepsy. These two associations show how neurological diseases might be caused by viral infections, but as suggested here, through completely different molecular mechanisms, in an autoimmune disease, such as MS, by triggering an overreaction of the immune system and in epilepsy by hampering internal cellular functions when the immune system fails to eliminate the virus. Understanding the viral mechanisms of primary infection and reactivation and their spectrum of associated symptoms will aid our ability to diagnose, treat and prevent these severe and chronic diseases. This review explores the role of HHV‐6A and HHV‐6B specifically in MS and epilepsy, the evidence to date and the future directions of this field.

Inherited Chromosomally Integrated Human Herpesvirus 6 Is a Risk Factor for Spontaneous Abortion

BACKGROUND

Human herpesvirus 6 (HHV-6) can be genetically transmitted from parent to child as inherited chromosomally integrated HHV-6 (iciHHV-6). HHV-6 reactivation occurs in pregnant women with iciHHV-6. We found no sex differences in the frequency of index cases with iciHHV-6 but inheritance from the father was more common. We evaluated the association between iciHHV-6 status and spontaneous abortion.

METHODS

iciHHV-6 was confirmed by high viral DNA copy numbers in whole blood and somatic cells. The origin of integrated viral genome, paternal or maternal, was examined using the same method. The pregnancy history of 23 mothers in families with iciHHV-6 and 285 mothers in families without iciHHV-6 was abstracted.

RESULTS

Of 23 iciHHV-6 index cases, 8 mothers and 15 fathers had iciHHV-6. Spontaneous abortion rates in mothers with and mothers without/fathers with iciHHV-6 and mothers in families without iciHHV-6 were 27.6%, 10.3%, and 14.8%, respectively (P = .012). Mothers with iciHHV-6 (odds ratio [OR], 6.41; 95% confidence interval [CI], 1.10-37.4) and maternal age at the most recent pregnancy ≥40 years (OR, 3.91; 95% CI, 1.30-11.8) were associated with 2 or more spontaneous abortions.

CONCLUSIONS

Mothers with iciHHV-6 is a risk factor for spontaneous abortion.

Using Direct RNA Nanopore Sequencing to Deconvolute Viral Transcriptomes

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 N6 -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.

Evolutionary analysis of exogenous and integrated HHV-6A/HHV-6B populations

Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6A/HHV-6B, iciHHV-6A/iciHHV-6B). We reconstructed the population structures of HHV-6A and HHV-6B, showing that HHV-6A diverged less than HHV-6B genomes from the projected common ancestral population. Thus, HHV-6B genomes experienced stronger drift, as also supported by calculation of nucleotide diversity and Tajima's D. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, ancestry proportions were similar in exogenous HHV-6B viruses and iciHHV-6B. In line with previous indications, this suggests the distinct exogenous viral populations that originated iciHHV-6B in subjects with European and Asian ancestry are still causing infections in the corresponding geographic areas. Notably, for both iciHHV-6A and iciHHV-6B, we found that European and American sequences tend to have high proportions of ancestry from viral populations that experienced considerable drift, suggesting that they underwent one or more bottlenecks followed by population expansion. Finally, analysis of HHV-6B exogenous viruses sampled in Japan indicated that proportions of ancestry components of most of these viruses are different from the majority of those sampled in the USA. More generally, we show that, in both viral species, both integrated and exogenous viral genomes have different ancestry components, partially depending on geographic location. It would be extremely important to determine whether such differences account for the diversity of HHV-6A/HHV-6B-associated clinical symptoms and epidemiology. Also, the sequencing of additional exogenous and integrated viral genomes will be instrumental to confirm and expand our conclusions, which are based on a relatively small number of genomes, sequenced with variable quality, and with unequal sampling in terms of geographic origin.

Pathogen or Bystander: Clinical Significance of Detecting Human Herpesvirus 6 in Pediatric Cerebrospinal Fluid

Human herpesvirus 6 (HHV-6) is an important cause of meningitis and meningoencephalitis. As testing for HHV-6 in cerebrospinal fluid (CSF) is more readily available using the FilmArray Meningitis/Encephalitis panel (FA-ME; BioFire Diagnostics, Salt Lake City, UT), we aimed to determine the clinical significance of detecting HHV-6 in order to identify true infections and to ensure appropriate antiviral initiation. Chart review on 25 patients positive for HHV-6 by FA-ME was performed to determine clinical presentation, comorbidity, treatment, and outcome. The presence of chromosomally integrated HHV-6 (ciHHV-6) DNA was also investigated. Of 1,005 children tested by FA-ME, HHV-6 was detected in 25 (2.5%). Five patients were diagnosed with either HHV-6 meningitis or meningoencephalitis based on HHV-6 detection in CSF, clinical presentation, and radiographic findings. Detection of HHV-6 by FA-ME led to discontinuation of acyclovir within 12.0 h in all 12 patients empirically treated with acyclovir. Six of the 12 patients were started on ganciclovir therapy within 6.8 h; 4 of these were treated specifically for HHV-6 infection, whereas therapy was discontinued in the remaining 2 patients. CSF parameters were not generally predictive of HHV-6 positivity. The presence of ciHHV-6 was confirmed in 3 of 18 patients who could be tested. Five of the 25 patients included in the study were diagnosed with HHV-6 meningitis/meningoencephalitis. FA-ME results led to discontinuation of empirical antiviral treatment in 12 patients and appropriate initiation of ganciclovir in 4 patients. In our institution, detection of HHV-6 using FA-ME led to faster establishment of disease etiology and optimization of antimicrobial therapy.

Cryptic transmission of SARS-CoV-2 in Washington State

Following its emergence in Wuhan, China, in late November or early December 2019, the SARS-CoV-2 virus has rapidly spread throughout the world. On March 11, 2020, the World Health Organization declared Coronavirus Disease 2019 (COVID-19) a pandemic. Genome sequencing of SARS-CoV-2 strains allows for the reconstruction of transmission history connecting these infections. Here, we analyze 346 SARS-CoV-2 genomes from samples collected between 20 February and 15 March 2020 from infected patients in Washington State, USA. We found that the large majority of SARS-CoV-2 infections sampled during this time frame appeared to have derived from a single introduction event into the state in late January or early February 2020 and subsequent local spread, strongly suggesting cryptic spread of COVID-19 during the months of January and February 2020, before active community surveillance was implemented. We estimate a common ancestor of this outbreak clade as occurring between 18 January and 9 February 2020. From genomic data, we estimate an exponential doubling between 2.4 and 5.1 days. These results highlight the need for large-scale community surveillance for SARS-CoV-2 introductions and spread and the power of pathogen genomics to inform epidemiological understanding.

Comprehensive annotations of human herpesvirus 6A and 6B genomes reveal novel and conserved genomic features

Human herpesvirus-6 (HHV-6) A and B are ubiquitous betaherpesviruses, infecting the majority of the human population. They encompass large genomes and our understanding of their protein coding potential is far from complete. Here, we employ ribosome-profiling and systematic transcript-analysis to experimentally define HHV-6 translation products. We identify hundreds of new open reading frames (ORFs), including upstream ORFs (uORFs) and internal ORFs (iORFs), generating a complete unbiased atlas of HHV-6 proteome. By integrating systematic data from the prototypic betaherpesvirus, human cytomegalovirus, we uncover numerous uORFs and iORFs conserved across betaherpesviruses and we show uORFs are enriched in late viral genes. We identified three highly abundant HHV-6 encoded long non-coding RNAs, one of which generates a non-polyadenylated stable intron appearing to be a conserved feature of betaherpesviruses. Overall, our work reveals the complexity of HHV-6 genomes and highlights novel features conserved between betaherpesviruses, providing a rich resource for future functional studies.

Betaherpesvirus assembly and egress: Recent advances illuminate the path

The human betaherpesviruses, human cytomegalovirus (HCMV; species Human betaherpesvirus 5) and human herpesviruses 6A, 6B, and 7 (HHV-6A, -6B, and -7; species Human betaherpesviruses 6A, 6B, and 7) are highly prevalent and can cause severe disease in immune-compromised and immune-naive populations in well- and under-developed communities. Herpesvirus virion assembly is an intricate process that requires viral orchestration of host systems. In this review, we describe recent advances in some of the many cellular events relevant to assembly and egress of betaherpesvirus virions. These include modifications of host metabolic, immune, and autophagic/recycling systems. In addition, we discuss unique aspects of betaherpesvirus virion structure, virion assembly, and the cellular pathways employed during virion egress.

Inherited Chromosomally Integrated Human Herpesvirus 6 Demonstrates Tissue-Specific RNA Expression In Vivo That Correlates with an Increased Antibody Immune Response

HHV-6A and -6B are human herpesviruses that have the unique property of being able to integrate into the telomeric regions of human chromosomes. Approximately 1% of the world’s population carries integrated HHV-6A/B genome in every cell of their body. Whether viral genes are transcriptionally active in these individuals is unclear. By taking advantage of a unique tissue-specific gene expression data set, we showed that the majority of tissues from iciHHV-6 individuals do not show HHV-6 gene expression. Brain and testes showed the highest tissue-specific expression of HHV-6 genes in two separate data sets. Two HHV-6 genes, U90 (immediate early 1 protein) and U100 (glycoproteins Q1 and Q2), were found to be selectively and consistently expressed across several human tissues. Expression of U90 translates into an increase in antigen-specific antibody response in iciHHV-6A/B⁺ subjects relative to controls. Future studies will be needed to determine the mechanism of gene expression, the effects of these genes on human gene transcription networks, and the pathophysiological impact of having increased viral protein expression in tissue in conjunction with increased antigen-specific antibody production.

Human herpesviruses 6A and 6B (HHV-6A and HHV-6B) are human viruses capable of chromosomal integration. Approximately 1% of the human population carries one copy of HHV-6A/B integrated into every cell in their body, referred to as inherited chromosomally integrated human herpesvirus 6A/B (iciHHV-6A/B). Whether iciHHV-6A/B is transcriptionally active in vivo and how it shapes the immunological response are still unclear. In this study, we screened DNA sequencing (DNA-seq) and transcriptome sequencing (RNA-seq) data for 650 individuals available through the Genotype-Tissue Expression (GTEx) project and identified 2 iciHHV-6A- and 4 iciHHV-6B-positive candidates. When corresponding tissue-specific gene expression signatures were analyzed, low levels HHV-6A/B gene expression was found across multiple tissues, with the highest levels of gene expression in the brain (specifically for HHV-6A), testis, esophagus, and adrenal gland. U90 and U100 were the most highly expressed HHV-6 genes in both iciHHV-6A- and iciHHV-6B-positive individuals. To assess whether tissue-specific gene expression from iciHHV-6A/B influences the immune response, a cohort of 15,498 subjects was screened and 85 iciHHV-6A/B⁺ subjects were identified. Plasma samples from iciHHV-6A/B⁺ and age- and sex-matched controls were analyzed for antibodies to control antigens (cytomegalovirus [CMV], Epstein-Barr virus [EBV], and influenza virus [FLU]) or HHV-6A/B antigens. Our results indicate that iciHHV-6A/B⁺ subjects have significantly more antibodies against the U90 gene product (IE1) than do non-iciHHV-6-positive individuals. Antibody responses against EBV and FLU antigens or HHV-6A/B gene products either not expressed or expressed at low levels, such as U47, U57, and U72, were identical between controls and iciHHV-6A/B⁺ subjects. CMV-seropositive individuals with iciHHV-6A/B⁺ have more antibodies against CMV pp150 than do CMV-seropositive controls. These results argue that spontaneous gene expression from integrated HHV-6A/B leads to an increase in antigenic burden that translates into a more robust HHV-6A/B-specific antibody response.

IMPORTANCE HHV-6A and -6B are human herpesviruses that have the unique property of being able to integrate into the telomeric regions of human chromosomes. Approximately 1% of the world’s population carries integrated HHV-6A/B genome in every cell of their body. Whether viral genes are transcriptionally active in these individuals is unclear. By taking advantage of a unique tissue-specific gene expression data set, we showed that the majority of tissues from iciHHV-6 individuals do not show HHV-6 gene expression. Brain and testes showed the highest tissue-specific expression of HHV-6 genes in two separate data sets. Two HHV-6 genes, U90 (immediate early 1 protein) and U100 (glycoproteins Q1 and Q2), were found to be selectively and consistently expressed across several human tissues. Expression of U90 translates into an increase in antigen-specific antibody response in iciHHV-6A/B⁺ subjects relative to controls. Future studies will be needed to determine the mechanism of gene expression, the effects of these genes on human gene transcription networks, and the pathophysiological impact of having increased viral protein expression in tissue in conjunction with increased antigen-specific antibody production.

Evaluation of liver failure in a pediatric transplant recipient of a liver allograft with inherited chromosomally integrated HHV-6B

BACKGROUND

Active infections of human herpesvirus 6B (HHV-6B) are frequent in immunocompromised recipients after transplantation. Nevertheless, they need to be distinguished from latent inherited chromosomally integrated genomes (iciHHV-6) present in about 1% of the population to avoid unnecessary administration of toxic antivirals.

METHODS

A 5-year-old child presented with acute liver allograft rejection associated with HHV-6 DNA in plasma, which led to an unfavorable outcome. We investigated the possibility of HHV-6 infection derived from an iciHHV-6 present in the donor's liver using molecular and histopathology studies in various tissues, including quantification of HHV-6 DNA, genotyping, sequencing for antiviral resistance genes, relative quantification of viral transcripts, and detection of gB and gH viral proteins.

RESULTS

The presence of iciHHV-6B was evidenced in the donor with signs of reactivation in the gallbladder and transplanted liver (detection of HHV-6B mRNA and late proteins). This localized expression could have played a role in liver rejection. Low viral loads in the recipient's plasma, with identical partial U39 sequences, were in favor of viral DNA released from the transplanted liver rather than a systemic infection.

CONCLUSIONS

Determination of iciHHV-6 status before transplantation should be considered to guide clinical decisions, such as antiviral prophylaxis, viral load monitoring, and antiviral therapy.

A complex evolutionary relationship between HHV-6A and HHV-6B

Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6, iciHHV6). We reconstructed the population structure of HHV-6 and we show that HHV-6A genomes diverged less than HHV-6B genomes from the ancestral common HHV-6A/B population. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, exogenous viral and iciHHV-6B populations were similar in terms of ancestry components, with no evident geographic structuring. Most HHV-6B genomes sampled to date derive from viral populations that experienced considerable drift. However, a population of HHV-6 exogenous viruses, currently classified as HHV-6B and sampled in New York state, formed a separate cluster (NY cluster) and harbored a considerable portion of HHV-6A-like ancestry. Recombination detection methods identified these viruses as interspecies recombinants, but phylogenetic reconstruction indicated that the recombination signals are due to shared ancestry. In analogy to iciHHV-6A, NY cluster viruses have high nucleotide diversity and constant population size. We propose that HHV-6A sequences and the NY cluster population diverged from an ancestral HHV-6A-like population. A relatively recent bottleneck of the NY (or a related) population with subsequent expansion originated most HHV-6B genomes currently sampled. Our findings indicate that the distinction between HHV-6A and -6B is not as clear-cut as previously thought. More generally, epidemiological and clinical surveys would benefit from taking HHV-6 genetic diversity into account.

Human Herpesvirus Sequencing in the Genomic Era: The Growing Ranks of the Herpetic Legion

The nine human herpesviruses are some of the most ubiquitous pathogens worldwide, causing life-long latent infection in a variety of different tissues. Human herpesviruses range from mild childhood infections to known tumour viruses and 'trolls of transplantation'. Epstein-Barr virus was the first human herpesvirus to have its whole genome sequenced; GenBank now includes thousands of herpesvirus genomes. This review will cover some of the recent advances in our understanding of herpesvirus diversity and disease that have come about as a result of new sequencing technologies, such as target enrichment and long-read sequencing. It will also look at the problem of resolving mixed-genotype infections, whether with short or long-read sequencing methods; and conclude with some thoughts on the future of the field as herpesvirus population genomics becomes a reality.

Genome-Wide Approach to the CD4 T-Cell Response to Human Herpesvirus 6B

Human herpesvirus 6 (HHV-6) and cytomegalovirus (CMV) are population-prevalent betaherpesviruses with intermittent lytic replication that can be pathogenic in immunocompromised hosts. Elucidation of the adaptive immune response is valuable for understanding pathogenesis and designing novel treatments. Knowledge of T-cell antigens has reached the genome-wide level for CMV and other human herpesviruses, but study of HHV-6 is at an earlier stage. Using rare-cell enrichment combined with an HLA-agnostic, proteome-wide approach, we queried HHV-6B-specific CD4 T cells from 18 healthy donors with each known HHV-6B protein. We detected a low abundance of HHV-6-specific CD4 T cells in blood; however, the within-person CD4 T-cell response is quite broad: the median number of open reading frame (ORF) products recognized was nine per person. Overall, the data expand the number of documented HHV-6B CD4 T-cell antigens from approximately 11 to 60. Epitopes in the proteins encoded by U14, U90, and U95 were mapped with synthetic peptides, and HLA restriction was defined for some responses. Intriguingly, CD4 T-cell antigens newly described in this report are among the most population prevalent, including U73, U72, U95, and U30. Our results indicate that selection of HHV-6B ORFs for immunotherapy should consider this expanded panel of HHV-6B antigens.IMPORTANCE Human herpesvirus 6 is highly prevalent and maintains chronic infection in immunocompetent individuals, with the potential to replicate widely in settings of immunosuppression, leading to clinical disease. Antiviral compounds may be ineffective and/or pose dose-limiting toxicity, and therefore, immune-based therapies have garnered increased interest in recent years. Attempts at addressing this unmet medical need begin with understanding the cellular response to HHV-6 at the individual and population levels. The present study provides a comprehensive assessment of HHV-6-specific T-cell responses that may inform the development of cell-based therapies directed at this virus.

Naturally processed HLA-DR3-restricted HHV-6B peptides are recognized broadly with polyfunctional and cytotoxic CD4 T-cell responses

Human herpes virus 6B (HHV-6B) is a widespread virus that infects most people early in infancy and establishes a chronic life-long infection with periodic reactivation. CD4 T cells have been implicated in control of HHV-6B, but antigenic targets and functional characteristics of the CD4 T-cell response are poorly understood. We identified 25 naturally processed MHC-II peptides, derived from six different HHV-6B proteins, and showed that they were recognized by CD4 T-cell responses in HLA-matched donors. The peptides were identified by mass spectrometry after elution from HLA-DR molecules isolated from HHV-6B-infected T cells. The peptides showed strong binding to matched HLA alleles and elicited recall T-cell responses in vitro. T-cell lines expanded in vitro were used for functional characterization of the response. Responding cells were mainly CD3⁺ CD4⁺ , produced IFN-γ, TNF-α, and low levels of IL-2, alone or in combination, highlighting the presence of polyfunctional T cells in the overall response. Many of the responding cells mobilized CD107a, stored granzyme B, and mediated specific killing of peptide-pulsed target cells. These results highlight a potential role for polyfunctional cytotoxic CD4 T cells in the long-term control of HHV-6B infection.

The Transcriptional Landscape of Marek's Disease Virus in Primary Chicken B Cells Reveals Novel Splice Variants and Genes

Marek's disease virus (MDV) is an oncogenic alphaherpesvirus that infects chickens and poses a serious threat to poultry health. In infected animals, MDV efficiently replicates in B cells in various lymphoid organs. Despite many years of research, the viral transcriptome in primary target cells of MDV remained unknown. In this study, we uncovered the transcriptional landscape of the very virulent RB1B strain and the attenuated CVI988/Rispens vaccine strain in primary chicken B cells using high-throughput RNA-sequencing. Our data confirmed the expression of known genes, but also identified a novel spliced MDV gene in the unique short region of the genome. Furthermore, de novo transcriptome assembly revealed extensive splicing of viral genes resulting in coding and non-coding RNA transcripts. A novel splicing isoform of MDV UL15 could also be confirmed by mass spectrometry and RT-PCR. In addition, we could demonstrate that the associated transcriptional motifs are highly conserved and closely resembled those of the host transcriptional machinery. Taken together, our data allow a comprehensive re-annotation of the MDV genome with novel genes and splice variants that could be targeted in further research on MDV replication and tumorigenesis.

RNA Sequencing of the In Vivo Human Herpesvirus 6B Transcriptome To Identify Targets for Clinical Assays Distinguishing between Latent and Active Infections

Human herpesvirus 6B (HHV-6B) DNA is frequently detected in human samples. Diagnostic assays distinguishing HHV-6B reactivation from latency are limited. This has impaired strategies to diagnose and treat HHV-6B-associated diseases. We used RNA sequencing to characterize and compare the HHV-6B transcriptome in multiple sample types, including (i) whole blood from hematopoietic cell transplant (HCT) recipients with and without HHV-6B plasma viremia, (ii) tumor tissue samples from subjects with large B cell lymphoma infected with HHV-6B, (iii) lymphoblastoid cell lines (LCLs) from subjects with inherited chromosomally integrated HHV-6B or latent infection with HHV-6B, and (iv) HHV-6B Z29 infected SupT1 CD4⁺ T cells. We demonstrated substantial overlap in the HHV-6B transcriptome observed in in vivo and in vitro samples, although there was variability in the breadth and quantity of gene expression across samples. The HHV-6B viral polymerase gene U38 was the only HHV-6B transcript detected in all next-generation RNA sequencing (RNA-seq) data sets and was one of the most highly expressed genes. We developed a novel reverse transcription-PCR assay targeting HHV-6B U38, which identified U38 mRNA in all tested whole-blood samples from patients with concurrent HHV-6B viremia. No HHV-6B U38 transcripts were detected by RNA-seq or reverse transcription-real-time quantitative PCR (RT-qPCR) in whole-blood samples from subjects without HHV-6B plasma detection or from latently infected LCLs. A RT-qPCR assay for HHV-6B U38 may be useful to identify lytic HHV-6B infection in nonplasma samples and samples from individuals with inherited chromosomally integrated HHV-6B. This study also demonstrates the feasibility of transcriptomic analyses for HCT recipients.IMPORTANCE Human herpesvirus 6B (HHV-6B) is a DNA virus that infects most children within the first few years of life. After primary infection, HHV-6B persists as a chronic, latent infection in many cell types. Additionally, HHV-6B can integrate into germ line chromosomes, resulting in individuals with viral DNA in every nucleated cell. Given that PCR to detect viral DNA is the mainstay for diagnosing HHV-6B infection, the characteristics of HHV-6B infection complicate efforts to distinguish between latent and active viral infection, particularly in immunocompromised patients who have frequent HHV-6B reactivation. In this study, we used RNA sequencing to characterize the HHV-6B gene expression profile in multiple sample types, and our findings identified evidence-based targets for diagnostic tests that distinguish between latent and active viral infection.

Ultrasensitive Capture of Human Herpes Simplex Virus Genomes Directly from Clinical Samples Reveals Extraordinarily Limited Evolution in Cell Culture

Herpes simplex viruses (HSVs) are difficult to sequence due to their large DNA genome, high GC content, and the presence of repeats. To date, most HSV genomes have been recovered from culture isolates, raising concern that these genomes may not accurately represent circulating clinical strains. We report the development and validation of a DNA oligonucleotide hybridization panel to recover nearly complete HSV genomes at abundances up to 50,000-fold lower than previously reported. Using copy number information on herpesvirus and host DNA background via quantitative PCR, we developed a protocol for pooling for cost-effective recovery of more than 50 HSV-1 or HSV-2 genomes per MiSeq run. We demonstrate the ability to recover >99% of the HSV genome at >100× coverage in 72 h at viral loads that allow whole-genome recovery from latently infected ganglia. We also report a new computational pipeline for rapid HSV genome assembly and annotation. Using the above tools and a series of 17 HSV-1-positive clinical swabs sent to our laboratory for viral isolation, we show limited evolution of HSV-1 during viral isolation in human fibroblast cells compared to the original clinical samples. Our data indicate that previous studies using low-passage-number clinical isolates of herpes simplex viruses are reflective of the viral sequences present in the lesion and thus can be used in phylogenetic analyses. We also detect superinfection within a single sample with unrelated HSV-1 strains recovered from separate oral lesions in an immunosuppressed patient during a 2.5-week period, illustrating the power of direct-from-specimen sequencing of HSV.IMPORTANCE Herpes simplex viruses affect more than 4 billion people across the globe, constituting a large burden of disease. Understanding the global diversity of herpes simplex viruses is important for diagnostics and therapeutics as well as cure research and tracking transmission among humans. To date, most HSV genomics has been performed on culture isolates and DNA swabs with high quantities of virus. We describe the development of wet-lab and computational tools that enable the accurate sequencing of near-complete genomes of HSV-1 and HSV-2 directly from clinical specimens at abundances >50,000-fold lower than previously sequenced and at significantly reduced cost. We use these tools to profile circulating HSV-1 strains in the community and illustrate limited changes to the viral genome during the viral isolation process. These techniques enable cost-effective, rapid sequencing of HSV-1 and HSV-2 genomes that will help enable improved detection, surveillance, and control of this human pathogen.

Copy Number Heterogeneity, Large Origin Tandem Repeats, and Interspecies Recombination in Human Herpesvirus 6A (HHV-6A) and HHV-6B Reference Strains

Quantitative PCR is a diagnostic pillar for clinical virology testing, and reference materials are necessary for accurate, comparable quantitation between clinical laboratories. Accurate quantitation of human herpesvirus 6A/B (HHV-6A/B) is important for detection of viral reactivation and inherited chromosomally integrated HHV-6A/B in immunocompromised patients. Reference materials in clinical virology commonly consist of laboratory-adapted viral strains that may be affected by the culture process. We performed next-generation sequencing to make relative copy number measurements at single nucleotide resolution of eight candidate HHV-6A and seven HHV-6B reference strains and DNA materials from the HHV-6 Foundation and Advanced Biotechnologies Inc. Eleven of 17 (65%) HHV-6A/B candidate reference materials showed multiple copies of the origin of replication upstream of the U41 gene by next-generation sequencing. These large tandem repeats arose independently in culture-adapted HHV-6A and HHV-6B strains, measuring 1,254 bp and 983 bp, respectively. The average copy number measured was between 5 and 10 times the number of copies of the rest of the genome. We also report the first interspecies recombinant HHV-6A/B strain with a HHV-6A backbone and a >5.5-kb region from HHV-6B, from U41 to U43, that covered the origin tandem repeat. Specific HHV-6A reference strains demonstrated duplication of regions at U1/U2, U87, and U89, as well as deletion in the U12-to-U24 region and the U94/U95 genes. HHV-6A/B strains derived from cord blood mononuclear cells from different laboratories on different continents with fewer passages revealed no copy number differences throughout the viral genome. These data indicate that large origin tandem duplications are an adaptation of both HHV-6A and HHV-6B in culture and show interspecies recombination is possible within the Betaherpesvirinae.IMPORTANCE Anything in science that needs to be quantitated requires a standard unit of measurement. This includes viruses, for which quantitation increasingly determines definitions of pathology and guidelines for treatment. However, the act of making standard or reference material in virology can alter its very accuracy through genomic duplications, insertions, and rearrangements. We used deep sequencing to examine candidate reference strains for HHV-6, a ubiquitous human virus that can reactivate in the immunocompromised population and is integrated into the human genome in every cell of the body for 1% of people worldwide. We found large tandem repeats in the origin of replication for both HHV-6A and HHV-6B that are selected for in culture. We also found the first interspecies recombinant between HHV-6A and HHV-6B, a phenomenon that is well known in alphaherpesviruses but to date has not been seen in betaherpesviruses. These data critically inform HHV-6A/B biology and the standard selection process.