Impacts of Genome-Wide Analyses on Our Understanding of Human Herpesvirus Diversity and Evolution

High-Fidelity Long-Read Sequencing of an Avian Herpesvirus Reveals Extensive Intrapopulation Diversity in Tandem Repeat Regions

Comparative genomic studies of Marek's disease virus (MDV) have suggested that attenuated and virulent strains share >98% sequence identity. However, these estimates fail to account for variation in regions of the MDV genome harboring tandem repeats. To resolve these loci and enable assessments of intrapopulation diversity, we used a PacBio Sequel II platform to sequence MDV strains CVI988/Rispens (attenuated), HPRS-B14 (virulent), Md5 (very virulent) and 675A (very virulent plus). This approach enabled us to identify patterns of variation in tandem repeat regions that are consistent with known phenotypic differences across these strains. We also found CVI988/Rispens variants showing a 4.3-kb deletion in the Unique Short (US) region, resulting in the loss of five genes. These findings support a potential link between MDV tandem repeats and phenotypic traits like virulence and attenuation, and demonstrate that DNA viruses can harbor high levels of intrapopulation diversity in tandem repeat regions.

N6-methyladenosine in inflammatory diseases: Important actors and regulatory targets

N6-methyladenosine (m6A) is one of the most prevalent epigenetic modifications in eukaryotic cells. It regulates RNA function and stability by modifying RNA methylation through writers, erasers, and readers. As a result, m6A plays a critical role in a wide range of biological processes. Inflammation is a common and fundamental pathological process. Numerous studies have investigated the role of m6A modifications in inflammatory diseases. This review highlights the mechanisms by which m6A contributes to inflammation, focusing on pathogen-induced infectious diseases, autoimmune disorders, allergic conditions, and metabolic disorder-related inflammatory diseases.

Update of the Genetic Variability of Monkeypox Virus Clade IIb Lineage B.1

From 1 January 2022 to 31 May 2024, the World Health Organization (WHO) reported 97,745 laboratory-confirmed Mpox cases, including 203 deaths, across 116 countries. Despite a 2.3% decrease in new cases in May 2024 compared to April 2024, significant regional variations persist. The African Region reported the highest proportion of new cases, while other regions experienced mixed trends. Phylogenomic analyses of the Mpox virus Clade IIb lineage B.1 reveal stable genetic variability with minimal diversification. The Bayesian Skyline Plot indicates a generally stable viral population size with a modest peak in late 2023, followed by a decline. In general, the data indicate that the MPXV outbreak is primarily localized within a few consistent geographic clusters. The virus's evolution is relatively slow, as indicated by its stable genetic variability, and Clade IIb lineage B.1 does not currently show signs of rapid genetic changes or population growth. The current low level of genetic diversity should not lead to complacency. Ongoing genomic surveillance is essential for effective outbreak management and understanding. This monitoring is crucial for identifying any shifts in the virus's behavior or transmission, allowing for prompt public health responses and adjustments. In addition, continued vigilance is necessary to detect any new variants that might influence the outbreak's trajectory.

Comparative analysis of multiple consensus genomes of the same strain of Marek's disease virus reveals intrastrain variation

Current strategies to understand the molecular basis of Marek's disease virus (MDV) virulence primarily consist of cataloging divergent nucleotides between strains with different phenotypes. However, most comparative genomic studies of MDV rely on previously published consensus genomes despite the confirmed existence of MDV strains as mixed viral populations. To assess the reliability of interstrain genomic comparisons relying on published consensus genomes of MDV, we obtained two additional consensus genomes of vaccine strain CVI988 (Rispens) and two additional consensus genomes of the very virulent strain Md5 by sequencing viral stocks and cultured field isolates. In conjunction with the published genomes of CVI988 and Md5, this allowed us to perform three-way comparisons between multiple consensus genomes of the same strain. We found that consensus genomes of CVI988 can vary in as many as 236 positions involving 13 open reading frames (ORFs). By contrast, we found that Md5 genomes varied only in 11 positions involving a single ORF. Notably, we were able to identify 3 single-nucleotide polymorphisms (SNPs) in the unique long region and 16 SNPs in the unique short (US) region of CVI988GenBank.BAC that were not present in either CVI988Pirbright.lab or CVI988USDA.PA.field. Recombination analyses of field strains previously described as natural recombinants of CVI988 yielded no evidence of crossover events in the US region when either CVI988Pirbright.lab or CVI988USDA.PA.field were used to represent CVI988 instead of CVI988GenBank.BAC. We were also able to confirm that both CVI988 and Md5 populations were mixed, exhibiting a total of 29 and 27 high-confidence minor variant positions, respectively. However, we did not find any evidence of minor variants in the positions corresponding to the 19 SNPs in the unique regions of CVI988GenBank.BAC. Taken together, our findings suggest that continued reliance on the same published consensus genome of CVI988 may have led to an overestimation of genomic divergence between CVI988 and virulent strains and that multiple consensus genomes per strain may be necessary to ensure the accuracy of interstrain genomic comparisons.

Whole Genome Sequence-Based Analysis of Bovine Gammaherpesvirus 4 Isolated from Bovine Abortions

Bovine gammaherpesvirus 4 (BoGHV4) is a member of the Gammaherspivirinae subfamily, Rhadinovirus genus. Its natural host is the bovine, and it is prevalent among the global cattle population. Although the complete genome of BoGHV4 has been successfully sequenced, the functions of most of its genes remain unknown. Currently, only six strains of BoGHV4, all belonging to Genotype 1, have been sequenced. This is the first report of the nearly complete genome of Argentinean BoGHV4 strains isolated from clinical cases of abortion, representing the first BoGHV4 Genotype 2 and 3 genomes described in the literature. Both Argentinean isolates presented the highest nt p-distance values, indicating a greater level of divergence. Overall, the considerable diversity observed in the complete genomes and open reading frames underscores the distinctiveness of both Argentinean isolates compared to the existing BoGHV4 genomes. These findings support previous studies that categorized the Argentinean BoGHV4 strains 07-435 and 10-154 as Genotypes 3 and 2, respectively. The inclusion of these sequences represents a significant expansion to the currently limited pool of BoGHV4 genomes while providing an important basis to increase the knowledge of local isolates.

Lytic and Latent Genetic Diversity of the Epstein-Barr Virus Reveals Raji-Related Variants from Southeastern Brazil Associated with Recombination Markers

Epstein-Barr virus (EBV) is a ubiquitous gammaherpesvirus etiologically associated with benign and malignant diseases. Since the pathogenic mechanisms of EBV are not fully understood, understanding EBV genetic diversity is an ongoing goal. Therefore, the present work describes the genetic diversity of the lytic gene BZLF1 in a sampling of 70 EBV-positive cases from southeastern Brazil. Additionally, together with the genetic regions previously characterized, the aim of the present study was to determine the impact of viral genetic factors that may influence EBV genetic diversity. Accordingly, the phylogenetic analysis of the BZLF1 indicated two main clades with high support, BZ-A and BZ-B (PP > 0.85). Thus, the BZ-A clade was the most diverse clade associated with the main polymorphisms investigated, including the haplotype Type 1 + V3 (p < 0.001). Furthermore, the multigene phylogenetic analysis (MLA) between BZLF1 and the oncogene LMP1 showed specific clusters, revealing haplotypic segregation that previous single-gene phylogenies from both genes failed to demonstrate. Surprisingly, the LMP1 Raji-related variant clusters were shown to be more diverse, associated with BZ-A/B and the Type 2/1 + V3 haplotypes. Finally, due to the high haplotypic diversity of the Raji-related variants, the number of DNA recombination-inducing motifs (DRIMs) was evaluated within the different clusters defined by the MLA. Similarly, the haplotype BZ-A + Raji was shown to harbor a greater number of DRIMs (p < 0.001). These results call attention to the high haplotype diversity of EBV in southeast Brazil and strengthen the hypothesis of the recombinant potential of South American Raji-related variants via the LMP1 oncogene.

Isolation of a recombinant simian adenovirus encoding the human adenovirus G52 hexon suggests a simian origin for human adenovirus G52

