Evolution and taxonomic classification of human papillomavirus 16 (HPV16)-related variant genomes: HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67

Background

Human papillomavirus 16 (HPV16) species group (alpha-9) of the Alphapapillomavirus genus contains HPV16, HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67. These HPVs account for 75% of invasive cervical cancers worldwide. Viral variants of these HPVs differ in evolutionary history and pathogenicity. Moreover, a comprehensive nomenclature system for HPV variants is lacking, limiting comparisons between studies.

Methods

DNA from cervical samples previously characterized for HPV type were obtained from multiple geographic regions to screen for novel variants. The complete 8 kb genomes of 120 variants representing the major and minor lineages of the HPV16-related alpha-9 HPV types were sequenced to capture maximum viral heterogeneity. Viral evolution was characterized by constructing phylogenic trees based on complete genomes using multiple algorithms. Maximal and viral region specific divergence was calculated by global and pairwise alignments. Variant lineages were classified and named using an alphanumeric system; the prototype genome was assigned to the A lineage for all types.

Results

The range of genome-genome sequence heterogeneity varied from 0.6% for HPV35 to 2.2% for HPV52 and included 1.4% for HPV31, 1.1% for HPV33, 1.7% for HPV58 and 1.1% for HPV67. Nucleotide differences of approximately 1.0% - 10.0% and 0.5%–1.0% of the complete genomes were used to define variant lineages and sublineages, respectively. Each gene/region differs in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) /noncoding region 2 (NCR2) >upstream regulatory region (URR)> E6/E7 > E2/L2 > E1/L1.

Conclusions

These data define maximum viral genomic heterogeneity of HPV16-related alpha-9 HPV variants. The proposed nomenclature system facilitates the comparison of variants across epidemiological studies. Sequence diversity and phylogenies of this clinically important group of HPVs provides the basis for further studies of discrete viral evolution, epidemiology, pathogenesis and preventative/therapeutic interventions.     

Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR

Naturally occurring genetic variants of human papillomavirus type 16 (HPV16) are common and have previously been classified into 4 major lineages; European-Asian (EAS), including the sublineages European (EUR) and Asian (As), African 1 (AFR1), African 2 (AFR2), and North-American/Asian-American (NA/AA). We aimed to improve the classification of HPV16 variant lineages by using a large resource of HPV16-positive cervical samples collected from geographically diverse populations in studies on HPV and/or cervical cancer undertaken by the International Agency for Research on Cancer. In total, we sequenced the entire E6 genes and long control regions (LCRs) of 953 HPV16 isolates from 27 different countries worldwide. Phylogenetic analyses confirmed previously described variant lineages and subclassifications. We characterized two new sublineages within each of the lineages AFR1 and AFR2 that are robustly classified using E6 and/or the LCR. We could differentiate previously identified AA1, AA2, and NA sublineages, although they could not be distinguished by E6 alone, requiring the LCR for correct phylogenetic classification. We thus provide a classification system for HPV16 genomes based on 13 and 32 phylogenetically distinguishing positions in E6 and the LCR, respectively, that distinguish nine HPV16 variant sublineages (EUR, As, AFR1a, AFR1b, AFR2a, AFR2b, NA, AA1, and AA2). Ninety-seven percent of all 953 samples fitted this classification perfectly. Other positions were frequently polymorphic within one or more lineages but did not define phylogenetic subgroups. Such a standardized classification of HPV16 variants is important for future epidemiological and biological studies of the carcinogenic potential of HPV16 variant lineages.     

Apoptosis induced by an antagonist peptide against HPV16 E7 in vitro and in vivo via restoration of p53

Human papilloma virus type 16 (HPV16) E7 is a viral oncoprotein that is believed to play a major role in cervical neoplasia. A novel antagonist peptide against HPV16 E7 was previously selected by phage display screening and the selected peptide was found to have anti-tumor efficacy against HPV16-positive cervical carcinoma through induction of cell cycle arrest. In the current study, to further elucidate the mechanisms of the antagonist peptide, the effects of the peptide on apoptosis are investigated by RT-PCR, Western blotting, MTT assay, TUNEL staining, Annexin V apoptosis assay, flow cytometry, and animal experiments. The antagonist peptide showed obvious anti-tumor efficacy through apoptosis induction, both in HPV16-positive cervical cancer cell lines and tumor xenografts. Our results also revealed that the peptide induced accumulation of cellular p53 and p21, and led to HPV16 E7 protein degradation. In the case of mRNA levels, it resulted in unaltered p53 and HPV16 E7 expression, but increased expression of p21. In contrast, the induction of apoptosis and p53 reactivation effects by the selected peptide were abolished after E7 knocked down with siRNA. These results demonstrate that the selected peptide can induce E7 degradation and lead to marked apoptosis in HPV16-related cancer cells by activating cellular p53 and its target genes, such as p21. Furthermore, the evident therapeutic efficacy obtained from the subcutaneous tumor model experiments in nude mice suggests a therapeutic potential for HPV16-related cancers of the selected peptide. Therefore, this specific peptide may be used to create specific biotherapies for the treatment of HPV 16-positive cervical cancers.     

