Islands of linkage in an ocean of pervasive recombination reveals two-speed evolution of human cytomegalovirus genomes

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.

Inhibiting HCMV pUL89-C Endonuclease with Metal-Binding Compounds

Human cytomegalovirus (HCMV) infects individuals of all ages and establishes a lifelong latency. Current antiviral drugs are suboptimal in efficacy and safety and ineffective against resistant/refractory HCMV. Therefore, there is an unmet clinical need for efficacious, safe, and mechanistically novel HCMV drugs. The recent Food and Drug Administration (FDA) approval of letermovir (LTV) validated the HCMV terminase complex as a new target for antiviral development. LTV targets terminase subunit pUL56 but not the main endonuclease enzymatic function housed in the C terminus of subunit pUL89 (pUL89-C). Structurally and mechanistically, pUL89-C is an RNase H-like viral endonuclease entailing two divalent metal ions at the active site. In recent years, numerous studies have extensively explored pUL89-C inhibition using metal-chelating chemotypes, an approach previously used for inhibiting HIV ribonuclease H (RNase H) and integrase strand transfer (INST). Collectively, the work summarized herein validates the use of metal-binding scaffolds for designing potent and specific pUL89-C inhibitors.

Direct Nanopore Sequencing of Human Cytomegalovirus Genomes from High-Viral-Load Clinical Samples

Nanopore sequencing is becoming increasingly commonplace in clinical settings, particularly for diagnostic assessments and outbreak investigations, due to its portability, low cost, and ability to operate in near real-time. Although high sequencing error rates initially hampered the wider implementation of this technology, improvements have been made continually with each iteration of the sequencing hardware and base-calling software. Here, we assess the feasibility of using nanopore sequencing to determine the complete genomes of human cytomegalovirus (HCMV) in high-viral-load clinical samples without viral DNA enrichment, PCR amplification, or prior knowledge of the sequences. We utilised a hybrid bioinformatic approach that involved assembling the reads de novo, improving the consensus sequence by aligning reads to the best-matching genome from a collated set of published sequences, and polishing the improved consensus sequence. The final genomes from a urine sample and a lung sample, the former with an HCMV to human DNA load approximately 50 times greater than the latter, achieved 99.97 and 99.93% identity, respectively, to the benchmark genomes obtained independently by Illumina sequencing. Thus, we demonstrated that nanopore sequencing is capable of determining HCMV genomes directly from high-viral-load clinical samples with a high accuracy.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

Identifying high-confidence variants in human cytomegalovirus genomes sequenced from clinical samples

Understanding the intrahost evolution of viral populations has implications in pathogenesis, diagnosis, and treatment and has recently made impressive advances from developments in high-throughput sequencing. However, the underlying analyses are very sensitive to sources of bias, error, and artefact in the data, and it is important that these are addressed adequately if robust conclusions are to be drawn. The key factors include (1) determining the number of viral strains present in the sample analysed; (2) monitoring the extent to which the data represent these strains and assessing the quality of these data; (3) dealing with the effects of cross-contamination; and (4) ensuring that the results are reproducible. We investigated these factors by generating sequence datasets, including biological and technical replicates, directly from clinical samples obtained from a small cohort of patients who had been infected congenitally with the herpesvirus human cytomegalovirus, with the aim of developing a strategy for identifying high-confidence intrahost variants. We found that such variants were few in number and typically present in low proportions and concluded that human cytomegalovirus exhibits a very low level of intrahost variability. In addition to clarifying the situation regarding human cytomegalovirus, our strategy has wider applicability to understanding the intrahost variability of other viruses.

Recent Advancements in Understanding Primary Cytomegalovirus Infection in a Mouse Model

Animal models that mimic human infections provide insights in virus–host interplay; knowledge that in vitro approaches cannot readily predict, nor easily reproduce. Human cytomegalovirus (HCMV) infections are acquired asymptomatically, and primary infections are difficult to capture. The gap in our knowledge of the early events of HCMV colonization and spread limits rational design of HCMV antivirals and vaccines. Studies of natural infection with mouse cytomegalovirus (MCMV) have demonstrated the olfactory epithelium as the site of natural colonization. Systemic spread from the olfactory epithelium is facilitated by infected dendritic cells (DC); tracking dissemination uncovered previously unappreciated DC trafficking pathways. The olfactory epithelium also provides a unique niche that supports efficient MCMV superinfection and virus recombination. In this review, we summarize recent advances to our understanding of MCMV infection and spread and the tissue-specific mechanisms utilized by MCMV to modulate DC trafficking. As these mechanisms are likely conserved with HCMV, they may inform new approaches for preventing HCMV infections in humans.

Human cytomegalovirus strain diversity and dynamics reveal the donor lung as a major contributor after transplantation

Mixed human cytomegalovirus (HCMV) strain infections are frequent in lung transplant recipients (LTRs). To date, the influence of the donor (D) and recipient (R) HCMV serostatus on intra-host HCMV strain composition and viral population dynamics after transplantation is only poorly understood. Here, we investigated ten pre-transplant lungs from HCMV-seropositive donors and 163 sequential HCMV-DNA-positive plasma and bronchoalveolar lavage samples from fifty LTRs with multiviremic episodes post-transplantation. The study cohort included D+R+ (38 per cent), D+R- (36 per cent), and D-R+ (26 per cent) patients. All samples were subjected to quantitative genotyping by short amplicon deep sequencing, and twenty-four of them were additionally PacBio long-read sequenced for genotype linkages. We find that D+R+ patients show a significantly elevated intra-host strain diversity compared to D+R- and D-R+ patients (P = 0.0089). Both D+ patient groups display significantly higher viral population dynamics than D- patients (P = 0.0061). Five out of ten pre-transplant donor lungs were HCMV DNA positive, whereof three multiple HCMV strains were detected, indicating that multi-strain transmission via lung transplantation is likely. Using long reads, we show that intra-host haplotypes can share distinctly linked genotypes, which limits overall intra-host diversity in mixed infections. Together, our findings demonstrate donor-derived strains as the main source of increased HCMV strain diversity and dynamics post-transplantation. These results foster strategies to mitigate the potential transmission of the donor strain reservoir to the allograft, such as ex vivo delivery of HCMV-selective immunotoxins prior to transplantation to reduce latent HCMV.

Host translesion polymerases are required for viral genome integrity

Human cells encode up to 15 DNA polymerases with specialized functions in chromosomal DNA synthesis and damage repair. In contrast, complex DNA viruses, such as those of the herpesviridae family, encode a single B-family DNA polymerase. This disparity raises the possibility that DNA viruses may rely on host polymerases for synthesis through complex DNA geometries. We tested the importance of error-prone Y-family polymerases involved in translesion synthesis (TLS) to human cytomegalovirus (HCMV) infection. We find most Y-family polymerases involved in the nucleotide insertion and bypass of lesions restrict HCMV genome synthesis and replication. In contrast, other TLS polymerases, such as the polymerase ζ complex, which extends past lesions, was required for optimal genome synthesis and replication. Depletion of either the polζ complex or the suite of insertion polymerases demonstrate that TLS polymerases suppress the frequency of viral genome rearrangements, particularly at GC-rich sites and repeat sequences. Moreover, while distinct from HCMV, replication of the related herpes simplex virus type 1 is impacted by host TLS polymerases, suggesting a broader requirement for host polymerases for DNA virus replication. These findings reveal an unexpected role for host DNA polymerases in ensuring viral genome stability.

