Variant upstream regulatory region sequences differentially regulate human papillomavirus type 16 DNA replication throughout the viral life cycle

Prognostic and therapeutic potential of STAT3: Opportunities and challenges in targeting HPV-mediated cervical carcinogenesis

Cervical cancer (CaCx) ranks as the fourth most prevalent cancer among women globally. Persistent infection of high-risk human papillomaviruses (HR-HPVs) is major etiological factor associated with CaCx. Signal Transducer and Activator of Transcription 3 (STAT3), a prominent member of the STAT family, has emerged as independent oncogenic driver. It is a target of many oncogenic viruses including HPV. How STAT3 influences HPV viral gene expression or gets affected by HPV is an area of active investigation. A better understanding of host-virus interaction will provide a prognostic and therapeutic window for CaCx control and management. In this comprehensive review, we delve into carcinogenic role of STAT3 in development of HPV-induced CaCx. With an emphasis on fascinating interplay between STAT3 and HPV genome, the review explores the diverse array of opportunities and challenges associated with this field to harness the prognostic and therapeutic potential of STAT3 in CaCx.

Investigation of molecular mechanisms underlying JAK/STAT signaling pathway in HPV-induced cervical carcinogenesis using 'omics' approach

The precise mechanism of action of Janus Kinases (JAK)/Signal Transducer and activator of Transcription (STAT) signaling in human papillomavirus (HPV)-associated cervical cancer (CaCx) is poorly defined. The present study dissected the underlying components of JAK/STAT signaling in HPV-positive cervical neoplasms. Whole transcriptome profile of CaCx cohort from TCGA database revealed elevated STAT3 and its impact on CaCx patients' survival. Using the RT² Profiler PCR Array, we analyzed 84 genes of interest associated with JAK/STAT signaling in mRNA derived from HPV-negative and HPV-positive cervical lesions which revealed 21 differentially expressed genes (DEGs). Analyses of DEGs using the Database for Annotation, Visualization and Integrated Discovery tool indicated maximum genes enriched in immune response and negative regulation of apoptotic process. Protein-protein network analysis indicated IL4, STAT5A, STAT4, and JAK3 to be the key genes in the interaction network. Further, 7 key DEGs (IL4R, IRF1, EGFR, OAS1, PIAS1, STAT4, and STAT5A) were validated in TCGA cohort using R2 platform. These genes were differentially expressed among HPV-positive cervical tissues and their correlation with STAT3 was established. EGFR and IL4R showed a comparatively strong correlation with STAT3 that supports their involvement in pathogenesis of CaCx. Finally, the Kaplan-Meier analysis established the prognostic association of the key DEGs, in CaCx cohort. The STAT3 and associated key genes discovered from our study establish a strong pathogenic role of JAK/STAT3 pathway in HPV-mediated cervical carcinogenesis.

Genetic variations in the long control region of human papillomavirus type 16 isolates from India: implications for cervical carcinogenesis

Introduction. Infection with high-risk human papillomavirus (HPV) types, specifically HPV type 16 (HPV16), is considered to be the most important risk factor in the development of cervical intraepithelial neoplasia and cancer. The long control region (LCR) is a noncoding region that comprises approximately 10 % of the HPV genome and contains regulatory elements for viral transcription and replication. Sequence variations in LCR may impact on the replication efficiency and oncogenic potential of the virus.Gap statement. Studies documenting variations in LCR of HPV16 isolates pertaining to cervical neoplastic status in India are limited.Aim. The present study was designed to characterize variations in the LCR of Indian isolates of HPV16 and study their association with cervical disease grades.Methodology. The LCR was amplified and sequenced from HPV16 positive cervical samples belonging to different cervical disease grades. Sequences were aligned to identify variations and potential transcription factor binding sites (TFbs) were predicted using the JASPAR database in addition to phylogenetic studies.Results. Among the 163 HPV16 isolates analysed, 47 different nucleotide variations were detected in the LCR, of which 25 are reported for first time in Indian isolates. Point mutations were detected in 35/54 (64.8 %) samples with normal cervical status, 44/50 (88 %) samples with low-grade cervical disease and 53/59 (89.8 %) samples with high-grade cervical disease. Variations T6586C, G6657A and T6850G were significantly associated with high-grade cervical status. Thirteen LCR variations were detected in the binding sites for CEBPB, ETS1, JUN, MYB, NFIL3, PHOX2A and SOX9 transcription factors.Conclusion. The present study helped to identify unique variations in the LCRs of HPV16 Indian isolates. The variations in the A4 sub-lineage were significantly associated with high-grade disease status. The isolates belonging to the A4 and D3 sub-lineages harboured mutations in putative TFbs, implying a potential impact on viral replication and progression to cervical cancer.

Complete genome and phylogenetic analysis of bovine papillomavirus type 15 in Southern Xinjiang dairy cow

Background

Bovine papilloma is a neoplastic disease caused by bovine papillomaviruses (BPVs), which were recently divided into 5 genera and at least 24 genotypes.

Objectives

The complete genome sequence of BPV type 15 (BPV Aks-02), a novel putative BPV type from skin samples from infected cows in Southern Xinjiang China, was determined by collecting warty lesions, followed by DNA extraction and amplicon sequencing.

Methods

DNA was analyzed initially by polymerase chain reaction (PCR) using the degenerate primers FAP59 and FAP64. The complete genome sequences of the BPV Aks-02 were amplified by PCR using the amplification primers and sequencing primers. Sequence analysis and phylogenetic analysis were performed using bio-informatic software.

Results

The nucleotide sequence of the L1 open reading frame (ORF) of BPV Aks-02 was 75% identity to the L1 ORF of BPV-9 reference strain from GenBank. The complete genome consisted of 7,189 base pairs (G + C content of 42.50%) that encoded 5 early (E8, E7, E1, E2, and E4) and 2 late (L1 and L2) genes. The E7 protein contained a consensus CX₂CX₂₉CX₂C zinc-binding domain and a LxCxE motif. Among the different members of this group, the percentages of the complete genome and ORFs (including 5 early and 2 late ORFs) sequence identity of BPV Aks-02 were closer to the genus Xipapillomavirus 1 of the Xipapillomavirus genus. Phylogenetic analysis and sequence similarities based on the L1 ORF of BPV Aks-02 revealed the same cluster.

Conclusions

The results suggest that BPV type (BPV Aks-02) clustered with members of the Xipapillomavirus genus as BPV 15 and were closely related to Xipapillomavirus 1.

Characterization of Bovine papillomavirus 28 (BPV28) and a novel genotype BPV29 associated with vulval papillomas in cattle

Papillomavirus (PV) infections are associated with the development of cutaneous and mucosal tumors in humans and various animal species. In humans, infection of high-risk human PVs (HPVs) causes anogenital cancers, while in animals, anogenital-associated PVs are not well understood. Among animal PVs, Bos taurus PVs (BPVs) have the most diverse genotypes, up to 28 of them. The present study will report two unique BPVs identified in vulval papilloma lesions from two Holstein Friesian cattle by conventional PCR and sequencing. In the first case, BPV28 harboring two L1 open reading frames (ORFs) due to a five-nucleotide deletion was identified. In the second case, histologically diagnosed as papilloma, an unclassified BPV genotype was detected. However, in both cases, the immunohistochemistry against PV antigen was negative. The full genome of the unclassified BPV was amplified by inverse PCR and analyzed by genome-walking sequencing. The L1 nucleotide sequence was most identical to BPV genotype 6 (BPV6), showing 78 % identity, indicating that this novel BPV should be classified as species Xipapillomavirus 1, genotype BPV29. The mRNA expression of three early genes (E1, E2, E10), but not L1, was confirmed in both BPV28- and BPV29-detected papilloma lesions. The present study suggests the involvement of novel types of BPV in vulval papilloma. The alteration of BPV28 pathogenicity due to the frameshift mutation of L1 needs to be elucidated in the future.

