Next Article in Journal
Correction: Zablocki, O.; et al. Diversity of dsDNA Viruses in a South African Hot Spring Assessed by Metagenomics and Microscopy. Viruses 2017, 9, 348
Next Article in Special Issue
The Inhibition of HIV-1 Entry Imposed by Interferon Inducible Transmembrane Proteins Is Independent of Co-Receptor Usage
Previous Article in Journal
A Model to Detect Autochthonous Group 1 and 2 Brazilian Vaccinia virus Coinfections: Development of a qPCR Tool for Diagnosis and Pathogenesis Studies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 10
Issue 1
10.3390/v10010016
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyses of Tissue Culture Adaptation of Human Herpesvirus-6A by Whole Genome Deep Sequencing Redefines the Reference Sequence and Identifies Virus Entry Complex Changes

by
Joshua G. Tweedy
1,
Eric Escriva
1,2,
Maya Topf
2 and
Ursula A. Gompels
1,*
1
Department of Pathogen Molecular Biology, London School of Hygiene & Tropical Medicine, University of London, London WC1E 7HT, UK
2
Institute for Structural and Molecular Biology, Department Biology, Birkbeck College University of London, London WC1E 7HX, UK
*
Author to whom correspondence should be addressed.
Viruses 2018, 10(1), 16; https://0-doi-org.brum.beds.ac.uk/10.3390/v10010016
Submission received: 30 November 2017 / Revised: 23 December 2017 / Accepted: 26 December 2017 / Published: 31 December 2017
(This article belongs to the Special Issue Breakthroughs in Viral Replication)
Browse Figures
Versions Notes

Abstract

:
Tissue-culture adaptation of viruses can modulate infection. Laboratory passage and bacterial artificial chromosome (BAC)mid cloning of human cytomegalovirus, HCMV, resulted in genomic deletions and rearrangements altering genes encoding the virus entry complex, which affected cellular tropism, virulence, and vaccine development. Here, we analyse these effects on the reference genome for related betaherpesviruses, Roseolovirus, human herpesvirus 6A (HHV-6A) strain U1102. This virus is also naturally “cloned” by germline subtelomeric chromosomal-integration in approximately 1% of human populations, and accurate references are key to understanding pathological relationships between exogenous and endogenous virus. Using whole genome next-generation deep-sequencing Illumina-based methods, we compared the original isolate to tissue-culture passaged and the BACmid-cloned virus. This re-defined the reference genome showing 32 corrections and 5 polymorphisms. Furthermore, minor variant analyses of passaged and BACmid virus identified emerging populations of a further 32 single nucleotide polymorphisms (SNPs) in 10 loci, half non-synonymous indicating cell-culture selection. Analyses of the BAC-virus genome showed deletion of the BAC cassette via loxP recombination removing green fluorescent protein (GFP)-based selection. As shown for HCMV culture effects, select HHV-6A SNPs mapped to genes encoding mediators of virus cellular entry, including virus envelope glycoprotein genes gB and the gH/gL complex. Comparative models suggest stabilisation of the post-fusion conformation. These SNPs are essential to consider in vaccine-design, antimicrobial-resistance, and pathogenesis.

1. Introduction

Establishing reference genomes are essential for characterisation of pathogens, their interaction with the host, development of vaccines, and responses to therapies including antimicrobial resistance. Betaherpesviruses are a classic example of effects of laboratory passage of viruses distorting understanding of natural transmission. In this case, betaherpesvirus human cytomegalovirus (HCMV) was originally cultured in human fibroblasts, the reference AD169 strain, with derivation of its genome being the largest at the time [1]. However, increasing evidence showed differences with strains grown in endothelial cells, the cell type relevant for characterising one of the main pathologies of the virus and placental infection causing congenital disease with cytomegalovirus (CMV) the main infectious cause of birth defects, including microcephaly even in the setting of emergent Zika virus. Derivation of the genomic sequence of early passage clinical strains showed that an almost 20kb deletion occurred in the AD169 genome [2]. Subsequent extensive analyses, using next generation sequencing technologies, have shown that clinical isolates, which are not cultured, do not have this deletion. Furthermore, genes located in the “deleted” region encoded essential elements dictating the cellular tropism of the virus. Alternate gH/gL complexes were established as main components of cellular tropism and fusogenic entry of the virus into cells. The trimer, gH/gL/gO, related to fibroblastic entry of cells while the competing pentamer form gH/gL/UL128/UL130/UL131 dictated endothelial, leukocyte, and epithelial cell entry [3]. The pentameric components were encoded by genes in the “deleted” region of the genomes of laboratory-passaged virus. bacterial artificial chromosome (BAC)mid “repair” of these genes rescued infection defects, but also could be mutagenic themselves, if there were single nucleotide polymorphisms (SNPs) affecting efficiency of tissue culture replication [4]. The resulting complexity of these components are essential to vaccine design and generation of a neutralising antibody, even to new vaccine vectors directed at other infections and modulating cellular immunity [5,6,7]. These could never have been identified without the derivation of the actual wild type genomes involved in natural transmission.
Similar analyses of polymorphisms in the reference human betaherpesvirus 6, HHV-6A, have not yet been conducted. However, isolated reports indicate homologous regions are targeted and may also occur in virus replicated in vivo. HHV-6A also encodes a gH/gL complex central to infectious entry, as well as alternate multimeric complexes affecting cellular tropism; these are a related trimeric complex gH/gL/gO and a quadrameric complex gH/gL/gQ1/gQ2 [8,9,10,11]. The gQ1 and gQ2 genes are in the homologous pathogenicity region to that in HCMV. Polymorphisms in gQ2 have similarly been detected in tissue culture passaged virus as well as in a BACmid clone of human herpesvirus 6A (HHV-6A) strain U1102 [12]. Furthermore, there is a positional homologue of HCMV chemokine like molecules UL128/130/131, which is HHV-6A chemokine U83A [13]. Polymorphisms in HHV-6 strain U1102 virus passaged in vitro as well as in vivo strains have been reported, which affect translated reading frames of this gene, alternating between an extended encoded mature secreted chemokine and a non-secreted version [13,14]. In addition, a few corrections have been reported, including a CG inversion in U83A [15], GG insertions in DR6 [16], and nucleotide insertion in U86 [17]. Two distinct but closely related genomes of HHV-6A strain GS have been reported (HHV-6A GS KJ123690.1 and KC465951.1), indicating potential for virus population polymorphisms; however, initial reports of this strain as a human B-lymphotropic virus (HBLV) describes isolation from six patients, therefore, there could be mixtures [18]. A systematic analysis has not been undertaken to evaluate in the type reference genome of HHV-6, HHV-6A strain U1102, populations of viruses generated from tissue culture passage or present from in vivo polymorphisms. This is essential to understanding the pathogenicity of HHV-6 as a HCMV.
Human herpesviruses cause a range of pathologies from oral diseases, incurable genital infections, to cancers. HHV-6A is linked with neurological and cardiac disorders in immune suppressed or naive patients [19], including fatal infant myocarditis [20,21]. Roseoloviruses, mainly HHV-6A and HHV-6B, are unique among human herpesviruses in integrating their genome in the human germline at the chromosomal telomere, termed CiHHV6A and CiHHV6B. This affects approximately 1% of people worldwide—there are upwards of 70 million people at risk of virus reactivation in every cell, and recent studies show links with impaired infant neurodevelopment and cardiac disease [14,22,23,24,25]. Genomic studies show the reference HHV-6A distinct from CiHHV6A, but the full extent of polymorphisms in the reference strain are required to understand relationships with CiHHV6A in order to interpret clinical conditions and potential treatment for virus reactivation.
Our studies using deep sequencing showed we could distinguish between known circulating HHV-6A viruses and those that are integrated in our genomes, CiHHV-6A [14,24]. With next generation sequencing technologies, we could amplify specifically the virus sequences separate from the human genome in order to characterise their differences. We combined this method with whole genome deep sequencing, where the short reads map to a reference template genome and are amplified thousands of times, which enables rare mutations in a population to be detected. We can significantly detect these “minor variants” down to 1% in a mixture and have used this method to detect minor variants from virus superinfection. Viruses superinfecting patients with integrated genomes were detected as 1–30% mixtures by using SNPs representative of superinfecting virus genomes [24].
This methodology was applied here to redefine the Roseolovirus type reference, HHV-6A strain U1102, using an early isolate and to characterise population polymorphisms and minor variants arising from effects of tissue culture passage of virus and in BACmid virus reconstruction. This was done in order to test whether the same mutagenic effects on the virus entry complex as shown for HCMV were also affected in HHV-6A. We further validate these polymorphisms by reference to in vivo virus “cloned” by integration in the germline. This not only defines the reference but also the changes critical to virus spread during cell culture that are relevant to general understanding of infection mechanisms.

2. Materials and Methods

2.1. Cells and Viruses

