Codon usage bias and A+T content variation in human papillomavirus genomes

The human papillomavirus late life cycle and links to keratinocyte differentiation

Regulation of human papillomavirus (HPV) gene expression is tightly linked to differentiation of the keratinocytes the virus infects. HPV late gene expression is confined to the cells in the upper layers of the epithelium where the virus capsid proteins are synthesized. As these proteins are highly immunogenic, and the upper epithelium is an immune-privileged site, this spatial restriction aids immune evasion. Many decades of work have contributed to the current understanding of how this restriction occurs at a molecular level. This review will examine what is known about late gene expression in HPV-infected lesions and will dissect the intricacies of late gene regulation. Future directions for novel antiviral approaches will be highlighted.

Nucleotide and codon usage biases involved in the evolution of African swine fever virus: A comparative genomics analysis

Since African swine fever virus (ASFV) replication is closely related to its host's machinery, codon usage of viral genome can be subject to selection pressures. A better understanding of codon usage can give new insights into viral evolution. We implemented information entropy and revealed that the nucleotide usage pattern of ASFV is significantly associated with viral isolation factors (region and time), especially the usages of thymine and cytosine. Despite the domination of adenine and thymine in the viral genome, we found that mutation pressure alters the overall codon usage pattern of ASFV, followed by selective forces from natural selection. Moreover, the nucleotide skew index at the gene level indicates that nucleotide usages influencing synonymous codon bias of ASFV are significantly correlated with viral protein hydropathy. Finally, evolutionary plasticity is proved to contribute to the weakness in synonymous codons with A- or T-end serving as optimal codons of ASFV, suggesting that fine-tuning translation selection plays a role in synonymous codon usages of ASFV for adapting host. Taken together, ASFV is subject to evolutionary dynamics on nucleotide selections and synonymous codon usage, and our detailed analysis offers deeper insights into the genetic characteristics of this newly emerging virus around the world.

Modeling HPV-Associated Disease and Cancer Using the Cottontail Rabbit Papillomavirus

Approximately 5% of all human cancers are attributable to human papillomavirus (HPV) infections. HPV-associated diseases and cancers remain a substantial public health and economic burden worldwide despite the availability of prophylactic HPV vaccines. Current diagnosis and treatments for HPV-associated diseases and cancers are predominantly based on cell/tissue morphological examination and/or testing for the presence of high-risk HPV types. There is a lack of robust targets/markers to improve the accuracy of diagnosis and treatments. Several naturally occurring animal papillomavirus models have been established as surrogates to study HPV pathogenesis. Among them, the Cottontail rabbit papillomavirus (CRPV) model has become known as the gold standard. This model has played a pivotal role in the successful development of vaccines now available to prevent HPV infections. Over the past eighty years, the CRPV model has been widely applied to study HPV carcinogenesis. Taking advantage of a large panel of functional mutant CRPV genomes with distinct, reproducible, and predictable phenotypes, we have gained a deeper understanding of viral–host interaction during tumor progression. In recent years, the application of genome-wide RNA-seq analysis to the CRPV model has allowed us to learn and validate changes that parallel those reported in HPV-associated cancers. In addition, we have established a selection of gene-modified rabbit lines to facilitate mechanistic studies and the development of novel therapeutic strategies. In the current review, we summarize some significant findings that have advanced our understanding of HPV pathogenesis and highlight the implication of the development of novel gene-modified rabbits to future mechanistic studies.

Compositional biases in RNA viruses: Causes, consequences and applications

If each of the four nucleotides were represented equally in the genomes of viruses and the hosts they infect, each base would occur at a frequency of 25%. However, this is not observed in nature. Similarly, the order of nucleotides is not random (e.g., in the human genome, guanine follows cytosine at a frequency of ~0.0125, or a quarter the number of times predicted by random representation). Codon usage and codon order are also nonrandom. Furthermore, nucleotide and codon biases vary between species. Such biases have various drivers, including cellular proteins that recognize specific patterns in nucleic acids, that once triggered, induce mutations or invoke intrinsic or innate immune responses. In this review we examine the types of compositional biases identified in viral genomes and current understanding of the evolutionary mechanisms underpinning these trends. Finally, we consider the potential for large scale synonymous recoding strategies to engineer RNA virus vaccines, including those with pandemic potential, such as influenza A virus and Severe Acute Respiratory Syndrome Coronavirus Virus 2. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.

Edging on Mutational Bias, Induced Natural Selection From Host and Natural Reservoirs Predominates Codon Usage Evolution in Hantaan Virus

