Abstract:Objective: Human papillomavirus (HPV) is the most common sexually transmitted agent in the world. HPV6b is a low-risk type of HPVs that causes benign verrucous hyperplastic lesions of the skin and anal genital mucosa. Previous research has indicated that HPV genotype 6 is sometimes associated with high-grade lesions and anal cancer. The pathogenesis of low-risk HPV infection and its relationship to high-risk HPV is not clear at present. The E7 protein, which is encoded by HPV early-expressing genes, plays an important role in HPV infection. The aim of this study is to investigate the human gene expression signature of the HPV6b E7-transfected HaCaT stable cell line. The identification of differentially expressed genes might provide a more comprehensive understanding of HPV6b infection and will allow us to explore the specific role of E7 protein. Methods: We established a stable cell line transfected with the HPV6b E7 gene and analyzed the line's genome-wide expression profile by microarray. Quantitative real-time PCR (qRT-PCR) was used to verify the differentially expressed genes. GO enrichment analysis was applied for gene annotation according to functions. KEGG analysis, a system for analyzing gene function and genome information, was used to help us integrate differentially expressed genes into pathways. Results: A total of 3519 genes were identified to be significantly differentially expressed between the HPV 6bE7-HaCaT stable cell line and a control cell line, among which 1884 genes were up-regulated and 1635 genes were down-regulated with a fold-change > 2.0 between the two groups. The expression profiles of the top 20 up regulated and the top 20 down-regulated genes in the HPV 6bE7-HaCaT stable cell line as analyzed by qRT-PCR were consistent with the microarray data. The most significantly enhanced genes HPV 6bE7-HaCaT cells were SIMC1, S100A8 and S100A9, whereas PXDN expression was markedly down-regulated. GO analysis showed that HPV 6bE7 primarily affected biological processes and that the most significant difference was in heart induction (GO:0003129). Many differentially expressed genes were linked to histone H4-K20 demethylation (GO:0035574). KEGG analysis showed that the most significant changes in gene expression were related to primary bile acid biosynthesis, and the most diverse biological processes were related to systemic lupus erythematosus pathogenesis. Conclusions: The global gene expression profile of the HPV 6bE7-HaCaT stable cell line was analyzed, revealing the genes regulated by E7 protein and providing insight into the pathogenesis of HPV6b infection.

Hpv-2 Isolates from Patients with Huge Verrucae Vulgaris Possess Stronger Promoter Activities

Detection of HPV-2 and identification of novel mutations by whole genome sequencing from biopsies of two patients with multiple cutaneous horns.

A Point Mutation in the DNA-binding Domain of HPV-2 E2 Protein Increases Its DNA-binding Capacity and Reverses Its Transcriptional Regulatory Activity on the Viral Early Promoter.

The enhancement of DNA binding ability of a mutated E2 (A338V) protein of HPV-2]

Digital RNA Sequencing of Human Epidermal Keratinocytes Carrying Human Papillomavirus Type 16 E7

[Mutations in Various Functional Domains of HPV2 E2 Protein Inhibit the Transcriptional Depression Activities].

Large Scale Genotype Comparison of Human Papillomavirus E2-Host Interaction Networks Provides New Insights for E2 Molecular Functions

Development of a Multiplex PCR Method for Detecting and Typing Human Papillomaviruses in Verrucae Vulgaris

HPV18 DNA Replication Inactivates the Early Promoter P 55 Activity and Prevents Viral E6 Expression

Genome-Wide Profiling of Cervical RNA-Binding Proteins Identifies Human Papillomavirus Regulation of RNASEH2A Expression by Viral E7 and E2F1

Molecular Epidemiological Study on Prevalence of Human Papillomaviruses in Patients with Common Warts in Beijing Area.

Identification of Potential Prognostic Biomarkers in Vulval Squamous Cell Carcinoma Based on Human Papillomavirus Infection Status-Analysis of GSE183454

Human papillomavirus type 58 genome variations and RNA expression in cervical lesions.

Functional effects of sequence variations in the E6 and E2 genes of human papillomavirus 16 European and Asian variants.

The possible impact of novel mutations in human papillomavirus 52 on the infection characteristics

The use of a human papillomavirus 18 promoter for tissue-specific expression in cervical carcinoma cells

Specific Up-Regulation of DNA Polymerase by Human Papillomavirus 16

Human Gene Expression Microarray Analysis of the HPV 6Be7-Hacat Stable Cell Line.

Human papillomaviruses 53 and 66: clinical aspects and genetic analysis

Genetic variability, phylogeny and functional implication of the long control region in human papillomavirus type 16, 18 and 58 in Chengdu, China

Transgenic HPV11-E2 protein modulates URR activity in vivo