Tissue-Specific Gene Expression during Productive Human Papillomavirus 16 Infection of Cervical, Foreskin, and Tonsil Epithelium

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.

Human Papillomavirus 16 Lineage A Variants Associated With Persistent Genital Infections in Men: The HPV Infection in Men (HIM) Study

BACKGROUND

Human papillomavirus (HPV) 16 non-A lineage variants have higher carcinogenic potential for cervical cancer. HPV-16 variants natural history among males is not established. We evaluated HPV-16 variants prevalence and persistence in the external genitalia of men enrolled in the prospective HPV Infection in Men (HIM) Study.

METHODS

The HIM Study included men from the United States, Brazil, and Mexico. HPV-16 variants were distinguished using polymerase chain reaction sequencing. The prevalence of HPV-16 variants was assessed, and associations with infection persistence were estimated.

RESULTS

We characterized the HPV-16 variants for 1700 genital swab samples from 753 men and 22 external genital lesions in 17 men. The prevalence of HPV-16 lineages differed by country and marital status (P < .001). Overall, 90.9% of participants harbored lineage A variants. The prevalence of non-A lineages was heterogenous among countries. HPV-16 lineage A variants were associated with a 2.69-fold increased risk of long-term persistent infections compared with non-A lineages. All high-grade penile intraepithelial neoplasia harbored lineage A variants and occurred in the context of long-term persistent infections with the same variants.

CONCLUSIONS

The prevalence and persistence of HPV-16 variants observed at the male external genitalia suggest differences in the natural history of these variants between men and women, which may be associated with intrinsic differences in the infected genital epithelia.

Characterization of 3D organotypic epithelial tissues reveals tonsil-specific differences in tonic interferon signaling

