tRNASer(CGA) differentially regulates expression of wild-type and codon-modified papillomavirus L1 genes

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.

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.

Codon usage variation of Zika virus: The potential roles of NS2B and NS4A in its global pandemic

A comprehensive demonstration of Zika virus (ZIKV) molecular evolution is essential for understanding its adaptation and expansion in its recent pandemics. Despite several studies on mutations and codon usage in ZIKVs, the variations in codon usage patterns across individual genes and their biological implication remains unclear. Here, we performed a gene-by-gene comparison of the codon usage variation in ZIKVs of the African and Asian lineages. We found that besides the evidence of positive selection (Ka/Ks >1) in the Asian lineage of the ZIKV genome, codon usage patterns were gene-specific and codon usage variation of ZIKV genes, was possibly constrained by their individual functional features, such as transmembrane domains, or antigenicity. In particular, the NS2B and NS4A genes showed distinct codon usage patterns, clearly separating them from the clusters of other genes in the correspondence analysis (CA). In the Asian lineage, the NS2B and NS4A genes showed the highest codon usage bias (ENC values: 51.01 ± 0.72 and 48.89 ± 0.99 respectively), and were subjected to the highest translation selection (ENC_obs/ENC_exp ratio: 0.847 ± 0.0297 and 0.828 ± 0.0233 respectively) in comparison to the African lineages of ZIKV. The CpG frequency of the NS2B showed a gradual ascending trend in the Asian ZIKV lineages, while in NS4A it was constrained along with the expansion of the Asian lineage. Furthermore, between the African and Asian lineages, differentiated and specific over-represented codons were more prominent in the NS2B and NS4A. Together, our study implies that ZIKVs are in the process of evolutionary fine tuning their codon as seen in the recent pandemics, and NS2B and NS4A could have played a potential role in the molecular evolution of the Asian lineage and their establishment.

Visualization and Modeling of the In Vivo Distribution of Mesenchymal Stem Cells

This unit describes a protocol for elucidating the in vivo disposition of administered mesenchymal stem cells (MSCs). Specifically, direct visualization of donor cell spatiotemporal distribution and assessment of donor cell quantity in recipient organs are described. Protocols for data analysis are suggested, with the goal of developing a model to characterize and predict the physiological kinetics of administered MSCs. The use of this model is described, suggesting that it can be applied to abnormal conditions and has potential interspecies and inter-route predictive capability. These universal methods can be employed, regardless of the type of stem cell or disease, to guide future experiments and design treatment protocols. © 2017 by John Wiley & Sons, Inc.

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.

Replication of human papillomavirus in culture

Human papillomaviruses (HPV) are the major factor in causing cervical cancer as well as being implicated in causing oral and anal cancers. The life cycle of HPV is tied to the epithelial differentiation system, as only native virus can be produced in stratified human skin. Initially, HPV research was only possible utilizing recombinant systems in monolayer culture. With new cell culture technology, systems using differentiated skin have allowed HPV to be studied in its native environment. Here, we describe current research studying native virions in differentiated skin including viral assembly, maturation, capsid protein interactions, and L2 cross-neutralizing epitopes. In doing so, we hope to show how differentiating skin systems have increased our knowledge of HPV biology and identify gaps in our knowledge about this important virus.

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.

Papillomavirus transcripts and posttranscriptional regulation

Papillomavirus gene expression is strictly linked to the differentiation state of the infected cell and is highly regulated at the level of transcription and RNA processing. All papillomaviruses make extensive use of alternative mRNA polyadenylation and splicing to control gene expression. This chapter contains a compilation of all known alternatively spliced papillomavirus mRNAs and it summarizes our current knowledge of viral RNA elements, and viral and cellular factors that control papillomavirus mRNA processing.

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.

Productive Lifecycle of Human Papillomaviruses that Depends Upon Squamous Epithelial Differentiation

Human papillomaviruses (HPVs) target the stratified epidermis, and can causes diseases ranging from benign condylomas to malignant tumors. Infections of HPVs in the genital tract are among the most common sexually transmitted diseases, and a major risk factor for cervical cancer. The virus targets epithelial cells in the basal layer of the epithelium, while progeny virions egress from terminally differentiated cells in the cornified layer, the surface layer of the epithelium. In infected basal cells, the virus maintains its genomic DNA at low-copy numbers, at which the viral productive lifecycle cannot proceed. Progression of the productive lifecycle requires differentiation of the host cell, indicating that there is tight crosstalk between viral replication and host differentiation programs. In this review, we discuss the regulation of the HPV lifecycle controlled by the differentiation program of the host cells.

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.

Synonymous codon usage influences the local protein structure observed

