Nucleotide polymorphisms of the human papillomavirus 16 E1 gene

Whole-genome analysis and evolutionary characterization of cervical and oral human papillomavirus 16

High-throughput genome-wide sequencing has revealed high genomic variability of HPV16 in different geographic regions which is the most predominant genotype in human papillomavirus (HPV)-associated malignancies. Analysis of the HPV16 by whole-genome sequence (WGS) is an advanced method for the identification of mutations in the genome. There is limited information about HPV16 diversity in Pakistan, especially at the genomic level. Till now, WGS for HPV16 has not been previously reported in Pakistan. The current study has sequenced three HPV16 viral genomes, from two cervical and one oral cavity positive sample of women presented with general gynecological problems without any evidence of precancerous or cancerous lesions using an ion ampliseq customized panel. Sequencing analysis detected 38 variations, including single-nucleotide polymorphisms (SNPs) and two Indels, across three samples with the highest number of SNPs present in E1, E2, and L2, respectively. A total of 20 non-synonymous and 11 synonymous mutations with amino acid substitutions (T1421C, G1515A, T2223C, T1389C, G1483A, and T2191C) were identified. The phylogenetic analysis revealed the genomes of HPV16 are closely associated with those reported from Thailand and the United States. These are the first HPV16 WGS from Pakistan. However, more research is needed with a large sample size from diversified areas to assess the carcinogenic consequences and impact of HPV vaccinations.

Mutation Profile of HPV16 L1 and L2 Genes in Different Geographic Areas

The causal relationship between HPV and cervical cancer in association with the high prevalence of high risk HPV genotypes led to the design of HPV vaccines based on the major capsid L1 protein. In recent years, capsid protein L2 has also become a focal point in the field of vaccine research. The present review focuses on the variability of HPV16 L1 and L2 genes, emphasizing the distribution of specific amino acid changes in the epitopes of capsid proteins. Moreover, a substantial bioinformatics analysis was conducted to describe the worldwide distribution of amino acid substitutions throughout HPV16 L1, L2 proteins. Five amino acid changes (T176N, N181T; EF loop), (T266A; FG loop), (T353P, T389S; HI loop) are frequently observed in the L1 hypervariable surface loops, while two amino acid substitutions (D43E, S122P) are adjacent to L2 specific epitopes. These changes have a high prevalence in certain geographic regions. The present review suggests that the extensive analysis of the amino acid substitutions in the HPV16 L1 immunodominant loops may provide insights concerning the ability of the virus in evading host immune response in certain populations. The genetic variability of the HPV16 L1 and L2 epitopes should be extensively analyzed in a given population.

Genetic variability of the HPV16 early genes and LCR. Present and future perspectives

Human papillomavirus 16 (HPV16) infection is the aetiologic factor for the development of cervical dysplasia and is regarded as highly carcinogen, because it is implicated in more than 50% of cervical cancer cases, worldwide. The tumourigenic potential of HPV16 has triggered the extensive sequence analysis of viral genome in order to identify nucleotide variations and amino acid substitutions that influence viral oncogenicity and subsequently the initiation and progression of cervical cancer. Nowadays, specific mutations of HPV16 DNA have been associated with an increased risk of high-grade squamous intraepithelial lesions and invasive cervical cancer (ICC) development, including E6: Q14H, H78Y, L83V, Ε7: N29S, S63F, E2: H35Q, P219S, T310K, E5: I65V, whereas highly conserved regions of viral DNA have been extensively characterised. In addition, numerous novel HPV16 mutations are observed among the studied populations from various geographic regions, hence advocating that different HPV16 strains seem to emerge with different tumourigenic capacities. The present review focuses on the variability of the early genes and the long control region, emphasising on the association of specific mutations with the development of severe dysplasia. Finally, it evaluates whether specific regions of HPV16 DNA are able to serve as valuable biomarkers for cervical cancer risk.

Epidemiology and Burden of Human Papillomavirus and Related Diseases, Molecular Pathogenesis, and Vaccine Evaluation

Diagnosed in more than 90% of cervical cancers, the fourth deadliest cancer in women, human papillomavirus (HPV) is currently the most common pathogen responsible for female cancers. Moreover, HPV infection is associated with many other diseases, including cutaneous and anogenital warts, and genital and upper aerodigestive tract cancers. The incidence and prevalence of these pathologies vary considerably depending on factors including HPV genotype, regional conditions, the study population, and the anatomical site sampled. Recently, features of the cervicovaginal microbiota are found to be associated with the incidence of HPV-related diseases, presenting a novel approach to identify high-risk women through both blood and cervical samples. Overall, the HPV repartition data show that HPV infection and related diseases are more prevalent in developing countries. Moreover, the available (2-, 4-, and 9-valent) vaccines based on virus-like particles, despite their proven effectiveness and safety, present some limitations in terms of system development cost, transport cold chain, and oncogenic HPV variants. In addition, vaccination programs face some challenges, leading to a considerable burden of HPV infection and related diseases. Therefore, even though the new (9-valent) vaccine seems promising, next-generation vaccines as well as awareness programs associated with HPV vaccination and budget reinforcements for immunization are needed.

