Presence and role of cytosine methylation in DNA viruses of animals

Host-encoded DNA methyltransferases modify the epigenome and host tropism of invading phages

Restriction modification (RM) systems are ubiquitous bacterial defense systems; however, some phages evade RM system and adapt to their bacterial hosts. In such cases, phages are thought to stochastically acquire DNA methylation from host-encoded DNA methyltransferases (MTases), facilitating host adaptation. However, no studies have directly compared the methylomes of host bacteria and their infecting phages. Here, we demonstrate the epigenetic landscape of adapted phages with diverse infection histories, focusing on the broad host-range phage KHP30T as its adapts to three Helicobacter pylori strains. Using a multistage infection system, we observed that the adapted phages displayed significantly high titers against the last infected H. pylori strain, suggesting an attendant change in host tropism. Single-molecule real-time sequencing revealed that methylated motifs were predominantly shared between the adapted phages and their most recent host. Our findings enhance our understanding of epigenetic phage-host interactions, which have significant implications for microbial ecology.

The ageing virus hypothesis: Epigenetic ageing beyond the Tree of Life

A recent thought-provoking theory argues that complex organisms using epigenetic information for their normal development and functioning must irreversibly age as a result of epigenetic signal loss. Importantly, the scope of this theory could be considerably expanded, with scientific benefits, by analyzing epigenetic ageing beyond the borders of the Tree of Life. Viruses that use epigenetic signals for their normal functioning may also age, that is, present an increasing risk of failing to complete their individual life cycle and to disappear with time. As viruses are ancient, abundant, and infect a considerable diversity of hosts, the ageing virus hypothesis, if verified, would have important consequences for many fields of the Life sciences. Uncovering ageing viruses would integrate the most abundant and biologically central entities on Earth into theories of ageing, enhance virology, gerontology, evolutionary biology, molecular ecology, genomics, and possibly medicine through the development of new therapies manipulating viral ageing.

Viral genomic methylation and the interspecies evolutionary relationships of ranavirus

Ranaviruses are capable of infecting both wild and farmed fish, amphibians, and reptiles, leading to significant economic losses and ecological risks. Currently, ranaviruses have been found in at least 175 species spanning six continents. Except for Singapore grouper iridovirus (SGIV), ranavirus genomes are generally regarded as highly methylated. Nevertheless, our comprehension of the methylation characteristics within ranaviruses remains limited. Despite the numerous genomes currently included in the GenBank database, a complete phylogenetic tree for ranaviruses has not yet been determined, and interspecific evolutionary relationships among ranaviruses have not been thoroughly investigated. In this study, the whole-genome methylation profile of mandarin fish ranavirus (MRV; a ranavirus) was investigated, revealing a methylation level of 16.04%, and hypomethylation of the MRV genome was detrimental to viral replication, speculating the genome methylation may play an important role in MRV replication. Furthermore, by combining with whole-genome DNA sequence phylogenetic analyses, we propose the possibility of an interspecies evolutionary relationship among ranaviruses, with the presence of four distinct evolutionary lineages within ranavirus evolution: "SGIV, SCRAV(MRV/LMBV), EHNV/ENARV/ATV, and CMTV/FV3", which might be also supported by the genomic collinearity, natural host range and host habitats. Furthermore, ranavirus genomic methylation levels may provide additional evidence for this hypothesis, but further proof is needed. Our work enhances the understanding of the role of genome methylation in ranaviruses and is beneficial for the prevention and control of ranavirus diseases; simultaneously, the proposed evolutionary hypothesis of ranavirus provides novel insights and ideas for exploring the evolutionary trajectory of viruses.

Are Viruses Taxonomic Units? A Protein Domain and Loop-Centric Phylogenomic Assessment

Virus taxonomy uses a Linnaean-like subsumption hierarchy to classify viruses into taxonomic units at species and higher rank levels. Virus species are considered monophyletic groups of mobile genetic elements (MGEs) often delimited by the phylogenetic analysis of aligned genomic or metagenomic sequences. Taxonomic units are assumed to be independent organizational, functional and evolutionary units that follow a 'natural history' rationale. Here, I use phylogenomic and other arguments to show that viruses are not self-standing genetically-driven systems acting as evolutionary units. Instead, they are crucial components of holobionts, which are units of biological organization that dynamically integrate the genetics, epigenetic, physiological and functional properties of their co-evolving members. Remarkably, phylogenomic analyses show that viruses share protein domains and loops with cells throughout history via massive processes of reticulate evolution, helping spread evolutionary innovations across a wider taxonomic spectrum. Thus, viruses are not merely MGEs or microbes. Instead, their genomes and proteomes conduct cellularly integrated processes akin to those cataloged by the GO Consortium. This prompts the generation of compositional hierarchies that replace the 'is-a-kind-of' by a 'is-a-part-of' logic to better describe the mereology of integrated cellular and viral makeup. My analysis demands a new paradigm that integrates virus taxonomy into a modern evolutionarily centered taxonomy of organisms.

Borg extrachromosomal elements of methane-oxidizing archaea have conserved and expressed genetic repertoires

Borgs are huge extrachromosomal elements (ECE) of anaerobic methane-consuming "Candidatus Methanoperedens" archaea. Here, we used nanopore sequencing to validate published complete genomes curated from short reads and to reconstruct new genomes. 13 complete and four near-complete linear genomes share 40 genes that define a largely syntenous genome backbone. We use these conserved genes to identify new Borgs from peatland soil and to delineate Borg phylogeny, revealing two major clades. Remarkably, Borg genes encoding nanowire-like electron-transferring cytochromes and cell surface proteins are more highly expressed than those of host Methanoperedens, indicating that Borgs augment the Methanoperedens activity in situ. We reconstructed the first complete 4.00 Mbp genome for a Methanoperedens that is inferred to be a Borg host and predicted its methylation motifs, which differ from pervasive TC and CC methylation motifs of the Borgs. Thus, methylation may enable Methanoperedens to distinguish their genomes from those of Borgs. Very high Borg to Methanoperedens ratios and structural predictions suggest that Borgs may be capable of encapsulation. The findings clearly define Borgs as a distinct class of ECE with shared genomic signatures, establish their diversification from a common ancestor with genetic inheritance, and raise the possibility of periodic existence outside of host cells.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Hypermethylated genome of a fish vertebrate iridovirus ISKNV plays important roles in viral infection

Iridoviruses are nucleocytoplasmic large dsDNA viruses that infect invertebrates and ectothermic vertebrates. The hypermethylated genome of vertebrate iridoviruses is unique among animal viruses. However, the map and function of iridovirus genomic methylation remain unknown. Herein, the methylated genome of Infectious spleen and kidney necrosis virus (ISKNV, a fish iridovirus), and its role in viral infection, are investigated. The methylation level of ISKNV is 23.44%. The hypermethylated genome is essential for ISKNV amplification, but there is no correlation between hypermethylation and viral gene expression. The hypomethylated ISKNV (obtained via 5-Azacytidine) activates a strong immunoreaction in vitro and reduces its pathogenicity in vivo. The unmethylated viral DNA can induce a stronger immunoreaction in vitro, whereas inactivated hypomethylated ISKNV can induce a stronger immunoreaction in vivo, suggesting ISKNV may evade from immune system by increasing its genome methylation level. Our work provides new insights into the role of genome methylation in viral infection.

Roving methyltransferases generate a mosaic epigenetic landscape and influence evolution in Bacteroides fragilis group