Human herpesvirus diversity is altered in HLA class I binding peptides

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

Neutralizing Antibodies Limit Cell-Associated Spread of Human Cytomegalovirus in Epithelial Cells and Fibroblasts

Human cytomegalovirus (HCMV) can cause severe clinical disease in immunocompromised individuals, such as allograft recipients and infants infected in utero. Neutralizing activity of antibodies, measured as the ability to prevent the entry of cell-free virus, has been correlated with the reduction in HCMV transmission and the severity of HCMV-associated disease. However, in vivo HCMV amplification may occur mainly via cell-to-cell spread. Thus, quantifying the inhibition of cell-to-cell transmission could be important in the evaluation of therapeutic antibodies and/or humoral responses to infection or immunization. Here, we established a quantitative plaque reduction assay, which allowed for the measurement of the capacity of antibodies to limit HCMV spread in vitro. Using an automated fluorescence spot reader, infection progression was assayed by the expansion of viral plaques during the course of infection with various GFP-expressing viruses. We found that in contrast to non-neutralizing monoclonal antibodies (mAbs), neutralizing mAbs against both glycoprotein B and H (gB and gH) could significantly inhibit viral plaque expansion of different HCMV strains and was equally efficient in fibroblasts as in epithelial cells. In contrast, an anti-pentamer mAb was active only in epithelial cells. Taken together, our data demonstrate that specific anti-HCMV mAbs can significantly limit cell-associated virus spread in vitro.

Long range PCR-based deep sequencing for haplotype determination in mixed HCMV infections

Background

Short read sequencing has been used extensively to decipher the genome diversity of human cytomegalovirus (HCMV) strains, but falls short to reveal individual genomes in mixed HCMV strain populations. Novel third-generation sequencing platforms offer an extended read length and promise to resolve how distant polymorphic sites along individual genomes are linked. In the present study, we established a long amplicon PacBio sequencing workflow to identify the absolute and relative quantities of unique HCMV haplotypes spanning over multiple hypervariable sites in mixtures. Initial validation of this approach was performed with defined HCMV DNA templates derived from cell-culture enriched viruses and was further tested for its suitability on patient samples carrying mixed HCMV infections.

Results

Total substitution and indel error rate of mapped reads ranged from 0.17 to 0.43% depending on the stringency of quality trimming. Artificial HCMV DNA mixtures were correctly determined down to 1% abundance of the minor DNA source when the total HCMV DNA input was 4 × 10⁴ copies/ml. PCR products of up to 7.7 kb and a GC content < 55% were efficiently generated when DNA was directly isolated from patient samples. In a single sample, up to three distinct haplotypes were identified showing varying relative frequencies. Alignments of distinct haplotype sequences within patient samples showed uneven distribution of sequence diversity, interspersed by long identical stretches. Moreover, diversity estimation at single polymorphic regions as assessed by short amplicon sequencing may markedly underestimate the overall diversity of mixed haplotype populations.

Conclusions

Quantitative haplotype determination by long amplicon sequencing provides a novel approach for HCMV strain characterisation in mixed infected samples which can be scaled up to cover the majority of the genome by multi-amplicon panels. This will substantially improve our understanding of intra-host HCMV strain diversity and its dynamic behaviour.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08272-z.

Optimization of a Lambda-RED Recombination Method for Rapid Gene Deletion in Human Cytomegalovirus

Human cytomegalovirus (HCMV) continues to be a major cause of morbidity in transplant patients and newborns. However, the functions of many of the more than 282 genes encoded in the HCMV genome remain unknown. The development of bacterial artificial chromosome (BAC) technology contributes to the genetic manipulation of several organisms including HCMV. The maintenance of the HCMV BAC in E. coli cells permits the rapid generation of recombinant viral genomes that can be used to produce viral progeny in cell cultures for the study of gene function. We optimized the Lambda-Red Recombination system to construct HCMV gene deletion mutants rapidly in the complete set of tested genes. This method constitutes a useful tool that allows for the quick generation of a high number of gene deletion mutants, allowing for the analysis of the whole genome to improve our understanding of HCMV gene function. This may also facilitate the development of novel vaccines and therapeutics.

Varicella-Zoster Virus-Genetics, Molecular Evolution and Recombination

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

Mutagenesis of Human Cytomegalovirus Glycoprotein L Disproportionately Disrupts gH/gL/gO over gH/gL/pUL128-131

Cell-free and cell-to-cell spread of herpesviruses involves a core fusion apparatus comprised of the fusion protein glycoprotein B (gB) and the regulatory factor gH/gL. The human cytomegalovirus (HCMV) gH/gL/gO and gH/gL/pUL128-131 facilitate spread in different cell types. The gO and pUL128-131 components bind distinct receptors, but how the gH/gL portions of the complexes functionally compare is not understood. We previously characterized a panel of gL mutants by transient expression and showed that many were impaired for gH/gL-gB-dependent cell-cell fusion but were still able to form gH/gL/pUL128-131 and induce receptor interference. Here, the gL mutants were engineered into the HCMV BAC clones TB40/e-BAC4 (TB), TR, and Merlin (ME), which differ in their utilization of the two complexes for entry and spread. Several of the gL mutations disproportionately impacted gH/gL/gO-dependent entry and spread over gH/gL/pUL128-131 processes. The effects of some mutants could be explained by impaired gH/gL/gO assembly, but other mutants impacted gH/gL/gO function. Soluble gH/gL/gO containing the L201 mutant failed to block HCMV infection despite unimpaired binding to PDGFRα, indicating the existence of other important gH/gL/gO receptors. Another mutant (L139) enhanced the gH/gL/gO-dependent cell-free spread of TR, suggesting a "hyperactive" gH/gL/gO. Recently published crystallography and cryo-electron microscopy studies suggest structural conservation of the gH/gL underlying gH/gL/gO and gH/gL/pUL128-131. However, our data suggest important differences in the gH/gL of the two complexes and support a model in which gH/gL/gO can provide an activation signal for gB. IMPORTANCE The endemic betaherpesvirus HCMV circulates in human populations as a complex mixture of genetically distinct variants, establishes lifelong persistent infections, and causes significant disease in neonates and immunocompromised adults. This study capitalizes on our recent characterizations of three genetically distinct HCMV BAC clones to discern the functions of the envelope glycoprotein complexes gH/gL/gO and gH/gL/pUL128-13, which are promising vaccine targets that share the herpesvirus core fusion apparatus component, gH/gL. Mutations in the shared gL subunit disproportionally affected gH/gL/gO, demonstrating mechanistic differences between the two complexes, and may provide a basis for more refined evaluations of neutralizing antibodies.

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.

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.

Human Cytomegalovirus Genome Diversity in Longitudinally Collected Breast Milk Samples