HHV-6A strain U1102 low passage from the original isolate, P3, was grown in JJhan cells as described [26]. The continuous cultured virus, termed pU1102 here, was grown on HSB2 cells (Kindly donated by B. Prusty from the University of Wurzberg, Würzburg, Germany) and a virus reconstructed from cloned BACmid version BAC-U1102 (kindly provided by Y. Mori, Kobe University, Japan and cultured by B. Prusty, University of Wurzberg, Würzburg, Germany) [12,27] were identically passaged in HSB2 cells, as described, for four serial weekly passages [28,29]. CiHHV-6A strains 5055, 2284, and 5814 were extracted from DNA from patients as described in [24].

2.2. Illumina Deep Sequencing

The five virus-infected cell DNA samples were extracted (Qiagen, Hilden, Germany) and CiHHV6A DNAs were prepared for target enrichment amplification and deep sequencing as described. These were HHV-6A U1102, bU1102, BAC-U1102, K21-HHV-6A-pU1102, and K21 BAC-U1102. For the specific amplification, 36 overlapping primer pairs were used which amplified overlapping PCR products across the complete genome as described, with all predicted sizes amplified (see primers in Table 1). These were separated by agarose gel electrophoresis, and the purified products pooled in equal quantities as measured using Qubit. Libraries were prepared by sonicating the DNA using Covaris, followed by purification and random PRC amplification using adaptor tags from the NEBNext DNA library kit. Tagged samples were run on a MiSeq and the resulting raw FASTQ files analysed after quality assessment and removal of adaptor tags and primer sequences as we described previously [24].
The reads were trimmed using trimmomatic with quality scores applied, then assembled using Samtools and Burrow-Wheeler Aligner (BWA)-mem alignment algorithm by mapping to reference genomes using the HHV-6A U1102 reference genome (updated from denovo assemblies) as described [24]. Coverage and qualities were assessed using BAC-U1102 complete coverage as a template. Contigs were ordered using the reference genome with manual adjustment using Artemis. Variant calling was applied with both Genome Analysis Toolkit (GATK) UnifiedGenotyper and SAMtools mpileup using BCFtools, vcfutils varFilter pipeline, quality scores >25, and minimum read depths of 10 for base call accuracy [24,29].
The rapid annotations tool Rapid Annotation Transfer Tool (RATT) was used to transfer annotations and syntenic coordinates from the reference genome [30] (HHV-6 U1102 NC_001664). All SNPs were compiled and compared in an excel database using a custom python script. Tables were tabulated for corrections and polymorphisms. All SNPs from K21 treated samples were compared using the unexposed HHV-6A pU1102 and BAC-U1102 samples as template. Average coverage was 2000–10,000 reads depth with a sensitivity variant cut-off of 1%. All results using target amplification were compared to those generated via SureSelect (Agilent, Santa Clara, CA, USA) enrichment and de novo assembly as described [31], and were in agreement with no evidence for gaps. The re-defined reference genome is deposited for European Bioinformatics Institute (EBI) (updated to the original deposition X83413) and Genbank (reference genome NC_001664).

2.3. Molecular Models

Comparative models were generated for HCMV gB from strain Merlin (UniProtKB/Swiss-Prot: F5HB53) and for HHV-6A U1102 gB from wt and passaged strains in the post-fusion conformation using Modeller9v7 [32]. The modelling was based on multiple templates from gB HCMV strain AD169 (PDB ID: 5cxf) and strain Towne (PDB ID: 4osn and 5c6t), gB Epstein-Barr virus (EBV) strain B95-8 (PDB ID: 3fvc), and gB Herpes Simplex virus (HSV)-1, strain KOS (PDB ID: 2gum and 3nwa). Closest alignments were shown with HCMV gB, 45% identity, 65% similarity (meta method M-Coffee, see Supplymentary Materials Figures S1–3). The models cover the ectodomain, residues 38–579, and domains I, II, III, and IV. The passage mutant, Thr193Ala, was modelled using swapaa command in Chimera [32]. The models of the HHV-6A gB trimer for the domain I of wt, and passage mutant, Thr193Ala, gB were obtained by structural alignments with the structure of the gB HCMV trimer strain AD169 (PDB ID: 5cxf) with each of its subunits. The free energy predictions of the Thr193Ala mutation on the subunit-subunit interaction stability was done using web server mCSM [33], which predicts stability changes of a wide range of mutations from graph-based signatures encoding distance patterns between atoms.

3. Results

3.1. Deep Sequence Analyses of Reference HHV-6A Strain U1102 Shows Corrections

Illumina deep sequencing was conducted deriving HHV-6A strain U1102 genomes between 2000–10,000 times per group. These included three versions of HHV-6A strain U1102: low passage (P3), continuous cell passage, and BAC reconstructed virus. In addition, sequences were derived for the continuous passage virus (pU1102) and the BAC virus (BAC-U1102). The Miseq run was derived at the same time as resequencing CiHHV6A strains 2284, 5055, and 5814 used to confirm SNPs. Consensus sequences were determined and minor variants tabulated and quantified.
The first analyses compared the NCBI Genbank reference genome HHV-6A strain U1102 genome to the deep sequenced early passage of the original isolate of strain U1102. In order to identify corrections to this reference sequence, all polymorphisms subsequently identified in the passaged strains were removed. Table 2 identifies 32 corrections to the reference sequence, including six previous citations, and one indel (U10). Nine of these included non-synonymous coding changes with five previously described changes (Table 2).
Newly identified changes were in genes U40, U57 (major capsid protein), U58, and U100 (gQ1 glycoprotein). It is possible these were mixtures of minor variants, below <1%, not detected in the original sequence. Therefore, these were also compared to genomes from cultured passaged virus, as well as from artificial cloned virus genomes using BACmid technology or naturally cloned virus genome present in the germline of people with chromosomally integrated virus genomes, CiHHV6. These were all re-sequenced at the same Illumina run, and none of these genomes showed evidence for these SNPs screened down to 1% sensitivity by deep sequencing.
Only five further SNPs were identified from a potential mixed original population. Two were previously described in U83 and gQ2 [13,14,27], but were not detected here down to 1% sensitivity. A further three SNPs comprised 88–97% of the population, but in non-coding repetitive sequences in the direct repeat (DR) region and the left end of the main unique region. The DR indels were also identified in the integrated genomes suggesting natural variants. This confirmed 32 corrections and 5 population variants out of 152,921 bp in the genome derived in the original isolate. Therefore, the original Sanger sequencing of the virus isolate was 99.98% accurate, which was not possible with Illumina next generation consensus sequencing methods, but resolved here with additional deep sequencing and analyses of population variants.

3.2. Changes with BACmid Cloning and Reconstruction

Having established a new reference genome sequence for HHV-6A strain U1102, this was then compared to the whole genome sequence of BACmid cloned virus. This had been cloned as BACmid using a BACmid cassette inserted in-between genes U53 and U54 [27], then the virus was reconstructed from tissue culture passage as described [28].
In total, 16 SNPs were identified including two indels (one in the U86 immendiate early 2 (IE2) repeat and one in an intron of U100 gQ1). Six of those were non-synonymous changes in structural coding genes U31 tegument, U33 capsid, U39 gB, U50 capsid, U86 IE2, and U100 gQ2 (Table 3). The SNP in gQ2 was distinct from a polymorphism previously reported which resulted in a premature stop codon giving a truncated gQ2 [12]. Instead, it was a single amino acid substitution (Table 3).
The only major rearrangement observed was at the site of the BAC cassette insertion. This cassette is derived from transfer cassette pHA2 [36] and includes loxP sites bounding a gpt, BAC related and green fluorescent protein (GFP) genes. This was rearranged, showing all the genes deleted with only one loxP site remaining indicative of recombination at this site. Moreover, the site of BAC insertion interrupts a newly identified distal poly A site for the essential U53/U53A proteinase/capsid scaffold genes (Figure 1). Therefore, mutagenic effects related to BAC reconstruction are demonstrated here.

3.3. Emerging Variant Genome Populations Generated after Tissue Culture Passage