Three types of DNA methyl modifications have been detected in bacterial genomes, and mechanistic studies have demonstrated roles for DNA methylation in physiological functions ranging from phage defense to transcriptional control of virulence and host-pathogen interactions. Despite the ubiquity of methyltransferases and the immense variety of possible methylation patterns, epigenomic diversity remains unexplored for most bacterial species. Members of the Bacteroides fragilis group (BFG) reside in the human gastrointestinal tract as key players in symbiotic communities but also can establish anaerobic infections that are increasingly multi-drug resistant. In this work, we utilize long-read sequencing technologies to perform pangenomic (n = 383) and panepigenomic (n = 268) analysis of clinical BFG isolates cultured from infections seen at the NIH Clinical Center over four decades. Our analysis reveals that single BFG species harbor hundreds of DNA methylation motifs, with most individual motif combinations occurring uniquely in single isolates, implying immense unsampled methylation diversity within BFG epigenomes. Mining of BFG genomes identified more than 6000 methyltransferase genes, approximately 1000 of which were associated with intact prophages. Network analysis revealed substantial gene flow among disparate phage genomes, implying a role for genetic exchange between BFG phages as one of the ultimate sources driving BFG epigenome diversity.

Selective Depletion of ZAP-Binding CpG Motifs in HCV Evolution

Hepatitis C virus (HCV) is a bloodborne pathogen that can cause chronic liver disease and hepatocellular carcinoma. The loss of CpGs from virus genomes allows escape from restriction by the host zinc-finger antiviral protein (ZAP). The evolution of HCV in the human host has not been explored in the context of CpG depletion. We analysed 2616 full-length HCV genomes from 1977 to 2021. During the four decades of evolution in humans, we found that HCV genomes have become significantly depleted in (a) CpG numbers, (b) CpG O/E ratios (i.e., relative abundance of CpGs), and (c) the number of ZAP-binding motifs. Interestingly, our data suggests that the loss of CpGs in HCV genomes over time is primarily driven by the loss of ZAP-binding motifs; thus suggesting a yet unknown role for ZAP-mediated selection pressures in HCV evolution. The HCV core gene is significantly enriched for the number of CpGs and ZAP-binding motifs. In contrast to the rest of the HCV genome, the loss of CpGs from the core gene does not appear to be driven by ZAP-mediated selection. This work highlights CpG depletion in HCV genomes during their evolution in humans and the role of ZAP-mediated selection in HCV evolution.

The Mutagenic Consequences of DNA Methylation within and across Generations

DNA methylation is an epigenetic modification with wide-ranging consequences across the life of an organism. This modification can be stable, persisting through development despite changing environmental conditions. However, in other contexts, DNA methylation can also be flexible, underlying organismal phenotypic plasticity. One underappreciated aspect of DNA methylation is that it is a potent mutagen; methylated cytosines mutate at a much faster rate than other genetic motifs. This mutagenic property of DNA methylation has been largely ignored in eco-evolutionary literature, despite its prevalence. Here, we explore how DNA methylation induced by environmental and other factors could promote mutation and lead to evolutionary change at a more rapid rate and in a more directed manner than through stochastic genetic mutations alone. We argue for future research on the evolutionary implications of DNA methylation driven mutations both within the lifetime of organisms, as well as across timescales.

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics

BACKGROUND

DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life.

RESULTS

We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark.

CONCLUSIONS

Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome. Video Abstract.

When does hepatitis B virus meet long-stranded noncoding RNAs?

Hepatitis B virus (HBV) infection in humans and its associated diseases are long-standing problems. HBV can produce a large number of non-self-molecules during its life cycle, which acts as targets for innate immune recognition and initiation. Among these, interferon and its large number of downstream interferon-stimulated gene molecules are important early antiviral factors. However, the development of an effective antiviral immune response is not simple and depends not only on the delicate regulation of the immune response but also on the various mechanisms of virus-related immune escape and immune tolerance. Therefore, despite there being a relatively well-established consensus on the major pathways of the antiviral response and their component molecules, the complete clearance of HBV remains a challenge in both basic and clinical research. Long-noncoding RNAs (lncRNAs) are generally >200 bp in length and perform different functions in the RNA strand encoding the protein. As an important part of the IFN-inducible genes, interferon-stimulated lncRNAs are involved in the regulation of several HBV infection-related pathways. This review traces the basic elements of such pathways and characterizes the various recent targets of lncRNAs, which not only complement the regulatory mechanisms of pathways related to chronic HBV infection, fibrosis, and cancer promotion but also present with new potential therapeutic targets for controlling HBV infection and the malignant transformation of hepatocytes.

Clinical sequencing uncovers the genomic characteristics and mutation spectrum of the 2018 African swine fever virus in Guangdong, China

African swine fever (ASF) outbreak have caused tremendous economic loss to the pig industry in China since its emergence in August 2018. Previous studies revealed that many published sequences are not suitable for detailed analyses due to the lack of data regarding quality parameters and methodology, and outdated annotations. Thus, high-quality genomes of highly pathogenic strains that can be used as references for early Chinese ASF outbreaks are still lacking, and little is known about the features of intra-host variants of ASF virus (ASFV). In this study, a full genome sequencing of clinical samples from the first ASF outbreak in Guangdong in 2018 was performed using MGI (MGI Tech Co., Ltd., Shenzhen, China) and Nanopore sequencing platforms, followed by Sanger sequencing to verify the variations. With 22 sequencing corrections, we obtained a high-quality genome of one of the earliest virulent isolates, GZ201801_2. After proofreading, we improved (add or modify) the annotations of this isolate using the whole genome alignment with Georgia 2007/1. Based on the complete genome sequence, we constructed the methylation profiles of early ASFV strains in China and predicted the potential 5mC and 6mA methylation sites, which are likely involved in metabolism, transcription, and replication. Additionally, the intra-host single nucleotide variant distribution and mutant allele frequency in the clinical samples of early strain were determined for the first time and found a strong preference for A and T substitution mutation, non-synonymous mutations, and mutations that resulted in amino acid substitutions into Lysine. In conclusion, this study provides a high-quality genome sequence, updated genome annotation, methylation profile, and mutation spectrum of early ASFV strains in China, thereby providing a reference basis for further studies on the evolution, transmission, and virulence of ASFV.

Integrative metabolomic and proteomic signatures define clinical outcomes in severe COVID-19

The coronavirus disease-19 (COVID-19) pandemic has ravaged global healthcare with previously unseen levels of morbidity and mortality. In this study, we performed large-scale integrative multi-omics analyses of serum obtained from COVID-19 patients with the goal of uncovering novel pathogenic complexities of this disease and identifying molecular signatures that predict clinical outcomes. We assembled a network of protein-metabolite interactions through targeted metabolomic and proteomic profiling in 330 COVID-19 patients compared to 97 non-COVID, hospitalized controls. Our network identified distinct protein-metabolite cross talk related to immune modulation, energy and nucleotide metabolism, vascular homeostasis, and collagen catabolism. Additionally, our data linked multiple proteins and metabolites to clinical indices associated with long-term mortality and morbidity. Finally, we developed a novel composite outcome measure for COVID-19 disease severity based on metabolomics data. The model predicts severe disease with a concordance index of around 0.69, and shows high predictive power of 0.83-0.93 in two independent datasets.

Metagenomic Sequencing for Microbial DNA in Human Samples: Emerging Technological Advances

Whole genome metagenomic sequencing is a powerful platform enabling the simultaneous identification of all genes from entirely different kingdoms of organisms in a complex sample. This technology has revolutionised multiple areas from microbiome research to clinical diagnoses. However, one of the major challenges of a metagenomic study is the overwhelming non-microbial DNA present in most of the host-derived specimens, which can inundate the microbial signals and reduce the sensitivity of microorganism detection. Various host DNA depletion methods to facilitate metagenomic sequencing have been developed and have received considerable attention in this context. In this review, we present an overview of current host DNA depletion approaches along with explanations of their underlying principles, advantages and disadvantages. We also discuss their applications in laboratory microbiome research and clinical diagnoses and, finally, we envisage the direction of the further perfection of metagenomic sequencing in samples with overabundant host DNA.

