Skip to main content
Download PDF

Unique LCR variations among lineages of HPV16, 18 and 45 isolates from women with normal cervical cytology in Ghana

Virology Journal volume 14, Article number: 85 (2017) Cite this article

Abstract

Background

In addition to being useful for classification, sequence variations of human Papillomavirus (HPV) genotypes have been implicated in differential oncogenic potential and a differential association with the different histological forms of invasive cervical cancer. These associations have also been indicated for HPV genotype lineages and sub-lineages. In order to better understand the potential implications of lineage variation in the occurrence of cervical cancers in Ghana, we studied the lineages of the three most prevalent HPV genotypes among women with normal cytology as baseline to further studies.

Methods

Of previously collected self- and health personnel-collected cervical specimen, 54, which were positive for HPV16, 18 and 45, were selected and the long control region (LCR) of each HPV genotype was separately amplified by a nested PCR. DNA sequences of 41 isolates obtained with the forward and reverse primers by Sanger sequencing were analysed.

Results

Nucleotide sequence variations of the HPV16 genotypes were observed at 30 positions within the LCR (7460 – 7840). Of these, 19 were the known variations for the lineages B and C (African lineages), while the other 11 positions had variations unique to the HPV16 isolates of this study. For the HPV18 isolates, the variations were at 35 positions, 22 of which were known variations of Africa lineages and the other 13 were unique variations observed for the isolates obtained in this study (at positions 7799 and 7813). HPV45 isolates had variations at 35 positions and 2 (positions 7114 and 97) were unique to the isolates of this study.

Conclusion

This study provides the first data on the lineages of HPV 16, 18 and 45 isolates from Ghana. Although the study did not obtain full genome sequence data for a comprehensive comparison with known lineages, these genotypes were predominately of the Africa lineages and had some unique sequence variations at positions that suggest potential oncogenic implications. These data will be useful for comparison with lineages of these genotypes from women with cervical lesion and all the forms of invasive cervical cancers.

Background

With an estimated 3052 new cases of cervical cancer in Ghana in the year 2012, cervical cancer is one of the two leading cause of cancer incidence and mortality in Ghana and in most sub-Saharan African countries. Cervical cancer is a significant health concern in Ghana in light of the facts that in spite of its high burden, cervical cancer screening programmes are still lacking. In the absence of empirical cancer registry based data, the estimated age-standardized cervical cancer incidence rate for Ghana (35.4 per 100,000) indicates a possible difference in the cervical cancer burden among close neighbours; Benin (27.6 per 100,000), Togo (21.5 per 100,000), Burkina Faso (23.3 per 100,000) and Cote d’Ivoire (21.7 per 100,000) [1].

These geographical variations in the incidence of cervical cancer are also seen in the prevalence of its etiological agent, which is the human Papillomavirus (HPV). In spite of these variations, HPV genotypes 16, 18 and 45 together are responsible for more than 80% of cervical cancer globally [2,3,4]. Molecular studies of each of these genotypes have shown sequence variations in some genes that have been the basis for their classification into phylogenetic linages and sub-lineages [5]. Conventionally, if the DNA sequence of the L1 gene of two HPV isolates defer by between 1 and 10%, then these two HPV isolates are variants (lineage) of a genotype. A difference of between 0.5 and 1% makes the HPVs sub-lineage (sub-variant) of a lineage [6, 7]. Furthermore, the level of nucleotide sequence identities or differences in other parts of the HPV genome (including the LCR, E6 and E7) and mostly in the full genome sequence have resulted in similar classification [5,6,7,8,9,10,11].

Specifically, variants of HPV16 have been classified most recently, as 4 phylogenetic lineages based on LCR and E6 open reading frames of its genome. These are, a) lineage A (also known as European-Asian, EAS lineage), consisting of at least the sub-lineages A1, A2, A3, (known previously as European, E lineage) and A4 (previously known as Asian, As lineage); b) Lineage B (known previously as African-1 (Afr-1) lineage), consisting of at least two sub-lineages B1 (Afr-1a) and B2 (Afr-1b); c) Lineage C (known previously as African-2 (Afr-2) lineage), previously consisting of the sub-lineages Afr-2a and Afr-2b, and d) lineage D (Asian-American/North American, AA/NA lineage), consisting of at least the sub-lineages, D1 (known previously as North-American (NA) lineage), D2 (known previously as Asian-American-1 (AA1) lineage), and D3 (Asian-American-2 (AA2) lineage). These groups of lineages have more than 95% sequence homology in the long control region (LCR) [5, 7, 12,13,14].

