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

Machine Learning Prediction of Adenovirus D8 Conjunctivitis Complications from Viral Whole-Genome Sequence

Objective

To obtain complete DNA sequences of adenoviral (AdV) D8 genome from patients with conjunctivitis and determine the relation of sequence variation to clinical outcomes.

Design

This study is a post hoc analysis of banked conjunctival swab samples from the BAYnovation Study, a previously conducted, randomized controlled clinical trial for AdV conjunctivitis.

Participants

Ninety-six patients with AdV D8-positive conjunctivitis who received placebo treatment in the BAYnovation Study were included in the study.

Methods

DNA from conjunctival swabs was purified and subjected to whole-genome viral DNA sequencing. Adenovirus D8 variants were identified and correlated with clinical outcomes, including 2 machine learning methods.

Main Outcome Measures

Viral DNA sequence and development of subepithelial infiltrates (SEIs) were the main outcome measures.

Results

From initial sequencing of 80 AdV D8-positive samples, full adenoviral genome reconstructions were obtained for 71. A total of 630 single-nucleotide variants were identified, including 156 missense mutations. Sequence clustering revealed 3 previously unappreciated viral clades within the AdV D8 type. The likelihood of SEI development differed significantly between clades, ranging from 83% for Clade 1 to 46% for Clade 3. Genome-wide analysis of viral single-nucleotide polymorphisms failed to identify single-gene determinants of outcome. Two machine learning models were independently trained to predict clinical outcome using polymorphic sequences. Both machine learning models correctly predicted development of SEI outcomes in a newly sequenced validation set of 16 cases (P = 1.5 × 10-5). Prediction was dependent on ensemble groups of polymorphisms across multiple genes.

Conclusions

Adenovirus D8 has ≥ 3 prevalent molecular substrains, which differ in propensity to result in SEIs. Development of SEIs can be accurately predicted from knowledge of full viral sequence. These results suggest that development of SEIs in AdV D8 conjunctivitis is largely attributable to pathologic viral sequence variants within the D8 type and establishes machine learning paradigms as a powerful technique for understanding viral pathogenicity.

Varicella Vaccine Meningitis as a Complication of Herpes Zoster in Twice-Immunized Immunocompetent Adolescents

Varicella-zoster virus vaccination is recommended for virtually all young children in the United States, Canada, and several other countries. Varicella vaccine is a live attenuated virus that retains some of its neurotropic properties. Herpes zoster caused by vaccine virus still occurs in immunized children, although the rate is much lower than in children who had wild-type varicella. It was commonly thought that 2 varicella vaccinations would protect children against the most serious complication of meningitis following herpes zoster; however, 2 meningitis cases have already been published. We now report a third case of varicella vaccine meningitis and define risk factors shared by all 3 immunized adolescents. The diagnosis in cerebrospinal fluid in this third case was verified by amplifying and sequencing portions of the viral genome, to document fixed alleles found only in the vaccine strain. Viral antibody was also detected in the cerebrospinal fluid by confocal microscopy. When compared with the other 2 cases, remarkably all 3 were 14 years old when meningitis occurred. All 3 were treated with intravenous acyclovir, with complete recovery. The adolescent in our case report also had recurrent asthma, which was treated with both prednisone tablets and beclomethasone inhaler before onset of meningitis. When the 3 cases were considered together, they suggested that immunity to varicella-zoster virus may be waning sufficiently in some twice-immunized adolescents to make them vulnerable to varicella vaccine virus reactivation and subsequent meningitis. This complication rarely happens in children after wild-type varicella.

Molecular Biology and Epidemiology of Neurotropic Viruses

Neurotropic viruses are those viruses that can cause central nervous system (CNS) diseases with both neuroinvasive and neurovirulence properties. It comprises a wide range of viruses, including herpes simplex virus, poliovirus, enteroviruses, parechovirus, West Nile virus, Japanese encephalitis virus, measles, and mumps viruses among others. Some of these viruses are highly neuroinvasive and neurovirulent, while others are weakly neuroinvasive and neurovirulent. Moreover, some of them, like herpes simplex viruses, are highly neuroinvasive but weakly neurovirulent for the peripheral nervous system and highly neurovirulent but weakly neuroinvasive for the central nervous system. All these disparities are a result of differences in their genomic constitution, associated vectors, geographical region, and environmental factors. Therefore, a successful intervention will be almost impossible without a clear understanding of the molecular biology and epidemiology of these viruses. Thus, we conducted a review of the published studies on the molecular biology and epidemiology of the common neurotropic viruses to make the viral genetic makeup more understandable for targeted intervention and provide the morbidity and mortality data of the different neurotropic viruses for more serious action.

