- Research Article
- Open Access
Genetic variability in E6 and E7 oncogenes of human papillomavirus Type 16 from Congolese cervical cancer isolates
Infectious Agents and Cancer volume 10, Article number: 15 (2015)
The molecular epidemiological studies showed that some variants of HPV-16, distributed geographically, would present a higher risk of causing cervical cancer. This study aimed to analyze nucleotide changes of HPV-16 E6 and E7 genomic regions from infected Southwestern Congolese women.
DNA of twenty HPV-16 isolates was analyzed by amplifying the E6 and E7 genes using type-specific primers PCR and direct sequencing. The sequences obtained were aligned with the HPV-16 GenBank reference sequences.
Thirteen (65.0%) out of 20 DNA-samples were successfully amplified. Genetic analysis revealed 18 and 4 nucleotide changes in E6 and E7 genomic regions respectively. The most frequently observed nucleotide variations were the missense C143G, G145T and C335T in E6 (100%), leading to the non-synonymous amino acid variation Q14D and H78Y. E7 genomic region was found to be highly conserved with two most common T789C and T795G (100%) silent variations. All HPV-16 variants identified belonged to the African lineage: 7 (53.8%) belonged to Af-1 lineage and 6 (46.1%) to Af-2 lineage. The missense mutation G622A (D21N) in the E7 region seems to be described for the first time in this study.
This study reported for the first time the distribution of HPV-16 E6 and E7 genetic variants in infected women from southwest Congo. The findings confirmed almost ascendancy of the African lineage in our study population.
The Republic of Congo has one of the highest rates of cervical cancer incidence in Sub-Sahara Africa. According to the World Health Organization (WHO), the reported age-standardized incidence rate and age-standardized mortality rate were 25.6 and 13.0 cases per 100 000 women in 2014 .
Epidemiological and molecular studies are showed that High-risk human papillomaviruses (HPV) types are the etiological agents of cervical cancer [2,3]. HPV-16 is the most common high-risk genotype involved in cervical cancer, representing for over 50% of all cases .
However, most HPV infections resolve spontaneously, suggesting that other co-factors are necessary in the persistence and lesion progression to cancer. These cofactors include some variants of high-risk types such as those of HPV16. Several molecular variants were described and classified according to geographical regions for the HPV-16 genotype [5,6]. Conventionally, a variant is defined when there is a difference of less than 2% between the complete genomes of the same HPV type. A difference from 1.0% has been used to define the lineage and the sub-lineages variants by a difference from 0.5 to 1.0% between the lineages [7-9]. The variants have been classified in six major phylogenetic lineages (European, E; Asian, As; African-1, Af-1; African-2, Af-2; Asian-American, AA and North-American, NA) showing more than 98% nucleotide level similarity in L1 or E6 genes [10,11]. However, a new phylogenetic nomenclature was proposed for the HPV-16 in four lineages as follows: A (previously called European-Asian, EAS), B (African 1, Af-1), C (African 2, Af-2), and D (North-American/Asian-American, NA/AA) [7,9].
HPV-16 variants studies have shown that certain intratypic-variants would contribute more than others to the persistence of HPV infection and cervical cancer progression [12,13]. In particular, non-European variants are more highly oncogenic them European variants [12,14,15]. Indeed, the importance of the E6 and E7 genes lies in their ability to transform and immortalize the host cell. A single nucleotide change in these genes could affect this ability, which would explain the difference in some variants to progress quickly them others to cervical cancer . E6 and E7 interacts with several cellular proteins in particular, they inactivate the p53 and the pRb respectively, which are two cancer suppressor major proteins .
In this perspective and with the advent of prophylactic vaccines against HPV-16/18 and cervical cancer, it becomes necessary to know the HPV genetic variants present in each population.
In Congo, our previous study showed that the most common genotype in invasive cervical cancer and its precancerous lesions was HPV-16 . Thus, to better formulate immunization strategies and other preventive measures in Congo, this preliminary study aimed to analyze the nucleotide changes of the HPV-16 E6 and E7 isolates in order to identify the different genetic variants prevalent among infected women with cervical carcinoma in southwest of the Congo.
Out of 20 HPV-16 single infected cervical cancer samples examined in this study, the E6 and E7 genes were able to be amplified and sequenced in 13 (13/20, 65.0%) DNA samples. Nucleotide (nt) changes were analyzed between nt 104 to 559 for E6 region and nt 562 to 858 for E7 region.
E6 nucleotide sequences analysis
By comparison with the HPV-16 reference sequence (European prototype, NC_001526), mutation analysis showed that 100% (13/13) of the analyzed sequences contained at least one nucleotide change in the E6 region. Eighteen single nucleotide changes were observed in HPV-16 E6 sequences, of which 13/18 (72.2%) were non-synonymous mutations (all missense) and 5/18 (27.8%) were synonymous mutations (silent). The most frequently observed variations were C143G (Q14D), G145T (Q14D), T286A, A289G and C335T (H78Y), which was found in all samples. All nucleotide changes, variants and their prevalence were reported in Table 1.
E7 nucleotide sequences analysis
The details of E7 mutations are summarized in Table 1. E7 gene was found to be highly conserved in all cases compared to E6 region. Only four single nucleotide changes were observed in HPV-16 E7 sequences, with 2/4 (50.0%) missense and 2/4 (50.0%) silent. The most frequently observed variations were the common silent mutations T789C and T795G at codons 76 and 78, which was found in all E7 sequences. However, the common non-synonymous mutation A to G transition at nt 647 and causing amino acid change N29S, was found in 3 (23.0%) out of 13 sequences. An undescribed mutation was found in three sequences in this study at nt 622, that was generated a non-synonymous nucleotide change G622A (D21N).
HPV-16 variants and phylogenetic analysis
All HPV-16 positive sequences analyzed belonged to the African lineage, of which 7 (53.8%) out of 13 belonged to Africain-1 lineage (Af-1) and 6 (46.1%) out of 13 to Africain-2 lineage (Af-2). All African lines showed a common pattern of five characteristic mutations in E6, namely, C143G, G145T, T286A, A289G, and C335T, which lead to two non-synonymous amino acid changes Q14D and H78Y. In E7, all HPV-16 isolates showed a common pattern of two silent mutations namely T789C and T795G at codons 76 and 78 respectively. No European (E), Asian (As), Asian-American (AA) or northern American (NA) lineages were detected in this study. Phylogenetic trees were respectively built from eight sequences E6 regions (Figure 1) and 5 of the E7 region (Figure 2) from this study and other HPV-16 published sequences available in GenBank.
The HPV-16 genotype was recently identified as the most prevalent genotype in cervical cancer cases and high grade precancerous lesions in the southwest of Congo .
To date, there are no molecular epidemiological studies on HPV-16 variants carried out in Congo. In the present study, we characterized the nucleotide changes of the HPV16 mono-infection isolates among infected Congolese women with cervical cancer based on the E6 and E7 viral oncogenes analysis.
Several studies showed that HPV-16 variants may influence viral persistence and progression to cervical cancer, which partly explains the high prevalence of this cancer in some populations compared to other . Also, their distribution is geographically related and associated with the ethnic group [19,20]. The HPV 16 E6 and E7 early genes are essential in cervical cancer pathogenesis . Any change in their sequences could induce an alteration of the biological function of the proteins encoded by these genes, which would justify the importance of these kinds of studies, especially in the poorly studied populations as ours.
In this study, phylogenetic analysis of HPV16 E6 and E7 genes was performed to analyze the local distribution of HPV16 lineages, according to the Yamada et al.  classification. A total of eighteen point mutations in E6 and four in E7 genomic regions were found. Most nucleotide changes reported in our study have been previously described. At least a nucleotide change was observed in all sequences analyzed. This result is consistent with the literature data indicating that more than 90% of the E6 sequences contains mutations in cervical carcinomas, whereas the E7 gene seems to be more conserved [5,6,22,23].
At the E6 genomic region, five characteristic mutations, C143G, G145T, T286A, A289G, and C335T defining the African lineage were found in our study as well as it’s was previously described [8,11,24]. However, some of the observed mutations are significant in the literature; in particular the two non-synonymous changes African lineage located in E6 N-terminus coding region at codon 10 and 14, which lead to amino acid changes R10T and Q14D respectively. These amino acid changes may lead to a difference in affinity with p53 and significantly altered its degradation rate by this oncogene [25,26]. Crook et al.  showed that the amino acid mutations in these positions increased by 180% the affinity of E6 with p53. It is also reported that the amino acid change R10I in addition to other mutations described between nucleotide position 106 and 113 could lead to a low affinity for binding of E6 and degrade p53 , although these findings remain controversial and poorly documented.
In contrast to E6, the E7 gene variability has been less studied. Our studies on E7 showed a highly conserved region of the HPV genome with some mutations (mainly point mutations) described and few lead to amino acid changes [23,28]. Out of four mutations identified in our study, two were the common silent variations found in all sequences. The two others were non-synonymous mutations, of which a common described N29S was found in three isolates (23%). This mutation often described in Asia, was located in the cervical carcinoma in Africa [23,29,30]. The importance of this amino acid change was suggested by Zehbe et al.  because of its location in an immunoreactive region. However, E7 protein with the N29S mutation present an identical power transforming has the prototype strain, after measuring the cooperation E7/activated ras gene in rat embryo fibroblast cells . In addition, N29S is involved in binding to pRB and could alter the affinity of the E7 protein for pRB, or modification of the oncogenic potential. But it was also proved that the variant A647G only binds with similar affinity to that of the prototype [32,33]. The D21N mutation (G622A) described in this study seems to be for the first time. Only one case showed this variation; future studies by sequencing after cloning in a vector are needed to confirm its description in our population.
All HPV-16 isolates examined in this study belonged to the African lineage, of which 53.8% belonged to Af-1 lineage and 46.1% to Af-2 lineage. Our results are in agreement with the literature. Indeed, the pioneering studies of Yamada and Wheeler showed that the African lineage was dominated in the continent. More than 90.7% of the variants in Africa region are African lineage, of which 61% belongs to the Af-1 lineage and 29.7% to the Af-2 lineage [6,11].
Tu et al.  in South Africa, Qmichou et al.  in Morocco and Buonaguro et al.  in Uganda reported have identified in their respective studies 41%, 64.5% and 100% of variants belonging to the African lineage. These results substantially similar show the importance of taking into account the African variants in the global effort to fight against cervical cancer. The findings obtained by Qmichou et al.  which reported 35.5% of variants belonged to the Af-1 lineage and 29% to the Af-2 lineage such as Yamada (61% Af-1 and 29.7% Af-2) are consistent with ours by showing that the Af-1 lineage is predominant on the continent . However, non-African lineages were also reported in the continent, especially in North Africa and South Africa regions. The introduction of these lineages could be explained by the diversity of their population with European populations. In our study, the absence of non-African lineages could be explained by the small size of our sample, but also by the fact that their frequency is very low in a homogeneous African black population. However, further investigation seems to be necessary to confirm this hypothesis.
Knowing that HPV variants data are important in developing HPV diagnostics, vaccines, and other therapeutic approaches to monitoring virus-induced diseases , the major limitation of this study was the small size of DNA-sample analyzed. However, this is the first study that provides information on HPV-16 E6 and E7 genetic diversity in Congolese women. A larger sample size of HPV 16 including E6 and E7, but also L1 and LCR genes will be necessary to better identify all the genomic changes and appreciate their impact on the various grades of cervical lesions.
The present study provides preliminary data on HPV16 E6 and E7 variants in the southwest of Congo. Findings confirmed the predominance of the African lineage in our samples, which could guide control strategies against cervical cancer in Congo. However, future studies with a large sample are needed to better assess the distribution of these variants depending on the grade and the risk of developing cervical cancer.
Materials and methods
DNA Biological samples
DNA of twenty HPV-16 single infection isolates of Congolese women with cervical cancer from our laboratory DNA bank (Laboratory of Virology, Microbiology and Quality/ETB, Faculty of Sciences and Techniques de Mohammedia, UH2C, Morocco) constituted from our previous study  was investigated. In brief, All the samples of biopsies were squamous cell carcinoma (SCC) selected based on the availability of the DNA sample. Samples were collected from the patients attending the pathology laboratory of the General Hospital of Loandjili (Pointe-Noire, Southwest Congo).
PCR amplification with type-specific (TS) primers and sequencing
Amplification of HPV-16 E6 and E7 genes was performed using TS primers flanking the encoding regions of HPV16 E6 ORF (nt: 41–576) and HPV16 E7 ORF (nt: 483–911), as previously described . PCR was performed with some minor modifications in 25 μL of reaction mixture containing 1X PCR buffer, 2.5 mM MgCl2, 0.8 mM of each dNTP, 0.4 μM/μL of sense and antisense primers, 0.025U/μl of Go Taq DNA polymerase (Promega) and 1 μL of template DNA. The thermal program (Perkin Elmer 2400 GeneAmp® PCR thermal Cycler, Scientific Support, Inc, Hayward, CA) started with a pre-heat of 94°C for 10 min, followed by 35 cycles according to the protocol: 45 sec denaturation at 94°C, 45 sec annealing at 49.5°C for E6 and 53.5°C for E7 and 45 sec extension at 72°C, with 7 min final elongation at 72°C. For each reaction, DNA from SiHa cell lines as positive controls and Ultra-pure PCR water nuclease-free (Bioline, UK) as negative PCR control were included.
PCR products were analyzed based on the expected specific bands in 2% agarose gel. Type-specific PCR products were directly sequenced using the “Big Dye Terminator v3.1 Cycle Sequencing kit” (Applied Biosystems) according to the manufacturer’s instructions in molecular and functional genomics platform (UATRS-CNRST, Rabat, Morocco).
Phylogenetic analysis of sequences
The sequences were subsequently analyzed by NCBI using online BLAST 2.0 software server (http://www.ncbi.nih.gov/BLAST/) to confirm genotype and amplified region. Multiple sequence alignments of the full-length of E6 and E7 genes were performed using BioEdit Sequence Alignment Software v220.127.116.11 and Clustal W program . HPV-16 DNA nucleotide positions were numbered according to the HPV-16 reference (European prototype published in GenBank with accession number: NC_001526.2) .
Phylogenetic trees were built by MEGA 5.1 package  using a distance-based criterion, Tamura 3-parameter as the substitution model  and Maximum Likelihood algorithm with bootstrap proportions were calculated with 1000 replicates to test the robustness of the major phylogenetic groups . The reference viral strains used for the construction of phylogenetic trees were obtained from the NCBI GenBank Database, which belong to Asian-American lineage (GenBank accession numbers: AF402678; HQ644247), African-1 (GenBank: AF472508; HQ644238), African-2 (GenBank: AF472509; HQ644237) lineage, and European lineage (GenBank: AF536179; HQ644236).
GenBank accession numbers
The representative nucleotide sequences from each region sequenced in this study were deposited into GenBank Database (http://www.ncbi.nlm.nih.gov/GenBank), under the accession numbers KP313765 to KP313772 for E6 and KP313773 to KP313777 for E7 genes.
During the collection of samples that have been the subject of this study, informed consent was obtained for all women and the study protocol was approved by the local ethical review board of research in health Sciences of Congo ("Comité d’éthique de la Recherche en Sciences de la Santé").
WHO/ICO: information centre on HPV and cervical cancer (HPV information center). Human Papillomavirus and related Deseases in Congo. Summary Report [17 June 2014] 2014, http://www.hpvcentre.net/statistics/reports/COG.pdf.
Zur Hausen H. Human papillomaviruses in the pathogenesis of anogenital cancer. Virology. 1991;184:9–13.
Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, et al. Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol. 1999;189:12–9.
de Sanjose S, Quint WGV, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, et al. Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study. Lancet Oncol. 2010;11:1048–56.
Yamada T, Wheeler CM, Halpern AL, Stewart AC, Hildesheim A, Jenison SA. Human papillomavirus type 16 variant lineages in United States populations characterized by nucleotide sequence analysis of the E6, L2, and L1 coding segments. J Virol. 1995;69:7743–53.
Wheeler CM, Yamada T, Hildesheim A, Jenison SA. Human papillomavirus type 16 sequence variants: identification by E6 and L1 lineage-specific hybridization. J Clin Microbiol. 1997;35:11–9.
Pérez S, Cid A, Iñarrea A, Pato M, Lamas MJ, Couso B, et al. Prevalence of HPV 16 and HPV 18 Lineages in Galicia. Spain PLoS ONE. 2014;9, e104678.
Cornet I, Gheit T, Franceschi S, Vignat J, Burk RD, Sylla BS, et al. Human Papillomavirus Type 16 Genetic Variants: Phylogeny and Classification Based on E6 and LCR. J Virol. 2012;86:6855–61.
Burk RD, Harari A, Chen Z. Human papillomavirus genome variants. Virology. 2013;445:232–43.
Cornet I, Gheit T, Iannacone MR, Vignat J, Sylla BS, Del Mistro A, et al. HPV16 genetic variation and the development of cervical cancer worldwide. Br J Cancer. 2013;108:240–4.
Yamada T, Manos MM, Peto J, Greer CE, Munoz N, Bosch FX, et al. Human papillomavirus type 16 sequence variation in cervical cancers: a worldwide perspective. J Virol Methods. 1997;71:2463–72.
Bernard HU, Calleja-Macias IE, Dunn ST. Genome variation of human papillomavirus types: phylogenetic and medical implications. Int J Cancer. 2006;118:1071–6.
de Villiers EM, Fauquet C, Broker TR, Bernard HU HU, zur Hausen H. Classification of papillomaviruses. Virology. 2004;324:17–27.
Freitas LB, Chen Z, Muqui EF, Boldrini NAT, Miranda AE, Spano LC, et al. Human Papillomavirus 16 Non-European variants are preferentially associated with high-grade cervical lesions. PLoS One. 2014;9, e100746.
Sichero L, Ferreira S, Trottier H, Duarte-Franco E, Ferenczy A, Franco EL, et al. High grade cervical lesions are caused preferentially by non-European variants of HPVs 16 and 18. Int J Cancer. 2007;120:1763–8.
Kim MK, Kim HS, Kim SH, Oh JM, Han JY, Lim JM, et al. Human papillomavirus type 16 E5 oncoprotein as a new target for cervical cancer treatment. Biochem Pharmacol. 2010;80:1930–5.
Lee K, Magalhaes I, Clavel C, Briolat J, Birembaut P, Tommasino M, et al. Human papillomavirus 16 E6, L1, L2 and E2 gene variants in cervical lesion progression. Virus Res. 2008;131:106–10.
Boumba A, Hilali L, Mouallif M, Moukassa D, Ennaji M. Specific genotypes of human papillomavirus in 125 high-grade squamous lesions and invasive cervical cancer cases from Congolese women. BMC Public Health. 2014;14:1320.
de Araujo Souza PS, Sichero L, Maciag PC. HPV variants and HLA polymorphisms: the role of variability on the risk of cervical cancer. Future Oncol. 2009;5:359–70.
Villa LL, Sichero L, Rahal P, Caballero O, Ferenczy A, Rohan T, et al. Molecular variants of human papillomavirus types 16 and 18 preferentially associated with cervical neoplasia. J Gen Virol. 2000;81:2959–68.
Psyrri A, DiMaio D. Human papillomavirus in cervical and head-and-neck cancer. Nat Clin Prac Oncol. 2008;5:24–31.
Zehbe I, Wilander E, Delius H, Tommasino M. Human papillomavirus 16 E6 variants are more prevalent in invasive cervical carcinoma than prototype. Cancer Res. 1998;58:829–33.
Eschle D, Durst M, ter Meulen J, Luande J, Eberhardt HC, Pawlita M, et al. Geographical dependence of sequence variation in the E7 gene of human papillomavirus type 16. J Gen Virol. 1992;73(Pt 7):1829–32.
Huertas-Salgado A, Martin-Gamez DC, Moreno P, Murillo R, Bravo MM, Villa L, et al. E6 molecular variants of human papillomavirus (HPV) type 16: an updated and unified criterion for clustering and nomenclature. Virology. 2011;410:201–15.
Stöppler MC, Ching K, Stöppler H, Clancy K, Schlegel R, Icenogle J. Natural variants of the human papillomavirus type 16 E6 protein differ in their abilities to alter keratinocyte differentiation and to induce p53 degradation. J Virol. 1996;70:6987–93.
Ellis JRM, Etherington I, Galloway D, Luesley D, Young LS. Antibody responses to HPV16 virus-like particles in women with cervical intraepithelial neoplasia infected with a variant HPV16. The Lancet. 1997;349:1069–70.
Crook T, Vousden KH, Tidy JA. Degradation of p53 can be targeted by HPV E6 sequences distinct from those required for p53 binding and trans-activation. Cell. 1991;67:547–56.
Pande S, Jain N, Prusty BK, Bhambhani S, Gupta S, Sharma R, et al. Human papillomavirus type 16 variant analysis of E6, E7, and L1 genes and long control region in biopsy samples from cervical cancer patients in north India. J Clin Microbiol. 2008;46:1060–6.
Song YS, Kee SH, Kim JW, Park NH, Kang SB, Chang WH, et al. Major sequence variants in E7 gene of human papillomavirus type 16 from cervical cancerous and noncancerous lesions of Korean women. Gynecol Oncol. 1997;66:275–81.
Qmichou Z, Ennaji MM, Amrani M, El Fahime El M, Meloul M, Meftah EL, et al. Molecular Characterization of HPV16 E6 and E7 Variants among Women with Cervical Cancer in Moroco. Britsh Microbiol Res J. 2013;3:692–705.
Fujinaga Y, Okazawa K, Nishikawa A, Yamakawa Y, Fukushima M, Kato I, et al. Sequence variation of human papillomavirus type 16 E7 in preinvasive and invasive cervical neoplasias. Virus Genes. 1994;9:85–92.
Chow VT, Loh E, Yeo WM, Tan SY, Chan R. Identification of multiple genital HPV types and sequence variants by consensus and nested type-specific PCR coupled with cycle sequencing. Pathology. 2000;32:204–8.
Wu Y, Chen Y, Li L, Yu G, He Y, Zhang Y. Analysis of mutations in the E6/E7 oncogenes and L1 gene of human papillomavirus 16 cervical cancer isolates from China. J Gen Virol. 2006;87:1181–8.
Tu JJ, Kuhn L, Denny L, Beattie KJ, Lorincz A, Wright TC. Molecular variants of human papillomavirus type 16 and risk for cervical neoplasia in South Africa. Int J Gynecol Cancer. 2006;16:736–42.
Buonaguro FM, Tornesello ML, Salatiello I, Okong P, Buonaguro L, Beth-Giraldo E, et al. The uganda study on HPV variants and genital cancers. J Clin Virol. 2000;19:31–41.
Garbuglia AR, Carletti F, Minosse C, Piselli P, Zaniratti MS, Serraino D, et al. Genetic variability in E6 and E7 genes of human papillomavirus -16, -18, -31 and -33 from HIV-1-positive women in Italy. New Microbiol. 2007;30:377–82.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Seedorf K, Krammer G, Durst M, Suhai S, Rowekamp WG. Human papillomavirus type 16 DNA sequence. Virology. 1985;145:181–5.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Tamura K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G + C-content biases. Mol Biol Evol. 1992;9:678–87.
Saitou N, Nei M. The neighbor-joiningmethod: a newmethod for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Authors gratefully acknowledge the Direction of General Hospital of Loandjili at Pointe-noire, Congo for having facilitated the transportation of samples to Morocco. We thank the functional genomics platform UATRS-CNRST Rabat, Morocco for sequencing. We also thank Dr. Siham Fellali for his critical reading of this article. This project was financially supported by the Moroccan Minister of Higher Education.
The authors declare that they have no competing interests.
LMAB was responsible for design and the experiment of the study. SZA and JVM participated in drafting the manuscript. DM and LMAB have conducted the recruitment process of the sample. LH and MME were responsible for project implementation and participated in the critical reading of the manuscript. All authors read and approved the final manuscript.