Human adenoviruses (HAdVs) are causative agents of morbidity and mortality throughout the world. These double-stranded DNA viruses are phylogenetically classified into seven different species (A-G). HAdV-G52, originally isolated in 2008 from a patient presenting with gastroenteritis, is the sole human-derived member of species G. Phylogenetic analysis previously suggested that HAdV-G52 may have a simian origin, indicating a potential zoonotic spillover into humans. However, evidence of HAdV-G52 in either human or simian populations has not been reported since. Here, we describe the isolation and in vitro characterization of rhesus (rh)AdV-69, a novel simian AdV with clear evidence of recombination with HAdV-G52, from the stool of a rhesus macaque. Specifically, the rhAdV-69 hexon capsid protein is 100% identical to that of HAdV-G52, whereas the remainder of the genome is most similar to rhAdV-55, sharing 95.36% nucleic acid identity. A second recombination event with an unknown adenovirus (AdV) is evident at the short fiber gene. From the same sample, we also isolated a second, highly related recombinant AdV (rhAdV-68) that harbors a distinct hexon gene but nearly identical backbone compared to rhAdV-69. In vitro, rhAdV-68 and rhAdV-69 demonstrate comparable growth kinetics and tropisms in human cell lines, nonhuman cell lines, and human enteroids. Furthermore, we show that coinfection of highly related AdVs is not unique to this sample since we also isolated coinfecting rhAdVs from two additional rhesus macaque stool samples. Our data collectively contribute to elucidating the origins of HAdV-G52 and provide insights into the frequency of coinfections and subsequent recombination in AdV evolution.IMPORTANCEUnderstanding the host origins of adenoviruses (AdVs) is critical for public health as transmission of viruses from animals to humans can lead to emergent viruses. Recombination between animal and human AdVs can also produce emergent viruses. HAdV-G52 is the only human-derived member of the HAdV G species. It has been suggested that HAdV-G52 has a simian origin. Here, we isolated from a rhesus macaque, a novel rhAdV, rhAdV-69, that encodes a hexon protein that is 100% identical to that of HAdV-G52. This observation suggests that HAdV-G52 may indeed have a simian origin. We also isolated a highly related rhAdV, differing only in the hexon gene, from the same rhesus macaque stool sample as rhAdV-69, illustrating the potential for co-infection of closely related AdVs and recombination at the hexon gene. Furthermore, our study highlights the critical role of whole-genome sequencing in understanding AdV evolution and monitoring the emergence of pathogenic AdVs.

Single-genome analysis reveals a heterogeneous association of the herpes simplex virus genome with H3K27me2 and the reader PHF20L1 following infection of human fibroblasts

The fate of herpesvirus genomes following entry into different cell types is thought to regulate the outcome of infection. For the Herpes simplex virus 1 (HSV-1), latent infection of neurons is characterized by association with repressive heterochromatin marked with Polycomb silencing-associated lysine 27 methylation on histone H3 (H3K27me). However, whether H3K27 methylation plays a role in repressing lytic gene expression in non-neuronal cells is unclear. To address this gap in knowledge, and with consideration that the fate of the viral genome and outcome of HSV-1 infection could be heterogeneous, we developed an assay to quantify the abundance of histone modifications within single viral genome foci of infected fibroblasts. Using this approach, combined with bulk epigenetic techniques, we were unable to detect any role for H3K27me3 during HSV-1 lytic infection of fibroblasts. By contrast, we could detect the lesser studied H3K27me2 on a subpopulation of viral genomes, which was consistent with a role for H3K27 demethylases in promoting lytic gene expression. In addition, viral genomes co-localized with the H3K27me2 reader protein PHF20L1, and this association was enhanced by inhibition of the H3K27 demethylases UTX and JMJD3. Notably, targeting of H3K27me2 to viral genomes was enhanced following infection with a transcriptionally defective virus in the absence of Promyelocytic leukemia nuclear bodies. Collectively, these studies implicate a role for H3K27me2 in fibroblast-associated HSV genome silencing in a manner dependent on genome sub-nuclear localization and transcriptional activity.

IMPORTANCE

Investigating the potential mechanisms of gene silencing for DNA viruses in different cell types is important to understand the differential outcomes of infection, particularly for viruses like herpesviruses that can undergo distinct types of infection in different cell types. In addition, investigating chromatin association with viral genomes informs on the mechanisms of epigenetic regulation of DNA processes. However, there is a growing appreciation for heterogeneity in the outcome of infection at the single cell, and even single viral genome, level. Here we describe a novel assay for quantifying viral genome foci with chromatin proteins and show that a portion of genomes are targeted for silencing by H3K27me2 and associate with the reader protein PHF20L1. This study raises important questions regarding the mechanism of H3K27me2-specific targeting to viral genomes, the contribution of epigenetic heterogeneity to herpesvirus infection, and the role of PHF20L1 in regulating the outcome of DNA virus infection.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Comparative Analysis of Multiple Consensus Genomes of the Same Strain of Marek's Disease Virus Reveals Intrastrain Variation

Current strategies to understand the molecular basis of Marek's disease virus (MDV) virulence primarily consist of cataloguing divergent nucleotides between strains with different phenotypes. However, each MDV strain is typically represented by a single consensus genome despite the confirmed existence of mixed viral populations. To assess the reliability of single-consensus interstrain genomic comparisons, we obtained two additional consensus genomes of vaccine strain CVI988 (Rispens) and two additional consensus genomes of the very virulent strain Md5 by sequencing viral stocks and cultured field isolates. In conjunction with the published genomes of CVI988 and Md5, this allowed us to perform 3-way comparisons between consensus genomes of the same strain. We found that consensus genomes of CVI988 can vary in as many as 236 positions involving 13 open reading frames (ORFs). In contrast, we found that Md5 genomes varied only in 11 positions involving a single ORF. Phylogenomic analyses showed all three Md5 consensus genomes clustered closely together, while also showing that CVI988 GenBank.BAC diverged from CVI988 Pirbright.lab and CVI988 USDA.PA.field . Comparison of CVI988 consensus genomes revealed 19 SNPs in the unique regions of CVI988 GenBank.BAC that were not present in either CVI988 Pirbright.lab or CVI988 USDA.PA.field . Finally, we evaluated the genomic heterogeneity of CVI988 and Md5 populations by identifying positions with >2% read support for alternative alleles in two ultra-deeply sequenced samples. We were able to confirm that both populations of CVI988 and Md5 were mixed, exhibiting a total of 29 and 27 high-confidence minor variant positions, respectively. We did not find any evidence of minor variants in the positions corresponding to the 19 SNPs in the unique regions of CVI988 GenBank.BAC . Taken together, our findings confirm that consensus genomes of the same strain of MDV can vary and suggest that multiple consensus genomes per strain are needed in order to maximize the accuracy of interstrain genomic comparisons.

Single-genome analysis reveals heterogeneous association of the Herpes Simplex Virus genome with H3K27me2 and the reader PHF20L1 following infection of human fibroblasts

The fate of herpesvirus genomes following entry into different cell types is thought to regulate the outcome of infection. For the Herpes simplex virus 1 (HSV-1), latent infection of neurons is characterized by association with repressive heterochromatin marked with Polycomb silencing-associated lysine 27 methylation on histone H3 (H3K27me). However, whether H3K27 methylation plays a role in repressing lytic gene expression in non-neuronal cells is unclear. To address this gap in knowledge, and with consideration that the fate of the viral genome and outcome of HSV-1 infection could be heterogeneous, we developed an assay to quantify the abundance of histone modifications within single viral genome foci of infected fibroblasts. Using this approach, combined with bulk epigenetic techniques, we were unable to detect any role for H3K27me3 during HSV-1 lytic infection of fibroblasts. In contrast, we could detect the lesser studied H3K27me2 on a subpopulation of viral genomes, which was consistent with a role for H3K27 demethylases in promoting lytic gene expression. This was consistent with a role for H3K27 demethylases in promoting lytic gene expression. In addition, viral genomes co-localized with the H3K27me2 reader protein PHF20L1, and this association was enhanced by inhibition of the H3K27 demethylases UTX and JMJD3. Notably, targeting of H3K27me2 to viral genomes was enhanced following infection with a transcriptionally defective virus in the absence of Promyelocytic leukemia nuclear bodies. Collectively, these studies implicate a role for H3K27me2 in fibroblast-associated HSV genome silencing in a manner dependent on genome sub-nuclear localization and transcriptional activity.

Viral Diversity and Epidemiology in Critically Endangered Yangtze Finless Porpoises (Neophocaena asiaeorientalis asiaeorientalis)

