Human Herpesvirus 6 Genome Research Articles

Molecular epidemiology of human herpesvirus 8 in patients with HHV-8-related diseases in Ireland.

Human Herpesvirus 8 (HHV-8) has been classified by sequence analysis of open reading frame (ORF) K1, ORF K15, and variable sequence loci within the central constant region. The purpose of this study was to examine the molecular epidemiology of HHV-8 in an Irish population. This retrospective study included 30 patients who had HHV-8 DNA detected in plasma. Nested end-point PCR was used to characterise four regions of the HHV-8 genome, K1, T0.7 (K12), ORF 75, and K15. Sequencing data were obtained for 23 specimens from 19 patients. Phylogenetic analysis of ORF K1 demonstrated that subtypes A, B, C and F were present in 37%, 11%, 47% and 5%, respectively. For T0.7 and ORF 75, sequencing data were obtained for 12 patients. For T0.7, subtypes A/C, J, B, R and Q were present in 58%, 17%, 8%, 8%, and 8%, respectively. For ORF 75, subtypes A, B, C and D were present in 58%, 8%, 25%, and 8%, respectively. K15 sequences were determined for 13 patients. 69% had the P allele and 31% had the M allele. The data generated by this study demonstrate that a broad variety of HHV-8 subtypes are represented in patients exhibiting HHV-8-related disease in Ireland, a low prevalence country. The predominance of C and A K1 subtypes was as expected for a Western European population. The 31% prevalence for K15 subtype M was higher than expected for a Western European population. This may represent the changing and evolving epidemiology in Ireland due to altered migration patterns.

Excision of Integrated Human Herpesvirus 6A Genomes Using CRISPR/Cas9 Technology.

Human herpesviruses 6A and 6B are betaherpesviruses that can integrate their genomes into the telomeres of latently infected cells. Integration can also occur in germ cells, resulting in individuals who harbor the integrated virus in every cell of their body and can pass it on to their offspring. This condition is termed inherited chromosomally integrated HHV-6 (iciHHV-6) and affects about 1% of the human population. The integrated HHV-6A/B genome can reactivate in iciHHV-6 patients and in rare cases can also cause severe diseases including encephalitis and graft-versus-host disease. Until now, it has remained impossible to prevent virus reactivation or remove the integrated virus genome. Therefore, we developed a system that allows the removal of HHV-6A from the host telomeres using the CRISPR/Cas9 system. We used specific guide RNAs (gRNAs) targeting the direct repeat region at the ends of the viral genome to remove the virus from latently infected cells generated in vitro and iciHHV-6A patient cells. Fluorescence-activated cell sorting (FACS), quantitative PCR (qPCR), and fluorescence in situ hybridization (FISH) analyses revealed that the virus genome was efficiently excised and lost in most cells. Efficient excision was achieved with both constitutive and transient expression of Cas9. In addition, reverse transcription-qPCR (RT-qPCR) revealed that the virus genome did not reactivate upon excision. Taken together, our data show that our CRISPR/Cas9 approach allows efficient removal of the integrated virus genome from host telomeres. IMPORTANCE Human herpesvirus 6 (HHV-6) infects almost all humans and integrates into the telomeres of latently infected cells to persist in the host for life. In addition, HHV-6 can also integrate into the telomeres of germ cells, which results in about 80 million individuals worldwide who carry the virus in every cell of their body and can pass it on to their offspring. In this study, we develop the first system that allows excision of the integrated HHV-6 genome from host telomeres using CRISPR/Cas9 technology. Our data revealed that the integrated HHV-6 genome can be efficiently removed from the telomeres of latently infected cells and cells of patients harboring the virus in their germ line. Virus removal could be achieved with both stable and transient Cas9 expression, without inducing viral reactivation.

A Panel of Kaposi's Sarcoma-Associated Herpesvirus Mutants in the Polycistronic Kaposin Locus for Precise Analysis of Individual Protein Products.

