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

Distribution of human papillomavirus among Vietnamese women with cervical cancer and unusual genetic variability of HPV16

HPV16, with typical mutations that differ in geographical distribution and carcinogenic potency, has implications for cervical cancer screening, clinical diagnosis, and treatment. DNASTAR and MEGA were used to identify HPV16 variants and construct a phylogenetic tree. The most prevalent HPV genotypes were HPV16 (63.9%), HPV18 (26.7%), and other HPV (6.9%). HPV16 alterations were found in all E6, E7, and L1 genes, including 15 missense and 18 synonymous mutations. Missense mutations include R10G, Q14H, D25E, H78Y, L83V (E6); M29V, R35K, L78R, L95P (E7); H73Y, T176 N, N178T, T317P, T386S, L472F/I (L1). HPV16 sublineages include A1 (17.2%), A2 (0.9%), A3 (56.0%), A4 (19.0%), D1 (4.3%), and D3 (2.6%). Although several mutations in the oncoproteins E6, E7, and L1 have been detected, mutations known to be associated with cervical cancer risk, such as D25E and L83V, occur at a relatively low frequency. This suggests that HPV16 mutations are associated with cervical cancer through a complicated mechanism.

The difference of transcriptome of HPV-infected patients contributes more to the occurrence of cervical cancer than the mutations of E6 and E7 genes in HPV16

Human papillomavirus (HPV) E6 and E7 genes are biomarkers and drivers of the progression of cervical cancer (CxCa). The aim of this study was to investigate the relationship between HPV16 E6, E7 gene mutations and the occurrence and development of CxCa. Cervical exfoliated cells and clinical data of patients with cervical diseases were collected. Sample DNA was extracted, the E6 and E7 gene fragments were amplified by PCR, and the mutations were detected by Sanger sequencing and compared with standard sequences. Microarray was used to sequence the transcriptome of cells. Data of transcriptome analyzed and visualized using R software and its packages. Analysis of clinical characteristics demonstrated the association of HPV16 infection with CxCa (P < .05). Sanger sequencing results showed that the mutation sites of E6 gene included T178G/A, T350G, A131C, and T241G; among these, A131C and T241G were synonymous mutations. The mutation sites of E7 gene included A647G, T846C, G666A, T843C, and T760C, and all of them were synonymous mutations except A647G. There was no significant difference in the distribution of HPV16 E6, E7 mutations among CxCa, cervical intraepithelial neoplasia, and infection groups (P > .05). Compared with the non- CxCa group, gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of differentially expressed genes (DEGs) showed more significant enrichment of DEGs in the biological processes, pathways, and diseases closely related to cancer. Compared with the non-mutation group, the DEGs in the E6, E7 gene mutation group were significantly enriched in the events related to infection and immunity. To summarize, HPV16 may be associated with the occurrence and development of CxCa, but HPV16 E6 and E7 gene mutations have little effect on the occurrence and development of CxCa. Individual differences may have a greater effect on the progression of CxCa.

Genetic variation of E6 and E7 genes of human papillomavirus type 16 from central China

BACKGROUND

Persistent high-risk human papillomavirus (HR-HPV) infection is an important factor in the development of cervical cancer, and human papillomavirus type 16 (HPV-16) is the most common HR-HPV type worldwide. The oncogenic potential of HPV-16 is closely related to viral sequence variation.

METHODS

In order to clarify the variant characteristics of HPV-16 E6 and E7 genes in central China, E6 and E7 sequences of 205 HPV-16 positive samples were amplified by polymerase chain reaction. PCR products of E6 and E7 genes were further sequenced and subjected to variation analysis, phylogenetic analysis, selective pressure analysis and B-cell epitope prediction.

RESULTS

Twenty-six single nucleotide variants were observed in E6 sequence, including 21 non-synonymous and 5 synonymous variants. Twelve single nucleotide variants were identified in E7 sequence, including 6 non-synonymous and 6 synonymous variants. Four new variants were found. Furthermore, nucleotide variation A647G (N29S) in E7 was significantly related to the higher risk of HSIL and cervical cancer. Phylogenetic analysis showed that the E6 and E7 sequences were all distributed in A lineage. No positively selected site was found in HPV-16 E6 and E7 sequences. Non-conservative substitutions in E6, H31Y, D32N, D32E, I34M, L35V, E36Q, L45P, N65S and K75T, affected multiple B-cell epitopes. However, the variation of E7 gene had little impact on the corresponding B-cell epitopes (score < 0.85).

