A full-genome phylogenetic analysis of varicella-zoster virus reveals a novel origin of replication-based genotyping scheme and evidence of recombination between major circulating clades

Genomic sequencing revealed recombination event between clade 1 and clade 2 occurs in circulating varicella-zoster virus in China

Varicella-zoster virus (VZV), a member of the Alphaherpesvirinae subfamily, causes varicella in primary infections and establishing a latent stage in sensory ganglia. Upon reactivation, VZV causes herpes zoster with severe neuralgia, especially in elderly patients. The mutation rate for VZV is comparatively lower than the other members of other alpha herpesviruses. Due to geographic isolation, different genotypes of VZV are circulating on separate continents. Here, we successfully isolated a VZV from the vesicular fluid of a youth zoster patient. Based on the single-nucleotide polymorphism profiles of different open reading frames that define the genotype, this newly isolated VZV primarily represents genotype clade 2 but also has characteristics of genotype clade 1. The next-generation sequencing provided a nearly full-length sequence, and further phylogenetic analysis revealed that this VZV isolate is distinct from clades 1 and 2. The Recombination Detection Program indicates that a possible recombinant event may occur between the VZV isolate and clade 1. In summary, we found that there is a circulating VZV isolate in China that may represent a recombinant between clade 1 and clade 2, providing new concerns that need to be considered in the future VZV vaccination program.

Human herpesvirus diversity is altered in HLA class I binding peptides

Herpesviruses are ubiquitous, genetically diverse DNA viruses, with long-term presence in humans associated with infrequent but significant pathology. Human leukocyte antigen (HLA) class I presents intracellularly derived peptide fragments from infected tissue cells to CD8+ T and natural killer cells, thereby directing antiviral immunity. Allotypes of highly polymorphic HLA class I are distinguished by their peptide binding repertoires. Because this HLA class I variation is a major determinant of herpesvirus disease, we examined if sequence diversity of virus proteins reflects evasion of HLA presentation. Using population genomic data from Epstein–Barr virus (EBV), human cytomegalovirus (HCMV), and Varicella–Zoster virus, we tested whether diversity differed between the regions of herpesvirus proteins that can be recognized, or not, by HLA class I. Herpesviruses exhibit lytic and latent infection stages, with the latter better enabling immune evasion. Whereas HLA binding peptides of lytic proteins are conserved, we found that EBV and HCMV proteins expressed during latency have increased peptide sequence diversity. Similarly, latent, but not lytic, herpesvirus proteins have greater population structure in HLA binding than nonbinding peptides. Finally, we found patterns consistent with EBV adaption to the local HLA environment, with less efficient recognition of EBV isolates by high-frequency HLA class I allotypes. Here, the frequency of CD8+ T cell epitopes inversely correlated with the frequency of HLA class I recognition. Previous analyses have shown that pathogen-mediated natural selection maintains exceptional polymorphism in HLA residues that determine peptide recognition. Here, we show that HLA class I peptide recognition impacts diversity of globally widespread pathogens.

Vaccine Strain and Wild-Type Clades of Varicella-Zoster Virus in Central Nervous System and Non-CNS Disease, New York State, 2004-2019

Since the introduction of the varicella-zoster virus (VZV) vaccine in the United States in 1995, there has been a dramatic decrease in both the number and severity of varicella cases. However, VZV surveillance data and information on the VZV clade distribution in central nervous system (CNS) disease and non-CNS disease in New York State is not available. To investigate this, cerebrospinal fluid (CSF) samples from patients with encephalitis or meningitis and non-CSF samples from patients with non-CNS disease manifestations consistent with VZV, collected from 2004 to 2019, were tested with molecular VZV assays. A total of 341 CSF and 1,398 non-CSF samples that tested positive by a VZV-specific real-time PCR assay were further characterized as wild-type or vaccine strain by 3 biallelic real-time PCR assays targeting single nucleotide polymorphism (SNP) markers in open reading frame (ORF) 62. Genotyping was then performed on wild-type strains by conventional PCR and sequencing of 500-bp regions in ORFs 21, 22, and 50. Sequence analysis identified clades 1 to 5 in both sample types with a virtually identical clade distribution between CSF and non-CSF samples. In addition, 19 clade 6 and 13 clade 9 samples were detected in non-CSF samples after implementation of an expanded genotyping scheme, including ORF 29, 38, and 67. These clades were not detected in any CSF samples. Finally, a total of 28 vaccine strains were detected, 25 in the non-CSF samples and 3 in the CSF samples. All three cases of vaccine strain with CNS involvement experienced relatively minor symptoms of aseptic meningitis and fully recovered. These results support the evidence that while the VZV vaccine is capable of causing CNS disease, it is still a rare event and symptoms are typically less severe than those caused by wild-type infection.

Whole-genome analysis of natural interspecific recombinant between bovine alphaherpesviruses 1 and 5

Bovine alphaherpesviruses 1 and 5 (BoHV-1 and BoHV-5) are closely related viruses that co-circulate in South America and recombine in the field. The complete genomes of three natural gB gene recombinant viruses between BoHV-1 and BoHV-5 were obtained by Illumina next-generation sequencing. Complete genome sequences of the three recombinant strains (RecA1, RecB2, and RecC2) have a similar size of approximately 138.3kb and a GC content of 75%. The genome structure corresponds to herpesvirus class D, with 69 open reading frames (ORFs) arranged in the same order as other bovine alphaherpesviruses related to BoHV-1. Their genomes were included in recombination network studies indicating statistically significant recombination evidence both based on the whole genome, as well as in the sub-regions. The novel recombinant region of 3074 nt of the RecB2 and RecC2 strains includes the complete genes of the myristylated tegument protein (UL11) and the glycoprotein M (UL10) and part of the helicase (UL9) gene, and it seems to have originated independently of the first recombinant event involving the gB gene. Phylogenetic analyzes performed with the amino acid sequences of UL9, UL 10, and UL11 indicated that RecB2 and RecC2 recombinants are closely related to the minor parental virus (BoHV-1.2b). On the contrary, RecA1 groups with the major parental (BoHV-5), thus confirming the absence of recombination in this region for this recombinant. One breakpoint in the second recombinant region lies in the middle of the UL9 reading frame, originating a chimeric enzyme half encoded by BoHV-5 and BoHV-1.2b parental strains. The chimeric helicases of both recombinants are identical and have 96.8 and 96.3% similarity with the BoHV-5 and BoHV-1 parents, respectively. In vitro characterization suggests that recombinants have delayed exit from the cell compared to parental strains. However, they produce the similar viral titer as their putative parents suggesting the accumulation of viral particles for the cell exit delayed on time. Despite in vitro different behavior, these natural recombinant viruses have been maintained in the bovine population for more than 30 years, indicating that recombination could be playing an important role in the biological diversity of these viral species. Our findings highlight the importance of studying whole genome diversity in the field and determining the role that homologous recombination plays in the structure of viral populations. A whole-genome recombinant characterization is a suitable tool to help understand the emergence of new viral forms with novel pathogenic features.

Varicella-Zoster Virus-Genetics, Molecular Evolution and Recombination

This chapter first details the structure, organization and coding content of the VZV genome to provide a foundation on which the molecular evolution of the virus can be projected. We subsequently describe the evolution of molecular profiling approaches from restriction fragment length polymorphisms to single nucleotide polymorphism profiling to modern day high-throughput sequencing approaches. We describe how the application of these methodologies led to our current model of VZV phylogeograpy including the number and structure of geographic clades and the role of recombination in reshaping these.

Possible European Origin of Circulating Varicella Zoster Virus Strains

Varicella zoster virus (VZV) is the causative agent of chickenpox and shingles. The geographic distribution of VZV clades was taken as evidence that VZV migrated out of Africa with human populations. We show that extant VZV strains most likely originated in Europe and not in Africa. Europe was also identified as the ancestral location for most internal nodes of the VZV phylogeny, including the ancestor of clade 5 strains. We also show that strains from clades 1, 2, 3, and 5 derived a major proportion of their ancestry from each of 4 ancestral populations. Conversely, viruses from other clades displayed variable levels of admixture. Some low-level admixture was also observed for clade 5 genomes, but only for non-African viruses. This pattern indicates that the clade 5 VZV strains do not represent recent introductions from Africa due to migratory fluxes. These data have also relevance for the definition and classification of VZV clades.

Alphaherpesvirus Genomics: Past, Present and Future