Microbial DNA Enrichment Promotes Adrenomedullary Inflammation, Catecholamine Secretion, and Hypertension in Obese Mice

Background

Obesity is an established risk factor for hypertension. Although obesity‐induced gut barrier breach leads to the leakage of various microbiota‐derived products into host circulation and distal organs, the roles of microbiota in mediating the development of obesity‐associated adrenomedullary disorders and hypertension have not been elucidated. We seek to explore the impacts of microbial DNA enrichment on inducing obesity‐related adrenomedullary abnormalities and hypertension.

Methods and Results

Obesity was accompanied by remarkable bacterial DNA accumulation and elevated inflammation in the adrenal glands. Gut microbial DNA containing extracellular vesicles (mEVs) were readily leaked into the bloodstream and infiltrated into the adrenal glands in obese mice, causing microbial DNA enrichment. In lean wild‐type mice, adrenal macrophages expressed CRIg (complement receptor of the immunoglobulin superfamily) that efficiently blocks the infiltration of gut mEVs. In contrast, the adrenal CRIg+ cell population was greatly decreased in obese mice. In lean CRIg^−/− or C3^−/− (complement component 3) mice intravenously injected with gut mEVs, adrenal microbial DNA accumulation elevated adrenal inflammation and norepinephrine secretion, concomitant with hypertension. In addition, microbial DNA promoted inflammatory responses and norepinephrine production in rat pheochromocytoma PC12 cells treated with gut mEVs. Depletion of microbial DNA cargo markedly blunted the effects of gut mEVs. We also validated that activation of cGAS (cyclic GMP‐AMP synthase)/STING (cyclic GMP–AMP receptor stimulator of interferon genes) signaling is required for the ability of microbial DNA to trigger adrenomedullary dysfunctions in both in vivo and in vitro experiments. Restoring CRIg+ cells in obese mice decreased microbial DNA abundance, inflammation, and hypertension.

Conclusions

The leakage of gut mEVs leads to adrenal enrichment of microbial DNA that are pathogenic to induce obesity‐associated adrenomedullary abnormalities and hypertension. Recovering the CRIg+ macrophage population attenuates obesity‐induced adrenomedullary disorders.

Innate Immune Modulation Induced by EBV Lytic Infection Promotes Endothelial Cell Inflammation and Vascular Injury in Scleroderma

Microvascular injury is considered an initial event in the pathogenesis of scleroderma and endothelial cells are suspected of being the target of the autoimmune process seen in the disease. EBV has long been proposed as a trigger for autoimmune diseases, including scleroderma. Nevertheless, its contribution to the pathogenic process remains poorly understood. In this study, we report that EBV lytic antigens are detected in scleroderma dermal vessels, suggesting that endothelial cells might represent a target for EBV infection in scleroderma skin. We show that EBV DNA load is remarkably increased in peripheral blood, plasma and circulating monocytes from scleroderma patients compared to healthy EBV carriers, and that monocytes represent the prominent subsets of EBV-infected cells in scleroderma. Given that monocytes have the capacity to adhere to the endothelium, we then investigated whether monocyte-associated EBV could infect primary human endothelial cells. We demonstrated that endothelial cells are infectable by EBV, using human monocytes bound to recombinant EBV as a shuttle, even though cell-free virus failed to infect them. We show that EBV induces activation of TLR9 innate immune response and markers of vascular injury in infected endothelial cells and that up-regulation is associated with the expression of EBV lytic genes in infected cells. EBV innate immune modulation suggests a novel mechanism mediating inflammation, by which EBV triggers endothelial cell and vascular injury in scleroderma. In addition, our data point to up-regulation of EBV DNA loads as potential biomarker in developing vasculopathy in scleroderma. These findings provide the framework for the development of novel therapeutic interventions to shift the scleroderma treatment paradigm towards antiviral therapies.

Comparative Genomics and Integrated Network Approach Unveiled Undirected Phylogeny Patterns, Co-mutational Hot Spots, Functional Cross Talk, and Regulatory Interactions in SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has resulted in 92 million cases in a span of 1 year. The study focuses on understanding population-specific variations attributing its high rate of infections in specific geographical regions particularly in the United States. Rigorous phylogenomic network analysis of complete SARS-CoV-2 genomes (245) inferred five central clades named a (ancestral), b, c, d, and e (subtypes e1 and e2). Clade d and subclade e2 were found exclusively comprised of U.S. strains. Clades were distinguished by 10 co-mutational combinations in Nsp3, ORF8, Nsp13, S, Nsp12, Nsp2, and Nsp6. Our analysis revealed that only 67.46% of single nucleotide polymorphism (SNP) mutations were at the amino acid level. T1103P mutation in Nsp3 was predicted to increase protein stability in 238 strains except for 6 strains which were marked as ancestral type, whereas co-mutation (P409L and Y446C) in Nsp13 were found in 64 genomes from the United States highlighting its 100% co-occurrence. Docking highlighted mutation (D614G) caused reduction in binding of spike proteins with angiotensin-converting enzyme 2 (ACE2), but it also showed better interaction with the TMPRSS2 receptor contributing to high transmissibility among U.S. strains. We also found host proteins, MYO5A, MYO5B, and MYO5C, that had maximum interaction with viral proteins (nucleocapsid [N], spike [S], and membrane [M] proteins). Thus, blocking the internalization pathway by inhibiting MYO5 proteins which could be an effective target for coronavirus disease 2019 (COVID-19) treatment. The functional annotations of the host-pathogen interaction (HPI) network were found to be closely associated with hypoxia and thrombotic conditions, confirming the vulnerability and severity of infection. We also screened CpG islands in Nsp1 and N conferring the ability of SARS-CoV-2 to enter and trigger zinc antiviral protein (ZAP) activity inside the host cell.

IMPORTANCE In the current study, we presented a global view of mutational pattern observed in SARS-CoV-2 virus transmission. This provided a who-infect-whom geographical model since the early pandemic. This is hitherto the most comprehensive comparative genomics analysis of full-length genomes for co-mutations at different geographical regions especially in U.S. strains. Compositional structural biology results suggested that mutations have a balance of opposing forces affecting pathogenicity suggesting that only a few mutations are effective at the translation level. Novel HPI analysis and CpG predictions elucidate the proof of concept of hypoxia and thrombotic conditions in several patients. Thus, the current study focuses the understanding of population-specific variations attributing a high rate of SARS-CoV-2 infections in specific geographical regions which may eventually be vital for the most severely affected countries and regions for sharp development of custom-made vindication strategies.

Dynamics of DNA Methylation and Its Functions in Plant Growth and Development

Epigenetic modifications in DNA bases and histone proteins play important roles in the regulation of gene expression and genome stability. Chemical modification of DNA base (e.g., addition of a methyl group at the fifth carbon of cytosine residue) switches on/off the gene expression during developmental process and environmental stresses. The dynamics of DNA base methylation depends mainly on the activities of the writer/eraser guided by non-coding RNA (ncRNA) and regulated by the developmental/environmental cues. De novo DNA methylation and active demethylation activities control the methylation level and regulate the gene expression. Identification of ncRNA involved in de novo DNA methylation, increased DNA methylation proteins guiding DNA demethylase, and methylation monitoring sequence that helps maintaining a balance between DNA methylation and demethylation is the recent developments that may resolve some of the enigmas. Such discoveries provide a better understanding of the dynamics/functions of DNA base methylation and epigenetic regulation of growth, development, and stress tolerance in crop plants. Identification of epigenetic pathways in animals, their existence/orthologs in plants, and functional validation might improve future strategies for epigenome editing toward climate-resilient, sustainable agriculture in this era of global climate change. The present review discusses the dynamics of DNA methylation (cytosine/adenine) in plants, its functions in regulating gene expression under abiotic/biotic stresses, developmental processes, and genome stability.