CONCLUSION

HPV-16 E6 and E7 sequences variation data may contribute to HR-HPV prevention and vaccine development in Jingzhou, central China.

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.

The polymorphism analysis and epitope predicted of Alphapapillomavirus 9 E6 in Sichuan, China

Background

The Alphapapillomavirus 9 (α-9 HPV) is a member of the Alphapapillomavirus genus and Papillomaviridae family. These viruses are almost all carcinogenic HPV, which is closely related to 75% of invasive cervical cancer worldwide, and has a high prevalence in Sichuan. The carcinogenic function is mainly realized by its E6 oncoprotein.

Methods

Cell samples were collected by cervical scraped for HPV detecting and typing. HPV-16, HPV-31, HPV-33, HPV-52, HPV-58 5 α-9 genus HPV subtype positive samples were selected, their E6 gene was sequenced and analyzed. The positive selection sites of HPV E6 genes were estimated by PAML 4.8 server. The secondary and tertiary structure of E6 protein were predicted by PSIPred and Swiss-model. The T-cell antigen epitopes of E6 protein were predicted by IEDB.

Results

α-9 HPV has a high prevalence in Sichuan, China. From 2012 to 2017, 18,067 cell cervical samples were collected, and 3135 were detected with α-9 HPV infection. Among which, 250 cases HPV-16 E6, 96 cases HPV-31 E6, 216 cases HPV-33 E6, 288 cases HPV-52 E6 and 405 cases HPV-58 E6 were successfully amplified, 17, 6, 6, 13, and 4 non-synonymous nucleotide mutations were respectively detected in HPV-16, 31, 33, 52, and 58 E6, 7 positive selection sites of α-9 HPV E6 were selected out (D32E of HPV-16 E6, K35N, K93N and R145I of HPV-33 E6, K93R of HPV-52 E6, K93N and R145K of HPV-58 E6). The structure and antigen epitopes of E6 protein with amino acid substitution differ from those of wild-type E6 protein, especially for the mutation located in the E6 positive selection site.

Conclusions

HPV E6 nucleotide non-synonymous mutation in the positive selection site influence the protein structure and decrease the antigen epitopes affinity of the E6 protein overall, making it more difficult for the HPV-infected cells to be detected by the immune system, and enhancing the HPV adaptability to the environment. Mutations influence the validity of HPV clinical diagnostic probes, the polymorphism analysis of α-9 HPV E6 enrich the data of HR-risk HPV in Sichuan China, and the detection probes designed with the polymorphism data in mind can improve the efficiency of clinical detection; Mutations influence epitopes affinity, the association of E6 polymorphism and epitope affinity can improve the design of therapeutic vaccine with good immunity and high generality antigen epitope; The above study all provide a good theoretical basis for the prevention and treatment of HPV-related diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01728-4.

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

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

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

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

Mutation Profiles, Glycosylation Site Distribution and Codon Usage Bias of Human Papillomavirus Type 16

Human papillomavirus type 16 (HPV16) is the most prevalent HPV type causing cervical cancers. Herein, using 1597 full genomes, we systemically investigated the mutation profiles, surface protein glycosylation sites and the codon usage bias (CUB) of HPV16 from different lineages and sublineages. Multiple lineage- or sublineage-conserved mutation sites were identified. Glycosylation analysis showed that HPV16 lineage D contained the highest number of different glycosylation sites from lineage A in both L1 and L2 capsid proteins, which might lead to their antigenic distances between the two lineages. CUB analysis showed that the HPV16 open reading frames (ORFs) preferred codons ending with A/T. The CUB of HPV16 ORFs was mainly affected by natural selection except for E1, E5 and L2. HPV16 only shared some of the preferred codons with humans, which might help reduce competition in translational resources. These findings increase our understanding of the heterogeneity between HPV16 lineages and sublineages, and the adaptation mechanism of HPV in human cells. In summary, this study might facilitate HPV classification and improve vaccine development and application.

