Depletion of CpG Dinucleotides in Papillomaviruses and Polyomaviruses: A Role for Divergent Evolutionary Pressures

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

Evolution and diversity of nucleotide and dinucleotide composition in poxviruses

Poxviruses (family Poxviridae) have long dsDNA genomes and infect a wide range of hosts, including insects, birds, reptiles and mammals. These viruses have substantial incidence, prevalence and disease burden in humans and in other animals. Nucleotide and dinucleotide composition, mostly CpG and TpA, have been largely studied in viral genomes because of their evolutionary and functional implications. We analysed here the nucleotide and dinucleotide composition, as well as codon usage bias, of a set of representative poxvirus genomes, with a very diverse host spectrum. After correcting for overall nucleotide composition, entomopoxviruses displayed low overall GC content, no enrichment in TpA and large variation in CpG enrichment, while chordopoxviruses showed large variation in nucleotide composition, no obvious depletion in CpG and a weak trend for TpA depletion in GC-rich genomes. Overall, intergenome variation in dinucleotide composition in poxviruses is largely accounted for by variation in overall genomic GC levels. Nonetheless, using vaccinia virus as a model, we found that genes expressed at the earliest times in infection are more CpG-depleted than genes expressed at later stages. This observation has parallels in betahepesviruses (also large dsDNA viruses) and suggests an antiviral role for the innate immune system (e.g. via the zinc-finger antiviral protein ZAP) in the early phases of poxvirus infection. We also analysed codon usage bias in poxviruses and we observed that it is mostly determined by genomic GC content, and that stratification after host taxonomy does not contribute to explaining codon usage bias diversity. By analysis of within-species diversity, we show that genomic GC content is the result of mutational biases. Poxvirus genomes that encode a DNA ligase are significantly AT-richer than those that do not, suggesting that DNA repair systems shape mutation biases. Our data shed light on the evolution of poxviruses and inform strategies for their genetic manipulation for therapeutic purposes.

Insights into Synonymous Codon Usage Bias in Hepatitis C Virus and Its Adaptation to Hosts

Hepatitis C virus (HCV) is enveloped RNA virus, encoding for a polyprotein that is processed by cellular proteases. The virus is responsible for liver cirrhosis, allograft rejection, and human hepatocellular carcinoma. Based on studies including compositional analysis, odds ratio analysis, parity analysis, skew analysis, relative synonymous codon usage, codon bias, and protein properties, it was evident that codon usage bias in HCV is dependent upon the nucleotide composition. Codon context analysis revealed CTC-CTG as a preferred codon pair. While CGA and CGT codons were rare, none of the codons were rare in HCV-like viruses envisaged in the present study. Many of the preferred codon pairs were valine amino acid-initiated, which possibly infers viral infectivity; hence the role of selection forces appears to act on the HCV genome, which was further validated by neutrality analysis where selection accounted for 87.28%, while mutation accounted for 12.72% force shaping codon usage. Furthermore, codon usage was correlated with the length of the genome. HCV viruses prefer valine-initiated codon pairs, while HCV-like viruses prefer alanine-initiated codon pairs. The HCV host range is very narrow and is confined to only humans and chimpanzees. Based on indices including codon usage correlation analysis, similarity index, and relative codon deoptimization index, it is evident in the study that the chimpanzee is the primary host of the virus. The present study helped elucidate the preferred host for HCV. The information presented in the study paved the way for generating an attenuated vaccine candidate through viral recoding, with finely tuned nucleotide composition and a perfect balance of preferred and rare codons.

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.

Synonymous nucleotide changes drive papillomavirus evolution

Papillomaviruses have been evolving alongside their hosts for at least 450 million years. This review will discuss some of the insights gained into the evolution of this diverse family of viruses. Papillomavirus evolution is constrained by pervasive purifying selection to maximize viral fitness. Yet these viruses need to adapt to changes in their environment, e.g., the host immune system. It has long been known that these viruses evolved a codon usage that doesn't match the infected host. Here we discuss how papillomavirus genomes evolve by acquiring synonymous changes that allow the virus to avoid detection by the host innate immune system without changing the encoded proteins and associated fitness loss. We discuss the implications of studying viral evolution, lifecycle, and cancer progression.