Reactivation and shedding of human cytomegalovirus (HCMV) in breast milk during lactation is highly frequent in HCMV-seropositive mothers. This represents a key transmission route for postnatal HCMV infection and can lead to severe disease in preterm neonates. Little is known about HCMV strain composition or longitudinal intrahost viral population dynamics in breast milk from immunocompetent women. We performed HCMV-specific target enrichment and high-throughput sequencing of 38 breast milk samples obtained in Germany between days 10 and 60 postpartum from 15 mothers with HCMV DNA lactia, and assembled HCMV consensus sequences de novo. The genotype distribution and number of HCMV strains present in each sample were determined by quantifying genotype-specific sequence motifs in 12 hypervariable viral genes, revealing a wide range of genotypes (82/109) for these genes in the cohort and a unique, longitudinally stable strain composition in each mother. Reactivation of up to three distinct HCMV strains was detected in 8/15 of mothers, indicating that a representative subset of the woman’s HCMV reservoir might be locally reactivated early during lactation. As described previously, nucleotide diversity of samples with multiple strains was much higher than that of samples with single strains. Breast milk as a main source of postnatal mother-to-infant transmission may serve as a repository for viral diversity and thus play an essential role in the natural epidemiology of HCMV.

Genetic Variability of Human Cytomegalovirus Clinical Isolates Correlates With Altered Expression of Natural Killer Cell-Activating Ligands and IFN-γ

Human cytomegalovirus (HCMV) infection often leads to systemic disease in immunodeficient patients and congenitally infected children. Despite its clinical significance, the exact mechanisms contributing to HCMV pathogenesis and clinical outcomes have yet to be determined. One of such mechanisms involves HCMV-mediated NK cell immune response, which favors viral immune evasion by hindering NK cell-mediated cytolysis. This process appears to be dependent on the extent of HCMV genetic variation as high levels of variability in viral genes involved in immune escape have an impact on viral pathogenesis. However, the link between viral genome variations and their functional effects has so far remained elusive. Thus, here we sought to determine whether inter-host genetic variability of HCMV influences its ability to modulate NK cell responses to infection. For this purpose, five HCMV clinical isolates from a previously characterized cohort of pediatric patients with confirmed HCMV congenital infection were evaluated by next-generation sequencing (NGS) for genetic polymorphisms, phylogenetic relationships, and multiple-strain infection. We report variable levels of genetic characteristics among the selected clinical strains, with moderate variations in genome regions associated with modulation of NK cell functions. Remarkably, we show that different HCMV clinical strains differentially modulate the expression of several ligands for the NK cell-activating receptors NKG2D, DNAM-1/CD226, and NKp30. Specifically, the DNAM-1/CD226 ligand PVR/CD155 appears to be predominantly upregulated by fast-replicating (“aggressive”) HCMV isolates. On the other hand, the NGK2D ligands ULBP2/5/6 are downregulated regardless of the strain used, while other NK cell ligands (i.e., MICA, MICB, ULBP3, Nectin-2/CD112, and B7-H6) are not significantly modulated. Furthermore, we show that IFN-γ; production by NK cells co-cultured with HCMV-infected fibroblasts is directly proportional to the aggressiveness of the HCMV clinical isolates employed. Interestingly, loss of NK cell-modulating genes directed against NK cell ligands appears to be a common feature among the “aggressive” HCMV strains, which also share several gene variants across their genomes. Overall, even though further studies based on a higher number of patients would offer a more definitive scenario, our findings provide novel mechanistic insights into the impact of HCMV genetic variability on NK cell-mediated immune responses.

Whole-Genome Sequencing of Mycobacterium tuberculosis Directly from Sputum Samples

Whole-genome sequencing is a powerful, high-resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography, and mutations associated with antimicrobial resistance. The ability to sequence pathogen genomes directly from clinical specimens, without the requirement for in vitro culturing, is attractive in terms of time- and labor-saving, especially in the case of slow growing pathogens, such as Mycobacterium tuberculosis. However, clinical samples typically contain too low levels of pathogen nucleic acid, plus relatively high levels of human and natural microbiota DNA/RNA, to make this a viable option. Using a combination of whole-genome enrichment and deep sequencing, which has been proven to be a nonmutagenic approach, we can capture all known variations found within M. tuberculosis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of M. tuberculosis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

Pathogenesis of human cytomegalovirus in the immunocompromised host

Human cytomegalovirus (HCMV) is a herpesvirus that infects ~60% of adults in developed countries and more than 90% in developing countries. Usually, it is controlled by a vigorous immune response so that infections are asymptomatic or symptoms are mild. However, if the immune system is compromised, HCMV can replicate to high levels and cause serious end organ disease. Substantial progress is being made in understanding the natural history and pathogenesis of HCMV infection and disease in the immunocompromised host. Serial measures of viral load defined the dynamics of HCMV replication and are now used routinely to allow intervention with antiviral drugs in individual patients. They are also used as pharmacodynamic read-outs to evaluate prototype vaccines that may protect against HCMV replication and to define immune correlates of this protection. This novel information is informing the design of randomized controlled trials of new antiviral drugs and vaccines currently under evaluation. In this Review, we discuss immune responses to HCMV and countermeasures deployed by the virus, the establishment of latency and reactivation from it, exogenous reinfection with additional strains, pathogenesis, development of end organ disease, indirect effects of infection, immune correlates of control of replication, current treatment strategies and the evaluation of novel vaccine candidates.

Human cytomegalovirus multiple-strain infections and viral population diversity in haematopoietic stem cell transplant recipients analysed by high-throughput sequencing

Human cytomegalovirus (HCMV) is an important opportunistic pathogen in allogeneic haematopoietic stem cell transplant (HSCT) recipients. High-throughput sequencing of target-enriched libraries was performed to characterise the diversity of HCMV strains present in this high-risk group. Forty-four HCMV-DNA-positive plasma specimens (median viral input load 321 IU per library) collected at defined time points from 23 HSCT recipients within 80 days of transplantation were sequenced. The genotype distribution for 12 hypervariable HCMV genes and the number of HCMV strains present (i.e. single- vs. multiple-strain infection) were determined for 29 samples from 16 recipients. Multiple-strain infection was observed in seven of these 16 recipients, and five of these seven recipients had the donor (D)/recipient (R) HCMV-serostatus combination D + R + . A very broad range of genotypes was detected, with an intrahost composition that was generally stable over time. Multiple-strain infection was not associated with particular virological or clinical features, such as altered levels or duration of antigenaemia, development of acute graft-versus-host disease or increased mortality. In conclusion, despite relatively low viral plasma loads, a high frequency of multiple-strain HCMV infection and a high strain complexity were demonstrated in systematically collected clinical samples from this cohort early after HSCT. However, robust evaluation of the pathogenic role of intrahost viral diversity and multiple-strain infection will require studies enrolling larger numbers of recipients.

Mixed cytomegalovirus genotypes in HIV-positive mothers show compartmentalization and distinct patterns of transmission to infants

Cytomegalovirus (CMV) is the commonest cause of congenital infection and particularly so among infants born to HIV-infected women. Studies of congenital CMV infection (cCMVi) pathogenesis are complicated by the presence of multiple infecting maternal CMV strains, especially in HIV-positive women, and the large, recombinant CMV genome. Using newly developed tools to reconstruct CMV haplotypes, we demonstrate anatomic CMV compartmentalization in five HIV-infected mothers and identify the possibility of congenitally transmitted genotypes in three of their infants. A single CMV strain was transmitted in each congenitally infected case, and all were closely related to those that predominate in the cognate maternal cervix. Compared to non-transmitted strains, these congenitally transmitted CMV strains showed statistically significant similarities in 19 genes associated with tissue tropism and immunomodulation. In all infants, incident superinfections with distinct strains from breast milk were captured during follow-up. The results represent potentially important new insights into the virologic determinants of early CMV infection.