E6/E7 Variants of Human Papillomavirus 16 Associated with Cervical Carcinoma in Women in Southern Mexico

Persistent infection with the human papillomavirus 16 (HPV 16) is the cause of half of all cervical carcinomas (CC) cases. Moreover, mutations in the oncoproteins E6 and E7 are associated with CC development. In this study, E6/E7 variants circulating in southern Mexico and their association with CC and its precursor lesions were evaluated. In total, 190 DNA samples were obtained from scrapes and cervical biopsies of women with HPV 16 out of which 61 are from patients with CC, 6 from patients with high-grade squamous intraepithelial lesions (HSIL), 68 from patients with low-grade squamous intraepithelial lesions (LSIL), and 55 from patients without intraepithelial lesions. For all E7 variants found, the E7-C732/C789/G795 variant (with three silent mutations) was associated with the highest risk of CC (odd ratio (OR) = 3.79, 95% confidence interval (CI) = 1.46–9.85). The analysis of E6/E7 bicistron conferred to AA-a*E7-C732/C789/G795 variants revealed the greatest increased risk of CC (OR = 110, 95% CI = 6.04–2001.3), followed by AA-c*E7-C732/C789/G795 and A176/G350*E7-p. These results highlight the importance of analyzing the combinations of E6/E7 variants in HPV 16 infection and suggest that AA-a*E7-C732/C789/G795, AA-c*E7-C732/C789/G795, and A176/G350*E7-p can be useful markers for predicting CC development.

Metabolic Reprogramming in Cancer: Role of HPV 16 Variants

Metabolic reprogramming is considered one of the hallmarks in cancer and is characterized by increased glycolysis and lactate production, even in the presence of oxygen, which leads the cancer cells to a process called “aerobic glycolysis” or “Warburg effect”. The E6 and E7 oncoproteins of human papillomavirus 16 (HPV 16) favor the Warburg effect through their interaction with a molecule that regulates cellular metabolism, such as p53, retinoblastoma protein (pRb), c-Myc, and hypoxia inducible factor 1α (HIF-1α). Besides, the impact of the E6 and E7 variants of HPV 16 on metabolic reprogramming through proteins such as HIF-1α may be related to their oncogenicity by favoring cellular metabolism modifications to satisfy the energy demands necessary for viral persistence and cancer development. This review will discuss the role of HPV 16 E6 and E7 variants in metabolic reprogramming and their contribution to developing and preserving the malignant phenotype of cancers associated with HPV 16 infection.

Prevalence of human papillomavirus in archival head and neck cancer in the eastern Inner Mongolian Autonomous Region, China

OBJECTIVE

To investigate the prevalence of human papillomavirus (HPV) in archival head and neck cancer (HNC) collected from the Tong-Liao area, which is located in east Inner Mongolia, China.

METHODS

The presence of HPV in 54 HNCs and 25 benign biopsies was detected and the sequence variation of the E6 gene in HPV-positive samples was analyzed to determine their lineage/sublineage classification.

RESULTS

HPV was detected in only 4 out of 54 HNCs and no benign biopsies were found to be HPV-positive. After further p16^INK4a immunostaining, only 3 cases of HNC were positive for both HPV and p16^INK4a. Phylogenetic analysis of the isolated E6 gene shows that the HPV 16, HPV 31, and HPV 58 isolated in this study belong to lineage A.

CONCLUSIONS

The prevalence of HPV in HNC from this area is very low. The lineage/sublineage classification of the 3 HPV types in HNC in this area is consistent with the previous reported data of HPV lineage distribution in cervical cancer within China.

Distribution of HPV 16 E6 gene variants in screening women and its associations with cervical lesions progression