Human genes with relative synonymous codon usage analogous to that of polyomaviruses are involved in the mechanism of polyomavirus nephropathy

Human polyomaviruses (HPyVs) can cause serious and deleterious infections in human. Yet, the molecular mechanism underlying these infections, particularly in polyomavirus nephropathy (PVAN), is not well-defined. In the present study, we aimed to identify human genes with codon usage bias (CUB) similar to that of HPyV genes and explore their potential involvement in the pathogenesis of PVAN. The relative synonymous codon usage (RSCU) values of genes of HPyVs and those of human genes were computed and used for Pearson correlation analysis. The involvement of the identified correlation genes in PVAN was analyzed by validating their differential expression in publicly available transcriptomics data. Functional enrichment was performed to uncover the role of sets of genes. The RSCU analysis indicated that the A- and T-ending codons are preferentially used in HPyV genes. In total, 5400 human genes were correlated to the HPyV genes. The protein-protein interaction (PPI) network indicated strong interactions between these proteins. Gene expression analysis indicated that 229 of these genes were consistently and differentially expressed between normal kidney tissues and kidney tissues from PVAN patients. Functional enrichment analysis indicated that these genes were involved in biological processes related to transcription and in pathways related to protein ubiquitination pathway, apoptosis, cellular response to stress, inflammation and immune system. The identified genes may serve as diagnostic biomarkers and potential therapeutic targets for HPyV associated diseases, especially PVAN.

Regulation of the Innate Immune Response during the Human Papillomavirus Life Cycle

High-risk human papillomaviruses (HR HPVs) are associated with multiple human cancers and comprise 5% of the human cancer burden. Although most infections are transient, persistent infections are a major risk factor for cancer development. The life cycle of HPV is intimately linked to epithelial differentiation. HPVs establish infection at a low copy number in the proliferating basal keratinocytes of the stratified epithelium. In contrast, the productive phase of the viral life cycle is activated upon epithelial differentiation, resulting in viral genome amplification, high levels of late gene expression, and the assembly of virions that are shed from the epithelial surface. Avoiding activation of an innate immune response during the course of infection plays a key role in promoting viral persistence as well as completion of the viral life cycle in differentiating epithelial cells. This review highlights the recent advances in our understanding of how HPVs manipulate the host cell environment, often in a type-specific manner, to suppress activation of an innate immune response to establish conditions supportive of viral replication.

Coevolutionary Analysis Implicates Toll-Like Receptor 9 in Papillomavirus Restriction

Upon infection, DNA viruses can be sensed by pattern recognition receptors (PRRs), leading to the activation of type I and III interferons to block infection. Therefore, viruses must inhibit these signaling pathways, avoid being detected, or both. Papillomavirus virions are trafficked from early endosomes to the Golgi apparatus and wait for the onset of mitosis to complete nuclear entry. This unique subcellular trafficking strategy avoids detection by cytoplasmic PRRs, a property that may contribute to the establishment of infection. However, as the capsid uncoats within acidic endosomal compartments, the viral DNA may be exposed to detection by Toll-like receptor 9 (TLR9). In this study, we characterized two new papillomaviruses from bats and used molecular archeology to demonstrate that their genomes altered their nucleotide compositions to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we showed that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for coevolution between papillomaviruses and their hosts. Finally, we demonstrated that the cancer-associated human papillomaviruses show a reduction in CpG dinucleotides within a TLR9 recognition complex.

Bias at the third nucleotide of codon pairs in virus and host genomes

