mRNA Instability Elements in the Human Papillomavirus Type 16 L2 Coding Region

Human papillomavirus capsid proteins L1 and L2 are detected only in terminally differentiated cells, indicating that expression of the L1 and L2 genes is blocked in dividing cells. The results presented here establish that the human papillomavirus type 16 L2 coding region contains cis-acting inhibitory sequences. When placed downstream of a reporter gene, the human papillomavirus type 16 L2 sequence reduced both mRNA and protein levels in an orientation-dependent manner. Deletion analysis revealed that the L2 sequence contains two cis-acting inhibitory RNA regions. We identified an inhibitory region in the 5′-most 845 nucleotides of L2 that acted by reducing cytoplasmic mRNA stability and a second, weaker inhibitory region in the 3′ end of L2. In contrast, human papillomavirus type 1 L1 and L2 genes did not encode strong inhibitory sequences. This result is consistent with observations of high virus production in human papillomavirus type 1-infected tissue, whereas only low levels of human papillomavirus type 16 virions are detectable in infected epithelium. The presence of inhibitory sequences in the L1 and L2 mRNAs may aid the virus in avoiding the host immunosurveillance and in establishing persistent infections.     

Specific inactivation of inhibitory sequences in the 5' end of the human papillomavirus type 16 L1 open reading frame results in production of high levels of L1 protein in human epithelial cells

The expression of human papillomavirus type 16 late genes encoding virus capsid proteins L1 and L2 is restricted to terminally differentiated epithelial cells in the superficial layers of the squamous epithelium. We wish to understand the molecular mechanisms that determine the levels of expression of the human papillomavirus type 16 late genes. We have previously shown that the L1 coding region contains inhibitory sequences. Here we extend previous findings to show that the 5' end of the L1 gene contains strong inhibitory sequences but that the 3' end does not. We show that the first 514 nucleotides of the L1 coding region contain multiple inhibitory elements that act independently of one another and that the major inhibitory element is located within the first 129 nucleotides of the L1 gene. Introduction of point mutations in the inhibitory elements in the 5' end of the L1 gene which altered the RNA sequence without affecting the protein sequence specifically inactivated the inhibitory elements and resulted in production of high levels of human papillomavirus type 16 L1 mRNA and protein in human epithelial cells. Furthermore, we show that inhibitory sequences are present in the L1 coding regions of multiple human papillomavirus types, demonstrating that these elements are conserved among the human papillomaviruses, and suggest that they have an important function in the viral life cycle.     

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

Exogenous transfer RNAs (tRNAs) favor translation of bovine papillomavirus 1 wild-type (wt) L1 mRNA in in vitro translation systems (Zhou et al. 1999, J. Virol., 73, 4972–4982). We, therefore, investigated whether papillomavirus (PV) wt L1 protein expression could be enhanced in eukaryotic cells following exogenous tRNA supplementation. Both Chinese hamster ovary (CHO) and Cos1 cells, transfected with PV1 wt L1 genes, effectively transcribed the genes but did not translate them. However, L1 protein translation was demonstrated following co-transfection with the L1 gene and a gene expressing tRNA^Ser(CGA). Cell lines, stably transfected with a bovine papillomavirus 1 (BPV1) wt L1 expression construct, produced L1 protein after the transfection of the tRNA^Ser(CGA) gene, but not following the transfection with basal vectors, suggesting that tRNA^Ser(CGA) gene enhanced wt L1 translation as a result of endogenous tRNA alterations and phosphorylation of translation initiation factors elF4E and elF2α in the tRNA^Ser(CGA) transfected L1 cell lines. The tRNA^Ser(CGA) gene expression significantly reduced translation of L1 proteins expressed from codon-modified (HB) PV L1 genes utilizing mammalian preferred codons, but had variable effects on translation of green fluorescent proteins (GFPs) expressed from six serine GFP variants. The changes of tRNA pools appear to match the codon composition of PV wt and HB L1 genes and serine GFP variants to regulate translation of their mRNAs. These findings demonstrate for the first time in eukaryotic cells that translation of the target genes can be differentially influenced by the provision of a single tRNA expression construct.     

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.     

HPV-16 E2 contributes to induction of HPV-16 late gene expression by inhibiting early polyadenylation

We provide evidence that the human papillomavirus (HPV) E2 protein regulates HPV late gene expression. High levels of E2 caused a read-through at the early polyadenylation signal pAE into the late region of the HPV genome, thereby inducing expression of L1 and L2 mRNAs. This is a conserved property of E2 of both mucosal and cutaneous HPV types. Induction could be reversed by high levels of HPV-16 E1 protein, or by the polyadenylation factor CPSF30. HPV-16 E2 inhibited polyadenylation in vitro by preventing the assembly of the CPSF complex. Both the N-terminal and hinge domains of E2 were required for induction of HPV late gene expression in transfected cells as well as for inhibition of polyadenylation in vitro. Finally, overexpression of HPV-16 E2 induced late gene expression from a full-length genomic clone of HPV-16. We speculate that the accumulation of high levels of E2 during the viral life cycle, not only turns off the expression of the pro-mitotic viral E6 and E7 genes, but also induces the expression of the late HPV genes L1 and L2.     

Genes expressed during sexual differentiation of Chlamydomonas reinhardtii