For HPV18 variants, two phylogenetic lineages (previously known as non-African and African) were identified based on full genome sequence analysis in an older classification. The non-African lineages had two sub-lineages namely, European (E) and Asian-American (AA) while the African lineage also had two sub-lineages, namely African 1 (Afr-1) and African-2 (Afr-2) [15]. However, with a more recent classification, three lineages have been identified; the previous non-Africa lineage was identified and designated as lineage A with five sub-lineages, designated A1 to A5. Part of the previous African lineage was identified as a separate lineage and designated lineage B with three sub-lineages, designated B1 to B3. The rest of the previous African lineage were identified as a new lineage and designated lineage C [6, 10, 15]. Also based on full genome sequence analysis, HPV 45 variants have been classified into two phylogenetic lineages, namely A and B. while lineage A has at least three sub-linages; A1, A2 and A3, lineage B has two sub-lineages, B1 and B2 [6, 8, 15].

The long control region (LCR) of the HPV genome is a 850 bp regulatory region that includes a section for regulating the transcription of HPV genes through interactions with many cellular and viral factors [13]. Therefore, sequence variation in the LCR in addition to being useful for classification has been implicated in the different extent of persistence and progression and oncogenic potential of HPV variants [15, 16]. It has also been shown that HPV genotypes (which vary in their LCR sequences) are differentially associated with different histological forms (squamous cell carcinoma (SCC); adenocarcinoma (ADC) and adenosquamous carcinoma (ADSC)) of invasive cervical cancer [17,18,19,20,21,22]. Specifically for invasive cervical cancer (ICC) reported in Ghana, HPV16 and closely related HPV types such HPV31 were common in SCC, while HPV18 and closely related types such HPV45 and HPV59 were common in ADC and ADSC [23, 24]. Such differential specificities have also been indicated for HPV variants/lineages [5, 14, 25,26,27,28,29].

It is therefore important for our population and in order to enhance cervical cancer screening, to study the variants of the reported most prevalent HPV genotypes in Ghana, to better understand the potential of developing cervical cancers. Although there are limited empirical data on the prevalence of human Papillomavirus (HPV) in the general population, the most prevalent HPV genotypes in Ghana are HPV 16, 18 59, and 45. Few hospital based studies on HPV and its distribution in cancers have reported, for instance, a 10% HPV prevalence among 75 women screened for cervical cancer in a tertiary hospital [30]. Three other studies, using formalin-fixed paraffin embedded (FFPE) cervical tissue samples of women diagnosed with high grade cervical lesion or invasive cervical cancer (at the same hospital) reported HPV prevalence of 93.9% [31], 98.0% [32] and 89.8% [23, 24]. Overall, only two of these studies reported the genotype specific prevalence, which indicated that HPV 18, 59, 16 and 45 were among the common genotypes in Ghana [23, 32]. However, there are no data on the variants of these common HPV genotypes in Ghana. With increasing evidence of HPV genotype, subtype and lineage specificity and association with different forms of invasive cervical cancer, subtype and lineage determination may become important in risk assessment. Additionally, kits development may depend on sequence data of HPV variants for specific targeting, particularly for the HPV 16, 18 and 45. Therefore, our study, focused on the determination of the lineages of HPV 16, 18 and 45 isolated from women participating in a community-based cervical cancer screening activity, all of whom were negative by Pap smear testing (normal cytology results).

Methods

Cervical specimen

Previously collected cervical specimen of women with normal cytology results, during a community based cervical cancer screening programme, were used for this study. The specimen positive for HPV16, 18 and 45 (either as a single infection or in a multiple infection) were selected and the long control region (LCR) of each HPV was separately amplified by a nested PCR to confirm the genotype present in the specimen. A total of 20 HPV16 positive specimens, 15 HPV18 positive specimens and 18 HPV45 positive specimens were used.

DNA extraction

DNA was extracted from these specimen on a MagNA Pure LC automated system (Roche Molecular Systems, Inc, Pleasanton, USA) using 200 μL of the cell suspension and the MagNA Pure DNA Extraction Kit (Roche Molecular Systems, Inc, Pleasanton, USA) according to the manufactures specification with a modification of an addition of 15 μL of RNase to the lysis buffer for each specimen. The negative controls included for each set of extraction were PBS, DNase/RNase free water, and DNAGard, while the positive control was a suspension of SiHa and HeLa cells with an integrated HPV16 and HPV18 genome respectively. The specimen/DNA quality was assessed by the amplification of the human RNase H gene in 10 μL of DNA extract and analysed with Light Cycler 480 and related software (Roche Molecular Systems, Inc, Pleasanton, USA).