ABSTRACTKaposi’s sarcoma-associated herpesvirus (KSHV) is the cause of several human cancers, including the endothelial cell (EC) malignancy, Kaposi’s sarcoma. Unique KSHV genes absent from other human herpesvirus genomes, the “K-genes,” are important for KSHV replication and pathogenesis. Among these, the kaposin transcript is highly expressed in all phases of infection, but its complex polycistronic nature has hindered functional analysis to date. At least three proteins are produced from the kaposin transcript: Kaposin A (KapA), B (KapB), and C (KapC). To determine the relative contributions of kaposin proteins during KSHV infection, we created a collection of mutant viruses unable to produce kaposin proteins individually or in combination. In previous work, we showed KapB alone recapitulated the elevated proinflammatory cytokine transcripts associated with KS via the disassembly of RNA granules called processing bodies (PBs). Using the new ΔKapB virus, we showed that KapB was necessary for this effect during latent KSHV infection. Moreover, we observed that despite the ability of all kaposin-deficient latent iSLK cell lines to produce virions, all displayed low viral episome copy number, a defect that became more pronounced after primary infection of naive ECs. For ΔKapB, provision of KapB in trans failed to complement the defect, suggesting a requirement for the kaposin locus in cis. These findings demonstrate that our panel of kaposin-deficient viruses enables precise analysis of the respective contributions of individual kaposin proteins to KSHV replication. Moreover, our mutagenesis approach serves as a guide for the functional analysis of other complex multicistronic viral loci.IMPORTANCE Kaposi’s sarcoma-associated herpesvirus (KSHV) expresses high levels of the kaposin transcript during both latent and lytic phases of replication. Due to its repetitive, GC-rich nature and polycistronic coding capacity, until now no reagents existed to permit a methodical analysis of the role of individual kaposin proteins in KSHV replication. We report the creation of a panel of recombinant viruses and matched producer cell lines that delete kaposin proteins individually or in combination. We demonstrate the utility of this panel by confirming the requirement of one kaposin translation product to a key KSHV latency phenotype. This study describes a new panel of molecular tools for the KSHV field to enable precise analysis of the roles of individual kaposin proteins during KSHV infection.

Metagenomic profiling of placental tissue suggests DNA virus infection of the placenta is rare.

It is widely recognized that pathogens can be transmitted across the placenta from mother to foetus. Recent re-evaluation of metagenomic studies indicates that the placenta has no unique microbiome of commensal bacteria. However, viral transmission across the placenta, including transmission of DNA viruses such as the human herpesviruses, is possible. A fuller understanding of which DNA virus sequence can be found in the placenta is required. We employed a metagenomic analysis to identify viral DNA sequences in placental metagenomes from full-term births (20 births), pre-term births (13 births), births from pregnancies associated with antenatal infections (12 births) or pre-term births with antenatal infections (three births). Our analysis found only a small number of DNA sequences corresponding to the genomes of human herpesviruses in four of the 48 metagenomes analysed. Therefore, our data suggest that DNA virus infection of the placenta is rare and support the concept that the placenta is largely free of pathogen infection.

Genome-wide analysis of codon usage pattern in herpesviruses and its relation to evolution

The preferential use of a specific codon, out of a group of synonymous codons encoding the same amino acid, in a gene transcript results from the bias in codon choice. Various evolutionary forces namely mutation pressure and natural selection influence the pattern of codon usage i.e. distinct for each gene/genome. We investigated the pattern of codon usage of eight human herpesvirus genomes and compared them with two other herpesvirus genomes namely murine herpesvirus 68 and bovine herpesvirus type 1.1 to elucidate its compositional features, pattern of codon usage across the genomes and report the differences of codon usage pattern of human herpesviruses from that of other two other viruses. We also identified the similarity of the codon usage of human herpesviruses with its host (human). The genes were found to be CG rich in HHV2, HHV3, HHV4, HHV6, HHV7 and BH genomes while TA rich in HHV1, HHV5, HHV8 and MH genomes. The codon usage bias (CUB) of genes was low. A highly significant correlation was found among compositional contents depicting the role of mutational pressure along with natural selection in framing CUB. Several more frequently used codons as well as less frequently used codons were identified to be similar between each human virus and its host (human), while murine herpesvirus 68 and bovine herpesvirus type 1.1 genomes did not possess similar adaptation strategy as human herpesviruses to human (host), thus we could conclude that viral CUB might have been shaped as per their host’s nature for better surveillance. Neutrality plot revealed mutational pressure mostly influenced the CUB of HHV1, HHV8 and MH viruses, while natural selection had a major impact in the CUB of HHV2, HHV3, HHV4, HHV5, HHV6, HHV7 and BH genomes.

Inherited chromosomally integrated human herpesvirus 6 and autoimmune connective tissue diseases.

