Human papillomavirus 16 non-European variants are preferentially associated with high-grade cervical lesions

Lineage and sublineage analysis of human papillomavirus types 51 and 59 in Iranian women

The current study aimed to investigate the sequence variations of HPV 51 and 59 in normal cervical cells and premalignant/malignant lesions of the cervix to know the common variants of HPV 51 and HPV 59 circulating in Iran. To do this, eighty-five samples that were infected by HPV 51 or HPV 59 were investigated using hemi-PCR to amplify the E6 gene followed by sequencing. Our findings indicated that lineages A and B were detected in 80.4% and 19.6% of HPV 51-positive cases, respectively. Among samples infected with HPV 59, 32.2% belonged to lineage A and 67.8% were classified with lineage B. In conclusion, our results showed that lineage A of HPV 51 and lineage B of HPV 59 are more prevalent and distributed in Iran.

Variant analysis of human papillomavirus type 52 in Iranian women during 2018-2020: A case-control study

Background and Aims

Knowing the regional variants of distinct human papillomavirus (HPV) types is valuable as it can be beneficial for studying their epidemiology, pathogenicity, and evolution. For this reason, the sequence variations of the E6 gene of HPV 52 were investigated among women with normal cervical cytology and premalignant/malignant cervical samples.

Methods

Sixty-four HPV 52-positive samples were analyzed using semi-nested PCR and sequencing.

Results

Our findings showed that all samples belonged to lineage A (61%) or B (39%). Among samples that were infected with the A lineage, sublineages A1 and A2 were detected and sublineage A1 was dominant. No association was found between lineages and stage of disease (p > 0.05).

Conclusion

Our results revealed that the A lineage, sublineage A1, and B lineage were common in Iranian women. Nevertheless, more studies with larger sample sizes are required to estimate the pathogenicity risk of HPV 52 lineages in Iranian women with cervical cancer.

HPV16 Phylogenetic Variants in Anogenital and Head and Neck Cancers: State of the Art and Perspectives

HPV16 is responsible for approximately 60% and 90% of global HPV-induced cervical and oropharyngeal cancers, respectively. HPV16 intratype variants have been identified by HPV genome sequencing and classified into four phylogenetic lineages (A-D). Our understanding of HPV16 variants mostly derives from epidemiological studies on cervical cancer (CC) in which HPV16 B, C, and D lineages (previously named "non-European" variants) were mainly associated with high-grade cervical lesions and cancer. Although a predominance of HPV16 lineage A (previously named "European variants") has been observed in head and neck squamous cell carcinoma (HNSCC), epidemiological and in vitro biological studies are still limited for this tumor site. Next Generation Sequencing (NGS) of the entire HPV genome has deepened our knowledge of the prevalence and distribution of HPV variants in CC and HNSCC. Research on cervical cancer has shown that certain HPV16 sublineages, such as D2, D3, A3, and A4, are associated with an increased risk of cervical cancer, and sublineages A4, D2, and D3 are linked to a higher risk of developing adenocarcinomas. Additionally, lineage C and sublineages D2 or D3 of HPV16 show an elevated risk of developing premalignant cervical lesions. However, it is still crucial to conduct large-scale studies on HPV16 variants in different HPV-related tumor sites to deeply evaluate their association with disease development and outcomes. This review discusses the current knowledge and updates on HPV16 phylogenetic variants distribution in HPV-driven anogenital and head and neck cancers.

The E6 gene polymorphism of Human papillomavirus 16 in relation to the risk of cervical cancer in Tunisian women

Human papillomavirus type 16 (HPV-16) is the most prevalent HPV type worldwide and in Tunisia and the major carcinogenic HPV type found in cervical precancers and cancers. Previous studies have reported that genetic diversity of HPV16-E6 oncoprotein might be associated with cervical intraepithelial neoplasia progression. In this study we aimed to investigate the prevalence of HPV-16 E6 variants in precancerous lesions in Tunisian population to assess potential correlation with disease severity. Positive HPV cervical samples were obtained from the Laboratory of Anatomy Pathology of Pasteur Institute of Tunis. Cytological study was performed to identify cervical precancerous lesions. HPVs were typed using Reverse Line Hybridization. Only samples with HPV-16 single infection were selected for HP16-E6 genetic diversity investigation. HPV-16 E6 gene amplification was performed by PCR using specific primers and sequenced by Sanger Sequencing. The multiple alignment of generated sequences was performed using MEGAX software. Phylogenetic tree was constructed using Maximum Likehood method. The ternary complex of E6, E6AP and p53 core domain was used to perform in silico point mutations and thermodynamic calculations to assess stability and binding affinity. Genetic analysis of Tunisian E6-HPV16 sequences showed the presence of three lineages: European (A), African (C) and Asian American (D). Interestingly, the EUR variants were identified as the dominant lineage of HPV-16 and HPV-16 E6 350 G (L83V) was the most detected mutation in precancerous lesions. Modelling data showed that African variants induced the largest destabilizing effect on E6 structure and decreasing thereby in the affinity toward E6AP. Therefore, women infected with European variants are associated with low and high intraepithelial lesions. The findings give useful information for personalized decision algorithms of intra-epithelial cervical neoplasia in Tunisian women.

HPV16 E7 oncoprotein test as a triage strategy for HPV16-positive women in cervical cancer screening: long-term follow-up outcome

Background

Colposcopy is recommended once human papillomavirus (HPV)16/18 infection is detected. However, not all HPV16/18-positive women will necessarily develop cervical lesions. Therefore, this study aimed to investigate the application of quantitative HPV16 E7 oncoprotein detection as a cervical cancer screening method for more efficient screening while minimizing unnecessary colposcopy.

Methods

E7 oncoprotein (HPV16) was quantitatively detected in cervical exfoliated cells of HPV16-positive women. The levels of HPV16 E7 oncoprotein in different degrees of cervical lesions were compared, and the optimal cut-off value for identifying HSIL+ was determined by receiver operating characteristic (ROC) curve analysis. With a pathological diagnosis as the gold standard, the sensitivity (SEN), specificity (SPE), positive predictive value (PPV), negative predictive value (NPV), and Kappa value were calculated to verify the diagnostic value of the method. Women diagnosed with low-grade squamous intraepithelial lesions (LSIL) and normal women were followed up for 5 years to evaluate the predictive value of HPV16 E7 protein for disease progression/persistent infection.

Results

The expression level of HPV16 E7 oncoprotein was positively correlated with the degree of the cervical lesion (r = 0.589, P < 0.01). The area under the ROC curve (AUC) was 0.817 (confidence interval: 0.729-0.904). The cut-off value of E7 oncoprotein for identifying HSIL+ was 8.68 ng/ml. The SEN, SPE, PPV, NPV, and Kappa values of HPV16 E7 oncoprotein for the identification of HSIL+ were 87.1%,70.0%, 87.1%, 70.0%, and 0.571, respectively, which were higher than those of ThinPrep cytology test (TCT). The SEN, SPE, PPV, and NPV of HPV16 E7 oncoprotein in predicting disease progression/persistent infection were 93.75%, 91.30%, 88.24%, and 95.45%, respectively.

Conclusion

The quantitative detection of HPV 16 E7 oncoprotein can not only accurately screen cervical lesions but also achieve efficient colposcopy referral. Additionally, HPV16 E7 oncoprotein can accurately predict the progression of cervical lesions and persistent HPV infection.

Human Papillomavirus 16 Lineage D is Associated with High Risk of Cervical Cancer in the Brazilian Northeast Region

OBJECTIVE

 Similar to Human Papillomavirus (HPV) genotypes, different lineages of a genotype also have different carcinogenic capabilities. Studies have shown that specific genotype lineages of oncogenic HPV are associated with variable risks for the development of cervical intraepithelial neoplasia (CIN2/CIN3) and cervical cancer. The present study aimed to analyze the genetic diversity of the HPV16 genotype in women with CIN2/CIN3 and cervical cancer, from the northeast region of Brazil.

METHODS

 A cross-sectional multicenter study was conducted in the northeast region of Brazil, from 2014 to 2016. This study included 196 cases of HPV16 variants (59 and 137 cases of CIN2/CIN3 and cervical cancer, respectively). The difference of proportion test was used to compare patients with CIN2/CIN3 and cervical cancer, based on the prevalent HPV16 lineage (p < 0.05).

RESULTS

 According to the histopathological diagnosis, the percentage of lineage frequencies revealed a marginal difference in the prevalence of lineage A in CIN2/CIN3, compared with that in cervical cancer (p = 0.053). For lineage D, the proportion was higher in cancer cases (32.8%), than in CIN2/CIN3 cases (16.9%), with p = 0.023.

CONCLUSION

 HPV16 lineage A was the most frequent lineage in both CIN2/CIN3 and cervical cancer samples, while lineage D was predominant in cervical cancer, suggesting a possible association between HPV16 lineage D and cervical cancer.

Impact of HPV-16 Lineages Infection in Response to Radio-Chemotherapy in Cervical Cancer

BACKGROUND

