Human alpha and beta papillomaviruses use different synonymous codon profiles

Why HPV16? Why, now, HPV42? How the discovery of HPV42 in rare cancers provides an opportunity to challenge our understanding about the transition between health and disease for common members of the healthy microbiota

In 2022, a bioinformatic, agnostic approach identified HPV42 as causative agent of a rare cancer, later confirmed experimentally. This unexpected association offers an opportunity to reconsider our understanding about papillomavirus infections and cancers. We have expanded our knowledge about the diversity of papillomaviruses and the diseases they cause. Yet, we still lack answers to fundamental questions, such as what makes HPV16 different from the closely related HPV31 or HPV33; or why the very divergent HPV13 and HPV32 cause focal epithelial hyperplasia, while HPV6 or HPV42 do not, despite their evolutionary relatedness. Certain members of the healthy skin microbiota are associated to rare clinical conditions. We propose that a focus on cellular phenotypes, most often transient and influenced by intrinsic and extrinsic factors, may help understand the continuum between health and disease. A conceptual switch is required towards an interpretation of biology as a diversity of states connected by transition probabilities, rather than quasi-deterministic programs. Under this perspective, papillomaviruses may only trigger malignant transformation when specific viral genotypes interact with precise cellular states. Drawing on Canguilhem's concepts of normal and pathological, we suggest that understanding the transition between fluid cellular states can illuminate how commensal-like infections transition from benign to malignant.

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.

Mutation Profiles, Glycosylation Site Distribution and Codon Usage Bias of Human Papillomavirus Type 16

Human papillomavirus type 16 (HPV16) is the most prevalent HPV type causing cervical cancers. Herein, using 1597 full genomes, we systemically investigated the mutation profiles, surface protein glycosylation sites and the codon usage bias (CUB) of HPV16 from different lineages and sublineages. Multiple lineage- or sublineage-conserved mutation sites were identified. Glycosylation analysis showed that HPV16 lineage D contained the highest number of different glycosylation sites from lineage A in both L1 and L2 capsid proteins, which might lead to their antigenic distances between the two lineages. CUB analysis showed that the HPV16 open reading frames (ORFs) preferred codons ending with A/T. The CUB of HPV16 ORFs was mainly affected by natural selection except for E1, E5 and L2. HPV16 only shared some of the preferred codons with humans, which might help reduce competition in translational resources. These findings increase our understanding of the heterogeneity between HPV16 lineages and sublineages, and the adaptation mechanism of HPV in human cells. In summary, this study might facilitate HPV classification and improve vaccine development and application.

Analysis of Codon Usage Patterns in the Human Papillomavirus Oncogenes

Background:

Persistent high-risk genital human papillomavirus (HPV) infection is a major cause of cervical cancer in women. The products of the viral transforming genes E6 and E7 in the high-risk HPVs are known to be similar in their amino acid composition and structure. We performed a comparative analysis of codon usage patterns in the E6 and E7 genes of HPVs.

Methods:

The E6 and E7 gene sequences of eight HPV subtypes were analyzed to determine their nucleotide composition, relative synonymous codon usage (RSCU), effective number of codons (ENC), neutrality, genetic variability, selection pressure, and codon adaptation index (CAI). Additionally, a correspondence analysis (CoA) was performed.

Results:

The analysis to determine the effects of differences in composition on the codon usage patterns revealed that there may be usage bias for ‘A’ nucleotides. This was consistent with the results of the RSCU analysis, which demonstrated that the selection of A/T-rich patterns and the preference for A/T-ended codons in HPVs are influenced by compositional constraints. Moreover, the results reveal that selection pressure plays an important role in the CoA results for the RSCU values, Tajima’s D tests, and neutrality tests.

Conclusion:

The results of this study are consistent with previous findings that most papillomavirus genes are under purifying selection pressure, which limits changes to the encoded proteins. Natural selection and mutation pressures resulting in changes in the nucleotide composition and codon usage bias in the two tumor genes of HPV act differently during the evolution of the HPV subtype; thus, throughout the viral life cycle, HPV can constantly evolve to adapt to a new environment.

The Biological Impact of Genomic Diversity in Cervical Cancer Development