Nested PCR amplification of HPV16, 18 and 45 LCR

Using 5 μL of DNA extracted from each of the selected specimen, first round PCRs were performed with the primer sets as indicated in Table 1 for HPV16, 18 and 45. The 50 μL final reaction volume also contained, 400 μM of each dNTP, 1X PCR buffer (200 mM Tris–HCl pH 8.4; 500 mM KCl), 2.5 mM of MgCl2, 1.25 μM of each of the two primers, and 5U of AmpliTaq DNA polymerase. The amplifications were achieved with the following conditions; an initial denaturation step at 95 °C for 5 min (10 min for HPV45) was followed by 60 cycles (40 cycles for HPV45) of denaturation at 95 °C for 30 s, primer annealing at 65 °C (55 °C for HPV45) for 30 s and an elongation at 72 °C for 30 s followed with a 7 min final extension at 72 °C and holding the temperature at 4 °C. Extracted DNA form Hela cells, SiHa cells with integrated HPV18, HPV16 genome respectively and confirmed clinical isolate of HPV45 were used as positive control and additional PCR water negative controls were used in each set of amplification.

Table 1 List of primer sequences use for first and second round PCR for each og the three HPV genotypes
Full size table

The second round PCR used 5 μL of the first round product with the second round primer sets indicated in Table 1 for the amplifications. The 100 μL final reaction volume contained 200 μM of each dNTP, 1 X PCR buffer (200 mM Tris–HCl pH 8.4; 500 mM KCl), 2.5 mM of MgCl2, 1.25 μM of each of the two primers, and 6.3 U of Taq DNA polymerase. The amplifications were achieved with the same cycling conditions as that of the first round PCR.

Visualization of second round PCR product

For all the second round PCR products of HPV16, 18 and 45, 20 μL were resolved on a 1.5% agarose gel in 1X TBE buffer with 1 μg/mL ethidium bromide at 120 V for an hour. The gels were documented with GelDoc™ XR+ System (Bio-Rad Laboratory, Mississauga, Canada).

Sequencing of second round PCR product

The second round PCR products were purified with the Amicon Ultra centrifuge filers (Millipore Ireland Ltd, Tullagreen, Ireland) with slight modification made to the manufactures instruction. Specifically, a filter column was place upright, into a micro-centrifuge tube for each sample, after which 200 μL of PCR water and 80 μL of a PCR product were added. The assembly was spun at 14000xg for 12 min. Additional 100 μL of PCR water were added to each filter column. The filter columns were thereafter transferred, inverted, to a new set of micro-centrifuge tubes and spun for 2 min at 1000xg. The columns were discarded and the eluted DNA in each tube was quantified at 230 nm using NanoDrop ND-100 spectrophotmeter (Thermo Fisher Scientific Inc, Wilmington, USA). The eluted DNA was diluted to 50 ng/μL and sequenced with ABI Applied Biosystems 3730xl DNA Analyzer BigDye® Terminator v 3.0 (Applied Biosystem, Canada). DNA Sequences obtained with the forward and reverse primers by Sanger sequencing were assembled with SeqMan Pro software (DNAStar). The sequences obtained for HPV16, HPV18 and HPV45 PCR products (listed in Additional file 1: Table S1) were aligned to HPV16, HPV18 and HPV45 reference sequences by the Cluster W method, using MEGA 5.2 software. Phylogenetic trees were constructed based on the Test Neighbour Joining using MEGA 6 software.

Ethical approval and consent to participants

The study that collected the specimen from which the HPV isolates were obtained received ethical approval form the Ethics Review Committee of the Ghana Health Service (ID No GHS-ERC: 06/11/10) and worked within the guidelines of the Ethical Review Committee and conducted the study in accordance to the ethical standards as declared in Declaration of Helsinki and the Belmont Report. Assurances were given and efforts were made to protect the privacy and confidentiality of the participants enrolled and their collected data. As such personal/demographic data as well as data that will be useful for identifying the participants were not included in this study.

Results