The molecular evolutionary dynamics that shape hantaviruses’ evolution are poorly understood even now, besides the contribution of virus-host interaction to their evolution remains an open question. Our study aimed to investigate these two aspects in Hantaan virus (HTNV)—the prototype of hantaviruses and an emerging zoonotic pathogen that infects humans, causing hemorrhagic fever with renal syndrome (HFRS): endemic in Far East Russia, China, and South Korea—via a comprehensive, phylogenetic-dependent codon usage analysis. We found that host- and natural reservoir-induced natural selection is the primary determinant of its biased codon choices, exceeding the mutational bias effect. The phylogenetic analysis of HTNV strains resulted in three distinct clades: South Korean, Russian, and Chinese. An effective number of codon (ENC) analysis showed a slightly biased codon usage in HTNV genomes. Nucleotide composition and RSCU analyses revealed a significant bias toward A/U nucleotides and A/U-ended codons, indicating the potential influence of mutational bias on the codon usage patterns of HTNV. Via ENC-plot, Parity Rule 2 (PR2), and neutrality plot analyses, we would conclude the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of HTNV; however, natural selection is the dominant factor influencing its codon usage bias. Codon adaptation index (CAI), Relative codon deoptimization index (RCDI), and Similarity Index (SiD) analyses uncovered the intense selection pressure from the host (Human) and natural reservoirs (Striped field mouse and Chinese white-bellied rat) in shaping HTNV biased codon choices. Our study clearly revealed the evolutionary processes in HTNV and the role of virus-host interaction in its evolution. Moreover, it opens the door for a more comprehensive codon usage analysis for all hantaviruses species to determine their molecular evolutionary dynamics and adaptability to several hosts and environments. We believe that our research will help in a better and deep understanding of HTNV evolution that will serve its future basic research and aid live attenuated vaccines design.

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.

HPV16 E7-impaired keratinocyte differentiation leads to tumorigenesis via cell cycle/pRb/involucrin/spectrin/adducin cascade

Here, we used codon usage technology to generate two codon-modified human papillomavirus (HPV)16 E7 genes and, together with wild-type E7, to construct three HPV16 E7 gene plasmids: Wt-E7, HB1-E7, and HB2-E7. The three HPV 16 E7 plasmids were used to investigate how HPV16 E7 protein was expressed in different cells and how this oncoprotein deregulated cellular and molecular events in human keratinocytes to induce carcinogenesis. We discovered that codon usage of HPV16 E7 gene played a key role in determining expression of E7 oncoprotein in all tested cells. HPV16 E7 inhibited significantly expression of pRb to impair keratinocyte differentiation and disrupted development of skin epidermis in mice. HPV16 E7 increased substantially the number of G0/G1 cells associated with upregulation of cyclin D2 and downregulation of cyclin B1 in keratinocytes. HPV16 E7 not only inhibited expression of involucrin and α-spectrin but also disrupted the organization of involucrin filaments and spectrin cytoskeleton. Furthermore, HPV16 E7 inhibited expression of β-adducin, destroyed its cytoskeletal structure and induced phosphorylation of β-adducin^(Ser662) in keratinocytes. Importantly, HPV16 E7 induced carcinogenesis in mice associated with expression of phosphorylated β-adducin^(Ser662) and its nucleus-translocation. In conclusion, we provided evidence that HPV16 E7 oncoprotein inhibited keratinocyte differentiation in vitro and in vivo leading to carcinogenesis through cell cycle arrest and disruption of pRb/involucrin/spectrin/adducin cascade.

The rabbit papillomavirus model: a valuable tool to study viral-host interactions

Cottontail rabbit papillomavirus (CRPV) was the first DNA virus shown to be tumorigenic. The virus has since been renamed and is officially known as Sylvilagus floridanus papillomavirus 1 (SfPV1). Since its inception as a surrogate preclinical model for high-risk human papillomavirus (HPV) infections, the SfPV1/rabbit model has been widely used to study viral-host interactions and has played a pivotal role in the successful development of three prophylactic virus-like particle vaccines. In this review, we will focus on the use of the model to gain a better understanding of viral pathogenesis, gene function and host immune responses to viral infections. We will discuss the application of the model in HPV-associated vaccine testing, in therapeutic vaccine development (using our novel HLA-A2.1 transgenic rabbits) and in the development and validation of novel anti-viral and anti-tumour compounds. Our goal is to demonstrate the role the SfPV1/rabbit model has played, and continues to play, in helping to unravel the intricacies of papillomavirus infections and to develop tools to thwart the disease. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

cRegions-a tool for detecting conserved cis-elements in multiple sequence alignment of diverged coding sequences

Identifying cis-acting elements and understanding regulatory mechanisms of a gene is crucial to fully understand the molecular biology of an organism. In general, it is difficult to identify previously uncharacterised cis-acting elements with an unknown consensus sequence. The task is especially problematic with viruses containing regions of limited or no similarity to other previously characterised sequences. Fortunately, the fast increase in the number of sequenced genomes allows us to detect some of these elusive cis-elements. In this work, we introduce a web-based tool called cRegions. It was developed to identify regions within a protein-coding sequence where the conservation in the amino acid sequence is caused by the conservation in the nucleotide sequence. The cRegion can be the first step in discovering novel cis-acting sequences from diverged protein-coding genes. The results can be used as a basis for future experimental analysis. We applied cRegions on the non-structural and structural polyproteins of alphaviruses as an example and successfully detected all known cis-acting elements. In this publication and in previous work, we have shown that cRegions is able to detect a wide variety of functional elements in DNA and RNA viruses. These functional elements include splice sites, stem-loops, overlapping reading frames, internal promoters, ribosome frameshifting signals and other embedded elements with yet unknown function. The cRegions web tool is available at http://bioinfo.ut.ee/cRegions/.