HPV is strongly related to cervical cancer. HPV lineages can contribute to a response to cervical cancer therapy. The aim of this research was to estimate the frequency of human papillomavirus (HPV)-16 lineages in specimens of cervical cancer, relate the pathological factors in these variants, and assess their response to treatment with radical chemoradiotherapy.

METHODS

Samples of cervical cancer were collected from women who were referred to a reference cancer hospital to test the presence of human papillomavirus-type DNA. The standard protocol of this service consisted of cisplatin-based chemotherapy of 40 mg/m², plus conventional pelvic irradiation in doses of 45-50.4 Gy and high dose-rate brachytherapy of 28-30 Gy to Point A. The response to chemotherapy was evaluated after three months in patients with the HPV-16 lineage.

RESULTS

HPV DNA was detected in 104 (88.1%) of the 118 patients. HPV-16 was present in 63 patients (53%). Lineages of HPV-16 were identified in 57 patients and comprised 33 instances of (57.8%) lineage A, 2 instances of lineage B (3.5%), 2 instances of lineage C (3.5%), and 20 instances of (35.0%) lineage D. The median age of the patients was 48.4 years (range 25-85 years). Squamous cell carcinoma was detected 48 times (84.2%). Adenocarcinoma was more likely to occur in lineage D, as three of the four cases occurred in this lineage. A total of 11 patients with the HPV-16 variant were treated with chemoradiotherapy. After three months, it was observed that nine of the eleven patients (81.8%) achieved a complete response, five with the lineage A type, two with the lineage C type, and two with the lineage D type. The two cases of partial response and disease progression, one of each, occurred in lineage A.

CONCLUSIONS

In addition to the small number of patients and HPV variants, we noticed a better response in patients with the HPV-16 lineage A. Increasing the sample size could be helpful to better assess the impact of HPV variants on cervical cancer treatment.

Human papillomavirus genomics: Understanding carcinogenicity

Human papillomavirus (HPV) causes virtually all cervical cancers and many cancers at other anatomical sites in both men and women. However, only 12 of 448 known HPV types are currently classified as carcinogens, and even the most carcinogenic type - HPV16 - only rarely leads to cancer. HPV is therefore necessary but insufficient for cervical cancer, with other contributing factors including host and viral genetics. Over the last decade, HPV whole genome sequencing has established that even fine-scale within-type HPV variation influences precancer/cancer risks, and that these risks vary by histology and host race/ethnicity. In this review, we place these findings in the context of the HPV life cycle and evolution at various levels of viral diversity: between-type, within-type, and within-host. We also discuss key concepts necessary for interpreting HPV genomic data, including features of the viral genome; events leading to carcinogenesis; the role of APOBEC3 in HPV infection and evolution; and methodologies that use deep (high-coverage) sequencing to characterize within-host variation, as opposed to relying on a single representative (consensus) sequence. Given the continued high burden of HPV-associated cancers, understanding HPV carcinogenicity remains important for better understanding, preventing, and treating cancers attributable to infection.

Characterization of the Human Papillomavirus 16 Oncogenes in K14HPV16 Mice: Sublineage A1 Drives Multi-Organ Carcinogenesis

The study of human papillomavirus (HPV)-induced carcinogenesis uses multiple in vivo mouse models, one of which relies on the cytokeratin 14 gene promoter to drive the expression of all HPV early oncogenes. This study aimed to determine the HPV16 variant and sublineage present in the K14HPV16 mouse model. This information can be considered of great importance to further enhance this K14HPV16 model as an essential research tool and optimize its use for basic and translational studies. Our study evaluated HPV DNA from 17 samples isolated from 4 animals, both wild-type (n = 2) and HPV16-transgenic mice (n = 2). Total DNA was extracted from tissues and the detection of HPV16 was performed using a qPCR multiplex. HPV16-positive samples were subsequently whole-genome sequenced by next-generation sequencing techniques. The phylogenetic positioning clearly shows K14HPV16 samples clustering together in the sub-lineage A1 (NC001526.4). A comparative genome analysis of K14HPV16 samples revealed three mutations to the human papillomaviruses type 16 sublineage A1 representative strain. Knowledge of the HPV 16 variant is fundamental, and these findings will allow the rational use of this animal model to explore the role of the A1 sublineage in HPV-driven cancer.

Epidemiology and Molecular Biology of HPV Variants in Cervical Cancer: The State of the Art in Mexico

Cervical cancer (CC) continues to be a major public health problem in Mexico, ranking second among cancers in women. A persistent infection with human papillomaviruses (HPV) is the main risk factor for CC development. In addition, a significant fraction of other cancers including those of the anus, oropharynx, and penis are also related to HPV infection. In CC, HPV-16 is the most prevalent high-risk HPV type, followed by HPV-18, both being responsible for 70% of cases. HPV intratype variant lineages differ in nucleotide sequences by 1–10%, while sublineages differ by 0.5–1%. Several studies have postulated that the nucleotide changes that occur between HPV intratype variants are reflected in functional differences and in pathogenicity. Moreover, it has been demonstrated that HPV-16 and -18 intratype variants differentially affect molecular processes in infected cells, changing their biological behavior that finally impacts in the clinical outcome of patients. Mexico has participated in providing knowledge on the geographical distribution of intratype variants of the most prevalent HPVs in premalignant lesions of the cervix and cervical cancer, as well as in other HPV-related tumors. In addition, functional studies have been carried out to assess the cellular effects of intratype variations in HPV proteins. This review addresses the state of the art on the epidemiology of HPV-16 and HPV-18 intratype variants in the Mexican population, as well as their association with persistence, precancer and cervical cancer, and functional aspects related to their biological behavior.

Deep Sequencing of HPV16 E6 Region Reveals Unique Mutation Pattern of HPV16 and Predicts Cervical Cancer

The genetic diversity of human papillomavirus (HPV) 16 within cervical cells and tissue is usually associated with persistent virus infection and precancerous lesions. To explore the HPV16 mutation patterns contributing to the cervical cancer (CC) progression, a total of 199 DNA samples from HPV16-positive cervical specimens were collected and divided into high‐grade squamous intraepithelial lesion (HSIL) and the non‐HSIL(NHSIL) groups. The HPV16 E6 region (nt 7125-7566) was sequenced using next-generation sequencing. Based on HPV16 E6 amino acid mutation features selected by Lasso algorithm, four machine learning approaches were used to establish HSIL prediction models. The receiver operating characteristic was used to evaluate the model performance in both training and validation cohorts. Western blot was used to detect the degradation of p53 by the E6 variants. Based on the 13 significant mutation features, the logistic regression (LR) model demonstrated the best predictive performance in the training cohort (AUC = 0.944, 95% CI: 0.913–0.976), and also achieved a high discriminative ability in the independent validation cohort (AUC = 0.802, 95% CI: 0.601–1.000). Among these features, the E6 D32E and H85Y variants have higher ability to degrade p53 compared to the E6 wildtype (P < 0.05). In conclusion, our study provides evidence for the first time that HPV16 E6 sequences contain vital mutation features in predicting HSIL. Moreover, the D32E and H85Y variants of E6 exhibited a significantly higher ability to degrade p53, which may play a vital role in the development of CC.

IMPORTANCE The study provides evidence for the first time that HPV16 E6 sequences contain vital mutation features in predicting the high‐grade squamous intraepithelial lesion and can reduce even more unneeded colposcopies without a loss of sensitivity to detect cervical cancer. Moreover, the D32E and H85Y variants of E6 exhibited a significantly higher ability to degrade p53, which may play a vital role in the development of cervical cancer.

Genetic variation within the human papillomavirus type 16 genome is associated with oropharyngeal cancer prognosis

PURPOSE

A significant barrier to adoption of de-escalated treatment protocols for human papillomavirus-driven oropharyngeal cancer (HPV-OPC) is that few predictors of poor prognosis exist. We conducted the first large whole-genome sequencing (WGS) study to characterize the genetic variation of the HPV type 16 (HPV16) genome and to evaluate its association with HPV-OPC patient survival.

PATIENTS AND METHODS

A total of 460 OPC tumor specimens from two large United States medical centers (1980-2017) underwent HPV16 whole-genome sequencing. Site-specific variable positions [single nucleotide polymorphisms (SNPs)] across the HPV16 genome were identified. Cox proportional hazards model estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall survival by HPV16 SNPs. Harrell C-index and time-dependent positive predictive value (PPV) curves and areas under the PPV curves were used to evaluate the predictive accuracy of HPV16 SNPs for overall survival.

RESULTS