Recommendations for the introduction of metagenomic high-throughput sequencing in clinical virology, part I: Wet lab procedure

Metagenomic high-throughput sequencing (mHTS) is a hypothesis-free, universal pathogen detection technique for determination of the DNA/RNA sequences in a variety of sample types and infectious syndromes. mHTS is still in its early stages of translating into clinical application. To support the development, implementation and standardization of mHTS procedures for virus diagnostics, the European Society for Clinical Virology (ESCV) Network on Next-Generation Sequencing (ENNGS) has been established. The aim of ENNGS is to bring together professionals involved in mHTS for viral diagnostics to share methodologies and experiences, and to develop application recommendations. This manuscript aims to provide practical recommendations for the wet lab procedures necessary for implementation of mHTS for virus diagnostics and to give recommendations for development and validation of laboratory methods, including mHTS quality assurance, control and quality assessment protocols.

Use of a small DNA virus model to investigate mechanisms of CpG dinucleotide-induced attenuation of virus replication

Suppression of the CpG dinucleotide is widespread in RNA viruses infecting vertebrates and plants, and in the genomes of retroviruses and small mammalian DNA viruses. The functional basis for CpG suppression in the latter was investigated through the construction of mutants of the parvovirus, minute virus of mice (MVM) with increased CpG or TpA dinucleotides in the VP gene. CpG-high mutants displayed extraordinary attenuation in A9 cells compared to wild-type MVM (>six logs), while TpA elevation showed no replication effect. Attenuation was independent of Toll-like receptor 9 and STING-mediated DNA recognition pathways and unrelated to effects on translation efficiency. While translation from codon-optimized VP RNA was enhanced in a cell-free assay, MVM containing this sequence was highly attenuated. Further mutational analysis indicated that this arose through its increased numbers of CpG dinucleotides (7→70) and separately from its increased G+C content (42.3→57.4 %), which independently attenuated replication. CpG-high viruses showed impaired NS mRNA expression by qPCR and reduced NS and particularly VP protein expression detected by immunofluorescence and replication in A549 cells, effects reversed in zinc antiviral protein (ZAP) knockout cells, even though nuclear relocalization of VP remained defective. The demonstrated functional basis for CpG suppression in MVM and potentially other small DNA viruses and the observed intolerance of CpGs in coding sequences, even after codon optimization, has implications for the use of small DNA virus vectors in gene therapy and immunization.

The Roles of Host 5-Methylcytosine RNA Methyltransferases during Viral Infections

Eukaryotic 5-methylcytosine RNA methyltransferases catalyze the transfer of a methyl group to the fifth carbon of a cytosine base in RNA sequences to produce 5-methylcytosine (m5C). m5C RNA methyltransferases play a crucial role in the maintenance of functionality and stability of RNA. Viruses have developed a number of strategies to suppress host innate immunity and ensure efficient transcription and translation for the replication of new virions. One such viral strategy is to use host m5C RNA methyltransferases to modify viral RNA and thus to affect antiviral host responses. Here, we summarize the latest findings concerning the roles of m5C RNA methyltransferases, namely, NOL1/NOP2/SUN domain (NSUN) proteins and DNA methyltransferase 2/tRNA methyltransferase 1 (DNMT2/TRDMT1) during viral infections. Moreover, the use of m5C RNA methyltransferase inhibitors as an antiviral therapy is discussed.

Evaluation of un-methylated DNA enrichment in sequencing of African swine fever virus complete genome

African swine fever is a febrile hemorrhagic fever disease that is caused by the African swine fever virus (ASFV) and is lethal for domestic pigs and wild boar. ASFV also infects soft ticks of the genus Ornithodoros, some species of which can act as a vector for ASFV. Whole genome sequencing of ASFV is a challenge because, due to the size difference of the host genome versus the viral genome, the higher proportion of host versus virus DNA fragments renders the virus sequencing poorly efficient. A novel approach of DNA enrichment, based on the separation of methylated and un-methylated DNA, has been reported but without an evaluation of its efficacy. In this study, the efficiency of the un-methylated DNA enrichment protocol was evaluated for pig and tick samples infected by ASFV. As expected, fewer reads corresponding to ASFV were found in the methylated fraction compared to the un-methylated fraction. However, the sequencing coverage of the un-methylated fraction was not improved compared to the untreated DNA. In our hands, the ASFV DNA enrichment was inefficient for tick samples and very limited for pig samples. This enrichment process represents extra work and cost without a significant improvement of ASFV genome coverage. The efficiency of this enrichment approach and the cost/benefit ratio are discussed.

The DNA methylation landscape of giant viruses

DNA methylation is an important epigenetic mark that contributes to various regulations in all domains of life. Giant viruses are widespread dsDNA viruses with gene contents overlapping the cellular world that also encode DNA methyltransferases. Yet, virtually nothing is known about the methylation of their DNA. Here, we use single-molecule real-time sequencing to study the complete methylome of a large spectrum of giant viruses. We show that DNA methylation is widespread, affecting 2/3 of the tested families, although unevenly distributed. We also identify the corresponding viral methyltransferases and show that they are subject to intricate gene transfers between bacteria, viruses and their eukaryotic host. Most methyltransferases are conserved, functional and under purifying selection, suggesting that they increase the viruses' fitness. Some virally encoded methyltransferases are also paired with restriction endonucleases forming Restriction-Modification systems. Our data suggest that giant viruses' methyltransferases are involved in diverse forms of virus-pathogens interactions during coinfections.

CpG-creating mutations are costly in many human viruses

Mutations can occur throughout the virus genome and may be beneficial, neutral or deleterious. We are interested in mutations that yield a C next to a G, producing CpG sites. CpG sites are rare in eukaryotic and viral genomes. For the eukaryotes, it is thought that CpG sites are rare because they are prone to mutation when methylated. In viruses, we know less about why CpG sites are rare. A previous study in HIV suggested that CpG-creating transition mutations are more costly than similar non-CpG-creating mutations. To determine if this is the case in other viruses, we analyzed the allele frequencies of CpG-creating and non-CpG-creating mutations across various strains, subtypes, and genes of viruses using existing data obtained from Genbank, HIV Databases, and Virus Pathogen Resource. Our results suggest that CpG sites are indeed costly for most viruses. By understanding the cost of CpG sites, we can obtain further insights into the evolution and adaptation of viruses.

Genomic and phylogenetic analysis of two guinea pig adenovirus strains recovered from archival lung tissue

Next generation sequencing was used to determine the whole genome sequence for two different strains of guinea pig adenovirus (GPAdV) detected in association with outbreaks of pneumonia in Australia in 1996, and in Germany in 1997 using total DNA extracted from infected archival frozen lung tissue as a template. The length of the determined genomic sequences was 37,031 bp and 37,070 bp, respectively. The nucleotide composition showed a relatively high content of guanine + cytosine (G + C) of 62 %. The 99.6 % nucleotide identity between the two sequenced viruses suggests that they may represent variants of the same genotype. The GPAdV genome exhibits the genomic features of a typical mastadenovirus with at least 32 open reading frames identified. Five novel open reading frames were found at the right end of the genomic sequence. One of them maps to the predicted E3 region and encodes a putative CR1 protein, two map to the E4 region, and two map to the l strand of L1 and L3, respectively. Our phylogenetic analysis of whole genome sequences showed that among the mammalian AdV species described to date, GPAdV is most closely related to MAdV-2 The characterization of this mastadenovirus species offers an opportunity to develop a new small animal model to study mammalian adenovirus pathogenesis.

Methylation of CpG 5962 in L1 of the human papillomavirus 16 genome as a potential predictive marker for viral persistence: A prospective large cohort study using cervical swab samples