Evolution of Synonymous Codon Usage Bias in West African and Central African Strains of Monkeypox Virus

The evolution of bias in synonymous codon usage in chosen monkeypox viral genomes and the factors influencing its diversification have not been reported so far. In this study, various trends associated with synonymous codon usage in chosen monkeypox viral genomes were investigated, and the results are reported. Identification of factors that influence codon usage in chosen monkeypox viral genomes was done using various codon usage indices, such as the relative synonymous codon usage, the effective number of codons, and the codon adaptation index. The Spearman rank correlation analysis and a correspondence analysis were used for correlating various factors with codon usage. The results revealed that mutational pressure due to compositional constraints, gene expression level, and selection at the codon level for utilization of putative optimal codons are major factors influencing synonymous codon usage bias in monkeypox viral genomes. A cluster analysis of relative synonymous codon usage values revealed a grouping of more virulent strains as one major cluster (Central African strains) and a grouping of less virulent strains (West African strains) as another major cluster, indicating a relationship between virulence and synonymous codon usage bias. This study concluded that a balance between the mutational pressure acting at the base composition level and the selection pressure acting at the amino acid level frames synonymous codon usage bias in the chosen monkeypox viruses. The natural selection from the host does not seem to have influenced the synonymous codon usage bias in the analyzed monkeypox viral genomes.

The full transcription map of mouse papillomavirus type 1 (MmuPV1) in mouse wart tissues

Mouse papillomavirus type 1 (MmuPV1) provides, for the first time, the opportunity to study infection and pathogenesis of papillomaviruses in the context of laboratory mice. In this report, we define the transcriptome of MmuPV1 genome present in papillomas arising in experimentally infected mice using a combination of RNA-seq, PacBio Iso-seq, 5’ RACE, 3’ RACE, primer-walking RT-PCR, RNase protection, Northern blot and in situ hybridization analyses. We demonstrate that the MmuPV1 genome is transcribed unidirectionally from five major promoters (P) or transcription start sites (TSS) and polyadenylates its transcripts at two major polyadenylation (pA) sites. We designate the P₇₅₀₃, P₃₆₀ and P₈₅₉ as “early” promoters because they give rise to transcripts mostly utilizing the polyadenylation signal at nt 3844 and therefore can only encode early genes, and P₇₁₀₇ and P₅₃₃ as “late” promoters because they give rise to transcripts utilizing polyadenylation signals at either nt 3844 or nt 7047, the latter being able to encode late, capsid proteins. MmuPV1 genome contains five splice donor sites and three acceptor sites that produce thirty-six RNA isoforms deduced to express seven predicted early gene products (E6, E7, E1, E1^M1, E1^M2, E2 and E8^E2) and three predicted late gene products (E1^E4, L2 and L1). The majority of the viral early transcripts are spliced once from nt 757 to 3139, while viral late transcripts, which are predicted to encode L1, are spliced twice, first from nt 7243 to either nt 3139 (P₇₁₀₇) or nt 757 to 3139 (P₅₃₃) and second from nt 3431 to nt 5372. Thirteen of these viral transcripts were detectable by Northern blot analysis, with the P₅₃₃-derived late E1^E4 transcripts being the most abundant. The late transcripts could be detected in highly differentiated keratinocytes of MmuPV1-infected tissues as early as ten days after MmuPV1 inoculation and correlated with detection of L1 protein and viral DNA amplification. In mature warts, detection of L1 was also found in more poorly differentiated cells, as previously reported. Subclinical infections were also observed. The comprehensive transcription map of MmuPV1 generated in this study provides further evidence that MmuPV1 is similar to high-risk cutaneous beta human papillomaviruses. The knowledge revealed will facilitate the use of MmuPV1 as an animal virus model for understanding of human papillomavirus gene expression, pathogenesis and immunology.

Papillomavirus (PV) infections lead to development of both benign warts and cancers. Because PVs are epitheliotropic and species specific, it has been extremely challenging to study PV infection in the context of a naturally occurring infection in a tractable laboratory animal. The recent discovery of the papillomavirus, MmuPV1, that infects laboratory mice, provides an important new animal model system for understanding the pathogenesis of papillomavirus-associated diseases. By using state of the art RNA-seq to provide deep sequencing analysis of what regions of the viral genome are transcribed and PacBio Iso-seq that produces longer reads to define the complete sequences of individual transcripts in combination with several conventional technologies to confirm transcription starts sites, splice sites, and polyadenylation sites, we provide the first detailed description of the MmuPV1 transcript map using RNA from MmuPV1-induced mouse warts. This study reveals the presence of mRNA transcripts capable of coding for ten protein products in the MmuPV1 genome and leads to correctly re-assigning the E1^E4, L2 and L1 coding regions. We were able to detect individual transcripts from the infected wart tissues by RT-PCR, Northern blot and RNA ISH, to define the temporal onset of productive viral infection and to ectopically express a predicted viral protein for functional studies. The constructed MmuPV1 transcript map provides a foundation to advance our understanding of papillomavirus biology and pathogenesis.