The phylogenetic approach for viral infectious disease evolution and epidemiology: An updating review

In the last decade, the phylogenetic approach is recurrent in molecular evolutionary analysis. On 12 May, 2019, about 2 296 213 papers are found, but typing "phylogeny" or "epidemiology AND phylogeny" only 199 804 and 20 133 are retrieved, respectively. Molecular epidemiology in infectious diseases is widely used to define the source of infection as so as the ancestral relationships of individuals sampled from a population. Coalescent theory and phylogeographic analysis have had scientific application in several, recent pandemic events, and nosocomial outbreaks. Hepatitis viruses and immunodeficiency virus (human immunodeficiency virus) have been largely studied. Phylogenetic analysis has been recently applied on Polyomaviruses so as in the more recent outbreaks due to different arboviruses type as Zika and chikungunya viruses discovering the source of infection and the geographic spread. Data on sequences isolated by the microorganism are essential to apply the phylogenetic tools and research in the field of infectious disease phylodinamics is growing up. There is the need to apply molecular phylogenetic and evolutionary methods in areas out of infectious diseases, as translational genomics and personalized medicine. Lastly, the application of these tools in vaccine strategy so as in antibiotic and antiviral researchers are encouraged.

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.

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.

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.

Approved Antiviral Drugs over the Past 50 Years

Since the first antiviral drug, idoxuridine, was approved in 1963, 90 antiviral drugs categorized into 13 functional groups have been formally approved for the treatment of the following 9 human infectious diseases: (i) HIV infections (protease inhibitors, integrase inhibitors, entry inhibitors, nucleoside reverse transcriptase inhibitors, nonnucleoside reverse transcriptase inhibitors, and acyclic nucleoside phosphonate analogues), (ii) hepatitis B virus (HBV) infections (lamivudine, interferons, nucleoside analogues, and acyclic nucleoside phosphonate analogues), (iii) hepatitis C virus (HCV) infections (ribavirin, interferons, NS3/4A protease inhibitors, NS5A inhibitors, and NS5B polymerase inhibitors), (iv) herpesvirus infections (5-substituted 2'-deoxyuridine analogues, entry inhibitors, nucleoside analogues, pyrophosphate analogues, and acyclic guanosine analogues), (v) influenza virus infections (ribavirin, matrix 2 protein inhibitors, RNA polymerase inhibitors, and neuraminidase inhibitors), (vi) human cytomegalovirus infections (acyclic guanosine analogues, acyclic nucleoside phosphonate analogues, pyrophosphate analogues, and oligonucleotides), (vii) varicella-zoster virus infections (acyclic guanosine analogues, nucleoside analogues, 5-substituted 2'-deoxyuridine analogues, and antibodies), (viii) respiratory syncytial virus infections (ribavirin and antibodies), and (ix) external anogenital warts caused by human papillomavirus infections (imiquimod, sinecatechins, and podofilox). Here, we present for the first time a comprehensive overview of antiviral drugs approved over the past 50 years, shedding light on the development of effective antiviral treatments against current and emerging infectious diseases worldwide.

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.

Varicella zoster virus genotyping in chickenpox patient's clinical isolates from Iran

Aim: The varicella zoster virus (VZV) causes chickenpox and zoster infections. This study aimed to investigate the distribution of VZV genotypes among Iranian patients. Materials & methods: From 2010 to 2015, 244 patients were enrolled in this cross-sectional study, 45 of whom were positive for VZV DNA. Both direct sequencing and restriction fragment length polymorphism assay were performed for 19 positive specimens. SPSS v.20 was used for statistics. Results: The predominant VZV genotype was M1 (84.2%) followed by genotype E (10.5%) and genotype J (5.3%). Restriction fragment length polymorphism demonstrated that 17 strains were PstI+ BglI+ (M1 and/or J genotypes) and 2 were PstI+ BglI- (E genotype). Conclusion: This research is a prelim study on VZV genotyping. Further investigations will help to confirm the VZV genotype prevalence reported here.

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.

Enzymatic Digestion Pattern of Varicella Zoster Virus ORF38 and ORF54 in Chickenpox Patients Using RFLP Technique

BACKGROUND