Four genes specifically expressed during gametogenesis of Chlamydomonas reinhardtii have been cloned and their expression patterns analyzed. mRNAs encoded by these gamete-specific genes (gas) were absent or present only at very low levels in vegetative cells and mature zygotes. In young zygotes 2 h after gamete fusion, the mRNAs of three gas genes still persisted. The gas mRNAs accumulated during gametic differentiation. The temporal patterns of accumulation of individual mRNAs differed; some started to increase early during gametogenesis, others accumulated in the late phase. The accumulation of one of the late mRNAs (gas28) was stricly light-dependent. To illustrate the utility of the genes cloned in the analysis of sexual differentiation in Chlamydomonas reinhardtii we show that in a gametogenesis-defective mutant, the expression of late genes is prevented while that of early genes is normal.     

RNA Regulatory Elements in the Genomes of Simple Retroviruses

The RNA genomes of simple retroviruses encode three genes (gag, pol,andenv) which are required for replication. In addition, there are at least three well-definedcis-acting structures which regulate important aspects of the viral life cycle. The packaging signal at the 5′ end of the RNA tags the genomic RNA for specific encapsidation into assembling virus. Since viral Env proteins are translated from spliced mRNAs,cis-acting splicing signals ensure that the proper ratio of spliced and unspliced viral RNAs is present in the infected cell. Finally,cis-acting elements at the 3′ end of the genome promote the export of unspliced RNAs from the nucleus for translation and encapsidation.     

Heterogeneous Nuclear Ribonucleoprotein C Proteins Interact with the Human Papillomavirus Type 16 (HPV16) Early 3′-Untranslated Region and Alleviate Suppression of HPV16 Late L1 mRNA Splicing

In order to identify cellular factors that regulate human papillomavirus type 16 (HPV16) gene expression, cervical cancer cells permissive for HPV16 late gene expression were identified and characterized. These cells either contained a novel spliced variant of the L1 mRNAs that bypassed the suppressed HPV16 late, 5′-splice site SD3632; produced elevated levels of RNA-binding proteins SRSF1 (ASF/SF2), SRSF9 (SRp30c), and HuR that are known to regulate HPV16 late gene expression; or were shown by a gene expression array analysis to overexpress the RALYL RNA-binding protein of the heterogeneous nuclear ribonucleoprotein C (hnRNP C) family. Overexpression of RALYL or hnRNP C1 induced HPV16 late gene expression from HPV16 subgenomic plasmids and from episomal forms of the full-length HPV16 genome. This induction was dependent on the HPV16 early untranslated region. Binding of hnRNP C1 to the HPV16 early, untranslated region activated HPV16 late 5′-splice site SD3632 and resulted in production of HPV16 L1 mRNAs. Our results suggested that hnRNP C1 controls HPV16 late gene expression.     

Characterization of Mus musculus Papillomavirus 1 Infection In Situ Reveals an Unusual Pattern of Late Gene Expression and Capsid Protein Localization

Full-length genomic DNA of the recently identified laboratory mouse papillomavirus 1 (MusPV1) was synthesized in vitro and was used to establish and characterize a mouse model of papillomavirus pathobiology. MusPV1 DNA, whether naked or encapsidated by MusPV1 or human papillomavirus 16 (HPV 16) capsids, efficiently induced the outgrowth of papillomas as early as 3 weeks after application to abraded skin on the muzzles and tails of athymic NCr nude mice. High concentrations of virions were extracted from homogenized papillomatous tissues and were serially passaged for >10 generations. Neutralization by L1 antisera confirmed that infectious transmission was capsid mediated. Unexpectedly, the skin of the murine back was much less susceptible to virion-induced papillomas than the muzzle or tail. Although reporter pseudovirions readily transduced the skin of the back, infection with native MusPV1 resulted in less viral genome amplification and gene expression on the back, including reduced expression of the L1 protein and very low expression of the L2 protein, results that imply skin region-specific control of postentry aspects of the viral life cycle. Unexpectedly, L1 protein on the back was predominantly cytoplasmic, while on the tail the abundant L1 was cytoplasmic in the lower epithelial layers and nuclear in the upper layers. Nuclear localization of L1 occurred only in cells that coexpressed the minor capsid protein, L2. The pattern of L1 protein staining in the infected epithelium suggests that L1 expression occurs earlier in the MusPV1 life cycle than in the life cycle of high-risk HPV and that virion assembly is regulated by a previously undescribed mechanism.     

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer cells.     

Single-cell analysis of Kaposi’s sarcoma-associated herpesvirus infection in three-dimensional air-liquid interface culture model

The oral cavity is the major site for transmission of Kaposi’s sarcoma-associated herpesvirus (KSHV), but how KSHV establishes infection and replication in the oral epithelia remains unclear. Here, we report a KSHV spontaneous lytic replication model using fully differentiated, three-dimensional (3D) oral epithelial organoids at an air-liquid interface (ALI). This model revealed that KSHV infected the oral epithelia when the basal epithelial cells were exposed by damage. Unlike two-dimensional (2D) cell culture, 3D oral epithelial organoid ALI culture allowed high levels of spontaneous KSHV lytic replication, where lytically replicating cells were enriched at the superficial layer of epithelial organoid. Single cell RNA sequencing (scRNAseq) showed that KSHV infection induced drastic changes of host gene expression in infected as well as uninfected cells at the different epithelial layers, resulting in altered keratinocyte differentiation and cell death. Moreover, we identified a unique population of infected cells containing lytic gene expression at the KSHV K2-K5 gene locus and distinct host gene expression compared to latent or lytic infected cells. This study demonstrates an in vitro 3D epithelial organoid ALI culture model that recapitulates KSHV infection in the oral cavity, where KSHV undergoes the epithelial differentiation-dependent spontaneous lytic replication with a unique cell population carrying distinct viral gene expression.