Keratinocyte Differentiation-Dependent Human Papillomavirus Gene Regulation

Human papillomaviruses (HPVs) cause diseases ranging from benign warts to invasive cancers. HPVs infect epithelial cells and their replication cycle is tightly linked with the differentiation process of the infected keratinocyte. The normal replication cycle involves an early and a late phase. The early phase encompasses viral entry and initial genome replication, stimulation of cell division and inhibition of apoptosis in the infected cell. Late events in the HPV life cycle include viral genome amplification, virion formation, and release into the environment from the surface of the epithelium. The main proteins required at the late stage of infection for viral genome amplification include E1, E2, E4 and E5. The late proteins L1 and L2 are structural proteins that form the viral capsid. Regulation of these late events involves both cellular and viral proteins. The late viral mRNAs are expressed from a specific late promoter but final late mRNA levels in the infected cell are controlled by splicing, polyadenylation, nuclear export and RNA stability. Viral late protein expression is also controlled at the level of translation. This review will discuss current knowledge of how HPV late gene expression is regulated.

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.

Proteomic analysis of the gamma human papillomavirus type 197 E6 and E7 associated cellular proteins

Gamma HPV197 was the most frequently identified HPV when human skin cancer specimens were analyzed by deep sequencing (Arroyo Muhr et al., Int. J. Cancer 136: 2546-55, 2015). To gain insight into the biological activities of HPV197, we investigated the cellular interactomes of HPV197 E6 and E7. HPV197 E6 protein interacts with a broad spectrum of cellular LXXLL domain proteins, including UBE3A and MAML1. HPV197 E6 also binds and inhibits the TP53 tumor suppressor and interacts with the CCR4-NOT ubiquitin ligase and deadenylation complex. Despite lacking a canonical retinoblastoma (RB1) tumor suppressor binding site, HPV197 E7 binds RB1 and activates E2F transcription. Hence, HPV197 E6 and E7 proteins interact with a similar set of cellular proteins as E6 and E7 proteins encoded by HPVs that have been linked to human carcinogenesis and/or have transforming activities in vitro.

Cloning and Expression of L1 Protein Human Papillomavirus Type 31 Isolated from Iranian Patients in Escherichia coli

Human papillomavirus (HPV), a major pathogen of human cervical cancer, contains a full-length L1 gene encoding its surface capsid protein. One group of potential vaccine candidates against this virus in Iranian patients is based on surface protein components such as HPV31 L1 protein that can make virus-like particles (VLPs). The high immunity response stimulation of this effecter VLP was observed in host, suggesting that the individual characteristics of a particular effecter may require empirical testing for vaccination. In the present study, we decided to clone and express HPV31 L1 protein to investigate its use as a subunit vaccine and furthermore to insert the gene into an Escherichia coli background so as to analyze production of this recombinant protein. We report the presentation of HPV31 in 100 cervical lesion tissue samples based on polymerase chain reaction (PCR). Type of lesion, age, and other characteristics were reviewed and confirmed by a pathologist. The sequence from L1 genes of HPV was selected using special primers. The gene encoding the major capsid protein L1 was used for subcloning in pTG19-T and pET-32a plasmid. The recombinant protein expression was confirmed by RT-PCR using L1 primers and detected by absorption sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and immunoblot testing. The results presented here offer new insights into the in vivo response of HPV31 in Iranian patients and European models. On the other hand, the use of recombinant L1 protein for Iranian patient protection as well as vaccination studies will permit testing of this antigen protection rate and open the way to the discovery of protein biomarkers for monitoring clinical and subclinical cervical cancers.

Revelation of Influencing Factors in Overall Codon Usage Bias of Equine Influenza Viruses

Equine influenza viruses (EIVs) of H3N8 subtype are culprits of severe acute respiratory infections in horses, and are still responsible for significant outbreaks worldwide. Adaptability of influenza viruses to a particular host is significantly influenced by their codon usage preference, due to an absolute dependence on the host cellular machinery for their replication. In the present study, we analyzed genome-wide codon usage patterns in 92 EIV strains, including both H3N8 and H7N7 subtypes by computing several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis disclosed bias of preferred synonymous codons towards A/U-ended codons. The overall codon usage bias in EIVs was slightly lower, and mainly affected by the nucleotide compositional constraints as inferred from the RSCU and effective number of codon (ENc) analysis. Our data suggested that codon usage pattern in EIVs is governed by the interplay of mutation pressure, natural selection from its hosts and undefined factors. The H7N7 subtype was found less fit to its host (horse) in comparison to H3N8, by possessing higher codon bias, lower mutation pressure and much less adaptation to tRNA pool of equine cells. To the best of our knowledge, this is the first report describing the codon usage analysis of the complete genomes of EIVs. The outcome of our study is likely to enhance our understanding of factors involved in viral adaptation, evolution, and fitness towards their hosts.

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

Background