Genomic characterization of a novel bovine papillomavirus type 28

Infection of bovine papillomavirus (BPV) has been associated with mucosal and/or cutaneous tumor development in bovids. To date, up to 27 genotypes of BPVs have been identified and classified based on the nucleotide sequence identity of L1 open reading frame. In the present study, the complete sequence of a novel BPV concurrently identified with BPV1 and BPV2 in the facial cutaneous papilloma lesion of a domestic cattle was characterized. The whole genome of the unclassified BPV was 7263 base pairs in full length with GC ratio of 42.9%. In comparison with published BPV sequences, L1 nucleotide sequence of the novel BPV shared 75% identity with BPV15, and was suggested to be classified in the genus, Xipapillomavirus. According to the criteria established by the International Committee on the Taxonomy of Viruses, the novel BPV was designated as BPV type 28.

Level of phospho-STAT3 (Tyr705) correlates with copy number and physical state of human papillomavirus 16 genome in cervical precancer and cancer lesions

Our earlier studies indicated an important role of inducible transcription factor STAT3 in the establishment of persistent infection of human papillomavirus (HPV) type 16 and promotion of cervical carcinogenesis. Since HPV load and its physical state are two potential determinants of this virally-induced carcinogensis, though with some exceptions, we extended our study to examine the role of active STAT3 level in cervical precancer and cancer lesions and it’s association with HPV viral load and physical state. An elevated level of active STAT3 was measured by assessing phospho-STAT3-Y705 (pSTAT3), in tumor tissues harboring higher viral load irrespective of the disease grade. Physical state analysis of HPV16 by assessing the degree of amplification of full length E2 and comparing it with E6 (E2:E6 ratio), which predominantly represent episomal form of HPV16, revealed low or undetectable pSTAT3. A strong pSTAT3 immunoreactivity was found in tissues those harbored either mixed or predominantly integrated form of viral genome. Cumulative analysis of pSTAT3 expression, viral load and physical state demonstrated a direct correlation between pSTAT3 expression, viral load and physical state of HPV. The study suggests that there exists a strong clinical correlation between level of active STAT3 expression and HPV genome copy number, and integrated state of the virus that may play a pivotal role in promotion/maintanence of tumorigenic phenotype.

Characterization of an HPV33 natural variant with enhanced transcriptional activity suggests a role for C/EBPβ in the regulation of the viral early promoter

The Long Control Region (LCR) of the human papillomavirus (HPV) genome encompasses the early promoter (EP) that drives expression of the viral oncogenes in infected cells and HPV-associated cancers. Here, we report on a natural variant of HPV33 that displays higher EP activity than the prototype in transfected C33A and HeLa cervical carcinoma cells, and in the osteosarcoma U2OS cell line which supports replication of HPV episomes. This increased promoter activity was ascribed to a single nucleotide variation in the LCR, T7791C, in a putative binding site for the transcription factor C/EBPβ. T7791C abrogated binding of recombinant C/EBPβ to this site in vitro and stimulated the EP in vivo, suggesting that it abrogates a negatively-acting regulatory element. A second C/EBPβ binding site was identified in vitro that activated the EP in vivo and whose function and location in the epithelial-specific enhancer is shown to be conserved in the highly prevalent HPV18. These results suggest that C/EBPβ is both an activator and a repressor of the HPV33 EP, acting via two distinct binding sites. Prediction of C/EBPβ sites in the LCR of 186 HPV types suggests that C/EBPβ regulation of the EP is common among high‐risk viruses from the α genus.

Phylogenetic analysis of Alphapapillomavirus based on L1, E6 and E7 regions suggests that carcinogenicity and tissue tropism have appeared multiple times during viral evolution

Members of the Alphapapillomavirus genus are causative agents for cervix cancer and benign lesions in humans. These viruses are classified according to sequence similarities in their L1 region. Yet, viral carcinogenicity has been associated with variations in the proteins encoded by the E6 and E7 genes. In order to relate evolutionary history with origin of carcinogenicity, we performed phylogenetic reconstructions using both nucleotide and predicted amino acid sequences of the L1, E6 and E7 genes. Whilst phylogenetic analysis of L1 reconstructed genus evolutionary history, phylogenies based on E6 and E7 proteins support the idea that mutations at amino acids S/Tx [V/L] (E6) and LxCxE (E7) might be responsible for carcinogenic potential. These findings indicate that virulence within Alphapapillomavirus have appeared multiple times during evolution. Our results reveal that oncogenic potential is not a monophyletic clade-specific adaptation but might be the result of positive selection on random mutations occurring on proteins involved in host infection during viral diversification.

Regulation of HPV transcription

Human papillomavirus infection is associated with the development of malignant and benign neoplasms. Approximately 40 viral types can infect the anogenital mucosa and are categorized into high- and low-risk oncogenic human papillomavirus, depending on their association with the development of cervical carcinoma. High-risk human papillomavirus 16 and 18 are detected in 55% and 15% of all invasive cervical squamous cell carcinomas worldwide, respectively. Low-risk human papillomavirus 6 and 11 are responsible for 90% of genital warts and are also associated with the development of recurrent respiratory papillomatosis. Human papillomavirus preferentially infects mitotic active cells of the basal layer from both mucosal and cutaneous epithelium through microabrasions. The viral life cycle synchronizes with the epithelial differentiation program, which may be due, in part, to the binding of differentially expressed cellular transcription factors to the long control region throughout the various epithelial layers. This review aimed to summarize the current knowledge regarding the mechanisms by which viral gene expression is regulated and the influence of human papillomavirus heterogeneity upon this phenomenon. A better understanding of the regulatory mechanisms may elucidate the particularities of human papillomavirus-associated pathogenesis and may provide new tools for antiviral therapy.

Naturally Occurring Variations Modulate the Activity of the HPV33 Early Promoter and its Affinity for the E2 Transcription Factor

The human papillomavirus (HPV) Long Control Region (LCR) encompasses the early promoter (EP) that drives transcription of the E6 and E7 oncogenes in keratinocytes and HPV-associated cancers. In this study, the transcriptional activities of the HPV33 EP from the prototype LCR and from eight variants representative of the worldwide diversity of the virus were examined in primary human keratinocytes (PHK) and in the HeLa cervical carcinoma cell line by luciferase reporter-gene assays. Remarkably, the two variations with the greatest effect on the EP in PHK were C7732G and a 79-bp deletion that were associated with high-grade cervical lesions and persistent infections, respectively, in epidemiological studies. In contrast, the three variations most active in HeLa cells were C7537A, A7874C and A7879G. A7874C, which lies within an E2-binding sequence, is also shown to increase the activity and binding of E2 at this site. Collectively, these results indicate that naturally-occurring variations affect the HPV33 EP differentially in PHK than in cancer cells and, furthermore, that they can also alter its regulation by E2. These findings provide a molecular basis for rationalizing the results of previous epidemiological studies and for understanding the contribution of LCR polymorphisms to the oncogenicity and persistence of HPV33 infections.