The Yangtze finless porpoise (YFP) (Neophocaena asiaeorientalis asiaeorientalis) is a critically endangered freshwater cetacean, with about 1,249 individuals thought to be left in the wild. However, viral entities and viral diseases of YFPs remain obscure. In this study, anal swabs for virome analysis were collected during the physical examination of YFPs in the Tian-E-Zhou Oxbow (TEO) ex situ reserve. A total of 19 eukaryotic viral species belonging to 9 families, including Papillomaviridae, Herpesviridae, Picornaviridae, Picobirnaviridae, Caliciviridae, Retroviridae, Parvoviridae, Virgaviridae, and Narnaviridae, and other unclassified viruses were identified based on metasequencing. Among these detected viruses, a novel herpesvirus (NaHV), two different kobuviruses (NaKV1-2), and six different papillomaviruses (NaPV1 to -6) were considered potential risks to YFPs and confirmed by PCR or reverse transcription-PCR (RT-PCR). Most YFPs sampled were found to harbor one or more kinds of detected viral genomes (52/58 [89.7%]). Surveillance results demonstrated that kobuvirus and herpesvirus displayed obvious age distribution and PVs showed significant gender difference in YFPs. According to species demarcation criteria in individual genera in Papillomaviridae, two novel species (referred to as Omikronpapillomavirus 2 and 3) and four novel isolates of PV were identified in YFPs. Further evolutionary analysis suggested that NaPVs would occupy the mucosal niche and that virus-host codivergence mixed with duplications and host-switching events drives the evolution of cetacean PVs. Divergence times of PVs in YFP and other cetacean reflect the incipient speciation of YFPs. In summary, our findings revealed the potential viral entities, their prevalence, and their evolutionary history in YFPs, which raises an important issue regarding effects of viral infection on the fitness of YFPs. IMPORTANCE The Yangtze finless porpoise (YFP) is the only cetacean species in freshwater following the functional extinction of the baiji (Lipotes vexillifer). Health management, disease treatment, and other special measures are important for maintaining the existing YFP populations, especially in in situ and ex situ reserves. The discovery of potential viral entities and their prevalence in YFPs raises an important issue regarding the effects of viral infection on the fitness of YFPs and may contribute to the conservation of YFPs. The evolutionary history of papillomaviruses in YFP and other cetaceans reflects the phylogeny of their hosts and supports the status of incipient species, opening a window to investigate the evolutionary adaptation of cetaceans to freshwater as well as their phylogeny to remedy the deficiency of fossil evidence.

High-throughput discovery and characterization of viral transcriptional effectors in human cells

Viruses encode transcriptional regulatory proteins critical for controlling viral and host gene expression. Given their multifunctional nature and high sequence divergence, it is unclear which viral proteins can affect transcription and which specific sequences contribute to this function. Using a high-throughput assay, we measured the transcriptional regulatory potential of over 60,000 protein tiles across ∼1,500 proteins from 11 coronaviruses and all nine human herpesviruses. We discovered hundreds of transcriptional effector domains, including a conserved repression domain in all coronavirus Spike homologs, dual activation-repression domains in viral interferon regulatory factors (VIRFs), and an activation domain in six herpesvirus homologs of the single-stranded DNA-binding protein that we show is important for viral replication and late gene expression in Kaposi's sarcoma-associated herpesvirus (KSHV). For the effector domains we identified, we investigated their mechanisms via high-throughput sequence and chemical perturbations, pinpointing sequence motifs essential for function. This work massively expands viral protein annotations, serving as a springboard for studying their biological and health implications and providing new candidates for compact gene regulation tools.

Surveillance for feline herpesvirus type 1 mutation and development of resistance in cats treated with antiviral medications

Feline herpesvirus type 1 (FHV-1) commonly causes ocular surface disease in cats and is treated with antiviral medications targeting viral DNA polymerase (UL30/42). Herein, we describe a method to assess the FHV-1 genome for mutation development and to assess the functional impact of mutations, if present. Fourteen shelter-housed domestic cats with FHV-1 ocular surface disease were assigned to one of four treatment groups: placebo (n = 3), cidofovir 0.5% ophthalmic solution (n = 3), famciclovir oral solution (n = 5), or ganciclovir 0.15% ophthalmic solution (n = 3). Swabs were collected before (day 1) and after (day 8) 1 week of twice-daily treatments to isolate viable FHV-1. Viral DNA was extracted for sequencing using Illumina MiSeq with subsequent genomic variant detection between paired day 1 and day 8 isolates. Plaque reduction assay was performed on paired isolates demonstrating non-synonymous variants. A total of 171 synonymous and 3 non-synonymous variants were identified in day 8 isolates. No variants were detected in viral UL23, UL30, or UL42 genes. Variant totals were not statistically different in animals receiving antiviral or placebo (p = 0.4997). A day 8 isolate from each antiviral treatment group contained a single non-synonymous variant in ICP4 (transcriptional regulator). These 3 isolates demonstrated no evidence of functional antiviral resistance when IC₅₀ was assessed. Most (10/14 pairs) day 1 and 8 viral isolate pairs from the same host animal were near-identical. While functional variants were not detected in this small sample, these techniques can be replicated to assess FHV-1 isolates suspected of having developed resistance to antiviral medications.

Cellular and Natural Viral Engineering in Cognition-Based Evolution

Neo-Darwinism conceptualizes evolution as the continuous succession of predominately random genetic variations disciplined by natural selection. In that frame, the primary interaction between cells and the virome is relegated to host-parasite dynamics governed by selective influences. Cognition-Based Evolution regards biological and evolutionary development as a reciprocating cognition-based informational interactome for the protection of self-referential cells. To sustain cellular homeorhesis, cognitive cells collaborate to assess the validity of ambiguous biological information. That collective interaction involves coordinate measurement, communication, and active deployment of resources as Natural Cellular Engineering. These coordinated activities drive multicellularity, biological development, and evolutionary change. The virome participates as the vital intercessory among the cellular domains to ensure their shared permanent perpetuation. The interactions between the virome and the cellular domains represent active virocellular cross-communications for the continual exchange of resources. Modular genetic transfers between viruses and cells carry bioactive potentials. Those exchanges are deployed as nonrandom flexible tools among the domains in their continuous confrontation with environmental stresses. This alternative framework fundamentally shifts our perspective on viral-cellular interactions, strengthening established principles of viral symbiogenesis. Pathogenesis can now be properly appraised as one expression of a range of outcomes between cells and viruses within a larger conceptual framework of Natural Viral Engineering as a co-engineering participant with cells. It is proposed that Natural Viral Engineering should be viewed as a co-existent facet of Natural Cellular Engineering within Cognition-Based Evolution.

Characterization of Pseudorabies Virus Associated with Severe Respiratory and Neuronal Signs in Old Pigs

Pseudorabies virus (PRV) represents a leading threat to the global pig industry. Generally, pigs exhibit a pronounced age resistance against PRV, and the virus generally does not cause severe clinical signs and even death in old pigs. However, we characterized two PRV strains (HeN21 and HuB20) associated with severe respiratory and neuronal signs in old pigs. Among these two strains, HeN21 was isolated from the tonsil of a 24-week-old pig that died from severe neuronal and respiratory signs in a PRV-outbreak farm where a commercial PRV attenuated vaccine developed based on a PRV variant was used; while, HuB20 was isolated from the lung and lymph node of a 20-week-old with symptoms in another farm where Bartha-K61 vaccine was used. In vitro evaluations in different cell models demonstrated that HeN21 and HuB20 led to similar cytotoxic effects to those caused by PRV variants on PK-15, Vero, and SK-N-SH cells after 30 hours of inoculation. However, HeN21 possessed a higher titer than the other PRV variants from the first to the fifth passage on PK-15 cells and induced plaques with larger size. In vivo assessments in mouse and fattening pig models showed that inoculations of HeN21 and HuB20 caused higher morbidity and mortality and severe pathological damages in tested animals. In particular, challenge of HeN21 led to severe respiratory and neuronal signs in 90-day-old pigs. These two strains displayed higher virus loads on the main organs of challenged mice and pigs. Phylogenetic analysis revealed that HeN21 and HuB20 belonged to genotype II. In addition, recombinant events were identified in the genomes of HeN21 and HuB20, and several events were located within genes associated with PRV virulence. Our data herein may suggest the emergence of novel PRV strains in China.

Virus-induced interference as a means for accelerating fitness-based selection of cyprinid herpesvirus 3 single-nucleotide variants in vitro and in vivo