The aim of this study is to investigate distribution of human papillomavirus (HPV) 16 variants in screening healthy women and the potential association between HPV 16 variants and progression of cervical lesions. For this study a total of 2000 healthy women in Tianjin urban area and 212 patients who were HPV 16 positive and underwent colposcopy were analyzed for HPV 16 variants by pyrosequencing. The results show that the HPV 16 was the most prevalent genotype in Tianjin healthy women and five HPV 16 variant types were detected. The HPV 16 variants were determined by sequencing parital E6 region and the detected variants were European prototype E-T350 (E-p), E-G350, E-C109 G, Asian (As) and Asian-American (AA), among which the E-p variant was the most prevalent (82.76%) followed by As variant. Interestingly, in patients with suspected cervical lesions the most prevalent variant was As variant (54.9%) by increasing significance with severity of cervical diseases (OR 4.337; 95% CI 1.248-15.067; P = 0.021), and followed by HPV 16 E-p variant while E-G350 variant only appeared in HSIL and cervical cancer. Our results show that HPV 16 E-p variant was more prevalent than As in Tianjin healthy screening women while As variant was the most frequently type in HSIL and cervical cancer. It is suggested that the mutation of HPV 16 Asian variants, comparing with HPV 16 E-p variants, might contribute to the transformation from HPV 16 persistent infection to cervical cancer.

The Prevalence of HPV Genotypes in Iranian Population: An Update

Background & Objective:

Human papillomavirus (HPV) is the main cause of genital warts and some anogenital cancers in male and female subjects which is commonly transmitted by sexual contacts. The objective of this cross-sectional study was to examine the prevalence of HPV genotypes in 10,266 Iranian male and female population, according to their age.

Methods:

Samples were collected from the penile and anal sites of male subjects and the vagina and cervix of female subjects in a time period between 2011 and 2016. HPV DNA was detected in PCR using the MY09 and MY11 primers, and the INNO-LiPA assay was applied for HPV genotyping. To investigate the relevance of HPV infection and age, the samples were classified into 4 age groups (13-29, 30-44, 45-59, and 60-74).

Results:

Totally, the most common low risk HPV genotypes detected in the studied male and female subjects were HPV-6 (77.7% and 43.3%) and HPV-11 (13.7% and 11.4%), and more frequent high risk HPV genotypes were HPV-16 (5.5% and 16.6%) and HPV-52 (3.2% and 9.6%), respectively. High burden of the HPV infection was observed at ranges of 30 and 44 years (51.8%) with a peak at ranges between 30 and 32 years. No considerable statistically significant correlation was found between HPV infection and age (P=1).

Conclusion:

This study gave an epidemiological overview of circulating HPV genotypes in Iranian population to develop future vaccination policies, though the findings of prevalent HPV genotypes in female subjects were inconsistent with the previous studies reported in Iran.

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.

Genetic variations in E6, E7 and the long control region of human papillomavirus type 16 among patients with cervical lesions in Xinjiang, China

Background

Xinjiang is one of the areas with the highest incidence of cervical cancer in China. Genetic variation in Human papillomavirus type 16 (HPV16) may increase the ability of the virus to mediate carcinogenesis and immune escape, which are risk factors for the progression of cervical cancer. We investigated polymorphism in HPV16 and the distribution of its sub-lineages in the region by analyzing the E6, E7 and long control region (LCR) gene sequences from women with HPV16-positive cervical samples in Xinjiang.

Methods

A total of 138 cases of cervical lesions and squamous cell carcinoma with infection of HPV16 virus were collected. The E6 and E7 genes and LCR of HPV16 virus were sequenced and compared with the HPV16 European prototype reference and other HPV16 mutants for single nucleotide polymorphisms. Neighbor-joining phylogenetic trees were constructed using E6, E7 and LCR sequences.

Results

Fourteen missense mutations were found in the E6 gene; the loci with the highest mutation frequency were T350G (36/75, 48%) and T178G (19/75, 25.3%). In the E7 gene, the locus with the highest mutation frequency was A647G (18/75, 24%). A total of 33 polymorphic sites were found in the LCR, of which T7447C (39/95, 40.1%) was the most frequent.

Conclusion

HPV16 in Xinjiang is mainly of the European variant, followed by the Asian variant type; no Africa 1, 2 or Asia–America variant types were found.

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

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

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.

Association of human papillomavirus type 16 and its genetic variants with cervical lesion in Korea

