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.

Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses

Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.

Papillomaviruses: Viral evolution, cancer and evolutionary medicine

Papillomaviruses (PVs) are a numerous family of small dsDNA viruses infecting virtually all mammals. PVs cause infections without triggering a strong immune response, and natural infection provides only limited protection against reinfection. Most PVs are part and parcel of the skin microbiota. In some cases, infections by certain PVs take diverse clinical presentations from highly productive self-limited warts to invasive cancers. We propose PVs as an excellent model system to study the evolutionary interactions between the immune system and pathogens causing chronic infections: genotypically, PVs are very diverse, with hundreds of different genotypes infecting skin and mucosa; phenotypically, they display extremely broad gradients and trade-offs between key phenotypic traits, namely productivity, immunogenicity, prevalence, oncogenicity and clinical presentation. Public health interventions have been launched to decrease the burden of PV-associated cancers, including massive vaccination against the most oncogenic human PVs, as well as systematic screening for PV chronic anogenital infections. Anti-PVs vaccines elicit protection against infection, induce cross-protection against closely related viruses and result in herd immunity. However, our knowledge on the ecological and intrapatient dynamics of PV infections remains fragmentary. We still need to understand how the novel anthropogenic selection pressures posed by vaccination and screening will affect viral circulation and epidemiology. We present here an overview of PV evolution and the connection between PV genotypes and the phenotypic, clinical manifestations of the diseases they cause. This differential link between viral evolution and the gradient cancer-warts-asymptomatic infections makes PVs a privileged playground for evolutionary medicine research.

Role of the host restriction factor APOBEC3 on papillomavirus evolution

More than 270 different types of papillomaviruses have been discovered in a wide array of animal species. Despite the great diversity of papillomaviruses, little is known about the evolutionary processes that drive host tropism and the emergence of oncogenic genotypes. Although host defense mechanisms have evolved to interfere with various aspects of a virus life cycle, viruses have also coevolved copious strategies to avoid host antiviral restriction. Our and other studies have shown that the cytidine deaminase APOBEC3 family members edit HPV genomes and restrict virus infectivity. Thus, we hypothesized that host restriction by APOBEC3 served as selective pressure during papillomavirus evolution. To test this hypothesis, we analyzed the relative abundance of all dinucleotide sequences in full-length genomes of 274 papillomavirus types documented in the Papillomavirus Episteme database (PaVE). Here, we report that TC dinucleotides, the preferred target sequence of several human APOBEC3 proteins (hA3A, hA3B, hA3F, and hA3H), are highly depleted in papillomavirus genomes. Given that HPV infection is highly tissue-specific, the expression levels of APOBEC3 family members were analyzed. The basal expression levels of all APOBEC3 isoforms, excluding hA3B, are significantly higher in mucosal skin compared with cutaneous skin. Interestingly, we reveal that Alphapapillomaviruses (alpha-PVs), a majority of which infects anogenital mucosa, display the most dramatic reduction in TC dinucleotide content. Computer modeling and reconstruction of ancestral alpha-PV genomes suggest that TC depletion occurred after the alpha-PVs diverged from their most recent common ancestor. In addition, we found that TC depletion in alpha-PVs is greatly affected by protein coding potential. Taken together, our results suggest that PVs replicating in tissues with high APOBEC3 levels may have evolved to evade restriction by selecting for variants that contain reduced APOBEC3 target sites in their genomes.

Evolution of the papillomaviridae

Viruses belonging to the Papillomaviridae family have been isolated from a variety of mammals, birds and non-avian reptiles. It is likely that most, if not all, amniotes carry a broad array of viral types. To date, the complete genomic sequence of more than 240 distinct viral types has been characterized at the nucleotide level. The analysis of this sequence information has begun to shed light on the evolutionary history of this important virus family. The available data suggests that many different evolutionary mechanisms have influenced the papillomavirus phylogenetic tree. Increasing evidence supports that the ancestral papillomavirus initially specialized to infect different ecological niches on the host. This episode of niche sorting was followed by extensive episodes of co-speciation with the host. This review attempts to summarize our current understanding of the papillomavirus evolution.

Human alpha and beta papillomaviruses use different synonymous codon profiles

Human papillomaviruses use rare codons relative to their hosts. It has been theorized that this is a mechanism to allow the virus to escape immune surveillance. In the present study, we examined the codings of four major genes of 21 human alpha (mucosatropic) viruses and 16 human beta (cutaneous-tropic) viruses. We compared the codon usage of different genes from a given papillomavirus and also the same genes from different papillomaviruses. Our data showed that codon usage was not always uniform between two genes of a given papillomavirus or between the same genes of papillomaviruses from different genera. We speculate as to why this might be and conclude that codon usage in the papillomaviruses may not only play a role in facilitating escape from immune surveillance but may also underlie some of the unanswered questions in the papillomavirus field.

RCDI/eRCDI: a web-server to estimate codon usage deoptimization

Background

The Relative Codon Deoptimization Index (RCDI) was developed by Mueller et al. (2006) as measure of codon deoptimization by comparing how similar is the codon usage of a gene and the codon usage of a reference genome.

Findings

RCDI/eRCDI is a web application server that calculates the Relative Codon Deoptimization Index and a new expected value for the RCDI (eRCDI). The RCDI is used to estimate the similarity of the codon frequencies of a specific gene in comparison to a given reference genome. The eRCDI is determined by generating random sequences with similar G+C and amino acid composition to the input sequences and may be used as an indicator of the significance of the RCDI values. RCDI/eRCDI is freely available at http://genomes.urv.cat/CAIcal/RCDI.

Conclusions

This web server will be a useful tool for genome analysis, to understand host-virus phylogenetic relationships or to infer the potential host range of a virus and its replication strategy, as well as in experimental virology to ease the step of gene design for heterologous protein expression.