Alphaherpesviruses, as large double-stranded DNA viruses, were long considered to be genetically stable and to exist in a homogeneous state. Recently, the proliferation of high-throughput sequencing (HTS) and bioinformatics analysis has expanded our understanding of herpesvirus genomes and the variations found therein. Recent data indicate that herpesviruses exist as diverse populations, both in culture and in vivo, in a manner reminiscent of RNA viruses. In this review, we discuss the past, present, and potential future of alphaherpesvirus genomics, including the technical challenges that face the field. We also review how recent data has enabled genome-wide comparisons of sequence diversity, recombination, allele frequency, and selective pressures, including those introduced by cell culture. While we focus on the human alphaherpesviruses, we draw key insights from related veterinary species and from the beta- and gamma-subfamilies of herpesviruses. Promising technologies and potential future directions for herpesvirus genomics are highlighted as well, including the potential to link viral genetic differences to phenotypic and disease outcomes.

Exocytosis of Progeny Infectious Varicella-Zoster Virus Particles via a Mannose-6-Phosphate Receptor Pathway without Xenophagy following Secondary Envelopment

The literature on the egress of different herpesviruses after secondary envelopment is contradictory. In this report, we investigated varicella-zoster virus (VZV) egress in a cell line from a child with Pompe disease, a glycogen storage disease caused by a defect in the enzyme required for glycogen digestion. In Pompe cells, both the late autophagy pathway and the mannose-6-phosphate receptor (M6PR) pathway are interrupted. We have postulated that intact autophagic flux is required for higher recoveries of VZV infectivity. To test that hypothesis, we infected Pompe cells and then assessed the VZV infectious cycle. We discovered that the infectious cycle in Pompe cells was remarkably different from that of either fibroblasts or melanoma cells. No large late endosomes filled with VZV particles were observed in Pompe cells; only individual viral particles in small vacuoles were seen. The distribution of the M6PR pathway (trans-Golgi network to late endosomes) was constrained in infected Pompe cells. When cells were analyzed with two different anti-M6PR antibodies, extensive colocalization of the major VZV glycoprotein gE (known to contain M6P residues) and the M6P receptor (M6PR) was documented in the viral highways at the surfaces of non-Pompe cells after maximum-intensity projection of confocal z-stacks, but neither gE nor the M6PR was seen in abundance at the surfaces of infected Pompe cells. Taken together, our results suggested that (i) Pompe cells lack a VZV trafficking pathway within M6PR-positive large endosomes and (ii) most infectious VZV particles in conventional cell substrates are transported via large M6PR-positive vacuoles without degradative xenophagy to the plasma membrane.IMPORTANCE The long-term goal of this research has been to determine why VZV, when grown in cultured cells, invariably is more cell associated and has a lower titer than other alphaherpesviruses, such as herpes simplex virus 1 (HSV1) or pseudorabies virus (PRV). Data from both HSV1 and PRV laboratories have identified a Rab6 secretory pathway for the transport of single enveloped viral particles from the trans-Golgi network within small vacuoles to the plasma membrane. In contrast, after secondary envelopment in fibroblasts or melanoma cells, multiple infectious VZV particles accumulated within large M6PR-positive late endosomes that were not degraded en route to the plasma membrane. We propose that this M6PR pathway is most utilized in VZV infection and least utilized in HSV1 infection, with PRV's usage being closer to HSV1's usage. Supportive data from other VZV, PRV, and HSV1 laboratories about evidence for two egress pathways are included.

Molecular Genetic Insights Into Varicella Zoster Virus (VZV), the vOka Vaccine Strain, and the Pathogenesis of Latency and Reactivation

Genetic tools for molecular typing of varicella zoster virus (VZV) have been used to understand the spread of virus, to differentiate wild-type and vaccine strains, and to understand the natural history of VZV infection in its cognate host. Molecular genetics has identified 7 clades of VZV (1-6 and 9), with 2 more mooted. Differences between the vOka vaccine strain and wild-type VZVs have been used to distinguish the cause of postimmunization events and to provide insight into the natural history of VZV infections. Importantly molecular genetics has shown that reinfection with establishment of latency by the reinfecting strain is common, that dual infections with different viruses can occur, and that reactivation of the superinfecting genotype can both occur. Whole-genome sequencing of the vOka vaccine has been used to show that vesicles form from a single virion, that latency is established within a few days of inoculation, and that all vaccine strains are capable of establishing latency and reactivating. Novel molecular tools have characterized the transcripts expressed during latent infection in vitro.

Progeny Varicella-Zoster Virus Capsids Exit the Nucleus but Never Undergo Secondary Envelopment during Autophagic Flux Inhibition by Bafilomycin A1

Varicella-zoster virus (VZV) is an alphaherpesvirus that lacks the herpesviral neurovirulence protein ICP34.5. The underlying hypothesis of this project was that inhibitors of autophagy reduce VZV infectivity. We selected the vacuolar proton ATPase inhibitor bafilomycin A1 for analysis because of its well-known antiautophagy property of impeding acidification during the late stage of autophagic flux. We documented that bafilomycin treatment from 48 to 72 h postinfection lowered VZV titers substantially (P ≤ 0.008). Because we were unable to define the site of the block in the infectious cycle by confocal microscopy, we turned to electron microscopy. Capsids were observed in the nucleus, in the perinuclear space, and in the cytoplasm adjacent to Golgi apparatus vesicles. Many of the capsids had an aberrant appearance, as has been observed previously in infections not treated with bafilomycin. In contrast to prior untreated infections, however, secondary envelopment of capsids was not seen in the trans-Golgi network, nor were prototypical enveloped particles with capsids (virions) seen in cytoplasmic vesicles after bafilomycin treatment. Instead, multiple particles with varying diameters without capsids (light particles) were seen in large virus assembly compartments near the disorganized Golgi apparatus. Bafilomycin treatment also led to increased numbers of multivesicular bodies in the cytoplasm, some of which contained remnants of the Golgi apparatus. In summary, we have defined a previously unrecognized property of bafilomycin whereby it disrupted the site of secondary envelopment of VZV capsids by altering the pH of the trans-Golgi network and thereby preventing the correct formation of virus assembly compartments.IMPORTANCE This study of VZV assembly in the presence of bafilomycin A1 emphasizes the importance of the Golgi apparatus/trans-Golgi network as a platform in the alphaherpesvirus life cycle. We have previously shown that VZV induces levels of autophagy far above the basal levels of autophagy in human skin, a major site of VZV assembly. The current study documented that bafilomycin treatment led to impaired assembly of VZV capsids after primary envelopment/de-envelopment but before secondary reenvelopment. This VZV study also complemented prior herpes simplex virus 1 and pseudorabies virus studies investigating two other inhibitors of endoplasmic reticulum (ER)/Golgi apparatus function: brefeldin A and monensin. Studies with porcine herpesvirus demonstrated that primary enveloped particles accumulated in the perinuclear space in the presence of brefeldin A, while studies with herpes simplex virus 1 documented an impaired secondary assembly of enveloped viral particles in the presence of monensin.

Autophagy Quantification and STAT3 Expression in a Human Skin Organ Culture Model for Innate Immunity to Herpes Zoster

The goal of this project was to document the autophagy response in human neonatal skin organ culture (SOC) after infection with varicella-zoster virus (VZV). The VZV-infected SOC model has attributes of herpes zoster, in that an injection of virus into the skin is analogous to exit of virus from the sensory nerve termini into skin during herpes zoster. Cultures were maintained for 28 days and periodically examined for an autophagy response by quantitation of autophagosomes with Imaris software. Expression of the STAT3 protein was plentiful in the VZV-infected SOC. Abundant autophagy was observed in VZV-infected SOC between 14 and 28 days after infection, while autophagy in mock-infected SOC was minimal (p = 0.0003). The autophagic response after infection of SOC with a recombinant VZV genome containing the herpes simplex virus ICP34.5 neurovirulence gene was similar to wild-type VZV (p = 0.3). These results suggested that the VZV-infected SOC system resembled biopsy data from herpes zoster infection of skin. An enhanced autophagy response has now been reported after infection with two additional alpha herpesviruses besides VZV, namely, pseudorabies virus and duck enteritis herpes virus; both lack the ICP34.5 protein.

Clinico-epidemiological investigation on Varicella Zoster Virus indicates multiple clade circulation in Maharashtra state, India

BACKGROUND

