Skip to main content

Fifteen new nucleotide substitutions in variants of human papillomavirus 18 in Korea

Korean HPV18 variants and clinical manifestation

Abstract

High-risk human papillomavirus (HPV) infection is an essential factor for the development of cervical cancer. HPV18 is the second most common carcinogenic HPV type following HPV16, but the lineages of HPV18 have been less well studied than those of HPV 16. The purpose of this study was to analyze the nucleotide variants in the E6, E7, and L1 genes of HPV18, to assess the prevalence of HPV18 variants in Korea and to explore the relationship between HPV18 genetic variants and the risk for cervical cancer.

A total of 170 DNA samples from HPV18-positive cervical specimens were collected from women admitted to a secondary referral hospital located in Seoul. Among them, the lineages of the 97 samples could be successfully determined by historical nomenclature.

All the studied HPV 18 variants were lineage A. Sublineages A1 and A4 comprised 91.7% (89/97) and 1.0% (1/97), respectively. Sublineages other than A1 or A4 comprised 7.2% (7/97). We identified 15 new nucleotide substitutions among 44 nucleotide substitutions: C158T, T317G, T443G, A560G, A5467G, A5560C, A5678C, A6155G, G6462A, T6650G, G6701A, T6809C, A6823G, T6941C and T6953C. Among them, 6 substitutions at positions 317, 443, 5467, 5560, 6462, and 6823 resulted in amino acid changes (E6: F71L and N113K; L1: H13R, H44P, A345T, and N465S, respectively). The pathologic results were classified as normal in 25.8% (25/97) of the women, atypical squamous cells of undermined significance (ASCUS) in 7.2% (7/97), cervical intraepithelial neoplasia (CIN) 1 in 36.1% (35/97), CIN2/3 in 19.6% (18/97), and carcinoma in 12.4% (12/97). There was no significant association between the HPV18 sublineages and the severity of pathologic lesion or the disease progression.

This study is the first to analyze the distribution of HPV18 variants in Korean and to associate the results with pathologic findings. Although the HPV18 variants had no significant effect on the degree and progression of the disease, the newly discovered nonsynonymous mutation in L1 might serve as a database to determine vaccine efficacy in Korean women.

Introduction

Cervical cancer is the fourth most common cancer among all malignancies in females worldwide and the seventh most common cancer in Korea. According to the World Health Organization’s (WHO’s) GLOBOCAN project in 2018, 569,847 new cases occur and 311,365 people die annually due to cervical cancer worldwide. In Korea, cervical cancer is the seventh most common cancer, with the development of 3348 new cases and 1029 deaths reported annually [1]. Epidemiologic, genetic, immunological and environmental factors are involved in carcinogenesis, and persistent and high-risk human papillomavirus (HPV) infection is an essential factor for the development of cervical cancer. The most deleterious type is HPV16, and the second is HPV18; these two infections are associated with approximately 70% of cervical cancers [2, 3].

HPV is a small double-stranded DNA virus with an 8-kb genome containing early expressed genes (E1, E2, E4, E5, E6, and E7), late genes (L1 and L2) and a long control region (LCR) [4]. The capsid proteins L1 and L2 play critical roles in viral structure formation and the infection process. In particular, purified L1, the major capsid protein, can form empty shells that resemble HPV, which are called virus-like particles (VLPs). These VLPs have hypervariable immunodominant loop structures on the surfaces of the virions that induce humoral immunity without oncogenic activity and are thus extensively used in HPV prophylactic vaccines [5,6,7,8]. E6 and E7 are major oncogenes that are highly expressed in tumors and are related to cellular immortalization, malignant transformation, and carcinogenesis. Based on these roles, proteins E6 and E7 are generally regarded as ideal targets for the development of therapeutic HPV vaccines [6, 9].