The obtained partial LCR nucleotide sequences for each of the HPV genotypes were of an average size of 375 bp. These have been deposited in NCBI GenBank and the accession number shown in the Additional file 1: Table S1. The obtained sequences were aligned with the LCR positions of the complete genome of the HPV16 prototype NC0015626.2 [33, 34], the HPV18 prototype AY262282_A1 [35] and the HPV45 prototype X74479.1 [36]. The alignments also included selected HPV variants with very similar sequences obtained by an NCBI BLAST of the isolates’ sequences. A neighbour joining phylogenetic tree was constructed for each HPV type (Figs. 1–3). All 20 HPV16 specimen were confirmed by positive results for the LCR HPV16 type specific PCR. Nucleotide sequence variations were observed at 30 positions within the LCR (Table 2). Of these 30 positions, 19 were the known variations often observed for African variants, while the other 11 positions had variations unique to the HPV16 isolates of this study (Table 2). The phylogenetic tree showed that about 85% of the HPV16 isolates were of the lineage C, 5% of lineage B, both of which were previously known as African variants (Fig. 1) and the lineages of two isolates (10%), isolate-866 and isolate 618 were not clearly identified, although they were most likely to belong to lineage B.

Fig. 1
figure 1

An evolutionary relationship of HPV16 isolates based on sequences between the genome positions 7469 and 7840 within the long control region (LCR.). The evolutionary history was inferred using the Neighbour-Joining method [37]. The optimal tree with the sum of branch length = 0.04170288 is shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method [38] and are in the units of the number of base substitutions per site. The analysis involved 16 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 876 positions in the final dataset. Evolutionary analyses were conducted in MEGA 6.0 [39]

Full size image
Fig. 2
figure 2

An evolutionary relationship of the taxa of HPV18 isolates based on sequences between the genome positions 7464 and 7839 within the long control region (LCR.). The evolutionary history was inferred using the Neighbor-Joining method [37]. The optimal tree with the sum of branch length = 0.04170288 is shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method [38] and are in the units of the number of base substitutions per site. The analysis involved 16 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 876 positions in the final dataset. Evolutionary analyses were conducted in MEGA 6.0 [39]

Full size image
Fig. 3
figure 3

An evolutionary relationship of HPV45 isolates based on sequences between the genome positions 7074 and 7858 within the long control region (LCR.). The evolutionary history was inferred using the Neighbor-Joining method [37]. The optimal tree with the sum of branch length = 0.04170288 is shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method [38] and are in the units of the number of base substitutions per site. The analysis involved 16 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 876 positions in the final dataset. Evolutionary analyses were conducted in MEGA [39]

Full size image
Table 2 Sequence variations within nucleotide positions 7469 to 7840 (LCR) of HPV16 variants
Full size table

Of the 15 specimen positive for HPV18 isolates, 13 were confirmed as positive by the LCR HPV18 type specific PCR, one of which gave poor sequence data. Nucleotide variations were determined in 35 positions, 22 of which were known variations of the lineages B and C (Africa variants) of HPV 18. The other 13 variations were unique to the isolates of this study (Table 3). The phylogenetic tree (Fig. 2) showed that most of the HPV18 isolate were of the lineages B2 (6 of 12), B3 (2 of 12), B1 (1 of 12) and C (1 of 12) (African variants) and only two isolates (isolates 613 and isolate 858) were of the non-African lineage (lineage A). Of the 18 specimen confirmed as HPV45 genotypes, 13 gave amplification products and 9 has good sequences. (Table 4). Compared to the HPV45 prototype, nucleotide sequence variations were determined in 35 positions, 2 of which were unique to the sequences of 8 isolates of this study. The phylogenetic tree showed that all the HPV45 were of the lineage A1 (Fig. 3).

Table 3 Sequence variations within nucleotide positions 7464 and 7839 (LCR) of HPV18 variants
Full size table
Table 4 Nucleotide sequence variations between positions 7074 and 7858 (LCR) of HPV45 variants
Full size table

Discussion