Varicella Zoster Virus (VZV) is consistently in circulation and shows an increase in disease burden during the spring season. Due to a wide range of clinical presentation from a vesicular rash to bleeding or neurological complications, it makes the clinical diagnosis difficult. The present study aims to understand whether the same strain of virus is responsible for the increase in the seasonal outbreaks occurring in different parts of the country with reference to the samples from Maharashtra, Rajasthan and Gujarat states of India.

MATERIALS AND METHODS

This study reports the clinico-epidemiological and laboratory findings of suspected Varicella cases. To understand the circulating clade few representative real-time Polymerase Chain Reaction (PCR) positive were analyzed by conventional PCR and partial Open Reading Frame (ORF) 22, partial ORF 38 and partial ORF 54 were sequenced to identify single nucleotide polymorphisms responsible for clade determination. Further partial glycoprotein B gene was sequenced, and a phylogenetic tree was generated.

RESULTS

A total of 50 cases from Maharashtra (Mumbai district) and referred clinical samples of Rajasthan (Barmer district; n = 12) and Gujarat States (Gandhi Nagar, Surat districts; n = 17) were tested for the presence of VZV. Vesicular rash with fever was a common clinical presentation with 82% cases having contact history with VZV positive cases, suggesting higher secondary attack rate. The vesicular fluid of all 50 cases from Mumbai revealed the presence of VZV by real-time PCR. Urine, serum and throat swab samples showed positivity by real-time PCR. Healthcare provider's samples from Rajasthan showed 36.4% [4/11] positivity. Clinical samples from Gujarat had positivity of 41.2% [7/17].

CONCLUSIONS

This study analyses the clade based circulation of VZV in three states in India and suggests different clades circulating in Maharashtra state. Health education amongst the general population is suggested to reduce the secondary cases by early diagnosis, effective isolation policies and vaccination to reduce the burden of disease.

Molecular Aspects of Varicella-Zoster Virus Latency

Primary varicella-zoster virus (VZV) infection causes varicella (chickenpox) and the establishment of a lifelong latent infection in ganglionic neurons. VZV reactivates in about one-third of infected individuals to cause herpes zoster, often accompanied by neurological complications. The restricted host range of VZV and, until recently, a lack of suitable in vitro models have seriously hampered molecular studies of VZV latency. Nevertheless, recent technological advances facilitated a series of exciting studies that resulted in the discovery of a VZV latency-associated transcript (VLT) and provide novel insights into our understanding of VZV latency and factors that may initiate reactivation. Deducing the function(s) of VLT and the molecular mechanisms involved should now be considered a priority to improve our understanding of factors that govern VZV latency and reactivation. In this review, we summarize the implications of recent discoveries in the VZV latency field from both a virus and host perspective and provide a roadmap for future studies.

Genomic, Recombinational and Phylogenetic Characterization of Global Feline Herpesvirus 1 Isolates

Feline herpes virus type 1 (FHV-1) is widely considered to be the leading cause of ocular disease in cats and has been implicated in upper respiratory tract infections. Little, however is known about interstrain phylogenetic relationships, and details of the genomic structure. For the present study, twenty-six FHV-1 isolates from different cats in animal shelters were collected from eight separate locations in the USA, and the genomes sequenced. Genomic characterization of these isolates includied short sequence repeat (SSR) detection, with fewer SSRs detected, compared to herpes simplex viruses type 1 and 2. For phylogenetic and recombination analysis, 27 previously sequenced isolates of FHV-1 were combined with the 26 strains sequenced for the present study. The overall genomic interstrain genetic distance between all available isolates was 0.093%. Phylogenetic analysis identified four main FHV-1 clades primarily corresponding to geographical collection site. Recombination analysis suggested that interclade recombination has occurred.

Varicella Virus Vaccination in the United States

Varicella zoster virus (VZV) is the cause of chickenpox (varicella) and shingles (zoster), and was once responsible for over 4 million infections in the United States annually. The development of a live attenuated VZV vaccine was initially viewed with extreme skepticism. Nonetheless, a VZV vaccine was developed in the 1970s by Takahashi and his colleagues in Japan and was eventually licensed in the US. It is now known to be one of the safest and most effective vaccines available and is administered worldwide. Here are described important factors that contributed to the successful research and licensure of the highly successful VZV vaccine.

Characterization and phylogenetic analysis of Varicella-zoster virus strains isolated from Korean patients

Varicella-zoster virus (VZV) is a causative agent of chickenpox in primary infection and shingles after its reactivation from latency. Complete or almost-complete genomic DNA sequences for various VZV strains have been reported. Recently, clinical VZV strains were isolated from Korean patients whose genome was sequenced using high-throughput sequencing technology. In this study, we analyzed single nucleotide polymorphism (SNP) of VZV strains to genetically characterize Korean clinical isolates. Phylogenetic analyses revealed that three Korean strains, YC01, YC02, and YC03, were linked to clade 2. Comprehensive SNP analysis identified 86 sites specific for the 5 VZV clades. VZV strains isolated from Korea did not form a phylogenetic cluster. Rather, YC02 and YC03 clustered strongly with Chinese strain 84-7 within clade 2, more specifically cluster 2a. Signature sequences for the cluster 2a were identified and found to play an important role in the separation of cluster 2a strains from other clade 2 strains, as shown in substitution studies. Further genetic analysis with additional strains isolated from Japan, China, and other Asian countries would provide a novel insight into the significance of two distinct subclades within clade 2.

Differentiation between wild-type and vaccines strains of varicella zoster virus (VZV) based on four single nucleotide polymorphisms

Varicella-zoster virus (VZV) infection (chickenpox) results in latency and subsequent reactivation manifests as shingles. Effective attenuated vaccines (vOka) are available for prevention of both illnesses. In this study, an amplicon-based sequencing method capable of differentiating between VZV wild-type (wt) strains and vOka vaccine is described. A total of 44 vesicular fluid specimens collected from 43 patients (16 from China and 27 from the UK) with either chickenpox or shingles were investigated, of which 10 had received previous vaccination. Four sets of polymerase chain reactions were set up simultaneously with primers amplifying regions encompassing four single nucleotide polymorphisms (SNPs), '69349-106262-107252-108111'. Nucleotide sequences were generated by Sanger sequencing. All samples except one had a wt SNP profile of 'A-T-T-T'. The sample collected from a patient who received vaccine 7-10 days ago, along with VZV vaccine preparations, Zostavax and Baike-varicella gave a SNP profile 'G-C-C-C'. The results show that this method can distinguish vaccine-derived virus from wt viruses from main four clades, (clades 1-4) and should be of utility worldwide.

Genetic analysis of varicella-zoster virus in the aqueous humor in uveitis with severe hyphema

BACKGROUND

Genetic variations have been identified in the genome of varicella-zoster virus (VZV) strains using vesicle fluid, varicella scabs and throat swab samples. We report a rare case of VZV-associated uveitis with severe hyphema, which was immediately diagnosed by polymerase chain reaction (PCR) using the aqueous humor, in which we were able to analyze the VZV genotype for the first time.

CASE PRESENTATION

A 16-year-old Japanese boy was referred to our hospital with a 20-day history of unilateral anterior uveitis and 11-day history of hyphema. At presentation, details of the iris, the iridocorneal angle, and the fundus were not visible due to the severe hyphema. Serum anti-VZV IgG and anti-VZV IgM were elevated, and 1.61 × 10⁹ copies/mL of VZV-DNA were detected by real-time PCR using the aqueous humor. As there were no eruptions on his face or body, we diagnosed zoster sine herpete and started intravenous administration of prednisolone and acyclovir. The hyphema completely disappeared 2 weeks after presentation, while sectorial iris atrophy and mild periphlebitis of the fundus became gradually apparent. Anterior inflammation and periphlebitis gradually improved and VZV-DNA in the aqueous humor was reduced to 1.02 × 10⁶ copies/mL at 4 weeks after presentation. Examination by slit lamp microscope revealed no inflammation after 5 months, and VZV-DNA could no longer be detected in the aqueous humor. Serum anti-VZV IgG and anti-VZV IgM also showed a gradual decrease along with improvement in ocular inflammation. The genetic analysis of multiple open reading frames and the R5 variable repeat region in the VZV genes, using DNA extracted from the aqueous humor at presentation, showed that the isolate was a wild-type clade 2 VZV strain (prevalent in Japan and surrounding countries) with R5A allele and one SNP unique to clade 1 (both are major types in Europe and North America).

CONCLUSIONS