Entire genome of human herpesvirus 6 (HHV-6) that integrates into human chromosomes is called chromosomally integrated HHV-6 (ciHHV-6). Several viral infections have been suggested to be involved in autoimmune connective tissue diseases (CTDs). Reactivated HHV-6 from the integrated viral genome can induce immune responses against the virus. Thus, it is plausible that ciHHV-6 is associated with autoimmune CTDs. We sought to determine whether the prevalence of ciHHV-6 was significantly higher in patients with autoimmune CTDs than in a healthy population. A total of 846 peripheral blood samples collected from autoimmune CTD patients were analyzed. Since there was a large number of samples, they were pooled into 24 samples per group. Copy numbers of HHV-6 DNA were measured by real-time PCR. The threshold level for distinguishing between ciHHV-6 and active viral infection and the reliability of pooled DNA analysis were examined as initial validation experiments. The threshold level was 1.6 × 10^6 copy/mL in whole blood. The reliability of pooled DNA analysis to identify one ciHHV-6 sample among 23 HHV-6 DNA-negative samples was high. No HHV-6 DNA was detected in any of the pooled DNA samples collected from the patients. The probability of the present study including the 846 autoimmune CTD patient's samples was statistically not different with a healthy Japanese population which was 0.2 % or 0.6 %. There was no significant difference in the prevalence of ciHHV-6 between a healthy population and patients with autoimmune CTDs.

A survey of RNA secondary structural propensity encoded within human herpesvirus genomes: global comparisons and local motifs.