This long control region of the HPV genome has been shown to be very informative and therefore very useful in the study of HPV variants [7, 12,13,14], as such it was used in this study. A close look at the phylogenetic tree (Fig. 1) suggests that Isolate-937, isolate-866 and isolate-618, may to belong to the HPV16 lineage B (Africa type 1 variants), although it showed that isolates-866, isolate-618 and 618_2 clustered slightly differently. A complete genome sequencing, which was not within the reach of this study would have clarified the lineages of these isolates [5]. Although not clearly understood, repeated PCR and subsequent sequencing of a single specimen produced the isolate-618 and isolate-618_2 with varied branch length and at two positions (Fig. 1). More importantly, the phylogenetic tree showed that the length of sequence obtained in this study were not enough to separate the sub-lineages of the HPV16 lineage B, since the nucleotide positions obtained in this study (Table 2) did not include the position (7438) which distinguishes the sub-lineages B1 and B2 (Af-1a and Af-1b respectively). This was evident by the close clustering of the previously described HPV16 sub-lineage B1 (KF466576) and HPV16 sub-lineage B2 (KF466592) in the phylogenetic tree (Fig. 1). These two variants were distinguished from each other by the substitution A7438C in an earlier study [7]. Therefore, this may imply that the clustering of isolates 866 and 618 may have been different with a longer sequence.

A more detail investigation of the variations in the sequences (Table 2) showed how very similar these isolates were and the bases for the differences seen in the phylogenetic tree. It also showed nucleotide variations unique to the isolates of this study. The occurrence of most of these unique variability specific to the isolates of this study, occurred between the genome positions 7801 and 7828. These may have implication for transcription activities in pathogenic pathways and therefore in oncogenesis, since this range of positions include the binding sites for transcription factors, activator protein 1 (AP1) and Octamer-binding protein 1 [13] which are involved the expression of the genome of the HPV.

The phylogenetic tree presented as Fig. 2 indicates that the sub-lineages B1 and B2 (African sub-lineage 1, old nomenclature) of the HPV18 isolates were the commonest. The sub-lineage B3 and lineage C (African sub-lineage 2, old nomenclature) were also common. Interestingly, only a few were of the A1, A2, and A3 sub-lineages (European sub-lineage (Eur), old nomenclature). A detailed comparison with sequences obtained by an NCBI BLAST and the HPV18 prototype AY262282_A1 [35] showed variations that were unique to the sequence of the isolates of this study, in addition to variations well known for the identified lineages. In all, 6 of the 11 isolates identified in this study had not been described by other studies while the other have been [6, 15]. Most of the nucleotide variations specific to isolates of this study were located within the genome positions 7769 to 7714 and 7795 to 7814; regions where the binding sites for the transcription factors NF1, YY1 and AP1 are located. These changes may have implication for transcription activities in pathogenic pathways [13] and therefore the oncogenic potential of the isolates of this study at the level of sub-lineages. However, it must be noted that these HPVs were isolated from women with normal cytology results. Therefore, a follow-up study may provide more information relating to the potential relationship between these sequence variations and oncogenicity.

All the 9 HPV45 genotypes sequenced were identified without difficulty and were determined to belong to the HPV45 lineage A1 (Table 4). Although the lineages of most of the isolates of this study were identified much clearly, the primary limitation of this study was that in the absence of a complete genome sequence, or the sequence of other genes, such as the E6 and E7 gene, the lineages of some isolate were not clearly identified and the few isolates that did not cluster within known lineage cannot be said to indicate new lineages.

Conclusion

This study provides the first DNA sequence data on the variants/lineages of HPV 16, 18 and 45 isolates from Ghana and has shown that predominately, these genotypes are of the Africa lineages. The observance of specific sequence variation within region associated with transcription binding, indicates a potential for differential oncogenicity of the isolate. The fact that these were obtained from the specimen of women with normal cytology, these data will be useful for comparison with variants of these genotypes from women with the different form of ICC and different grades of cervical lesions to improve our understanding of relationship between lineage and oncogenicity. However, there is need for additional gene or a full genome sequencing to better clarify the classification of a few of the isolates that did not cluster with most of the isolates within a clear lineage.

Abbreviations

ADC:

Adenocarcinoma

ADSC:

Adenosquamous carcinoma

EAS:

European-Asian

HPV:

Human Papillomavirus

ICC:

Invasive cervical cancer

LCR:

Long control region

SCC:

Squamous cell carcinoma