Translation of mRNA into protein is a unidirectional information flow process. Analysing the input (mRNA) and output (protein) of translation, we find that local protein structure information is encoded in the mRNA nucleotide sequence. The Coding Sequence and Structure (CSandS) database developed in this work provides a detailed mapping between over 4000 solved protein structures and their mRNA. CSandS facilitates a comprehensive analysis of codon usage over many organisms. In assigning translation speed, we find that relative codon usage is less informative than tRNA concentration. For all speed measures, no evidence was found that domain boundaries are enriched with slow codons. In fact, genes seemingly avoid slow codons around structurally defined domain boundaries. Translation speed, however, does decrease at the transition into secondary structure. Codons are identified that have structural preferences significantly different from the amino acid they encode. However, each organism has its own set of ‘significant codons’. Our results support the premise that codons encode more information than merely amino acids and give insight into the role of translation in protein folding.

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.

Expression of papillomavirus L1 proteins regulated by authentic gene codon usage is favoured in G2/M-like cells in differentiating keratinocytes

We investigated whether differentiation-dependent expression of papillomavirus (PV) L1 genes is influenced by the cell cycle state in keratinocytes (KCs) grown in vitro or in vivo. In primary keratinocytes, flow cytometry revealed a clear shift from predominantly G0/G1 to G2/M cells from day 1 to day 7, with a three-fold increase in G2/M-like cells in day 7 keratinocytes that showed approximately 50% of the cells expressed a terminal differentiation marker involucrin. The correlation between the levels of the L1 proteins expressed from authentic (Nat) L1 genes of HPV6b and BPV1 and the frequencies of the G2/M-like KCs was significantly positive, while in contrast, a significantly negative correlation in the levels of L1 proteins expressed from codon-modified (Mod) L1 genes of HPV6b and BPV1 with the frequencies of the G2/M-like KCs was observed. Experiments using cell cycle arrest reagents (all-trans retinoic acid (RA) and colchicine) confirmed that L1 proteins expressed from PV Nat L1 genes were facilitated in G2/M-like KCs upon differentiation. Using immunofluorescence microscopy, it appears that L1 proteins from PV Nat L1 genes were co-expressed with cyclin B1, while the L1 proteins expressed from PV Mod L1 genes were preferentially associated with cyclin D2 in KCs in vitro and in mouse skin. Our results demonstrate that (1) expression of the L1 proteins from Nat L1 genes of HPV6b and BPV1 that have strong codon usage bias with A or T at codon third position dependent on KC differentiation is favoured by the G2/M-like environment and (2) codon modifications can alter the cell differentiation-dependent and cell cycle-associated patterns of expression of the PV L1 proteins in KCs.

Replication and assembly of human papillomaviruses

Human papillomaviruses (HPVs) are small dsDNA tumor viruses, which are the etiologic agents of most cervical cancers and are associated with a growing percentage of oropharyngeal cancers. The HPV capsid is non-enveloped, having a T=7 icosahedral symmetry formed via the interaction among 72 pentamers of the major capsid protein, L1. The minor capsid protein L2 associates with L1 pentamers, although it is not known if each L1 pentamer contains a single L2 protein. The HPV life cycle strictly adheres to the host cell differentiation program, and as such, native HPV virions are only produced in vivo or in organotypic "raft" culture. Research producing synthetic papillomavirus particles--such as virus-like particles (VLPs), papillomavirus-based gene transfer vectors, known as pseudovirions (PsV), and papillomavirus genome-containing quasivirions (QV)--has bypassed the need for stratifying and differentiating host tissue in viral assembly and has allowed for the rapid analysis of HPV infectivity pathways, transmission, immunogenicity, and viral structure.

Adenovirus E4orf4 induces HPV-16 late L1 mRNA production

The adenovirus E4orf4 protein regulates the switch from early to late gene expression during the adenoviral replication cycle. Here we report that overexpression of adenovirus E4orf4 induces human papillomavirus type 16 (HPV-16) late gene expression from subgenomic expression plasmids. E4orf4 specifically overcomes the negative effects of two splicing silencers at the two late HPV-16 splice sites SD3632 and SA5639. This results in the production of HPV-16 spliced L1 mRNAs. We show that the interaction of E4orf4 with protein phosphatase 2A (PP2A) is necessary for induction of HPV-16 late gene expression. Also an E4orf4 mutant that fails to bind the cellular splicing factor ASF/SF2 fails to induce L1 mRNA production. Collectively, these results suggest that dephosphorylation of SR proteins by E4orf4 activates HPV-16 late gene expression. Indeed, a mutant ASF/SF2 protein in which the RS-domain had been deleted could itself induce HPV-16 late gene expression, whereas wild type ASF/SF2 could not.

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.

Generalized substitution of isoencoding codons shortens the duration of papillomavirus L1 protein expression in transiently gene-transfected keratinocytes due to cell differentiation

Recently we reported that gene codon composition determines differentiation-dependent expression of the PV L1 genes in mouse primary keratinocytes (KCs) in vitro and in vivo (Zhao et al. 2005, Mol. Cell Biol. 25:8643–8655). Here, we investigated whether generalized substitution of isoencoding codons affects the duration of expression of PV L1 genes in mouse and human KCs in day 1 culture transiently transfected with native (Nat) and codon modified (Mod) L1 genes. Following transient transfection, KC continuously transcribed both Nat and Mod PV L1 genes for at least 12 days, with the levels of L1 mRNAs from the Mod L1 genes significantly higher than those from the Nat L1 genes. However, continuous L1 protein expression at day 9 post-transfection was observed for both mouse and human KCs transfected with the Nat L1 genes only. Further, aa-tRNAs prepared from D8 KC cultures enhanced translation of two PV Nat L1 DNAs in RRL lysate and PV Nat L1 mRNAs in D0 cell-free lysate, whereas aa-tRNAs from D0 KCs enhanced translation of PV Mod L1 mRNAs in D8 cell-free lysate. It appears that aa-tRNAs in less-differentiated and differentiated KCs differentially match the PV Nat and Mod L1 mRNAs to regulate their translations in vitro.