The association of human papillomavirus type 16 E2 variations with cervical cancer in a Han Chinese population

Human papillomavirus type 16 (HPV16) is considered to be the primary pathogen related to cervical cancer. The HPV16 E2 protein plays an important role in tumourigenicity of cervical carcinoma. In the current study, we enrolled 121 HPV16-positive cervical cancer patients in the case group and 130 HPV16-positive asymptomatic individuals in the control group, and we investigated the association between HPV16 E2 gene variations and cervical cancer. The HPV16 E2 DNA was amplified and sequenced. We identified two HPV variants (EUR and As) in the control group; the As variant was predominant (68.5%), followed by the EUR variant (31.5%). In the case group, three HPV variants (EUR, As and AA) were observed; the As variant was predominant (72.7%), followed by the EUR variant (22.3%) and the AA variant (5.0%). Our results showed a significant difference in the distribution of the HPV16 variants between the case and control groups (P < 0.05). Moreover, in the HPV16 E2 gene variation analysis, the distribution of sixteen variations was significantly different between the case and control groups (P < 0.05), and all of these variations were present in the AA variant. In the subgroup analysis, the frequency of the T3575G (S274A) variation in the EUR variant was significantly different between the case and control groups (P = 0.029); however, there was no significant difference in the frequency of the variations in the As variant between the case and control groups. Our findings in the current study could provide a better understanding of the relationship between HPV16 variants, E2 gene variations and cervical cancer.

Phylogeny and polymorphism in the long control regions E6, E7, and L1 of HPV Type 56 in women from southwest China

Globally, human papillomavirus (HPV)‑56 accounts for a small proportion of all high‑risk HPV types; however, HPV‑56 is detected at a higher rate in Asia, particularly in southwest China. The present study analyzed polymorphisms, intratypic variants, and genetic variability in the long control regions (LCR), E6, E7, and L1 of HPV‑56 (n=75). The LCRs, E6, E7 and L1 were sequenced using a polymerase chain reaction and the sequences were submitted to GenBank. Maximum‑likelihood trees were constructed using Kimura's two‑parameter model, followed by secondary structure analysis and protein damaging prediction. Additionally, in order to assess the effect of variations in the LCR on putative binding sites for cellular proteins, MATCH server was used. Finally, the selection pressures of the E6‑E7 and L1 genes were estimated. A total of 18 point substitutions, a 42‑bp deletion and a 19‑bp deletion of LCR were identified. Some of those mutations are embedded in the putative binding sites for transcription factors. 18 single nucleotide changes occurred in the E6‑E7 sequence, 11/18 were non‑synonymous substitutions and 7/18 were synonymous mutations. A total 24 single nucleotide changes were identified in the L1 sequence, 6/24 being non‑synonymous mutations and 18/24 synonymous mutations. Selective pressure analysis predicted that the majority of mutations of HPV‑56 E6, E7 and L1 were of positive selection. The phylogenetic tree demonstrated that the isolates distributed in two lineages. Data on the prevalence and genetic variation of HPV‑56 types in southwest China may aid future studies on viral molecular mechanisms and contribute to future investigations of diagnostic probes and therapeutic vaccines.

The polymorphisms of LCR, E6, and E7 of HPV-58 isolates in Yunnan, Southwest China

Backgroud

Variations in HPV LCR/E6/E7 have been shown to be associated with the viral persistence and cervical cancer development. So far, there are few reports about the polymorphisms of the HPV-58 LCR/E6/E7 sequences in Southwest China. This study aims to characterize the gene polymorphisms of the HPV-58 LCR/E6/E7 sequences in women of Southwest China, and assess the effects of variations on the immune recognition of viral E6 and E7 antigens.

Methods

Twelve LCR/E6/E7 of the HPV-58 isolates were amplified and sequenced. A neighbor-joining phylogenetic tree was constructed by MEGA 7.0, followed by the secondary structure prediction of the related proteins using PSIPRED v3.3. The selection pressure acting on the HPV-58 E6 and E7 coding regions was estimated by Bayes empirical Bayes analysis of PAML 4.8. Meanwhile, the MHC class-I and II binding peptides were predicted by the ProPred-I server and ProPred server. The transcription factor binding sites in the HPV-58 LCR were analyzed using the JASPAR database.

Results

Twenty nine SNPs (20 in the LCR, 3 in the E6, 6 in the E7) were identified at 27 nucleotide sites across the HPV-58 LCR/E6/E7. From the most variable to the least variable, the nucleotide variations were LCR > E7 > E6. The combinations of all the SNPs resulted in 11 unique sequences, which were clustered into the A lineage (7 belong to A1, 2 belong to A2, and 2 belong to A3). An insertion (TGTCAGTTTCCT) was found between the nucleotide sites 7280 and 7281 in 2 variants, and a deletion (TTTAT) was found between 7429 and 7433 in 1 variant. The most common non-synonymous substitution V77A in the E7 was observed in the sequences encoding the α-helix. 63G in the E7 was determined to be the only one positively selected site in the HPV-58 E6/E7 sequences. Six non-synonymous amino acid substitutions (including S71F and K93 N in the E6, and T20I, G41R, G63S/D, and V77A in the E7) were affecting multiple putative epitopes for both CD4⁺ and CD8⁺ T-cells. In the LCR, C7265G and C7266T were the most variable sites and were the potential binding sites for the transcription factor SOX10.

Conclusion

These results provide an insight into the intrinsic geographical relatedness and biological differences of the HPV-58 variants, and contribute to further research on the HPV-58 epidemiology, carcinogenesis, and therapeutic vaccine development.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-0986-7) contains supplementary material, which is available to authorized users.

Different Association of Human Papillomavirus 16 Variants with Early and Late Presentation of Cervical Cancer

The median age of cervical cancer (CC) presentation coincides with the mean age of menopause presentation (49 years) in Mexico. Here, we investigated the association between different HPV16 variants and early (≤ 49 years) or delayed (≥ 50 years) CC presentation. We conducted a case-case study that included 462 CCs, 386 squamous cell carcinomas (SCC), 63 adenocarcinomas (ACC), and 13 additional cell types. Variants were identified by PCR and DNA sequencing. The risk conferred by each variant for developing CC earlier than 50 years was analyzed using a univariate logistic regression model considering old-aged patients (≥ 50 years) and non-HPV16 cases as the reference variables. Overall, the frequency of HPV16 was 50.9%, and the only identified variants were the European A1/2 (31.2%) and the Asian-American D2 (10.8%), and D3 (8.9%). D2 was mainly associated with ≤ 49-year-old patients (15.9%); A1/2 was uniformly distributed between the two age groups (~31%), whereas D3 increased with age to a frequency of 11.8% in the older group. Only the D2 variant conferred a 3.3-fold increase in the risk of developing CC before 50 years of age (OR = 3.3, 95% CI = 1.7-6.6, p < 0.001) in relation with non-HPV16 cases. Remarkably, this risk was higher for ACC (OR = 6.0, 95% CI = 1.1-33, p < 0.05) than for SCC (OR = 2.8, 95% CI = 1.3-5.9, p < 0.01). Interestingly, when analyzing only the HPV16-positive CC, D2 increases (OR = 2.5, 95% CI = 1.2-5, p < 0.05) and D3 decreases (OR = 0.45, 95% CI 0.2-0.9, p < 0.05) the risk to develop CC before 50 years old in relation with A1/2 variant. These results indicated that D2 variant is associated with early and D3 with delayed CC presentation, whereas A1/2 variant was uniformly distributed between the two age groups.

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.