GRACy: A tool for analysing human cytomegalovirus sequence data

Modern DNA sequencing has instituted a new era in human cytomegalovirus (HCMV) genomics. A key development has been the ability to determine the genome sequences of HCMV strains directly from clinical material. This involves the application of complex and often non-standardized bioinformatics approaches to analysing data of variable quality in a process that requires substantial manual intervention. To relieve this bottleneck, we have developed GRACy (Genome Reconstruction and Annotation of Cytomegalovirus), an easy-to-use toolkit for analysing HCMV sequence data. GRACy automates and integrates modules for read filtering, genotyping, genome assembly, genome annotation, variant analysis, and data submission. These modules were tested extensively on simulated and experimental data and outperformed generic approaches. GRACy is written in Python and is embedded in a graphical user interface with all required dependencies installed by a single command. It runs on the Linux operating system and is designed to allow the future implementation of a cross-platform version. GRACy is distributed under a GPL 3.0 license and is freely available at https://bioinformatics.cvr.ac.uk/software/ with the manual and a test dataset.

[RETRACTED: High-throughput sequencing in diagnostics and prevention of herpes simplex virus infection (Herpesviridae, Alphaherpesvirinae, Simplexvirus, Human alphaherpesvirus 1)]

RETRACTEDHerpes simplex viruses types 1 (HSV-1) and 2 (HSV-2) are among the most common viruses in the human population. The clinical manifestations of HSV infection vary widely, which necessitates reliable molecular methods for the timely diagnosis of herpes virus infection, as well as for detection of mutations in the genes responsible for drug resistance. PCR is often unable to detect HSV isolates with nucleotide substitutions at the primer binding site. Sanger sequencing of the whole genome reveals mutations mainly at the consensus level, which accumulate at advanced stages of viral infection. High-throughput sequencing (HTS, next generation sequencing) offers an obvious advantage both in early diagnosis of herpes virus infection and identification of HSV variants.

Contribution of Viral Genomic Diversity to Oyster Susceptibility in the Pacific Oyster Mortality Syndrome

Juvenile Pacific oysters (Crassostrea gigas) are subjected to recurrent episodes of mass mortalities that constitute a threat for the oyster industry. This mortality syndrome named “Pacific Oyster Mortality Syndrome” (POMS) is a polymicrobial disease whose pathogenesis is initiated by a primary infection by a variant of an Ostreid herpes virus named OsHV-1 μVar. The characterization of the OsHV-1 genome during different disease outbreaks occurring in different geographic areas has revealed the existence of a genomic diversity for OsHV-1 μVar. However, the biological significance of this diversity is still poorly understood. To go further in understanding the consequences of OsHV-1 diversity on POMS, we challenged five biparental families of oysters to two different infectious environments on the French coasts (Atlantic and Mediterranean). We observed that the susceptibility to POMS can be different among families within the same environment but also for the same family between the two environments. Viral diversity analysis revealed that Atlantic and Mediterranean POMS are caused by two distinct viral populations. Moreover, we observed that different oyster families are infected by distinct viral populations within a same infectious environment. Altogether these results suggest that the co-evolutionary processes at play between OsHV-1 μVar and oyster populations have selected a viral diversity that could facilitate the infection process and the transmission in oyster populations. These new data must be taken into account in the development of novel selective breeding programs better adapted to the oyster culture environment.

Influence of Human Cytomegalovirus Glycoprotein O Polymorphism on the Inhibitory Effect of Soluble Forms of Trimer- and Pentamer-Specific Entry Receptors

Human cytomegalovirus (HCMV) envelope glycoprotein complexes, gH/gL/gO trimer and gH/gL/UL128-131 pentamer, are important for cell-free HCMV entry. While soluble NRP2-Fc (sNRP2-Fc) interferes with epithelial/endothelial cell entry through UL128, soluble platelet-derived growth factor receptor α-Fc (sPDGFRα-Fc) interacts with gO, thereby inhibiting infection of all cell types. Since gO is the most variable subunit, we investigated the influence of gO polymorphism on the inhibitory capacities of sPDGFRα-Fc and sNRP2-Fc. Accordingly, gO genotype 1c (GT1c) sequence was fully or partially replaced by gO GT2b, GT3, and GT5 sequences in the bacterial artificial chromosome (BAC) TB40-BAC4-luc background (where luc is luciferase). All mutants were tested for fibroblast and epithelial cell infectivity, for virion content of gB, gH, and gO, and for infection inhibition by sPDGFRα-Fc and sNRP2-Fc. Full-length and partial gO GT swapping may increase epithelial-to-fibroblast ratios due to subtle alterations in fibroblast and/or epithelial infectivity but without substantial changes in gB and gH levels in mutant virions. All gO GT mutants except recombinant gO GT1c/3 displayed a nearly complete inhibition at 1.25 μg/ml sPDGFRα-Fc on epithelial cells (98% versus 91%), and all experienced complete inhibition on fibroblasts (≥99%). While gO GT replacement did not influence sNRP2-Fc inhibition at 1.25 μg/ml on epithelial cells (97% to 99%), it rendered recombinant mutant GT1c/3 moderately accessible to fibroblast inhibition (40%). In contrast to the steep sPDGFRα-Fc inhibition curves (slope of >1.0), sNRP2-Fc dose-response curves on epithelial cells displayed slopes of ∼1.0, suggesting functional differences between these entry inhibitors. Our findings demonstrate that artificially generated gO recombinants rather than the major gO genotypic forms may affect the inhibitory capacities of sPDGFRα and sNRP2 in a cell type-dependent manner.IMPORTANCE Human cytomegalovirus (HCMV) is known for its broad cell tropism, as reflected by the different organs and tissues affected by HCMV infection. Hence, inhibition of HCMV entry into distinct cell types could be considered a promising therapeutic option to limit cell-free HCMV infection. Soluble forms of cellular entry receptor PDGFRα rather than those of entry receptor neuropilin-2 inhibit infection of multiple cell types. sPDGFRα specifically interacts with gO of the trimeric gH/gL/gO envelope glycoprotein complex. HCMV strains may differ with respect to the amounts of trimer in virions and the highly polymorphic gO sequence. In this study, we show that the major gO genotypes of HCMV that are also found in vivo are similarly well inhibited by sPDGFRα. Novel gO genotypic forms potentially emerging through recombination, however, may evade sPDGFRα inhibition on epithelial cells. These findings provide useful additional information for the future development of anti-HCMV therapeutic compounds based on sPDGFRα.

Specialization for Cell-Free or Cell-to-Cell Spread of BAC-Cloned Human Cytomegalovirus Strains Is Determined by Factors beyond the UL128-131 and RL13 Loci

It is widely held that clinical isolates of human cytomegalovirus (HCMV) are highly cell associated, and mutations affecting the UL128-131 and RL13 loci that arise in culture lead to the appearance of a cell-free spread phenotype. The bacterial artificial chromosome (BAC) clone Merlin (ME) expresses abundant UL128-131, is RL13 impaired, and produces low infectivity virions in fibroblasts, whereas TB40/e (TB) and TR are low in UL128-131, are RL13 intact, and produce virions of much higher infectivity. Despite these differences, quantification of spread by flow cytometry revealed remarkably similar spread efficiencies in fibroblasts. In epithelial cells, ME spread more efficiently, consistent with robust UL128-131 expression. Strikingly, ME spread far better than did TB or TR in the presence of neutralizing antibodies on both cell types, indicating that ME is not simply deficient at cell-free spread but is particularly efficient at cell-to-cell spread, whereas TB and TR cell-to-cell spread is poor. Sonically disrupted ME-infected cells contained scant infectivity, suggesting that the efficient cell-to-cell spread mechanism of ME depends on features of the intact cells such as junctions or intracellular trafficking processes. Even when UL128-131 was transcriptionally repressed, cell-to-cell spread of ME was still more efficient than that of TB or TR. Moreover, RL13 expression comparably reduced both cell-free and cell-to-cell spread of all three strains, suggesting that it acts at a stage of assembly and/or egress common to both routes of spread. Thus, HCMV strains can be highly specialized for either for cell-free or cell-to-cell spread, and these phenotypes are determined by factors beyond the UL128-131 or RL13 loci.IMPORTANCE Both cell-free and cell-to-cell spread are likely important for the natural biology of HCMV. In culture, strains clearly differ in their capacity for cell-free spread as a result of differences in the quantity and infectivity of extracellular released progeny. However, it has been unclear whether "cell-associated" phenotypes are simply the result of poor cell-free spread or are indicative of particularly efficient cell-to-cell spread mechanisms. By measuring the kinetics of spread at early time points, we were able to show that HCMV strains can be highly specialized to either cell-free or cell-to-cell mechanisms, and this was not strictly linked the efficiency of cell-free spread. Our results provide a conceptual approach to evaluating intervention strategies for their ability to limit cell-free or cell-to-cell spread as independent processes.

Within-Host Viral Diversity: A Window into Viral Evolution

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

Human Cytomegalovirus Genomes Sequenced Directly From Clinical Material: Variation, Multiple-Strain Infection, Recombination, and Gene Loss

The genomic characteristics of human cytomegalovirus (HCMV) strains sequenced directly from clinical pathology samples were investigated, focusing on variation, multiple-strain infection, recombination, and gene loss. A total of 207 datasets generated in this and previous studies using target enrichment and high-throughput sequencing were analyzed, in the process enabling the determination of genome sequences for 91 strains. Key findings were that (i) it is important to monitor the quality of sequencing libraries in investigating variation; (ii) many recombinant strains have been transmitted during HCMV evolution, and some have apparently survived for thousands of years without further recombination; (iii) mutants with nonfunctional genes (pseudogenes) have been circulating and recombining for long periods and can cause congenital infection and resulting clinical sequelae; and (iv) intrahost variation in single-strain infections is much less than that in multiple-strain infections. Future population-based studies are likely to continue illuminating the evolution, epidemiology, and pathogenesis of HCMV.

Where do we Stand after Decades of Studying Human Cytomegalovirus?

Human cytomegalovirus (HCMV), a linear double-stranded DNA betaherpesvirus belonging to the family of Herpesviridae, is characterized by widespread seroprevalence, ranging between 56% and 94%, strictly dependent on the socioeconomic background of the country being considered. Typically, HCMV causes asymptomatic infection in the immunocompetent population, while in immunocompromised individuals or when transmitted vertically from the mother to the fetus it leads to systemic disease with severe complications and high mortality rate. Following primary infection, HCMV establishes a state of latency primarily in myeloid cells, from which it can be reactivated by various inflammatory stimuli. Several studies have shown that HCMV, despite being a DNA virus, is highly prone to genetic variability that strongly influences its replication and dissemination rates as well as cellular tropism. In this scenario, the few currently available drugs for the treatment of HCMV infections are characterized by high toxicity, poor oral bioavailability, and emerging resistance. Here, we review past and current literature that has greatly advanced our understanding of the biology and genetics of HCMV, stressing the urgent need for innovative and safe anti-HCMV therapies and effective vaccines to treat and prevent HCMV infections, particularly in vulnerable populations.

Evolution and Genetic Diversity of Primate Cytomegaloviruses

Cytomegaloviruses (CMVs) infect many mammals, including humans and non-human primates (NHPs). Human cytomegalovirus (HCMV) is an important opportunistic pathogen among immunocompromised patients and represents the most common infectious cause of birth defects. HCMV possesses a large genome and very high genetic diversity. NHP-infecting CMVs share with HCMV a similar genomic organization and coding content, as well as the course of viral infection. Recent technological advances have allowed the sequencing of several HCMV strains from clinical samples and provided insight into the diversity of NHP-infecting CMVs. The emerging picture indicates that, with the exclusion of core genes (genes that have orthologs in all herpesviruses), CMV genomes are relatively plastic and diverse in terms of gene content, both at the inter- and at the intra-species level. Such variability most likely underlies the strict species-specificity of these viruses, as well as their ability to persist lifelong and with relatively little damage to their hosts. However, core genes, despite their strong conservation, also represented a target of adaptive evolution and subtle changes in their coding sequence contributed to CMV adaptation to different hosts. Indubitably, important knowledge gaps remain, the most relevant of which concerns the role of viral genetics in HCMV-associated human disease.

Polymorphisms in Human Cytomegalovirus Glycoprotein O (gO) Exert Epistatic Influences on Cell-Free and Cell-to-Cell Spread and Antibody Neutralization on gH Epitopes

Human cytomegalovirus (HCMV) glycoproteins H and L (gH/gL) can be bound by either gO or the UL128 to UL131 proteins (referred to here as UL128-131) to form complexes that facilitate entry and spread, and the complexes formed are important targets of neutralizing antibodies. Strains of HCMV vary considerably in the levels of gH/gL/gO and gH/gL/UL128-131, and this can impact infectivity and cell tropism. In this study, we investigated how natural interstrain variation in the amino acid sequence of gO influences the biology of HCMV. Heterologous gO recombinants were constructed in which 6 of the 8 alleles or genotypes (GT) of gO were analyzed in the backgrounds of strains TR and Merlin (ME). The levels of gH/gL complexes were not affected, but there were impacts on entry, spread, and neutralization by anti-gH antibodies. AD169 (AD) gO (GT1a) [referred to here as ADgO(GT1a)] drastically reduced cell-free infectivity of both strains on fibroblasts and epithelial cells. PHgO(GT2a) increased cell-free infectivity of TR in both cell types, but spread in fibroblasts was impaired. In contrast, spread of ME in both cell types was enhanced by Towne (TN) gO (GT4), despite similar cell-free infectivity. TR expressing TNgO(GT4) was resistant to neutralization by anti-gH antibodies AP86 and 14-4b, whereas ADgO(GT1a) conferred resistance to 14-4b but enhanced neutralization by AP86. Conversely, ME expressing ADgO(GT1a) was more resistant to 14-4b. These results suggest that (i) there are mechanistically distinct roles for gH/gL/gO in cell-free and cell-to-cell spread, (ii) gO isoforms can differentially shield the virus from neutralizing antibodies, and (iii) effects of gO polymorphisms are epistatically dependent on other variable loci.IMPORTANCE Advances in HCMV population genetics have greatly outpaced understanding of the links between genetic diversity and phenotypic variation. Moreover, recombination between genotypes may shuffle variable loci into various combinations with unknown outcomes. UL74(gO) is an important determinant of HCMV infectivity and one of the most diverse loci in the viral genome. By analyzing interstrain heterologous UL74(gO) recombinants, we showed that gO diversity can have dramatic impacts on cell-free and cell-to-cell spread as well as on antibody neutralization and that the manifestation of these impacts can be subject to epistatic influences of the global genetic background. These results highlight the potential limitations of laboratory studies of HCMV biology that use single, isolated genotypes or strains.

Intrahost cytomegalovirus population genetics following antibody pretreatment in a monkey model of congenital transmission

Human cytomegalovirus (HCMV) infection is the leading non-genetic cause of congenital birth defects worldwide. While several studies have addressed the genetic composition of viral populations in newborns diagnosed with HCMV, little is known regarding mother-to-child viral transmission dynamics and how therapeutic interventions may impact within-host viral populations. Here, we investigate how preexisting CMV-specific antibodies shape the maternal viral population and intrauterine virus transmission. Specifically, we characterize the genetic composition of CMV populations in a monkey model of congenital CMV infection to examine the effects of passively-infused hyperimmune globulin (HIG) on viral population genetics in both maternal and fetal compartments. In this study, 11 seronegative, pregnant monkeys were challenged with rhesus CMV (RhCMV), including a group pretreated with a standard potency HIG preparation (n = 3), a group pretreated with a high-neutralizing potency HIG preparation (n = 3), and an untreated control group (n = 5). Targeted amplicon deep sequencing of RhCMV glycoprotein B and L genes revealed that one of the three strains present in the viral inoculum (UCD52) dominated maternal and fetal viral populations. We identified minor haplotypes of this strain and characterized their dynamics. Many of the identified haplotypes were consistently detected at multiple timepoints within sampled maternal tissues, as well as across tissue compartments, indicating haplotype persistence over time and transmission between maternal compartments. However, haplotype numbers and diversity levels were not appreciably different between control, standard-potency, and high-potency pretreatment groups. We found that while the presence of maternal antibodies reduced viral load and congenital infection, it had no apparent impact on intrahost viral genetic diversity at the investigated loci. Interestingly, some minor haplotypes present in fetal and maternal-fetal interface tissues were also identified as minor haplotypes in corresponding maternal tissues, providing evidence for a loose RhCMV mother-to-fetus transmission bottleneck even in the presence of preexisting antibodies.

Whole genome sequencing of Herpes Simplex Virus 1 directly from human cerebrospinal fluid reveals selective constraints in neurotropic viruses

Herpes Simplex Virus type 1 (HSV-1) chronically infects over 70 per cent of the global population. Clinical manifestations are largely restricted to recurrent epidermal vesicles. However, HSV-1 also leads to encephalitis, the infection of the brain parenchyma, with high associated rates of mortality and morbidity. In this study, we performed target enrichment followed by direct sequencing of HSV-1 genomes, using target enrichment methods on the cerebrospinal fluid (CSF) of clinical encephalitis patients and from skin swabs of epidermal vesicles on non-encephalopathic patients. Phylogenetic analysis revealed high inter-host diversity and little population structure. In contrast, samples from different lesions in the same patient clustered with similar patterns of allelic variants. Comparison of consensus genome sequences shows HSV-1 has been freely recombining, except for distinct islands of linkage disequilibrium (LD). This suggests functional constraints prevent recombination between certain genes, notably those encoding pairs of interacting proteins. Distinct LD patterns characterised subsets of viruses recovered from CSF and skin lesions, which may reflect different evolutionary constraints in different body compartments. Functions of genes under differential constraint related to immunity or tropism and provide new hypotheses on tissue-specific mechanisms of viral infection and latency.

In vitro evolution of herpes simplex virus 1 (HSV-1) reveals selection for syncytia and other minor variants in cell culture

The large dsDNA virus herpes simplex virus 1 (HSV-1) is considered to be genetically stable, yet it can rapidly evolve in response to strong selective pressures such as antiviral treatment. Deep sequencing has revealed that clinical and laboratory isolates of this virus exist as populations that contain a mixture of minor alleles or variants, similar to many RNA viruses. The classic virology approach of plaque purifying virus creates a genetically homogenous population, but it is not clear how closely this represents the mixed virus populations found in nature. We sought to study the evolution of mixed versus highly purified HSV-1 populations in controlled cell culture conditions, to examine the impact of this genetic diversity on evolution. We found that a mixed population of HSV-1 acquired more genetic diversity and underwent a more dramatic phenotypic shift than a plaque-purified population, producing a viral population that was almost entirely syncytial after just ten passages. At the genomic level, adaptation and genetic diversification occurred at the level of minor alleles or variants in the viral population. Certain genetic variants in the mixed viral population appeared to be positively selected in cell culture, and this shift was also observed in clinical samples during their first passages in vitro. In contrast, the plaque-purified viral population did not appear to change substantially in phenotype or overall quantity of minor allele diversity. These data indicate that HSV-1 is capable of evolving rapidly in a given environment, and that this evolution is facilitated by diversity in the viral population.

Human Herpesvirus Sequencing in the Genomic Era: The Growing Ranks of the Herpetic Legion

The nine human herpesviruses are some of the most ubiquitous pathogens worldwide, causing life-long latent infection in a variety of different tissues. Human herpesviruses range from mild childhood infections to known tumour viruses and 'trolls of transplantation'. Epstein-Barr virus was the first human herpesvirus to have its whole genome sequenced; GenBank now includes thousands of herpesvirus genomes. This review will cover some of the recent advances in our understanding of herpesvirus diversity and disease that have come about as a result of new sequencing technologies, such as target enrichment and long-read sequencing. It will also look at the problem of resolving mixed-genotype infections, whether with short or long-read sequencing methods; and conclude with some thoughts on the future of the field as herpesvirus population genomics becomes a reality.

Cytomegalovirus distribution and evolution in hominines

Herpesviruses are thought to have evolved in very close association with their hosts. This is notably the case for cytomegaloviruses (CMVs; genus Cytomegalovirus) infecting primates, which exhibit a strong signal of co-divergence with their hosts. Some herpesviruses are however known to have crossed species barriers. Based on a limited sampling of CMV diversity in the hominine (African great ape and human) lineage, we hypothesized that chimpanzees and gorillas might have mutually exchanged CMVs in the past. Here, we performed a comprehensive molecular screening of all 9 African great ape species/subspecies, using 675 fecal samples collected from wild animals. We identified CMVs in eight species/subspecies, notably generating the first CMV sequences from bonobos. We used this extended dataset to test competing hypotheses with various degrees of co-divergence/number of host switches while simultaneously estimating the dates of these events in a Bayesian framework. The model best supported by the data involved the transmission of a gorilla CMV to the panine (chimpanzee and bonobo) lineage and the transmission of a panine CMV to the gorilla lineage prior to the divergence of chimpanzees and bonobos, more than 800,000 years ago. Panine CMVs then co-diverged with their hosts. These results add to a growing body of evidence suggesting that viruses with a double-stranded DNA genome (including other herpesviruses, adenoviruses, and papillomaviruses) often jumped between hominine lineages over the last few million years.

Multiple-Strain Infections of Human Cytomegalovirus With High Genomic Diversity Are Common in Breast Milk From Human Immunodeficiency Virus-Infected Women in Zambia

BACKGROUND

In developed countries, human cytomegalovirus (HCMV) is a major pathogen in congenitally infected and immunocompromised individuals, where multiple-strain infection appears linked to disease severity. The situation is less documented in developing countries. In Zambia, breast milk is a key route for transmitting HCMV and carries higher viral loads in human immunodeficiency virus (HIV)-infected women. We investigated HCMV strain diversity.