Genomic characterization of the first insectivoran papillomavirus reveals an unusually long, second non-coding region and indicates a close relationship to Betapapillomavirus

Knowledge about biological diversity is the prerequisite to reliably reconstruct the evolution of pathogens such as papillomaviruses (PV). However, complete genomes of non-human PV have only been cloned and sequenced from 8 out of 18 orders within the Placentalia, although the host-specific variety of PV is considered much larger. We isolated and sequenced the complete genome of the first insectivoran PV type from hair follicle cells of the European hedgehog (Erinaceus europaeus), designated EHPV. We conducted phylogenetic analyses (maximum-likelihood criterion and Bayesian inference) with the genomic information of a systematically representative set of 67 PV types including EHPV. As inferred from amino acid sequence data of the separate genes E1, E2 and L1 as well as of the gene combination E6–E7–E1–E2–L1, EHPV clustered within the β-γ-π-Ξ-PV supertaxon and constituted the closest relative of genus Betapapillomavirus infecting primates. Beside the typical organization of the PV genome, EHPV exhibited a 1172 bp, non-coding region between the E2 and the L2 open reading frames. This trait has been previously described for the only distantly related Lambdapapillomavirus, but a common evolutionary origin of both non-coding regions is unlikely. Our results underscore the modular organization of the PV genome and the complex natural history of PV.

The first hydrophobic region of the HPV16 E5 protein determines protein cellular location and facilitates anchorage-independent growth

The human papillomavirus type 16 E5 protein (HPV16 E5) is 83 amino acids in length and contains three well-defined hydrophobic regions. The protein is expressed at very limited amounts in transfected cells and the absence of specific antibodies has strongly hampered functional analyses. To investigate the relationship between structure and function we have synthesized a codon-adapted version of the gene (hE5) and prepared a series of N-terminal and C-terminal deletions. Immunofluorescence analyses show colocaliation of the protein with calnexin, an ER marker, EEA-1, an early endosomes marker, and Lamp-2, a lysosomal marker. No major colocalization was found between hE5 and the Golgi marker 58 K. Whereas deletions at the C-terminal end of the protein do not greatly alter the localisation pattern, deletion of the first hydrophobic region results in loss of colocalisation with the ER, early endosomes and lysosomes. Further, we show that while the complete E5 protein confers to HaCaT cells the property to grow in an anchorage-independent manner, deletion of the first hydrophobic region results in loss of growth in soft agar. We conclude that the first hydrophobic region of the E5 protein largely determines the biological properties of the viral protein.

Expression of the HPV11 E2 gene in transgenic mice does not result in alterations of the phenotypic pattern

The E2 early protein of human papillomaviruses (HPV) has been found associated with the mitotic spindle therefore being implicated in the partition of the replicated viral DNA to daughter cells. In addition, E2 proteins bind to the upstream regulatory region of the virus and to cellular promoters modulating thereby cellular transcription and differentiation. In many cervical cancers, the E2 reading frame is interrupted upon incorporation of the viral genome into the host DNA. This results in the loss of the E2 mediated transcriptional repression and uncontrolled expression of the viral oncogenes. All these results have been obtained in transfected cells but no information is available on the E2 effects in the context of the entire organism. Transgenic mice were generated expressing the E2 protein of HPV11 under the control of the Ubiquitin C promoter. E2 mRNA is present in all mice tissues analysed and the E2 protein expressed in the skin (the target tissue of HPV11) was shown by Western blotting, albeit at a very low level. Analysis of the transgenic mice shows no major histological changes in the skin or all other tissues investigated. These data indicate that in transgenic mice the human papillomavirus type 11 E2 does not grossly modulate cellular proliferation or differentiation events.

Phylogeny and evolution of papillomaviruses based on the E1 and E2 proteins

Papillomaviridae are a family of small double-stranded DNA viruses that infect stratified squamous epithelia in vertebrates. Members of this family are causative agents of malignant tumours, such as cervical cancer while others are associated with benign proliferative lesions. So far, Papillomaviruses (PVs) are classified according to the sequence identity in the capsid gene L1. However, evidence has accumulated indicating a discontinuity in the evolutionary history of the L1 and L2 genes of many PVs, giving rise to differences in the phylogenetic reconstructions of the early and of the late genes. Neither the oncogenes E5, E6 and E7 nor the upstream regulatory region are suitable for phylogenetic inference due to the poor conservation along the Papillomaviridae family. We have analysed here the evolutionary relationships of the PVs with respect to the E1 and E2 proteins, and the results provide both phylogeny and biologic behaviour of the viruses. The hierarchical taxonomic relationships can be structured as an alternative classification system in which mucosal high-risk viruses, mucosal low-risk viruses and viruses associated with cutaneous lesions are grouped separately and do not appear intermingled. Some important trends are also observed: first, evolution of the PVs has not been homogeneous, even in viruses that infect the same host, and second mucosal human PVs have evolved faster than their cutaneous counterparts. The evolutionary analysis based on the E1 and E2 proteins will allow us to better understand the generation of the diversity of the PVs and the development of malignancy associated with these viruses.

Papillomaviruses: different genes have different histories

Papillomaviruses (PVs) infect stratified squamous epithelia in vertebrates. Some PVs are associated with different types of cancer and with certain benign lesions. It has been assumed that PVs coevolved with their hosts. However, recently it has been shown that different regions of the genome have different evolutionary histories. The PV genome has a modular nature and appeared after the addition of pre-existent blocks. This order of appearance in the PV genome is evident today in the different evolutionary rates of the different genes, with new genes--E5, E6 and E7--diverging faster than old genes--E1, E2, L2 and L1. Here, we propose an evolutionary framework aiming to integrate genome evolution, PV biology and epidemiology of PV infections.