Next, using deep sequencing, the new reference HHV-6A strain U1102 genome sequence was compared to genomes derived from virus, which had been continually passaged (four passages) as infected cells. The results showed that the passaged virus consisted now of an expanded population of minor variants with the now expanded major variants emerging to dominate (Table 3). Interestingly, the variants, which arose by culture, included all the five non-synonymous SNPs in structural genes identified in the BAC virus reconstructions. In addition, there were three further non-synonymous mutations in replication genes. In total, there were 24 SNPs, including six indels, all non-coding, four in DR, and one in the left end repeats of the Unique region, one in the R2 repeat and one in the intron within U100 gQ1. The 24 SNPs ranged from 52% to 99% of the mixed population.
Almost half of the SNPs, 11/24, were shared with the SNPs in the BACmid population of reconstructed virus genomes, therefore, indicating selection for tissue culture. These included the non-synonymous coding changes in U31, U33, U39 gB, U50 capsid, and U100 gQ2, which are mainly structural or envelope genes affecting infectivity. The new SNPs in this population included non-synonymous coding changes in U2, U41 DNA binding protein, U73 origin binding protein, U86 IE2, and three in U95. There were also two SNPs within the origin for replication. Of these, only three have been reported in another strain (GS). This strain is also a mixture from six original patients; two SNPs were reported for GS-KJ123690.1 (U41 and origin of replication, ORI, repeat) and two SNPs for GS-KC465951.1 (U50 capsid and the ORI repeat). The other non-coding changes are uniquely reported here, U2, U73 origin binding protein, U86 IE2, and the three in U95, which are mainly non-structural genes affecting replication. These may represent low-level variants present in the original virus isolate expanded here or the effects of adaptation to tissue culture and improved replication. Notably, all the BACmid virus SNP coding changes were present as variants in the continually passaged virus population, identifying these as most likely tissue culture adaptations.

3.4. Virus Envelope Genes Mutations from Tissue Culture Passage Affect Fusion Complex

Both the BACmid and cultured virus accumulated mutations in envelope glycoproteins of the entry complex, which mediate cellular fusion. This included glycoprotein gB, which has been characterised as a fusogen conserved in herpesvirus and also gQ1 and gQ2 which are part of the gH/gL/gQ1/gQ2 complex. Like gB, the gH/gL complexes are also conserved in herpesviruses, but there are additional trimer or quadrameric versions in Roseolovirus. The gH/gL/gQ1/gQ2 components are specific to Roseoloviruses; this complex interacts with a cellular receptor CD46, and is a target for a neutralising antibody, which inhibits cellular spread [11,41,42]. These genes are positional homologues of the UL128/130/131 pentameric complexes in HCMV which are also mutated by virus passage in cell culture [43]. The remaining complex gH/gL/gO is unaltered from tissue culture of both viruses. Although at present there is not a structural model for gQ1 or gQ2, the complex interacting with gH/gL can form the fusion complex as for other herpesviruses [12].
There are crystal structures available for homologous gB molecules from HCMV, HSV, and EBV (see Material and Methods) in the post-fusion conformation [40,44]. Using these structures, we modelled the tertiary structure of HHV-6A U1102 gB affected by the mutation Thr193Ala from cell passage (Table 3, Figure 2 and Figure 3). This mutation is not found in any homologous gB sequence analysed to date (BLAST Genbank scans). Alignments show that this mutation is located in the Beta 12 pleated sheet of Domain I (DI) of gB (based on the HCMV homolog) in the vicinity of the second fusion loop, FL2 (residues 184–191; Figure 2 and see Figures S1–3). Thr193 aligns to the unpaired Cys246 residue in HCMV gB. Further, DI not only is the fusion domain, but also the site suggested of interaction with gH/gL in potential triggering of cell fusion [40]. Thr193 is conserved in HSV and EBV gB next to a conserved motif central to the beta sheet (Figure 2 and Figure 3).
The passage mutation Thr193Ala is unlikely to affect the beta sheet prediction adjacent to FL2 (beta sheet 12 in HCMV domain I), but the reduced polar residues away from FL2 and increased hydrophobicity may promote membrane association via the Domain I structure as the mutation is within the FL2 adjacent beta-strand (Figure 2 and Figure 3). However, this residue seems to be interacting with residues at the trimer interfaces in Domain I (Leu93, Phe95, Ala147, Asp151, Thr191, and Leu195) (Figure 3 and Figure S2). Alternatively, the Thr193Ala mutation may modify the angle of inclination of the beta-sheet, increasing the interaction between the fusion loops and the membrane, although more likely stabilising the subunit interactions in the hydrophobic pocket as shown by increases in stability predictions (Figure S2).

4. Discussion

We used whole genome deep sequencing to re-define the reference genome of the Roseolovirus betaherpesvirus with unprecedented precision and further demonstrated emerging populations of minor variants from the effects of tissue culture passaging of virus and BACmid reconstructed virus. The results have implications for understanding cellular spread as well as clinical interpretation, analyses of drug resistance, and vaccine designs to inhibit virus entry. We analysed the original HHV-6A U1102 isolate in comparison to cell culture adapted virus and cell culture reconstructed BACmid virus. We also compared it to natural “clones” of HHV-6A in the integrated CiHHV-6A genomes we previously derived and re-sequenced here. Specifically, 32 SNP corrections were made to the reference genome sequence, which included nine corrections previously reported (Table 2) giving an overall original accuracy of 99.98%, higher than that of consensus sequencing by Illumina which is 99%. However, combined with deep sequencing of the original isolate, this allowed enhanced resolution as well as definition of emerging minor variants in the cell culture passaged virus, and many of these dominated the mixed population and were detected as “major” variants. There were mixtures of variants in both the passaged HHV-6A U1102 and BAC-U1102 virus stocks with 34–99% mixtures compared to the redefined reference genome. In the passaged stocks, HHV-6A pU1102 and BAC-U1102, there were 10 and 6 genes with non-synonymous SNPs, respectively. Therefore, both the HHV-6 pU1102 and BAC-HHV-6-U1102 stock viruses had polymorphic mixtures. This also is the first derivation of the genome of the reconstructed BAC-HHV-6-U1102 and the polymorphisms shown here comprise a subset of those in the passaged virus, identifying cell culture specific effects.
The passaged stocks were shown comprised of complex polymorphic mixtures, with SNPs averaging 50% mixtures with the wild type reference virus genome. We used screens to the 85 defined coding genes as described for both virus and integrated genomes [24,26]. There were SNP mixtures in both the passaged virus stock as well as the cloned BACmid virus stocks. The passaged virus had 24 SNPs identified.
Mixtures and variants of the strains of HHV-6A as well as related HHV-6B have been described previously. In these cases, the dominant SNP created the consensus sequence. The initial virus isolates from patients were not clonal and since the virus is cell associated, the virus strains were not clone purified previously. For example, there are at least two consensus genomes for HHV-6A strain GS [47] (HHV-6A GS KJ123,690.1 and KC465,951.1).
Using BACmid technology, which “clones” the virus genome using a bacterial artificial chromosome sequence, BAC, with plasmid elements, this should clonally purify the virus stock. However, recovered virus has to replicate in cell culture, and studies from another herpesvirus, HCMV, show that complex mixtures of SNPs can accumulate during selection for replication in tissue culture. In the case of the only described HHV-6A BACmid, BAC-U1102, we also show here there are SNPs compared to the wild type genome. Some could be present in the virus stock used to generate this BAC, but the mixtures shown here also indicate there is some selection during growth in culture.
In the original description of the derivation of BAC HHV-6A U1102, four clones showed varying BAMH1 restriction endonuclease digestion profiles, indicating either mixtures present in the parent U1102 stock used or mutagenesis during the derivation of the BACmid clones. Notably, only one of these clones, G1, was able to reconstitute virus after transfection and co-culture with cord blood mononuclear cells (CBMC). This BAMH1 restriction enzyme profile most closely resembles that reported for the original isolate [19,48]. However, reports using this BAC cloned virus have observed low replication efficiency, so it is possible that the SNPs identified here contribute to this phenotype. It is notable that there are fewer SNPs in the parent BACmid stock, 16 SNPs compared to 24 in the passaged HHV-6 U1102 stock, which could indicate that continuous passage of virus in culture permits accumulation of further SNPs.
Rearrangements in the BACmid may also affect virus replication efficiency. The BAC virus genomes had deletion of the GFP gene. Only the recombined loxP site and associated sequences from the original vectors were maintained in both the parent BAC and the K21 treated BAC virus. Although this could be from homologous recombination mediated by the host, the efficient rearrangement in the entire population suggests a virus gene product is mediating this event. This would confound any selection process dependent on GFP, since mutations or rearrangements would have to occur to reduce this activity. Since GFP expression has been used to characterise genome integration, it is likely that selection for GFP may be artificially enhancing integration as demonstrated [49,50,51]. Analyses showed that the GFP insertion disrupts utilisation of the downstream polyadenylation site, which could be a more efficient transcript. The upstream polyadenylation signal is prior to the stop codon, and with cleavage 10–30 nucleotides downstream, this could result in transcript cleavage and polyadenylation before the stop codon, resulting in an unstable transcript. This could affect the interpretation of studies using BACmids to study replication or integrations. For example, if GFP is selected, this could select for virus with inefficient virus replication, and promote virus genome integration. Further, the remaining bacterial origin sequences may affect replication and integration. The GFP insertion may interrupt the alternate downstream polyadenylation site, giving inefficient expression of the essential scaffold protein p53A and lower virions. Therefore, use of this BACmid for replication or integration studies can give confounding results and either new non-rearranging BACmid or different approaches are required to understand these pathways.
The SNPs in BACmid reconstructed virus are shared with those in the passaged virus, showing effects of cell culture passage, and all are in structural genes. In comparison, the virus-passaged genomes had only three additional SNPs, all related to replication, in the origin binding protein, and the immediate early genes IE2 and U95 “IE3”. There were five shared structural changes, two in capsid and tegument genes, U31 and U33, with the remaining three in envelope glycoprotein genes affecting virus cellular entry by cell fusion. These included gB and the gH/gL complex components gQ1 and gQ2. These may be culture adapted changes for improved growth in leukemic cell lines JJhan, SupT1, and HSB2, which may affect the selection. In contrast, clinical isolates were prepared on primary cord blood stem cells. The population of BACmid SNPs were in structural genes and the virus envelope and these, including gB Thr193Ala, were all contained in the passaged population of SNPs, which also included SNPs in essential replication genes. Of particular relevance, we showed treatment of BAC virus with a novel antiviral affected membranes, reverted the population to wild type, and eliminated the cell culture passage SNPs in the envelope glycoproteins [29].
In HCMV, cell culture causes changes in the virus entry complex, specifically in components of the pentameric complex of gH/gL, in UL128, UL130, and UL131. In HHV-6A, the genes in the homologous genomic position also encode a fusion complex alternate to the trimeric gH/gL/gO. This is the quadrameric complex of gH/gL including components gQ1 and gQ2. In addition, there is a chemokine gene, a virokine, U83A, which is the positional homologue of UL128, but appears to independently direct cellular recruitment for dissemination or immune evasion. There is polymeric DNA editing polymorphism in this gene disrupting expression, the full-length gene is only rarely identified in culture, but has been identified in “cloned” integrated genomes [14]. Therefore, effects of cell passage of both virus and BAC virus show similarities to the mutations arising in HCMV after cell culture passage. The gH/gL/gQ1/gQ2 complex binds the receptor CD46, and it appears that in HHV-6A cell culture, passage favours mutation of the alternate complex to gH/gL/gO, which may also distort cellular tropism and immune responses as shown for CMV. Distinct from HCMV, we show virus and BAC-virus laboratory passage also results in mutation of the fusogen gB, in addition to the fusion regulator, the gH/gL complexes.
HHV-6A gB is closely related to HCMV gB, which enabled modelling of the ectodomain in the post-fusion conformation with high confidence. This demonstrated the mutation, Thr193Ala, was in Domain I, which mediates cell fusion, specifically in proximity to the fusion loop 2 (FL2) adjacent region within a conserved beta sheet. Interestingly, recent studies identified a cell culture induced mutation in the Ebola virus GP1,2 glycoprotein, which increases virus cellular spread. This is also in the fusion loop adjacent region, a Thr544Ile mutation, which also increases hydrophobicity [52]. The modelling of HHV-6 gB mutation showed that the mutation is located at the predicted trimer interface region, which may increase the stability of the post-fusion conformation through increased hydrophobicity (Figure S2). Interestingly, in treatment with a membrane directed antiviral, the treated virus no longer retained the gB mutation, and comparisons with the mutated passaged virus in immune precipitation analyses did show increased stability of gB [28,29].
The cell culture passage mutation in the FL2 proximal region is within domain I, the site of effects by gH/gL based on experimental results, although no stable intermediate has been identified [40]. In HSV gH, the N-terminal region proposed to affect gB overlaps with the site of gL interaction and a complex conformational epitope, interacting with a neutralising LP11 antibody [53,54]. A similar structure has been demonstrated at the N-terminus of HHV-6A gH, with mutational analyses demonstrating that a conformational change between interactions of the N and C terminal domains is required for cell fusion [8]. Campadelli-Fiume and colleagues have demonstrated that for the HSV entry complex, gL removal is required for fusion [55]. Therefore, the cell passage mutation in the HHV-6A FL2 proximal region in gB domain I could potentially affect this cell fusion regulation favouring cell-associated virus spread. Recent studies on the gB fusion loops indicate the FL2 proximal region may influence the angle of membrane association [56].
Taken together, our overall studies have redefined the reference genome and identified SNPs associated with cell culture and BAC virus passage, which can affect virus infection by cell fusion as similarly demonstrated in the related betaherpesvirus HCMV. Characterisation of host interactions, genome integration, immune interactions, and antiviral resistance all need to take account of the effects of a diverse mixed genome population present in laboratory-passaged virus and highlight the components of the virus entry complex affecting cellular spread.