VZV-associated uveitis may develop hyphema that obscures ocular inflammation, thus PCR analysis using the aqueous humor is the key investigation necessary for the diagnosis. The measurement of VZV-DNA copies by real-time PCR would be useful for evaluation of therapeutic effects. We could amplify and analyze VZV genotype using the aqueous humor including a very large number of VZV-DNA copies (1.61 × 10⁹ copies/mL).

Nationwide distribution of varicella-zoster virus clades in China

Background

In 2010, a universal nomenclature for varicella-zoster virus (VZV) clades was established, which is very useful in the monitoring of viral evolution, recombination, spread and genetic diversity. Currently, information about VZV clades has been disclosed worldwide, however, there are limited data regarding the characterization of circulating VZV clades in China, even where varicella remains widely epidemic.

Methods

From 2008 to 2012, clinical samples with varicella or zoster were collected in General Hospital in eight provinces and analyzed by PCR, restriction endonuclease digestion and sequencing. The viral clades were determined by analysis of five single nucleotide polymorphisms (SNPs) within the 447-bp fragment of open reading frame (ORF) 22, and the restriction fragment length polymorphisms (RFLPs) of ORF 38 (PstI), ORF 54 (BglI) and ORF 62 (SmaI) were evaluated to understand genetic diversity of VZV and determinate varicella vaccine adverse event (VVAE).

Results

Seventy-seven varicella and 11 zoster samples were identified as being positive for VZV. The five SNPs profile showed that the majority of VZV strains belonged to clade 2, but clade 5 and clade 4 strains were also found in Guangdong. The RFLPs analysis of the DNA fragments of ORF 38, 54 and 62 showed that 85 of these samples were characterized as PstI + BglI + SamI-, and the remaining three VZV strains from varicella patients were characterized as PstI-BglI + SamI+ which is the genetic profile of VVAEs.

Conclusions

The study suggested that the predominant clade 2 VZVs had been continually circulating since at least the 1950s in China. Nearly all VZV strains except VVAEs possessed the genetic profile of PstI + BglI + Sam-. However, the other clades were also found to be co-circulating with clade 2, especially in the border regions. These results highlighted the need for the constant and broad use of virologic surveillance to provide an important genetic baseline for varicella control and vaccination programs in China.

DNA from Dust: Comparative Genomics of Large DNA Viruses in Field Surveillance Samples

The intensification of the poultry industry over the last 60 years facilitated the evolution of increased virulence and vaccine breaks in Marek's disease virus (MDV-1). Full-genome sequences are essential for understanding why and how this evolution occurred, but what is known about genome-wide variation in MDV comes from laboratory culture. To rectify this, we developed methods for obtaining high-quality genome sequences directly from field samples without the need for sequence-based enrichment strategies prior to sequencing. We applied this to the first characterization of MDV-1 genomes from the field, without prior culture. These viruses were collected from vaccinated hosts that acquired naturally circulating field strains of MDV-1, in the absence of a disease outbreak. This reflects the current issue afflicting the poultry industry, where virulent field strains continue to circulate despite vaccination and can remain undetected due to the lack of overt disease symptoms. We found that viral genomes from adjacent field sites had high levels of overall DNA identity, and despite strong evidence of purifying selection, had coding variations in proteins associated with virulence and manipulation of host immunity. Our methods empower ecological field surveillance, make it possible to determine the basis of viral virulence and vaccine breaks, and can be used to obtain full genomes from clinical samples of other large DNA viruses, known and unknown. IMPORTANCE Despite both clinical and laboratory data that show increased virulence in field isolates of MDV-1 over the last half century, we do not yet understand the genetic basis of its pathogenicity. Our knowledge of genome-wide variation between strains of this virus comes exclusively from isolates that have been cultured in the laboratory. MDV-1 isolates tend to lose virulence during repeated cycles of replication in the laboratory, raising concerns about the ability of cultured isolates to accurately reflect virus in the field. The ability to directly sequence and compare field isolates of this virus is critical to understanding the genetic basis of rising virulence in the wild. Our approaches remove the prior requirement for cell culture and allow direct measurement of viral genomic variation within and between hosts, over time, and during adaptation to changing conditions.

Microbiology laboratory and the management of mother-child varicella-zoster virus infection

Varicella-zoster virus, which is responsible for varicella (chickenpox) and herpes zoster (shingles), is ubiquitous and causes an acute infection among children, especially those aged less than six years. As 90% of adults have had varicella in childhood, it is unusual to encounter an infected pregnant woman but, if the disease does appear, it can lead to complications for both the mother and fetus or newborn. The major maternal complications include pneumonia, which can lead to death if not treated. If the virus passes to the fetus, congenital varicella syndrome, neonatal varicella (particularly serious if maternal rash appears in the days immediately before or after childbirth) or herpes zoster in the early years of life may occur depending on the time of infection. A Microbiology laboratory can help in the diagnosis and management of mother-child infection at four main times: (1) when a pregnant woman has been exposed to varicella or herpes zoster, a prompt search for specific antibodies can determine whether she is susceptible to, or protected against infection; (2) when a pregnant woman develops clinical symptoms consistent with varicella, the diagnosis is usually clinical, but a laboratory can be crucial if the symptoms are doubtful or otherwise unclear (atypical patterns in immunocompromised subjects, patients with post-vaccination varicella, or subjects who have received immunoglobulins), or if there is a need for a differential diagnosis between varicella and other types of dermatoses with vesicle formation; (3) when a prenatal diagnosis of uterine infection is required in order to detect cases of congenital varicella syndrome after the onset of varicella in the mother; and (4) when the baby is born and it is necessary to confirm a diagnosis of varicella (and its complications), make a differential diagnosis between varicella and other diseases with similar symptoms, or confirm a causal relationship between maternal varicella and malformations in a newborn.

Analysis of single nucleotide polymorphism among Varicella-Zoster Virus and identification of vaccine-specific sites

Varicella-zoster virus (VZV) is a causative agent for chickenpox and zoster. Live attenuated vaccines have been developed based on Oka and MAV/06 strains. In order to understand the molecular mechanisms of attenuation, complete genome sequences of vaccine and wild-type strains were compared and single nucleotide polymorphism (SNP) was analyzed. ORF22 and ORF62 contained the highest number of SNPs. The detailed analysis of the SNPs suggested 24 potential vaccine-specific sites. All the mutational events found in vaccine-specific sites were transitional, and most of them were substitution of AT to GC pair. Interestingly, 18 of the vaccine-specific sites of the vaccine strains appeared to be genetically heterogeneous. The probability of a single genome of vaccine strain to contain all 24 vaccine-type sequences was calculated to be less than 4%. The average codon adaptation index (CAI) value of the vaccine strains was significantly lower than the CAI value of the clinical strains.

Genotyping of clinical varicella-zoster virus isolates collected from Yunnan in Southwestern China

Varicella-zoster virus (VZV) belongs to the α-herpesvirus family. Genetically, it is stable and is divided into several genotypes based upon the genetic variations. The genotypes of VZV are rarely studied in the Southwestern region of China. In the present study, the common genetic variations in the VZV genes were examined in 42 VZV isolates collected from the patients with herpes zoster in the Yunnan province (Southwestern China). The restriction fragment length polymorphism analysis of open reading frames (ORFs) 38, 54 and 62 in the VZV genes showed that all the collected VZV isolates were PstI, BglI and SmaI positive. The R5 variable-repeat region in these isolates was variable (R5A: 46.4%; R5B: 53.6%). The sequencing data of ORFs 1, 21, 22 and 54 indicated that 41 of the 42 collected VZV isolates could be grouped into genotype J or J1. Only one VZV isolate was identified as genotype A1 or M2. No new substitutions in the sequenced fragments were found in the collected VZV isolates. The results of the present study provided a preliminary genetic characterization of the VZV strains in the Yunnan province of Southwestern China.

Evidence of widespread natural recombination among field isolates of equine herpesvirus 4 but not among field isolates of equine herpesvirus 1

Recombination in alphaherpesviruses allows evolution to occur in viruses that have an otherwise stable DNA genome with a low rate of nucleotide substitution. High-throughput sequencing of complete viral genomes has recently allowed natural (field) recombination to be studied in a number of different alphaherpesviruses, however, such studies have not been applied to equine herpesvirus 1 (EHV-1) or equine herpesvirus 4 (EHV-4). These two equine alphaherpesviruses are genetically similar, but differ in their pathogenesis and epidemiology. Both cause economically significant disease in horse populations worldwide. This study used high-throughput sequencing to determine the full genome sequences of EHV-1 and EHV-4 isolates (11 and 14 isolates, respectively) from Australian or New Zealand horses. These sequences were then analysed and examined for evidence of recombination. Evidence of widespread recombination was detected in the genomes of the EHV-4 isolates. Only one potential recombination event was detected in the genomes of the EHV-1 isolates, even when the genomes from an additional 11 international EHV-1 isolates were analysed. The results from this study reveal another fundamental difference between the biology of EHV-1 and EHV-4. The results may also be used to help inform the future safe use of attenuated equine herpesvirus vaccines.