References

  1. GLOBOCAN, IARC. GLOBOCAN. Estimated Cancer incidence, mortality and Prevalence Worldwide in 2012; Fact Sheets. 2012. http://globocan.iarc.fr/Default.aspx. Accessed 14 Jan 2013.

    Google Scholar 

  2. Muñoz N, Bosch FX, de Sanjosé S, Herrero R, Castellsagué X, Shah KV, et al. Epidemiologic Classification of Human Papillomavirus Types Associated with Cervical Cancer. N Engl J Med. 2003;348:518–27.

    Article  PubMed  Google Scholar 

  3. Bosch FX, Qiao Y-L, Castellsagué X. CHAPTER 2 The epidemiology of human papillomavirus infection and its association with cervical cancer. Int J Gynecol Obstet. 2006;94:S8–21.

    Article  Google Scholar 

  4. Bosch FX, Burchell AN, Schiffman M, Giuliano AR, de Sanjose S, Bruni L, Supplement 10, et al. Epidemiology and Natural History of Human Papillomavirus Infections and Type-Specific Implications in Cervical Neoplasia. Vaccine. 2008;26(Supplement 10):K1–16.

    Article  PubMed  Google Scholar 

  5. Burk RD, Harari A, Chen Z. Human papillomavirus genome variants. Virology. 2013;445:232–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Chen Z, Schiffman M, Herrero R, DeSalle R, Anastos K, Segondy M, et al. Evolution and Taxonomic Classification of Alphapapillomavirus 7 Complete Genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70. PLoS ONE. 2013;8:e72565.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Chen AA, Heideman DAM, Boon D, Gheit T, Snijders PJF, Tommasino M, et al. Human Papillomavirus 45 Genetic Variation and Cervical Cancer Risk Worldwide. J Virol. 2014;88:4514–21.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Jelen MM, Chen Z, Kocjan BJ, Hošnjak L, Burt FJ, Chan PKS, et al. Global Genomic Diversity of Human Papillomavirus 11 Based on 433 Isolates and 78 Complete Genome Sequences. J Virol. 2016;90:5503–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Chen AA, Gheit T, Franceschi S, Tommasino M, Clifford GM. Human Papillomavirus 18 Genetic Variation and Cervical Cancer Risk Worldwide. J Virol. 2015;89:10680–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. de Villiers E-M, Fauquet C, Broker TR, Bernard H-U, zur Hausen H. Classification of papillomaviruses. Virology. 2004;324:17–27.

    Article  PubMed  Google Scholar 

  12. Al-Awadhi R, Chehadeh W, Al-Jassar W, Al-Harmi J, Al-Saleh E, Kapila K. Phylogenetic analysis of partial L1 gene of 10 human papillomavirus types isolated most commonly from women with normal and abnormal cervical cytology in Kuwait. Arch Virol. 2013;158:1687–99.

    Article  CAS  PubMed  Google Scholar 

  13. Mendoza L, Picconi MA, Mirazo S, Mongelós P, Giménez G, Basiletti J, et al. Distribution of HPV-16 variants among isolates from Paraguayan women with different grades of cervical lesion. Int J Gynecol Obstet. 2013;122:44–7.

    Article  Google Scholar 

  14. Tornesello ML, Losito S, Benincasa G, Fulciniti F, Botti G, Greggi S, et al. Human papillomavirus (HPV) genotypes and HPV16 variants and risk of adenocarcinoma and squamous cell carcinoma of the cervix. Gynecol Oncol. 2011;121:32–42.

    Article  CAS  PubMed  Google Scholar 

  15. Chen Z, DeSalle R, Schiffman M, Herrero R, Burk RD. Evolutionary Dynamics of Variant Genomes of Human Papillomavirus Types 18, 45, and 97. J Virol. 2009;83:1443–55.

    Article  CAS  PubMed  Google Scholar 

  16. Schmidt M, Kedzia W, Goździcka-Józefiak A. Intratype HPV16 sequence variation within LCR of isolates from asymptomatic carriers and cervical cancers. J Clin Virol. 2001;23:65–77.

    Article  CAS  PubMed  Google Scholar 

  17. Li N, Franceschi S, Howell-Jones R, Snijders PJF, Clifford GM. Human papillomavirus type distribution in 30,848 invasive cervical cancers worldwide: Variation by geographical region, histological type and year of publication. Int J Cancer. 2011;128:927–35.

    Article  CAS  PubMed  Google Scholar 

  18. Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R, et al. Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int J Cancer. 2007;121:621–32.

    Article  CAS  PubMed  Google Scholar 

  19. Altekruse SF, Lacey JV, Brinton LA, Gravitt PE, Silverberg SG, Barnes WA, et al. Comparison of human papillomavirus genotypes, sexual, and reproductive risk factors of cervical adenocarcinoma and squamous cell carcinoma: Northeastern United States. Am J Obstet Gynecol. 2003;188:657–63.

    Article  PubMed  Google Scholar 

  20. Bosch FX, Manos MM, Muñoz N, Sherman M, Jansen AM, Peto J, et al. Prevalence of human papillomavirus in cervical cancer: a worldwide perspective. International biological study on cervical cancer (IBSCC) Study Group. J Natl Cancer Inst. 1995;87:796–802.

    Article  CAS  PubMed  Google Scholar 

  21. Bosch FX, Lorincz A, Muñoz N, Meijer CJLM, Shah KV. The causal relation between human papillomavirus and cervical cancer. J Clin Pathol. 2002;55:244–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Lacey JV, Frisch M, Brinton LA, Abbas FM, Barnes WA, Gravitt PE, et al. Associations between smoking and adenocarcinomas and squamous cell carcinomas of the uterine cervix (United States). Cancer Causes Control. 2001;12:153–61.

    Article  PubMed  Google Scholar 

  23. Awua AK, Sackey ST, Osei YD, Asmah RH, Wiredu EK. Prevalence of human papillomavirus genotypes among women with cervical cancer in Ghana. Infect Agents Cancer. 2016;11:4. doi:10.1186/s13027-016-0050-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Awua AK, Wiredu EK, Osei YD, Sackey ST, Asmah RH, Tettey Y, et al. Oncogenic genotypes of Human Papillomavirus associated with cervical cancer in Ghanaian women. In: BooK of Abstracts. College of Health Sciences, University of Ghana: Accra; 2007.

    Google Scholar 

  25. Mirabello L, Yeager M, Cullen M, Boland JF, Chen Z, Wentzensen N, et al. HPV16 Sublineage Associations With Histology-Specific Cancer Risk Using HPV Whole-Genome Sequences in 3200 Women. J Natl Cancer Inst. 2016;108(9):djw100.

  26. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Quint KD, de Koning MNC, van Doorn L-J, Quint WGV, Pirog EC. HPV genotyping and HPV16 variant analysis in glandular and squamous neoplastic lesions of the uterine cervix. Gynecol Oncol. 2010;117:297–301.

    Article  CAS  PubMed  Google Scholar 

  28. Burk RD, Terai M, Gravitt PE, Brinton LA, Kurman RJ, Barnes WA, et al. Distribution of human papillomavirus types 16 and 18 variants in squamous cell carcinomas and adenocarcinomas of the cervix. Cancer Res. 2003;63:7215–20.

    CAS  PubMed  Google Scholar 

  29. Berumen J, Ordoñez RM, Lazcano E, Salmeron J, Galvan SC, Estrada RA, et al. Asian-American variants of human papillomavirus 16 and risk for cervical cancer: a case–control study. J Natl Cancer Inst. 2001;93:1325–30.

    Article  CAS  PubMed  Google Scholar 

  30. Domfeh A, Wiredu E, Adjei A, Ayeh-Kumi P, Adiku T, Tettey Y, et al. Cervical human papillomavirus infection in Accra, Ghana. Ghana Med J. 2008;42:71–8.

    PubMed  PubMed Central  Google Scholar 

  31. Denny L, Adewole I, Anorlu R, Dreyer G, Moodley M, Smith T, et al. Human papillomavirus prevalence and type distribution in invasive cervical cancer in sub-Saharan Africa: Cervical Cancer in sub-Saharan Africa. Int J Cancer. 2014. doi:10.1002/ijc.28425.

  32. Attoh S, Asmah R, Wiredu EK, Gyasi R, Tettey Y. Human papilloma virus genotypes in Ghanaian women with cervical carcinoma. East Afr Med J. 2010;87:345–9.

    CAS  PubMed  Google Scholar 

  33. Kennedy IM, Haddow JK, Clements JB. A negative regulatory element in the human papillomavirus type 16 genome acts at the level of late mRNA stability. J Virol. 1991;65:2093–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Seedorf K, Krämmer G, Dürst M, Suhai S, Röwekamp WG. Human papillomavirus type 16 DNA sequence. Virology. 1985;145:181–5.

    Article  CAS  PubMed  Google Scholar 

  35. Narechania A, Chen Z, DeSalle R, Burk RD. Phylogenetic incongruence among oncogenic genital alpha human papillomaviruses. J Virol. 2005;79:15503–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Delius H, Hofmann B. Primer-directed sequencing of human papillomavirus types. Curr Top Microbiol Immunol. 1994;186:13–31.

    CAS  PubMed  Google Scholar 

  37. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.

    CAS  PubMed  Google Scholar 

  38. Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci U S A. 2004;101:11030–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors are grateful to the Ghanaian-German Centre for Health Research Ph.D programme and its counterparts (Prof. Fred Binka, Prof. Thomas Junghanss, Prof. Oliver Razum), coordinators (Prof. Richard M K. Adanu, Dr. Micheal Kaeser and Dr. Fenna Veeltmann) and the support staff of this programme. Also, we are grateful to the German Academic Exchange Services (DAAD) for its funding and sponsorship of this programme and to the staff of Department of Epidemiology and Disease Control, School of Public Health, University of Ghana. We sincerely acknowledge the varying contributions of Staff of the Viral Sexually Transmitted Disease Section of the Public Health Agency of Canada and the National Microbiology Laboratory (NML) of Canada for funding and hosting the molecular laboratory analysis of this study.