5. Conclusions

The reference virus genome was re-defined by deep sequencing of the early clinical isolate. Tissue culture passaged virus accumulated a population of genetic variants in genes affecting infection and replication. Virus culture of BACmid virus accumulated a subset of these SNPs in genes affecting infection, including envelope genes gB and gQ1 and gQ2 of the gH/gL complex. Efforts to study virus properties need to focus on the newly defined reference strain. New BACs or different culturing protocols may be required to do this. Additional study of these cell associated mutations from BAC and passaged virus would provide further understanding of cell associated infection and regulation of virus production versus integration.

Supplementary Materials

The following are available online at www.mdpi.com/1999-4915/10/1/16/s1, Figure S1: Alignment HHV-6A strain U1102 and HCMV strain Merlin gB ectodomains; Figure S2: The passage mutation in HHV-6 gB incrases interface hydrophobicity affecting subunit interactions; Figure S3: Multiple sequence alignment of HHV-6A gB with model templates from HCMV, HSV and EBV.

Acknowledgments

We thank the HHV-6 Foundation, LSHTM GTA and Bloomsbury PhD Studentships (Joshua G. Tweedy and Eric Escriva), and MRC project grant MR/M019292/1 (Maya Topf) for support.

Author Contributions

Joshua G. Tweedy and Ursula A. Gompels designed the study, Joshua G. Tweedy performed experimentation, Eric Escriva, Maya Topf, and Ursula A. Gompels constructed protein alignments and models, Joshua G. Tweedy and Ursula A. Gompels analysed the genomics, all contributed to manuscript preparation/review, and Ursula A. Gompels wrote the final paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chee, M.S.; Bankier, A.T.; Beck, S.; Bohni, R.; Brown, C.M.; Cerny, R.; Horsnell, T.; Hutchison, C.A., 3rd; Kouzarides, T.; Martignetti, J.A.; et al. Analysis of the protein-coding content of the sequence of human cytomegalovirus strain AD169. Curr. Top. Microbiol. Immunol. 1990, 154, 125–169. [Google Scholar] [PubMed]
  2. Cha, T.A.; Tom, E.; Kemble, G.W.; Duke, G.M.; Mocarski, E.S.; Spaete, R.R. Human cytomegalovirus clinical isolates carry at least 19 genes not found in laboratory strains. J. Virol. 1996, 70, 78–83. [Google Scholar] [PubMed]
  3. Ciferri, C.; Chandramouli, S.; Donnarumma, D.; Nikitin, P.A.; Cianfrocco, M.A.; Gerrein, R.; Feire, A.L.; Barnett, S.W.; Lilja, A.E.; Rappuoli, R.; et al. Structural and biochemical studies of HCMV gH/gL/gO and Pentamer reveal mutually exclusive cell entry complexes. Proc. Natl. Acad. Sci. USA 2015, 112, 1767–1772. [Google Scholar] [CrossRef] [PubMed]
  4. Murrell, I.; Wilkie, G.S.; Davison, A.J.; Statkute, E.; Fielding, C.A.; Tomasec, P.; Wilkinson, G.W.; Stanton, R.J. Genetic Stability of Bacterial Artificial Chromosome-Derived Human Cytomegalovirus during Culture In Vitro. J. Virol. 2016, 90, 3929–3943. [Google Scholar] [CrossRef] [PubMed]
  5. Ciferri, C.; Chandramouli, S.; Leitner, A.; Donnarumma, D.; Cianfrocco, M.A.; Gerrein, R.; Friedrich, K.; Aggarwal, Y.; Palladino, G.; Aebersold, R.; et al. Antigenic Characterization of the HCMV gH/gL/gO and Pentamer Cell Entry Complexes Reveals Binding Sites for Potently Neutralizing Human Antibodies. PLoS Pathog. 2015, 11, e1005230. [Google Scholar] [CrossRef] [PubMed]
  6. Hansen, S.G.; Sacha, J.B.; Hughes, C.M.; Ford, J.C.; Burwitz, B.J.; Scholz, I.; Gilbride, R.M.; Lewis, M.S.; Gilliam, A.N.; Ventura, A.B.; et al. Cytomegalovirus vectors violate CD8+ T cell epitope recognition paradigms. Science 2013, 340, 1237874. [Google Scholar] [CrossRef] [PubMed]
  7. Paterson, D.A.; Dyer, A.P.; Milne, R.S.; Sevilla-Reyes, E.; Gompels, U.A. A role for human cytomegalovirus glycoprotein O (gO) in cell fusion and a new hypervariable locus. Virology 2002, 293, 281–294. [Google Scholar] [CrossRef] [PubMed]
  8. Anderson, R.A.; Gompels, U.A. N- and C-terminal external domains of human herpesvirus-6 glycoprotein H affect a fusion-associated conformation mediated by glycoprotein L binding the N terminus. J. Gen. Virol. 1999, 80, 1485–1494. [Google Scholar] [CrossRef] [PubMed]
  9. Anderson, R.A.; Liu, D.X.; Gompels, U.A. Definition of a human herpesvirus-6 betaherpesvirus-specific domain in glycoprotein gH that governs interaction with glycoprotein gL: Substitution of human cytomegalovirus glycoproteins permits group-specific complex formation. Virology 1996, 217, 517–526. [Google Scholar] [CrossRef] [PubMed]
  10. Jasirwan, C.; Furusawa, Y.; Tang, H.; Maeki, T.; Mori, Y. Human herpesvirus-6A gQ1 and gQ2 are critical for human CD46 usage. Microbiol. Immunol. 2014, 58, 22–30. [Google Scholar] [CrossRef] [PubMed]
  11. Mori, Y.; Yang, X.; Akkapaiboon, P.; Okuno, T.; Yamanishi, K. Human herpesvirus 6 variant A glycoprotein H-glycoprotein L-glycoprotein Q complex associates with human CD46. J. Virol. 2003, 77, 4992–4999. [Google Scholar] [CrossRef] [PubMed]
  12. Tang, H.; Hayashi, M.; Maeki, T.; Yamanishi, K.; Mori, Y. Human herpesvirus 6 glycoprotein complex formation is required for folding and trafficking of the gH/gL/gQ1/gQ2 complex and its cellular receptor binding. J. Virol. 2011, 85, 11121–11130. [Google Scholar] [CrossRef] [PubMed]
  13. Dewin, D.R.; Catusse, J.; Gompels, U.A. Identification and characterization of U83A viral chemokine, a broad and potent beta-chemokine agonist for human CCRs with unique selectivity and inhibition by spliced isoform. J. Immunol. 2006, 176, 544–556. [Google Scholar] [CrossRef] [PubMed]
  14. Tweedy, J.; Spyrou, M.A.; Hubacek, P.; Kuhl, U.; Lassner, D.; Gompels, U.A. Analyses of germline, chromosomally integrated human herpesvirus 6A and B genomes indicate emergent infection and new inflammatory mediators. J. Gen. Virol. 2015, 96, 370–389. [Google Scholar] [CrossRef] [PubMed]
  15. French, C.; Menegazzi, P.; Nicholson, L.; Macaulay, H.; DiLuca, D.; Gompels, U.A. Novel, nonconsensus cellular splicing regulates expression of a gene encoding a chemokine-like protein that shows high variation and is specific for human herpesvirus 6. Virology 1999, 262, 139–151. [Google Scholar] [CrossRef] [PubMed]
  16. Schleimann, M.H.; Hoberg, S.; Solhoj Hansen, A.; Bundgaard, B.; Witt, C.T.; Kofod-Olsen, E.; Hollsberg, P. The DR6 protein from human herpesvirus-6B induces p53-independent cell cycle arrest in G2/M. Virology 2014, 452–453, 254–263. [Google Scholar] [CrossRef] [PubMed]
  17. Papanikolaou, E.; Kouvatsis, V.; Dimitriadis, G.; Inoue, N.; Arsenakis, M. Identification and characterization of the gene products of open reading frame U86/87 of human herpesvirus 6. Virus Res. 2002, 89, 89–101. [Google Scholar] [CrossRef]
  18. Salahuddin, S.Z.; Ablashi, D.V.; Markham, P.D.; Josephs, S.F.; Sturzenegger, S.; Kaplan, M.; Halligan, G.; Biberfeld, P.; Wong-Staal, F.; Kramarsky, B.; et al. Isolation of a new virus, HBLV, in patients with lymphoproliferative disorders. Science 1986, 234, 596–601. [Google Scholar] [CrossRef] [PubMed]
  19. Gompels, U.A. Human herpesvirus 6A, 6B and 7 (Herpesviridae). In Reference Module in Biomedical Sciences, Encyclopedia of Virology, 3rd ed.; Mahy, B., Ed.; Elsevier: Amsterdam, The Netherlands, 2016; pp. 498–505. [Google Scholar]
  20. Simpson, K.E.; Storch, G.A.; Lee, C.K.; Ward, K.E.; Danon, S.; Simon, C.M.; Delaney, J.W.; Tong, A.; Canter, C.E. High Frequency of Detection by PCR of Viral Nucleic Acid in The Blood of Infants Presenting with Clinical Myocarditis. Pediatr. Cardiol. 2016, 37, 399–404. [Google Scholar] [CrossRef] [PubMed]
  21. Stefanski, H.E.; Thibert, K.A.; Pritchett, J.; Prusty, B.K.; Wagner, J.E.; Lund, T.C. Fatal Myocarditis Associated With HHV-6 Following Immunosuppression in Two Children. Pediatrics 2016, 137. [Google Scholar] [CrossRef] [PubMed]
  22. Gravel, A.; Dubuc, I.; Morissette, G.; Sedlak, R.H.; Jerome, K.R.; Flamand, L. Inherited chromosomally integrated human herpesvirus 6 as a predisposing risk factor for the development of angina pectoris. Proc. Natl. Acad. Sci. USA 2015, 112, 8058–8063. [Google Scholar] [CrossRef] [PubMed]
  23. Kuhl, U.; Lassner, D.; Wallaschek, N.; Gross, U.M.; Krueger, G.R.; Seeberg, B.; Kaufer, B.B.; Escher, F.; Poller, W.; Schultheiss, H.P. Chromosomally integrated human herpesvirus 6 in heart failure: Prevalence and treatment. Eur. J. Heart Fail. 2015, 17, 9–19. [Google Scholar] [CrossRef] [PubMed]
  24. Tweedy, J.; Spyrou, M.A.; Pearson, M.; Lassner, D.; Kuhl, U.; Gompels, U.A. Complete Genome Sequence of Germline Chromosomally Integrated Human Herpesvirus 6A and Analyses Integration Sites Define a New Human Endogenous Virus with Potential to Reactivate as an Emerging Infection. Viruses 2016, 8. [Google Scholar] [CrossRef] [PubMed]
  25. Caserta, M.T.; Hall, C.B.; Canfield, R.L.; Davidson, P.; Lofthus, G.; Schnabel, K.; Carnahan, J.; Shelley, L.; Wang, H. Early developmental outcomes of children with congenital HHV-6 infection. Pediatrics 2014, 134, 1111–1118. [Google Scholar] [CrossRef] [PubMed]
  26. Gompels, U.A.; Nicholas, J.; Lawrence, G.; Jones, M.; Thomson, B.J.; Martin, M.E.; Efstathiou, S.; Craxton, M.; Macaulay, H.A. The DNA sequence of human herpesvirus-6: Structure, coding content, and genome evolution. Virology 1995, 209, 29–51. [Google Scholar] [CrossRef] [PubMed]
  27. Tang, H.; Kawabata, A.; Yoshida, M.; Oyaizu, H.; Maeki, T.; Yamanishi, K.; Mori, Y. Human herpesvirus 6 encoded glycoprotein Q1 gene is essential for virus growth. Virology 2010, 407, 360–367. [Google Scholar] [CrossRef] [PubMed]
  28. Gulve, N.; Kimmerling, K.; Johnston, A.D.; Krueger, G.R.; Ablashi, D.V.; Prusty, B.K. Anti-herpesviral effects of a novel broad range anti-microbial quaternary ammonium silane, K21. Antivir. Res. 2016, 131, 166–173. [Google Scholar] [CrossRef] [PubMed]
  29. Tweedy, J.G.; Prusty, B.K.; Gompels, U.A. Use of whole genome deep sequencing to define emerging minority variants in virus envelope genes in herpesvirus treated with novel antimicrobial K21. Antivir. Res. 2017, 146, 201–204. [Google Scholar] [CrossRef] [PubMed]
  30. Otto, T.D.; Dillon, G.P.; Degrave, W.S.; Berriman, M. RATT: Rapid Annotation Transfer Tool. Nucleic Acids Res. 2011, 39, e57. [Google Scholar] [CrossRef] [PubMed]
  31. Tweedy, J.; Spyrou, M.A.; Donaldson, C.D.; Depledge, D.; Breuer, J.; Gompels, U.A. Complete Genome Sequence of the Human Herpesvirus 6A Strain AJ from Africa Resembles Strain GS from North America. Genome Announc. 2015, 3. [Google Scholar] [CrossRef] [PubMed]
  32. Yang, Z.; Lasker, K.; Schneidman-Duhovny, D.; Webb, B.; Huang, C.C.; Pettersen, E.F.; Goddard, T.D.; Meng, E.C.; Sali, A.; Ferrin, T.E. UCSF Chimera, MODELLER, and IMP: An integrated modeling system. J. Struct. Biol. 2012, 179, 269–278. [Google Scholar] [CrossRef] [PubMed]
  33. Pires, D.E.V.; Ascher, D.B.; Blundell, T.L. mCSM: Predicting the effects of mutations in proteins using graph-based signatures. Bioinformatics 2014, 30, 335–342. [Google Scholar] [CrossRef] [PubMed]
  34. Neipel, F.; Ellinger, K.; Fleckenstein, B. Gene for the major antigenic structural protein (p100) of human herpesvirus 6. J. Virol. 1992, 66, 3918–3924. [Google Scholar] [PubMed]
  35. Nicholas, J.; Martin, M.E. Nucleotide sequence analysis of a 38.5-kilobase-pair region of the genome of human herpesvirus 6 encoding human cytomegalovirus immediate-early gene homologs and transactivating functions. J. Virol. 1994, 68, 597–610. [Google Scholar] [PubMed]
  36. Adler, H.; Messerle, M.; Wagner, M.; Koszinowski, U.H. Cloning and mutagenesis of the murine gammaherpesvirus 68 genome as an infectious bacterial artificial chromosome. J. Virol. 2000, 74, 6964–6974. [Google Scholar] [CrossRef] [PubMed]
  37. Thompson, J.; Choudhury, S.; Kashanchi, F.; Doniger, J.; Berneman, Z.; Frenkel, N.; Rosenthal, L.J. A transforming fragment within the direct repeat region of human herpesvirus type 6 that transactivates HIV-1. Oncogene 1994, 9, 1167–1175. [Google Scholar] [PubMed]
  38. Kondo, K.; Shimada, K.; Sashihara, J.; Tanaka-Taya, K.; Yamanishi, K. Identification of human herpesvirus 6 latency-associated transcripts. J. Virol. 2002, 76, 4145–4151. [Google Scholar] [CrossRef] [PubMed]
  39. Akkapaiboon, P.; Mori, Y.; Sadaoka, T.; Yonemoto, S.; Yamanishi, K. Intracellular processing of human herpesvirus 6 glycoproteins Q1 and Q2 into tetrameric complexes expressed on the viral envelope. J. Virol. 2004, 78, 7969–7983. [Google Scholar] [CrossRef] [PubMed]
  40. Burke, H.G.; Heldwein, E.E. Crystal Structure of the Human Cytomegalovirus Glycoprotein B. PLoS Pathog. 2015, 11, e1005227. [Google Scholar] [CrossRef] [PubMed]
  41. Santoro, F.; Kennedy, P.E.; Locatelli, G.; Malnati, M.S.; Berger, E.A.; Lusso, P. CD46 is a cellular receptor for human herpesvirus 6. Cell 1999, 99, 817–827. [Google Scholar] [CrossRef]
  42. Hansen, A.S.; Bundgaard, B.B.; Biltoft, M.; Rossen, L.S.; Hollsberg, P. Divergent tropism of HHV-6AGS and HHV-6BPL1 in T cells expressing different CD46 isoform patterns. Virology 2017, 502, 160–170. [Google Scholar] [CrossRef] [PubMed]
  43. Dolan, A.; Cunningham, C.; Hector, R.D.; Hassan-Walker, A.F.; Lee, L.; Addison, C.; Dargan, D.J.; McGeoch, D.J.; Gatherer, D.; Emery, V.C.; et al. Genetic content of wild-type human cytomegalovirus. J. Gen. Virol. 2004, 85, 1301–1312. [Google Scholar] [CrossRef] [PubMed]
  44. Zeev-Ben-Mordehai, T.; Vasishtan, D.; Hernandez Duran, A.; Vollmer, B.; White, P.; Prasad Pandurangan, A.; Siebert, C.A.; Topf, M.; Grunewald, K. Two distinct trimeric conformations of natively membrane-anchored full-length herpes simplex virus 1 glycoprotein B. Proc. Natl. Acad. Sci. USA 2016, 113, 4176–4181. [Google Scholar] [CrossRef] [PubMed]
  45. Achour, A.; Malet, I.; Le Gal, F.; Dehee, A.; Gautheret-Dejean, A.; Bonnafous, P.; Agut, H. Variability of gB and gH genes of human herpesvirus-6 among clinical specimens. J. Med. Virol. 2008, 80, 1211–1221. [Google Scholar] [CrossRef] [PubMed]
  46. Franti, M.; Aubin, J.T.; Gautheret-Dejean, A.; Malet, I.; Cahour, A.; Huraux, J.M.; Agut, H. Preferential associations of alleles of three distinct genes argue for the existence of two prototype variants of human herpesvirus 7. J. Virol. 1999, 73, 9655–9658. [Google Scholar] [PubMed]
  47. Gravel, A.; Ablashi, D.; Flamand, L. Complete Genome Sequence of Early Passaged Human Herpesvirus 6A (GS Strain) Isolated from North America. Genome Announc. 2013, 1. [Google Scholar] [CrossRef] [PubMed]
  48. Martin, M.E.; Thomson, B.J.; Honess, R.W.; Craxton, M.A.; Gompels, U.A.; Liu, M.Y.; Littler, E.; Arrand, J.R.; Teo, I.; Jones, M.D. The genome of human herpesvirus 6: Maps of unit-length and concatemeric genomes for nine restriction endonucleases. J. Gen. Virol. 1991, 72, 157–168. [Google Scholar] [CrossRef] [PubMed]
  49. Gulve, N.; Frank, C.; Klepsch, M.; Prusty, B.K. Chromosomal integration of HHV-6A during non-productive viral infection. Sci. Rep. 2017, 7, 512. [Google Scholar] [CrossRef] [PubMed]
  50. Wallaschek, N.; Gravel, A.; Flamand, L.; Kaufer, B.B. The putative U94 integrase is dispensable for human herpesvirus 6 (HHV-6) chromosomal integration. J. Gen. Virol. 2016, 97, 1899–1903. [Google Scholar] [CrossRef] [PubMed]
  51. Wallaschek, N.; Sanyal, A.; Pirzer, F.; Gravel, A.; Mori, Y.; Flamand, L.; Kaufer, B.B. The Telomeric Repeats of Human Herpesvirus 6A (HHV-6A) Are Required for Efficient Virus Integration. PLoS Pathog. 2016, 12, e1005666. [Google Scholar] [CrossRef] [PubMed]
  52. Ruedas, J.B.; Ladner, J.T.; Ettinger, C.R.; Gummuluru, S.; Palacios, G.; Connor, J.H. Spontaneous Mutation at Amino Acid 544 of the Ebola Virus Glycoprotein Potentiates Virus Entry and Selection in Tissue Culture. J. Virol. 2017, 91. [Google Scholar] [CrossRef] [PubMed]
  53. Gompels, U.A.; Carss, A.L.; Saxby, C.; Hancock, D.C.; Forrester, A.; Minson, A.C. Characterization and sequence analyses of antibody-selected antigenic variants of herpes simplex virus show a conformationally complex epitope on glycoprotein H. J. Virol. 1991, 65, 2393–2401. [Google Scholar] [PubMed]
  54. Chowdary, T.K.; Cairns, T.M.; Atanasiu, D.; Cohen, G.H.; Eisenberg, R.J.; Heldwein, E.E. Crystal structure of the conserved herpesvirus fusion regulator complex gH-gL. Nat. Struct. Mol. Biol. 2010, 17, 882–888. [Google Scholar] [CrossRef] [PubMed]
  55. Gianni, T.; Massaro, R.; Campadelli-Fiume, G. Dissociation of HSV gL from gH by alphavbeta6- or alphavbeta8-integrin promotes gH activation and virus entry. Proc. Natl. Acad. Sci. USA 2015, 112, E3901–E3910. [Google Scholar] [CrossRef] [PubMed]
  56. Vallbracht, M.; Brun, D.; Tassinari, M.; Vaney, M.C.; Pehau-Arnaudet, G.; Guardado-Calvo, P.; Haouz, A.; Klupp, B.G.; Mettenleiter, T.C.; Rey, F.A.; et al. Structure-function dissection of the Pseudorabies virus glycoprotein B fusion loops. J. Virol. 2017. [CrossRef] [PubMed]