Cyprinid herpesvirus 3 (CyHV-3) is the archetype of fish alloherpesviruses and is advantageous to research because, unlike many herpesviruses, it can be studied in the laboratory by infection of the natural host (common and koi carp). Previous studies have reported a negative correlation among CyHV-3 strains between viral growth in vitro (in cell culture) and virulence in vivo (in fish). This suggests the existence of genovariants conferring enhanced fitness in vitro but reduced fitness in vivo and vice versa. Here, we identified the syncytial plaque formation in vitro as a common trait of CyHV-3 strains adapted to cell culture. A comparison of the sequences of virion transmembrane protein genes in CyHV-3 strains, and the use of various recombinant viruses, demonstrated that this trait is linked to a single-nucleotide polymorphism (SNP) in the open reading frame (ORF) 131 coding sequence (C225791T mutation) that results in codon 183 encoding either an alanine (183A) or a threonine (183T) residue. In experiments involving infections with recombinant viruses differing only by this SNP, the 183A genovariant associated with syncytial plaque formation was the more fit in vitro but the less fit in vivo. In experiments involving coinfection with both viruses, the more fit genovariant contributed to the purifying selection of the less fit genovariant by outcompeting it. In addition, this process appeared to be accelerated by viral stimulation of interference at a cellular level and stimulation of resistance to superinfection at a host level. Collectively, this study illustrates how the fundamental biological properties of some viruses and their hosts may have a profound impact on the degree of diversity that arises within viral populations.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

Performance Evaluation of the Fully Automated NeuMoDx RT-PCR Platform for the Quantification of CMV and EBV DNA in EDTA Plasma: Implications for Clinical Management and Establishment of a Conversion Formula

The NeuMoDx96 platform is a fully automated real-time PCR (RT-PCR) system. To provide continued testing quality with the introduction of new assays, the primary aim of this study was to evaluate the analytical and clinical performance of the NeuMoDx platform for the detection and quantification of CMV and EBV DNA in EDTA plasma. As no conversion from log₁₀ international units per milliliter to copies per milliliter was provided, the secondary aim was to calculate and establish a conversion factor for the output of results in copies per milliliter for CMV and EBV. Archived ETDA plasma samples (cytomegalovirus [CMV], n = 290; Ebstein-Barr virus [EBV], n = 254) were used to evaluate the analytical performance of the NeuMoDx96 platform against the routine real-time quantitative PCR (qPCR) assays. Additionally, the first WHO international standards (WHO-IS) for CMV (n = 70) and EBV (n = 72) were used for the calculation of the intra- and interassay variation. WHO-IS qualitative agreement between the assays was 100%. Intra-assay variability was low for both CMV assays (coefficient of variation [CV], phosphate-buffered saline [PBS], 3 log₁₀ IU/mL NeuMoDx, 3.67%; Abbott RealTime, CMV, 3.35%) and NeuMoDx EBV assay (CV, PBS, 3 log₁₀ IU/mL, 3.05%) but high for the Altona EBV assay (CV, PBS, 3 log₁₀ IU/mL, 26.13%). The overall qualitative concordance in clinical samples was 96.8% (270/279) for CMV and 96.7% (237/245) for EBV. The mean difference between the assays was -0.2 log₁₀ IU/mL (CMV) and -0.18 log₁₀ IU/mL (EBV). High qualitative concordance and a significant correlation of quantitative values for both assays make NeuMoDx CMV and EBV assays suitable for routine diagnostic testing. The new RT-PCR system and conversion formulas to report results in copies per milliliter are now applied in clinical routine testing. IMPORTANCE Clinical management of solid organ transplant (SOT) patients requires the careful monitoring of immunosuppression and viral infection or reactivation. qPCR is the gold standard for the detection and quantification of very small amounts of viral DNA and allows for an early assessment of viral load kinetics. The tested NeuMoDx 96 platform provides faster results than the previously used RT-PCR workflows for CMV (Abbott m2000 and RealTime CMV assay) and EBV (LightCycler 480 II, Roche high pure extraction, and Altona RealStar EBV assay) DNA detection. The implemented conversion formulas allow the continued reporting in clinically established copies per milliliter, important for long-term care of SOT patients.

Two-Color CRISPR Imaging Reveals Dynamics of Herpes Simplex Virus 1 Replication Compartments and Virus-Host Interactions

Real-time imaging tools for single-virus tracking provide spatially resolved, quantitative measurements of viral replication and virus-host interactions. However, efficiently labeling both parental and progeny viruses in living host cells remains challenging. Here, we developed a novel strategy using the CRISPR-Tag system to detect herpes simplex virus 1 (HSV-1) DNA in host cells. We created recombinant HSV-1 harboring an ~600-bp CRISPR-Tag sequence which can be sufficiently recognized by dCas9-fluorescent protein (FP) fusion proteins. CRISPR-assisted single viral genome tracking (CASVIT) allows us to assess the temporal and spatial information of viral replication at the single-cell level. Combining the advantages of SunTag and tandem split green fluorescent protein (GFP) in amplifying fluorescent signals, dSaCas9-tdTomato_10x and dSpCas9-GFP_14x were constructed to enable efficient two-color CASVIT detection. Real-time two-color imaging indicates that replication compartments (RCs) frequently come into contact with each other but do not mix, suggesting that RC territory is highly stable. Last, two-color CASVIT enables simultaneous tracking of viral DNA and host chromatin, which reveals that a dramatic loss of telomeric and centromeric DNA occurs in host cells at the early stage of viral replication. Overall, our work has established a framework for developing CRISPR-Cas9-based imaging tools to study DNA viruses in living cells. IMPORTANCE Herpes simplex virus 1 (HSV-1), a representative of the family Herpesviridae, is a ubiquitous pathogen that can establish lifelong infections and widely affects human health. Viral infection is a dynamic process that involves many steps and interactions with various cellular structures, including host chromatin. A common viral replication strategy is to form RCs that concentrate factors required for viral replication. Efficient strategies for imaging the dynamics of viral genomes, RC formation, and the interaction between the virus and host offer the opportunity to dissect the steps of the infection process and determine the mechanism underlying each step. We have developed an efficient two-color imaging system based on CRISPR-Cas9 technology to detect HSV-1 genomes quantitatively in living cells. Our results shed light on novel aspects of RC dynamics and virus-host interactions.

Evaluation of tangential flow filtration coupled to long-read sequencing for ostreid herpesvirus type 1 genome assembly

Whole-genome sequencing is widely used to better understand the transmission dynamics, the evolution and the emergence of new variants of viral pathogens. This can bring crucial information to stakeholders for disease management. Unfortunately, aquatic virus genomes are usually difficult to characterize because most of these viruses cannot be easily propagated in vitro. Developing methodologies for routine genome sequencing of aquatic viruses is timely given the ongoing threat of disease emergence. This is particularly true for pathogenic viruses infecting species of commercial interest that are widely exchanged between production basins or countries. For example, the ostreid herpesvirus type 1 (OsHV-1) is a Herpesvirus widely associated with mass mortality events of juvenile Pacific oyster Crassostrea gigas. Genomes of Herpesviruses are large and complex with long direct and inverted terminal repeats. In addition, OsHV-1 is unculturable. It therefore accumulates several features that make its genome sequencing and assembly challenging. To overcome these difficulties, we developed a tangential flow filtration (TFF) method to enrich OsHV-1 infective particles from infected host tissues. This virus purification allowed us to extract high molecular weight and high-quality viral DNA that was subjected to Illumina short-read and Nanopore long-read sequencing. Dedicated bioinformatic pipelines were developed to assemble complete OsHV-1 genomes with reads from both sequencing technologies. Nanopore sequencing allowed characterization of new structural variations and major viral isomers while having 99,98 % of nucleotide identity with the Illumina assembled genome. Our study shows that TFF-based purification method, coupled with Nanopore sequencing, is a promising approach to enable in field sequencing of unculturable aquatic DNA virus.

Genetic Variability of the Monkeypox Virus Clade IIb B.1

Monkeypox is caused by a sylvatic, double-stranded DNA zoonotic virus. Since 1 January 2022, monkeypox cases have been reported to WHO from 106 Member States across six WHO regions, and as of 2 October 2022, a total of 68,900 confirmed cases, including 25 deaths, occurred. Here, by using a whole genome approach, we perform a genetic and phylodynamic survey of the monkeypox virus Clade IIb B.1, which is the lineage causing the current multi-country outbreak. Results suggest that outbreaks seem to be isolated and localized in several epidemic clusters with geographic consistency. Currently, monkeypox appears to be a virus with a flattened genetic variability in terms of evolutionary path, with a very slow rate of growth in the population size. This scenario confirms that the monkeypox virus lacks the evolutionary advantage, given by the high level of mutation rate, which is very strong in RNA viruses. Of course, constant genome-based monitoring must be performed over time in order to detect the change in its genetic composition, if any.

Multi-phenotype analysis for enhanced classification of 11 herpes simplex virus 1 strains