Genetic variation of high-risk HPV-16 from a cervical cancer biopsy collected in Pakistan

Aim: HPV-16 is one of the causative agents of cervical cancer. The HPV-16 genome is characterized widely to assist mutagenesis studies, phylogenetic analysis and protein annotations. This is the first report of HPV-16 sequence analysis from Pakistan. Materials & methods: The gene and long control region specific primers were used to generate the targeted amplicons that were then cloned and sequenced. Results: The complete genome analysis revealed nine changes both at amino acid and nucleotide level including seven novel and two previously reported mutations. The evolutionary relationship of HPV from Pakistan and reported sequences shows that strain reported in this study segregates with European taxa. Conclusion: This report lays foundation to conduct HPV genome mutation analysis in Pakistan.

Differences in Transcriptional Activity of Human Papillomavirus Type 6 Molecular Variants in Recurrent Respiratory Papillomatosis

A significant proportion of recurrent respiratory papillomatosis (RRP) is caused by human papillomavirus type 6 (HPV-6). The long control region (LCR) contains cis-elements for regulation of transcription. Our aim was to characterize LCR HPV-6 variants in RRP cases, compare promoter activity of these isolates and search for cellular transcription factors (TFs) that could explain the differences observed. The complete LCR from 13 RRP was analyzed. Transcriptional activity of 5 variants was compared using luciferase assays. Differences in putative TFs binding sites among variants were revealed using the TRANSFAC database. Chromatin immunoprecipation (CHIP) and luciferase assays were used to evaluate TF binding and impact upon transcription, respectively. Juvenile-onset RRP cases harbored exclusively HPV-6vc related variants, whereas among adult-onset cases HPV-6a variants were more prevalent. The HPV-6vc reference was more transcriptionally active than the HPV-6a reference. Active FOXA1, ELF1 and GATA1 binding sites overlap variable nucleotide positions among isolates and influenced LCR activity. Furthermore, our results support a crucial role for ELF1 on transcriptional downregulation. We identified TFs implicated in the regulation of HPV-6 early gene expression. Many of these factors are mutated in cancer or are putative cancer biomarkers, and must be further studied.

Human papillomavirus type 56 polymorphism in Canadian women with and without cervical lesions

BACKGROUND

The genomic diversity of high-risk human papillomaviruses (HPV) has been associated with viral persistence and HPV-induced lesions. Studies on HPV56 persistence are still pending.

OBJECTIVE

To assess the association between HPV56 polymorphism and HPV56 persistence and presence of high-grade cervical intraepithelial neoplasia (CIN2,3) or cancer.

STUDY DESIGN

HPV56-positive cervical specimens from 204 women selected from a total of 4669 participants recruited in 5 epidemiological studies (parent studies) were further analyzed by PCR-sequencing of the long control region (LCR).

RESULTS

Of the 81 women followed prospectively in cohort studies who could be classified, 34 had persistent and 47 had transient HPV56 infections. Variant HPV56-LCR-MTL-21 was detected more frequently in persistent infections (52.9%, 95% CI: 36.7-68.6%) than in transient infections (25.5%, 95% CI: 15.1-39.4). Considering only women recruited in a cohort of women infected or at high risk for HIV infection, infection with variant HPV56-LCR-MTL-21 (OR=4.4, 95% CI: 1.3-14.5) was significantly associated with HPV56 persistence controlling in multivariate analysis for high risk HPV detection and HIV infection. A variation at nucleotide 7800 in HPV56-LCR-MTL-21 resulted in the loss of a binding site for Elf-1 embedded in one of the E2 binding sites, a potential activator or repressor of expression of the HPV genome. HPV56 polymorphism was not associated with CIN2,3 or cancer in women enrolled in cross-sectional and case-control studies.

CONCLUSION

Polymorphism in HPV56 may influence the risk that infections with this type will persist.

Genetic variations of E6 and long control region of human papillomavirus type 16 from patients with cervical lesion in Liaoning, China

Background

High-risk human papillomavirus type 16 (HPV16) is a risk factor for cervical cancer. Previous studies suggest that polymorphisms in the E6 gene or the long control region(LCR)of HPV16 may alter the oncogenic potential of the virus. The aims of this study were to investigate the genetic variations of HPV16 E6 gene and LCR in isolates from Chinese population and correlation of the E6 and LCR polymorphisms with disease status of infected patients.

Methods

HPV16 positive endocervical specimens were collected from 304 women living in Northeast of China. Sequences of E6 gene and LCR were analyzed by PCR-sequencing.

Results

Two lineages were found in the populations, including EUR lineage and As lineage. Based on the HPV16 prototype, the most frequent variation in the E6 gene was T178A/G (48.7%), followed by mutations of G94A (12.2%) and T350G (9.9%). The rank orders of incidence of E6 variations in amino acid were as follows: D25E (46.3%), L83V (9.9%) and H78Y (4.3%). Nucleotide variations in LCR were found in all the 304 isolates from HPV16 positive cervical samples. The most commonly observed LCR variations were the transition replacement G7193T, 7434CIns, G7521A and 7863ADel (100%). The As lineage was associated with HPV persistent infections and with disease status of ≥CIN2,3. The EUR lineage variants showed a negative trend of association with the severity of ≥CIN2,3. Among 41 variations found in LCR, 25 (61.0%) were located at the binding sites for transcription factors. Occurrence of ≥CIN2,3 was significantly associated with the mutations of R10G/L83V in E6 and the C7294T co-variation in LCR, after adjusting for ages of infected patients.

Conclusions

Associations between As lineage and HPV persistent infections, and with disease status of ≥CIN2,3, and an association between the EUR lineage and negative trend of association with the severity of ≥CIN2,3 were found in this study. An association between a co-variation of R10G/L83V in E6 and C7294T in LCR and an increased risk for developing CIN-2,3 was found in a HPV16 infected population of Chinese women. These findings indicate that HPV16 polymorphism influences development of CIN-2,3.

Whole genome analysis of human papillomavirus genotype 11 from cervix, larynx and lung

The prevalence of human papillomavirus genotypes differs in various target organs. HPV16 is the most prevalent genotype in the cervix while genotypes 6 and 11 are highly prevalent in skin and aero-digestive tract infections. In this study HPV11 positive specimens were selected from cervix, larynx and lung biopsy tissue to analyze the whole genome by PCR and direct sequencing. Five HPV11 whole genomes were characterized, consisting of two cervical specimens, two laryngeal specimens and one lung specimen. The results showed high homology of HPV11 in these organs. Phylogenetic analysis showed that all HPV11 derived from various organs belonged to the same lineage. Molecular characterization and functional studies can further our understanding of virulence, expression or transmission. Additional studies on functional protein expression at different organ sites will also contribute to our knowledge of HPV infection in various organs.