Papillomaviruses and polyomaviruses are small ds-DNA viruses infecting a wide-range of vertebrate hosts. Evidence supporting co-evolution of the virus with the host does not fully explain the evolutionary path of papillomaviruses and polyomaviruses. Studies analyzing CpG dinucleotide frequencies in virus genomes have provided interesting insights on virus evolution. CpG dinucleotide depletion has not been extensively studied among papillomaviruses and polyomaviruses. We sought to analyze the relative abundance of dinucleotides and the relative roles of evolutionary pressures in papillomaviruses and polyomaviruses.

Methods

We studied 127 full-length sequences from papillomaviruses and 56 full-length sequences from polyomaviruses. We analyzed the relative abundance of dinucleotides, effective codon number (ENC), differences in synonymous codon usage. We examined the association, if any, between the extent of CpG dinucleotide depletion and the evolutionary lineage of the infected host. We also investigated the contribution of mutational pressure and translational selection to the evolution of papillomaviruses and polyomaviruses.

Results

All papillomaviruses and polyomaviruses are CpG depleted. Interestingly, the evolutionary lineage of the infected host determines the extent of CpG depletion among papillomaviruses and polyomaviruses. CpG dinucleotide depletion was more pronounced among papillomaviruses and polyomaviruses infecting human and other mammals as compared to those infecting birds. Our findings demonstrate that CpG depletion among papillomaviruses is linked to mutational pressure; while CpG depletion among polyomaviruses is linked to translational selection. We also present evidence that suggests methylation of CpG dinucleotides may explain, at least in part, the depletion of CpG dinucleotides among papillomaviruses but not polyomaviruses.

Conclusions

The extent of CpG depletion among papillomaviruses and polyomaviruses is linked to the evolutionary lineage of the infected host. Our results highlight the existence of divergent evolutionary pressures leading to CpG dinucleotide depletion among small ds-DNA viruses infecting vertebrate hosts.

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.

Codon usage bias in human cytomegalovirus and its biological implication

Human cytomegalovirus (HCMV) infection, a worldwide contagion, causes a serious disorder in infected individuals. Analysis of codon usage can reveal much molecular information about this virus. The effective number of codon (ENC) values, relative synonymous codon usage (RSCU) values, codon adaptation index (CAI), and nucleotide contents was investigated in approximately 160 coding sequences (CDS) among 17 human cytomegalovirus genomes using the software CodonW. Linear regression analysis and logistic regression were performed to explore the preliminary data. The results showed that, overall, HCMV genomes had low codon usage bias (mean ENC=47.619). However, the ENC of individual CDS varied widely and was distributed unevenly between host-related genes and viral-self-function genes (P=0.002, odds ratio (OR)=3.194), as did the GC content (P=0.016, OR=2.178). The ENC values correlated with CAI, GC content, and the nucleotide composing at the 3rd codon position (GC3s) (P<0.001). There was a significant variation in the codon preference that depended on the RSCU data. The predicted ENC curve suggested that mutational pressure, rather than natural selection, was one of the main factors that determined the codon usage bias in HCMV. Among 123 genes with known function, the genes related to viral self-replication and viral-host interaction showed different ENC and CAI values, and GC and GC3s contents. In conclusion, the detailed codon usage bias theoretically revealed information concerning HCMV evolution and could be a valuable additional parameter for HCMV gene function research.

System analysis of synonymous codon usage biases in archaeal virus genomes

Recent studies of geothermally heated aquatic ecosystems have found widely divergent viruses with unusual morphotypes. Archaeal viruses isolated from these hot habitats usually have double-stranded DNA genomes, linear or circular, and can infect members of the Archaea domain. In this study, the synonymous codon usage bias (SCUB) and dinucleotide composition in the available complete archaeal virus genome sequences have been investigated. It was found that there is a significant variation in SCUB among different Archaeal virus species, which is mainly determined by the base composition. The outcome of correspondence analysis (COA) and Spearman׳s rank correlation analysis shows that codon usage of selected archaeal virus genes depends mainly on GC richness of genome, and the gene׳s function, albeit with smaller effects, also contributes to codon usage in this virus. Furthermore, this investigation reveals that aromaticity of each protein is also critical in affecting SCUB of these viral genes although it was less important than that of the mutational bias. Especially, mutational pressure may influence SCUB in SIRV1, SIRV2, ARV1, AFV1, and PhiCh1 viruses, whereas translational selection could play a leading role in HRPV1׳s SCUB. These conclusions not only can offer an insight into the codon usage biases of archaeal virus and subsequently the possible relationship between archaeal viruses and their host, but also may help in understanding the evolution of archaeal viruses and their gene classification, and more helpful to explore the origin of life and the evolution of biology.

•

The SCUB of archaeal virus genes depends mainly on GC richness of genome.

•

The mutational pressure is the main factor that influences SCUB.

•

The aromaticity of each protein is also critical in affecting SCUB.

•

The translational selection could play a leading role in HRPV1׳s SCUB.

•

The mode is helpful to explore the origin of life and the evolution of biology.

In vitro translation of papillomavirus authentic and codon-modified L1 capsid gene mRNAs in mouse keratinocyte cell-free lysate