Several studies have demonstrated that the viral genome can be methylated by the host cell during progression from persistent infection to cervical cancer. The aim of this study was to investigate whether methylation at a specific site could predict the development of viral persistence and whether viral load shows a correlation with specific methylation patterns. HPV16‐positive samples from women aged 20–29 years (n = 99) with a follow‐up time of 13 years, were included from a Danish cohort comprising 11 088 women. Viral load was measured by real‐time PCR and methylation status was determined for 39 CpG sites in the upstream regulatory region (URR), E6/E7, and L1 region of HPV16 by next‐generation sequencing. Participants were divided into two groups according to whether they were persistently (≥ 24 months) or transiently HPV16 infected. The general methylation status was significantly different between women with a persistent and women with a transient infection outcome (P = .025). One site located in L1 (nt. 5962) was statistically significantly (P = .00048) different in the methylation status after correction using the Holm‐Sidak method (alpha = 0.05). Correlation analyses of samples from HPV16 persistently infected women suggest that methylation is higher although viral load is lower. This study indicates that methylation at position 5962 of the HPV16 genome within the L1 gene might be a predictive marker for the development of a persistent HPV16 infection.

Methylation Status of the Adeno-Associated Virus Type 2 (AAV2)

To analyze the methylation status of wild-type adeno-associated virus type 2 (AAV2), bisulfite PCR sequencing (BPS) of the packaged viral genome and its integrated form was performed and 262 of the total 266 CG dinucleotides (CpG) were mapped. In virion-packaged DNA, the ratio of the methylated cytosines ranged between 0⁻1.7%. In contrast, the chromosomally integrated AAV2 genome was hypermethylated with an average of 76% methylation per CpG site. The methylation level showed local minimums around the four known AAV2 promoters. To study the effect of methylation on viral rescue and replication, the replication initiation capability of CpG methylated and non-CpG methylated AAV DNA was compared. The in vitro hypermethylation of the viral genome does not inhibit its rescue and replication from a plasmid transfected into cells. This insensitivity of the viral replicative machinery to methylation may permit the rescue of the integrated heavily methylated AAV genome from the host's chromosomes.

The effect of different dietary levels of DL-methionine and DL-hydroxy analogue on the antioxidant status of young turkeys infected with the haemorrhagic enteritis virus

BACKGROUND

The results of experiments involving broiler chickens and turkeys indicate that increased dietary methionine (Met) levels may improve the antioxidant protection of tissues in fast-growing birds. This is an important consideration since viral infections induce oxidative stress. The aim of this study was to verify the hypothesis that turkey diets with increased Met content can suppress oxidation processes induced by infection caused by the haemorrhagic enteritis virus (HEV), and that the noted effect is determined by the chemical form of this amino acid: DL-methionine (DLM) or DL-hydroxy analogue of Met (MHA).

RESULTS

Dietary Met content above 40% higher than the level recommended by the NRC (1994) intensified lipid peroxidation in the small intestine, leading to an increase in malondialdehyde (MDA) and lipid peroxide (LOOH) levels, but it also stimulated antioxidant mechanisms in the blood and liver of turkeys infected with HEV. In comparison with DLM, MHA contributed to more severe symptoms of oxidative stress, such as elevated MDA levels in the intestines, and a decrease in glutathione peroxidase (GPx) activity and ferric-reducing ability of plasma (FRAP).

CONCLUSIONS

In HEV-infected turkeys, diets with increased Met content did not exert a clear antioxidant effect, which was noted in uninfected birds. The prooxidant activity of Met observed in the small intestinal wall was suppressed in the blood and liver of turkeys, most likely due to intensified synthesis of uric acid and glutathione. In comparison with MHA, DLM had a more beneficial influence on the analysed parameters of the redox status in the small intestine, blood and liver of turkeys.

Methylation in host and viral genes as marker of aggressiveness in cervical lesions: Analysis in 543 unscreened women

OBJECTIVE

The present study aimed to evaluate the association between altered methylation and histologically confirmed high grade cervical intraepithelial neoplasia (hgCIN).

METHODS

Methylation levels in selected host (CADM1, MAL, DAPK1) and HPV (L1_I, L1_II, L2) genes were measured by pyrosequencing in DNA samples obtained from 543 women recruited in Curitiba (Brazil), 249 with hgCIN and 294 without cervical lesions. Association of methylation status with hgCIN was estimated by Odds Ratio (OR) with 95% confidence interval (CI).

RESULTS

The mean methylation level increased with severity of the lesion in the host and viral genes (p-trend < 0.05), with the exception of L1_II region (p-trend = 0.075). Positive association was found between methylation levels for host genes and CIN2 and CIN3 lesions respectively [CADM1: OR 4.17 (95%CI 2.03-8.56) and OR 9.54 (95%CI 4.80-18.97); MAL: OR 5.98 (95%CI 2.26-15.78) and OR 22.66 (95%CI 9.21-55.76); DAPK1: OR 3.37 (95%CI 0.93-12.13) and OR 6.74 (95%CI 1.92-23.64)]. Stronger risk estimates were found for viral genes [L1_I: OR 10.74 (95%CI 2.66-43.31) and OR 15.00 (95%CI 3.00-74.98); L1_II: OR 73.18 (95%CI 4.07-1315.94) and OR 32.50 (95%CI 3.86-273.65); L2: OR 4.73 (95%CI 1.55-14.44) and OR 10.62 (95%CI 2.60-43.39)]. The cumulative effect of the increasing number of host and viral methylated genes was associated with the risk of CIN2 and CIN3 lesions (p-trend < 0.001).

CONCLUSIONS

Our results, empowered by a wide cervical sample series with a large number of hgCIN, supported the role of methylation as marker of aggressiveness.

Cellular Antisilencing Elements Support Transgene Expression from Herpes Simplex Virus Vectors in the Absence of Immediate Early Gene Expression

Inactivation of all herpes simplex virus (HSV) immediate early (IE) genes to eliminate vector cytotoxicity results in rapid silencing of the viral genome, similar to the establishment of HSV latency. We recently reported that silencing of a nonviral reporter cassette could be overcome in nonneuronal cells by positioning the cassette in the viral latency (LAT) locus between resident chromatin boundary elements. Here, we tested the abilities of the chicken hypersensitive site 4 insulator and the human ubiquitous chromatin opening element A2UCOE to promote transgene expression from an IE-gene-inactivated HSV vector. We found that A2UCOE was particularly active in nonneuronal cells and reduced reporter promoter occupancy by a repressive histone mark. We determined whether multiple transgenes could be expressed under the control of different promoters from different loci of the same virus. The results showed abundant coexpression of LAT-embedded and A2UCOE-flanked genes in nonneuronal cells. In addition, a third reporter gene without known protective elements was active in cultured rat sensory neurons. These findings indicate that cellular antisilencing sequences can contribute to the expression of multiple genes from separate promoters in fully IE gene-disabled HSV vectors, providing an opportunity for therapeutic applications requiring mutually independent expression of different gene products from a single vector.IMPORTANCE Gene therapy has now entered a phase of development in which a growing number of recessive single gene defects can be successfully treated by vector-mediated introduction of a wild-type copy of the gene into the appropriate tissue. However, many disease conditions, such as neurodegeneration, cancer, and inflammatory processes, are more complex, requiring either multiple gene corrections or provision of coordinated gene activities to achieve a therapeutic outcome. Although herpes simplex virus (HSV) vectors have the capacity to meet this need, the challenge has been to genetically engineer the HSV genome in a manner to prevent expression of any viral genes while retaining the ability to express multiple therapeutic transgenes under independent transcriptional control. Here, we show that non-HSV insulator elements can be applied to retain at least transient transgene activity from multiple viral loci, thereby opening the door for more complex gene therapy applications in the future.