Cell-type-dependent activities of regulatory regions and E2 proteins derived from carcinogenic and non-carcinogenic human alphapapillomaviruses

A large number of studies have revealed that persistent infections with certain human papillomavirus (HPV) types are necessary for the development of invasive cancer of the cervix. Recent studies have shown that not only do the major carcinogenic HPV types 16 and 18 encode E6 and E7 oncoproteins with immortalizing activity but also the very weakly or non-carcinogenic types 53, 66, 70 and 82. Currently, it is unknown whether transcriptional differences exist between these viruses that account for carcinogenicity in vivo. Therefore, we compared for the first time the activities of the upstream regulatory regions (URRs) that drive E6 and E7 expression derived from HPV16, -18, -31, -53, -66, -70 and -82 in the absence and presence of the viral E2 transcriptional regulator. URR activities in the absence of E2 varied widely and were further modulated by the cellular background. The co-expression of homologous E2 proteins resulted in repression of the URRs of only some HPV types and this varied with cell type. Activation by E2 proteins was less cell-type dependent but differed in an HPV-type-dependent manner. However, basal URR activity, repression of the URR by E2 and transcriptional activation by E2 did not correlate with HPV carcinogenicity in vivo. In summary, our data do not support the model that the transcriptional activity of human alphapapillomavirus types correlates with epidemiological risk classification.

Human papillomavirus type 52 polymorphism and high-grade lesions of the uterine cervix

The association between polymorphism of human papillomavirus type 52 (HPV52) and high-grade cervical intraepithelial neoplasia (CIN2,3) was investigated in Canadian women. HPV-52-positive endocervical specimens collected from 216 women selected from a total of 3,614 participants recruited in two case-control and two cohort studies conducted in Canada, were further analyzed by PCR-sequencing of the LCR and E6 gene. Overall, the HPV52 LCR prototype was detected more frequently in Caucasian women (69 of 132, 52.3%, 95% confidence interval (CI): 43.8%-60.6%) than in non-Caucasian women (15 of 48, 31.3%, 95% CI 19.9%-45.4%). In two cohort studies, HPV52 prototype was detected in seven of 15 (46.7%, 95% CI 24.8-69.9) HPV52 persistent infections and 14 of 35 (40.0%, 95% CI 25.5-56.5) transient infections (p = 0.76). In two case-control studies, 30 participants did not have CIN, 18 had low-grade CIN (CIN1), 64 had CIN2,3, seven had cervical cancer and the diagnosis was undefined for 27 women. Variant MTL-52-LCR-02 was detected more frequently in women with cancer (28.6%, 95% CI 7.6%-64.8%) than in women without cancer or CIN2,3 (0%, 95% CI 0.0%-9.2%; p = 0.015). CIN2,3 risk was significantly associated with a deletion at nucleotide position 7695 in the LCR (OR 4.9, 95% CI 1.2-20.8), the T7744C variation in the LCR (OR 5.7, 95% CI 1.1-32.0), and the K93R variation in E6 (OR 6.9, 95% CI 1.3-36.8), after adjusting for age, detection of HPV16 or 18 and study site. These findings indicate that HPV52 polymorphism influences risk of CIN-2,3 and possibly invasive cancer.

Genomic polymorphism of human papillomavirus type 52 in women from Northeast China

Human papillomavirus (HPV) 52 is an oncogenic HPV type prevalent in Asia. The aim of the study was to analyze HPV 52 genetic variations in women from Northeast China. To explore the intratypic variants of HPV 52, the genomic regions of L1, E6, E7 and long control region (LCR) of HPV 52, which have been identified in women from Northeast China by HPV GenoArray test, were analyzed. Twenty-five mutations were identified in the regions examined. Of the mutations found in the L1 gene, three novel nonsynonymous mutations of C5640T, A5641T and G5642A were located within the region that encodes the binding domain of neutralizing antibodies against HPV 52. Although four variations were identified in HPV 52 E6 and E7 genes, no significant association was found between the mutations and the cytological lesion of the patients. Eight mutations, including a novel CTT7681−7683 deletion, found in the LCR of HPV 52 encompassed the known transcription binding sites, which may possibly affect the transcription of the oncogenic genes of E6 and E7. The most prevalent HPV 52 variant in women from northeastern China belongs to clade L1-LN-A. The genetic variations of HPV 52, including three novel nonsynonymous mutations of C5640T, A5641T and G5642A in the L1 gene and a novel CTT7681−7683 deletion in the LCR, were first documented in strains from women in Northeast China. The statistical result showed no associations between the variants and the severities of the infected women. These findings provide new data regarding gene variations of HPV 52.

Intratypic changes of the E1 gene and the long control region affect ori function of human papillomavirus type 18 variants

A persistent infection with high-risk human papillomavirus (HPV) constitutes the main aetiological factor for cervical cancer development. HPV16 and 18 are the most prevalent types found in cervical cancer worldwide. It has been proposed that HPV intratype variations may result in differences in biological behaviour. Three different HPV18 variants belonging to the Asian Amerindian (AsAi), European (E) and African (Af) branches have been associated with specific histological types of cervical cancer with different relative prognoses, suggesting that HPV18 genomic variations might participate in disease evolution. The E1 viral protein plays a critical role in controlling viral replication and load, requiring interaction with the E2 protein to bind to the long control region (LCR). In this work, we analysed if intratype variations in the LCR and E1 and E2 genes of HPV18 impact ori replication. While the changes found in E2 genes of the tested variants were irrelevant in replication, we found that variations in E1 and LCR in fact affect ori function. It was demonstrated that nucleotide differences in the LCR variants impact ori function. Nevertheless, HPV18 E1 Af gene was mainly involved in the highest ori replication, compared with the E and AsAi E1 variants. Immunofluorescence analysis showed increased levels of Af E1 in the nucleus, correlating with the enhanced ori function. Site-directed mutagenesis revealed that at least two positions in the N-terminal domain of E1 could impact its nuclear accumulation.

Subtypes of associated protein-DNA (Transcription Factor-Transcription Factor Binding Site) patterns

In protein–DNA interactions, particularly transcription factor (TF) and transcription factor binding site (TFBS) bindings, associated residue variations form patterns denoted as subtypes. Subtypes may lead to changed binding preferences, distinguish conserved from flexible binding residues and reveal novel binding mechanisms. However, subtypes must be studied in the context of core bindings. While solving 3D structures would require huge experimental efforts, recent sequence-based associated TF-TFBS pattern discovery has shown to be promising, upon which a large-scale subtype study is possible and desirable. In this article, we investigate residue-varying subtypes based on associated TF-TFBS patterns. By re-categorizing the patterns with respect to varying TF amino acids, statistically significant (P values ≤ 0.005) subtypes leading to varying TFBS patterns are discovered without using TF family or domain annotations. Resultant subtypes have various biological meanings. The subtypes reflect familial and functional properties and exhibit changed binding preferences supported by 3D structures. Conserved residues critical for maintaining TF-TFBS bindings are revealed by analyzing the subtypes. In-depth analysis on the subtype pair PKVVIL-CACGTG versus PKVEIL-CAGCTG shows the V/E variation is indicative for distinguishing Myc from MRF families. Discovered from sequences only, the TF-TFBS subtypes are informative and promising for more biological findings, complementing and extending recent one-sided subtype and familial studies with comprehensive evidence.

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.