Genomes of different sizes and complexity can be compared using common features. Most genomes contain open reading frames, and most genomes use the same genetic code. Redundancy in the genetic code means that different biases in the third nucleotide position of a codon exist in different genomes. However, the nucleotide composition of viruses can be quite different from host nucleotide composition making it difficult to assess the relevance of these biases. Here we show that grouping codons of a codon-pair according to the GC content of the first two nucleotide positions of each codon reveals patterns in nucleotide usage at the third position of the 1st codon. Differences between the observed and expected biases occur predominantly when the first two nucleotides of the 2nd codon are both S (strong, G or C) or both W (weak, A or T), not a mixture of strong and weak. The data indicates that some codon pairs are preferred because of the strength of the interactions between the codon and anticodon, the adjacent tRNAs and the ribosome. Using base-pairing strength and third position bias facilitates the comparison of genomes of different size and nucleotide composition and reveals patterns not previously described.

The slowing rate of CpG depletion in SARS-CoV-2 genomes is consistent with adaptations to the human host

Depletion of CpG dinucleotides in SARS-CoV-2 genomes has been linked to virus evolution, host-switching, virus replication, and innate immune responses. Temporal variations, if any, in the rate of CpG depletion during virus evolution in the host remain poorly understood. Here, we analysed the CpG content of over 1.4 million full-length SARS-CoV-2 genomes representing over 170 million documented infections during the first 17 months of the pandemic. Our findings suggest that the extent of CpG depletion in SARS-CoV-2 genomes is modest. Interestingly, the rate of CpG depletion is highest during early evolution in humans and it gradually tapers off almost reaching an equilibrium; this is consistent with adaptations to the human host. Furthermore, within the coding regions, CpG depletion occurs predominantly at codon positions 2-3 and 3-1. Loss of ZAP-binding motifs in SARS-CoV-2 genomes is primarily driven by the loss of the terminal CpG in the motifs. Nonetheless, majority of the CpG depletion in SARS-CoV-2 genomes occurs outside ZAP-binding motifs. SARS-CoV-2 genomes selectively lose CpGs-motifs from a U-rich context; this may help avoid immune recognition by TLR7. SARS-CoV-2 alpha-, beta- and delta-variants of concern have reduced CpG content compared to sequences from the beginning of the pandemic. In sum, we provide evidence that the rate of CpG depletion in virus genomes is not uniform and it greatly varies over time and during adaptations to the host. This work highlights how temporal variations in selection pressures during virus adaption may impact the rate and the extent of CpG depletion in virus genomes.

K-Mer Analyses Reveal Different Evolutionary Histories of Alpha, Beta, and Gamma Papillomaviruses

Papillomaviruses (PVs) are a heterogeneous group of DNA viruses that can infect fish, birds, reptiles, and mammals. PVs infecting humans (HPVs) phylogenetically cluster into five genera (Alpha-, Beta-, Gamma-, Mu- and Nu-PV), with differences in tissue tropism and carcinogenicity. The evolutionary features associated with the divergence of Papillomaviridae are not well understood. Using a combination of k-mer distributions, genetic metrics, and phylogenetic algorithms, we sought to evaluate the characteristics and differences of Alpha-, Beta- and Gamma-PVs constituting the majority of HPV genomes. A total of 640 PVs including 442 HPV types, 27 non-human primate PV types, and 171 non-primate animal PV types were evaluated. Our analyses revealed the highest genetic diversity amongst Gamma-PVs compared to the Alpha and Beta PVs, suggesting reduced selective pressures on Gamma-PVs. Using a sequence alignment-free trimer (k = 3) phylogeny algorithm, we reconstructed a phylogeny that grouped most HPV types into a monophyletic clade that was further split into three branches similar to alignment-based classifications. Interestingly, a subset of low-risk Alpha HPVs (the species Alpha-2, 3, 4, and 14) split from other HPVs and were clustered with non-human primate PVs. Surprisingly, the trimer-constructed phylogeny grouped the Gamma-6 species types originally isolated from the cervicovaginal region with the main Alpha-HPV clade. These data indicate that characterization of papillomavirus heterogeneity via orthogonal approaches reveals novel insights into the biological understanding of HPV genomes.

Limited expression of non-integrating CpG-free plasmid is associated with increased nucleosome enrichment