Epstein-Barr virus DNA modulates regulatory T-cell programming in addition to enhancing interleukin-17A production via Toll-like receptor 9

Infection with the Epstein-Barr virus (EBV) has been associated with several autoimmune diseases including rheumatoid arthritis (RA). We have previously reported that DNA from this virus enhances production of the pro-autoimmune interleukin 17A (IL-17A) in mice. In this study we assessed the effect of EBV DNA on regulatory T cell programming and examined whether it mediated its effects via Toll-like receptor 9 (TLR9) in mice; moreover, we evaluated whether EBV DNA in humans had similar effects to those seen in mice. For this purpose, we assessed the linearity of the correlation between EBV DNA and IL-17A levels in RA subjects and matched controls. A modulatory effect for the viral DNA was observed for regulatory T cell markers with an inhibitory effect observed for CTLA4 expression in the EBV DNA-treated mice. To examine whether TLR9 mediated the detection of EBV DNA and enhancement of IL-17A production, mouse peripheral blood mononuclear cells were treated with the DNA in the presence or absence of the TLR9 inhibitor ODN 2088. Subsequently, IL-17A production from these cells was assessed. Treatment with the TLR9 inhibitor resulted in a significant decrease in IL-17A production indicating that TLR9 is involved in this pathway. In human subjects, examining the linearity of the correlation between EBV DNA and IL-17A levels in RA subjects showed a propensity for linearity that was not observed in controls. Our data thus indicates that EBV DNA itself acts as a modulator of the Th17 compartment as well as that of regulatory T cell mechanisms. The involvement of TLR9 in the EBV DNA-triggered induction of IL-17A suggests therapeutic targeting of this endosomal receptor in EBV positive subjects with an autoimmune flare-up or possibly for prophylactic purposes.

The analysis of genome composition and codon bias reveals distinctive patterns between avian and mammalian circoviruses which suggest a potential recombinant origin for Porcine circovirus 3

Members of the genus Circovirus are host-specific viruses, which are totally dependent on cell machinery for their replication. Consequently, certain mimicry of the host genome features is expected to maximize cellular replicative system exploitation and minimize the recognition by the innate immune system. In the present study, the analysis of several genome composition and codon bias parameters of circoviruses infecting avian and mammalian species demonstrated the presence of quite distinctive patterns between the two groups. Remarkably, a higher deviation from the expected values based only on mutational patterns was observed for mammalian circoviruses both at dinucleotide and codon levels. Accordingly, a stronger selective pressure was estimated to shape the genome of mammalian circoviruses, particularly in the Cap encoding gene, compared to avian circoviruses. These differences could be attributed to different physiological and immunological features of the two host classes and suggest a trade-off between a tendency to optimize the capsid protein translation while minimizing the recognition of the genome and the transcript molecules. Interestingly, the recently identified Porcine circovirus 3 (PCV-3) had an intermediate pattern in terms of genome composition and codon bias. Particularly, its Rep gene appeared closely related to other mammalian circoviruses (especially bat circoviruses) while the Cap gene more closely resembled avian circoviruses. These evidences, coupled with the high selective forces apparently modelling the PCV-3 Cap gene composition, suggest the potential recombinant origin, followed or preceded by a host jump, of this virus.

UBXN3B positively regulates STING-mediated antiviral immune responses

The ubiquitin regulatory X domain-containing proteins (UBXNs) are likely involved in diverse biological processes. Their physiological functions, however, remain largely unknown. Here we present physiological evidence that UBXN3B positively regulates stimulator-of-interferon genes (STING) signaling. We employ a tamoxifen-inducible Cre-LoxP approach to generate systemic Ubxn3b knockout in adult mice as the Ubxn3b-null mutation is embryonically lethal. Ubxn3b^-/-, like Sting^-/- mice, are highly susceptible to lethal herpes simplex virus 1 (HSV-1) and vesicular stomatitis virus (VSV) infection, which is correlated with deficient immune responses when compared to Ubxn3b^+/+ littermates. HSV-1 and STING agonist-induced immune responses are also reduced in several mouse and human Ubxn3b^-/- primary cells. Mechanistic studies demonstrate that UBXN3B interacts with both STING and its E3 ligase TRIM56, and facilitates STING ubiquitination, dimerization, trafficking, and consequent recruitment and phosphorylation of TBK1. These results provide physiological evidence that links the UBXN family with antiviral immune responses.

Assessment of viral methylation levels for high risk HPV types by newly designed consensus primers PCR and pyrosequencing

Background

Measuring viral DNA methylation in human papillomavirus (HPV) infected women showed promise for accurate detection of high-grade cervical lesions and cancer. Methylation status has been widely investigated for HPV16, sporadically for other HPV types.

Methods

Objective of this methodological study was to set up molecular methods to test the methylation levels in the twelve oncogenic HPV types by pyrosequencing, minimizing the number of HPV type-specific PCR protocols. Target CpGs were selected on the HPV L1 (two regions, L1 I and L1 II) and L2 genes.

Study samples included DNA stored at Turin, Italy, purified by cervical cells collected in Standard Transport Medium or PreservCyt from women who participated in two studies (N = 126 and 140) nested within the regional organized screening programme.

PCR consensus primers were designed by PyroMark Assay Design software to be suitable for amplification of many different oncogenic HPV types.

Results

Generation of consensus primers was successful for L1 I and II regions, unsuccessful for L2 region, for which HPV type-specific primers remained necessary. The difference between replicated tests on the same sample was ≤4% in 88%, 77% and 91% of cases when targeting the L1 I, L1 II and L2 regions, respectively. The corresponding intra-class correlation coefficients (ICC) were 0.94, 0.87 and 0.97 respectively. When comparing methylation measures based on consensus and type-specific primers, ICC was 0.97 for the L1 I region and 0.99 the for L1 II region.

Conclusions

The proposed protocols, applying consensus primers suitable to amplify the oncogenic HPV types and minimize the number of PCR reactions, represent a promising tool to quantify viral methylation in women positive for any high risk HPV type.

Impact

Potential application of these methylation protocols in screening settings can be explored to identify women with high probability of progression to high grade lesions.

DNA Tumor Virus Regulation of Host DNA Methylation and Its Implications for Immune Evasion and Oncogenesis

Viruses have evolved various mechanisms to evade host immunity and ensure efficient viral replication and persistence. Several DNA tumor viruses modulate host DNA methyltransferases for epigenetic dysregulation of immune-related gene expression in host cells. The host immune responses suppressed by virus-induced aberrant DNA methylation are also frequently involved in antitumor immune responses. Here, we describe viral mechanisms and virus–host interactions by which DNA tumor viruses regulate host DNA methylation to evade antiviral immunity, which may contribute to the generation of an immunosuppressive microenvironment during cancer development. Recent trials of immunotherapies have shown promising results to treat multiple cancers; however, a significant number of non-responders necessitate identifying additional targets for cancer immunotherapies. Thus, understanding immune evasion mechanisms of cancer-causing viruses may provide great insights for reversing immune suppression to prevent and treat associated cancers.

Biology of Porcine Parvovirus (Ungulate parvovirus 1)

Porcine parvovirus (PPV) is among the most important infectious agents causing infertility in pigs. Until recently, it was thought that the virus had low genetic variance, and that prevention of its harmful effect on pig fertility could be well-controlled by vaccination. However, at the beginning of the third millennium, field observations raised concerns about the effectiveness of the available vaccines against newly emerging strains. Subsequent investigations radically changed our view on the evolution and immunology of PPV, revealing that the virus is much more diverse than it was earlier anticipated, and that some of the “new” highly virulent isolates cannot be neutralized effectively by antisera raised against “old” PPV vaccine strains. These findings revitalized PPV research that led to significant advancements in the understanding of early and late viral processes during PPV infection. Our review summarizes the recent results of PPV research and aims to give a comprehensive update on the present understanding of PPV biology.