Herpes simplex virus 1 (HSV1) is best known for causing oral lesions and mild clinical symptoms, but it can produce a significant range of disease severities and rates of reactivation. To better understand this phenotypic variation, we characterized 11 HSV1 strains that were isolated from individuals with diverse infection outcomes. We provide new data on genomic and in vitro plaque phenotype analysis for these isolates and compare these data to previously reported quantitation of the disease phenotype of each strain in a murine animal model. We show that integration of these three types of data permitted clustering of these HSV1 strains into four groups that were not distinguishable by any single dataset alone, highlighting the benefits of combinatorial multi-parameter phenotyping. Two strains (group 1) produced a partially or largely syncytial plaque phenotype and attenuated disease phenotypes in mice. Three strains of intermediate plaque size, causing severe disease in mice, were genetically clustered to a second group (group 2). Six strains with the smallest average plaque sizes were separated into two subgroups (groups 3 and 4) based on their different genetic clustering and disease severity in mice. Comparative genomics and network graph analysis suggested a separation of HSV1 isolates with attenuated vs. virulent phenotypes. These observations imply that virulence phenotypes of these strains may be traceable to genetic variation within the HSV1 population.

Whole genome sequence analysis of equid gammaherpesvirus -2 field isolates reveals high levels of genomic diversity and recombination

Background

Equid gammaherpesvirus 2 (EHV2) is a gammaherpesvirus with a widespread distribution in horse populations globally. Although its pathogenic significance can be unclear in most cases of infection, EHV2 infection can cause upper respiratory tract disease in foals. Co-infection of different strains of EHV2 in an individual horse is common. Small regions of the EHV2 genome have shown considerable genetic heterogeneity. This could suggest genomic recombination between different strains of EHV2, similar to the extensive recombination networks that have been demonstrated for some alphaherpesviruses. This study examined natural recombination and genome diversity of EHV2 field isolates.

Results

Whole genome sequencing analysis of 18 EHV2 isolates, along with analysis of two publicly available EHV2 genomes, revealed variation in genomes sizes (from 173.7 to 184.8 kbp), guanine plus cytosine content (from 56.7 to 57.8%) and the size of the terminal repeat regions (from 17,196 to 17,551 bp). The nucleotide sequence identity between the genomes ranged from 86.2 to 99.7%. The estimated average inter-strain nucleotide diversity between the 20 EHV2 genomes was 2.9%. Individual gene sequences showed varying levels of nucleotide diversity and ranged between 0 and 38.1%. The ratio of nonsynonymous substitutions, Ka, to synonymous substitutions, Ks, (Ka/Ks) suggests that over 50% of EHV2 genes are undergoing diversifying selection. Recombination analyses of the 20 EHV2 genome sequences using the recombination detection program (RDP4) and SplitsTree revealed evidence of viral recombination.

Conclusions

Analysis of the 18 new EHV2 genomes alongside the 2 previously sequenced genomes revealed a high degree of genetic diversity and extensive recombination networks. Herpesvirus genome diversification and virus evolution can be driven by recombination, and our findings are consistent with recombination being a key mechanism by which EHV2 genomes may vary and evolve.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08789-x.

Herpes Simplex Virus 1 Entry Glycoproteins Form Complexes before and during Membrane Fusion

Herpesviruses—ubiquitous pathogens that cause persistent infections—have some of the most complex cell entry mechanisms. Entry of the prototypical herpes simplex virus 1 (HSV-1) requires coordinated efforts of 4 glycoproteins, gB, gD, gH, and gL. The current model posits that the glycoproteins do not interact before receptor engagement and that binding of gD to its receptor causes a “cascade” of sequential pairwise interactions, first activating the gH/gL complex and subsequently activating gB, the viral fusogen. But how these glycoproteins interact remains unresolved. Here, using a quantitative split-luciferase approach, we show that pairwise HSV-1 glycoprotein complexes form before fusion, interact at a steady level throughout fusion, and do not depend on the presence of the cellular receptor. Based on our findings, we propose a revised “conformational cascade” model of HSV-1 entry. We hypothesize that all 4 glycoproteins assemble into a complex before fusion, with gH/gL positioned between gD and gB. Once gD binds to a cognate receptor, the proximity of the glycoproteins within this complex allows for efficient transmission of the activating signal from the receptor-activated gD to gH/gL to gB through sequential conformational changes, ultimately triggering the fusogenic refolding of gB. Our results also highlight previously unappreciated contributions of the transmembrane and cytoplasmic domains to glycoprotein interactions and fusion. Similar principles could be at play in other multicomponent viral entry systems, and the split-luciferase approach used here is a powerful tool for investigating protein-protein interactions in these and a variety of other systems.

A Single Herpes Simplex Virus 1 Genome Reactivates from Individual Cells

Latent infection is a characteristic feature of herpesviruses’ life cycle. Herpes simplex virus 1 is a common human pathogen that establishes lifelong latency in peripheral neurons. Symptomatic or asymptomatic periodic reactivations from the latent state allow the virus to replicate and spread among individuals. The latent viral genomes are found as several quiescent episomes inside the infected nuclei; however, it is not clear if and how many latent genomes are able to reactivate together. To address this question, we developed a quiescent infection assay, which provides a quantitative analysis of the number of genomes reactivating per cell, in cultured immortalized fibroblasts. We found that, almost always, only one viral genome reactivates per cell. We showed that different timing of entry to quiescence did not result in a significant change in the probability of reactivating. Reactivation from this quiescent state allowed only limited intergenomic recombination between two viral strains compared to lytic infection. Following coinfection with a mutant that is unable to reactivate, only coreactivation with a reactivation-proficient recombinant can provide the opportunity for the mutant to reactivate. We speculate that each individual quiescent viral genome has a low and stochastic chance to reactivate in each cell, an assumption that can explain the limited number of genomes reactivating per cell.

IMPORTANCE Herpesviruses are highly prevalent and cause significant morbidity in the human and animal populations. Most individuals who are infected with herpes simplex virus (HSV-1), a common human pathogen, will become lifelong carriers of the virus, as HSV-1 establishes latent (quiescent) infections in the host cells. Reactivation from the latent state leads to many of the viral symptoms and to the spread of the virus among individuals. While many triggers for reactivation were identified, how many genomes reactivate from an individual cell and how are these genomes selected remain understudied. Here, we identify that, in most cases, only one genome per cell reactivates. Mutated HSV-1 genomes require coinfection with another strain to allow coreactivation. Our findings suggest that the decision to reactivate is determined for each quiescent genome separately and support the notion that reactivation preferences occur at the single-genome level.

A surface pocket in the cytoplasmic domain of the herpes simplex virus fusogen gB controls membrane fusion

Membrane fusion during the entry of herpesviruses is carried out by the viral fusogen gB that is activated by its partner protein gH in some manner. The fusogenic activity of gB is controlled by its cytoplasmic (or intraviral) domain (gBCTD) and, according to the current model, the gBCTD is a trimeric, inhibitory clamp that restrains gB in the prefusion conformation. But how the gBCTD clamp is released by gH is unclear. Here, we identified two new regulatory elements within gB and gH from the prototypical herpes simplex virus 1: a surface pocket within the gBCTD and residue V831 within the gH cytoplasmic tail. Mutagenesis and structural modeling suggest that gH V831 interacts with the gB pocket. The gB pocket is located above the interface between adjacent protomers, and we hypothesize that insertion of the gH V831 wedge into the pocket serves to push the protomers apart, which releases the inhibitory clamp. In this manner, gH activates the fusogenic activity of gB. Both gB and gH are conserved across all herpesviruses, and this activation mechanism could be used by other gB homologs. Our proposed mechanism emphasizes a central role for the cytoplasmic regions in regulating the activity of a viral fusogen.

Comparison of herpes simplex virus 1 genomic diversity between adult sexual transmission partners with genital infection

Herpes simplex virus (HSV) causes chronic infection in the human host, characterized by self-limited episodes of mucosal shedding and lesional disease, with latent infection of neuronal ganglia. The epidemiology of genital herpes has undergone a significant transformation over the past two decades, with the emergence of HSV-1 as a leading cause of first-episode genital herpes in many countries. Though dsDNA viruses are not expected to mutate quickly, it is not yet known to what degree the HSV-1 viral population in a natural host adapts over time, or how often viral population variants are transmitted between hosts. This study provides a comparative genomics analysis for 33 temporally-sampled oral and genital HSV-1 genomes derived from five adult sexual transmission pairs. We found that transmission pairs harbored consensus-level viral genomes with near-complete conservation of nucleotide identity. Examination of within-host minor variants in the viral population revealed both shared and unique patterns of genetic diversity between partners, and between anatomical niches. Additionally, genetic drift was detected from spatiotemporally separated samples in as little as three days. These data expand our prior understanding of the complex interaction between HSV-1 genomics and population dynamics after transmission to new infected persons.