CpG-free pDNA was reported to facilitate sustained transgene expression with minimal inflammation in vivo as compared to CpG-containing pDNA. However, the expression potential and impact of CpG-free pDNA in in vitro model have never been described. Hence, in this study, we analyzed the transgene expression profiles of CpG-free pDNA in vitro to determine the influence of CpG depletion from the transgene. We found that in contrast to the published in vivo studies, CpG-free pDNA expressed a significantly lower level of luciferase than CpG-rich pDNA in several human cell lines. By comparing novel CpG-free pDNA carrying CpG-free GFP (pZGFP: 0 CpG) to CpG-rich GFP (pRGFP: 60 CpGs), we further showed that the discrepancy was not influenced by external factors such as gene transfer agent, cell species, cell type, and cytotoxicity. Moreover, pZGFP exhibited reduced expression despite having equal gene dosage as pRGFP. Analysis of mRNA distribution revealed that the mRNA export of pZGFP and pRGFP was similar; however, the steady state mRNA level of pZGFP was significantly lower. Upon further investigation, we found that the CpG-free transgene in non-integrating CpG-free pDNA backbone acquired increased nucleosome enrichment as compared with CpG-rich transgene, which may explain the observed reduced level of steady state mRNA. Our findings suggest that nucleosome enrichment could regulate non-integrating CpG-free pDNA expression and has implications on pDNA design.

Mutation Patterns of Human SARS-CoV-2 and Bat RaTG13 Coronavirus Genomes Are Strongly Biased Towards C>U Transitions, Indicating Rapid Evolution in Their Hosts

The pandemic caused by the spread of SARS-CoV-2 has led to considerable interest in its evolutionary origin and genome structure. Here, we analyzed mutation patterns in 34 human SARS-CoV-2 isolates and a closely related RaTG13 isolated from Rhinolophus affinis (a horseshoe bat). We also evaluated the CpG dinucleotide contents in SARS-CoV-2 and other human and animal coronavirus genomes. Out of 1136 single nucleotide variations (~4% divergence) between human SARS-CoV-2 and bat RaTG13, 682 (60%) can be attributed to C>U and U>C substitutions, far exceeding other types of substitutions. An accumulation of C>U mutations was also observed in SARS-CoV2 variants that arose within the human population. Globally, the C>U substitutions increased the frequency of codons for hydrophobic amino acids in SARS-CoV-2 peptides, while U>C substitutions decreased it. In contrast to most other coronaviruses, both SARS-CoV-2 and RaTG13 exhibited CpG depletion in their genomes. The data suggest that C-to-U conversion mediated by C deamination played a significant role in the evolution of the SARS-CoV-2 coronavirus. We hypothesize that the high frequency C>U transitions reflect virus adaptation processes in their hosts, and that SARS-CoV-2 could have been evolving for a relatively long period in humans following the transfer from animals before spreading worldwide.

Genetic and codon usage bias analyses of major capsid protein gene in Ranavirus

The Ranavirus (one genus of Iridovidae family) is an emerging pathogen that infects fish, amphibian, and reptiles, and causes great economical loss and ecological threat to farmed and wild animals globally. The major capsid protein (MCP) has been used as genetic typing marker and as target to design vaccines. Herein, the codon usage pattern of 73 MCP genes of Ranavirus and Lymphocystivirus are studied by calculating effective number of codons (ENC), relative synonymous codon usage (RSCU), codon adaptation index (CAI), and relative codon deoptimization index (RCDI), and similarity index (SiD). The Ranavirus are confirmed to be classified into five groups by using phylogenetic analysis, and varied nucleotide compositions and hierarchical cluster analysis based on RSCU. The results revealed different codon usage patterns among Lymphocystivirus and five groups of Ranavirus. Ranavirus had six over-represented codons ended with G/C nucleotide, while Lymphocystivirus had six over-represented codons ended with A/T nucleotide. A comparative analysis of parameters that define virus and host relatedness in terms of codon usage were analyzed indicated that Amphibian-like ranaviruses (ALRVs) seem to possess lower ENC values and higher CAIs in contrast to other ranaviruses isolated from fishes, and two groups (FV3-like and CMTV-like group) of them had received higher selection pressure from their hosts as having higher relative codon deoptimization index (RCDI) and similarity index (SiD). The correspondence analysis (COA) and Spearman's rank correlation analyses revealed that nucleotide compositions, relative dinucleotide frequency, mutation pressure, and natural translational selection shape the codon usage pattern in MCP genes and the ENC-GC_3S and neutrality plots indicated that the natural selection is the predominant factor. These results contribute to understanding the evolution of Ranavirus and their adaptions to their hosts.