Keratinocytes are the major cell type in the epidermis responsible for constructing the protective barrier of mammalian skin. Primary keratinocytes cultured in vitro can mimic the in vivo differentiation process, providing an abundant source of pure keratinocytes used for investigating the regulatory mechanisms of gene expression associated with cell proliferation and differentiation. We developed a primary mouse keratinocyte cell-free lysate system for translation of papillomavirus authentic and codon-modified L1 capsid gene mRNAs in vitro. We demonstrated that the viral gene codon usage matched available aminoacyl-tRNAs to determine their translational efficiencies, which were associated with differentiation status of the keratinocytes used for preparing cell-free lysate. We revealed that a novel regulatory mechanism of gene expression is utilized by papillomavirus to direct viral capsid protein expression to the site of virion assembly in mature keratinocytes. Here, I describe the methods in details of how to establish primary mouse keratinocyte culture, to prepare cell-free lysate, to carry out in vitro translation of the viral gene mRNAs, and to detect the translated products using Western blotting analysis.

An optimized yeast cell-free lysate system for in vitro translation of human virus mRNA

Yeast (Saccharomyces cerevisiae) as a model organism has long been established to study various aspects of eukaryotic cellular and molecular biology. Cell-free lysates prepared from different yeast strains have been used as a powerful tool to study eukaryotic protein expression in vitro. Recently, we established a yeast cell-free lysate system for in vitro translation long and short L1 capsid gene mRNAs of human papillomavirus type 58. We were able to significantly improve the translation efficiencies of the viral mRNAs in the established system by optimizing the concentrations of potassium and magnesium and controlling the physiological status of the yeast cells used for lysate preparation. We proved that a single specific amino acid can be rate limiting for translation of a target mRNA. Here, we describe the materials, method, and technique used for the development of an efficient yeast cell-free translation system.

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.

Nucleotide and phylogenetic analysis of human papillomavirus types 6 and 11 isolated from recurrent respiratory papillomatosis in Brazil

There are few studies about the distribution of natural molecular variants of low-risk HPVs. Our aim was to evaluate the E6 early gene variability among HPV-6 and HPV-11 isolates detected in recurrent respiratory papillomatosis (RRP) samples obtained in a cohort of Brazilian patients. We also performed a phylogenetic analysis in order to compare nucleotide sequences identified in our study with previously reported isolates from different anatomic sites (laryngeal papillomas, genital warts, cervical cancer and anal swabs) obtained from other parts of the world to determine the phylogenetic relationships of variants detected in Brazil. The complete coding region of the E6 gene of 25 samples was cloned and sequenced: 18 isolates of HPV-6 (72%) and 7 isolates of HPV-11 (28%). A total of four different HPV-6 genomic variants and two HPV-11 genomic variants was identified. It was not possible to correlate specific variants with disease severity. Phylogenetic trees for both HPV types were constructed enclosing both E6 sequences detected in our study and formerly published sequences. In both phylogenetic trees, the sequences from Brazil did not group together. We could not establish a geographical association between HPV-6 or HPV-11 variants, unlike HPV-16 and HPV-18.

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.

Analysis of synonymous codon usage patterns in torque teno sus virus 1 (TTSuV1)

Torque teno sus virus 1 (TTSuV1) is a novel virus that has been found widely distributed in the swine population in recent years. Analysis of codon usage can reveal much about the molecular evolution of TTSuV1. In this study, synonymous codon usage patterns and the key determinants in the coding region of 29 available complete TTSuV1 genome sequences were examined. By calculating the nucleotide content and relative synonymous codon usage (RSCU) of TTSuV1 coding sequences, we found that the preferentially used codons were mostly those ending with A or C nucleotides; less-used codons were mostly codons ending with U or G nucleotides, and these were mainly affected by composition constraints. Although there was a variation in codon usage bias among different TTSuV1 genomes, the codon usage bias and GC content in the TTSuV1 coding region was lower, which was mainly determined by the base composition in the third codon position and the effective number of codons (ENC) value. Moreover, the results of correspondence analysis (COA) indicated that the codon usage patterns of TTSuV1 isolated from different countries varied greatly and had significant differences. In addition, Spearman’s rank correlation analysis and an ENC plot revealed that apart from mutation pressure, which was critical in determining the codon usage pattern, other factors were involved in shaping the evolution of codon usage bias in TTSuV1, such as natural selection. Those results suggested that synonymous codon usage patterns of TTSuV1 genomes were the result of interaction between mutation pressure and natural selection. The information from this study not only provides important insights into the synonymous codon usage pattern of TTSuV1, but also helps to identify the main factors affecting codon usage by this virus.

Characterization and whole genome analysis of human papillomavirus type 16 e1-1374^63nt variants

BACKGROUND

The variation of the most common Human papillomavirus (HPV) type found in cervical cancer, the HPV16, has been extensively investigated in almost all viral genes. The E1 gene variation, however, has been rarely studied. The main objective of the present investigation was to analyze the variability of the E6 and E1 genes, focusing on the recently identified E1-1374^63nt variant.

METHODOLOGY/PRINCIPAL FINDINGS