A total of 384 OPC tumor specimens (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. Some 284 HPV16 SNPs with a minor allele frequency ≥1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate correction (individual prevalence: 1.0%-5.5%; combined prevalence: 15.10%); E1 gene position 1053 [HR for overall survival (HRos): 3.75, 95% CI 1.77-7.95; Pfdr = 0.0099]; L2 gene positions 4410 (HRos: 5.32, 95% CI 1.91-14.81; Pfdr = 0.0120), 4539 (HRos: 6.54, 95% CI 2.03-21.08; Pfdr = 0.0117); 5050 (HRos: 6.53, 95% CI 2.34-18.24; Pfdr = 0.0030), and 5254 (HRos: 7.76, 95% CI 2.41-24.98; Pfdr = 0.0030); and L1 gene positions 5962 (HRos: 4.40, 95% CI 1.88-10.31; Pfdr = 0.0110) and 6025 (HRos: 5.71, 95% CI 2.43-13.41; Pfdr = 0.0008) and position 7173 within the upstream regulatory region (HRos: 9.90, 95% CI 3.05-32.12; Pfdr = 0.0007). Median survival time for patients with ≥1 high-risk HPV16 SNPs was 3.96 years compared with 18.67 years for patients without a high-risk SNP; log-rank test P < 0.001. HPV16 SNPs significantly improved the predictive accuracy for overall survival above traditional factors (age, smoking, stage, treatment); increase in C-index was 0.069 (95% CI 0.019-0.119, P < 0.001); increase in area under the PPV curve for predicting 5-year survival was 0.068 (95% CI 0.015-0.111, P = 0.008).

CONCLUSIONS

HPV16 genetic variation is associated with HPV-OPC prognosis and can improve prognostic accuracy.

Prevalence and distribution of human papillomavirus (HPV) in Luoyang city of Henan province during 2015-2021 and the genetic variability of HPV16 and 52

Background

Persistent high-risk Human papillomavirus (HPV) subtypes infection has been implicated as a causative of cervical cancer. Distribution and genotypes of HPV infection among females and their variations would assist in the formulation of preventive strategy for cervical cancer. The purpose of the present study is to investigate the prevalence of HPV among females in central China.

Methods

The distribution and genotypes of HPV among 9943 females attending the gynecological examinations in central of China during 2015–2021 were investigated. HPV genotypes were detected using a commercial kit. Nucleotides sequences of L1, E6 and E7 genes in HPV16 or HPV52 positive samples collected in 2021 were amplified by polymerase chain reaction (PCR). Variations of L1, E6 and E7 in HPV16 and HPV52 were gained by sequencing and compared with the reference sequence. Sublineages of HPV16 and HPV52 were determined by the construction of phylogenetic tree based on L1 gene.

Results

The overall prevalence of HPV infection was 22.81%, with the infection rate of high-risk human papillomavirus (HR-HPV) was 19.02% and low-risk human papillomavirus (LR-HPV) was 6.40%. The most top five genotypes of HPV infection were HPV16 (7.49%), HPV52 (3.04%), HPV58 (2.36%), HPV18 (1.65%) and HPV51 (1.61%). Plots of the age-infection rate showed that the single HPV, multiple HPV, HR-HPV, LR-HPV infection revealed the same tendency with two peaks of HPV infection were observed among females aged ≤ 20 year-old and 60–65 year-old. The predominant sublineage of HPV16 was A1 and B2 for HPV52. For HPV16, The most prevalent mutations were T266A (27/27) and N181T (7/27) for L1, D32E for E6 and S63F for E7 in HPV16. For HPV52, all of the nucleotide changes were synonymous mutation in L1 (except L5S) and E7 genes. The K93R mutation was observed in most HPV52 E6 protein.

Conclusions

The present study provides basic information about the distribution, genotypes and variations of HPV among females population in Henan province, which would assist in the formulation of preventive strategies and improvements of diagnostic probe and vaccine for HPV in this region.

HPV16 E6 gene polymorphisms and the functions of the mutation site in cervical cancer among Uygur ethnic and Han nationality women in Xinjiang, China

BACKGROUND

To investigate the genotype distribution of human papillomavirus (HPV) in infected Uygur and Han women in Xinjiang, China; analyze the HPV16 E6 gene polymorphism site and relationship with the development of cervical cancer.

METHODS

The HPV16 E6 sequence was analyzed using the European standard prototype to perform an evolutionary tree. HPV16 E6-T295/T350, G295/G350, and T295/G350 GV230 vectors were stably transfected into cervical cancer C33A cells to analyze the cell proliferation, migration and invasion, apoptosis by CCK8 and clonogenic assays, transwell and cell scratch assays, FACS experiments.

RESULTS

The total HPV infection rate was 26.390% (760/2879), whereas the Uygur 22.87% (196/857) and the Han was 27.89% (564/2022) (P < 0.05). Among 110 mutations, 65 cases of E6 genes were mutated at nucleotide 350 (T350G) with the leucine changing to valine (L83V). Moreover, there were 7 cases of E6 gene mutated at nucleotide 295 (T295G) with aspartic changing to glutamic (D64E). When E6 vector(s) of mutations sites were transfected into C33A cells, they were found to promote cellular proliferation, migration, invasion, and inhibit apoptosis. T295/G350-E6 was significantly stronger than G295/G350 and T295/T350, G295/G350 was significantly stronger than T295/T350 (P < 0.05). The T295/G350 had the strongest effect on C33A cells and G295/G350 was significantly stronger than T295/T350 (P < 0.05).

CONCLUSIONS

The positive HPV infection rates differed between the Uygur and Han in Xinjiang, China, and the genotype distribution of infection was different. After transfecting C33A cells with different eukaryotic expression vectors, the T295/G350 mutation site promoted the proliferation, migration, and invasion of C33A cells to a greater extent than G295/G350; however, G295/G350 had a stronger effect than T295/T350.

Clinical and Molecular Profile of Patients With Condyloma Acuminatum Treated in the Brazilian Public Healthcare System

Condyloma acuminatum is a common clinical outcome of human papillomavirus (HPV) in men. A prospective investigation was performed of the clinical and molecular profile of 122 patients with condyloma acuminatum treated by the Brazilian public healthcare system. The patients were evaluated clinically. The fragments were submitted to molecular analysis for the identification of HPV genotypes. A total of 104 (85.2%) patients presented multiple warts and 18 (14.8%) presented a single wart. The predominant location was the body of the penis (48.4% of cases of multiple warts and 7.4% of cases of single warts), and 49 (40.2%) cases were recurrences and 73 (59.8%) were initial occurrences. Regarding sexual activity, 56 patients (45.9%) had multiple partners and 65 (53.3%) had a single partner. The most frequent genotype was HPV6 (70%). In conclusion, the most frequent anatomic location of condyloma acuminatum was the body of the penis in the present sample. The present findings suggest that the natural history of infection by HPV is not yet completely known and that greater care is needed to ensure clinical safety in the follow-up of these patients due to the oncogenic potential.

APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma

APOBEC is a mutagenic source in human papillomavirus (HPV)-mediated malignancies, including HPV+ oropharyngeal squamous cell carcinoma (HPV + OPSCC), and in HPV genomes. It is unknown why APOBEC mutations predominate in HPV + OPSCC, or if the APOBEC-induced mutations observed in both human cancers and HPV genomes are directly linked. We performed sequencing of host somatic exomes, transcriptomes, and HPV16 genomes from 79 HPV + OPSCC samples, quantifying APOBEC mutational burden and activity in both host and virus. APOBEC was the dominant mutational signature in somatic exomes. In viral genomes, there was a mean of five (range 0-29) mutations per genome. The mean of APOBEC mutations in viral genomes was one (range 0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele fraction mutations, suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection. HPV16 APOBEC-induced mutation patterns in OPSCC were similar to those previously observed in cervical samples. Paired host and viral analyses revealed that APOBEC-enriched tumor samples had higher viral APOBEC mutation rates (p = 0.028), and APOBEC-associated RNA editing (p = 0.008), supporting the concept that APOBEC mutagenesis in host and viral genomes is directly linked and occurrs during infection. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes, we demonstrated for the first-time definitive evidence of concordance between tumor and viral APOBEC mutagenesis. This finding provides a missing link connecting APOBEC mutagenesis in host and virus and supports a common mechanism driving APOBEC dysregulation.

Lineage and sublineage analysis of human papillomavirus type 56 in cervical samples of Iranian women

Understanding the regional lineages and sublineages of human papillomavirus type 56 (HPV 56) would be of great importance for further evolutionary, epidemiological, and biological investigations. To identify the distribution of lineages and sublineages of HPV 56 in Iran, the sequence variations of the E6 gene were analyzed in normal, premalignant, and malignant samples obtained from the cervix. In total, 58 HPV 56-positive samples were investigated by nested-PCR and followed by bidirectional direct nucleotide sequencing analysis. Both lineages A and B were identified in the studied samples. Lineage B was dominant as it was detected in 88.4% of all samples and the remaining samples belonged to lineage A (11.6%). Sublineages A1 and A2 were detected in 3.3% and 8.3% of all samples, respectively. With regard to the pathological stages of cervical specimens, no statistically significant differences were found in the three studied groups (p > 0.05). In conclusion, our findings showed that lineage B of HPV 56 was prevalent in Iran. However, further studies with a larger sample size are warranted to estimate the pathogenicity risk of HPV 56 lineages/sublineages to the progression of cervical cancer among Iranian women.

Lineage analysis of human papillomavirus type 39 in cervical samples of Iranian women

Background

The data with regards to the regional variants of distinct HPV types is of great value. Accordance with this, this study aimed to investigate the sequence variations of E6 gene and long control region of HPV 39 among normal, premalignant and malignant cervical samples in order to characterize the frequent HPV 39 variants circulating in Tehran, Iran.

Methods

In total, 70 cervical samples (45 normal, 16 premalignant, and 9 malignant samples) infected with HPV 39 were analyzed by nested-PCR and sequencing.

Results

Our results revealed that all samples belonged to A lineage. Almost all sequences (98.6%) were classified in A1 sublineage and only one sample (1.4%) was A2 sub lineage.

Conclusions

Our findings showed that lineages A, sublineage A1, is dominant in Tehran, Iran. However, the small sample size was the most important limitations of this study. Further studies with larger sample size from different geographical regions of Iran are necessary to estimate the pathogenicity risk of HPV 39 variants in this population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01619-8.

Association of cervical carcinogenesis risk with HPV16 E6 and E7 variants in the Taizhou area, China

Background

Human papillomavirus (HPV) type 16 accounts for a larger share of cervical cancer and has been a major health problem worldwide for decades. The progression of initial infection to cervical cancer has been linked to viral sequence properties; however, the role of HPV16 variants in the risk of cervical carcinogenesis, especially with longitudinal follow-up, is not fully understood in China.

Methods

We aimed to investigate the genetic variability of HPV16 E6 and E7 oncogenes in isolates from cervical exfoliated cells. Between December 2012 and December 2014, a total of 310 single HPV16-positive samples were selected from women living in the Taizhou area, China. Sequences of all E6 and E7 oncogenes were analysed by PCR-sequencing assay. Detailed sequence comparison, genetic heterogeneity analyses and maximum-likelihood phylogenetic tree construction were performed with BioEdit Sequence Alignment Editor and MEGA X software. Data for cytology tests and histological diagnoses were obtained from our Taizhou Area Study with longitudinal follow-up for at least 5 years. The relationship between HPV16 variants and cervical carcinogenesis risk was analysed by the chi-square test or Fisher’s exact test.

Results

In this study, we obtained 64 distinct variation patterns with the accession GenBank numbers MT681266-MT681329. Phylogenetic analysis revealed that 98.3% of HPV16 variants belong to lineage A, in which the A4 (Asian) sublineage was dominant (64.8%), followed by A2 (12.1%), A1 (11.4%), and A3 (10.0%). The A4 (Asian) sublineage had a higher risk of CIN2+ than the A1–3 (European) sublineages (OR = 2.69, 95% CI = 1.04–6.97, P < 0.05). Furthermore, nucleotide variation in HPV16 E6 T178G is associated with the development of cervical cancer.

Conclusion

These data could provide novel insights into the role of HPV16 variants in cervical carcinogenesis risk in China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08531-y.

Lineage analysis of human papillomavirus types 31 and 45 in cervical samples of Iranian women

Knowing the regional lineages/sublineages of human papillomavirus 31 (HPV 31) and 45 would be of great importance for further evolutionary, epidemiological, and biological analysis. In this regard, to characterize more common lineages and sublineages of HPV 31 and 45, the sequence variations of E6 gene were investigated in normal, premalignant, and malignant samples collected from the cervix in Iran. In total, 54 HPV 31- and 24 HPV 45-positive samples were analyzed by hemi-nested polymerase chain reaction (PCR) and nested-PCR, respectively. All PCR products were subjected to direct sequencing analysis. The results indicated that all three lineages A, B, and C were detected in HPV 31-positive samples; among which HPV 31 lineage A was dominant as it was found in 66.7% of all samples. HPV 31 lineages B and C were identified in 5.5% and 27.8% of samples, respectively. In HPV 45-infected samples, lineage B comprised of 62.5% of all samples and the remaining 37.5%  belonged to lineage A. In conclusion, our findings showed that lineage A of HPV 31 was predominant in Iran. Lineage B of HPV 45 was also dominant among Iranian women. However, further studies with larger sample size should be addressed to estimate the pathogenicity risk of HPV 31 or HPV 45 lineages/sublineages in the development of cervical cancer among Iranian women.

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.

The frequency of high-risk human papillomavirus types, HPV16 lineages, and their relationship with p16INK4a and NF-κB expression in head and neck squamous cell carcinomas in Southwestern Iran

High-risk human papillomaviruses (hr-HPVs) are the key risk factors implicated in the development of a significant proportion of head and neck squamous cell carcinomas (HNSCCs). We aimed to investigate the distribution of hr-HPV types and HPV16 lineages in a sample of patients with HNSCC and the possible association between HPV status and the expression of P16INK4A and NF-κB in Iranian HNSCC patients. We examined 108 formalin-fixed, paraffin-embedded (FFPE) histologically confirmed primary SCC tissue specimens of different head and neck anatomical sites. HPV types and HPV16 lineages were determined by nested PCR and overlapping nested PCR assays, respectively, followed by gene sequencing and phylogenetic analysis. The expression of p16INK4a and NF-κB was evaluated by immunohistochemistry. Twenty-five (23.1%) HNSCC tissue specimens were tested positive for HPV infection. The most prevalent HPV type was HPV-16, followed by HPV18 and HPV11. HPV16 variants belonged to the lineage A and lineage D which were further sorted into sublineages A1, A2, and D2. A significant association between HPV status and p16INK4a immunoreactivity was observed in more than 76% of the HPV-related HNSCCs (P < 0.0001). The overexpression of p16INK4a and cytoplasmic NF-κB was more common in low-grade HNSCC tumors. Our data highlights that HPV16, in particular the A2 sublineage, followed by A1 and D2 sublineages are the major agents associated with HNSCCs in Iran. Based on HPV16 predominance and its lineage distribution pattern, it seems that the prophylactic vaccines developed for cervical cancer prevention could also be applicable for the prevention of HPV-related HNSCCs in our population.

Effects of HPV16 E6 protein on Daxx-induced apoptosis in C33A cells

Aims

Daxx is a highly conserved nuclear protein with an important role in transcription, apoptosis and other cell processes. We investigated the role of HPV16 E6 in Daxx-induced apoptosis through their interactions in C33A cells.

Methods

The binding of HPV16 E6 and Daxx was confirmed in C33A cells using co-immunoprecipitation and indirect immunofluorescence assays. Quantitative PCR and western blotting were performed to determine the RNA and protein expressions of Daxx, respectively. Automatic cell count and MTT assays were performed to investigate the proliferation of C33A cells. The apoptosis rate of C33A cells was determined via flow cytometry using Annexin V-FITC/PI staining. The relative activity of caspase-8 was tested using ELISA.

Results

HPV16 E6 can bind with Daxx and cause its translocation in C33A cells. The transfected HPV16 E6 can cause a decrease in relative quantification for Daxx in Daxx-overexpressing cells. After Daxx transfection, cell proliferation was found to decrease sharply and cell apoptosis to increase sharply. However, when HPV16 E6 was co-transfected with Daxx, this decrease and increase both became gentle. Similarly, HPV16 E6 made the Daxx-induced increase in caspase-8 activity milder.

Conclusions

HPV16 E6 is involved in inhibiting apoptosis through deregulation of Daxx-induced caspase-8 activities.

Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance

HPV16 causes half of cervical cancers worldwide; for unknown reasons, most infections resolve within two years. Here, we analyze the viral genomes of 5,328 HPV16-positive case-control samples to investigate mutational signatures and the role of human APOBEC3-induced mutations in viral clearance and cervical carcinogenesis. We identify four de novo mutational signatures, one of which matches the COSMIC APOBEC-associated signature 2. The viral genomes of the precancer/cancer cases are less likely to contain within-host somatic HPV16 APOBEC3-induced mutations (Fisher's exact test, P = 6.2 x 10-14), and have a 30% lower nonsynonymous APOBEC3 mutation burden compared to controls. We replicate the low prevalence of HPV16 APOBEC3-induced mutations in 1,749 additional cases. APOBEC3 mutations also historically contribute to the evolution of HPV16 lineages. We demonstrate that cervical infections with a greater burden of somatic HPV16 APOBEC3-induced mutations are more likely to be benign or subsequently clear, suggesting they may reduce persistence, and thus progression, within the host.

Analysis of Human Papillomavirus (HPV) 16 Variants Associated with Cervical Infection in Italian Women

This study aims to evaluate HPV16 variants distribution in a population of Italian women living in two different regions (Lombardy and Sardinia) by sequence analyses of HPV16-positive cervical samples, in order to reconstruct the phylogenetic relationship among variants to identify the currently circulating lineages. Analyses were conducted starting from DNA isolated from 67 HPV16-positive cervical samples collected from two different Italian centres (31 from Lombardy and 36 from Sardinia) of women with normal and abnormal cervical cytology. The entire long control region (LCR) and 300 nt of the E6 gene was sequenced to identify intra-type variants. Sequence comparison and phylogenetic analysis were made using a distance-based neighbour joining method (NJ) and Kimura two-parameter model. Data obtained reported that Italian sequences mainly belonged to the European lineage, in particular sublineage A2. Only five sequences clustered in non-European branches: two in North American lineage (sublineage D1), two in African-1 (sublineage B1) and one in African-2. A new 27 nucleotide duplication in the central segment of the LCR region was found in a sequence obtained from a sample isolated in Sardinia. A predominance of European variants was detected, with some degree of variability among the studied HPV16 strains. This study contributes to the implementation of data regarding the molecular epidemiology of HPV16 variants.

Phylogeny and polymorphism in the E6 and E7 of human papillomavirus: alpha-9 (HPV16, 31, 33, 52, 58), alpha-5 (HPV51), alpha-6 (HPV53, 66), alpha-7 (HPV18, 39, 59, 68) and alpha-10 (HPV6, 44) in women from Shanghai

Background

Persistent infection with human papillomaviruses (HPVs) has been associated with cervical intraepithelial neoplasia (CIN) and cervical cancer. However, why only a fraction of HPV cases progress to cancer is still unclear.

Methods

We focused on the heterogeneity, classification, evolution and dispersal of variants for 14 common HPV types in 262 HPV-positive patients with cervical lesions. The E6 and E7 genes of HPV were sequenced and compared with the HPV reference for sequence analysis. Phylogenetic trees were constructed using the neighbour-joining tree method with MEGA 7.0.

Results

In this study, 233 E6 and 212 E7 sequences were successfully amplified by PCR, and these sequences were divided into 5 species groups: alpha-9 (HPV16, 31, 33, 52, 58), alpha-5 (HPV51), alpha-6 (HPV53, 66), alpha-7 (HPV18, 39, 59, 68) and alpha-10 (HPV6, 44). The incidence of high-grade squamous intraepithelial lesion (HSIL) in patients infected with alpha-9 HPV was significantly increased compared with other groups (P < 0.0001), especially HPV16 (P < 0.0001). Strikingly, E7 had significantly fewer nonsynonymous variants in the HSIL compared to <HSIL groups (P = 3.17× 10^- 4). The A388C (K93 N) variation in HPV58 E6 can significantly reduce the risk of HSIL (P = 0.015). However, T7220G (D32E) variation in HPV16 E6 and A7689G (N29S) in HPV16 E7 increased the incidence of HSIL compared to the <HSIL group (P = 0.036 and 0.022).

Conclusions

Strict conservation of E7 is important for HPV carcinogenicity, especially N29 of HPV16. The findings in this work provide preventative/therapeutic interventions for HPV infections and CIN.

Genetic variability and functional implication of HPV16 from cervical intraepithelial neoplasia in Shanghai women

Human papillomavirus (HPV)16 gene mutation is usually associated with persistent HPV infection and cervical intraepithelial neoplasia (CIN). However, the functional implications of HPV16 mutations remain poorly understood.145 LCR/E6/E7 of the HPV16 isolates were amplified and sequenced, and HPV16 integration status was detected. In total, 89 SNPs (68 in the LCR, 13 in E6, 8 in E7) were discovered, 11 of which were nonsynonymous mutations (8 in E6, 3 in E7). The H85Y and E120D variants in E6 were significantly reduced in the high-grade squamous intraepithelial lesion (HSIL) group compared to the <HSIL group (P = .046 and .005), conversely the N29S in E7(P = .01). Amino acid substitutions (D32N/E, E36Q, H85Y, and E120D in E6 and N29H/S and R77C in E7) were predicted to have an effect on conserved structural and functional residues, and five amino acid substitutions (H85Y, E36Q, I34L, and D32E in E6; R77C in E7) would potentially change the secondary structure. "6329G>T," a potential binding site for TATA-binding protein, is the most common in LCR variants. A4 (Asian) was associated with an increased risk of HSIL compared to A1-3(P = .009). The H85/E120 in E6 and N29 in HPV16 E7 might play a critical role in carcinogenesis by disrupting p53 and Rb degradation due to affecting their interaction, respectively. In a word, the findings in this study provide preventative and therapeutic interventions of HPV16 -related cervical lesions/cancer.

Molecular characterisation of genital human papillomavirus among women in Southwestern, Nigeria

Background

Persistent infections with high-risk genital Human papillomavirus (HPV) especially types 16 and 18, are associated with cervical cancer. However, distribution of HPV types varies greatly across geographical regions and the available vaccines target only few types. This study was designed to determine the HPV types circulating in Southwestern Nigeria, thereby providing necessary information for effective control of the virus.

Methods

Endocervical swab samples were collected from a total of 295 consenting women attending routine cervical cancer screening, STI clinics and community-based outreach programme. Viral DNA was extracted from the samples and the consensus region of the HPV DNA was amplified by PCR using GP-E6/E7 primers. Type-specific nested multiplex PCR and Sanger sequencing were used to genotype the HPV isolates.

Results

In this study, 51 (17.3%) individuals were positive for HPV DNA using consensus primers that target the E6/E7 genes but only 48 (16.3%) were genotyped. A total of 15 HPV types (HPV-6, 16, 18, 31, 33, 35, 42, 43, 44, 52, 58, 66, 74, 81, 86) were detected, with HPV-31 being the most predominant (32.8%), followed by HPV-35 (17.2%) and HPV-16 (15.5%). Two rare HPV types; 74 and 86 were also detected. The HPV-74 isolate had three nucleotide (CCT) insertions at E7 gene that translated into amino acid proline. Highest nucleotide substitutions (n = 32) were found in HPV-44 genotype. Among positive individuals, 20.8% had dual infections and 86.2% had High-risk HPV types.

Conclusions

Multiple Human papillomavirus types co-circulated in the study. Most of the circulating Human papillomavirus are high-risk type with type 31 being the most predominant. Although the implication of HPV-74 with proline insertion detected for the first time is unknown, it may have effect on the transformation potential of the virus. Polyvalent HPV vaccine will be more effective for the infection control in Nigeria.

Whole-Genome Analysis of Human Papillomavirus Type 16 Prevalent in Japanese Women with or without Cervical Lesions

Recent large-scale genomics studies of human papillomaviruses (HPVs) have shown a high level of genomic variability of HPV16, the most prevalent genotype in HPV-associated malignancies, and provided new insights into the biological and clinical relevance of its genetic variations in cervical cancer development. Here, we performed deep sequencing analyses of the viral genome to explore genetic variations of HPV16 that are prevalent in Japan. A total of 100 complete genome sequences of HPV16 were determined from cervical specimens collected from Japanese women with cervical intraepithelial neoplasia and invasive cervical cancer, or without cervical malignancies. Phylogenetic analyses revealed the variant distribution in the Japanese HPV16 isolates; overall, lineage A was the most prevalent (94.0%), in which sublineage A4 was dominant (52.0%), followed by sublineage A1 (21.0%). The relative risk of sublineage A4 for cervical cancer development was significantly higher compared to sublineages A1/A2/A3 (odds ratio = 6.72, 95% confidence interval = 1.78–28.9). Interestingly, a novel cluster of variants that branched from A1/A2/A3 was observed for the Japanese HPV16 isolates, indicating that unique HPV16 variants are prevalent among Japanese women.

Human papillomavirus type 16 E6 and E7 gene variations associated with cervical cancer in a Han Chinese population

BACKGROUD

Human papillomavirus type 16 (HPV16) is a high-risk HPV subtype and a potent carcinogen. The HPV16 E6 and E7 genes are considered oncogenes that play a core role in the development of cervical cancer.

METHODS

In the current study, we enrolled 97 HPV16-positive cervical cancer patients (case group) and 136 HPV16-positive asymptomatic individuals (control group) in a study to analyse the association between HPV16 E6 and E7 gene variations and cervical cancer.

RESULTS

Our results showed that three HPV16 sub-lineages (A1-A3, A4 and D3) were present; the distribution of these variants between the case and control group was not significantly different (P = 0.178). When the distribution of the HPV16 E6 and E7 gene variations was compared, the distribution of only A131C (R10R) in the E6 gene showed a different trend between the case and control groups and C749T (S63F) in the E7 gene was significantly different between the case and control groups (P = 0.071 and P = 4.861 × 10^-10, respectively). Regarding the sub-lineages, no variations in the E6 gene were significantly different between the case and control group for the A4 (As) and A1-A3 (EUR) sub-lineages. However, the distribution of C749T (S63F) in the E7 gene was significantly different between the case and control groups for the A4 (As) and A1-A3 (EUR) sub-lineages (P = 1.815 × 10^-8 and P = 0.008). In the current study, we found that the C749T (S63F) variation in the HPV16 E7 gene was associated with cervical cancer not only in the A4 (As) sub-lineage but also in the A1-A3 (EUR) sub-lineage.

CONCLUSION

Our study will provide a good reference for further functional studies of the relationship between cervical cancer carcinogenesis and the HPV16 E6 and E7 genes.

Prevalence of oncogenic Human papillomavirus and genetic diversity in the L1 gene of HPV16 HPV 18 HPV31 and HPV33 found in women from Vojvodina Province Serbia

The main purpose of this paper is to estimate the pre-vaccination prevalence of 12 hrHPV types among 564 women from Vojvodina province (Serbia). The corrected contingency coefficient (Ccorr) was used to estimate the importance of association of examined HPV types and cytological diagnosis. The highest association with the abnormal cytology was observed for HPV 16 (Ccorr = 0.493) in all age groups of participants. The effect of HPV 16 was especially clear within the group of women older than 35 years (Ccorr = 0.691), compared with women younger than 35 (Ccorr = 0.333). The molecular characterization at the level of L1 gene of HPV 16, 18, 31 and 33 variants was for the first time assessed in our region. Nearly all HPV 16 isolates cluster with variant lineage A (96.4%) the remaining isolates clustering with variant lineage D. All of HPV 18 and HPV 33 isolates are clustering within the lineage A while isolates of HPV 31 group with lineages A and C. This contributes to understanding of intrinsic geographical and biological differences of examined HPV types and could be useful for development of cervical cancer screening strategies in Vojvodina (Serbia) and diagnosis of HPV related cervical cancer in general.

Identification of human papillomavirus type 16 variants circulating in the Calabria region by sequencing and phylogenetic analysis of HPV16 from cervical smears

Sequence analysis of HPV16 isolates reveals the presence of genome variants with characteristic mutations. The HPV16 variants have different geographical distribution and diverge into four phylogenetic lineages (A, B, C and D) and 16 sub-lineages: A1, A2, A3 (previously known as European variants), A4 (Asian variant), B1, B2, B3, B4, C1, C2, C3, and C4 (African variants), D1 (North-American variant), D2, D3 (Asian-American variants) and D4. Population studies showed that infections with viruses belonging to specific HPV16 sublineages confer different risks of viral persistence and cancer. In this study, 39 HPV16-positive cervical smears from European women living in Calabria (Italy) were analyzed for the presence of HPV16 variants. Cervical DNA extracts were processed by PCR to amplify L1, the Long Control Region (LCR), E6 and E7, which were sequenced. The sequences were concatenated and the 3169 nucleotides long fragments were characterized by BLAST and phylogenetic analysis. A total of 96 Single Nucleotide Polymorphism (SNPs) were detected, 29 of which mapping in the L1, 45 in the LCR, 15 in the E6 and 7 in the E7. The most common SNP was the T350G (29/39 samples, 74.4%), causing the L83 V amino acid change in the E6. Most of the HPV16 isolates (89.7%) had 99% of nucleotide (nt) identity to members of the A1 and A2 sublineages, while 4 isolates had 99% nt identity to members of the B2, B4, C1 and D4 sublineages. In conclusion, viruses belonging to the A1, A2, B2, B4, C1 and D4 HPV16 sublineages were found to circulate in the Calabria region.

The role of eukaryotic translation initiation factor 5A-1 (eIF5A-1) gene in HPV 16 E6 induces cell growth in human cervical squamous carcinoma cells

Human Papillomavirus (HPV) is considered as the major risk factor for the development and progression of cervical cancer. The high expression of HPV 16 E6 may be the causative factor for induction and maintenance of the transformed phenotype. These oncoproteins would interact with several intracellular proteins, such as eukaryotic translation initiation factor 5A-1 (eIF5A-1) that is essential for proliferation of eukaryotic cells. Our study explored the expression level of HPV 16 E6 and eIF5A-1 in human cervical squamous carcinoma samples and the adjacent non-cancerous cervix samples. Both C33a cells and SiHa cells transfected with a vector encoding HPV 16 E6 resulted in increase of eIF5A-1 expression level and enhancement of viability, migration and proliferation of these cells, furthermore, these effects in both C33a cells and SiHa cells could be abrogated by treatment with eIF5A-1 small-interfering RNA (siRNA) or the specific inhibitors ciclopirox (CPX) that was used to inhibit the function of eIF5A-1 via blocking the main enzymes deoxyhypusine hydroxylase (DOHH). Our results support that the silencing the eIF5A-1 gene or blocking the DOHH could induce the apoptosis of HPV 16 E6-infected cervical carcinoma cells. Thus might provide a new approach to preventing and treating cervical cancer.

Eurogin roadmap 2017: Triage strategies for the management of HPV-positive women in cervical screening programs

Cervical cancer screening will rely, increasingly, on HPV testing as a primary screen. The requirement for triage tests which can delineate clinically significant infection is thus prescient. In this EUROGIN 2017 roadmap, justification behind the most evidenced triages is outlined, as are challenges for implementation. Cytology is the triage with the most follow-up data; the existence of an HR-HPV-positive, cytology-negative group presents a challenge and retesting intervals for this group (and choice of retest) require careful consideration. Furthermore, cytology relies on subjective skills and while adjunctive dual-staining with p16/Ki67 can mitigate inter-operator/-site disparities, clinician-taken samples are required. Comparatively, genotyping and methylation markers are objective and are applicable to self-taken samples, offering logistical advantages including in low and middle income settings. However, genotyping may have diminishing returns in immunised populations and type(s) included must balance absolute risk for disease to avoid low specificity. While viral and cellular methylation markers show promise, more prospective data are needed in addition to refinements in automation. Looking forward, systems that detect multiple targets concurrently such as next generation sequencing platforms will inform the development of triage tools. Additionally, multistep triage strategies may be beneficial provided they do not create complex, unmanageable pathways. Inevitably, the balance of risk to cost(s) will be key in decision making, although defining an acceptable risk will likely differ between settings. Finally, given the significant changes to cervical screening and the variety of triage strategies, appropriate education of both health care providers and the public is essential.

Anal HPV infection and correlates in HIV-infected patients attending a Sexually Transmitted Infection clinic in Brazil

Objective

To estimate the prevalence of anal HPV infection, genotype distribution, intraepithelial neoplasia (AIN) and correlates in a cohort of HIV-infected patients attending at Sexually Transmitted Infections (STI) clinic in Brazil.

Study design

A descriptive analysis was performed which includes, demographic, behavioral and clinical data. Anal specimens from HIV-positive men and women were collected during a regular visit and they were used for cytology and histopathology tests, as well as for HPV molecular identification.

Results

A total of 223 patients (143 females and 80 males) were enrolled in the study and, HPV was identified in 68.6% of the sample. The frequency of HR-HPV, HPV16/18 and multiple HPV infection were similar in both groups. The upstream regulatory region (URR) sequencing was carried out in 38 samples identified as HPV16-positive, and European variants were the most frequent (69.2%), followed by Africans (25.6%) and Asiatic-Americans (5.1%). Having more than 20 sexual partners was associated with multiple HPV infection (p = 0.000) while, anal sex and the first intercourse before 15 years of age was a risk factor for any HPV infection (p = 0.001). Being MSM (men who have sex with men) was a risk factor for any HPV and multiple infections (p = 0.002). The CD4 count >500 cells/mm³ was a protective factor for the HPV16/18 (p = 0.048) and multiple infections (p = 0.023), and the undetectable viral load and HAART treatment were both protective for any HPV (p = 0.010), HR-HPV (p = 0.091) and multiple infections (p = 0.006). Abnormal anoscopy was found in 23.7% (53/223) of the total number of patients, and this was significantly associated with all types of investigated HPV infections (p<0.0001).

Conclusions

In this study, anal HPV infection was common among young HIV-positive men and women, particularly in MSM. Anal cancer screening in patients at risk, such as those who are HIV-positive, and mainly those with anal HPV infection and a history of STI, will increase the likelihood of detecting anal intraepithelial neoplasia.

The Intersection of HPV Epidemiology, Genomics and Mechanistic Studies of HPV-Mediated Carcinogenesis

Of the ~60 human papillomavirus (HPV) genotypes that infect the cervicovaginal epithelium, only 12–13 “high-risk” types are well-established as causing cervical cancer, with HPV16 accounting for over half of all cases worldwide. While HPV16 is the most important carcinogenic type, variants of HPV16 can differ in their carcinogenicity by 10-fold or more in epidemiologic studies. Strong genotype-phenotype associations embedded in the small 8-kb HPV16 genome motivate molecular studies to understand the underlying molecular mechanisms. Understanding the mechanisms of HPV genomic findings is complicated by the linkage of HPV genome variants. A panel of experts in various disciplines gathered on 21 November 2016 to discuss the interdisciplinary science of HPV oncogenesis. Here, we summarize the discussion of the complexity of the viral–host interaction and highlight important next steps for selected applied basic laboratory studies guided by epidemiological genomic findings.

The interaction between mitochondria and oncoviruses

Mitochondria play important roles in multiple aspects of viral tumorigenesis. Mitochondrial genomes contribute to the host's genetic background. After viruses enter the cell, they modulate mitochondrial function and thus alter bioenergetics and retrograde signaling pathways. At the same time, mitochondria also regulate and mediate viral oncogenesis. In this context, oncogenesis by oncoviruses like Hepatitis B virus (HBV), Hepatitis C virus (HCV), Human papilloma virus (HPV), Human Immunodeficiency virus (HIV) and Epstein-Barr virus (EBV) will be discussed.

Human papillomavirus type 16 lineage analysis based on E6 region in cervical samples of Iranian women

It is suggested that distinct HPV 16 variants differ in oncogenic potential and geographic distribution. As such, understanding the regional variants of HPV 16 would be of great importance for evolutionary, epidemiological and biological analysis. In this regard, the sequence variations of E6 gene were investigated to characterize more common variants of HPV 16 in normal cells, premalignant and malignant lesions of the cervix. In total, 106 isolates of HPV 16 were analyzed by PCR and sequencing. Overall, two different lineages (A and D) were identified. Lineage D comprised 70.7% of samples and the remaining 29.3% belonged to lineage A. Regarding to cytology/histology, lineage D was dominant in both normal+CIN I-II and CIN III+ICC groups as it was detected in 80% and 66.2% of cases, respectively. The comparison of the lineages between different groups (35 normal+CIN I-II samples and 71 CIN III+ICC samples) revealed that lineage A is more prevalent in cervical cancer cases (7 (20%) vs. 24 (33.8%)) although the difference observed did not reach statistical significance (p=0.07). In conclusion, our findings confirm that HPV lineages A and D are more prevalent in Iran, with the lineage D as the most dominant in all studied groups.

Genetic diversity of human papillomavirus types 35, 45 and 58 in cervical cancer in Brazil

In Brazil, most studies of intra-type variants of human papillomavirus (HPV) have focused on HPV16 and HPV18, but other high-risk HPV types have not been studied. Here, we report the prevalence of lineages and variants of HPV35, HPV45 and HPV58 in cervical cancers from the Amazonian and Southeast Brazilian regions. The most frequent sublineages were A1 for HPV35, B2 for HPV45, and A2 for HPV58. The Southeast region had a higher frequency of the B2 sublineage of HPV45, and for HPV35, the genetic and nucleotide sequence diversity were higher in the Southeast region, suggesting that regional factors are influencing the diversity and lineage prevalence.

The high prevalence of HPV and HPV16 European variants in cervical and anal samples of HIV-seropositive women with normal Pap test results

Human Immunodeficiency Virus (HIV)-seropositive women are more likely to have anogenital cancer, and high risk-HPV (HR-HPV) infection is the main associated factor. Between August 2013 and December 2015, we conducted a descriptive study to determine the HPV genotypes and HPV16 variants in cervical and anal samples of HIV-seropositive women with a normal Pap test. The viral DNA was amplified by PCR using the PGMY09/11 set of primers. Reverse line blot (RLB), restriction fragment length polymorphism (RFLP) and sequencing assays were used to determine the HPV genotypes. HPV16 variants were identified by gene sequencing. We found a high frequency of HR-HPV (60.3%; 76/126) at the anogenital site among HIV-seropositive women and without association with anal intercourse. HPV16 and European variant predominated among the HR-HPV. Mixed infections with at least three different HPV types were common, particularly at the anal site. CD4+ T-cell counts below 500 cells/mm³, a HIV viral load above 50 copies/mL and an age of 18 to 35 years old were all related to HPV anal infection. Our study showed a high frequency of HR-HPV in both cervical and anal sites of women with negative cytology belonging to a risk group for the development of anogenital cancer.

Risk Factors for Cervical Cancer and CIN3 in Jewish Women in Israel - Two Case Control Studies

PURPOSE

The aim of the study was to identify risk and protective factors/markers for cervical cancer and cervical intraepithelial neoplasia 3 (CIN3) in Israeli Jewish women in order to settle the discrepancy of low incidence rate of cervical cancer and relatively high incidence rate of CIN3.

MATERIALS AND METHODS

We conducted two case control studies, which examined the association between potential risk and protective factors/markers for cervical cancer or CIN3 using self administered detailed questionnaires.

RESULTS

For studying cervical cancer, 40 cases and 40 matched controls were interviewed. In the univariable and multivariable analyses older age, depression or anxiety and ever smoking seemed to act as independent risk factors/markers, while older age at first intercourse was protective. For studying CIN3, 99 cases and 79 controls were interviewed. Multivariable analysis has demonstrated that being born in Israel, depression or anxiety and ever smoking were independent risk factors/markers for CIN3.

CONCLUSIONS

The risk factors/markers studied, that were associated with cervical cancer or CIN3 among Jewish women in Israel, are similar to those reported in other parts of the world, and do not explain the observed discrepancy of high in-situ cervical cancer rates and low invasive cervical cancer incidence in Israel.

Transmission between Archaic and Modern Human Ancestors during the Evolution of the Oncogenic Human Papillomavirus 16

Every human suffers through life a number of papillomaviruses (PVs) infections, most of them asymptomatic. A notable exception are persistent infections by Human papillomavirus 16 (HPV16), the most oncogenic infectious agent for humans and responsible for most infection-driven anogenital cancers. Oncogenic potential is not homogeneous among HPV16 lineages, and genetic variation within HPV16 exhibits some geographic structure. However, an in-depth analysis of the HPV16 evolutionary history was still wanting. We have analyzed extant HPV16 diversity and compared the evolutionary and phylogeographical patterns of humans and of HPV16. We show that codivergence with modern humans explains at most 30% of the present viral geographical distribution. The most explanatory scenario suggests that ancestral HPV16 already infected ancestral human populations and that viral lineages co-diverged with the hosts in parallel with the split between archaic Neanderthal-Denisovans and ancestral modern human populations, generating the ancestral HPV16A and HPV16BCD viral lineages, respectively. We propose that after out-of-Africa migration of modern human ancestors, sexual transmission between human populations introduced HPV16A into modern human ancestor populations. We hypothesize that differential coevolution of HPV16 lineages with different but closely related ancestral human populations and subsequent host-switch events in parallel with introgression of archaic alleles into the genomes of modern human ancestors may be largely responsible for the present-day differential prevalence and association with cancers for HPV16 variants.

The differential role of HTRA1 in HPV-positive and HPV-negative cervical cell line proliferation

BACKGROUND

High-risk human papillomaviruses (HPVs) are strongly associated with the development of some malignancies. The E6 and E7 viral oncoproteins are the primary proteins responsible for cell homeostasis alteration and immortalization. Furthermore, the E6 protein from high-risk HPVs can interact with the PDZ (PSD-90/Dlg/ZO-1) domains of cellular proteins, triggering cell transformation. One protein that is associated with pathological conditions and has a PDZ domain is the protease HTRA1 (high temperature requirement 1). This protein is poorly expressed in some cancers, suggesting a tumor suppressor role. The aim of this study was to evaluate the effect of HTRA1 overexpression in HPV16-positive (CasKi) and HPV-negative (C33) cervical cell lines.

METHODS

The cells were transfected with a vector containing the HTRA1 ORF or an empty vector. HTRA1 overexpression was confirmed by qRT-PCR. The cells were subjected to cell proliferation, colony formation, apoptosis and cell cycle assays.

RESULTS

C33 cells expressing HTRA1 grew significantly fewer colonies and showed less proliferation than cells without HTRA1 expression. In contrast, in the CasKi cells overexpressing HTRA1, there was an increase in the cell growth rate and in the colonies density compared to cells expressing low levels of HTRA1. An apoptosis assay showed that HTRA1 does not interfere with the apoptosis rate in these cells. A cell cycle immunofluorescence assay revealed more CasKi cells overexpressing HTRA1 in the S phase and more C33 HTRA1-transfected cells in the G0/G1 phase, suggesting that HTRA1 plays different roles in the cell cycle progression of these cells.

CONCLUSIONS

HTRA1 overexpression prevents cell proliferation in the HPV-negative cell line and increases cell proliferation in the HPV-positive cell line. Although the E6/HTRA1 interaction has already been described in the literature, more studies are required to confirm whether the present functional findings are a result of this interaction.

Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology

BACKGROUND

Human papillomaviruses (HPVs) are a worldwide burden as they are a widespread group of tumour viruses in humans. Having a tropism for mucosal tissues, high-risk HPVs are detected in nearly all cervical cancers. HPV16 is the most common high-risk type but not all women infected with high-risk HPV develop a malignant tumour. Likely relevant, HPV genomes are polymorphic and some HPV16 single nucleotide polymorphisms (SNPs) are under evolutionary constraint instigating variable oncogenicity and immunogenicity in the infected host.

RESULTS

To investigate the tumourigenicity of two common HPV16 variants, we used our recently developed, three-dimensional organotypic model reminiscent of the natural HPV infectious cycle and conducted various "omics" and bioinformatics approaches. Based on epidemiological studies we chose to examine the HPV16 Asian-American (AA) and HPV16 European Prototype (EP) variants. They differ by three non-synonymous SNPs in the transforming and virus-encoded E6 oncogene where AAE6 is classified as a high- and EPE6 as a low-risk variant. Remarkably, the high-risk AAE6 variant genome integrated into the host DNA, while the low-risk EPE6 variant genome remained episomal as evidenced by highly sensitive Capt-HPV sequencing. RNA-seq experiments showed that the truncated form of AAE6, integrated in chromosome 5q32, produced a local gene over-expression and a large variety of viral-human fusion transcripts, including long distance spliced transcripts. In addition, differential enrichment of host cell pathways was observed between both HPV16 E6 variant-containing epithelia. Finally, in the high-risk variant, we detected a molecular signature of host chromosomal instability, a common property of cancer cells.

CONCLUSIONS

We show how naturally occurring SNPs in the HPV16 E6 oncogene cause significant changes in the outcome of HPV infections and subsequent viral and host transcriptome alterations prone to drive carcinogenesis. Host genome instability is closely linked to viral integration into the host genome of HPV-infected cells, which is a key phenomenon for malignant cellular transformation and the reason for uncontrolled E6 oncogene expression. In particular, the finding of variant-specific integration potential represents a new paradigm in HPV variant biology.

Genomic diversity and phylogenetic relationships of human papillomavirus 16 (HPV16) in Nepal

OBJECTIVES/BACKGROUND

Sequence variants in HPV16 confer differences in oncogenic potential; however, to date there have not been any HPV sequence studies performed in Nepal. The objective of this study was to characterize HPV16 viral genome sequences from Nepal compared to a reference sequence in order to determine their lineages. Additionally, we sought to determine if five High-grade Squamous Intraepithelial Lesion (HSIL) subjects were genetically distinct from the non-HSIL subjects.

METHODS

DNA was isolated from exfoliated cervical cells from 17 individuals in Nepal who were previously identified to be HPV16-positive. A custom HPV16 Ion Ampliseq panel of multiplexed degenerate primers was designed that generated 47 overlapping amplicons and covered 99% of the viral genome for all known HPV16 variant lineages. All sequence data were processed through a custom quality control and analysis pipeline of sequence comparisons and phylogenetic analysis.

RESULTS

There were high similarities across the genomes, with two major indels observed in the non-coding region between E5 and L2. Compared to the PAVE reference HPV16 genome, there were up to 9, 4, 38, 27, 8, 7, 52, and 32 nucleotide variants in the E6, E7, E1, E2, E4, E5, L2, and L1 genes in the Nepalese samples, respectively. Based on sequence variation, HPV16 from Nepal falls across the A, C, and D lineages in this study. We found no evidence of genetic distinctness between HSIL and non-HSIL subjects.

CONCLUSIONS

The evolutionary and pathological characteristics of the representative HPV16 genomes from Nepal seem similar to results from other parts of the world and provide the basis for further studies.

HPV16 variants distribution in invasive cancers of the cervix, vulva, vagina, penis, and anus

Human papillomavirus (HPV)16 is the most oncogenic human papillomavirus, responsible for most papillomavirus‐induced anogenital cancers. We have explored by sequencing and phylogenetic analysis the viral variant lineages present in 692 HPV16‐monoinfected invasive anogenital cancers from Europe, Asia, and Central/South America. We have assessed the contribution of geography and anatomy to the differential prevalence of HPV16 variants and to the nonsynonymous E6 T350G polymorphism. Most (68%) of the variance in the distribution of HPV16 variants was accounted for by the differential abundance of the different viral lineages. The most prevalent variant (above 70% prevalence) in all regions and in all locations was HPV16_A1‐3, except in Asia, where HPV16_A4 predominated in anal cancers. The differential prevalence of variants as a function of geographical origin explained 9% of the variance, and the differential prevalence of variants as a function of anatomical location accounted for less than 3% of the variance. Despite containing similar repertoires of HPV16 variants, we confirm the worldwide trend of cervical cancers being diagnosed significantly earlier than other anogenital cancers (early fifties vs. early sixties). Frequencies for alleles in the HPV16 E6 T350G polymorphism were similar across anogenital cancers from the same geographical origin. Interestingly, anogenital cancers from Central/South America displayed higher 350G allele frequencies also within HPV16_A1‐3 lineage compared with Europe. Our results demonstrate ample variation in HPV16 variants prevalence in anogenital cancers, which is partly explained by the geographical origin of the sample and only marginally explained by the anatomical location of the lesion, suggesting that tissue specialization is not essential evolutionary forces shaping HPV16 diversity in anogenital cancers.

HPV16 Sublineage Associations With Histology-Specific Cancer Risk Using HPV Whole-Genome Sequences in 3200 Women

BACKGROUND

HPV16 is a common sexually transmitted infection although few infections lead to cervical precancer/cancer; we cannot distinguish nor mechanistically explain why only certain infections progress. HPV16 can be classified into four main evolutionary-derived variant lineages (A, B, C, D) that have been previously suggested to have varying disease risks.

METHODS

We used a high-throughput HPV16 whole-genome sequencing assay to investigate variant lineage risk among 3215 HPV16-infected women. Using sublineages A1/A2 as the reference, we assessed all variant lineage associations with infection outcome over three or more years of follow-up: 1107 control subjects (

RESULTS

A4 sublineage was associated with an increased risk of cancer, specifically adenocarcinoma (OR = 9.81, 95% CI = 2.02 to 47.69, P = 4.7x10(-03)). Lineage B had a lower risk of CIN3 (OR = 0.51, 95% CI = 0. 28 to 0.91, P = 02) while lineage C showed increased risk (OR = 2.06, 95% CI = 1.09 to 3.89, P = 03). D2/D3 sublineages were strongly associated with an increased risk of CIN3 and cancer, particularly D2 (OR for cancer = 28.48, 95% CI = 9.27 to 87.55, P = 5.0x10(-09)). D2 had the strongest increased risk of glandular lesions, AIS (OR = 29.22, 95% CI = 8.94 to 95.51, P = 2.3x10(-08)), and adenocarcinomas (OR = 137.34, 95% CI = 37.21 to 506.88, P = 1.5x10(-13)). Moreover, the risk of precancer and cancer for specific variant lineages varied by a women's race/ethnicity; those women whose race/ethnicity matched that of the infecting HPV16 variant had an increased risk of CIN3 + (P < 001).

CONCLUSIONS

Specific HPV16 variant sublineages strongly influence risk of histologic types of precancer and cancer, and viral genetic variation may help explain its unique carcinogenic properties.

Collapse

Key Words Collapse

MESH Headings

• Adenocarcinoma/ethnology
• Adenocarcinoma/virology
• Adenocarcinoma in Situ/ethnology
• Adenocarcinoma in Situ/virology
• Adult
• Aftercare
• California/epidemiology
• Carcinoma/ethnology
• Carcinoma/virology
• Carcinoma, Squamous Cell/ethnology
• Carcinoma, Squamous Cell/virology
• Female
• Genome
• High-Throughput Nucleotide Sequencing
• Human papillomavirus 16/classification
• Human papillomavirus 16/genetics
• Human papillomavirus 16/pathogenicity
• Humans
• Middle Aged
• Papillomavirus Infections/ethnology
• Papillomavirus Infections/virology
• Phylogeny
• Precancerous Conditions/ethnology
• Precancerous Conditions/virology
• Risk Factors
• Uterine Cervical Neoplasms/ethnology
• Uterine Cervical Neoplasms/virology
• Young Adult
• Uterine Cervical Dysplasia/ethnology
• Uterine Cervical Dysplasia/virology

Collapse

Grants

• U01 CA078527 NCI NIH HHS

Collapse

Affiliation(s)

Lisa Mirabello Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Meredith Yeager Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Michael Cullen Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Joseph F Boland Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Zigui Chen Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Nicolas Wentzensen Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Xijun Zhang Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Kai Yu Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Qi Yang Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Jason Mitchell Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
David Roberson Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Sara Bass Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Yanzi Xiao Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Laurie Burdett Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Tina Raine-Bennett Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Thomas Lorey Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Philip E Castle Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Robert D Burk Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Mark Schiffman Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)

Collapse

Variants of human papillomaviruses 16 (HPV16) in Uigur women in Xinjiang, China

BACKGROUND

Persistent infection of high-risk human papillomaviruses 16 (HPV16) has been considered as the leading cause of cervical cancer. In this study we assessed HPV16 sequence variation and genetic diversity of HPV16 variants in cervical cancer in Uigur women in Xinjiang, China. We analyzed the nucleotide sequences of the open reading frames of E6 and E7, and part of the open reading frames of L1 of HPV16 in Uigur women.

METHODS

Biopsies of histologically confirmed HPV16 infections with cervical cancer were obtained from 43 Uigur women in Xinjiang, China. E6, E7 and L1 genes of HPV16 of all samples were amplified and sequenced; the sequences were used in phylogenetic analysis of HPV16 variants.

RESULTS

Our analysis revealed nine nucleotide changes in E6 (five changes), E7 (one change) and L1 (three changes) gene. The most frequently observed variations were T350G (79.1 %). One variation T295G (D64E) at E6 were detected in 6 cases (KT959536, KT959542, KT959546, KT959550, KT959553, KT959558). Deletion (464Asp) along with insertion (448Ser) were observed in L1 (100 %). Most variants were European lineage (97.7 %); only one belongs to Asia variants with common T178G (D25E) in E6 and A647G (N29S) in E7.

CONCLUSION

The most prevalent HPV16 variants in the Uigur women we studied were of the European lineage. Our results indicate that HPV16 European lineage may serve as a harmful factor associated with the development and progression of cervical cancer.