Varicella zoster virus (VZV) causes chickenpox in children and zoster (zona) in the elderly. Using RFLP-PCR method for detection of VZV specific SNPs ORF38, 54 and 62 could distinguish the profile of VZV isolates. The aim of this study was to investigate enzymatic digestion pattern of VZV ORF38 and ORF54 in chickenpox patients using RFLP technique.

METHODS

Thirty-eight chickenpox patients, who referred to the hospitals of Iran University of Medical Sciences in Tehran from May 2010 to June 2015 were enrolled in this cross sectional study. After the DNA extraction, PCR amplification of 38 VZV isolates performed by specific primers of ORFs 38 and 54, then RFLP assay and digestion carried out by PstI (for ORF38) and BglI (for ORF54) restriction enzymes.

RESULTS

Of 38 positive VZV DNA, the mean age (yr)±SD was 34.4±23.3 (range: 7-89). 22 (57.9%) were female and 16 (42.1%) were male. The predominant VZV profile of BglI(+) PstI(+) were 89.5% (34/38) followed by 10.5% (4/38) PstI(+) BglI‾. Statistical analysis showed that there was no significant relationship between genotype, age, sex, and year of infection variables (P value> 0.05). The common VZV genotype among Iranian patients with chickenpox and zona infection is genotype BglI(+) PstI(+) followed by PstI(+) BglI‾.

CONCLUSION

There are different VZV circulating genotypes that call for for more research on this field by widely population and other methods such as nucleotide sequencing to justify the accurate VZV genotype prevalence in Iran.

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.

Rates of vaccine evolution show strong effects of latency: implications for varicella zoster virus epidemiology

Varicella-zoster virus (VZV) causes chickenpox and shingles, and is found in human populations worldwide. The lack of temporal signal in the diversity of VZV makes substitution rate estimates unreliable, which is a barrier to understanding the context of its global spread. Here, we estimate rates of evolution by studying live attenuated vaccines, which evolved in 22 vaccinated patients for known periods of time, sometimes, but not always undergoing latency. We show that the attenuated virus evolves rapidly (∼ 10(-6) substitutions/site/day), but that rates decrease dramatically when the virus undergoes latency. These data are best explained by a model in which viral populations evolve for around 13 days before becoming latent, but then undergo no replication during latency. This implies that rates of viral evolution will depend strongly on transmission patterns. Nevertheless, we show that implausibly long latency periods are required to date the most recent common ancestor of extant VZV to an "out-of-Africa" migration with humans, as has been previously suggested.

Varicella-zoster virus and herpes simplex virus 1 can infect and replicate in the same neurons whether co- or superinfected

The two human neurotropic alphaherpesviruses varicella-zoster virus (VZV) and herpes simplex virus type 1 (HSV1) both establish latency in sensory ganglia. Human trigeminal ganglia are known to frequently harbor both viruses, and there is evidence to suggest the presence of both VZV and HSV1 DNA in the same neuron. We ask here whether VZV and HSV1 can exclude themselves and each other and whether they can productively infect the same cells in human neurons and human foreskin fibroblasts (HFF). Simultaneous infection (coinfection) or consecutive infection (superinfection) was assessed using cell-free HSV1 and VZV expressing fluorescent reporter proteins. Automated analysis was carried out to detect singly and dually infected cells. We demonstrate that VZV and HSV1 both display efficient superinfection exclusion (SE) in HFF, with each virus excluding either itself or the other virus. While SE also occurred in neurons, it was with much lower efficiency. Both alphaherpesviruses productively infected the same neurons, whether applied simultaneously or even consecutively, albeit at lower frequencies.

IMPORTANCE

Superinfection exclusion by VZV for itself or the related neurotropic alphaherpesvirus HSV1 has been studied here for the first time. We find that while these viruses display classic SE in fibroblasts, SE is less efficient for both HSV1 and VZV in human neurons. The ability of multiple VZV strains to productively infect the same neurons has important implications in terms of recombination of both wild-type and vaccine strains in patients.

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.

Breakthrough varicella in a cancer patient with persistent varicella antibody after one varicella vaccination

A boy with Hodgkin disease contracted breakthrough varicella from his father, who had chickenpox. The boy had received a single varicella vaccination and was seropositive by enzyme-linked immunosorbent assay before being diagnosed with breakthrough varicella. Seropositivity after a single varicella vaccination does not guarantee complete protection in an immunocompromised child.