Oncogenic comparison of human papillomavirus type 58 E7 variants

Human papillomavirus 58 (HPV58) ranks the second or third in East Asian cervical cancers. Current studies on HPV58 are scarce and focus on the prototype. Previously, we identified the three most common circulating HPV58 E7 strains contained amino acid alterations: G41R/G63D (51%), T20I/G63S (22%) and T74A/D76E (14%) respectively. Among them, the T20I/G63S variant (V1) had a stronger epidemiological association with cervical cancer. We therefore suggested that V1 possessed stronger oncogenicity than the other two variants. Here, we performed phenotypic assays to characterize and compare their oncogenicities with HPV58 E7 prototype. Our results showed that overexpression of V1 conferred a higher colony‐forming ability to primary murine epithelial cells than prototype (P < 0.05) and other variants, implicating its higher immortalising potential. Further experiments showed that both V1 and prototype enhanced the anchorage‐independent growth of NIH/3T3 cells (P < 0.001), implicating their stronger transforming power than the two other variants. Moreover, they possessed an increased ability to degrade pRb (P < 0.001), which is a major effector pathway of E7‐driven oncogenesis. Our work represents the first study to compare the oncogenicities of HPV58 E7 prototype and variants. These findings deepened our understanding of HPV58 and might inform clinical screening and follow‐up strategy.     

Human Papillomavirus 16 Non-European Variants Are Preferentially Associated with High-Grade Cervical Lesions

HPV16 accounts for 50–70% of cervical cancer cases worldwide. Characterization of HPV16 variants previously indicated that they differ in risks for viral persistence, progression to cervical precancer and malignant cancer. The aim of this study was to examine the association of severity of disease with HPV16 variants identified in specimens (n = 281) obtained from a Cervical Pathology and Colposcopy outpatient clinic in the University Hospital of Espírito Santo State, Southeastern Brazil, from April 2010 to November 2011. All cytologic and histologic diagnoses were determined prior to definitive treatment. The DNA was isolated using QIAamp DNA Mini Kit and HPV was detected by amplification with PGMY09/11 primers and positive samples were genotyped by RFLP analyses and reverse line blot. The genomes of the HPV16 positive samples were sequenced, from which variant lineages were determined. Chi² statistics was performed to test the association of HPV16 variants between case and control groups. The prevalence of HR-HPV types in <CIN1, CIN2 and CIN3+ were 33.7%, 84.4% and 91.6%, respectively. Thirty-eight of 49 (78%) HPV16 positive samples yielded HPV16 sequence information; of which, 32 complete genomes were sequenced and an additional 6 samples were partially sequenced. Phylogenetic analysis and patterns of variations identified 65.8% (n = 25) as HPV16 European (E) and 34.2% (n = 13) as non-European (NE) variants. Classification of disease into CIN3+ vs. <CIN3 indicated that NE types were associated with high-grade disease with an OR = 4.6 (1.07–20.2, p = 0.05). The association of HPV16 NE variants with an increased risk of CIN3+ is consistent with an HPV16 genetically determined enhanced oncogenicity. The prevalence of genetic variants of HPV16 is distributed across different geographical areas and with recent population admixture, only empiric data will provide information on the highest risk HPV16 variants within a given population.     

Genetic Variation of Human Papillomavirus Type 16 in Individual Clinical Specimens Revealed by Deep Sequencing