Genetic diversity and connectivity of the Ostreid herpesvirus 1 populations in France: A first attempt to phylogeographic inference for a marine mollusc disease

The genetic diversity of viral populations is a key driver of the spatial and temporal diffusion of viruses; yet, studying the diversity of whole genomes from natural populations still remains a challenge. Phylodynamic approaches are commonly used for RNA viruses harboring small genomes but have only rarely been applied to DNA viruses with larger genomes. Here, we used the Pacific oyster mortality syndrome (a disease that affects oyster farms around the world) as a model to study the genetic diversity of its causative agent, the Ostreid herpesvirus 1 (OsHV-1) in the three main French oyster-farming areas. Using ultra-deep sequencing on individual moribund oysters and an innovative combination of bioinformatics tools, we de novo assembled twenty-one OsHV-1 new genomes. Combining quantification of major and minor genetic variations, phylogenetic analysis, and ancestral state reconstruction of discrete traits approaches, we assessed the connectivity of OsHV-1 viral populations between the three oyster-farming areas. Our results suggest that the Marennes-Oléron Bay represents the main source of OsHV-1 diversity, from where the virus has dispersed to other farming areas, a scenario consistent with current practices of oyster transfers in France. We demonstrate that phylodynamic approaches can be applied to aquatic DNA viruses to determine how epidemiological, immunological, and evolutionary processes act and potentially interact to shape their diversity patterns.

Fraternal Twins: The Enigmatic Role of the Immune System in Alphaherpesvirus Pathogenesis and Latency and Its Impacts on Vaccine Efficacy

Although the establishment, maintenance and reactivation from alphaherpesvirus latency is far from fully understood, some things are now manifestly clear: Alphaherpesvirus latency occurs in neurons of the peripheral nervous system and control of the process is multifactorial and complex. This includes components of the immune system, contributions from non-neuronal cells surrounding neurons in ganglia, specialized nucleic acids and modifications to the viral DNA to name some of the most important. Efficacious vaccines have been developed to control both acute varicella and zoster, the outcome of reactivation, but despite considerable effort vaccines for acute herpes simplex virus (HSV) infection or reactivated lesions have thus far failed to materialize despite considerable effort. Given the relevance of the immune system to establish and maintain HSV latency, a vaccine designed to tailor the HSV response to maximize the activity of components most critical for controlling reactivated infection might limit the severity of recurrences and hence reduce viral transmission. In this review, we discuss the current understanding of immunological factors that contribute to HSV and VZV latency, identify differences between varicella-zoster virus (VZV) and HSV that could explain why vaccines have been valuable at controlling VZV disease but not HSV, and finish by outlining possible strategies for developing effective HSV vaccines.

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

Rapid typing of infectious laryngotracheitis virus directly from tracheal tissues based on next-generation sequencing

Infectious laryngotracheitis virus (ILTV) is the causative agent of an economically important disease of chickens causing upper respiratory tract infection. Strains of ILTV are commonly identified by polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) and/or PCR high resolution melt (PCR-HRM) curve analysis targeting several genes. However, these techniques examine only a limited number of mutations present inside the target regions and may generate unreliable results when the sample contains more than one strain. Here, we attempted to sequence the whole genome of ILTV with known identity (class 9) directly from tracheal scrapings to circumvent in vitro culturing, which can potentially introduce variations into the genome. Despite the large number of quality reads, mapping was compromised by poor overlapping and gaps, and assembly of the complete genome sequence was not possible. In a map-to-reference alignment, the regions with low coverage were deleted, those with high coverage were concatenated and a genome sequence of 139,465 bp was obtained, which covered 91% of the ILTV genome. Sixteen single-nucleotide polymorphisms (SNPs) were found between the ILTV isolate examined and ILTV class 9 (JN804827). Despite only 91% genome coverage, using sequence analysis and comparison with previously sequenced ILTVs, we were able to classify the isolate as class 9. Therefore, this technique has the potential to replace the current PCR-HRM technique, as it provides detailed information about the ILTV isolates.

HSV-1 DNA Replication-Coordinated Regulation by Viral and Cellular Factors

DNA replication is an integral step in the herpes simplex virus type 1 (HSV-1) life cycle that is coordinated with the cellular DNA damage response, repair and recombination of the viral genome, and viral gene transcription. HSV-1 encodes its own DNA replication machinery, including an origin binding protein (UL9), single-stranded DNA binding protein (ICP8), DNA polymerase (UL30), processivity factor (UL42), and a helicase/primase complex (UL5/UL8/UL52). In addition, HSV-1 utilizes a combination of accessory viral and cellular factors to coordinate viral DNA replication with other viral and cellular processes. The purpose of this review is to outline the roles of viral and cellular proteins in HSV-1 DNA replication and replication-coupled processes, and to highlight how HSV-1 may modify and adapt cellular proteins to facilitate productive infection.

Herpes simplex virus 2 (HSV-2) evolves faster in cell culture than HSV-1 by generating greater genetic diversity

Herpes simplex virus type 1 and 2 (HSV-1 and HSV-2, respectively) are prevalent human pathogens of clinical relevance that establish long-life latency in the nervous system. They have been considered, along with the Herpesviridae family, to exhibit a low level of genetic diversity during viral replication. However, the high ability shown by these viruses to rapidly evolve under different selective pressures does not correlates with that presumed genetic stability. High-throughput sequencing has revealed that heterogeneous or plaque-purified populations of both serotypes contain a broad range of genetic diversity, in terms of number and frequency of minor genetic variants, both in vivo and in vitro. This is reminiscent of the quasispecies phenomenon traditionally associated with RNA viruses. Here, by plaque-purification of two selected viral clones of each viral subtype, we reduced the high level of genetic variability found in the original viral stocks, to more genetically homogeneous populations. After having deeply characterized the genetic diversity present in the purified viral clones as a high confidence baseline, we examined the generation of de novo genetic diversity under culture conditions. We found that both serotypes gradually increased the number of de novo minor variants, as well as their frequency, in two different cell types after just five and ten passages. Remarkably, HSV-2 populations displayed a much higher raise of nonconservative de novo minor variants than the HSV-1 counterparts. Most of these minor variants exhibited a very low frequency in the population, increasing their frequency over sequential passages. These new appeared minor variants largely impacted the coding diversity of HSV-2, and we found some genes more prone to harbor higher variability. These data show that herpesviruses generate de novo genetic diversity differentially under equal in vitro culture conditions. This might have contributed to the evolutionary divergence of HSV-1 and HSV-2 adapting to different anatomical niche, boosted by selective pressures found at each epithelial and neuronal tissue.

Herpesviruses are highly human pathogens that establish latency in neurons of the peripheral nervous system. Colonization of nerve endings is required for herpes simplex virus (HSV) persistence and pathogenesis. HSV-1 global prevalence is much higher than HSV-2, in addition to their preferential tendency to infect the oronasal and genital areas, respectively. How these closely related viruses have been adapting and evolving to replicate and colonize these two different anatomical areas remains unclear. Herpesviruses were presumed to mutate much less than viruses with RNA genomes, due to the higher fidelity of the DNA polymerase and proofreading mechanisms when replicating. However, the worldwide accessibility and development of high-throughput sequencing technologies have revealed the heterogenicity and high diversity present in viral populations clinically isolated. Here we show that HSV-2 mutates much faster than HSV-1, when compared under similar and controlled cell culture conditions. This high mutation rate is translated into an increase in coding diversity, since the great majority of these new mutations lead to nonconservative changes in viral proteins. Understanding how herpesviruses differentially mutate under similar selective pressures is critical to prevent resistance to anti-viral drugs.