There are nine herpesviruses known to infect humans, of which Epstein–Barr virus (EBV) is the most widely distributed (>90% of adults infected). This ubiquitous virus is implicated in a variety of cancers and autoimmune diseases. Previous analyses of the EBV genome revealed numerous regions with evidence of generating unusually stable and conserved RNA secondary structures and led to the discovery of a novel class of EBV non-coding (nc)RNAs: the stable intronic sequence (sis)RNAs. To gain a better understanding of the roles of RNA structure in EBV biology and pathogenicity, we revisit EBV using recently developed tools for genome-wide motif discovery and RNA structural characterization. This corroborated previous results and revealed novel motifs with potential functionality; one of which has been experimentally validated. Additionally, since many herpesviruses increasingly rival the seroprevalence of EBV (VZV, HHV-6 and HHV-7 being the most notable), analyses were expanded to include all sequenced human Herpesvirus RefSeq genomes, allowing for genomic comparisons. In total 10 genomes were analyzed, for EBV (types 1 and 2), HCMV, HHV-6A, HHV-6B, HHV-7, HSV-1, HSV-2, KSHV, and VZV. All resulting data were archived in the RNAStructuromeDB (https://structurome.bb.iastate.edu/herpesvirus) to make them available to a wide array of researchers.

Endogenization and excision of human herpesvirus 6 in human genomes.

Sequences homologous to human herpesvirus 6 (HHV-6) are integrated within the nuclear genome of about 1% of humans, but it is not clear how this came about. It is also uncertain whether integrated HHV-6 can reactivate into an infectious virus. HHV-6 integrates into telomeres, and this has recently been associated with polymorphisms affecting MOV10L1. MOV10L1 is located on the subtelomere of chromosome 22q (chr22q) and is required to make PIWI-interacting RNAs (piRNAs). As piRNAs block germline integration of transposons, piRNA-mediated repression of HHV-6 integration has been proposed to explain this association. In vitro, recombination of the HHV-6 genome along its terminal direct repeats (DRs) leads to excision from the telomere and viral reactivation, but the expected “solo-DR scar” has not been described in vivo. Here we screened for integrated HHV-6 in 7,485 Japanese subjects using whole-genome sequencing (WGS). Integrated HHV-6 was associated with polymorphisms on chr22q. However, in contrast to prior work, we find that the reported MOV10L1 polymorphism is physically linked to an ancient endogenous HHV-6A variant integrated into the telomere of chr22q in East Asians. Unexpectedly, an HHV-6B variant has also endogenized in chr22q; two endogenous HHV-6 variants at this locus thus account for 72% of all integrated HHV-6 in Japan. We also report human genomes carrying only one portion of the HHV-6B genome, a solo-DR, supporting in vivo excision and possible viral reactivation. Together these results explain the recently-reported association between integrated HHV-6 and MOV10L1/piRNAs, suggest potential exaptation of HHV-6 in its coevolution with human chr22q, and clarify the evolution and risk of reactivation of the only intact (non-retro)viral genome known to be present in human germlines.

Revisiting the genomes of herpesviruses.

Combining integrative genomics and systems biology approaches has revealed new and conserved features in the genome of human herpesvirus 6.

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.

Human herpesvirus 8 DNA detection and variant analysis in patients with multiple sclerosis

Several studies reported a complex interplay between viral infections and neural cells leading to multiple sclerosis. A role for some viral infections has been proposed in MS. In this study, DNA sequences of human herpesvirus 8 (HHV-8) were searched in the peripheral blood of 54 patients with multiple sclerosis and 130 healthy subjects using nested-PCR assay to amplify ORF26 locus. Furthermore, HHV-8 positive samples were subjected to a nested-PCR to amplify K1 gene of HHV-8 followed by direct nucleotide sequencing. HHV-8 genome was detected in 18.5% (10/54) and 3% (4/130) of MS patients and controls, respectively, and the difference reached statistically significant level (P = 0.0017). Genotyping analysis revealed that genotype C was common (88.9%) in all study subjects, followed by genotype A. Our results showed higher detection of HHV-8 DNA in MS patients than control group.

High Prevalence of Distinct Human Herpesvirus 8 Contributes to the High Incidence of Non-acquired Immune Deficiency Syndrome-Associated Kaposi's Sarcoma in Isolated Japanese Islands.

Non-acquired immune deficiency syndrome (AIDS) Kaposi's sarcoma (KS) is extremely rare in Japan but highly endemic in Okinawa, especially in Miyako Islands. We aimed to elucidate the exact incidence and cause of this high prevalence. Non-AIDS KS cases in Okinawa Prefecture over the past 31 years were reviewed, and human herpesvirus 8 (HHV8) seroprevalence in Miyako Islands was determined. We examined whole-genome sequences of 3 HHV8 strains and performed whole-exome sequencing of 4 male patients from Miyako Islands. Approximately half of the non-AIDS KS cases in Okinawa Prefecture were from Miyako Islands. The age-adjusted incidence rate was 0.87/105 per year for Miyako Islands and 0.056/105 per year for the rest of Okinawa. Human herpesvirus 8 seroprevalence was 15.4% in Miyako Islands. The 3 HHV8 genomes isolated from Miyako islanders formed a phylogenetically branch distinct from those of previously sequenced HHV8 strains and shared specific mutations in 9 proteins. These mutations were verified in Okinawan patients other than those from Miyako Islands. Whole-exome sequencing of the 4 male Miyako Islanders did not reveal shared pathogenic mutations. Miyako Islands are an endemic area of non-AIDS KS. The high rate of a distinct HHV8 may contribute to the high incidence of KS in the region.

K-mer Content, Correlation, and Position Analysis of Genome DNA Sequences for the Identification of Function and Evolutionary Features.

In genome analysis, k-mer-based comparison methods have become standard tools. However, even though they are able to deliver reliable results, other algorithms seem to work better in some cases. To improve k-mer-based DNA sequence analysis and comparison, we successfully checked whether adding positional resolution is beneficial for finding and/or comparing interesting organizational structures. A simple but efficient algorithm for extracting and saving local k-mer spectra (frequency distribution of k-mers) was developed and used. The results were analyzed by including positional information based on visualizations as genomic maps and by applying basic vector correlation methods. This analysis was concentrated on small word lengths (1 ≤ k ≤ 4) on relatively small viral genomes of Papillomaviridae and Herpesviridae, while also checking its usability for larger sequences, namely human chromosome 2 and the homologous chromosomes (2A, 2B) of a chimpanzee. Using this alignment-free analysis, several regions with specific characteristics in Papillomaviridae and Herpesviridae formerly identified by independent, mostly alignment-based methods, were confirmed. Correlations between the k-mer content and several genes in these genomes have been found, showing similarities between classified and unclassified viruses, which may be potentially useful for further taxonomic research. Furthermore, unknown k-mer correlations in the genomes of Human Herpesviruses (HHVs), which are probably of major biological function, are found and described. Using the chromosomes of a chimpanzee and human that are currently known, identities between the species on every analyzed chromosome were reproduced. This demonstrates the feasibility of our approach for large data sets of complex genomes. Based on these results, we suggest k-mer analysis with positional resolution as a method for closing a gap between the effectiveness of alignment-based methods (like NCBI BLAST) and the high pace of standard k-mer analysis.

Possible chromosomal and germline integration of human herpesvirus 7.

Human herpesvirus 7 (HHV-7) is a betaherpesvirus, and is phylogenetically related to both HHV-6A and HHV-6B. The presence of telomeric repeat sequences at both ends of its genome should make it equally likely to integrate into the human telomere as HHV-6. However, numerous studies have failed to detect germline integration of HHV-7, suggesting an important difference between the HHV-6A/-6B and HHV-7 genomes. In search of possible germline integrated HHV-7, we developed a sensitive and quantitative real-time PCR assay and discovered that primers designed against some parts of the HHV-7 genome can frequently miss HHV-7 positive clinical samples even though they work efficiently in cell-culture-derived HHV-7 positive materials. Using a primer pair against the U90 ORF of HHV-7, we identified a possible case of germline integration of HHV-7 with one copy of viral genome per cell in both peripheral blood cells and hair follicles. Chromosomal integration of HHV-7 in these individuals was confirmed by fluorescence in situ hybridization analysis. Germline integration of HHV-7 was further confirmed by detection of ~2.6 copies of HHV-7 in the hair follicles of one of the parents. Our results shed light on the complex nature of the HHV-7 genome in human-derived materials in comparison to cell-culture-derived materials and show the need for stringent criteria in the selection of primers for epidemiological HHV-7 studies.

Prevalence of human herpesvirus-8 among HIV-infected patients, intravenous drug users and the general population in Iran.

Studies looking at the frequency of human herpesvirus-8 (HHV-8) among Iranian blood donors have produced conflicting results. The aim of this study was to investigate the prevalence of HHV-8 DNA by using polymerase chain reaction methods among 168 healthy individuals, 60 intravenous drug users and 100 HIV-infected patients from Iran. The prevalence of HHV-8 was significantly higher among intravenous drug users (13.3%) compared with the general population (3.6%; P=0.017). The HHV-8 genome was mostly detected among intravenous drug users who displayed high-risk sexual behaviours. Moreover, the HHV-8 genome was also detected in 8% of HIV-infected patients. The present study findings support the likelihood that the transmission of HHV-8 is via a sexual route in the Iranian population.

Comparative Genomics of Herpesviridae Family to Look for Potential Signatures of Human Infecting Strains

Herpesviridae family is one of the significant viral families which comprises major pathogens of a wide range of hosts. This family includes at least eight species of viruses which are known to infect humans. This family has evolved 180–220 million years ago and the present study highlights that it is still evolving and more genes can be added to the repertoire of this family. In addition, its core-genome includes important viral proteins including glycoprotein B and helicase. Most of the infections caused by human herpesviruses have no definitive cure; thus, search for new therapeutic strategies is necessary. The present study finds core-genome of human herpesviruses that differs from that of Herpesviridae family and nonhuman herpes strains of this family and might be a putative target for vaccine development. The phylogenetic reconstruction based upon the protein sequences of core gene set of Herpesviridae family reveals the sharp splits of its different subfamilies and supports the hypothesis of coevolution of viruses with their hosts. In addition, data mining for cis-elements in the genomes of human herpesviruses results in the prediction of numerous regulatory elements which can be used for regulating the expression of viral based vectors implicated in gene therapies.

Association of Mycosis Fungoides and Large Plaque Parapsoriasis with Human Herpes Virus 8

Background: The association of the human herpes virus 8 (HHV-8) with Mycosis fungoides (MF) and large plaque parapsoriasis (LPP) has been assessed previously and contradictory results were reported. Although it was investigated in different countries and different ethnic groups, there was no study concerning Iranian population with these lymphoproliferative diseases. Objectives: We aimed to assess the association of HHV-8 with MF and LPP in biopsy samples of Iranian patients with these lymphoproliferative conditions. Patients and Methods: Paraffin-embedded samples from patients who were diagnosed by histopathological examination as MF or LPP were assessed by real-time polymerase chain reaction (RT-PCR) in a molecular epidemiology case-control study. Primerdesign Genesig kit for HHV-8 Genomes was used to extract total DNA. Samples from chronic dermatitis were used as control; five samples from patients with kaposi sarcoma were used as positive control in PCR. Results: The investigated samples consisted of 13 patient with MF, 10 patients with LPP, and 23 patients with chronic dermatitis; the diagnosis of these skin conditions were confirmed histopathologically. The RT-PCR could not detect HHV-8 DNA sequence in any of the examined samples, ie. MF, LPP, and chronic dermatitis. The PCR results for all kaposi sarcoma samples were positive for HHV-8 genome. Conclusion: There was no association between HHV-8 and MF or LPP in Iranian population. This is in accordance with most of the reports from different countries.

Presence of HHV-6 genome in spermatozoa in a context of couples with low fertility: what type of infection?

Human herpesvirus-6 (HHV-6) is a betaherpesvirus whose genome may integrate into human chromosomes. Chromosomally integrated HHV-6 (ciHHV-6) may be transmitted vertically from parents to children. HHV-6 DNA has been detected in semen, but its integrated or extrachromosomal status has not yet been characterised. The aim of this study was to determine the prevalence of HHV-6 DNA and to search for ciHHV-6 forms in spermatozoa purified from semen obtained from subjects explored for low fertility. A total of 184 sperm samples were purified using PureSperm(®) . HHV-6 viral load and species identification were performed by real-time polymerase chain reaction. Of 179 sperm specimens analysed, three were positive for HHV-6 (1.7%). Two samples (1.1%) had viral loads of 680232 and 2834075copies per million spermatozoa, compatible with loads expected for a ciHHV-6 form. The viral load of the third positive sample (73684copies per million spermatozoa) was lower than would be expected for ciHHV-6 infection, implying that the HHV-6 DNA detected in spermatozoa corresponds mainly to ciHHV-6. However, viral DNA may also be detected at a low level that is not in favour of the presence of ciHHV-6. Further studies are necessary to determine the origin of detected viral genomes.

Reactivation of Chromosomally Integrated Human Herpesvirus-6 by Telomeric Circle Formation

More than 95% of the human population is infected with human herpesvirus-6 (HHV-6) during early childhood and maintains latent HHV-6 genomes either in an extra-chromosomal form or as a chromosomally integrated HHV-6 (ciHHV-6). In addition, approximately 1% of humans are born with an inheritable form of ciHHV-6 integrated into the telomeres of chromosomes. Immunosuppression and stress conditions can reactivate latent HHV-6 replication, which is associated with clinical complications and even death. We have previously shown that Chlamydia trachomatis infection reactivates ciHHV-6 and induces the formation of extra-chromosomal viral DNA in ciHHV-6 cells. Here, we propose a model and provide experimental evidence for the mechanism of ciHHV-6 reactivation. Infection with Chlamydia induced a transient shortening of telomeric ends, which subsequently led to increased telomeric circle (t-circle) formation and incomplete reconstitution of circular viral genomes containing single viral direct repeat (DR). Correspondingly, short t-circles containing parts of the HHV-6 DR were detected in cells from individuals with genetically inherited ciHHV-6. Furthermore, telomere shortening induced in the absence of Chlamydia infection also caused circularization of ciHHV-6, supporting a t-circle based mechanism for ciHHV-6 reactivation.

Genome-wide prediction of vaccine targets for human herpes simplex viruses using Vaxign reverse vaccinology

Herpes simplex virus (HSV) types 1 and 2 (HSV-1 and HSV-2) are the most common infectious agents of humans. No safe and effective HSV vaccines have been licensed. Reverse vaccinology is an emerging and revolutionary vaccine development strategy that starts with the prediction of vaccine targets by informatics analysis of genome sequences. Vaxign (http://www.violinet.org/vaxign) is the first web-based vaccine design program based on reverse vaccinology. In this study, we used Vaxign to analyze 52 herpesvirus genomes, including 3 HSV-1 genomes, one HSV-2 genome, 8 other human herpesvirus genomes, and 40 non-human herpesvirus genomes. The HSV-1 strain 17 genome that contains 77 proteins was used as the seed genome. These 77 proteins are conserved in two other HSV-1 strains (strain F and strain H129). Two envelope glycoproteins gJ and gG do not have orthologs in HSV-2 or 8 other human herpesviruses. Seven HSV-1 proteins (including gJ and gG) do not have orthologs in all 40 non-human herpesviruses. Nineteen proteins are conserved in all human herpesviruses, including capsid scaffold protein UL26.5 (NP_044628.1). As the only HSV-1 protein predicted to be an adhesin, UL26.5 is a promising vaccine target. The MHC Class I and II epitopes were predicted by the Vaxign Vaxitop prediction program and IEDB prediction programs recently installed and incorporated in Vaxign. Our comparative analysis found that the two programs identified largely the same top epitopes but also some positive results predicted from one program might not be positive from another program. Overall, our Vaxign computational prediction provides many promising candidates for rational HSV vaccine development. The method is generic and can also be used to predict other viral vaccine targets.