Viral genetic diversity within infected cells or tissues, called viral quasispecies, has been mostly studied for RNA viruses, but has also been described among DNA viruses, including human papillomavirus type 16 (HPV16) present in cervical precancerous lesions. However, the extent of HPV genetic variation in cervical specimens, and its involvement in HPV-induced carcinogenesis, remains unclear. Here, we employ deep sequencing to comprehensively analyze genetic variation in the HPV16 genome isolated from individual clinical specimens. Through overlapping full-circle PCR, approximately 8-kb DNA fragments covering the whole HPV16 genome were amplified from HPV16-positive cervical exfoliated cells collected from patients with either low-grade squamous intraepithelial lesion (LSIL) or invasive cervical cancer (ICC). Deep sequencing of the amplified HPV16 DNA enabled de novo assembly of the full-length HPV16 genome sequence for each of 7 specimens (5 LSIL and 2 ICC samples). Subsequent alignment of read sequences to the assembled HPV16 sequence revealed that 2 LSILs and 1 ICC contained nucleotide variations within E6, E1 and the non-coding region between E5 and L2 with mutation frequencies of 0.60% to 5.42%. In transient replication assays, a novel E1 mutant found in ICC, E1 Q381E, showed reduced ability to support HPV16 origin-dependent replication. In addition, partially deleted E2 genes were detected in 1 LSIL sample in a mixed state with the intact E2 gene. Thus, the methods used in this study provide a fundamental framework for investigating the influence of HPV somatic genetic variation on cervical carcinogenesis.     

Properties of p53 mutations detected in primary and secondary cervical cancers suggest mechanisms of metastasis and involvement of environmental carcinogens.

Primary human papillomavirus (HPV) positive anogenital cancers normally develop without somatic mutation within the p53 gene. In this study, however, we have identified p53 point mutations in metastases arising from HPV positive cervical carcinomas, suggesting that acquisition of p53 mutation may play a role in the progression of some HPV associated primary cancers. p53 mutants identified in anogenital cancers exhibit a dominant transforming phenotype and increased resistance to HPV16 E6 directed degradation. The association of p53 mutation with metastases may explain the poor prognosis reported for HPV negative primary cancers, many of which already contain mutant p53. A high proportion of p53 mutations detected in both primary and metastatic cancers are GC-->TA transversions, strongly suggesting a role for external carcinogens in the development of these cancers.     

Tumourigenesis Driven by the Human Papillomavirus Type 16 Asian-American E6 Variant in a Three-Dimensional Keratinocyte Model

Infection with a transforming human papillomavirus (HPV) such as type 16 (of species Alphapapillomavirus 9) causes ano-genital and oral tumours via viral persistence in human squamous cell epithelia. Epidemiological studies showed that the naturally occurring HPV16 Asian-American (AA) variant (sublineage D2/D3) is found more often than the European Prototype (EP) (sublineage A1) in high-grade cervical neoplasia and tumours compared to non-cancer controls. Just three amino acid changes within the early gene, E6, of HPV16 AA have been linked to this augmented tumourigenicity. The AAE6 variant''s greater immortalizing and transforming potential over EPE6 has recently been confirmed in retrovirally-transduced keratinocytes expressing the E6 gene only. However, the tumourigenic role of the full-length viral genome of HPV16 has not yet been addressed with regard to these E6 variants. To investigate this process in the context of these two HPV16 E6 genotypes, an organotypic tissue culture model was used to simulate the HPV infectious life cycle. The AAE6 variant demonstrated an enhanced ability over EPE6 to drive the viral life cycle toward tumourigenesis, as evidenced phenotypically—by a more severe grade of epithelial dysplasia with higher proliferation and deregulated differentiation, and molecularly—by high viral oncogene E6 and E7 expression, but lack of productive viral life cycle markers. In contrast, EPE6 had low E6 and E7 but high E1∧E4 expression, indicative of a productive life cycle. We suggest increased viral integration into the host genome for AAE6 as one possible mechanism for these observed differences from EPE6. Additionally, we found downstream effects on immortalization and host innate immune evasion. This study highlights how minor genomic variations in transforming viruses can have a significant affect on their tumourigenic ability.     

Genetic diversity and phylogenetic analysis of HPV 16 & 18 variants isolated from cervical specimens of women in Saudi Arabia