Roles of APOBEC3A and APOBEC3B in Human Papillomavirus Infection and Disease Progression

The apolipoprotein B messenger RNA-editing, enzyme-catalytic, polypeptide-like 3 (APOBEC3) family of cytidine deaminases plays an important role in the innate immune response to viral infections by editing viral genomes. However, the cytidine deaminase activity of APOBEC3 enzymes also induces somatic mutations in host genomes, which may drive cancer progression. Recent studies of human papillomavirus (HPV) infection and disease outcome highlight this duality. HPV infection is potently inhibited by one family member, APOBEC3A. Expression of APOBEC3A and APOBEC3B is highly elevated by the HPV oncoproteins E6 and E7 during persistent virus infection and disease progression. Furthermore, there is a high prevalence of APOBEC3A and APOBEC3B mutation signatures in HPV-associated cancers. These findings suggest that induction of an APOBEC3-mediated antiviral response during HPV infection may inadvertently contribute to cancer mutagenesis and virus evolution. Here, we discuss current understanding of APOBEC3A and APOBEC3B biology in HPV restriction, evolution, and associated cancer mutagenesis.

Novel bat adenoviruses with low G+C content shed new light on the evolution of adenoviruses

Bats have been reported to carry diverse adenoviruses. However, most bat adenoviruses have been identified on the basis of partial genome sequences, and knowledge on the evolution of bat adenoviruses remains limited. In this study, we isolated and characterized four novel adenoviruses from two distinct bat species, and their full-length genomes were sequenced. Sequence analysis revealed that these isolates represented three distinct species of the genus Mastadenovirus. However, all isolates had an exceptionally low G+C content and relatively short genomes compared with other known mastadenoviruses. We further analysed the relationships among the G+C content, 5'-C-phosphate-G-3' (CpG) representation and genome size in the family Adenoviridae. Our results revealed that the CpG representation in adenoviral genomes depends primarily on the level of methylation, and the genome size displayed significant positive correlations with both G+C content and CpG representation. Since ancestral adenoviruses are believed to have contained short genomes, those probably had a low G+C content, similar to the genomes of these bat strains. Our results suggest that bats are important natural reservoirs for adenoviruses and play important roles in the evolution of adenoviruses.

Human Herpesvirus 6B Induces Hypomethylation on Chromosome 17p13.3, Correlating with Increased Gene Expression and Virus Integration

Human herpesvirus 6B (HHV-6B) is a neurotropic betaherpesvirus that achieves latency by integrating its genome into host cell chromosomes. Several viruses can induce epigenetic modifications in their host cells, but no study has investigated the epigenetic modifications induced by HHV-6B. This study analyzed methylation with an Illumina 450K array, comparing HHV-6B-infected and uninfected Molt-3 T cells 3 days postinfection. Bisulfite pyrosequencing was used to validate the Illumina results and to investigate methylation over time in vitro Expression of genes was investigated using quantitative PCR (qPCR), and virus integration was investigated with PCR. A total of 406 CpG sites showed a significant HHV-6B-induced change in methylation in vitro Remarkably, 86% (351/406) of these CpGs were located <1 Mb from chromosomal ends and were all hypomethylated in virus-infected cells. This was most evident at chromosome 17p13.3, where HHV-6B had induced CpG hypomethylation after 2 days of infection, possibly through TET2, which was found to be upregulated by the virus. In addition, virus-induced cytosine hydroxymethylation was observed. Genes located in the hypomethylated region at 17p13.3 showed significantly upregulated expression in HHV-6B-infected cells. A temporal experiment revealed HHV-6B integration in Molt-3 cell DNA 3 days after infection. The telomere at 17p has repeatedly been described as an integration site for HHV-6B, and we show for the first time that HHV-6B induces hypomethylation in this region during acute infection, which may play a role in the integration process, possibly by making the DNA more accessible.IMPORTANCE The ability to establish latency in the host is a hallmark of herpesviruses, but the mechanisms differ. Human herpesvirus 6B (HHV-6B) is known to establish latency through integration of its genome into the telomeric regions of host cells, with the ability to reactivate. Our study is the first to show that HHV-6B specifically induces hypomethylated regions close to the telomeres and that integrating viruses may use the host methylation machinery to facilitate their integration process. The results from this study contribute to knowledge of HHV-6B biology and virus-host interaction. This in turn will lead to further progress in our understanding of the underlying mechanisms by which HHV-6B contributes to pathological processes and may have important implications in both disease prevention and treatment.

K-mer Content, Correlation, and Position Analysis of Genome DNA Sequences for the Identification of Function and Evolutionary Features

In genome analysis, k-mer-based comparison methods have become standard tools. However, even though they are able to deliver reliable results, other algorithms seem to work better in some cases. To improve k-mer-based DNA sequence analysis and comparison, we successfully checked whether adding positional resolution is beneficial for finding and/or comparing interesting organizational structures. A simple but efficient algorithm for extracting and saving local k-mer spectra (frequency distribution of k-mers) was developed and used. The results were analyzed by including positional information based on visualizations as genomic maps and by applying basic vector correlation methods. This analysis was concentrated on small word lengths (1 ≤ k ≤ 4) on relatively small viral genomes of Papillomaviridae and Herpesviridae, while also checking its usability for larger sequences, namely human chromosome 2 and the homologous chromosomes (2A, 2B) of a chimpanzee. Using this alignment-free analysis, several regions with specific characteristics in Papillomaviridae and Herpesviridae formerly identified by independent, mostly alignment-based methods, were confirmed. Correlations between the k-mer content and several genes in these genomes have been found, showing similarities between classified and unclassified viruses, which may be potentially useful for further taxonomic research. Furthermore, unknown k-mer correlations in the genomes of Human Herpesviruses (HHVs), which are probably of major biological function, are found and described. Using the chromosomes of a chimpanzee and human that are currently known, identities between the species on every analyzed chromosome were reproduced. This demonstrates the feasibility of our approach for large data sets of complex genomes. Based on these results, we suggest k-mer analysis with positional resolution as a method for closing a gap between the effectiveness of alignment-based methods (like NCBI BLAST) and the high pace of standard k-mer analysis.

Immunoassay and molecular methods to investigate DNA methylation changes in peripheral blood mononuclear cells in HIV infected patients on cART

This study aimed to investigate the influence of antiretroviral therapy on methylation markers, in a group of HIV infected, heavily treated patients. Immune and molecular methods were used to investigate potential changes in methylation profile in DNA isolated from peripheral blood mononuclear cells collected from antiretroviral-experienced HIV infected patients and healthy controls. The percentage of 5-methylcytosine was inversely correlated with proviral DNA and active replication while DNMT1 (p = 0.01) and DNMT3A (p = 0.004) independently correlated with active viral replication. DNMT3A expression increased with total treatment duration (p = 0.03), number of antiretroviral drugs ever used (p = 0.003), and cumulative exposure to protease inhibitors (p = 0.02) even in currently HIV undetectable patients.

Genome-Wide Estimation of the Spontaneous Mutation Rate of Human Adenovirus 5 by High-Fidelity Deep Sequencing