Genotype analysis of varicella-zoster virus isolates from suburban Shanghai Municipal Province, China

To determine the predominant genotype of the varicella-zoster virus (VZV) in suburban Shanghai Municipal Province, specimens were collected from the lesions of 95 outpatients clinically diagnosed with varicella or herpes zoster. Of these, 69 patients (72.6%) were positive for VZV DNA. The 69 isolates were all genotyped as the genotype J1/clade 2. Based on sequencing of the 447 bp sequence in ORF22, 66 isolates were identified as genotype J/clade 2 strains and three were identified as type M2/clade 4 strains. To confirm the classification of these three strains, we determined the presence of 27 single-nucleotide polymorphisms (SNPs) and found that isolates 1270/1450 shared seven SNPs that differed from those of clade 2, in which three SNPs were unique to clade 3 and another three were unique to clade 4. Isolate 1456 had two markers of clade 4 that differed from clade 2. The phylogenetic tree showed that our isolates clustered primarily with clade 2 and that the three M2/J1 strains clustered between clades 2 and 4. It is likely that isolates 1270/1450/1446 may represent a new subclade of either clade 2 or 4, or some recombinant events. In addition, our isolates were WT strains. We also observed significant inter-strain variations.

Recombination of Globally Circulating Varicella-Zoster Virus

Varicella-zoster virus (VZV) is a human herpesvirus, which during primary infection typically causes varicella (chicken pox) and establishes lifelong latency in sensory and autonomic ganglia. Later in life, the virus may reactivate to cause herpes zoster (HZ; also known as shingles). To prevent these diseases, a live-attenuated heterogeneous vaccine preparation, vOka, is used routinely in many countries worldwide. Recent studies of another alphaherpes virus, infectious laryngotracheitis virus, demonstrate that live-attenuated vaccine strains can recombine in vivo, creating virulent progeny. These findings raised concerns about using attenuated herpesvirus vaccines under conditions that favor recombination. To investigate whether VZV may undergo recombination, which is a prerequisite for VZV vaccination to create such conditions, we here analyzed 115 complete VZV genomes. Our results demonstrate that recombination occurs frequently for VZV. It thus seems that VZV is fully capable of recombination if given the opportunity, which may have important implications for continued VZV vaccination. Although no interclade vaccine wild-type recombinant strains were found, intraclade recombinants were frequently detected in clade 2, which harbors the vaccine strains, suggesting that the vaccine strains have already been involved in recombination events, either in vivo or in vitro during passages in cell culture. Finally, previous partial and complete genomic studies have described strains that do not cluster phylogenetically to any of the five established clades. The additional VZV strains sequenced here, in combination with those previously published, have enabled us to formally define a novel sixth VZV clade.

IMPORTANCE

Although genetic recombination has been demonstrated to frequently occur for other human alphaherpesviruses, herpes simplex viruses 1 and 2, only a few ancient and isolated recent recombination events have hitherto been demonstrated for VZV. In the present study, we demonstrate that VZV also frequently undergoes genetic recombination, including strains belonging to the clade containing the vOKA strain.

Full-Genome Sequence of a Novel Varicella-Zoster Virus Clade Isolated in Mexico

Varicella-zoster virus (VZV) is a member of the Herpesviridae family, which causes varicella (chicken pox) and herpes zoster (shingles) in humans. Here, we report the complete genome sequence of varicella-zoster virus, isolated from a vesicular fluid sample, revealing the circulation of VZV clade VIII in Mexico.

Autophagic flux without a block differentiates varicella-zoster virus infection from herpes simplex virus infection

Autophagy is a process by which misfolded and damaged proteins are sequestered into autophagosomes, before degradation in and recycling from lysosomes. We have extensively studied the role of autophagy in varicella-zoster virus (VZV) infection, and have observed that vesicular cells are filled with >100 autophagosomes that are easily detectable after immunolabeling for the LC3 protein. To confirm our hypothesis that increased autophagosome formation was not secondary to a block, we examined all conditions of VZV infection as well as carrying out two assessments of autophagic flux. We first investigated autophagy in human skin xenografts in the severe combined immunodeficiency (SCID) mouse model of VZV pathogenesis, and observed that autophagosomes were abundant in infected human skin tissues. We next investigated autophagy following infection with sonically prepared cell-free virus in cultured cells. Under these conditions, autophagy was detected in a majority of infected cells, but was much less than that seen after an infected-cell inoculum. In other words, inoculation with lower-titered cell-free virus did not reflect the level of stress to the VZV-infected cell that was seen after inoculation of human skin in the SCID mouse model or monolayers with higher-titered infected cells. Finally, we investigated VZV-induced autophagic flux by two different methods (radiolabeling proteins and a dual-colored LC3 plasmid); both showed no evidence of a block in autophagy. Overall, therefore, autophagy within a VZV-infected cell was remarkably different from autophagy within an HSV-infected cell, whose genome contains two modifiers of autophagy, ICP34.5 and US11, not present in VZV.

Evolution of cocirculating varicella-zoster virus genotypes during a chickenpox outbreak in Guinea-Bissau

Varicella-zoster virus (VZV), a double-stranded DNA alphaherpesvirus, is associated with seasonal outbreaks of varicella in nonimmunized populations. Little is known about whether these outbreaks are associated with a single or multiple viral genotypes and whether new mutations rapidly accumulate during transmission. Here, we take advantage of a well-characterized population cohort in Guinea-Bissau and produce a unique set of 23 full-length genome sequences, collected over 7 months from eight households. Comparative sequence analysis reveals that four distinct genotypes cocirculated among the population, three of which were present during the first week of the outbreak, although no patients were coinfected, which indicates that exposure to infectious virus from multiple sources is common during VZV outbreaks. Transmission of VZV was associated with length polymorphisms in the R1 repeat region and the origin of DNA replication. In two cases, these were associated with the formation of distinct lineages and point to the possible coevolution of these loci, despite the lack of any known functional link in VZV or related herpesviruses. We show that these and all other sequenced clade 5 viruses possess a distinct R1 repeat motif that increases the acidity of an ORF11p protein domain and postulate that this has either arisen or been lost following divergence of the major clades. Thus, sequencing of whole VZV genomes collected during an outbreak has provided novel insights into VZV biology, transmission patterns, and (recent) natural history.

IMPORTANCE VZV is a highly infectious virus and the causative agent of chickenpox and shingles, the latter being particularly associated with the risk of painful complications. Seasonal outbreaks of chickenpox are very common among young children, yet little is known about the dynamics of the virus during person-to-person to transmission or whether multiple distinct viruses seed and/or cocirculate during an outbreak. In this study, we have sequenced chickenpox viruses from an outbreak in Guinea-Bissau that are supported by detailed epidemiological data. Our data show that multiple different virus strains seeded and were maintained throughout the 6-month outbreak period and that viruses transmitted between individuals accumulated new mutations in specific genomic regions. Of particular interest is the potential coevolution of two distinct parts of the genomes and our calculations of the rate of viral mutation, both of which increase our understanding of how VZV evolves over short periods of time in human populations.

Varicella-zoster virus glycoprotein expression differentially induces the unfolded protein response in infected cells

Varicella-zoster virus (VZV) is a human herpesvirus that spreads to children as varicella or chicken pox. The virus then establishes latency in the nervous system and re-emerges, typically decades later, as zoster or shingles. We have reported previously that VZV induces autophagy in infected cells as well as exhibiting evidence of the Unfolded Protein Response (UPR): XBP1 splicing, a greatly expanded Endoplasmic Reticulum (ER) and CHOP expression. Herein we report the results of a UPR specific PCR array that measures the levels of mRNA of 84 different components of the UPR in VZV infected cells as compared to tunicamycin treated cells as a positive control and uninfected, untreated cells as a negative control. Tunicamycin is a mixture of chemicals that inhibits N-linked glycosylation in the ER with resultant protein misfolding and the UPR. We found that VZV differentially induces the UPR when compared to tunicamycin treatment. For example, tunicamycin treatment moderately increased (8-fold) roughly half of the array elements while downregulating only three (one ERAD and two FOLD components). VZV infection on the other hand upregulated 33 components including a little described stress sensor CREB-H (64-fold) as well as ER membrane components INSIG and gp78, which modulate cholesterol synthesis while downregulating over 20 components mostly associated with ERAD and FOLD. We hypothesize that this expression pattern is associated with an expanding ER with downregulation of active degradation by ERAD and apoptosis as the cell attempts to handle abundant viral glycoprotein synthesis.