Human papillomavirus (HPV) are well known to be associated with the development of cervical cancer. HPV16 and HPV 18 are known as high-risk types and reported to be predominantly associated with cervical cancer. The prevalence and genetic diversity of HPV have been well documented globally but, in the Kingdom of Saudi Arabia, data on HPV genetic diversity are lacking. In this study, we have analyzed the genetic diversity of both HPV16 and HPV18 based on their L1 gene sequence because L1 gene is a major capsid protein gene and has been utilized to develop a prophylactic vaccine. In January 2011–2012, a total of forty samples from cervical specimens of women in Saudi Arabia were collected. The association of HPV16, HPV18 was detected by polymerase chain reaction, sequenced and submitted to GenBank. The sequences identity matrix and the phylogenetic relationship were analyzed with selected HPVs. The highest sequence identity (99.5%) for HPV16 and (99.3%) for HPV was observed with selected HPVs. The phylogenetic analysis results showed that HPVs from Saudi Arabia formed a closed cluster with African, Asian, East Asian as well as American HPVs distributed into multiple linages from various geographical locations. The results provided the valuable information about genetic diversity, but there is an urgent need to generate full genome sequence information which will provide a clearer picture of the genetic diversity and evolution of HPVs in Saudi Arabia. In conclusion, the generated data will be highly beneficial for developing molecular diagnostic tools, analyzing and correlating the epidemiological data to determine the risk of cervical cancer and finally to develop a vaccine for Saudi Arabian population.     

Prevalence of human papillomavirus type 16 variants in the Federal District, Central Brazil

We report the prevalence of human papillomavirus type 16 (HPV-16) variants in women with cervical lesions from the Federal District, Central Brazil. We analyzed 34 HPV-16 samples, identifying the sequence variations of E6 and L1 genes and correlating variant frequency with disease status. The most prevalent HPV-16 variant was the European (50%), followed by Asian-American (41.2%), African-1 (5.9%), and African-2 (2.9%). European and non-European variants appeared in equal frequencies among the cytological types of lesions - atypical squamous or glandular cells of undetermined significance, cytological alterations suggesting HPV infection, cervical intraepithelial neoplasias, squamous cell carcinoma, and adenocarcinoma.     

Human Papillomavirus Type 6 and 11 Genetic Variants Found in 71 Oral and Anogenital Epithelial Samples from Australia

Genetic variation of 49 human papillomavirus (HPV) 6 and 22 HPV11 isolates from recurrent respiratory papillomatosis (RRP) (n = 17), genital warts (n = 43), anal cancer (n = 6) and cervical neoplasia cells (n = 5), was determined by sequencing the long control region (LCR) and the E6 and E7 genes. Comparative analysis of genetic variability was examined to determine whether different disease states resulting from HPV6 or HPV11 infection cluster into distinct variant groups. Sequence variation analysis of HPV6 revealed that isolates cluster into variants within previously described HPV6 lineages, with the majority (65%) clustering to HPV6 sublineage B1 across the three genomic regions examined. Overall 72 HPV6 and 25 HPV11 single nucleotide variations, insertions and deletions were observed within samples examined. In addition, missense alterations were observed in the E6/E7 genes for 6 HPV6 and 5 HPV11 variants. No nucleotide variations were identified in any isolates at the four E2 binding sites for HPV6 or HPV11, nor were any isolates found to be identical to the HPV6 lineage A or HPV11 sublineage A1 reference genomes. Overall, a high degree of sequence conservation was observed between isolates across each of the regions investigated for both HPV6 and HPV11. Genetic variants identified a slight association with HPV6 and anogenital lesions (p = 0.04). This study provides important information on the genetic diversity of circulating HPV 6 and HPV11 variants within the Australian population and supports the observation that the majority of HPV6 isolates cluster to the HPV6 sublineage B1 with anogenital lesions demonstrating an association with this sublineage (p = 0.02). Comparative analysis of Australian isolates for both HPV6 and HPV11 to those from other geographical regions based on the LCR revealed a high degree of sequence similarity throughout the world, confirming previous observations that there are no geographically specific variants for these HPV types.     

Human papillomavirus type 16 E7 oncoprotein associates with the cullin 2 ubiquitin ligase complex, which contributes to degradation of the retinoblastoma tumor suppressor