Viruses 10 00016 g001 550
Figure 1. Passaged bacterial artificial chromosome (BAC)-virus deletes the green fluorescent protein (GFP) cassette. In the BAC virus, the GFP cassette disrupts the U53/U53a transcriptional control elements, and after tissue culture passage it is deleted. Shaded and underlined are palindrome1, palindrome 2, and LoxP sites; note the insertion is in an end repaired EcoRI site. Underlined are the original and new poly A regulatory sequences disrupted by the BAC cassette.
Figure 1. Passaged bacterial artificial chromosome (BAC)-virus deletes the green fluorescent protein (GFP) cassette. In the BAC virus, the GFP cassette disrupts the U53/U53a transcriptional control elements, and after tissue culture passage it is deleted. Shaded and underlined are palindrome1, palindrome 2, and LoxP sites; note the insertion is in an end repaired EcoRI site. Underlined are the original and new poly A regulatory sequences disrupted by the BAC cassette.
Viruses 10 00016 g001
Viruses 10 00016 g002 550
Figure 2. Human herpesvirus gB fusion loop 2, FL2, region sequence comparisons show unique SNP substitution in passaged HHV-6A virus. Comparison of the FL2 (in bold) and associated sequences of HHV-6A gB (residues 184–191 in FL2 in bold) from passaged virus (and cultured BAC virus) compared to gB from 51 clinical samples of HHV-6A/B [45], CiHHV-6A/B [14,24], and 50 clinical samples of HHV-7 gB [46] demonstrating the Thr193Ala mutation is only found in passaged virus. This is compared to reference sequences from alpha, beta, and gamma herpesvirus gB where the structure has been determined, as shown for HCMV gB structure, identifying the beta sheet 12 in domain I [40]. Alignments were made using Clustal as described [40] with gB from HCMV (shown for strains Merlin, Towne, and TB40E and 74 clinical strains; strain AD169 and six clinical strains have Val252Leu), EBV (strain B958), and HSV-1 (strain KOS) (UniProtKB P06473, Q777B0, P06437, respectively). The FL2 in HCMV, EBV, and HSV-1 gB are in bold, and predicted FL2 for HHV-6A/B and HHV-7 are similarly indicated. FL2 in HHV-6A/B gB is amino acids 184–191 and in HCMV, strain Merlin, 237–244. This is adjacent to a conserved beta sheet shown in the structures of EBV, HSV, and cytomegalovirus (CMV) gB as underlined, and that predicted for HHV-6A/B and HHV-7 indicated (See Methods). The Thr193Ala passage SNP is shaded in the FL2 adjacent region, disrupting the polar patch and increasing hydrophobicity. This Thr is conserved in HHV-6A/B, HHV-7, EBV, and HSV.
Figure 2. Human herpesvirus gB fusion loop 2, FL2, region sequence comparisons show unique SNP substitution in passaged HHV-6A virus. Comparison of the FL2 (in bold) and associated sequences of HHV-6A gB (residues 184–191 in FL2 in bold) from passaged virus (and cultured BAC virus) compared to gB from 51 clinical samples of HHV-6A/B [45], CiHHV-6A/B [14,24], and 50 clinical samples of HHV-7 gB [46] demonstrating the Thr193Ala mutation is only found in passaged virus. This is compared to reference sequences from alpha, beta, and gamma herpesvirus gB where the structure has been determined, as shown for HCMV gB structure, identifying the beta sheet 12 in domain I [40]. Alignments were made using Clustal as described [40] with gB from HCMV (shown for strains Merlin, Towne, and TB40E and 74 clinical strains; strain AD169 and six clinical strains have Val252Leu), EBV (strain B958), and HSV-1 (strain KOS) (UniProtKB P06473, Q777B0, P06437, respectively). The FL2 in HCMV, EBV, and HSV-1 gB are in bold, and predicted FL2 for HHV-6A/B and HHV-7 are similarly indicated. FL2 in HHV-6A/B gB is amino acids 184–191 and in HCMV, strain Merlin, 237–244. This is adjacent to a conserved beta sheet shown in the structures of EBV, HSV, and cytomegalovirus (CMV) gB as underlined, and that predicted for HHV-6A/B and HHV-7 indicated (See Methods). The Thr193Ala passage SNP is shaded in the FL2 adjacent region, disrupting the polar patch and increasing hydrophobicity. This Thr is conserved in HHV-6A/B, HHV-7, EBV, and HSV.
Viruses 10 00016 g002
Viruses 10 00016 g003 550
Figure 3. Structural models of HHV-6A gB show conservation with HCMV gB structure identifying the passaged virus has a mutation in the fusion associated domain I affecting the trimer interface. (a) Model of the HCMV reference strain Merlin (ectodomain in residues 91–642 from a total of 907), alignment with corresponding region of HHV-6A gB is shown in the Figures S1–3. (b) Model of HHV-6A gB ectodomain, in the corresponding residues 38 to 579 (total 839 residues) using multiple templates of pdb (RCSB protein database) folds from gB HCMV strain AD169 (PDB ID: 5cxf) and strain Towne (PDB IDs: 4osn, 5c6t), gB EBV strain B95-8 (PDB ID: 3fvc), and gB HSV-1 strain KOS (PDB IDs: 2gum and 3nwa). Using the structure of most closely related HCMV gB, as a model, shows Domains I (blue), II (green), III (yellow), and IV (orange) in the post-fusion conformation. Domain I (blue) shows conserved beta sheets connected by the two fusion loops, FL1 and FL2, at the base of the model (purple). (c) Longitudinal view of the predicted trimer of HHV-6A gB showing the Thr193 position in red within three subunits, indicated in purple, blue, and orange. (d) Domain I (residues 79–291) is enlarged as indicated by the grey T bar, showing the position of Thr193, in red, in the beta sheet at the trimer interfaces with the fusion loops FL1 (residues 100–105) and FL2 (residues 184–191) exposed at the base of the structure.
Figure 3. Structural models of HHV-6A gB show conservation with HCMV gB structure identifying the passaged virus has a mutation in the fusion associated domain I affecting the trimer interface. (a) Model of the HCMV reference strain Merlin (ectodomain in residues 91–642 from a total of 907), alignment with corresponding region of HHV-6A gB is shown in the Figures S1–3. (b) Model of HHV-6A gB ectodomain, in the corresponding residues 38 to 579 (total 839 residues) using multiple templates of pdb (RCSB protein database) folds from gB HCMV strain AD169 (PDB ID: 5cxf) and strain Towne (PDB IDs: 4osn, 5c6t), gB EBV strain B95-8 (PDB ID: 3fvc), and gB HSV-1 strain KOS (PDB IDs: 2gum and 3nwa). Using the structure of most closely related HCMV gB, as a model, shows Domains I (blue), II (green), III (yellow), and IV (orange) in the post-fusion conformation. Domain I (blue) shows conserved beta sheets connected by the two fusion loops, FL1 and FL2, at the base of the model (purple). (c) Longitudinal view of the predicted trimer of HHV-6A gB showing the Thr193 position in red within three subunits, indicated in purple, blue, and orange. (d) Domain I (residues 79–291) is enlarged as indicated by the grey T bar, showing the position of Thr193, in red, in the beta sheet at the trimer interfaces with the fusion loops FL1 (residues 100–105) and FL2 (residues 184–191) exposed at the base of the structure.
Viruses 10 00016 g003
Table
Table 1. Long range PCR genome amplicon primers.
Table 1. Long range PCR genome amplicon primers.
AmpliconForward Primer (5′-3′)Reverse Primer (5′-3′)Size (bp)
DR1ACAACCGCATCTTCTTCAGATGAACGGAGATCTGGGAG1750
DR1DR6CCTCATCTGTTATTCCCTTCCACTCACCGCAGTCTTGAAC2872
DR6TCCAAACTGTACGGCTATGTGGCGTCTAATCCAGGTAC1858
DR_ULCATGAGAGGACACTGGACGTACTCGACATCTGTGGACCATGCA2339
0203ACCATGTGCGGAGAAAGTTGACGAGAATGAATATGCCGATG2174
0306ACTCCTATCCGTCTAAACTGCTATTTCTCACGCCGGTATTC4417
0611ATGCCAAATCTAGCTGCTGGAATCGGAGGAAGCTAACGT4757
1114TCTTCAGGTGTCTCAACCTTCTGATGATCTGCTGTTGGAAC5737
1419TTGAGATGCAGAGCCTCTCACGTGGTGGATAATTCTTCAG5370
1923CTACATCATCGGTTGCCATATACAGAAGTTTCGGCAACTG5216
2328TAGTCGCTGGATGTCGATACGAGGAATTGTGGACAACGCT4721
2830CATCCAACTCCTTGTCCGTTCGTTCGACGTGTCTGAC5753
3031GACACTGTCCAACGCATCTCTCTTGGGCATAAGTCCAGTA6552
3134ACAGAGCGACACCCAATATTCAGCTGTGTCGATGTCTTGCTA2787
3439TGCATTCTCTCACAGCATGTGTTTCCTTATGCCACCAATC6309
3941TGAGGCATGAGACTCAGAAGAGCGGCATCCATAATGCTTA4796
4142AGCATCGTTTGAAGACACCGAGTTGTTGTTCGCTGGTCT2810
4142INTTTGTTTCAGGCCAGCATCATACCGAAAACTCAGGCAAC2983
4245AGACCGCATTCGTATCGCTCACGATCAAGTTCGAATCGTA5092
4549GAGTGCTCGGATAAGTTCATTACGATACAATCCGGAATGC5922
4953TGGAACTCGGAATTGTTTCTGTTTGAGACTGTATCTGCGCA4282
5357TATGGTGAAGAACCGTCGTCCGCAGATAACAGAATGG4707
5758CAGTGCCATCTGTTATGAATAGAGTCGCAGTAAGGTCCACGT6792
5864TGCTGAAATTACCTCAGTGAGCCATTCAGTATAATGCATCC4229
6469CTGTATCAATATCAAGGCGGACCACTGAATACGAACGCTC5111
6974GCGATTGTTCTGCGATAGAGCTTATATCCGAGCCTTGCAGT6472
7476GAGAACATGCTAGACAATTGGCAAGAACTGCGACTCAATCC2668
7679GTTGCTAGTTGTATGACTTGGCAATATCCACCGTTAGAGAAC7681
7984AGATGTATGCTGAAGAACGTGGGATCGTCAACCGTTAGTG5109
8486CACAGTGTGTTCGCCGGAAGACAGACAATGCACATCCTCTG4773
8690AGGTTGATACGGCAACGAGGCAATCATTAGCATACAGATG3326
8690INTTGAGGAATCACGTGTTTGATCGAATCTATCCATGAAGATG1746
90R3CATCATTGTTATCGCTTTCACGCAACCGCAGTTCCTGTT4912
R395GCGGTACCCACTGATCTTAGTCTACCAGGCATTCCGT6250
95100GAGGAGGGTCTGTCTAGATGTTCGGAGATATCATAATCTGCGT4031
100DRTTATAGTTGCTCCCGAAAGCAAGAAGATGCGGTTGTCTTG2358
Table
Table 2. Comparison of Illumina resequence human betaherpesvirus 6 (HHV-6A) U1102 to Reference sequence NC_001664.
Table 2. Comparison of Illumina resequence human betaherpesvirus 6 (HHV-6A) U1102 to Reference sequence NC_001664.
ORF or Locus Refseq bpRefseq:SNPNon-Synonymous SNP Amino Acid Change, Codons and (Reverse Complement Codons)This Study or as Cited
DR-L DR66663G:GG
“6665G:GGFrameshift extension[16]
DR-L repeat7623G:C““
Repeats left end U8262AGGAG:A
“8362C:CGCG
“8518T:C
“8559ACAACA:CAC
“8591AT:A
“8602GG:G
U1018884C:Indel1Frameshift extension[34,35]
U4063871G:TThr:His (ACG:CAG)
“63872T:G“
U4270446A:ATCHisCys:GlnLeuAsp (CACTGT:CAACTGCATG)Jones & Teo 1993 in Genbank L20954
“70449G:GT“Jones, M. in Genbank X92436
U57 mcp92149C:GLys:Asn ( AAG:AAC)
U58 tbp94068A:GMet:Val ATG:GTG
U83 chemokine123740C:GAspGlu:GluGln GACGAG:GAGCAG[15]
“123741G:C““
U86 IE2(128132)(Indel G)Frameshift U86/U87In RefSeq; [17]
R1 repeats128325C:T
“128443C:T
“128914G:A
“129917A:ATCC
R2 repeats132312GG:G
R3 repeats137831C:G
“137832G:C
“137906C:G
“137907G:C
U100148154C:G
“148155G:C
U100 gQ1149553C:GArg:Ala (CGT: GCT)
“149554G:C“
DR-R DR6157898G:GGFrameshift extension[16]
“157900G:GG“
DR-R repeat158857G:C
ORF = open reading frame; SNP = single nucleotide polymorphisms; Repetitive regions are shaded; Indel1: CACACCTATCGCGACCCTGTTTTATGTAGAAAGAATGACAAGCCAGAACAATGA; (“) indicates same as above.
Table
Table 3. Minor variant polymorphism populations in HHV-6A U1102, passaged, and BACmid virus.
Table 3. Minor variant polymorphism populations in HHV-6A U1102, passaged, and BACmid virus.
ORF or LociRefSeq
bp		RefSeq
:SNP		Amino Acid Change, Codons and (Reverse Complement Codons)% New Reference% BAC Virus% Passaged VirusComment, Citations, and Other Strains with SNPs
DR8016A:C -87%-End T2
“8077C:CGA 96%--HHV-6A AJ KP257,584.1, GS KC465,951.1; CiHHV6A 2284 KT895,199.1, 1501 KT355,575.1; [37]; Sites cleavages after pac2; [26,31]
“8080G:GAC -86%-CiHHV6A 2284 KT895,199.1, 5055, 5814 [16,24]; “
Uleft repeat8157G:GC --78%
“8558A:AC 88%AAC:A 99%90%AC repeat
“8561AAC:A -100%98%HHV-6A AJ [15,31]
“8601CG:C 97%87%86%
U28876G:TAla314:Glu(GCG:GAG)--53%
U19 UTR28,371C:CA -52%72%HHV-6A AJ
U31 tegument45,667A:TAsn173:Ile AAC:ATC-73%96%
U33 capsid52,148A:GTyr330:His(TAT:CAT)-54%69%
U39 gB*61,162T:CThr193:Ala(ACG:GCG)-58%99%
U41 dbp64,707C:TAla972:Thr (GCT:ACT)--90%HHV-6A GS KJ123,690.1
ORI repeat67,674T:C --54%
“67,860T:C --34%HHV-6A GS KJ123,690.1 & KC465,951.1
U50 capsid81,583G:AAla258:Thr(GCA:ACA)-59%78%HHV-6A GS KC465,951.1
“82,158G:C -67%-
U73 obp109,214G:ASer297:Asn(AGT:AAT)--99%
U83 chemokine123,504G:GTTf/s N-terminal extension <1%--cDNA clone [13,14,15]; HHV-6A AJ KP257584;
CiHHV6A: 5055 AIX09949, 23,090 AIX09956, 1501 KT355,575
U86 IE2126,126T:CAsn395:Asp(AAT:GAT)--99%
“127,092C:CTGA Ser976:SerSer
(TCA repeat ×9:×10)		-45%-2 Ser insertion CiHHV6A 5055; 5814; [24]
R2 repeat131,092C:T --56%IE2 intron
“131,743G:C --52%“
“131,745A:C --62%“
R2 region 132,303TG:T -64%56%“ TGG:T HHV-6A AJ [31]
R3 repeat140,872G:A -96%-End 30× repeats from 137,775; latency IE1/2 intron [38]
U95142,969T:CTrp10:Arg TGG:CGG--54%
“142,979T:AIle13:Asn ATT:AAT--53%
“143,591C:ASer217:Tyr TCC:TAC--76%
U100 gQ2146,729C:TC-terminal truncation; Trp186:stop (TGG:TGA)<1%--[12]1 (in 7/8 gQ2 cDNA clones HHV-6A U1102)
“146,862G:ASer42:Leu (TCG:TTG) -60%99%Not in HHV-6A HHV-6A U1102 from [39]
U100 gQ1149,251Indel2:Aintron-44%83%HHV-6A AJ; CiHHV6A 2284, 5055, 5814 [24]
Indel2: AGGGGGGGG. Repetitive regions are shaded, DR repeats listed once here. * Adjacent to fusion loop 2 of Domain I (DI) homologue of human cytomegalovirus (HCMV) gB C246 [40], see text. (“) indicates same as above. (–) indicates no change in comparison to the reference genome sequence.

Share and Cite

MDPI and ACS Style

Tweedy, J.G.; Escriva, E.; Topf, M.; Gompels, U.A. Analyses of Tissue Culture Adaptation of Human Herpesvirus-6A by Whole Genome Deep Sequencing Redefines the Reference Sequence and Identifies Virus Entry Complex Changes. Viruses 2018, 10, 16. https://0-doi-org.brum.beds.ac.uk/10.3390/v10010016

AMA Style

Tweedy JG, Escriva E, Topf M, Gompels UA. Analyses of Tissue Culture Adaptation of Human Herpesvirus-6A by Whole Genome Deep Sequencing Redefines the Reference Sequence and Identifies Virus Entry Complex Changes. Viruses. 2018; 10(1):16. https://0-doi-org.brum.beds.ac.uk/10.3390/v10010016

Chicago/Turabian Style

Tweedy, Joshua G., Eric Escriva, Maya Topf, and Ursula A. Gompels. 2018. "Analyses of Tissue Culture Adaptation of Human Herpesvirus-6A by Whole Genome Deep Sequencing Redefines the Reference Sequence and Identifies Virus Entry Complex Changes" Viruses 10, no. 1: 16. https://0-doi-org.brum.beds.ac.uk/10.3390/v10010016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop