The role of DNA repair in herpesvirus pathogenesis

Transcriptomic profiling of thymic dysregulation and viral tropism after neonatal roseolovirus infection

Introduction

Herpesviruses, including the roseoloviruses, have been linked to autoimmune disease. The ubiquitous and chronic nature of these infections have made it difficult to establish a causal relationship between acute infection and subsequent development of autoimmunity. We have shown that murine roseolovirus (MRV), which is highly related to human roseoloviruses, induces thymic atrophy and disruption of central tolerance after neonatal infection. Moreover, neonatal MRV infection results in development of autoimmunity in adult mice, long after resolution of acute infection. This suggests that MRV induces durable immune dysregulation.

Methods

In the current studies, we utilized single-cell RNA sequencing (scRNAseq) to study the tropism of MRV in the thymus and determine cellular processes in the thymus that were disrupted by neonatal MRV infection. We then utilized tropism data to establish a cell culture system.

Results

Herein, we describe how MRV alters the thymic transcriptome during acute neonatal infection. We found that MRV infection resulted in major shifts in inflammatory, differentiation and cell cycle pathways in the infected thymus. We also observed shifts in the relative number of specific cell populations. Moreover, utilizing expression of late viral transcripts as a proxy of viral replication, we identified the cellular tropism of MRV in the thymus. This approach demonstrated that double negative, double positive, and CD4 single positive thymocytes, as well as medullary thymic epithelial cells were infected by MRV in vivo. Finally, by applying pseudotime analysis to viral transcripts, which we refer to as "pseudokinetics," we identified viral gene transcription patterns associated with specific cell types and infection status. We utilized this information to establish the first cell culture systems susceptible to MRV infection in vitro.

Conclusion

Our research provides the first complete picture of roseolovirus tropism in the thymus after neonatal infection. Additionally, we identified major transcriptomic alterations in cell populations in the thymus during acute neonatal MRV infection. These studies offer important insight into the early events that occur after neonatal MRV infection that disrupt central tolerance and promote autoimmune disease.

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.

Contribution of the Epstein-Barr virus to the oncogenesis of mature T-cell lymphoproliferative neoplasms

EBV is a lymphotropic virus, member of the Herpesviridae family that asymptomatically infects more than 90% of the human population, establishing a latent infection in memory B cells. EBV exhibits complex survival and persistence dynamics, replicating its genome through the proliferation of infected B cells or production of the lytic virions. Many studies have documented the infection of T/NK cells by EBV in healthy individuals during and after primary infection. This feature has been confirmed in humanized mouse models. Together these results have challenged the hypothesis that the infection of T/NK cells per se by EBV could be a triggering event for lymphomagenesis. Extranodal NK/T-cell lymphoma (ENKTCL) and Epstein-Barr virus (EBV)-positive nodal T- and NK-cell lymphoma (NKTCL) are two EBV-associated lymphomas of T/NK cells. These two lymphomas display different clinical, histological and molecular features. However, they share two intriguing characteristics: the association with EBV and a geographical prevalence in East Asia and Latin America. In this review we will discuss the genetic characteristics of EBV in order to understand the possible role of this virus in the oncogenesis of ENKTCL and NKTCL. In addition, the main immunohistological, molecular, cytogenetic and epigenetic differences between ENKTCL and NKTCL will be discussed, as well as EBV differences in latency patterns and other viral molecular characteristics.

Comprehensive analysis of HHV-6 and HHV-7-related gene signature in prognosis and response to temozolomide of glioma

Human herpesvirus (HHV)-6 and HHV-7 have been detected in central nervous system and glioma tissue, while their exact role in glioma remains uncertain. Omics profiles and clinical information were downloaded from public databases, including The Cancer Genome Atlas cohort for training set and the Chinese Glioma Genome Atlas cohorts for validation sets. Differentially expressed genes between HHV-6 and HHV-7 infected or noninfected glioma patients were screened for establishing the HHV-6 and HHV-7 infection (HI) model through Lasso regression analysis. Bioinformatics methods were used to analyze the correlation between HI scores and prognosis, metastasis in glioma patients. Predictable efficacy of HI in temozolomide-resistance and HI-related genetic signatures were also explored. The HI model was constructed as: Risk score = (0.014709*DIRAS3) + (0.029787*TEX26) + (0.223492*FBXO39) + (0.074951*MYBL1) + (0.060202*HILS1). The five gene signature showed good performance in predicting survival time for glioma patients, while higher HI score is correlated with malignant features. Moreover, DNA mismatch repair genes were augmented in glioma patients with higher HI score as well as nonresponse to temozolomide treatment, which was in parallel with the transcriptomic result of temozolomide-resistant glioma cell. Targeting the five gene signature is beneficial for prognosis of glioma patients, especially in glioma patients underwent temozolomide treatment.

Genomic changes in Kaposi Sarcoma-associated Herpesvirus and their clinical correlates

Kaposi sarcoma (KS), a common HIV-associated malignancy, presents a range of clinicopathological features. Kaposi sarcoma-associated herpesvirus (KSHV) is its etiologic agent, but the contribution of viral genomic variation to KS development is poorly understood. To identify potentially influential viral polymorphisms, we characterized KSHV genetic variation in 67 tumors from 1-4 distinct sites from 29 adults with advanced KS in Kampala, Uganda. Whole KSHV genomes were sequenced from 20 tumors with the highest viral load, whereas only polymorphic genes were screened by PCR and sequenced from 47 other tumors. Nine individuals harbored ≥1 tumors with a median 6-fold over-coverage of a region centering on K5 and K6 genes. K8.1 gene was inactivated in 8 individuals, while 5 had mutations in the miR-K10 microRNA coding sequence. Recurring inter-host polymorphisms were detected in K4.2 and K11.2. The K5-K6 region rearrangement breakpoints and K8.1 mutations were all unique, indicating that they arise frequently de novo. Rearrangement breakpoints were associated with potential G-quadruplex and Z-DNA forming sequences. Exploratory evaluations of viral mutations with clinical and tumor traits were conducted by logistic regression without multiple test corrections. K5-K6 over-coverage and K8.1 inactivation were tentatively correlated (p<0.001 and p = 0.005, respectively) with nodular rather than macular tumors, and with individuals that had lesions in ≤4 anatomic areas (both p≤0.01). Additionally, a trend was noted for miR-K10 point mutations and lower survival rates (HR = 4.11, p = 0.053). Two instances were found of distinct tumors within an individual sharing the same viral mutation, suggesting metastases or transmission of the aberrant viruses within the host. To summarize, KSHV genomes in tumors frequently have over-representation of the K5-K6 region, as well as K8.1 and miR-K10 mutations, and each might be associated with clinical phenotypes. Studying their possible effects may be useful for understanding KS tumorigenesis and disease progression.

Pathways involved in viral oncogenesis: New perspectives from virus-host protein interactomics

Oncogenic viruses are among the apparent causes of cancer-associated mortality. It was estimated that 12% to 15% of human malignancies are linked to oncoviruses. Although modernist strategies and traditional genetic studies have defined host-pathogen interactions of the oncoviruses, their host functions which are critical for the establishment of infection still remain mysterious. However, over the last few years, it has become clear that infections hijack and modify cellular pathways for their benefit. In this context, we constructed the virus-host protein interaction networks of seven oncoviruses (EBV, HBV, HCV, HTLV-1, HHV8, HPV16, and HPV18), and revealed cellular pathways hijacking as a result of oncogenic virus infection. Several signaling pathways/processes such as TGF-β signaling, cell cycle, retinoblastoma tumor suppressor protein, and androgen receptor signaling were mutually targeted by viruses to induce oncogenesis. Besides, cellular pathways specific to a certain virus were detected. By this study, we believe that we improve the understanding of the molecular pathogenesis of viral oncogenesis and provide information in setting new targets for treatment, prognosis, and diagnosis.

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.

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.

Genetic diversity of equine herpesvirus 1 isolated from neurological, abortigenic and respiratory disease outbreaks

Equine herpesvirus 1 (EHV‐1) causes respiratory disease, abortion, neonatal death and neurological disease in equines and is endemic in most countries. The viral factors that influence EHV‐1 disease severity are poorly understood, and this has hampered vaccine development. However, the N752D substitution in the viral DNA polymerase catalytic subunit has been shown statistically to be associated with neurological disease. This has given rise to the term “neuropathic strain,” even though strains lacking the polymorphism have been recovered from cases of neurological disease. To broaden understanding of EHV‐1 diversity in the field, 78 EHV‐1 strains isolated over a period of 35 years were sequenced. The great majority of isolates originated from the United Kingdom and included in the collection were low passage isolates from respiratory, abortigenic and neurological outbreaks. Phylogenetic analysis of regions spanning 80% of the genome showed that up to 13 viral clades have been circulating in the United Kingdom and that most of these are continuing to circulate. Abortion isolates grouped into nine clades, and neurological isolates grouped into five. Most neurological isolates had the N752D substitution, whereas most abortion isolates did not, although three of the neurological isolates from linked outbreaks had a different polymorphism. Finally, bioinformatic analysis suggested that recombination has occurred between EHV‐1 clades, between EHV‐1 and equine herpesvirus 4, and between EHV‐1 and equine herpesvirus 8.

Heterogeneity of the Epstein-Barr Virus (EBV) Major Internal Repeat Reveals Evolutionary Mechanisms of EBV and a Functional Defect in the Prototype EBV Strain B95-8

Epstein-Barr virus (EBV) is a ubiquitous pathogen of humans that can cause several types of lymphoma and carcinoma. Like other herpesviruses, EBV has diversified through both coevolution with its host and genetic exchange between virus strains. Sequence analysis of the EBV genome is unusually challenging because of the large number and lengths of repeat regions within the virus. Here we describe the sequence assembly and analysis of the large internal repeat 1 of EBV (IR1; also known as the BamW repeats) for more than 70 strains. The diversity of the latency protein EBV nuclear antigen leader protein (EBNA-LP) resides predominantly within the exons downstream of IR1. The integrity of the putative BWRF1 open reading frame (ORF) is retained in over 80% of strains, and deletions truncating IR1 always spare BWRF1. Conserved regions include the IR1 latency promoter (Wp) and one zone upstream of and two within BWRF1. IR1 is heterogeneous in 70% of strains, and this heterogeneity arises from sequence exchange between strains as well as from spontaneous mutation, with interstrain recombination being more common in tumor-derived viruses. This genetic exchange often incorporates regions of <1 kb, and allelic gene conversion changes the frequency of small regions within the repeat but not close to the flanks. These observations suggest that IR1-and, by extension, EBV-diversifies through both recombination and breakpoint repair, while concerted evolution of IR1 is driven by gene conversion of small regions. Finally, the prototype EBV strain B95-8 contains four nonconsensus variants within a single IR1 repeat unit, including a stop codon in the EBNA-LP gene. Repairing IR1 improves EBNA-LP levels and the quality of transformation by the B95-8 bacterial artificial chromosome (BAC).IMPORTANCE Epstein-Barr virus (EBV) infects the majority of the world population but causes illness in only a small minority of people. Nevertheless, over 1% of cancers worldwide are attributable to EBV. Recent sequencing projects investigating virus diversity to see if different strains have different disease impacts have excluded regions of repeating sequence, as they are more technically challenging. Here we analyze the sequence of the largest repeat in EBV (IR1). We first characterized the variations in protein sequences encoded across IR1. In studying variations within the repeat of each strain, we identified a mutation in the main laboratory strain of EBV that impairs virus function, and we suggest that tumor-associated viruses may be more likely to contain DNA mixed from two strains. The patterns of this mixing suggest that sequences can spread between strains (and also within the repeat) by copying sequence from another strain (or repeat unit) to repair DNA damage.

Oyster RNA-seq Data Support the Development of Malacoherpesviridae Genomics

The family of double-stranded DNA (dsDNA) Malacoherpesviridae includes viruses able to infect marine mollusks and detrimental for worldwide aquaculture production. Due to fast-occurring mortality and a lack of permissive cell lines, the available data on the few known Malacoherpesviridae provide only partial support for the study of molecular virus features, life cycle, and evolutionary history. Following thorough data mining of bivalve and gastropod RNA-seq experiments, we used more than five million Malacoherpesviridae reads to improve the annotation of viral genomes and to characterize viral InDels, nucleotide stretches, and SNPs. Both genome and protein domain analyses confirmed the evolutionary diversification and gene uniqueness of known Malacoherpesviridae. However, the presence of Malacoherpesviridae-like sequences integrated within genomes of phylogenetically distant invertebrates indicates broad diffusion of these viruses and indicates the need for confirmatory investigations. The manifest co-occurrence of OsHV-1 genotype variants in single RNA-seq samples of Crassostrea gigas provide further support for the Malacoherpesviridae diversification. In addition to simple sequence motifs inter-punctuating viral ORFs, recombination-inducing sequences were found to be enriched in the OsHV-1 and AbHV1-AUS genomes. Finally, the highly correlated expression of most viral ORFs in multiple oyster samples is consistent with the burst of viral proteins during the lytic phase.