A single viral gene determines lethal cross-species neurovirulence of baboon herpesvirus HVP2

Alpha-herpesviruses can produce more severe infections in non-natural host species than in their natural host. Isolates of the baboon alpha-herpesvirus Papiine herpesvirus 2 (HVP2) are either very neurovirulent in mice (subtype nv) or non-virulent (subtype ap), but no such difference is evident in the natural baboon host. Comparative genome sequencing was used to identify subtype-specific sequence differences (SSDs) between HVP2nv and HVP2ap isolates. Some genes were identified that despite exhibiting sequence variation among isolates did not have any SSDs, while other genes had comparatively high levels of SSDs. Construction of genomic recombinants between HVP2nv and HVP2ap isolates mapped the mouse neurovirulence determinant to within three genes. Construction of gene-specific recombinants demonstrated that the UL39 ORF is responsible for determining the lethal neurovirulence phenotype of HVP2 in mice. These results demonstrate that differences in a single viral gene can determine the severity of herpesvirus infection in a non-natural host species.

Intratypic variants of human papillomavirus type 16 and risk of cervical neoplasia in Taiwan

The associations between variants of human papillomavirus (HPV) 16 and risk of cervical neoplasia have been reported, but nucleotide variations of HPV 16 in Asian populations and their association with cervical neoplasia have not been evaluated extensively. During 1991-1992, 11,923 women from seven townships in Taiwan were enrolled. The HPV DNA in cervical cells was detected and genotyped using EasyChip HPV blot. Nucleotide variations in the long control region (LCR), E6, and E7 genes were determined using DNA sequencing for 170 HPV 16-positive cervical samples. The Asian variant was the most prevalent variant (81.8%) of HPV 16 in Taiwan, and was also associated with increased prevalence of histologically confirmed cervical intraepithelial neoplasia grade 3 or worse, showing an age-adjusted odds ratio (exact confidence limits) of 10.70 (1.62-451.05; P = 0.0049) compared to the HPV 16 European variant. Similar significant associations with cervical intraepithelial neoplasia grade 3 or worse were also observed for distinct nucleotide substitutions, including T178A/G, A647G, A7730C/G, T7781C, G7842A, and C24T/G. These results demonstrate that non-European variants (non-E) of HPV 16, predominantly Asian variants, are associated with increased risk for severe cervical neoplasia, compared with European variants. Molecular mechanisms accounting for varied cervical neoplasia risk among different HPV 16 variants warrant further investigation.

Evolution and diversity in human herpes simplex virus genomes

Herpes simplex virus 1 (HSV-1) causes a chronic, lifelong infection in >60% of adults. Multiple recent vaccine trials have failed, with viral diversity likely contributing to these failures. To understand HSV-1 diversity better, we comprehensively compared 20 newly sequenced viral genomes from China, Japan, Kenya, and South Korea with six previously sequenced genomes from the United States, Europe, and Japan. In this diverse collection of passaged strains, we found that one-fifth of the newly sequenced members share a gene deletion and one-third exhibit homopolymeric frameshift mutations (HFMs). Individual strains exhibit genotypic and potential phenotypic variation via HFMs, deletions, short sequence repeats, and single-nucleotide polymorphisms, although the protein sequence identity between strains exceeds 90% on average. In the first genome-scale analysis of positive selection in HSV-1, we found signs of selection in specific proteins and residues, including the fusion protein glycoprotein H. We also confirmed previous results suggesting that recombination has occurred with high frequency throughout the HSV-1 genome. Despite this, the HSV-1 strains analyzed clustered by geographic origin during whole-genome distance analysis. These data shed light on likely routes of HSV-1 adaptation to changing environments and will aid in the selection of vaccine antigens that are invariant worldwide.

Autophagy and the effects of its inhibition on varicella-zoster virus glycoprotein biosynthesis and infectivity

Autophagy and the effects of its inhibition or induction were investigated during the entire infectious cycle of varicella-zoster virus (VZV), a human herpesvirus. As a baseline, we first enumerated the number of autophagosomes per cell after VZV infection compared with the number after induction of autophagy following serum starvation or treatment with tunicamycin or trehalose. Punctum induction by VZV was similar in degree to punctum induction by trehalose in uninfected cells. Treatment of infected cells with the autophagy inhibitor 3-methyladenine (3-MA) markedly reduced the viral titer, as determined by assays measuring both cell-free virus and infectious foci (P < 0.0001). We next examined a virion-enriched band purified by density gradient sedimentation and observed that treatment with 3-MA decreased the amount of VZV gE, while treatment with trehalose increased the amount of gE in the same band. Because VZV gE is the most abundant glycoprotein, we selected gE as a representative viral glycoprotein. To further investigate the role of autophagy in VZV glycoprotein biosynthesis as well as confirm the results obtained with 3-MA inhibition, we transfected cells with ATG5 small interfering RNA to block autophagosome formation. VZV-induced syncytium formation was markedly reduced by ATG5 knockdown (P < 0.0001). Further, we found that both expression and glycan processing of VZV gE were decreased after ATG5 knockdown, while expression of the nonglycosylated IE62 tegument protein was unchanged. Taken together, our cumulative results not only documented abundant autophagy within VZV-infected cells throughout the infectious cycle but also demonstrated that VZV-induced autophagy facilitated VZV glycoprotein biosynthesis and processing.

Deep sequencing of viral genomes provides insight into the evolution and pathogenesis of varicella zoster virus and its vaccine in humans

Immunization with the vOka vaccine prevents varicella (chickenpox) in children and susceptible adults. The vOka vaccine strain comprises a mixture of genotypes and, despite attenuation, causes rashes in small numbers of recipients. Like wild-type virus, the vaccine establishes latency in neuronal tissue and can later reactivate to cause Herpes zoster (shingles). Using hybridization-based methodologies, we have purified and sequenced vOka directly from skin lesions. We show that alleles present in the vaccine can be recovered from the lesions and demonstrate the presence of a severe bottleneck between inoculation and lesion formation. Genotypes in any one lesion appear to be descended from one to three vaccine-genotypes with a low frequency of novel mutations. No single vOka haplotype and no novel mutations are consistently present in rashes, indicating that neither new mutations nor recombination with wild type are critical to the evolution of vOka rashes. Instead, alleles arising from attenuation (i.e., not derived from free-living virus) are present at lower frequencies in rash genotypes. We identify 11 loci at which the ancestral allele is selected for in vOka rash formation and show genotypes in rashes that have reactivated from latency cannot be distinguished from rashes occurring immediately after inoculation. We conclude that the vOka vaccine, although heterogeneous, has not evolved to form rashes through positive selection in the mode of a quasispecies, but rather alleles that were essentially neutral during the vaccine production have been selected against in the human subjects, allowing us to identify key loci for rash formation.

Phylogenetic and molecular epidemiological studies reveal evidence of multiple past recombination events between infectious laryngotracheitis viruses

In contrast to the RNA viruses, the genome of large DNA viruses such as herpesviruses have been considered to be relatively stable. Intra-specific recombination has been proposed as an important, but underestimated, driving force in herpesvirus evolution. Recently, two distinct field strains of infectious laryngotracheitis virus (ILTV) have been shown to have arisen from independent recombination events between different commercial ILTV vaccines. In this study we sequenced the genomes of additional ILTV strains and also utilized other recently updated complete genome sequences of ILTV to confirm the existence of a number of ILTV recombinants in nature. Multiple recombination events were detected in the unique long and repeat regions of the genome, but not in the unique short region. Most recombinants contained a pair of crossover points between two distinct lineages of ILTV, corresponding to the European origin and the Australian origin vaccine strains of ILTV. These results suggest that there are two distinct genotypic lineages of ILTV and that these commonly recombine in the field.

Molecular characterization of varicella-zoster virus clinical isolates from 2006 to 2008 in a tertiary care hospital, Dublin, Ireland, using different genotyping methods