Over 200 HPV types have been identified based on L1 sequences. HPV18 variants were originally grouped into European (E), Asian-Amerindian (AA) or African (AFR) lineages according to E6-E7, L1, and/or LCR sequences [10,11,12,13,14]. This classification has been superseded by a whole viral genome sequencing approach that has defined three major lineages (A, B, and C) and additional sublineages (A1 to A5 and B1 to B3) [15] that can be translated from the historical nomenclature (A1 and A2 are AA, A3 to A5 are E and B/C are AFR) [16, 17]. In addition, a recent study published in China proposed new A6 to A8 sublineages and classified them as the E lineage [18].

HPV genetic variants have been reported to differ in their viral assembly, pathogenic potential, and host immune response depending on geographic and ethnic features [2, 17,18,19,20]; however, no data on HPV18 variation among strains from Korean women have been reported thus far. Furthermore, there is some debate regarding the carcinogenic properties of HPV18 lineages [16,17,18,19,20,21,22,23,24].

Methods

From 2010 to 2017, 7992 women admitted to the Seoul National University Boramae Medical Center were tested for cervical HPV genotype more than once. Among them, 3926 (3926/7992 = 49.1%) were positive for HPV, and 170 (170/3926 = 4.3%) were positive for type 18 and negative for other types. HPV detection and typing were performed using a liquid bead microarray, namely, the GeneFinder HPV PCR Kit (Infopia, Seoul, Korea).

Amplification and sequencing of HPV18 E6, E7, and L1 genes were performed using type-specific primers, which are shown in Table 1 [11, 24, 25]. The cycling conditions were as follows: 5 min at 95 °C for initial denaturation; 45 s at 94 °C, 45 s at 55 °C, and 60 s at 72 °C for 35 cycles; and 10 min at 72 °C for final elongation. Amplicons were visualized on 2.0% agarose gels stained with ethidium bromide under UV transillumination. PCR products were automatically sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) and an ABI 3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA) according to the manufacturer’s instructions. E7 sequencing had a high success rate (165/170), so the sequencing was completed with only 1 primer set (=1 trial). However, E6 and L1 sequencing had low success rates, so we attempted 3 trials each (Table 1). All data were confirmed by repeating the PCR amplification and sequence analysis at least twice.

Table 1 Primers used for PCR amplification and sequence analysis