Hytrosaviruses: current status and perspective

Salivary gland hytrosaviruses (SGHVs) are entomopathogenic dsDNA, enveloped viruses that replicate in the salivary glands (SGs) of the adult dipterans, Glossina spp (GpSGHV) and Musca domestica (MdSGHV). Although belonging to the same virus family (Hytrosaviridae), SGHVs have distinct morphologies and pathobiologies. Two GpSGHV strains potentially account for the differential pathologies in lab-bred tsetse. New data suggest incorporation of host-derived cellular proteins and lipids into mature SGHVs. In addition to within the SGs, MdSGHV undergoes limited replication in the corpora allata, potentially disrupting hormone biosynthesis, and GpSGHV replicates in the milk glands providing a transmission conduit to progeny tsetse. Whereas MdSGHV is a potential biocontrol agent, the vertically transmitted GpSGHV is unsuitable for tsetse vector control but does jeopardize tsetse mass rearing.

Herpes Simplex Virus Latency: The DNA Repair-Centered Pathway

Like all herpesviruses, herpes simplex virus 1 (HSV1) is able to produce lytic or latent infections depending on the host cell type. Lytic infections occur in a broad range of cells while latency is highly specific for neurons. Although latency suggests itself as an attractive target for novel anti-HSV1 therapies, progress in their development has been slowed due in part to a lack of agreement about the basic biochemical mechanisms involved. Among the possibilities being considered is a pathway in which DNA repair mechanisms play a central role. Repair is suggested to be involved in both HSV1 entry into latency and reactivation from it. Here I describe the basic features of the DNA repair-centered pathway and discuss some of the experimental evidence supporting it. The pathway is particularly attractive because it is able to account for important features of the latent response, including the specificity for neurons, the specificity for neurons of the peripheral compared to the central nervous system, the high rate of genetic recombination in HSV1-infected cells, and the genetic identity of infecting and reactivated virus.

Possible Mechanisms of Acquisition of Herpesvirus Virokines

The genomes of certain types of human and primate herpesviruses contain functional homologs of important host cytokines (IL-6, IL-17, and IL-10), or so-called virokines. Virokines can interact with immune cell receptors, transmit a signal to them, and thus switch the type of immune response that facilitates viral infection development. In this work, we have summarized possible ways of virokine origin and proposed an evolutionary scenario of virokine acquisition with involvement of retroviral coinfection of the host. This scenario is probably valid for vIL-6 of HHV-8 and MRV-5 viruses, vIL-17 of HVS virus, and vIL-10 of HHV-4, Bonobo-HV, RhLCV, and BaLCV viruses. The ability to acquire cytokine genes allows herpesviruses to implement unique strategies of avoiding the immune response and provides them an evolutionary advantage: more than 90% of the host population can be chronically infected with different herpesviruses. It is possible that the biological success of herpesviruses can be partially due to their cooperation with another group of viruses. This hypothesis emphasizes the importance of studies on the reciprocal influence of pathogens on their coinfection, as well as their impact on the host organism.

Recombination in viruses: mechanisms, methods of study, and evolutionary consequences

Recombination is a pervasive process generating diversity in most viruses. It joins variants that arise independently within the same molecule, creating new opportunities for viruses to overcome selective pressures and to adapt to new environments and hosts. Consequently, the analysis of viral recombination attracts the interest of clinicians, epidemiologists, molecular biologists and evolutionary biologists. In this review we present an overview of three major areas related to viral recombination: (i) the molecular mechanisms that underlie recombination in model viruses, including DNA-viruses (Herpesvirus) and RNA-viruses (Human Influenza Virus and Human Immunodeficiency Virus), (ii) the analytical procedures to detect recombination in viral sequences and to determine the recombination breakpoints, along with the conceptual and methodological tools currently used and a brief overview of the impact of new sequencing technologies on the detection of recombination, and (iii) the major areas in the evolutionary analysis of viral populations on which recombination has an impact. These include the evaluation of selective pressures acting on viral populations, the application of evolutionary reconstructions in the characterization of centralized genes for vaccine design, and the evaluation of linkage disequilibrium and population structure.