Genomic characterization of the first oral avian papillomavirus in a colony of breeding canaries (Serinus canaria)

Papillomaviruses are non-enveloped, DNA viruses that infect skin and mucosa of a wide variety of vertebrates, causing neoplasias or simply persisting asymptomatically. Avian papillomaviruses, with six fully sequenced genomes, are the second most studied group after mammalian papillomaviruses. In this study, we describe the first oral avian papillomavirus, detected in the tongue of a dead Yorkshire canary (Serinus canaria) and in oral swabs of the same bird and other two live canaries from an aviary in Madrid, Spain. Its genome is 8,071 bp and presents the canonical papillomavirus architecture with six early (E6, E7, E1, E9, E2, E4) and two late open reading frames (L1 and L2) and a long control region between L1 and E6. This new avian papillomavirus L1 gene shares a 64% pairwise identity with FcPV1 L1, so it has been classified as a new species (ScPV1) within the Ethapapillomavirus genus. Although the canary died after showing breathing problems, there is no evidence that the papillomavirus caused those symptoms so it could be part of the oral microbiota of the birds. Hence, future investigations are needed to evaluate the clinical relevance of the virus.

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.

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.

The CpG dinucleotide content of the HIV-1 envelope gene may predict disease progression

The clinical course of HIV-1 varies greatly among infected individuals. Despite extensive research, virus factors associated with slow-progression remain poorly understood. Identification of unique HIV-1 genomic signatures linked to slow-progression remains elusive. We investigated CpG dinucleotide content in HIV-1 envelope gene as a potential virus factor in disease progression. We analysed 1808 HIV-1 envelope gene sequences from three independent longitudinal studies; this included 1280 sequences from twelve typical-progressors and 528 sequences from six slow-progressors. Relative abundance of CpG dinucleotides and relative synonymous codon usage (RSCU) for CpG-containing codons among HIV-1 envelope gene sequences from typical-progressors and slow-progressors were analysed. HIV-1 envelope gene sequences from slow-progressors have high-CpG dinucleotide content and increased number of CpG-containing codons as compared to typical-progressors. Our findings suggest that observed differences in CpG-content between typical-progressors and slow-progressors is not explained by differences in the mononucleotide content. Our results also highlight that the high-CpG content in HIV-1 envelope gene from slow-progressors is observed immediately after seroconversion. Thus CpG dinucleotide content of HIV-1 envelope gene is a potential virus-related factor that is linked to disease progression. The CpG dinucleotide content of HIV-1 envelope gene may help predict HIV-1 disease progression at early stages after seroconversion.

Evasion of host immune defenses by human papillomavirus

A majority of human papillomavirus (HPV) infections are asymptomatic and self-resolving in the absence of medical interventions. Various innate and adaptive immune responses, as well as physical barriers, have been implicated in controlling early HPV infections. However, if HPV overcomes these host immune defenses and establishes persistence in basal keratinocytes, it becomes very difficult for the host to eliminate the infection. The HPV oncoproteins E5, E6, and E7 are important in regulating host immune responses. These oncoproteins dysregulate gene expression, protein-protein interactions, posttranslational modifications, and cellular trafficking of critical host immune modulators. In addition to the HPV oncoproteins, sequence variation and dinucleotide depletion in papillomavirus genomes has been suggested as an alternative strategy for evasion of host immune defenses. Since anti-HPV host immune responses are also considered to be important for antitumor immunity, immune dysregulation by HPV during virus persistence may contribute to immune suppression essential for HPV-associated cancer progression. Here, we discuss cellular pathways dysregulated by HPV that allow the virus to evade various host immune defenses.