Novel HPV-6 variants of human papillomavirus causing recurrent respiratory papillomatosis in southern Africa

There is currently no information regarding the genetic diversity of HPV-6 variants circulating in South Africa. The aim of this study was to determine the HPV-6 variants affecting patients with recurrent respiratory papillomatosis, to determine whether mutations correlate with disease severity and identify molecular determinants of virulence with prognostic relevance. HPV-6 variants were identified based on genome changes within the 712-991 bp region encompassing the non-coding region (URR) of the genome, with variations in length resulting from insertions and duplications, and the 453-bp gene encoding the E6 protein. Based on manual comparison of sequence data from the URR, the isolates were identified as HPV-6a and HPV-6vc variants. Three novel HPV-6 variants were identified: one based on a mutation in the E6 region; two based on changes in the URR including a unique substitution detected in three isolates and an insertion and 170-bp duplication in the URR genome in one patient, who had clinical features of severe disease.

Human papillomavirus type 16 variant analysis of E6, E7, and L1 [corrected] genes and long control region in [corrected] cervical carcinomas in patients in northeast China

Human papillomavirus type 16 (HPV 16) plays a cardinal role in the pathogenesis of cervical cancer. HPV 16 has intratypic variants which show different geographical distributions and different oncogenic potentials. To analyze the presence of sequence variations of HPV 16 variants in northeast China, 71 cervical carcinomas were identified by HPV typing. HPV 16-positive specimens were analyzed by PCR-directed sequencing in the E6, E7, and L1 genes and the LCR (long control region). The variation data were compared with those of neighboring districts. In this hospital-based study, HPV 16 was the most common type (73.24%). In HPV 16-positive specimens, 67.31% belonged to the European (E) lineage, while 32.69% were Asian (As) variants. The Asian-American (AA), African-1 (Af-1), African-2 (Af-2), and northern American (NA) lineages were not detected. The most frequently observed variation sites were T178G (32.69%) in E6; A647G (34.62%), G666A (38.46%), and T846C (32.69%) in E7; C6826T (36.17%) and G7060A (61.70%) in L1; and G7521A (98.08%) in the LCR. The most prevalent amino acid variations were D25E in E6 and N29S in E7. In addition, 28 novel variations of HPV 16 were reported. Some covariations between different genes were obtained. In this study, HPV 16 variants belonged to the European lineage and the Asian lineage. Compared with neighboring districts, the distribution of HPV 16 variants in northeast China had a typical pattern. As the first report on HPV 16 variants in northeast China, it should be helpful for designing a HPV vaccine and HPV vaccination program in China.

Human papillomavirus 16 variants from Zambian women with normal pap smears