Variation within the E6 of 786 HPV16 positive cervical samples was analyzed using high-resolution melting, while the E1-1374^63nt duplication was assayed by PCR. Both techniques were supplemented with sequencing. The E1-1374^63nt duplication was linked with the E-G350 and the E-C109/G350 variants. In comparison to the referent HPV16, the E1-1374^63nt E-G350 variant was significantly associated with lower grade cervical lesions (p = 0.029), while the E1-1374^63nt E-C109/G350 variant was equally distributed between high and low grade lesions. The E1-1374^63nt variants were phylogenetically closest to E-G350 variant lineage (A2 sub-lineage based on full genome classification). The major differences between E1-1374^63nt variants were within the LCR and the E6 region. On the other hand, changes within the E1 region were the major differences from the A2 sub-lineage, which has been historically but inconclusively associated with high grade cervical disease. Thus, the shared variations cannot explain the particular association of the E1-1374^63nt variant with lower grade cervical lesions.

CONCLUSIONS/SIGNIFICANCE

The E1 region has been thus far considered to be well conserved among all HPVs and therefore uninteresting for variability studies. However, this study shows that the variations within the E1 region could possibly affect cervical disease, since the E1-1374^63nt E-G350 variant is significantly associated with lower grade cervical lesions, in comparison to the A1 and A2 sub-lineage variants. Furthermore, it appears that the silent variation 109T>C of the E-C109/G350 variant might have a significant role in the viral life cycle and warrants further study.

Analysis of CpG methylation sites and CGI among human papillomavirus DNA genomes

Background

The Human Papillomavirus (HPV) genome is divided into early and late coding sequences, including 8 open reading frames (ORFs) and a regulatory region (LCR). Viral gene expression may be regulated through epigenetic mechanisms, including cytosine methylation at CpG dinucleotides. We have analyzed the distribution of CpG sites and CpG islands/clusters (CGI) among 92 different HPV genomes grouped in function of their preferential tropism: cutaneous or mucosal. We calculated the proportion of CpG sites (PCS) for each ORF and calculated the expected CpG values for each viral type.

Results

CpGs are underrepresented in viral genomes. We found a positive correlation between CpG observed and expected values, with mucosal high-risk (HR) virus types showing the smallest O/E ratios. The ranges of the PCS were similar for most genomic regions except E4, where the majority of CpGs are found within islands/clusters. At least one CGI belongs to each E2/E4 region. We found positive correlations between PCS for each viral ORF when compared with the others, except for the LCR against four ORFs and E6 against three other ORFs. The distribution of CpG islands/clusters among HPV groups is heterogeneous and mucosal HR-HPV types exhibit both lower number and shorter island sizes compared to cutaneous and mucosal Low-risk (LR) HPVs (all of them significantly different).

Conclusions

There is a difference between viral and cellular CpG underrepresentation. There are significant correlations between complete genome PCS and a lack of correlations between several genomic region pairs, especially those involving LCR and E6. L2 and L1 ORF behavior is opposite to that of oncogenes E6 and E7. The first pair possesses relatively low numbers of CpG sites clustered in CGIs while the oncogenes possess a relatively high number of CpG sites not associated to CGIs. In all HPVs, E2/E4 is the only region with at least one CGI and shows a higher content of CpG sites in every HPV type with an identified E4. The mucosal HR-HPVs show either the shortest CGI size, followed by the mucosal LR-HPVs and lastly by the cutaneous viral subgroup, and a trend to the lowest CGI number, followed by the cutaneous viral subgroup and lastly by the mucosal LR-HPVs.

Bioinformatics analysis of rabbit haemorrhagic disease virus genome

Background

Rabbit haemorrhagic disease virus (RHDV), as the pathogeny of Rabbit haemorrhagic disease, can cause a highly infectious and often fatal disease only affecting wild and domestic rabbits. Recent researches revealed that it, as one number of the Caliciviridae, has some specialties in its genome, its reproduction and so on.

Results

In this report, we firstly analyzed its genome and two open reading frameworks (ORFs) from this aspect of codon usage bias. Our researches indicated that mutation pressure rather than natural is the most important determinant in RHDV with high codon bias, and the codon usage bias is nearly contrary between ORF1 and ORF2, which is maybe one of factors regulating the expression of VP60 (encoding by ORF1) and VP10 (encoding by ORF2). Furthermore, negative selective constraints on the RHDV whole genome implied that VP10 played an important role in RHDV lifecycle.

Conclusions

We conjectured that VP10 might be beneficial for the replication, release or both of virus by inducing infected cell apoptosis initiate by RHDV. According to the results of the principal component analysis for ORF2 of RSCU, we firstly separated 30 RHDV into two genotypes, and the ENC values indicated ORF1 and ORF2 were independent among the evolution of RHDV.

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.

Human papillomavirus: gene expression, regulation and prospects for novel diagnostic methods and antiviral therapies