Varicella-zoster virus (VZV), a herpesvirus, is a ubiquitous organism that causes considerable morbidity worldwide and can cause severe complications on reactivation. Phylogenetic analysis was performed on 19 clinical VZV isolates (16 zoster and 3 varicella) found in Ireland, between December 2006 and November 2008, in order to determine whether previously reported viral heterogeneity was still present and whether viral recombination was evident. Open reading-frames (ORFs) from genes 1, 21, 50, and 54, were sequenced. Clades 1, 2, 3, and 5 were identified. Four putative recombinant isolates were detected (three clade 3/1 and one clade 5/3/1). Further sequencing and examination of ORF 22 and 21/50, did not elucidate the putative recombinant genotypes further. These two previously published genotyping schemes were examined in light of the new consensus genotyping scheme proposed in 2010. Remarkable VZV heterogeneity remains prevalent in Ireland. This is the first evidence of putative VZV recombination found in Ireland.

Analysis of enzymatic digestion pattern of two open reading frames of Varciella-Zoster genome from Kuwaiti patients using the RFLP technique

BACKGROUND AND OBJECTIVES

Varicella-Zoster virus (VZV) is a human herpes virus that usually attacks young children and commonly causes chicken pox (Varicella). Following primary infection, a lifelong latent infection is established. The virus often reactivates during adulthood or senesces to cause shingles (Zoster). Little is known regarding the genotypes of Varicella in Kuwait. The aim of this study was to genotype Varicella samples collected from patients in Kuwait.

MATERIALS AND METHODS

Samples from 60 cases of chicken pox were typed. The DNA extraction was performed using the commercially available DNA extraction kit. Two sets of oligonucleotide primers were used to amplify the intervening sequences with polymerase chain reaction to identify VZV DNA in clinical samples. The BglI and PstI endonucleases were used to digest. The PCR amplicons for PCR-RFLP typing.

RESULTS

Relatively consistent restriction enzyme digestion profiles for different VZV strains were observed. Limited genetic differences between VZV samples were found. Three VZV strains were identified (A, B and C) with type B representing 86.6%, type A 11.7% and type C being 1.7%. We found that distinct restriction fragment length polymorphism isolates from the same origin or nationality were very similar.

CONCLUSION

Varicella strains with cutting sites for both enzyme PstI and BglI (typeB) were more prevalent. Molecular amplification of viral DNA by PCR and restriction digestion could be used for VZV typing as an alternative method to serological assays.

Bioinformatics of varicella-zoster virus: single nucleotide polymorphisms define clades and attenuated vaccine genotypes

Varicella zoster virus (VZV) is one of the human herpesviruses. To date, over 40 complete VZV genomes have been sequenced and analyzed. The VZV genome contains around 125,000 base pairs including 70 open reading frames (ORFs). Enumeration of single nucleotide polymorphisms (SNPs) has determined that the following ORFs are the most variable (in descending order): 62, 22, 29, 28, 37, 21, 54, 31, 1 and 55. ORF 62 is the major immediate early regulatory VZV gene. Further SNP analysis across the entire genome has led to the observation that VZV strains can be broadly grouped into clades within a phylogenetic tree. VZV strains collected in Singapore provided important sequence data for construction of the phylogenetic tree. Currently five VZV clades are recognized; they have been designated clades 1 through 5. Clades 1 and 3 include European/North American strains; clade 2 includes Asian strains, especially from Japan; and clade 5 includes strains from India. Clade 4 includes some strains from Europe, but its geographic origins need further documentation. Within clade 1, five variant viruses have been isolated with a missense mutation in the gE (ORF 68) glycoprotein; these strains have an altered increased cell spread phenotype. Bioinformatics analyses of the attenuated vaccine strains have also been performed, with a subsequent discovery of a stop-codon SNP in ORFO as a likely attenuation determinant. Taken together, these VZV bioinformatics analyses have provided enormous insights into VZV phylogenetics as well as VZV SNPs associated with attenuation.

The attenuated genotype of varicella-zoster virus includes an ORF0 transitional stop codon mutation

Varicella-zoster virus (VZV) is the first of the human herpesviruses to be attenuated and subsequently approved as a live vaccine to prevent varicella and herpes zoster. Both the attenuated VZV vaccine, called vaccine Oka or vOka, and the parental strain pOka have been completely sequenced. Yet the specific determinants of attenuation are uncertain. The open reading frame (ORF) with the most single nucleotide polymorphisms (SNPs), ORF62, encodes the regulatory protein IE62, but IE62 studies have failed to define a specific SNP associated with attenuation. We have completed next-generation sequencing of the VZV Ellen genome, a strain known to be highly attenuated by its very limited replication in human skin xenografts in the SCID mouse model of VZV pathogenesis. A comparative analysis of the Ellen sequence with all other complete VZV sequences was extremely informative. In particular, an unexpected finding was a stop codon mutation in Ellen ORF0 (herpes simplex virus UL56 homolog) identical to one found in vOka, combined with the absence of polymorphisms in most Ellen ORFs that were known to be mutated in vOka. The mutated ORF0 protein was also imaged in both two dimensions and three dimensions by confocal microscopy. The probability of two VZV strains not connected by a recent common ancestor having an identical ORF0 SNP by chance would be 1 × 10(-8), in other words, extremely unlikely. Taken together, these bioinformatics analyses strongly suggest that the stop codon ORF0 SNP is one of the determinants of the attenuation genotype of live VZV vaccines.

Pangaea and the Out-of-Africa Model of Varicella-Zoster Virus Evolution and Phylogeography

The goal of this minireview is to provide an overview of varicella-zoster virus (VZV) phylogenetics and phylogeography when placed in the broad context of geologic time. Planet Earth was formed over 4 billion years ago, and the supercontinent Pangaea coalesced around 400 million years ago (mya). Based on detailed tree-building models, the base of the phylogenetic tree of the Herpesviridae family has been estimated at 400 mya. Subsequently, Pangaea split into Laurasia and Gondwanaland; in turn, Africa rifted from Gondwanaland. Based on available data, the hypothesis of this minireview is that the ancestral alphaherpesvirus VZV coevolved in simians, apes, and hominins in Africa. When anatomically modern humans first crossed over the Red Sea 60,000 years ago, VZV was carried along in their dorsal root ganglia. Currently, there are five VZV clades, distinguishable by single nucleotide polymorphisms. These clades likely represent continued VZV coevolution, as humans with latent VZV infection left Arabia and dispersed into Asia (clades 2 and 5) and Europe (clades 1, 3, and 4). The prototype VZV sequence contains nearly 125,000 bp, divided into 70 open reading frames. Generally, isolates within a clade display >99.9% identity to one another, while members of one clade compared to a second clade show 99.8% identity to one another. Recently, four different VZV genotypes that do not segregate into the previously defined five clades have been identified, a result indicating a wider than anticipated diversity among newly collected VZV strains around the world.

Novel genetic variation identified at fixed loci in ORF62 of the Oka varicella vaccine and in a case of vaccine-associated herpes zoster

The live attenuated Oka varicella vaccine (vOka), derived from clade 2 wild-type (wt) virus pOka, is used for routine childhood immunization in several countries, including the United States, which has caused dramatic declines in the incidence of varicella. vOka can cause varicella, establish latency, and reactivate to cause herpes zoster (HZ). Three loci in varicella-zoster virus (VZV) open reading frame 62 (ORF62) (106262, 107252, and 108111) are used to distinguish vOka from wt VZV. A fourth position (105705) is also fixed for the vOka allele in nearly all vaccine batches. These 4 positions and two vOka mutations (106710 and 107599) reportedly absent from Varivax were analyzed on Varivax-derived ORF62 TOPO TA clones. The wt allele was detected at positions 105705 and 107252 on 3% and 2% of clones, respectively, but was absent at positions 106262 and 108111. Position 106710 was fixed for the wt allele, whereas the vOka allele was present on 18.4% of clones at position 107599. We also evaluated the 4 vOka markers in an isolate obtained from a case of vaccine-caused HZ. The isolate carried the vOka allele at positions 105705, 106262, and 108111. However, at position 107252, the wt allele was present. Thus, all of the ORF62 vOka markers previously regarded as fixed occur as the wt allele in a small percentage of vOka strains. Characterization of all four vOka markers in ORF62 and of the clade 2 subtype marker in ORF38 is now necessary to confirm vOka adverse events.

Sequencing and characterization of Varicella-zoster virus vaccine strain SuduVax

BACKGROUND