Genomic Diversity of the Ostreid Herpesvirus Type 1 Across Time and Location and Among Host Species

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. This is particularly true for pathogens with low per-site mutation rates, such as DNA viruses, that do not exhibit a large amount of evolutionary change among genetic sequences sampled at different time points. However, whole-genome sequencing can reveal the accumulation of novel genetic variation between samples, promising to render most, if not all, microbial pathogens measurably evolving and suitable for analytical techniques derived from population genetic theory. Here, we aim to assess the measurability of evolution on epidemiological time scales of the Ostreid herpesvirus 1 (OsHV-1), a double stranded DNA virus of which a new variant, OsHV-1 μVar, emerged in France in 2008, spreading across Europe and causing dramatic economic and ecological damage. We performed phylogenetic analyses of heterochronous (n = 21) OsHV-1 genomes sampled worldwide. Results show sufficient temporal signal in the viral sequences to proceed with phylogenetic molecular clock analyses and they indicate that the genetic diversity seen in these OsHV-1 isolates has arisen within the past three decades. OsHV-1 samples from France and New Zealand did not cluster together suggesting a spatial structuration of the viral populations. The genome-wide study of simple and complex polymorphisms shows that specific genomic regions are deleted in several isolates or accumulate a high number of substitutions. These contrasting and non-random patterns of polymorphism suggest that some genomic regions are affected by strong selective pressures. Interestingly, we also found variant genotypes within all infected individuals. Altogether, these results provide baseline evidence that whole genome sequencing could be used to study population dynamic processes of OsHV-1, and more broadly herpesviruses.

Whole-Genome Sequence Analysis of Pseudorabies Virus Clinical Isolates from Pigs in China between 2012 and 2017 in China

Pseudorabies virus (PRV) is an economically significant swine infectious agent. A PRV outbreak took place in China in 2011 with novel virulent variants. Although the association of viral genomic variability with pathogenicity is not fully confirmed, the knowledge concerning PRV genomic diversity and evolution is still limited. Here, we sequenced 54 genomes of novel PRV variants isolated in China from 2012 to 2017. Phylogenetic analysis revealed that China strains and US/Europe strains were classified into two separate genotypes. PRV strains isolated from 2012 to 2017 in China are highly related to each other and genetically close to classic China strains such as Ea, Fa, and SC. RDP analysis revealed 23 recombination events within novel PRV variants, indicating that recombination contributes significantly to the viral evolution. The selection pressure analysis indicated that most ORFs were under evolutionary constraint, and 19 amino acid residue sites in 15 ORFs were identified under positive selection. Additionally, 37 unique mutations were identified in 19 ORFs, which distinguish the novel variants from classic strains. Overall, our study suggested that novel PRV variants might evolve from classical PRV strains through point mutation and recombination mechanisms.

A holistic perspective on herpes simplex virus (HSV) ecology and evolution

Herpes simplex viruses (HSV) cause chronic infection in humans that are characterized by periodic episodes of mucosal shedding and ulcerative disease. HSV causes millions of infections world-wide, with lifelong bouts of viral reactivation from latency in neuronal ganglia. Infected individuals experience different levels of disease severity and frequency of reactivation. There are two distantly related HSV species, with HSV-1 infections historically found most often in the oral niche and HSV-2 infections in the genital niche. Over the last two decades, HSV-1 has emerged as the leading cause of first-episode genital herpes in multiple countries. While HSV-1 has the highest level of genetic diversity among human alpha-herpesviruses, it is not yet known how quickly the HSV-1 viral population in a human host adapts over time, or if there are population bottlenecks associated with viral reactivation and/or transmission. It is also unknown how the ecological environments in which HSV infections occur influence their evolutionary trajectory, or that of co-occurring viruses and microbes. In this review, we explore how HSV accrues genetic diversity within each new infection, and yet maintains its ability to successfully infect most of the human population. A holistic examination of the ecological context of natural human infections can expand our awareness of how HSV adapts as it moves within and between human hosts, and reveal the complexity of these lifelong human-virus interactions. These insights may in turn suggest new areas of exploration for other chronic pathogens that successfully evolve and persist among their hosts.

Current Drugs to Treat Infections with Herpes Simplex Viruses-1 and -2

Herpes simplex viruses-1 and -2 (HSV-1 and -2) are two of the three human alphaherpesviruses that cause infections worldwide. Since both viruses can be acquired in the absence of visible signs and symptoms, yet still result in lifelong infection, it is imperative that we provide interventions to keep them at bay, especially in immunocompromised patients. While numerous experimental vaccines are under consideration, current intervention consists solely of antiviral chemotherapeutic agents. This review explores all of the clinically approved drugs used to prevent the worst sequelae of recurrent outbreaks by these viruses.

Combination of long- and short-read sequencing fully resolves complex repeats of herpes simplex virus 2 strain MS complete genome

Herpes simplex virus serotype 2 (HSV-2) is a ubiquitous human pathogen that causes recurrent genital infections and ulcerations. Many HSV-2 strains with different biological properties have been identified, but only the genomes of HSV-2 strains HG52, SD90e and 333 have been reported as complete and fully characterized sequences. We de novo assembled, annotated and manually curated the complete genome sequence of HSV-2 strain MS, a highly neurovirulent strain, originally isolated from a multiple sclerosis patient. We resolved both DNA ends, as well as the complex inverted repeats regions present in HSV genomes, usually undisclosed in previous published partial herpesvirus genomes, using long reads from Pacific Biosciences (PacBio) technology. Additionally, we identified isomeric genomes by determining the alternative relative orientation of unique fragments in the genome of the sequenced viral population. Illumina short-read sequencing was crucial to examine genetic variability, such as nucleotide polymorphisms, insertion/deletions and sequence determinants of strain-specific virulence factors. We used Illumina data to fix two disrupted open reading frames found in coding homopolymers after PacBio assembly. These results support the combination of long- and short-read sequencing technologies as a precise and effective approach for the accurate de novo assembly and curation of complex microbial genomes.

Multifaceted Roles of ICP22/ORF63 Proteins in the Life Cycle of Human Herpesviruses

Herpesviruses are extremely successful parasites that have evolved over millions of years to develop a variety of mechanisms to coexist with their hosts and to maintain host-to-host transmission and lifelong infection by regulating their life cycles. The life cycle of herpesviruses consists of two phases: lytic infection and latent infection. During lytic infection, active replication and the production of numerous progeny virions occur. Subsequent suppression of the host immune response leads to a lifetime latent infection of the host. During latent infection, the viral genome remains in an inactive state in the host cell to avoid host immune surveillance, but the virus can be reactivated and reenter the lytic cycle. The balance between these two phases of the herpesvirus life cycle is controlled by broad interactions among numerous viral and cellular factors. ICP22/ORF63 proteins are among these factors and are involved in transcription, nuclear budding, latency establishment, and reactivation. In this review, we summarized the various roles and complex mechanisms by which ICP22/ORF63 proteins regulate the life cycle of human herpesviruses and the complex relationships among host and viral factors. Elucidating the role and mechanism of ICP22/ORF63 in virus-host interactions will deepen our understanding of the viral life cycle. In addition, it will also help us to understand the pathogenesis of herpesvirus infections and provide new strategies for combating these infections.

Genetic changes in plaque-purified varicella vaccine strain Suduvax during in vitro propagation in cell culture

Infection by varicella-zoster virus (VZV) can be prevented by using live attenuated vaccines. VZV vaccine strains are known to evolve rapidly in vivo, however, their genetic and biological effects are not known. In this study, the plaque-purified vaccine strain Suduvax (PPS) was used to understand the genetic changes that occur during the process of propagation in in vitro cell culture. Full genome sequences of three different passages (p4, p30, and p60) of PPS were determined and compared for genetic changes. Mutations were found at 59 positions. The number of genetically polymorphic sites (GPS) and the average of minor allele frequency (MAF) at GPSs were not significantly altered after passaging in cell culture up to p60. The number of variant nucleotide positions (VNPs), wherein GPS was found in at least one passage of PPS, was 149. Overall, MAF changed by less than 5% at 52 VNPs, increased by more than 5% at 42 VNPs, and decreased by more than 5% at 55 VNPs in p60, compared with that seen in p4. More complicated patterns of changes in MAF were observed when genetic polymorphism at 149 VNPs was analyzed among the three passages. However, MAF decreased and mixed genotypes became unequivocally fixed to vaccine type in 23 vaccine-specific positions in higher passages of PPS. Plaque-purified Suduvax appeared to adapt to better replication during in vitro cell culture. Further studies with other vaccine strains and in vivo studies will help to understand the evolution of the VZV vaccine.

An MHV-68 Mutator Phenotype Mutant Virus, Confirmed by CRISPR/Cas9-Mediated Gene Editing of the Viral DNA Polymerase Gene, Shows Reduced Viral Fitness