Human papillomavirus (HPV) type 16 is the most prevalent high-risk viral genotype associated with cervical cancer. Six distinct phylogenetic clusters of HPVs have been identified and are distributed differently across five continents. HPV16 DNA was extracted from cervicolavage samples from women with normal pap smears. The LCR regions were amplified in triplicate, cloned, sequenced, and analyzed from a total of 11 recovered HPV16 positive samples [Ng'andwe et al. (2007): BMC Infect Dis 7:77] were analyzed for sequence variation. The HPV16 LCR variants were assessed for promoter activity by use of a luciferase reporter gene. Six novel HPV16 variants with nucleotide exchanges in the LCR region were identified. Five clones were classified as European group HPV16 variants and one as an African group variant. Two of these variants had relatively lower promoter activity, 30% of that of the wild-type strain. The decreased promoter activity of some HPV16 variants may decrease expression of viral oncogenes and may be linked with the development, phenotype and severity of the cervical lesions in women infected with these across HPV16 variants.

HPV-related carcinogenesis: basic concepts, viral types and variants

Human papillomavirus (HPV) constitutes a diverse group of small DNA virus, some extensively studied during the last three decades due to their carcinogenic potential. Persistence of viral infections and uncontrolled expression of E6 and E7 viral oncogenes are critical events in transformation. A surprisingly large number of different HPV types have been identified and classified (>100) and it has been anticipated that almost 200 may exist. HPV types are thought to have originated very early during human evolution and are now defined by their L1 genomic sequence, differing by >10% among them. Importantly, viral types are cell-type specific and usually produce different kinds of lesions, benign or malignant. In addition, these types have co-evolved with their hosts and have generated what we call now intratype variants. Variants of HPV types are found associated with the ethnicity of the populations and have been grouped geographically. It is believed that HPV intratype variants may differ in biological behavior. Recognition of the crucial role that some specific HPV types play in cervical cancer development is highly important for their prevention and implementation of public health strategies to control cervical cancer, still the leading cause of death among cancer patients in many developing nations. Here we review basic concepts of HPV-induced carcinogenesis and molecular differences found among HPV types and intratype variants and discuss their clinical and functional implications.

Cellular factor YY1 downregulates the human papillomavirus 16 E6/E7 promoter, P97, in vivo and in vitro from a negative element overlapping the transcription-initiation site

Cellular factors that bind to cis sequences in the human papillomavirus 16 (HPV-16) upstream regulatory region (URR) positively and negatively regulate the viral E6 and E7 oncogene promoter, P97. DNase I footprinting has revealed the binding of cellular proteins to two previously undetected cis elements overlapping and 3′ of the transcription-initiation site of the P97 promoter. Mutations within homologous motifs found in both of these cis elements abolished their negative function in vivo and the binding of the same cellular complex in vitro. This factor was identified as YY1 by complex mobility and binding specificity in comparison with vaccinia virus-expressed, purified recombinant YY1 protein and by antigenic reactivity with YY1 antisera. Cis mutations in the ‘initiator’ YY1 site activated the P97 promoter in vivo and in vitro. P97 was also activated threefold in vitro by depletion of endogenous YY1 with wild-type, but not mutant, YY1 oligonucleotides from the IgH kappa E3′ enhancer. Furthermore, increasing concentrations of exogenous, purified recombinant YY1 repressed wild-type P97 transcript levels by up to threefold, but did not influence the P97 promoter mutated in the ‘initiator’ YY1 site. Thus, the promoter-proximal YY1 site was not necessary for correct transcription initiation at the P97 promoter, but was found to be required for downregulation of P97 transcription in vivo and in vitro. In contrast to other viral and cellular promoters, where YY1 is thought to function as a positive transcription-‘initiator’ factor, HPV-16 P97 transcription is downregulated by YY1 from a critical motif overlapping the transcription start site.

Human papillomavirus (HPV) type 18 induces extended growth in primary human cervical, tonsillar, or foreskin keratinocytes more effectively than other high-risk mucosal HPVs

Mucosal high-risk (HR) human papillomaviruses (HPVs) that cause cervical and other anogenital cancers also are found in approximately 25% of head and neck carcinomas (HNCs), especially those arising in the oropharynx and the tonsils. While many HR HPV types are common in anogenital cancer, over 90% of HPV-positive HNCs harbor HPV type 16 (HPV-16). Using a quantitative colony-forming assay, we compared the ability of full-length mucosal HPV genomes, i.e., the low-risk HPV-11 and HR HPV-16, -18, and -31, to persist in and alter the growth of primary human keratinocytes from the foreskin, cervix, and tonsils. The HR HPV types led to the formation of growing keratinocyte colonies in culture independent of the site of epithelial origin. However, HPV-18 induced colony growth in all keratinocytes >4-fold more effectively than HPV-16 or HPV-31 and >20-fold more efficiently than HPV-11 or controls. HPV-11-transfected or control colonies failed to expand beyond 32 to 36 population doublings postexplantation. In contrast, individual HR HPV-transfected clones exhibited no apparent slowdown of growth or "crisis," and many maintained HPV plasmid persistence beyond 60 population doublings. Keratinocyte clones harboring extrachromosomal HR HPV genomes had shorter population doubling times and formed dysplastic stratified epithelia in organotypic raft cultures, mirroring the pathological features of higher-grade intraepithelial lesions, yet did not exhibit chromosomal instability. We conclude that, in culture, the HR HPV type, rather than the site of epithelial origin of the cells, determines the efficacy of inducing continued growth of individual keratinocytes, with HPV-18 being the most aggressive mucosal HR HPV type tested.

Upstream regulatory region alterations found in human papillomavirus type 16 (HPV-16) isolates from cervical carcinomas increase transcription, ori function, and HPV immortalization capacity in culture

Human papillomavirus (HPV) DNAs isolated from cervical and head and neck carcinomas frequently contain nucleotide sequence alterations in the viral upstream regulatory region (URR). Our study has addressed the role such sequence changes may play in the efficiency of establishing HPV persistence and altered keratinocyte growth. Genomic mapping of integrated HPV type 16 (HPV-16) genomes from 32 cervical cancers revealed that the viral E6 and E7 oncogenes, as well as the L1 region/URR, were intact in all of them. The URR sequences from integrated and unintegrated viral DNA were found to harbor distinct sets of nucleotide substitutions. A subset of the altered URRs increased the potential of HPV-16 to establish persistent, cell growth-altering viral-genome replication in the cell. This aggressive phenotype in culture was not solely due to increased viral early gene transcription, but also to augmented initial amplification of the viral genome. As revealed in a novel ori-dependent HPV-16 plasmid amplification assay, the altered motifs that led to increased viral transcription from the intact genome also greatly augmented HPV-16 ori function. The nucleotide sequence changes correlate with those previously described in the distinct geographical North American type 1 and Asian-American variants that are associated with more aggressive disease in epidemiologic studies and encompass, but are not limited to, alterations in previously characterized sites for the negative regulatory protein YY1. Our results thus provide evidence that nucleotide alterations in HPV regulatory sequences could serve as potential prognostic markers of HPV-associated carcinogenesis.

HPV variants and HLA polymorphisms: the role of variability on the risk of cervical cancer

Human papillomavirus (HPV) infection is linked to the development of cervical cancer, and several cofactors contribute to the risk of disease. Research on the intratypic variability of HPVs has defined variants that are associated with persistent infections and are potentially more oncogenic, translating to a higher risk of malignant disease. The genetic variability of the host also plays a role in the risk of cervical cancer, especially genes controlling the immune response, such as HLA class I and II. These highly polymorphic genes are important risk determinants of HPV persistence and disease progression. The interaction between host and viral factors is complex and needs to be further investigated, paving the way to better define the patients at the highest risk of developing malignant diseases linked to HPV infection.

Functional mapping of the human papillomavirus type 16 E1 cistron

Replication of the double-stranded, circular human papillomavirus (HPV) genomes requires the viral DNA replicase E1. Here, we report an initial characterization of the E1 cistron of HPV type 16 (HPV-16), the most common oncogenic mucosal HPV type found in cervical and some head and neck cancers. The first step in HPV DNA replication is an initial burst of plasmid viral DNA amplification. Complementation assays between HPV-16 genomes carrying mutations in the early genes confirmed that the expression of E1 was necessary for initial HPV-16 plasmid synthesis. The major early HPV-16 promoter, P97, was dispensable for E1 production in the initial amplification because cis mutations inactivating P97 did not affect the trans complementation of E1- mutants. In contrast, E1 expression was abolished by cis mutations in the splice donor site at nucleotide (nt) 226, the splice acceptor site at nt 409, or a TATAA box at nt 7890. The mapping of 5' mRNA ends using rapid amplification of cDNA ends defined a promoter with a transcription start site at HPV-16 nt 14, P14. P14-initiated mRNA levels were low and required intact TATAA (7890). E1 expression required the HPV-16 keratinocyte-dependent enhancer, since cis mutations in its AP-2 and TEF-1 motifs abolished the ability of the mutant genomes to complement E1- genomes, and it was further modulated by origin-proximal and -distal binding sites for the viral E2 gene products. We conclude that P14-initiated E1 expression is critical for and limiting in the initial amplification of the HPV-16 genome.

Phylogeny and evolution of papillomaviruses based on the E1 and E2 proteins

Papillomaviridae are a family of small double-stranded DNA viruses that infect stratified squamous epithelia in vertebrates. Members of this family are causative agents of malignant tumours, such as cervical cancer while others are associated with benign proliferative lesions. So far, Papillomaviruses (PVs) are classified according to the sequence identity in the capsid gene L1. However, evidence has accumulated indicating a discontinuity in the evolutionary history of the L1 and L2 genes of many PVs, giving rise to differences in the phylogenetic reconstructions of the early and of the late genes. Neither the oncogenes E5, E6 and E7 nor the upstream regulatory region are suitable for phylogenetic inference due to the poor conservation along the Papillomaviridae family. We have analysed here the evolutionary relationships of the PVs with respect to the E1 and E2 proteins, and the results provide both phylogeny and biologic behaviour of the viruses. The hierarchical taxonomic relationships can be structured as an alternative classification system in which mucosal high-risk viruses, mucosal low-risk viruses and viruses associated with cutaneous lesions are grouped separately and do not appear intermingled. Some important trends are also observed: first, evolution of the PVs has not been homogeneous, even in viruses that infect the same host, and second mucosal human PVs have evolved faster than their cutaneous counterparts. The evolutionary analysis based on the E1 and E2 proteins will allow us to better understand the generation of the diversity of the PVs and the development of malignancy associated with these viruses.

High grade cervical lesions are caused preferentially by non-European variants of HPVs 16 and 18

The intratypic variability of HPVs 16 and 18 has been extensively studied and has been used as an important tool in epidemiological studies of viral transmission, persistence and progression to clinically relevant cervical lesions. Infections by non-European variants of HPVs 16 and 18 are associated with an increased risk for the development of high grade squamous intraepithelial lesions (HSIL). Our aim was to correlate the intratypic molecular variability of both HPV types and risk of persistent infection and lesion outcome in a cohort study conducted in Brazil. We characterized molecular variants of HPV types 16 and 18 by sequencing a fragment of the LCR, and of the E6 and L1 genes, for HPV-16 variants only. For both types, European variants composed the most prevalent and diverse group. Persistent infections with HPV-18 were associated with continuous detection of European variants. However, risk for simultaneous detection of HSIL and HPV DNA was higher in women harboring non-European variants of HPV-16. The same trend was observed with HSIL detected during follow-up. Our study confirms the association between non-European variants and risk of cervical neoplasia, and highlights the importance of their geographic distribution for cervical cancer risk assessment.

Epidemiological and functional implications of molecular variants of human papillomavirus

Human papillomavirus genomes are classified into molecular variants when they present more than 98% of similarity to the prototype sequence within the L1 gene. Comparative nucleotide sequence analyses of these viruses have elucidated some features of their phylogenetic relationship. In addition, human papillomavirus intratype variability has also been used as an important tool in epidemiological studies of viral transmission, persistence and progression to clinically relevant cervical lesions. Until the present, little has been published concerning the functional significance of molecular variants. It has been shown that nucleotide variability within the long control region leads to differences in the binding affinity of some cellular transcriptional factors and to the enhancement of the expression of E6 and E7 oncogenes. Furthermore, in vivo and in vitro studies revealed differences in E6 and E7 biochemical and biological properties among molecular variants. Nevertheless, further correlation with additional functional information is needed to evaluate the significance of genome intratypic variability. These results are also important for the development of vaccines and to determine the extent to which immunization with L1 virus-like particles of one variant could induce antibodies that cross-neutralize other variants.

Different papillomaviruses have different repertoires of transcription factor binding sites: convergence and divergence in the upstream regulatory region

Background

Papillomaviruses (PVs) infect stratified squamous epithelia in warm-blooded vertebrates and have undergone a complex evolutionary process. The control of the expression of the early ORFs in PVs depends on the binding of cellular and viral transcription factors to the upstream regulatory region (URR) of the virus. It is believed that there is a core of transcription factor binding sites (TFBS) common to all PVs, with additional individual differences, although most of the available information focuses only on a handful of viruses.

Results

We have studied the URR of sixty-one PVs, covering twenty different hosts. We have predicted the TFBS present in the URR and analysed these results by principal component analysis and genetic algorithms. The number and nature of TFBS in the URR might be much broader than thus far described, and different PVs have different repertoires of TFBS.

Conclusion

There are common fingerprints in the URR in PVs that infect primates, although the ancestors of these viruses diverged a long time ago. Additionally, there are obvious differences between the URR of alpha and beta PVs, despite these PVs infect similar histological cell types in the same host, i.e. human. A thorough analysis of the TFBS in the URR might provide crucial information about the differential biology of cancer-associated PVs.

Genome variation of human papillomavirus types: phylogenetic and medical implications

Human papillomaviruses (HPVs) are described as "types" based on their genome sequences and identified by a number. For example, HPV-6 is associated with genital warts, and HPV-16 with anogenital cancers. The genomes of many HPV types have been reisolated, sequenced and compared to reference "prototypes" countless times by laboratories throughout the world. It was found that each HPV type occurs in the form of "variants", identified by about 2% nucleotide differences in most genes and 5% in less conserved regions. Less than 100 variants of any HPV type have been detected, a scenario that is very different from the quasi-species formed by many RNA viruses. The variants of each HPV type form phylogenetic trees, and variants from specific branches are often unique to specific ethnic groups. Immigrant populations contain, depending on their respective ethnic origins, mixtures of variants. The absence of HPV genomes intermediate to specific types show that all HPV types existed already when humans became a species. Consequently, humans had always suffered from lesions like anogenital cancer, genital warts and common warts. A growing number of epidemiological, etiological and molecular data suggest that variants of the same HPV type are biologically distinct and may confer differential pathogenic risks. Since the distribution of some variants of HPV-16 and 18 correlates with the distribution of human populations that have an increased risk to develop anogenital cancer, the study of HPV type variation may point to one of the reasons for the higher incidence rates of these lesions in specific cohorts.

Worldwide genomic diversity of the high-risk human papillomavirus types 31, 35, 52, and 58, four close relatives of human papillomavirus type 16

Among the more than one hundred formally described human papillomavirus (HPV) types, 18 are referred to as high-risk HPV types due to their association with anogenital cancer. Despite pathogenic similarities, these types form three remotely related taxonomic groups. One of these groups is called HPV species 9 and is formed by HPV-16, the most common and best-studied type, together with HPV-31, -33, -35, -52, -58, and -67. Previous worldwide comparisons of HPV-16 samples showed about 2% nucleotide diversity between isolates, which were subsequently termed variants. The distribution of divergent variants has been found to correlate frequently with the geographic origin and the ethnicity of the infected patients and led to the concept of unique African, European, Asian, and Native American HPV-16 variants. In the current study, we address the question of whether geography and ethnicity also correlate with sequence variations found for HPV-31, -35, -52, and -58. This was done by sequencing the long control region in samples derived from Europe, Asia, and Africa, and from immigrant populations in North and South America. We observed maximal divergence between any two variants within each of these four HPV types ranging from 1.8 to 3.6% based on nucleotide exchanges and, occasionally, on insertions and deletions. Similar to the case with HPV-16, these mutations are not random but indicate a relationship between the variants in form of phylogenetic trees. An interesting example is presented by a 16-bp insert in select variants of HPV-35, which appears to have given rise to additional variants by nucleotide exchanges within the insert. All trees showed distinct phylogenetic topologies, ranging from dichotomic branching in the case of HPV-31 to star phylogenies of the other three types. No clear similarities between these types or between these types and HPV-16 exist. While variant branches in some types were specific for Europe, Africa, or East Asia, none of the four trees reflected human evolution and spread to the extent illustrated by HPV-16. One possible explanation is that the rare HPV types that we studied spread and thereby diversified more slowly than the more abundant HPV-16 and may have established much of today's variant diversity already before the worldwide spread of humans 100,000 years ago. Most variants had prototypic amino acid sequences within the E6 oncoprotein and a segment of the L1 capsid protein. Some had one, two, or three amino acid substitutions in these regions, which might indicate biological and pathogenic diversity between the variants of each HPV type.

Human papillomavirus type 16 E1circumflexE4 contributes to multiple facets of the papillomavirus life cycle

The life cycle of human papillomaviruses (HPVs) is tightly linked to the differentiation program of the host's stratified epithelia that it infects. E1(circumflex)E4 is a viral protein that has been ascribed multiple biochemical properties of potential biological relevance to the viral life cycle. To identify the role(s) of the viral E1(circumflex)E4 protein in the HPV life cycle, we characterized the properties of HPV type 16 (HPV16) genomes harboring mutations in the E4 gene in NIKS cells, a spontaneously immortalized keratinocyte cell line that when grown in organotypic raft cultures supports the HPV life cycle. We learned that E1(circumflex)E4 contributes to the replication of the viral plasmid genome as a nuclear plasmid in basal cells, in which we also found E1(circumflex)E4 protein to be expressed at low levels. In the suprabasal compartment of organotypic raft cultures harboring E1(circumflex)E4 mutant HPV16 genomes there were alterations in the frequency of suprabasal cells supporting DNA synthesis, the levels of viral DNA amplification, and the degree to which the virus perturbs differentiation. Interestingly, the comparison of the phenotypes of various mutations in E4 indicated that the E1(circumflex)E4 protein-encoding requirements for these various processes differed. These data support the hypothesis that E1(circumflex)E4 is a multifunctional protein and that the different properties of E1(circumflex)E4 contribute to different processes in both the early and late stages of the virus life cycle.