Nucleotide sequences were translated by the translate-tool of ExPASy (http://web.expasy.org/translate/) for the determination of amino acid changes. PSIPRED v.4.0 (http://bioinf.cs.ucl.ac.uk/psipred/) was used for secondary structure prediction, as it provides a simple and accurate secondary structure prediction method.

All sequence data were assembled against the HPV18 reference strain, GenBank sequence NC_001357, using Sequencher version 5.2.3 (Gene Codes Corporation, Ann Arbor, MI, USA). Twenty-eight other sequences were downloaded from GenBank, as in the latest Taizhou area study [18]. The accession numbers used in this analysis were as follows. A1 sublineage: EF202143-EF202145, A2 sublineage: EF202146, A3 sublineage: EF202147-EF202149, A4 sublineage: EF202150-EF202151, A5 sublineage: GQ180787, A6 sublineage: KY457833-KY457836, A7 sublineage: KY45737-KY45840, A8 sublineage: KY457826-KY457827; B1 sublineage: EF202153-EF202155, B2 sublineage: KC470224-KC470225, B3 sublineage: EF202152; and C lineage: KC470229-KC470230. Evolutionary analyses were conducted in MEGA X (http://megasoftware.net) [22, 23].

Based on cytological and histological evaluations of fresh specimens, the cervical lesions were graded according to their severity as follows: normal, atypical squamous cells of undetermined significance (ASCUS), low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), cervical intraepithelial neoplasia grade 1, 2 or 3 (CIN1, 2 or 3) and cervical cancer. The histological diagnosis of each case was reviewed by an experienced pathologist who was unaware of the HPV testing results.

Mann-Whitney, Fisher exact and linear by linear association tests were used for comparisons between AA and E lineages. Variables affecting cervical cancer risk were analyzed by a logistic regression model. All statistical analyses were carried out with SPSS, version 22.0 (IBM, Armonk, NY, USA).

Results

HPV18 variants

Ninety-seven (97/170 = 57.1%) E6, 165 (165/170 = 97.1%) E7, 74 (74/170 = 43.5%) entire L1 and 90 (90/170 = 52.9%) partial L1 sequences from HPV18 isolates were successfully amplified by PCR and sequenced. In total, only 97 samples (E6-E7-L1 sequences) could be classified as lineages and sublineages according to E6-E7 sequence-based historical nomenclature rules of previous studies [9, 10].

The genetic analysis of the regions of the 97 samples is shown in Fig. 1. They were composed of 26 intact and 10 partial E6-E7-L1 sequences (BRM01 ~ 36 variants). In total, we found 44 nucleotide substitutions. Among them, 29 nucleotide substitutions were already reported: T104C, A171G, T173G, C287G, A482C, T485C, C549A, C554T, C751T, C860T, C5478T, A5497G, G5609A, T5619A, C5701G, C5875A, C5920T, A6059G, T6131G, T6146G, A6401G, A6430C, A6441C, C6460G, C6625G, C6842G, A6970G, and G7032A [11,12,13,14, 17, 18, 24, 26,27,28,29], and 15 nucleotide substitutions were newly identified: C158T of BRM05, T317G of BRM07, T443G of BRM08, A560G of BRM10, A5467G of BRM12, A5560C of BRM16, A5678C of BRM17, A6155G of BRM21, G6462A of BRM22, T6650G of BRM23, G6701A of BRM24, T6809C of BRM26, A6823G of BRM27, T6941C of BRM28, and T6953C of BRM36. Among the 15 new nucleotide substitutions, 6 nonsynonymous amino acid substitutions were at positions 317, 443, 5467, 5560, 6462, and 6823 (F71L and N113K in E6; H13R, H44P, A345T, and N465S in L1, respectively). In addition, all the 15 newly identified sequences in this study were submitted to GenBank under the accession numbers MK813921-MK813935 (Fig. 1).

Fig. 1
figure1

Nucleotide sequence variations in the HPV E6-E7-L1 isolates in Korean women. *Previously reported: N, No Grade of cervical lesions: Normal, ASCUS, LSIL/CIN1, HSIL/CIN2 and 3, Cancer L1 secondary structure: H, helix; C, coil; S, strand Abbreviations: AA, Asian-Amerindian; E, European; ASCUS, Atypical squamous cells of undetermined significance; LSIL, Low-grade squamous intraepithelial lesion; HSIL, High-grade squamous intraepithelial lesion; CIN, Cervical intraepithelial neoplasia; na, not assessed; aa, amino acid. E6, E7 and L1 nucleotide positions at which variations were observed are written vertically across the top. Amino acid translations and immunodominant loop structures of L1 are shown at the bottom. The GenBank accession numbers of the newly discovered sequences are on the left. The phylogenetic groupings based on the analysis of E6-E7 are indicated. For each variant sequence, the positions that do not vary relative to the HPV reference are marked with a blank. The regions where sequencing failed are marked with an X. The cervical lesion grade for each variant is shown on the right

The phylogenetic tree analysis was performed using the E6-E7-L1 genes of the 97 Korean HPV18 isolates and 28 already reported variants (Fig. 2). According to the E6-E7 sequence-based previous nomenclature rule [10, 11], all the previous AA lineages (BRM01 ~ 28, total 89 samples) were matched updated A1 sublineage sequences. However, except for one A4 sublineage (BRM32, 1 sample), the other seven variants (BRM29 ~ 31 and 33 ~ 36, total 8 samples) of the previous E lineage did not match the reported sublineages (Figs. 1, 2).

Fig. 2
figure2

Phylogenetic tree of the HPV18 variants by the Maximum Likelihood method. The evolutionary history was inferred using the Maximum Likelihood method with 1000 bootstraps in a Tamura-Nei model. All positions with less than 95% site coverage were eliminated; i.e., fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any position (partial deletion option). Numbers near the line indicate bootstrap values. Evolutionary analyses were conducted in MEGA X

Association between HPV18 variants and cervical lesion severity

All 97 of the women underwent pathological examinations, pap smears (cytological) and/or cervical biopsies (histological) after HPV genotyping within 3 months. The cervical pathologic results were classified as normal in 25 women (25/97 = 25.8%); ASCUS in 7 (7/97 = 7.2%); LSIL/CIN 1 in 35 (35/97 = 36.1%); HSIL/CIN2 and 3 in 18 (18/97 = 18.6%); and carcinoma in 12 (12/97 = 12.4%). There was no statistically significant association between HPV18 lineage and cervical lesion severity (Table 2).

Table 2 Distribution and statistical comparison of HPV18 lineages between the AA and E lineages

Most women underwent HPV genotyping several times; 23 were found to have a different types of HPV infection, aside from type 18, in their lifetime. In particular, 17 women had an infectious history of other high-risk HPVs, including types 16, 31, 33, 35, 39, 45, 51, 52, 58, or 59, as diagnosed according to the International Agency for Research on Cancer (IARC) [30]. There was no association between HPV18 lineages and other types of HPV co-infection (Table 2). Additional infections by other types are not related to the development of cancer (data not shown).

Of the 97 women involved in this study, 54 women underwent an additional follow-up pathological evaluation. Twenty-two women had worsening lesions confirmed by serial pathological tests. There was no association between HPV18 lineages and disease progression (P = .773) (Table 2).

Discussion

From 2010 to 2017, almost half (49.1%) of requested cervical HPV genotyping tests were positive for any type of HPV in a secondary referral hospital in Seoul, Korea. Among the positive results, the frequency of HPV 18 single positive was 4.2%. The frequency of the present study was comparable to other reports in Korea (total HPV prevalence was 16.7% ~ 40.7%, and the HPV 18 prevalence among HPV-positive women was 0.5% ~ 3.6%) [31,32,33].

The distributions of HPV variants differ among geographic origins, evolutionary dynamics, and pathogenicity. In our population, 6 substitutions, namely, C287G in E6 and G5503A, C5701G, C6460G, C6625G, and C6842G in L1, were found in all HPV18 variants. These 6 substitutions were also found in all HPV18 isolates in southeastern and northeastern China [11, 18], which Korea is located next to, but found in 40% of HPV 18 in southwest and central China [24]. These findings support the geographical distribution of HPV lineages.

We previously discovered that the HPV 16 variants were composed of A1 ~ 3, A4, C, and D sublineages (54.1, 37.8, 0.7, and 7.4% of samples, respectively) in Korea [21]. Upon comparing HPV 16, little diversity was evident, indicating that previous lineages AA and E comprised 91.7 and 8.3%, respectively, in Korean HPV 18.

Previous studies reported that the risk of developing high-grade CIN is significantly increased with the non-European variants [13, 34]. One study [35] reported that the AA and European variants had significantly higher associations with pre-invasive lesions than the African variants. In contrast, other studies showed that no significant difference in pre-invasive lesion risk was observed between the variant lineages (A, B, and C) [14, 17, 26]. Our results are in line with the latter conclusion; there was no statistically significant association between HPV18 lineages and cervical pathologic lesions in Korea.

In the aspect of public health, many countries have implemented national policies of HPV vaccination. In Korea, HPV vaccination has been free for 12-year-old girls since 2016, and the government intends to expand the vaccination targets. Three prophylactic vaccines have received licensure from Korea Food and Drug Administration: the AS04-adjuvanted bivalent (HPV16/18) vaccine (Cervarix®, GlaxoSmithKline, Belgium), which was licensed in 2008; the aluminum hydroxyphosphate sulfate (AAHS) adjuvant quadrivalent (HPV6/11/16/18) vaccine (Gardasil®, Merck, US), which was licensed in 2007; and the AAHS adjuvant 9-valent (HPV6/11/16/18/31/33/45/52/58) vaccine (Gardasil®9, Merck, US), which was licensed in 2016. These prophylactic HPV vaccines are composed of L1 proteins of multiple HPV combinations [5, 7]. The loop structures of the HPV L1 major capsid protein contribute to the epitopes of vaccine-induced cross-neutralizing antibodies. Therefore, amino acid changes in the L1 loop region could be a critical issue for vaccination development. Our data on the genetic diversity of the HPV18 variants in Korea show two nonsynonymous substitutions in the loop structures, L64M within the BC loop (of BRM34) and T149N within the DE loop (of BRM20, BRM31 ~ 36). It may be helpful to design second-generation prophylactic HPV vaccines and implement feasible nationwide vaccination programs.

With the development of next-generation sequencing (NGS), the whole-genome sequencing (WGS) of 8 kb HPV became easier, making the analysis of lineages and single-nucleotide polymorphisms (SNPs) relatively faster and more accurate. However, WGS analysis pipelines in microbiological fields have not yet been established systematically, and it is difficult to analyze multiple samples with limited resources. In addition, although WGS identifies more variants and contributes to the construction of more accurate phylogenetic trees than partial sequencing, the data composed of E6-E7-L1 sequences over the past 20 years does not have significantly reduced reliability compared to WGS. Thus far, many studies have been selectively conducted on the oncogenic proteins E6 and E7, and the major capsid protein L1 plays an important role in the prophylactic vaccine, as our study suggests.

Conclusions

In summary, this research identified 15 new nucleotide substitutions and is the first to analyze the distribution of HPV18 variants in Korean women and to match the results to pathologic results. Data regarding the geographic/ethnic HPV18 genetic diversity distributions may be helpful for designing diagnostic probes, correlating epidemiologic cervical cancer risks, and analyzing the efficacy of HPV vaccines for targeted populations in Korea.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

AA:

Asian-Amerindian

AAHS:

Aluminum hydroxyphosphate sulfate

AFR:

African

ASCUS:

Atypical squamous cells of undermined significance

CIN:

Cervical intraepithelial neoplasia

E:

European

E :

Early gene

HPV:

Human papillomavirus

HSIL:

High-grade squamous intraepithelial lesion

IARC:

International Agency for Research on Cancer

LSIL:

Low-grade squamous intraepithelial lesion

L :

Late gene

LCR :

Long control region

NGS:

Next-generation sequencing

SNP:

Single-nucleotide polymorphism

VLP:

Virus-like protein

WGS:

Whole-genome sequencing

WHO:

World Health Organization

References

  1. 1.

    WHO Global Cancer Observatory [Available from: https://gco.iarc.fr/].

  2. 2.

    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(3):621–32.

    CAS  Article  Google Scholar 

  3. 3.

    Munoz N, Bosch FX, Castellsague X, Diaz M, de Sanjose S, Hammouda D, et al. Against which human papillomavirus types shall we vaccinate and screen? The international perspective. Int J Cancer. 2004;111(2):278–85.

    CAS  Article  Google Scholar 

  4. 4.

    Bernard HU. Taxonomy and phylogeny of papillomaviruses: an overview and recent developments. Infect Genet Evol. 2013;18:357–61.

    Article  Google Scholar 

  5. 5.

    Kirnbauer R, Booy F, Cheng N, Lowy DR, Schiller JT. Papillomavirus L1 major capsid protein self-assembles into virus-like particles that are highly immunogenic. Proc Natl Acad Sci. 1992;89(24):12180–4.

    CAS  Article  Google Scholar 

  6. 6.

    Doorbar J, Quint W, Banks L, Bravo IG, Stoler M, Broker TR, et al. The biology and life-cycle of human papillomaviruses. Vaccine. 2012;30(Suppl 5):F55–70.

    CAS  Article  Google Scholar 

  7. 7.

    Bissett SL, Godi A, Beddows S. The DE and FG loops of the HPV major capsid protein contribute to the epitopes of vaccine-induced cross-neutralising antibodies. Sci Rep. 2016;6:39730.

    CAS  Article  Google Scholar 

  8. 8.

    Shen-Gunther J, Cai H, Zhang H, Wang Y. Abundance of HPV L1 Intra-Genotype Variants With Capsid Epitopic Modifications Found Within Low- and High-Grade Pap Smears With Potential Implications for Vaccinology. Front Genet. 2019;10:489.

    CAS  Article  Google Scholar 

  9. 9.

    Chabeda A, Yanez RJR, Lamprecht R, Meyers AE, Rybicki EP, Hitzeroth II. Therapeutic vaccines for high-risk HPV-associated diseases. Papillomavirus Res(Amsterdam, Netherlands). 2018;5:46–58.

    Article  Google Scholar 

  10. 10.

    Ong CK, Chan SY, Campo MS, Fujinaga K, Mavromara-Nazos P, Labropoulou V, et al. Evolution of human papillomavirus type 18: an ancient phylogenetic root in Africa and intratype diversity reflect coevolution with human ethnic groups. J Virol. 1993;67(11):6424–31.

    CAS  Article  Google Scholar 

  11. 11.

    Sun Z, Liu J, Wang G, Zhou W, Liu C, Ruan Q. Variant lineages of human papillomavirus type 18 in Northeast China populations characterized by sequence analysis of E6, E7, and L1 regions. Int J Gynecol Cancer. 2012;22(6):930–6.

    Article  Google Scholar 

  12. 12.

    Schlecht NF, Burk RD, Palefsky JM, Minkoff H, Xue X, Massad LS, et al. Variants of human papillomaviruses 16 and 18 and their natural history in human immunodeficiency virus-positive women. J Gen Virol. 2005;86(10):2709–20.

    CAS  Article  Google Scholar 

  13. 13.

    Arroyo SL, Basaras M, Arrese E, Hernáez S, Andía D, Esteban V, et al. Human papillomavirus (HPV) genotype 18 variants in patients with clinical manifestations of HPV related infections in Bilbao, Spain. Virology J. 2012;9:258.

    Article  Google Scholar 

  14. 14.

    Chen AA, Gheit T, Franceschi S, Tommasino M, Clifford GM. Human papillomavirus 18 genetic variation and cervical cancer risk worldwide. J Virol. 2015;89(20):10680–7.

    CAS  Article  Google Scholar 

  15. 15.

    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(8):e72565.

    CAS  Article  Google Scholar 

  16. 16.

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

    CAS  Article  Google Scholar 

  17. 17.

    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(3):1443–55.

    CAS  Article  Google Scholar 

  18. 18.

    Xu HH, Zheng LZ, Lin AF, Dong SS, Chai ZY, Yan WH. Human papillomavirus (HPV) 18 genetic variants and cervical cancer risk in Taizhou area. China Gene. 2018;647:192–7.

    CAS  Article  Google Scholar 

  19. 19.

    Sait K, Turki R, Abuzenadah AM, Jiffiri OH, Bohmaidah A, Sohrab SS. Genetic diversity and phylogenetic analysis of HPV 16 & 18 variants isolated from cervical specimens of women in Saudi Arabia. Saudi J Biol Sci. 2019;26(2):317–24.

    CAS  Article  Google Scholar 

  20. 20.

    Frati E, Bianchi S, Amendola A, Colzani D, Petrelli F, Zehender G, et al. Genetic characterization of variants of HPV-16, HPV-18 and HPV-52 circulating in Italy among general and high-risk populations. Mol Med Rep. 2019;21:894.

    PubMed  Google Scholar 

  21. 21.

    Park JS, Shin S, Kim EC, Kim JE, Kim YB, Oh S, et al. Association of human papillomavirus type 16 and its genetic variants with cervical lesion in Korea. APMIS. 2016;124(11):950–7.

    CAS  Article  Google Scholar 

  22. 22.

    Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30(5):1229–35.

    CAS  Article  Google Scholar 

  23. 23.

    Kumar S, Stecher G, Li M, Knyaz C, Tamura K, MEGA X. Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.

    CAS  Article  Google Scholar 

  24. 24.

    Shen M, Ding X, Li T, Chen G, Zhou X. Sequence variation analysis of HPV-18 isolates in southwest China. PloS one. 2013;8(2):e56614.

    CAS  Article  Google Scholar 

  25. 25.

    Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3--new capabilities and interfaces. Nucleic acids research. 2012;40(15):e115-e. [Available from: http://bioinfo.ut.ee/primer3-0.4.0/].

    CAS  Article  Google Scholar 

  26. 26.

    Arias-Pulido H, Peyton CL, Torrez-Martinez N, Anderson DN, Wheeler CM. Human papillomavirus type 18 variant lineages in United States populations characterized by sequence analysis of LCR-E6, E2, and L1 regions. Virology. 2005;338(1):22–34.

    CAS  Article  Google Scholar 

  27. 27.

    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(8):1763–8.

    CAS  Article  Google Scholar 

  28. 28.

    Liu Y, Pan Y, Gao W, Ke Y, Lu Z. Whole-genome analysis of human papillomavirus types 16, 18, and 58 isolated from cervical Precancer and Cancer samples in Chinese women. Sci Rep. 2017;7(1):263.

    Article  Google Scholar 

  29. 29.

    van der Weele P, Meijer CJLM, King AJ. High whole-genome sequence diversity of human papillomavirus type 18 isolates. Viruses. 2018;10(2):68.

    Article  Google Scholar 

  30. 30.

    IARC 2019 [Available from: https://monographs.iarc.fr/agents-classified-by-the-iarc/].

  31. 31.

    So KA, Hong JH, Lee JK. Human papillomavirus prevalence and type distribution among 968 women in South Korea. J Cancer Prev. 2016;21(2):104–9.

    Article  Google Scholar 

  32. 32.

    Lee EH, Um TH, Chi H-S, Hong Y-J, Cha YJ. Prevalence and distribution of human papillomavirus infection in Korean women as determined by restriction fragment mass polymorphism assay. J Korean Med Sci. 2012;27(9):1091–7.

    Article  Google Scholar 

  33. 33.

    Nah EH, Cho S, Kim S, Cho HI. Human papillomavirus genotype distribution among 18,815 women in 13 Korean cities and relationship with cervical cytology findings. Ann Laboratory Med. 2017;37(5):426–33.

    Article  Google Scholar 

  34. 34.

    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 General Virol. 2000;81(Pt 12):2959–68.

    CAS  Article  Google Scholar 

  35. 35.

    Xi LF, Koutsky LA, Hildesheim A, Galloway DA, Wheeler CM, Winer RL, et al. Risk for high-grade cervical intraepithelial neoplasia associated with variants of human papillomavirus types 16 and 18. Cancer Epidemiol Biomarkers Prev. 2007;16(1):4–10.

    CAS  Article  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

Not applicable.

Author information

Affiliations

Authors

Contributions

NK and SS: Study concept and design; NK and JEK: Performed the experiments; NK, JSP, and JEK: Acquisition of data; NK and SS: Analysis and interpretation of data; NK and SS: Drafting of the manuscript; JHP, HP, EYR, and JHY: Critical revision of the manuscript for important intellectual content; and SS: Study supervision. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Sue Shin.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of Seoul National University Boramae hospital (30–2018-40).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kim, N., Park, J.S., Kim, J.E. et al. Fifteen new nucleotide substitutions in variants of human papillomavirus 18 in Korea. Virol J 17, 70 (2020). https://doi.org/10.1186/s12985-020-01337-7

Download citation

Keywords

  • Human papillomavirus (HPV) 18
  • E6, E7 and L1 genes
  • Variants
  • Lineage
  • Cervical cancer