Three-dimensional (3D) culturing techniques can recapitulate the stratified nature of multicellular epithelial tissues. Organotypic 3D epithelial tissue culture methods have several applications, including the study of tissue development and function, drug discovery and toxicity testing, host-pathogen interactions, and the development of tissue-engineered constructs for use in regenerative medicine. We grew 3D organotypic epithelial tissues from foreskin, cervix, and tonsil-derived primary cells and characterized the transcriptome of these in vitro tissue equivalents. Using the same 3D culturing method, all three tissues yielded stratified squamous epithelium, validated histologically using basal and superficial epithelial cell markers. The goal of this study was to use RNA-seq to compare gene expression patterns in these three types of epithelial tissues to gain a better understanding of the molecular mechanisms underlying their function and identify potential therapeutic targets for various diseases. Functional profiling by over-representation and gene set enrichment analysis revealed tissue-specific differences: i.e., cutaneous homeostasis and lipid metabolism in foreskin, extracellular matrix remodeling in cervix, and baseline innate immune differences in tonsil. Specifically, tonsillar epithelia may play an active role in shaping the immune microenvironment of the tonsil balancing inflammation and immune responses in the face of constant exposure to microbial insults. Overall, these data serve as a resource, with gene sets made available for the research community to explore, and as a foundation for understanding the epithelial heterogeneity and how it may impact their in vitro use. An online resource is available to investigate these data (https://viz.datascience.arizona.edu/3DEpiEx/).

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

Deregulation of host gene expression by HPV16 E8^E2 knock-out genomes is due to increased productive replication

Productive replication of human papillomaviruses (HPV) only takes place in differentiating keratinocytes. The HPV16 E8^E2 protein acts as a repressor of viral gene expression and genome replication and HPV16 E8^E2 knock-out (E8-) genomes display enhanced viral late protein expression in differentiated cells. Global transcriptome analysis of differentiated HPV16 wild-type and E8-cell lines revealed a small number of differentially expressed genes which are not related to cell cycle, DNA metabolism or keratinocyte differentiation. The analysis of selected genes suggested that deregulation requires cell differentiation and positively correlated with the expression of viral late, not early transcripts. Consistent with this, the additional knock-out of the viral E4 and E5 genes, which are known to enhance productive replication, attenuated the deregulation of these host cell genes. In summary, these data reveal that productive HPV16 replication modulates host cell transcription.

Characterization of 3D organotypic epithelial tissues reveals tonsil-specific differences in tonic interferon signaling

Three-dimensional (3D) culturing techniques can recapitulate the stratified nature of multicellular epithelial tissues. Organotypic 3D epithelial tissue culture methods have several applications, including the study of tissue development and function, drug discovery and toxicity testing, host-pathogen interactions, and the development of tissue-engineered constructs for use in regenerative medicine. We grew 3D organotypic epithelial tissues from foreskin, cervix, and tonsil-derived primary cells and characterized the transcriptome of these in vitro tissue equivalents. Using the same 3D culturing method, all three tissues yielded stratified squamous epithelium, validated histologically using basal and superficial epithelial cell markers. The goal of this study was to use RNA-seq to compare gene expression patterns in these three types of epithelial tissues to gain a better understanding of the molecular mechanisms underlying their function and identify potential therapeutic targets for various diseases. Functional profiling by over-representation and gene set enrichment analysis revealed tissue-specific differences: i.e. , cutaneous homeostasis and lipid metabolism in foreskin, extracellular matrix remodeling in cervix, and baseline innate immune differences in tonsil. Specifically, tonsillar epithelia may play an active role in shaping the immune microenvironment of the tonsil balancing inflammation and immune responses in the face of constant exposure to microbial insults. Overall, these data serve as a resource, with gene sets made available for the research community to explore, and as a foundation for understanding the epithelial heterogeneity and how it may impact their in vitro use. An online resource is available to investigate these data ( https://viz.datascience.arizona.edu/3DEpiEx/ ).

The HPV Induced Cancer Resource (THInCR): a Suite of Tools for Investigating HPV-Dependent Human Carcinogenesis

Human papillomaviruses (HPVs) are highly infectious and cause the most common sexually transmitted viral infections. They induce hyperproliferation of squamous epithelial tissue, often forming warts. Virally encoded proteins reprogram gene expression and cell growth to create an optimal environment for viral replication. In addition to their normal roles in infection, functional alterations induced by viral proteins establish conditions that frequently contribute to human carcinogenesis. In fact, ~5% of human cancers are caused by HPVs, with virtually all cervical squamous cell carcinomas (CESC) and an increasing number of head and neck squamous cell carcinomas (HNSC) attributed to HPV infection. The Cancer Genome Atlas (TCGA) molecularly characterized thousands of primary human cancer samples in many cancer types, including CESC and HNSC, and created a comprehensive atlas of genomic, epigenomic, and transcriptomic data. This publicly available genome-wide information provides an unprecedented opportunity to expand the knowledge of the role that HPV plays in human carcinogenesis. While many tools exist to mine these data, few, if any, focus on the comparison of HPV-positive cancers with their HPV-negative counterparts or adjacent normal control tissue. We have constructed a suite of web-based tools, The HPV Induced Cancer Resource (THInCR), to utilize TCGA data for research related to HPV-induced CESC and HNSC. These tools allow investigators to gain greater biological and medical insights by exploring the impacts of HPV on cellular gene expression (mRNA and microRNA), altered gene methylation, and associations with patient survival and immune landscape features. These tools are accessible at https://thincr.ca/. IMPORTANCE The suite of analytical tools of THInCR provides the opportunity to investigate the roles that candidate target genes identified in cell lines or other model systems contribute to in actual HPV-dependent human cancers and is based on large-scale TCGA data sets. Expression of target genes, including both mRNA and microRNA, can be correlated with HPV gene expression, epigenetic changes in DNA methylation, patient survival, and numerous immune features, like leukocyte infiltration, interferon gamma response, T cell response, etc. Data from these analyses may immediately provide evidence to validate in vitro observations, reveal insights into mechanisms of virus-mediated alterations in cell growth, behavior, gene expression, and innate and adaptive immunity and may help hypothesis generation for further investigations.

Generation of human tonsil epithelial organoids as an ex vivo model for SARS-CoV-2 infection

The palatine tonsils (hereinafter referred to as “tonsils”) serve as a reservoir for viral infections and play roles in the immune system's first line of defense. The aims of this study were to establish tonsil epithelial cell–derived organoids and examine their feasibility as an ex vivo model for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. The tonsil organoids successfully recapitulated the key characteristics of the tonsil epithelium, including cellular composition, histologic properties, and biomarker distribution. Notably, the basal layer cells of the organoids express molecules essential for SARS-CoV-2 entry, such as angiotensin-converting enzyme 2 (ACE2), transmembrane serine protease 2 (TMPRSS2) and furin, being susceptible to the viral infection. Changes in the gene expression profile in tonsil organoids revealed that 395 genes associated with oncostatin M signaling and lipid metabolism were highly upregulated within 72 h after SARS-CoV-2 infection. Notably, remdesivir suppressed the viral RNA copy number in organoid culture supernatants and intracellular viral protein levels in a dose-dependent manner. Here, we suggest that tonsil epithelial organoids could provide a preclinical and translational research platform for investigating SARS-CoV-2 infectivity and transmissibility or for evaluating antiviral candidates.

The Key Differences between Human Papillomavirus-Positive and -Negative Head and Neck Cancers: Biological and Clinical Implications

Head and neck squamous cell carcinoma (HNSCC) is a unique malignancy associated with two distinct risk factors: exposure to typical carcinogens and infection of human papillomavirus (HPV). HPV encodes the potent oncoproteins E6 and E7, which bypass many important oncogenic processes and result in cancer development. In contrast, HPV-negative HNSCC is developed through multiple mutations in diverse oncogenic driver genes. While the risk factors associated with HPV-positive and HPV-negative HNSCCs are discrete, HNSCC patients still show highly complex molecular signatures, immune infiltrations, and treatment responses even within the same anatomical subtypes. Here, we summarize the current understanding of biological mechanisms, treatment approaches, and clinical outcomes in comparison between HPV-positive and -negative HNSCCs.

Coronavirus Porcine Epidemic Diarrhea Virus Nucleocapsid Protein Interacts with p53 To Induce Cell Cycle Arrest in S-Phase and Promotes Viral Replication

Subversion of the host cell cycle to facilitate viral replication is a common feature of coronavirus infections. Coronavirus nucleocapsid (N) protein can modulate the host cell cycle, but the mechanistic details remain largely unknown. Here, we investigated the effects of manipulation of porcine epidemic diarrhea virus (PEDV) N protein on the cell cycle and the influence on viral replication. Results indicated that PEDV N induced Vero E6 cell cycle arrest at S-phase, which promoted viral replication (P < 0.05). S-phase arrest was dependent on the N protein nuclear localization signal S⁷¹NWHFYYLGTGPHADLRYRT⁹⁰ and the interaction between N protein and p53. In the nucleus, the binding of N protein to p53 maintained consistently high-level expression of p53, which activated the p53-DREAM pathway. The key domain of the N protein interacting with p53 was revealed to be S¹⁷¹RGNSQNRGNNQGRGASQNRGGNN¹⁹⁴ (N^S171-N194), in which G¹⁸³RG¹⁸⁵ are core residues. N^S171-N194 and G¹⁸³RG¹⁸⁵ were essential for N-induced S-phase arrest. Moreover, small molecular drugs targeting the N^S171-N194 domain of the PEDV N protein were screened through molecular docking. Hyperoside could antagonize N protein-induced S-phase arrest by interfering with interaction between N protein and p53 and inhibit viral replication (P < 0.05). The above-described experiments were also validated in porcine intestinal cells, and data were in line with results in Vero E6 cells. Therefore, these results reveal the PEDV N protein interacts with p53 to activate the p53-DREAM pathway, and subsequently induces S-phase arrest to create a favorable environment for virus replication. These findings provide new insight into the PEDV-host interaction and the design of novel antiviral strategies against PEDV. IMPORTANCE Many viruses subvert the host cell cycle to create a cellular environment that promotes viral growth. PEDV, an emerging and reemerging coronavirus, has led to substantial economic loss in the global swine industry. Our study is the first to demonstrate that PEDV N-induced cell cycle arrest during the S-phase promotes viral replication. We identified a novel mechanism of PEDV N-induced S-phase arrest, where the binding of PEDV N protein to p53 maintains consistently high levels of p53 expression in the nucleus to mediate S-phase arrest by activating the p53-DREAM pathway. Furthermore, a small molecular compound, hyperoside, targeted the PEDV N protein, interfering with the interaction between the N protein and p53 and, importantly, inhibited PEDV replication by antagonizing cell cycle arrest. This study reveals a new mechanism of PEDV-host interaction and also provides a novel antiviral strategy for PEDV. These data provide a foundation for further research into coronavirus-host interactions.

Cytokeratin 7 and 19 expression in oropharyngeal and oral squamous cell carcinoma

PURPOSE

The precise etiopathogenesis of human papillomavirus (HPV)-related oropharyngeal squamous cell carcinoma (OPSCC), and reasons for predilection for crypt epithelium, remain uncertain. The purpose of this study is to investigate the interaction between HPV and specific cytokeratins 7 (CK7) and 19 (CK19) in crypt epithelium.

METHODS

This is a retrospective cohort study of patients presenting between 1999 and 2015 at a tertiary referral center. CK7 and CK19 positivity and H Scores were determined by immunohistochemistry. Disease-specific and overall survival rates were analyzed.

RESULTS

There were 253 patients presenting with OPSCC (134), squamous cell carcinoma (SCC) of unknown primary site (22), and oral tongue SCC (97). Primary tumor CK7 and CK19 positivity and H Scores were significantly higher in HPV-positive OPSCC than HPV-negative OPSCC and oral tongue SCC. Higher CK19 Scores, but not CK7 Scores, were also seen in regional metastases from HPV-positive OPSCC than other sites. No impact on disease-specific or overall survival was identified on multivariate analysis.

CONCLUSION

The increased expression of CK7 and CK19 in HPV-positive OPSCC compared to HPV-negative disease supports the theory for a role for these cytokeratins in the etiopathogenesis of HPV-related OPSCC.

Identification of significant genes signatures and prognostic biomarkers in cervical squamous carcinoma via bioinformatic data

Background

Cervical squamous cancer (CESC) is an intractable gynecological malignancy because of its high mortality rate and difficulty in early diagnosis. Several biomarkers have been found to predict the prognose of CESC using bioinformatics methods, but they still lack clinical effectiveness. Most of the existing bioinformatic studies only focus on the changes of oncogenes but neglect the differences on the protein level and molecular biology validation are rarely conducted.

Methods

Gene set data from the NCBI-GEO database were used in this study to compare the differences of gene and protein levels between normal and cancer tissues through significant pathway selection and core gene signature analysis to screen potential clinical biomarkers of CESC. Subsequently, the molecular and protein levels of clinical samples were verified by quantitative transcription PCR, western blot and immunohistochemistry.

Results

Three differentially expressed genes (RFC4, MCM2, TOP2A) were found to have a significant survival (P < 0.05) and highly expressed in CESC tissues. Molecular biological verification using quantitative reverse transcribed PCR, western blotting and immunohistochemistry assays exhibited significant differences in the expression of RFC4 between CESC and para-cancerous tissues (P < 0.05).

Conclusion

This study identified three potential biomarkers (RFC4, MCM2, TOP2A) of CESC which may be useful to clarify the underlying mechanisms of CESC and predict the prognosis of CESC patients.

Genome-Wide Transcriptome Analysis of Human Papillomavirus 16-Infected Primary Keratinocytes Reveals Subtle Perturbations Mostly due to E7 Protein Expression

It is established that the host cell transcriptomes of natural lesions, organotypic rafts, and human papillomavirus (HPV)-immortalized keratinocytes are altered in the presence of HPV genomes. However, the establishment of HPV-harboring cell lines requires selection and immortalization, which makes it impossible to distinguish between alterations directly induced by HPV or indirectly by the need for immortalization or selection. To address direct effects of HPV infection on the host cell transcriptome, we have used our recently established infection model that allows efficient infection of primary keratinocytes with HPV16 virions. We observed only a small set of genes to be deregulated at the transcriptional level at 7 days postinfection (dpi), most of which fall into the category regulated by pocket proteins pRb, p107, and p130. Furthermore, cell cycle genes were not deregulated in cells infected with a virus lacking E7 despite the presence of episomal genome and viral transcripts. These findings imply that the majority of transcriptional changes are due to the E7 protein impairing pocket protein function. Additional pathways, such as the Fanconi anemia-BRCA pathway, became perturbed only after long-term culturing of infected cells. When grown as organotypic raft cultures, keratinocytes infected with wild-type but not E7 mutant virus had perturbed transcriptional regulation of pathways previously identified in natural lesions and in rafts derived from immortalized keratinocytes. We conclude that the HPV infection model provides a valuable tool to distinguish immediate transcriptional alterations from those induced by persistent infection and the need for selection and immortalization.IMPORTANCE To establish infection and complete the viral life cycle, human papillomavirus (HPV) needs to alter the transcriptional program of host cells. Until recently, studies were restricted to keratinocyte-derived cell lines immortalized by HPV due to the lack of experimental systems to efficiently infect primary keratinocytes. Need for selection and immortalization made it impossible to distinguish between alterations induced by HPV and secondary adaptation due to selection and immortalization. With our recent establishment of an extracellular matrix (ECM)-to-cell transfer system allowing efficient infection of primary keratinocytes, we were able to identify transcriptional changes attributable to HPV16 infection. Most perturbed genes fall into the class of S-phase genes, which are regulated by pocket proteins. Indeed, infection with viruses lacking E7 abrogated most transcriptional changes. It is important to note that many transcriptional alterations thought to be important for the HPV life cycle are actually late events that may reflect immortalization and, possibly, disease progression.