METHODS

High-throughput sequence datasets were generated from 28 HCMV-positive breast milk samples donated by 22 mothers (15 HIV-infected and 7 HIV-negative) at 4-16 weeks postpartum, then analyzed by genome assembly and novel motif-based genotyping in 12 hypervariable HCMV genes.

RESULTS

Among the 20 samples from 14 donors (13 HIV-infected and one HIV-negative) who yielded data meeting quality thresholds, 89 of the possible 109 genotypes were detected, and multiple-strain infections involving up to 5 strains per person were apparent in 9 HIV-infected women. Strain diversity was extensive among individuals but conserved compartmentally and longitudinally within them. Genotypic linkage was maintained within hypervariable UL73/UL74 and RL12/RL13/UL1 loci for virus entry and immunomodulation, but not between genes more distant from each other.

CONCLUSIONS

Breast milk from HIV-infected women contains multiple HCMV strains of high genotypic complexity and thus constitutes a major source for transmitting viral diversity.

Co-evolution of sites under immune selection shapes Epstein-Barr Virus population structure

Epstein–Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host’s innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of worldwide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main subpopulations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV.

Human cytomegalovirus genomics and transcriptomics through the lens of next-generation sequencing: revision and future challenges

The human cytomegalovirus (HCMV) genome was sequenced by hierarchical shotgun almost 30 years ago. Over these years, low and high passaged strains have been sequenced, improving, albeit still far from complete, the understanding of the coding potential, expression dynamics and diversity of wild-type HCMV strains. Next-generation sequencing (NGS) platforms have enabled a huge advancement, facilitating the comparison of differentially passaged strains, challenging diagnostics and research based on a single or reduced gene set genotyping. In addition, it allowed to link genetic features to different viral phenotypes as for example, correlating large genomic re-arrangements to viral attenuation or different mutations to antiviral resistance and cell tropism. NGS platforms provided the first high-resolution experiments to HCMV dynamics, allowing the study of intra-host viral population structures and the description of rare transcriptional events. Long-read sequencing has recently become available, helping to identify new genomic re-arrangements, partially accounting for the genetic variability displayed in clinical isolates, as well as, in changing the understanding of the HCMV transcriptome. Better knowledge of the transcriptome resulted in a vast number of new splicing events and alternative transcripts, although most of them still need additional validation. This review summarizes the sequencing efforts reached so far, discussing its approaches and providing a revision and new nuances on HCMV sequence variability in the sequencing field.

Human cytomegalovirus haplotype reconstruction reveals high diversity due to superinfection and evidence of within-host recombination

Human cytomegalovirus (HCMV) is a major global cause of congenital disability and transplant-related morbidity. Excessive levels of within-host HCMV nucleotide diversity are attributed to unexpectedly high mutation rates. Here, we show that high HCMV diversity is due to the frequent presence of mixed infections with genetically distinct strains, whereas HCMV in nonmixed infections is no more diverse than other DNA viruses. Using serial patient samples, we reconstruct viral strain haplotypes to pinpoint the timing of HCMV superinfections occurring within the study sampling time frame and uncover within-host viral recombination. From these results, we identify likely sources of infection and demonstrate probable selection for recombinant viruses. These results generate new, yet testable, insights into putative viral and host drivers of HCMV evolution and pathogenesis.

Recent sequencing efforts have led to estimates of human cytomegalovirus (HCMV) genome-wide intrahost diversity that rival those of persistent RNA viruses [Renzette N, Bhattacharjee B, Jensen JD, Gibson L, Kowalik TF (2011) PLoS Pathog 7:e1001344]. Here, we deep sequence HCMV genomes recovered from single and longitudinally collected blood samples from immunocompromised children to show that the observations of high within-host HCMV nucleotide diversity are explained by the frequent occurrence of mixed infections caused by genetically distant strains. To confirm this finding, we reconstructed within-host viral haplotypes from short-read sequence data. We verify that within-host HCMV nucleotide diversity in unmixed infections is no greater than that of other DNA viruses analyzed by the same sequencing and bioinformatic methods and considerably less than that of human immunodeficiency and hepatitis C viruses. By resolving individual viral haplotypes within patients, we reconstruct the timing, likely origins, and natural history of superinfecting strains. We uncover evidence for within-host recombination between genetically distinct HCMV strains, observing the loss of the parental virus containing the nonrecombinant fragment. The data suggest selection for strains containing the recombinant fragment, generating testable hypotheses about HCMV evolution and pathogenesis. These results highlight that high HCMV diversity present in some samples is caused by coinfection with multiple distinct strains and provide reassurance that within the host diversity for single-strain HCMV infections is no greater than for other herpesviruses.

Genotypic and Phenotypic Diversity of Herpes Simplex Virus 2 within the Infected Neonatal Population

Herpes simplex virus (HSV) causes invasive disease in half of infected neonates, resulting in significant mortality and permanent cognitive morbidity. The factors that contribute to invasive disease are not understood. This study revealed diversity among HSV isolates from infected neonates and detected the first associations between viral genetic variations and clinical disease manifestations. We found that viruses isolated from newborns with encephalitis showed enhanced spread in culture. These viruses contained protein-coding variations not found in viruses causing noninvasive disease. Many of these variations were found in proteins known to impact neurovirulence and viral spread between cells. This work advances our understanding of HSV diversity in the neonatal population and how it may impact disease outcome.

More than 14,000 neonates are infected with herpes simplex virus (HSV) annually. Approximately half display manifestations limited to the skin, eyes, or mouth (SEM disease). The rest develop invasive infections that spread to the central nervous system (CNS disease or encephalitis) or throughout the infected neonate (disseminated disease). Invasive HSV disease is associated with significant morbidity and mortality, but the viral and host factors that predispose neonates to these forms are unknown. To define viral diversity within the infected neonatal population, we evaluated 10 HSV-2 isolates from newborns with a range of clinical presentations. To assess viral fitness independently of host immune factors, we measured viral growth characteristics in cultured cells and found diverse in vitro phenotypes. Isolates from neonates with CNS disease were associated with larger plaque size and enhanced spread, with the isolates from cerebrospinal fluid (CSF) exhibiting the most robust growth. We sequenced complete viral genomes of all 10 neonatal viruses, providing new insights into HSV-2 genomic diversity in this clinical setting. We found extensive interhost and intrahost genomic diversity throughout the viral genome, including amino acid differences in more than 90% of the viral proteome. The genes encoding glycoprotein G (gG; US4), glycoprotein I (gI; US7), and glycoprotein K (gK; UL53) and viral proteins UL8, UL20, UL24, and US2 contained variants that were found in association with CNS isolates. Many of these viral proteins are known to contribute to cell spread and neurovirulence in mouse models of CNS disease. This report represents the first application of comparative pathogen genomics to neonatal HSV disease.

IMPORTANCE Herpes simplex virus (HSV) causes invasive disease in half of infected neonates, resulting in significant mortality and permanent cognitive morbidity. The factors that contribute to invasive disease are not understood. This study revealed diversity among HSV isolates from infected neonates and detected the first associations between viral genetic variations and clinical disease manifestations. We found that viruses isolated from newborns with encephalitis showed enhanced spread in culture. These viruses contained protein-coding variations not found in viruses causing noninvasive disease. Many of these variations were found in proteins known to impact neurovirulence and viral spread between cells. This work advances our understanding of HSV diversity in the neonatal population and how it may impact disease outcome.