Varicella-zoster virus (VZV) causes chickenpox in children and shingles in older people. Currently, live attenuated vaccines based on the Oka strain are available worldwide. In Korea, an attenuated VZV vaccine has been developed from a Korean isolate and has been commercially available since 1994. Despite this long history of use, the mechanism for the attenuation of the vaccine strain is still elusive. We attempted to understand the molecular basis of attenuation mechanism by full genome sequencing and comparative genomic analyses of the Korean vaccine strain SuduVax.

RESULTS

SuduVax was found to contain a genome that was 124,759 bp and possessed 74 open reading frames (ORFs). SuduVax was genetically most close to Oka strains and these Korean-Japanese strains formed a strong clade in phylogenetic trees. SuduVax, similar to the Oka vaccine strains, underwent T- > C substitution at the stop codon of ORF0, resulting in a read-through mutation to code for an extended form of ORF0 protein. SuduVax also shared certain deletion and insertion mutations in ORFs 17, 29, 56 and 60 with Oka vaccine strains and some clinical strains.

CONCLUSIONS

The Korean VZV vaccine strain SuduVax is genetically similar to the Oka vaccine strains. Further comparative genomic and bioinformatics analyses will help to elucidate the molecular basis of the attenuation of the VZV vaccine strains.

Sequencing of 21 varicella-zoster virus genomes reveals two novel genotypes and evidence of recombination

Genotyping of 21 varicella-zoster virus (VZV) strains using a scattered single nucleotide polymorphism (SNP) method revealed ambiguous SNPs and two nontypeable isolates. For a further genetic characterization, the genomes of all strains were sequenced using the 454 technology. Almost-complete genome sequences were assembled, and most remaining gaps were closed with Sanger sequencing. Phylogenetic analysis of 42 genomes revealed five established and two novel VZV genotypes, provisionally termed VIII and IX. Genotypes VIII and IX are distinct from the previously reported provisional genotypes VI and VII as judged from the SNP pattern. The alignments showed evidence of ancient recombination events in the phylogeny of clade 4 and recent recombinations within single strains: 3/2005 (clade 1), 11 and 405/2007 (clade 3), 8 and DR (clade 4), CA123 and 413/2000 (clade 5), and strains of the novel genotypes VIII and IX. Bayesian tree inference of the thymidine kinase and the polymerase genes of the VZV clades and other varicelloviruses revealed that VZV radiation began some 110,000 years ago, which correlates with the out-of-Africa dispersal of modern humans. The split of ancestral clades 2/4 and 1/3/5/VIII/IX shows the greatest node height.

Specific capture and whole-genome sequencing of viruses from clinical samples

Whole genome sequencing of viruses directly from clinical samples is integral for understanding the genetics of host-virus interactions. Here, we report the use of sample sparing target enrichment (by hybridisation) for viral nucleic acid separation and deep-sequencing of herpesvirus genomes directly from a range of clinical samples including saliva, blood, virus vesicles, cerebrospinal fluid, and tumour cell lines. We demonstrate the effectiveness of the method by deep-sequencing 13 highly cell-associated human herpesvirus genomes and generating full length genome alignments at high read depth. Moreover, we show the specificity of the method enables the study of viral population structures and their diversity within a range of clinical samples types.

Comparative analysis of the complete genome sequences of two Australian origin live attenuated vaccines of infectious laryngotracheitis virus

Infectious laryngotracheitis virus (ILTV) is an alphaherpesvirus that causes acute respiratory disease in poultry. Live attenuated ILTV vaccines have been used extensively to help control outbreaks of disease. Two Australian-origin attenuated vaccine strains, SA2 and A20 ILTV, are commercially available and are in frequent use in Australia. Both these vaccines are of chicken embryo origin (CEO). The A20 ILTV strain was developed from the SA2 ILTV strain by sequential passage of SA2 ILTV in tissue culture in order to reduce its residual virulence. Previous studies in our laboratories have demonstrated the greater attenuation of A20 ILTV under controlled experimental conditions, but the genetic basis of the in vivo phenotypes of A20 and SA2 ILTV has not been elucidated. In this study, the genetic differences between A20 and SA2 ILTV were examined by performing complete genome sequencing and comparative analysis. The genome sequences were also compared to a reference sequence from another CEO ILTV vaccine (Serva ILTV: GenBank accession number HQ_630064) of European-origin. Additional in ovo studies to assess cell to cell spread were performed in order to allow further comparisons of the pathogenicity of SA2 and A20 ILTV. The sequencing results showed that the genome sizes of SA2 and A20 ILTV were 152,975 and 152,978bp, respectively, while Serva ILTV had a genome size of 152,630bp. The genomes of SA2 and A20 ILTV shared 99.9% nucleotide sequence identity with each other, but only 99.2% identity with Serva ILTV. In complete genome alignments between SA2 and A20 ILTV, a total of 24 single nucleotide polymorphisms (SNPs) were identified, but only two of these were non-synonymous. These were located in the ORF B and UL15 genes. Four indels were detected in non-coding regions. The findings from this study demonstrate the general genetic stability of ILTV, but also show that non-synonymous changes in the ORF B and UL15 genes have arisen following tissue culture passage of SA2 ILTV to produce the A20 vaccine. It is likely that these non-synonymous changes are related to the greater attenuation of A20 ILTV compared to SA2 ILTV, and to the reduced ability of A20 ILTV to spread from cell to cell, as observed in this study. The results from this study also demonstrate the divergence between the genomes of the Australian-origin ILTV vaccine strains and the Serva vaccine strain.

A wide extent of inter-strain diversity in virulent and vaccine strains of alphaherpesviruses

Alphaherpesviruses are widespread in the human population, and include herpes simplex virus 1 (HSV-1) and 2, and varicella zoster virus (VZV). These viral pathogens cause epithelial lesions, and then infect the nervous system to cause lifelong latency, reactivation, and spread. A related veterinary herpesvirus, pseudorabies (PRV), causes similar disease in livestock that result in significant economic losses. Vaccines developed for VZV and PRV serve as useful models for the development of an HSV-1 vaccine. We present full genome sequence comparisons of the PRV vaccine strain Bartha, and two virulent PRV isolates, Kaplan and Becker. These genome sequences were determined by high-throughput sequencing and assembly, and present new insights into the attenuation of a mammalian alphaherpesvirus vaccine strain. We find many previously unknown coding differences between PRV Bartha and the virulent strains, including changes to the fusion proteins gH and gB, and over forty other viral proteins. Inter-strain variation in PRV protein sequences is much closer to levels previously observed for HSV-1 than for the highly stable VZV proteome. Almost 20% of the PRV genome contains tandem short sequence repeats (SSRs), a class of nucleic acids motifs whose length-variation has been associated with changes in DNA binding site efficiency, transcriptional regulation, and protein interactions. We find SSRs throughout the herpesvirus family, and provide the first global characterization of SSRs in viruses, both within and between strains. We find SSR length variation between different isolates of PRV and HSV-1, which may provide a new mechanism for phenotypic variation between strains. Finally, we detected a small number of polymorphic bases within each plaque-purified PRV strain, and we characterize the effect of passage and plaque-purification on these polymorphisms. These data add to growing evidence that even plaque-purified stocks of stable DNA viruses exhibit limited sequence heterogeneity, which likely seeds future strain evolution.

Alphaherpesviruses such as herpes simplex virus (HSV) are ubiquitous in the human population. HSV causes oral and genital lesions, and has co-morbidities in acquisition and spread of human immunodeficiency virus (HIV). The lack of a vaccine for HSV hinders medical progress for both of these infections. A related veterinary alphaherpesvirus, pseudorabies virus (PRV), has long served as a model for HSV vaccine development, because of their similar pathogenesis, neuronal spread, and infectious cycle. We present here the first full genome characterization of a live PRV vaccine strain, Bartha, and reveal a spectrum of unique mutations that are absent from two divergent wild-type PRV strains. These mutations can now be examined individually for their contribution to vaccine strain attenuation and for potential use in HSV vaccine development. These inter-strain comparisons also revealed an abundance of short repetitive elements in the PRV genome, a pattern which is repeated in other herpesvirus genomes and even the unrelated Mimivirus. We provide the first global characterization of repeats in viruses, comparing both their presence and their variation among different viral strains and species. Repetitive elements such as these have been shown to serve as hotspots of variation between individuals or strains of other organisms, generating adaptations or even disease states through changes in length of DNA-binding sites, protein folding motifs, and other structural elements. These data suggest for the first time that similar mechanisms could be widely distributed in viral biology as well.