Human papillomaviruses (HPVs) are the etiologic agents of cervical cancer, the unique human neoplasia that has one single necessary cause. The diversity of HPVs is well described, with 200 HPV types existing as distinct taxonomic units and each receiving an Arabic number. On a clinical basis, they are usually grouped by their site of occurrence and disease associations. Those types inhabiting the anogenital mucosa are more intensively studied and further divided into cancer-associated HPVs, which are termed 'high risk', while those linked to benign proliferative lesions are assigned as 'low risk'. HPV16 is responsible for approximately 50% of all ICC cases, and paradoxically is one of the most prevalent types among healthy women. Longitudinal studies have shown that when an incidental HPV16 infection becomes persistent it will result in an enhanced risk for the development of high-grade lesions. However, it is unknown why some persistent, HPV16 infections (or infections by other HR-HPV types) progress to CIN3+ while most clear spontaneously. Several epidemiological investigations have focused on cofactors, from the most obvious such as cigarette and other carcinogenic exposures, to coinfections by other STDs such as chlamydia, with no significant findings. Thus, the current focus is on genomic variation from both virus and host. Such studies have been potentialized by the enormous technical advances in nucleic acid sequencing, allowing this relationship to be broadly interrogated. Corroborating subgenomic data from decades ago, an association between HPV16 lineages and carcinogenesis is being revealed. However, this effect does not seem to apply across female populations from different continents/ethnicities, again highlighting a role played by HPV16 adaptation and evasion from the host over time.

HPV disease transmission protection and control

Human papillomaviruses (HPVs) represent a large collection of viral types associated with significant clinical disease of cutaneous and mucosal epithelium. HPV-associated cancers are found in anogenital and oral mucosa, and at various cutaneous sites. Papillomaviruses are highly species and tissue restricted, and these viruses display both mucosotropic, cutaneotropic or dual tropism for epithelial tissues. A subset of HPV types, predominantly mucosal, are also oncogenic and cancers with these HPV types account for more than 200,000 deaths world-wide. Host control of HPV infections requires both innate and adaptive immunity, but the viruses have developed strategies to escape immune detection. Viral proteins can disrupt both innate pathogen-sensing pathways and T-cell based recognition and subsequent destruction of infected tissues. Current treatments to manage HPV infections include mostly ablative strategies in which recurrences are common and only active disease is treated. Although much is known about the papillomavirus life cycle, viral protein functions, and immune responsiveness, we still lack knowledge in a number of key areas of PV biology including tissue tropism, site-specific cancer progression, codon usage profiles, and what are the best strategies to mount an effective immune response to the carcinogenic stages of PV disease. In this review, disease transmission, protection and control are discussed together with questions related to areas in PV biology that will continue to provide productive opportunities of discovery and to further our understanding of this diverse set of human viral pathogens.

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.

Role of genomic and proteomic tools in the study of host-virus interactions and virus evolution

Viruses have short replication cycles and produce genomic variants within a host, a process that seems to adapt to their specific host and also enable them to infect new hosts. The recent emergence of viral genomic variants from the circulating pool within the host population and re-emergence of the old ones are posing serious threat to agriculture, animal husbandry and humanity as a whole. This review assesses the potential role of genomic and proteomic tools that can monitor not only the course of infection and pathogenesis, but also predict the pandemic or zoonotic epidemic potential of a virus in a previously exposed or immunologically naive biological population.

Synonymous codon changes in the oncogenes of the cottontail rabbit papillomavirus lead to increased oncogenicity and immunogenicity of the virus

Papillomaviruses use rare codons with respect to the host. The reasons for this are incompletely understood but among the hypotheses is the concept that rare codons result in low protein production and this allows the virus to escape immune surveillance. We changed rare codons in the oncogenes E6 and E7 of the cottontail rabbit papillomavirus to make them more mammalian-like and tested the mutant genomes in our in vivo animal model. While the amino acid sequences of the proteins remained unchanged, the oncogenic potential of some of the altered genomes increased dramatically. In addition, increased immunogenicity, as measured by spontaneous regression, was observed as the numbers of codon changes increased. This work suggests that codon usage may modify protein production in ways that influence disease outcome and that evaluation of synonymous codons should be included in the analysis of genetic variants of infectious agents and their association with disease.

Codon usage roles in human papillomavirus

Human papillomavirus (HPV) genomes, similar to other virus genomes, frequently have a G + C content significantly different from their host species. The HPV genomes show a strong codon usage bias to 18 codons, with 14 showing T at the third position amongst degenerately encoded amino acids. The codon usage pattern in HPV genome plays an important role, which regulates low or non-translational expression of the viral capsid genes and results in very weak protein expression of oncogenes in a wide range of mammalian cells. Codon modification has been proved to be a powerful technology to overcome the translational blockage and weak expression of both HPV capsid genes and oncogenes in different expression systems. Furthermore, keratinocytes are the host cells of HPV infection; the codon usage in HPV capsid genes matches available aminoacyl-tRNAs in differentiated keratinocytes to modulate their protein expression. HPV DNA vaccines with codon optimization have been shown to have higher immunogenicity and induce both strong cellular and humoral responses in animal models, which may be a promising form of therapeutic HPV vaccines.