Rates of spontaneous mutation determine the ability of viruses to evolve, infect new hosts, evade immunity and undergo drug resistance. Contrarily to RNA viruses, few mutation rate estimates have been obtained for DNA viruses, because their high replication fidelity implies that new mutations typically fall below the detection limits of Sanger and standard next-generation sequencing. Here, we have used a recently developed high-fidelity deep sequencing technique (Duplex Sequencing) to score spontaneous mutations in human adenovirus 5 under conditions of minimal selection. Based on >200 single-base spontaneous mutations detected throughout the entire viral genome, we infer an average mutation rate of 1.3 × 10⁻⁷ per base per cell infection cycle. This value is similar to those of other, large double-stranded DNA viruses, but an order of magnitude lower than those of single-stranded DNA viruses, consistent with the possible action of post-replicative repair. Although the mutation rate did not vary strongly along the adenovirus genome, we found several sources of mutation rate heterogeneity. First, two regions mapping to transcription units L3 and E1B-IVa2 were significantly depleted for mutations. Second, several point insertions/deletions located within low-complexity sequence contexts appeared recurrently, suggesting mutational hotspots. Third, mutation probability increased at GpC dinucleotides. Our findings suggest that host factors may influence the distribution of spontaneous mutations in human adenoviruses and potentially other nuclear DNA viruses.

Next-generation sequencing has provided a powerful tool for studying microbial genetic diversity but suffers from relatively low per-base accuracy, limiting our ability to detect low-frequency polymorphisms and spontaneous mutations. However, this limitation has been solved recently by the development of high-fidelity deep sequencing techniques. Taking advantage of these advancements, here we provide the first unbiased genome-wide characterization of the rate of spontaneous mutation of a human DNA virus (adenovirus 5) under controlled laboratory conditions. The adenovirus genome shows a relatively low mutation rate, consistent with high replication fidelity and the action of post-replicative repair. We also found evidence for mutation rate heterogeneities and regions of genetic instability in the viral genome. Together with previous reports, our findings indicate that DNA viruses with large double-stranded genomes mutate significantly slower than those with small single-stranded genomes.

Functions and Malfunctions of Mammalian DNA-Cytosine Deaminases

The AID/APOBEC family enzymes convert cytosines in single-stranded DNA to uracils, causing base substitutions and strand breaks. They are induced by cytokines produced during the body's inflammatory response to infections, and they help combat infections through diverse mechanisms. AID is essential for the maturation of antibodies and causes mutations and deletions in antibody genes through somatic hypermutation (SHM) and class-switch recombination (CSR) processes. One member of the APOBEC family, APOBEC1, edits mRNA for a protein involved in lipid transport. Members of the APOBEC3 subfamily in humans (APOBEC3A, APOBEC3B, APOBEC3C, APOBEC3D, APOBEC3F, APOBEC3G, and APOBEC3H) inhibit infections of viruses such as HIV-1, HBV, and HCV, and retrotransposition of endogenous retroelements through mutagenic and nonmutagenic mechanisms. There is emerging consensus that these enzymes can cause mutations in the cellular genome at replication forks or within transcription bubbles depending on the physiological state of the cell and the phase of the cell cycle during which they are expressed. We describe here the state of knowledge about the structures of these enzymes, regulation of their expression, and both the advantageous and deleterious consequences of their expression, including carcinogenesis. We highlight similarities among them and present a holistic view of their regulation and function.

Functional Upregulation of the DNA Cytosine Deaminase APOBEC3B by Polyomaviruses

The APOBEC3 family of DNA cytosine deaminases has important roles in innate immunity and cancer. It is unclear how DNA tumor viruses regulate these enzymes and how these interactions, in turn, impact the integrity of both the viral and cellular genomes. Polyomavirus (PyVs) are small DNA pathogens that contain oncogenic potentials. In this study, we examined the effects of PyV infection on APOBEC3 expression and activity. We demonstrate that APOBEC3B is specifically upregulated by BK polyomavirus (BKPyV) infection in primary kidney cells and that the upregulated enzyme is active. We further show that the BKPyV large T antigen, as well as large T antigens from related polyomaviruses, is alone capable of upregulating APOBEC3B expression and activity. Furthermore, we assessed the impact of A3B on productive BKPyV infection and viral genome evolution. Although the specific knockdown of APOBEC3B has little short-term effect on productive BKPyV infection, our informatics analyses indicate that the preferred target sequences of APOBEC3B are depleted in BKPyV genomes and that this motif underrepresentation is enriched on the nontranscribed stand of the viral genome, which is also the lagging strand during viral DNA replication. Our results suggest that PyV infection upregulates APOBEC3B activity to influence virus sequence composition over longer evolutionary periods. These findings also imply that the increased activity of APOBEC3B may contribute to PyV-mediated tumorigenesis.

IMPORTANCE

Polyomaviruses (PyVs) are a group of emerging pathogens that can cause severe diseases, including cancers in immunosuppressed individuals. Here we describe the finding that PyV infection specifically induces the innate immune DNA cytosine deaminase APOBEC3B. The induced APOBEC3B enzyme is fully functional and therefore may exert mutational effects on both viral and host cell DNA. We provide bioinformatic evidence that, consistent with this idea, BK polyomavirus genomes are depleted of APOBEC3B-preferred target motifs and enriched for the corresponding predicted reaction products. These data imply that the interplay between PyV infection and APOBEC proteins may have significant impact on both viral evolution and virus-induced tumorigenesis.

Poly(ADP-ribose) polymerase-1 silences retroviruses independently of viral DNA integration or heterochromatin formation

PARP-1 silences retrotransposons in Drosophila, through heterochromatin maintenance, and integrated retroviruses in chicken. Here, we determined the role of viral DNA integration and cellular heterochromatin in PARP-1-mediated retroviral silencing using HIV-1-derived lentiviral vectors and Rous-associated virus type 1 (RAV-1) as models. Analysis of the infection of PARP-1 knockout and control cells with HIV-1 harbouring WT integrase, in the presence or absence of an integrase inhibitor, or catalytic-dead mutant integrase indicated that silencing does not require viral DNA integration. The mechanism involves the catalytic activity of histone deacetylases but not that of PARP-1. In contrast to Drosophila, lack of PARP-1 in avian cells did not affect chromatin compaction globally or at the RAV-1 provirus, or the cellular levels of histone H3 N-terminal acetylated or Lys27 trimethylated, as indicated by micrococcal nuclease accessibility and immunoblot assays. Therefore, PARP-1 represses retroviruses prior to viral DNA integration by mechanisms involving histone deacetylases but not heterochromatin formation.

Epigenetic Segregation of Microbial Genomes from Complex Samples Using Restriction Endonucleases HpaII and McrB

We describe continuing work to develop restriction endonucleases as tools to enrich targeted genomes of interest from diverse populations. Two approaches were developed in parallel to segregate genomic DNA based on cytosine methylation. First, the methyl-sensitive endonuclease HpaII was used to bind non-CG methylated DNA. Second, a truncated fragment of McrB was used to bind CpG methylated DNA. Enrichment levels of microbial genomes can exceed 100-fold with HpaII allowing improved genomic detection and coverage of otherwise trace microbial genomes from sputum. Additionally, we observe interesting enrichment results that correlate with the methylation states not only of bacteria, but of fungi, viruses, a protist and plants. The methods presented here offer promise for testing biological samples for pathogens and global analysis of population methylomes.

ViralEpi v1.0: a high-throughput spectrum of viral epigenomic methylation profiles from diverse diseases

AIMS

To develop a computational resource for viral epigenomic methylation profiles from diverse diseases.

MATERIALS & METHODS

Methylation patterns of Epstein-Barr virus and hepatitis B virus genomic regions are provided as web platform developed using open source Linux-Apache-MySQL-PHP (LAMP) bundle: programming and scripting languages, that is, HTML, JavaScript and PERL.

RESULTS

A comprehensive and integrated web resource ViralEpi v1.0 is developed providing well-organized compendium of methylation events and statistical analysis associated with several diseases. Additionally, it also facilitates 'Viral EpiGenome Browser' for user-affable browsing experience using JavaScript-based JBrowse.

CONCLUSION

This web resource would be helpful for research community engaged in studying epigenetic biomarkers for appropriate prognosis and diagnosis of diseases and its various stages.