Human papillomaviruses (HPVs) cause diseases ranging from benign warts to invasive tumors. A subset of these viruses termed 'high risk' infect the cervix where persistent infection can lead to cervical cancer. Although many HPV genomes have been sequenced, knowledge of virus gene expression and its regulation is still incomplete. This is due in part to the lack, until recently, of suitable systems for virus propagation in the laboratory. HPV gene expression is polycistronic initiating from multiple promoters. Gene regulation occurs at transcriptional, but particularly post-transcriptional levels, including RNA processing, nuclear export, mRNA stability and translation. A close association between the virus replication cycle and epithelial differentiation adds a further layer of complexity. Understanding HPV mRNA expression and its regulation in the different diseases associated with infection may lead to development of novel diagnostic approaches and will reveal key viral and cellular targets for development of novel antiviral therapies.

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.

Viral adaptation to host: a proteome-based analysis of codon usage and amino acid preferences

Viruses differ markedly in their specificity toward host organisms. Here, we test the level of general sequence adaptation that viruses display toward their hosts. We compiled a representative data set of viruses that infect hosts ranging from bacteria to humans. We consider their respective amino acid and codon usages and compare them among the viruses and their hosts. We show that bacteria-infecting viruses are strongly adapted to their specific hosts, but that they differ from other unrelated bacterial hosts. Viruses that infect humans, but not those that infect other mammals or aves, show a strong resemblance to most mammalian and avian hosts, in terms of both amino acid and codon preferences. In groups of viruses that infect humans or other mammals, the highest observed level of adaptation of viral proteins to host codon usages is for those proteins that appear abundantly in the virion. In contrast, proteins that are known to participate in host-specific recognition do not necessarily adapt to their respective hosts. The implication for the potential of viral infectivity is discussed.

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.

Analysis of synonymous codon usage in classical swine fever virus

Using the complete genome sequences of 35 classical swine fever viruses (CSFV) representing all three genotypes and all three kinds of virulence, we analyzed synonymous codon usage and the relative dinucleotide abundance in CSFV. The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in CSFV. Furthermore, we observed that the relative abundance of dinucleotides in CSFV is independent of the overall base composition but is still the result of differential mutational pressure, which also shapes codon usage. In addition, other factors, such as the subgenotypes and aromaticity, also influence the codon usage variation among the genomes of CSFV. This study represents the most comprehensive analysis to date of CSFV codon usage patterns and provides a basic understanding of the mechanisms for codon usage bias.

CAIcal: a combined set of tools to assess codon usage adaptation

BACKGROUND

The Codon Adaptation Index (CAI) was first developed to measure the synonymous codon usage bias for a DNA or RNA sequence. The CAI quantifies the similarity between the synonymous codon usage of a gene and the synonymous codon frequency of a reference set.

RESULTS

We describe here CAIcal, a web-server available at http://genomes.urv.es/CAIcal that includes a complete set of utilities related with the CAI. The server provides useful important features, such as the calculation and graphical representation of the CAI along either an individual sequence or a protein multiple sequence alignment translated to DNA. The automated calculation of CAI and its expected value is also included as one of the CAIcal tools. The software is also free to be downloaded as a stand alone application for local use.

CONCLUSION

The CAIcal server provides a complete set of tools to assess codon usage adaptation and to help in genome annotation.

CRPV genomes with synonymous codon optimizations in the CRPV E7 gene show phenotypic differences in growth and altered immunity upon E7 vaccination

Papillomaviruses use rare codons relative to their hosts. Recent studies have demonstrated that synonymous codon changes in viral genes can lead to increased protein production when the codons are matched to those of cells in which the protein is being expressed. We theorized that the immunogenicity of the virus would be enhanced by matching codons of selected viral genes to those of the host. We report here that synonymous codon changes in the E7 oncogene are tolerated in the context of the cottontail rabbit papillomavirus (CRPV) genome. Papilloma growth rates differ depending upon the changes made indicating that synonymous codons are not necessarily neutral. Immunization with wild type E7 DNA yielded significant protection from subsequent challenge by both wild type and codon-modified genomes. The reduction in growth was most dramatic with the genome containing the greatest number of synonymous codon changes.

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.

Analysis of synonymous codon usage in 11 human bocavirus isolates

Human Bocavirus (HBoV) is a novel virus which can cause respiratory tract disease in infants or children. In this study, the codon usage bias and the base composition variations in the available 11 complete HBoV genome sequences have been investigated. Although, there is a significant variation in codon usage bias among different HBoV genes, codon usage bias in HBoV is a little slight, which is mainly determined by the base compositions on the third codon position and the effective number of codons (ENC) value. The results of correspondence analysis (COA) and Spearman's rank correlation analysis reveals that the G + C compositional constraint is the main factor that determines the codon usage bias in HBoV and the gene's function also contributes to the codon usage in this virus. Moreover, it was found that the hydrophobicity of each protein and the gene length are also critical in affecting these viruses’ codon usage, although they were less important than that of the mutational bias and the genes’ function. At last, the relative synonymous codon usage (RSCU) of 44 genes from these 11 HBoV isolates is analyzed using a hierarchical cluster method. The result suggests that genes with same function yet from different isolates are classified into the same lineage and it does not depend on geographical location. These conclusions not only can offer an insight into the codon usage patterns and gene classification of HBoV, but also may help in increasing the efficiency of gene delivery/expression systems.