Human papillomavirus type 16 (HPV16) and other high-risk HPVs are etiologically linked to the development of cervical carcinomas and contribute to a number of other tumors of the anogenital tract, as well as oral cancers. The high-risk HPV E6 and E7 oncoproteins are consistently expressed in cervical cancer cells and are necessary for the induction and maintenance of the transformed phenotype. An important aspect of HPV16 E7's oncogenic activities is destabilization of the retinoblastoma tumor suppressor (pRB) through a ubiquitin/proteasome-dependent mechanism, although the exact molecular mechanism is unknown. Here, we report that HPV16 E7 is associated with an enzymatically active cullin 2 ubiquitin ligase complex and that the HPV16 E7/pRB complex contains cullin 2. Depletion of cullin 2 by RNA interference causes increased steady-state levels and stability of pRB in HPV16 E7-expressing cells, and ectopic expression of HPV16 E7 and the cullin 2 complex leads to pRB ubiquitination in vivo. Hence, we propose that the HPV16 E7-associated cullin 2 ubiquitin ligase complex contributes to aberrant degradation of the pRB tumor suppressor in HPV16 E7-expressing cells.     

水稻一多拷贝微卫星DNA多态性分析

用一个多拷贝微卫星 DNA标记分析了238份栽培水稻的遗传多样性及遗传多样性从农家品种到现代栽培品种的动态变化。共检测出16种长度变异类型和32种表现型,但长度变异类型数、表现型数和表现型多样性水平在农家品种和现代栽培品种间没有显著差别。该标记在遗传资源评价和DNA指纹图谱建设中有重要的应用价值,也是研究环境或人工对多拷贝基因选择作用的一个良好模型。 Abstract:The genetic diversity of 238 cultivated rice and its difference between landraces and modern cultivars have been assayed using a multiple copy microsatellite DNA marker.Although 16 length variants and 32 phenotypes were found in all accessions,no significant differences in number of length variants,number of phenotypes and phenotypic diversity were detected between landraces and modern cultivars.This marker is of great value in evaluating genetic germplasm and constructing DNA fingerprints of rice and provides an elite model for tracing effectiveness of natural and artificial selection on multiple copy genes.     

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

While the central role of the viral upstream regulatory region (URR) in the human papillomavirus (HPV) life cycle has been well established, its effects on viral replication factor expression and plasmid replication of HPV type 16 (HPV16) remain unclear. Some nonprototypic variants of HPV16 contain altered URR sequences and are considered to increase the oncogenic risk of infections. To determine the relationship between viral replication and variant URRs, hybrid viral genomes were constructed with the replication-competent HPV16 prototype W12 and analyzed in assays which recapitulate the different phases of normal viral replication. The establishment efficiencies of hybrid HPV16 genomes differed about 20-fold among European prototypes and variants from Africa and America. Generally, European and African genomes exhibited the lowest replication efficiencies. The high replication levels observed with American variants were primarily attributable to their efficient expression of the replication factors E1 and E2. The maintenance levels of these viral genomes varied about fivefold, which correlated with their respective establishment phenotypes and published P(97) activities. Vegetative DNA amplification could also be observed with replicating HPV16 genomes. These results indicate that efficient E1/E2 expression and elevated plasmid replication levels during the persistent stage of infection may comprise a risk factor in HPV16-mediated oncogenesis.     

From loss to gain: role for SUN1 in laminopathies

Background

About 15 types of human papillomavirus (HPV) are classified as high-risk based on their epidemiological link with cervical cancer. These HPV types have deferent degrees of oncogenicity and their distribution among cervical precancers and cancers varies ethnogeographically. HPV58 is rare worldwide but being found more commonly in East Asia.

Findings

A high prevalence of HPV58 among squamous cell carcinoma has been reported from China (28% in Shanghai, 10% in Hong Kong and 10% in Taiwan) and other countries in East Asia including Korea (16%) and Japan (8%). HPV58 ranks the third in Asia overall, but contributes to only 3.3% of cervical cancers globally. The reasons for a difference in disease attribution may lie on the host as well as the virus itself. HLA-DQB1*06 was found to associate with a higher risk of developing HPV58-positive cervical neoplasia in Hong Kong women, but not neoplasia caused by other HPV types. An HPV58 variant (E7 T20I, G63S) commonly detected in Hong Kong was found to confer a 6.9-fold higher risk of developing cervical cancer compared to other variants. A study involving 15 countries/cities has shown a predilection in the distribution of HPV58 variant lineages. Sublineage A1, the prototype derived from a cancer patient in Japan, was rare worldwide except in Asia.

Conclusions

HPV58 accounts for a larger share of disease burden in East Asia, which may be a result of differences in host genetics as well as the oncogenicity of circulating variants. These unique characteristics of HPV58 should be considered in the development of next generation vaccines and diagnostic assays.     

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.     

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.