Calcium enhances mouse keratinocyte differentiation in vitro to differentially regulate expression of papillomavirus authentic and codon modified L1 genes

Here, we first wished to establish for mouse primary keratinocytes (KCs) the Ca(2+) concentrations that were associated with KC differentiation in vitro. Using the range of Ca(2+) concentrations (0-6 mM) to differentiate primary KCs in culture to varying extents for 2 days, we then examined how KC differentiation impacted on expression of papillomavirus (PV) native (Nat) and codon modified (Mod) L1 genes. L1 mRNAs transcribed from either Nat or Mod L1 genes were present in similar amounts in KCs exposed to six Ca(2+) concentrations. However, expression of the L1 proteins from two Mod L1 genes were down-regulated, with no L1 signal detected in KCs exposed to 6 mM Ca(2+). In contrast, L1 proteins expressed from the two Nat L1 genes were not detectable in KCs without Ca(2+), but dramatically up-regulated as the KC cultures exposed to Ca(2+) from 0.5 to 2 mM, then down-regulated in KCs exposed to Ca(2+) from 4 to 6 mM. The different regulatory roles of the Ca(2+) in L1 protein expression from Nat and Mod L1 genes in cultured KCs were confirmed by TGF-beta1 experiments. We observed that aminoacyl-tRNAs (aa-tRNAs) from the 2 mM Ca(2+)-treated KCs only significantly enhanced the Nat L1 mRNAs translation in vitro, suggesting that aa-tRNAs play a differentially regulatory role in translations of the PV Nat and Mod L1 mRNAs. Importantly, the Ca(2+) experimental model provides evidence that mouse primary KCs could be transiently infected by BPV1 virus to express L1 mRNA and protein, which is very useful for future HPV virus infection study.

Papillomavirus genome structure, expression, and post-transcriptional regulation

Papillomaviruses are a group of small non-enveloped DNA tumor viruses whose infection usually causes benign epithelial lesions (warts). Certain types of HPVs, such as HPV-16, HPV-18, and HPV-31, have been recognized as causative agents of cervical cancer and anal cancer and their infections, which arise via sexual transmission, are associated with more than 95% of cervical cancer. Papillomaviruses infect keratinocytes in the basal layer of stratified squamous epithelia and replicate in the nucleus of infected keratinocytes in a differentiation-dependent manner. Viral gene expression in infected cells depends on cell differentiation and is tightly regulated at the transcriptional and post-transcriptional levels. A noteworthy feature of all papillomavirus transcripts is that they are transcribed as a bicistronic or polycistronic form containing two or more ORFs and are polyadenylated at either an early or late poly(A) site. In the past ten years, remarkable progress has been made in understanding how this complex viral gene expression is regulated at the level of transcription (such as via DNA methylation) and particularly post-transcription (including RNA splicing, polyadenylation, and translation). Current knowledge of papillomavirus mRNA structure and RNA processing has provided some clues on how to control viral oncogene expression. However, we still have little knowledge about which mRNAs are used to translate each viral protein. Continuing research on post-transcriptional regulation of papillomavirus infection will remain as a future focus to provide more insights into papillomavirus-host interactions, the virus life-cycle, and viral oncogenesis.

Gene codon composition determines differentiation-dependent expression of a viral capsid gene in keratinocytes in vitro and in vivo

By establishing mouse primary keratinocytes (KCs) in culture, we were able, for the first time, to express papillomavirus major capsid (L1) proteins by transient transfection of authentic or codon-modified L1 gene expression plasmids. We demonstrate in vitro and in vivo that gene codon composition is in part responsible for differentiation-dependent expression of L1 protein in KCs. L1 mRNA was present in similar amounts in differentiated and undifferentiated KCs transfected with authentic or codon-modified L1 genes and had a similar half-life, demonstrating that L1 protein production is posttranscriptionally regulated. We demonstrate further that KCs substantially change their tRNA profiles upon differentiation. Aminoacyl-tRNAs from differentiated KCs but not undifferentiated KCs enhanced the translation of authentic L1 mRNA, suggesting that differentiation-associated change to tRNA profiles enhances L1 expression in differentiated KCs. Thus, our data reveal a novel mechanism for regulation of gene expression utilized by a virus to direct viral capsid protein expression to the site of virion assembly in mature KCs. Analysis of two structural proteins of KCs, involucrin and keratin 14, suggests that translation of their mRNAs is also regulated, in association with KC differentiation in vitro, by a similar mechanism.