Drug resistance studies on human γ-herpesviruses are hampered by the absence of an in vitro system that allows efficient lytic viral replication. Therefore, we employed murine γ-herpesvirus-68 (MHV-68) that efficiently replicates in vitro as a model to study the antiviral resistance of γ-herpesviruses. In this study, we investigated the mechanism of resistance to nucleoside (ganciclovir (GCV)), nucleotide (cidofovir (CDV), HPMP-5azaC, HPMPO-DAPy) and pyrophosphate (foscarnet (PFA)) analogues and the impact of these drug resistance mutations on viral fitness. Viral fitness was determined by dual infection competition assays, where MHV-68 drug-resistant viral clones competed with the wild-type virus in the absence and presence of antivirals. Using next-generation sequencing, the composition of the viral populations was determined at the time of infection and after 5 days of growth. Antiviral drug resistance selection resulted in clones harboring mutations in the viral DNA polymerase (DP), denoted Y383S^GCV, Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA, G302W+K442T^PFA, C297W^HPMPO-DAPy and C981Y^CDV. Without antiviral pressure, viral clones Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA and G302W+K442T^PFA grew equal to the wild-type virus. However, in the presence of antivirals, these mutants had a growth advantage over the wild-type virus that was moderately to very strongly correlated with antiviral resistance. The Y383S^GCV mutant was more fit than the wild-type virus with and without antivirals, except in the presence of brivudin. The C297W and C981Y changes were associated with a mutator phenotype and had a severely impaired viral fitness in the absence and presence of antivirals. The mutator phenotype caused by C297W in MHV-68 DP was validated by using a CRISPR/Cas9 genome editing approach.

Genetic diversity through social heterosis can increase virulence in RNA viral infections and cancer progression

In viral infections and cancer tumours, negative health outcomes often correlate with increasing genetic diversity. Possible evolutionary processes for such relationships include mutant lineages escaping host control or diversity, per se, creating too many immune system targets. Another possibility is social heterosis where mutations and replicative errors create clonal lineages varying in intrinsic capability for successful dispersal; improved environmental buffering; resource extraction or effective defence against immune systems. Rather than these capabilities existing in one genome, social heterosis proposes complementary synergies occur across lineages in close proximity. Diverse groups overcome host defences as interacting 'social genomes' with group genetic tool kits exceeding limited individual plasticity. To assess the possibility of social heterosis in viral infections and cancer progression, we conducted extensive literature searches for examples consistent with general and specific predictions from the social heterosis hypothesis. Numerous studies found supportive patterns in cancers across multiple tissues and in several families of RNA viruses. In viruses, social heterosis mechanisms probably result from long coevolutionary histories of competition between pathogen and host. Conversely, in cancers, social heterosis is a by-product of recent mutations. Investigating how social genomes arise and function in viral quasi-species swarms and cancer tumours may lead to new therapeutic approaches.

Mutational pressure by host APOBEC3s more strongly affects genes expressed early in the lytic phase of herpes simplex virus-1 (HSV-1) and human polyomavirus (HPyV) infection

Herpes-Simplex Virus 1 (HSV-1) infects most humans when they are young, sometimes with fatal consequences. Gene expression occurs in a temporal order upon lytic HSV-1 infection: immediate early (IE) genes are expressed, then early (E) genes, followed by late (L) genes. During this infection cycle, the HSV-1 genome has the potential for exposure to APOBEC3 (A3) proteins, a family of cytidine deaminases that cause C>U mutations on single-stranded DNA (ssDNA), often resulting in a C>T transition. We developed a computational model for the mutational pressure of A3 on the lytic cycle of HSV-1 to determine which viral kinetic gene class is most vulnerable to A3 mutations. Using in silico stochastic methods, we simulated the infectious cycle under varying intensities of A3 mutational pressure. We found that the IE and E genes are more vulnerable to A3 than L genes. We validated this model by analyzing the A3 evolutionary footprints in 25 HSV-1 isolates. We find that IE and E genes have evolved to underrepresent A3 hotspot motifs more so than L genes, consistent with greater selection pressure on IE and E genes. We extend this model to two-step infections, such as those of polyomavirus, and find that the same pattern holds for over 25 human Polyomavirus (HPyVs) genomes. Genes expressed earlier during infection are more vulnerable to mutations than those expressed later.

The Pacific Oyster Mortality Syndrome, a Polymicrobial and Multifactorial Disease: State of Knowledge and Future Directions

The Pacific oyster (Crassostreae gigas) has been introduced from Asia to numerous countries around the world during the 20th century. C. gigas is the main oyster species farmed worldwide and represents more than 98% of oyster production. The severity of disease outbreaks that affect C. gigas, which primarily impact juvenile oysters, has increased dramatically since 2008. The most prevalent disease, Pacific oyster mortality syndrome (POMS), has become panzootic and represents a threat to the oyster industry. Recently, major steps towards understanding POMS have been achieved through integrative molecular approaches. These studies demonstrated that infection by Ostreid herpesvirus type 1 µVar (OsHV-1 µvar) is the first critical step in the infectious process and leads to an immunocompromised state by altering hemocyte physiology. This is followed by dysbiosis of the microbiota, which leads to a secondary colonization by opportunistic bacterial pathogens, which in turn results in oyster death. Host and environmental factors (e.g. oyster genetics and age, temperature, food availability, and microbiota) have been shown to influence POMS permissiveness. However, we still do not understand the mechanisms by which these different factors control disease expression. The present review discusses current knowledge of this polymicrobial and multifactorial disease process and explores the research avenues that must be investigated to fully elucidate the complexity of POMS. These discoveries will help in decision-making and will facilitate the development of tools and applied innovations for the sustainable and integrated management of oyster aquaculture.

Genomic diversity in a population of Spodoptera frugiperda nucleopolyhedrovirus

Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV) represents a strong candidate to develop environmental-friendly pesticides against the fall armyworm (Spodoptera frugiperda), a widespread pest that poses a severe threat to different crops around the world. To date, SfMNPV genomic diversity of different isolates has been mainly studied by means of restriction pattern analyses and by sequencing of the egt region. Here, the genomic diversity present inside an isolate of SfMNPV was explored using high-throughput sequencing for the first time. We identified 704 intrahost single nucleotide variants, from which 184 are nonsynonymous mutations distributed among 82 different coding sequences. We detected several structural variants affecting SfMNPV genome, including two previously reported deletions inside the egt region. A comparative analysis between polymorphisms present in different SfMNPV isolates and our intraisolate diversity data suggests that coding regions with higher genetic diversity are associated with oral infectivity or unknown functions. In this context, through molecular evolution studies we provide evidence of diversifying selection acting on sf29, a putative collagenase which could contribute to the oral infectivity of SfMNPV. Overall, our results contribute to deepen our understanding of the coevolution between SfMNPV and the fall armyworm and will be useful to improve the applicability of this virus as a biological control agent.

Mechanisms of DNA Virus Evolution

The mechanisms by which DNA viruses adapt and evolve over time include minor accumulated changes associated with genetic drift – such as single nucleotide changes and small insertions or deletions – as well as more substantial changes equivalent to genetic shift. The latter case is most often driven by recombination between co-replicating genomes of the same or different species, but it can also include large segmental duplications and functionally tolerated deletions. Finally, instances of genetic exchange between virus and host or between virus species – i.e., horizontal gene transfer (HGT) – have driven the evolution of DNA viruses as well as their host organisms.

Interrogating Host Antiviral Environments Driven by Nuclear DNA Sensing: A Multiomic Perspective

Nuclear DNA sensors are critical components of the mammalian innate immune system, recognizing the presence of pathogens and initiating immune signaling. These proteins act in the nuclei of infected cells by binding to foreign DNA, such as the viral genomes of nuclear-replicating DNA viruses herpes simplex virus type 1 (HSV-1) and human cytomegalovirus (HCMV). Upon binding to pathogenic DNA, the nuclear DNA sensors were shown to initiate antiviral cytokines, as well as to suppress viral gene expression. These host defense responses involve complex signaling processes that, through protein–protein interactions (PPIs) and post-translational modifications (PTMs), drive extensive remodeling of the cellular transcriptome, proteome, and secretome to generate an antiviral environment. As such, a holistic understanding of these changes is required to understand the mechanisms through which nuclear DNA sensors act. The advent of omics techniques has revolutionized the speed and scale at which biological research is conducted and has been used to make great strides in uncovering the molecular underpinnings of DNA sensing. Here, we review the contribution of proteomics approaches to characterizing nuclear DNA sensors via the discovery of functional PPIs and PTMs, as well as proteome and secretome changes that define a host antiviral environment. We also highlight the value of and future need for integrative multiomic efforts to gain a systems-level understanding of DNA sensors and their influence on epigenetic and transcriptomic alterations during infection.

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.