Persistent human papillomavirus type 16 (HPV16) is the major risk factor for cervical cancer. HPV16 intratypic variants differ in their geographical distribution and oncogenic potential. This study aimed to analyze the distribution of HPV16 variants and their association with cervical lesion histopathology in Korean women. In total, 133 HPV16-positive cervical samples from women admitted to Seoul National University Boramae Hospital were analyzed by sequencing E6, E7, and L1 genes and the long control region (LCR), and the variant distribution according to cervical lesion grade was determined. Isolates were grouped into a phylogenetic lineage, and A1-3, A4, C, and D sublineages were detected in 54.1, 37.8, 0.7, and 7.4% of samples, respectively. The most commonly observed LCR variations were 7521G>A (91.5%), 7730A>C (59.6%), and 7842G>A (59.6%). Furthermore, A4 or D sublineage-positive women had a higher risk for cervical cancer than women who were positive for A1-3. Among HPV phylogenetic clusters, A1-3 was the predominant sublineage, and within A1-3, the 350G polymorphism was highly frequent. These results differed from those of previous studies in Korea and other Asian countries. The findings suggest that cervical neoplasia incidence in HPV16-infected patients could be affected by the distribution of HPV16 variants in the population.

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.

Distribution of human papillomavirus type 16 variants in Lithuanian women with cervical cancer

BACKGROUND AND OBJECTIVE

Cervical cancer usually is caused by HPV 16. However, HPV 16 varies within type; different genotypes are described as prototype or variants. Prevalence of different variants differ according the geographic regions and has an unequal impact for cervical cancer development. Our study aimed to identify which variant of HPV 16 was most prevalent in biological samples taken from Lithuanian women with cervical cancer.

MATERIALS AND METHODS

A total of 122 HPV 16 positive cervical samples (invasive cancer and cervical intraepithelial neoplasia) were investigated and sequenced to identify different variants. HPV 16 was detected using type specific PCR, exact sequence of the virus was obtained by viral DNA sequencing.

RESULTS

Adequate HPV sequence was detected in 106 cases from 122 (86.9% of all cases). After histological confirmation, 96 cases were included in the final analysis. In 33 cases (34.4%) HPV 16 prototype was detected; in 50 cases (52.1%), L83V variant; and in remaining 13 cases (13.5%), multivariant of HPV 16. The frequency of L83V variant in invasive cancer and carcinoma in situ samples was the same (66.7% and 62.0%, respectively; P=0.696). Of analyzed multivariants, 10 were attributed to the European phylogenetic line; 1, to the North American, and 1, to the Asian-American. One sample was not attributed to any of the known phylogenetic lines.

CONCLUSIONS

The European HPV 16 L83V variant is usually associated with high risk of cervical cancer among women. However, statistically significant difference was not achieved when comparing difference of L83V variants between investigated groups and in HPV 16 L83V variant and prototype distribution in CIN3/Ca in situ and cancer.

Lineage distribution and E2 sequence variation of high-risk human papillomavirus types isolated from patients with cervical cancer in Sichuan province, China

To explore the nucleotide sequence variability of the E2 gene in high-risk HPV types in cervical cancer patients from Sichuan province, China, the E2 genes of eight high-risk HPV types were amplified and sequenced. Several novel nucleotide substitutions and deletions were observed. The lineages to which the isolates belonged were determined by phylogenetic analysis, employing the sequence of the representative lineages/sublineages in the coherent classification and nomenclature system as references. This study updates the lineage distribution data on high-risk HPV types in Southwest China and helps broaden understanding of the polymorphism of the E2 gene.

Human Papillomavirus 16 Variants May Be Identified by E6 Gene Analysis

Aims: The aims of the study were (1) to characterize the genetic variability of human papillomavirus (HPV) genotype 16 in the E6 region when this genotype is present in multiple infection samples, (2) to assess the prevalence of variants in our region and (3) to analyze the relationship between variants, patients' ages and pathology. Methods: The Clinical Microbiology and Infection Control Department analyzed samples which were positive for genotype 16 and other genotypes from 2007 to 2013. Variants were assigned to European, Euro-German, Asian, Asian-American or African lineage by sequence analysis. The relationship among variants, age and different types of lesion was studied. Results: In HPV-16 sequence analysis, the European variant was detected in 85.10% of samples, the Asian-American in 7.80%, the African in 4.25% and the Euro-German in 2.83% of specimens. Sequence genetic variability showed 16 nucleotide substitutions. Moreover, non-European variants were mainly found in old women and in isolates from high-grade squamous intraepithelial lesions since European variants were mainly detected in negative cytologies. Conclusion: Multiple infections may take effect on nucleotide substitution and the appearance of recombinant samples. Single gene analysis makes it impossible to detect recombination which has a great influence on drug response and vaccine strategies. Thus, E6 gene analysis would be enough to identify HPV-16 intratypic variants but not to confirm the results.