Human Papillomavirus Type 16 E1 Mutations Associated with Cervical Cancer in a Han Chinese Population

Human papillomavirus type 16 (HPV16) is a high-risk HPV type and a potent carcinogen. HPV E1 is one of the most highly conserved proteins and it plays a central role in initiating HPV DNA replication. In current study, we enrolled 161 HPV16-positive cervical cancer patients (case group) and 171 HPV16-positive asymptomatic individuals (control group) in a study to analyse the association between HPV16 E1 genetic mutations and cervical cancer. The samples of case group were cervical cancer tissues and the samples of control group were cervical exfoliated cells. Three variants (A4, A1-A3 and D3) were found in the case group, 68.3% of the HPV16 E1 sequences belonged to the A4 (As) sub-lineage, 29.2% belonged to the A1-A3 (EUR) sub-lineage, and 2.5% belonged to the D3 (AA1) sub-lineage. Two variants (A4 and A1-A3) occurred in the control group. The A4 (As) sub-lineage was predominant in this group as well (66.1%), followed by the A1-A3 (EUR) sub-lineage (33.9%), but the D3 (AA1) sub-lineage was not found in the control group. The distribution of the HPV16 variants between the case and control groups was significantly different (P<0.05). When the distribution of the HPV16 E1 gene mutations was compared, the distribution of twenty-seven mutations was significantly different between the case and control groups (P<0.05), and twenty-two mutations occurred only in the D3 (AA1) sub-lineage, two were found only in the A4 (As) sub-lineage, one was found in the A1-A3 (EUR) sub-lineage, two was found in both the A4 (As) and A1-A3 (EUR) sub-lineages. In the sub-lineage analysis, the differences in the T933A (A23A), T1014G (D50E) and G2160A (R432R) mutations were statistically significant between the case and control groups for the A4 (As) sub-lineage (P<0.05), and the differences in the T2232C (F456F), G2337A (M491I) and A2547G (P561P) mutations were statistically significant between the case and control groups for the A1-A3 (EUR) sub-lineage (P<0.05). In the current study, we describe specific mutations in the HPV16 E1 gene associated with cervical cancer, and our study will provide a good reference for further functional studies of the relationship between cervical cancer carcinogenesis and HPV genes.

HPV16 viral characteristics in primary, recurrent and metastatic vulvar carcinoma

Vulvar carcinoma is the fourth most common gynecological malignancy. Two separate carcinogenic pathways are suggested, where one is associated with the human papillomavirus (HPV) and HPV16 the most common genotype. The aim of this study was to evaluate HPV-markers in a set of primary tumors, metastases and recurrent lesions of vulvar squamous cell carcinomas (VSCC). Ten HPV16-positive VSCC with metastatic regional lymph nodes, distant lymphoid/hematogenous metastases or local recurrent lesions were investigated for HPV genotype, HPV16 variant, HPV16 viral load, HPV16 integration and HPV16 E2BS3 and 4 methylation. In all 10 analyzed case series, the same HPV genotype (HPV16), HPV16 variant and level of viral load were detected in all lesions within a patient case. Primary tumors with a high E2/E6 ratio were found to have fewer vulvar recurrences and/or metastases after diagnosis and treatment. Also, a significantly lower viral load was evident in regional lymph nodes compared to primary tumors. The data presented strengthens the evidence for a clonal HPV-induced pathway for vulvar carcinoma.

Determination of human papillomavirus 16 physical status through E1/E6 and E2/E6 ratio analysis

Human papillomavirus (HPV) 16 genome integration into the host chromosome is a crucial event during the life cycle of the virus and a major step towards carcinogenesis. The integration of HPV16 DNA promotes a constitutive high expression level of E6 and E7 oncoproteins, resulting in the extensive proliferation of the infected epithelial cells. In the present report the physical status of the HPV16 genome was studied, through determination of E1/E6 and E2/E6 DNA copy number ratios in 61 cervical samples of low- and high-grade malignancy and 8 cervical cancer samples, all of them associated with HPV16 infection. The selection of E1, E2 and E6 amplification target regions was performed according to the most prevalent deleted/disrupted sites of E1 and E2 genes. For this target selection we also considered the most conserved regions of E1, E2 and E6 genes among the same HPV16 isolates that were recently reported by our group. The analysis of HPV16 DNA form revealed a significant association among the mixed DNA forms in low-grade and high-grade malignancies, (χ2, P<0.01). The comparative analysis of E1/E6 and E2/E6 in the same cervical samples provides an accurate picture of HPV16 DNA form and may reveal whether different HPV16 DNA integrants coexist in the same cervical sample or not. This study proposes that E1/E6 and E2/E6 ratios determine with accuracy the HPV16 DNA integration pattern and may predict multiple integration events in the examined sample, thus providing significant information about the progression of cervical dysplasia.

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.