Podcast: A podcast concerning this article is available.

Herpes Simplex Virus Disease Management and Diagnostics in the Era of High-Throughput Sequencing

Herpes simplex virus (HSV) serotypes 1 and 2 are among the most widespread human viruses. HSV disease has a complex phenotype, with symptoms that can range from mild lesions to encephalitis. In the clinical setting, this diversity of outcomes poses a major challenge, making timely disease diagnosis and treatment challenging. High-throughput sequencing (HTS) has been one of the breakthrough technologies in the modern era of molecular biology, and it is revolutionizing the study of pathogen biology and clinical diagnostics. Here, we review recent studies that have used HTS to answer questions related to the evolution of drug resistance, transmission and spread, virulence marker identification, and the design of better antiviral therapeutics for HSV. We also highlight practical considerations for handling computational analysis of HSV genomes and adoption of HTS as a routine diagnostic procedure in the clinical laboratories.

Recombination rates along the entire Epstein Barr virus genome display a highly heterogeneous landscape

Epstein Barr virus (EBV) has a large DNA genome assumed to be stable, but also subject to mutational processes such as nucleotide substitution and recombination, the latter explored to a lesser extent. Moreover, differences in the extent of recombination events across herpes sub-families were recently reported. Given the relevance of recombination in viral evolution and its possible impact in pathogenesis, we aimed to fully characterize and quantify its extension in all available EBV complete genome by assessing global and local recombination rate values (⍴/bp). Our results provide the first EBV recombination map based on recombination rates assessment, both at a global and gene by gene level, where the mean value for the entire genome was 0.035 (HPDI 0.020-0.062) ⍴/bp. We quantified how this evolutionary process changes along the EBV genome, and proved it to be non-homogeneous, since regulatory regions depicted the lowest recombination rate values while repetitive regions the highest signal. Moreover, GC content rich regions seem not to be linked to high recombination rates as previously reported. At an intragenic level, four genes (EBNA3C, EBNA3B, BRRF2 and BBLF2-BBLF3) presented a recombination rate above genome average. We specifically quantified the signal strength among different recombination-initiators previously described features and concluded that those which elicited the greatest amount of changes in ⍴/bp, TGGAG and CCCAG, were two well characterized recombination inducing motifs in eukaryotic cells. Strikingly, although TGGAG was not the most frequently detected DNA motif across the EBV genome (697 hits), it still induced a significantly greater proportion of initiation events (0.025 events/hits) than other more represented motifs, p-value = 0.04; one tailed proportion test. Present results support the idea that diversity and evolution of herpesviruses are impacted by mechanisms, such as recombination, which extends beyond the usual consideration of point mutations.

Nonprimary Maternal Cytomegalovirus Infection After Viral Shedding in Infants

BACKGROUND

Most infants with congenital Cytomegalovirus (CMV) infection are born to seropositive women as a result of maternal CMV nonprimary infection (reinfection or reactivation). Although infected children are known to transmit CMV to their seronegative mothers, the frequency and magnitude of nonprimary maternal CMV infection after exposure to viral shedding by children in their household have not been characterized.

METHODS

A cohort of Ugandan newborns and their mothers were tested weekly for CMV by quantitative polymerase chain reaction of oral swabs. Infant primary infection and maternal nonprimary infection were defined by the onset of persistent high-level oral CMV shedding. Strain-specific antibody testing was used to assess maternal reinfection. Cox regression models with time-dependent covariates were used to evaluate risk factors for nonprimary maternal infection.

RESULTS

Nonprimary CMV infection occurred in 15 of 30 mothers, all after primary infection of their infants by a median of 6 weeks (range: 1-10) in contrast to none of the mothers of uninfected infants. The median duration of maternal oral shedding lasted 18 weeks (range: 4-42) reaching a median maximum viral load of 4.69 log copies/mL (range: 3.22-5.64). Previous-week infant CMV oral quantities strongly predicted maternal nonprimary infection (hazard ratio: 2.32 per log10 DNA copies/swab increase; 95% confidence interval: 1.63-3.31). Maternal nonprimary infections were not associated with changes in strain-specific antibody responses.

CONCLUSIONS

Nonprimary CMV infection was common in mothers after primary infection in their infants, consistent with infant-to-mother transmission. Because infants frequently acquire CMV from their mothers, for example, through breast milk, this suggests the possibility of "ping-pong" infections. Additional research is needed to characterize the antigenic and genotypic strains transmitted among children and their mothers.

Host subspecific viral strains in European house mice: Murine cytomegalovirus in the Eastern (Mus musculus musculus) and Western house mouse (Mus musculus domesticus)

Murine cytomegalovirus (MCMV) has been reported from house mice (Mus musculus) worldwide, but only recently from Eastern house mice (M. m. musculus), of particular interest because they form a semi-permeable species barrier in Europe with Western house mice, M. m. domesticus. Here we report genome sequences of EastMCMV (from Eastern mice), and set these in the context of MCMV genomes from genus Mus hosts. We show EastMCMV and WestMCMV are genetically distinct. Phylogeny splitting analyses show a genome wide (94%) pattern consistent with no West-East introgression, the major exception (3.8%) being a genome-terminal region of duplicated genes involved in host immune system evasion. As expected from its function, this is a region of maintenance of ancestral polymorphism: The lack of clear splitting signal cannot be interpreted as evidence of introgression. The EastMCMV genome sequences reported here can therefore serve as a well-described resource for exploration of murid MCMV diversity.

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

Until fairly recently, genome-wide evolutionary dynamics and within-host diversity were more commonly examined in the context of small viruses than in the context of large double-stranded DNA viruses such as herpesviruses. The high mutation rates and more compact genomes of RNA viruses have inspired the investigation of population dynamics for these species, and recent data now suggest that herpesviruses might also be considered candidates for population modeling. High-throughput sequencing (HTS) and bioinformatics have expanded our understanding of herpesviruses through genome-wide comparisons of sequence diversity, recombination, allele frequency, and selective pressures. Here we discuss recent data on the mechanisms that generate herpesvirus genomic diversity and underlie the evolution of these virus families. We focus on human herpesviruses, with key insights drawn from veterinary herpesviruses and other large DNA virus families. We consider the impacts of cell culture on herpesvirus genomes and how to accurately describe the viral populations under study. The need for a strong foundation of high-quality genomes is also discussed, since it underlies all secondary genomic analyses such as RNA sequencing (RNA-Seq), chromatin immunoprecipitation, and ribosome profiling. Areas where we foresee future progress, such as the linking of viral genetic differences to phenotypic or clinical outcomes, are highlighted as well.

Clinical and biological insights from viral genome sequencing

Whole-genome sequencing (WGS) of pathogens is becoming increasingly important not only for basic research but also for clinical science and practice. In virology, WGS is important for the development of novel treatments and vaccines, and for increasing the power of molecular epidemiology and evolutionary genomics. In this Opinion article, we suggest that WGS of viruses in a clinical setting will become increasingly important for patient care. We give an overview of different WGS methods that are used in virology and summarize their advantages and disadvantages. Although there are only partially addressed technical, financial and ethical issues in regard to the clinical application of viral WGS, this technique provides important insights into virus transmission, evolution and pathogenesis.