Prevalence of human papillomavirus variants and genetic diversity in the L1 gene and long control region of HPV16, HPV31, and HPV58 found in North-East Brazil

This study showed the prevalence of human papillomavirus (HPV) variants as well as nucleotide changes within L1 gene and LCR of the HPV16, HPV31, and HPV58 found in cervical lesions of women from North-East Brazil.

Detection and analysis of human papillomavirus (HPV) DNA in breast cancer patients by an effective method of HPV capture

Despite an increase in the number of molecular epidemiological studies conducted in recent years to evaluate the association between human papillomavirus (HPV) and the risk of breast carcinoma, these studies remain inconclusive. Here we aim to detect HPV DNA in various tissues from patients with breast carcinoma using the method of HPV capture combined with massive paralleled sequencing (MPS). To validate the confidence of our methods, 15 cervical cancer samples were tested by PCR and the new method. Results showed that there was 100% consistence between the two methods.DNA from peripheral blood, tumor tissue, adjacent lymph nodes and adjacent normal tissue were collected from seven malignant breast cancer patients, and HPV type 16(HPV16) was detected in 1/7, 1/7, 1/7and 1/7 of patients respectively. Peripheral blood, tumor tissue and adjacent normal tissue were also collected from two patients with benign breast tumor, and 1/2, 2/2 and 2/2 was detected to have HPV16 DNA respectively. MPS metrics including mapping ratio, coverage, depth and SNVs were provided to characterize HPV in samples. The average coverage was 69% and 61.2% for malignant and benign samples respectively. 126 SNVs were identified in all 9 samples. The maximum number of SNVs was located in the gene of E2 and E4 among all samples. Our study not only provided an efficient method to capture HPV DNA, but detected the SNVS, coverage, SNV type and depth. The finding has provided further clue of association between HPV16 and breast cancer.

Genome sequencing accuracy by RCA-seq versus long PCR template cloning and sequencing in identification of human papillomavirus type 58

Background

Genome variations in human papillomaviruses (HPVs) are common and have been widely investigated in the past two decades. HPV genotyping depends on the finding of the viral genome variations in the L1 ORF. Other parts of the viral genome variations have also been implicated as a possible genetic factor in viral pathogenesis and/or oncogenicity.

Results

In this study, the HPV58 genome in cervical lesions was completely sequenced both by rolling-circle amplification of total cell DNA and deep sequencing (RCA-seq) and by long PCR template cloning and sequencing. By comparison of three HPV58 genome sequences decoded from three clinical samples to reference HPV-58, we demonstrated that RCA-seq is much more accurate than long-PCR template cloning and sequencing in decoding HPV58 genome. Three HPV58 genomes decoded by RCA-seq displayed a total of 52 nucleotide substitutions from reference HPV58, which could be verified by long PCR template cloning and sequencing. However, the long PCR template cloning and sequencing led to additional nucleotide substitutions, insertions, and deletions from an authentic HPV58 genome in a clinical sample, which vary from one cloned sequence to another. Because the inherited error-prone nature of Tgo DNA polymerase used in preparation of the long PCR templates of HPV58 genome from the clinical samples, the measurable error rate in incorporation of nucleotide into an elongating DNA template was about 0.149% ±0.038% in our studies.

Conclusions

Since PCR template cloning and sequencing is widely used in identification of single nucleotide polymorphism (SNP), our data indicate that a serious caution should be taken in finding of true SNPs in various genetic studies.

Nucleotide and phylogenetic analysis of human papillomavirus types 6 and 11 isolated from recurrent respiratory papillomatosis in Brazil