Funding

German Academic Exchange Services (DAAD) through the Ghanaian-German Centre for Health Research Ph.D programme provided scholarship and support funding for travel, accommodation and living expense. Additionally, the National Microbiology Laboratory of the Public Health Agency of Canada, Winnipeg, hosted, provided technical support and all the inputs for the molecular laboratory analyses on the collected samples.

Availability of data and materials

The DNA sequence data generated and analysed in this study is available in NCBI GenBank with the Accession numbers ranging from KM226842.1 to KM226874.1 (https://www.ncbi.nlm.nih.gov/nuccore/?term=AWUA) and are available from the corresponding author on reasonable request.

Authors’ contributions

Conception of the study and the drafting of proposal: AKA. Review of and contribution to draft proposal and supervising the progress of study: RMKA, EKW, EAA, AS. Conduct of specimen collection: AKA, RMKA, EAA, RHA. Conduct of molecular laboratory analyses: AKA, VAZ, AS. Draft of manuscript: AKA. Review of and contribution to the manuscript: All the authors. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethical approval and consent to participants

The study that collected the specimen from which the HPV isolates were obtained received ethical approval form the Ethics Review Committee of the Ghana Health Service and worked within the guidelines of the Ethical Review Committee and conducted the study in accordance to the ethical standards as declared in Declaration of Helsinki and the Belmont Report. Assurances were given and efforts were made to protect the privacy and confidentiality of the participants enrolled and their collected data. As such personal/demographic data as well as data that will be useful for identifying the participants were not included in this study.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

  1. Department of Epidemiology and Disease Control, School of Public Health, University of Ghana, Accra, Ghana

    Adolf K. Awua & Edwin A. Afari

  2. Cellular and Clinical Research Centre, Radiological and Medical Sciences Research Institute, Ghana Atomic Energy Commission, Accra, Ghana

    Adolf K. Awua

  3. Population, Family and Reproductive Health, School of Public Health, University of Ghana, Accra, Ghana

    Richard M. K. Adanu

  4. Department of Pathology, School of Biomedical and Allied Health Science, University of Ghana, Korle-Bu, Accra, Ghana

    Edwin K. Wiredu

  5. National Microbiology Laboratory, Public Health Agency of Canada, Winnipeg, MB, Canada

    Vanessa A. Zubuch & Alberto Severini

  6. Department of Medical Laboratory Sciences, School of Biomedical and Allied Health Science, University of Ghana, Korle-Bu, Accra, Ghana

    Richard H. Asmah

  7. University of Manitoba, Winnipeg, MB, Canada

    Alberto Severini

Authors
  1. Adolf K. Awua
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Richard M. K. Adanu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Edwin K. Wiredu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Edwin A. Afari
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Vanessa A. Zubuch
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Richard H. Asmah
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Alberto Severini
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Adolf K. Awua.

Additional file

Additional file 1: Table S1.

Isolates’ GenBank accession numbers and lineage classifications. (DOCX 12 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Awua, A.K., Adanu, R.M.K., Wiredu, E.K. et al. Unique LCR variations among lineages of HPV16, 18 and 45 isolates from women with normal cervical cytology in Ghana. Virol J 14, 85 (2017). https://doi.org/10.1186/s12985-017-0755-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-017-0755-z

Keywords

Virology Journal

ISSN: 1743-422X