There are few studies about the distribution of natural molecular variants of low-risk HPVs. Our aim was to evaluate the E6 early gene variability among HPV-6 and HPV-11 isolates detected in recurrent respiratory papillomatosis (RRP) samples obtained in a cohort of Brazilian patients. We also performed a phylogenetic analysis in order to compare nucleotide sequences identified in our study with previously reported isolates from different anatomic sites (laryngeal papillomas, genital warts, cervical cancer and anal swabs) obtained from other parts of the world to determine the phylogenetic relationships of variants detected in Brazil. The complete coding region of the E6 gene of 25 samples was cloned and sequenced: 18 isolates of HPV-6 (72%) and 7 isolates of HPV-11 (28%). A total of four different HPV-6 genomic variants and two HPV-11 genomic variants was identified. It was not possible to correlate specific variants with disease severity. Phylogenetic trees for both HPV types were constructed enclosing both E6 sequences detected in our study and formerly published sequences. In both phylogenetic trees, the sequences from Brazil did not group together. We could not establish a geographical association between HPV-6 or HPV-11 variants, unlike HPV-16 and HPV-18.

Whole genome analysis of human papillomavirus type 16 multiple infection in cervical cancer patients

The characterization of the whole genome of human papillomavirus type 16 (HPV16) from cervical cancer specimens with multiple infections in comparison with single infection samples as the oncogenic potential of the virus may differ. Cervical carcinoma specimens positive for HPV16 by PCR and INNO-LiPA were randomly selected for whole genome characterization. Two HPV16 single infection and six HPV16 multiple infection specimens were subjected to whole genome analysis by using conserved primers and subsequent sequencing. All HPV16 whole genomes from single infection samples clustered in the European (E) lineage while all multiple infection specimens belonged to the non-European lineage. The variations in nucleotide sequences in E6, E7, E2, L1 and Long control region (LCR) were evaluated. In the E6 region, amino acid changes at L83V were related to increased cancer progression. An amino acid variation N29S within the E7 oncoprotein significantly associated with severity of lesion was also discovered. In all three domains of the E2 gene non synonymous mutations were found. The L1 region showed various mutations which may be related to conformation changes of viral epitopes. Some transcription factor binding sites in the LCR region correlated to virulence were shown on GRE/1, TEF- 1, YY14 and Oct-1. HPV16 European variant prone to single infection may harbor a major variation at L83V which significantly increases the risk for developing cervical carcinoma. HPV16 non-European variants prone to multiple infections may require many polymorphisms to enhance the risk of cervical cancer development.

Whole genome sequencing and evolutionary analysis of human papillomavirus type 16 in central China

Human papillomavirus type 16 plays a critical role in the neoplastic transformation of cervical cancers. Molecular variants of HPV16 existing in different ethnic groups have shown substantial phenotypic differences in pathogenicity, immunogenicity and tumorigenicity. In this study, we sequenced the entire HPV16 genome of 76 isolates originated from Anyang, central China. Phylogenetic analysis of these sequences identified two major variants of HPV16 in the Anyang area, namely the European prototype (E(p)) and the European Asian type (E(As)). These two variants show a high degree of divergence between groups, and the E(p) comprised higher genetic diversity than the E(As). Analysis with two measurements of genetic diversity indicated that viral population size was relatively stable in this area in the past. Codon based likelihood models revealed strong statistical support for adaptive evolution acting on the E6 gene. Bayesian analysis identified several important amino acid positions that may be driving adaptive selection in the HPV 16 population, including R10G, D25E, L83V, and E113D in the E6 gene. We hypothesize that the positive selection at these codons might be a contributing factor responsible for the phenotypic differences in carcinogenesis and immunogenicity among cervical cancers in China based on the potential roles of these molecular variants reported in other studies.

Susceptibility to cervical cancer: an overview

Cervical cancer is the second most common cancer in females worldwide. It is well-established that Human Papillomavirus (HPV) infections play a critical role in the development of cervical cancer. However, a large number of women infected with oncogenic HPV types will never develop cervical cancer. Thus, there are several external environment and genetic factors involved in the progression of a precancerous lesion to invasive cancer. In this review article, we addressed possible susceptible phenotypes to cervical cancer, focusing on host genome and HPV DNA variability, multiple HPV infections, co-infection with other agents, circulating HPV DNA and lifestyle.