Skip to main content

Human papillomavirus genotyping by Linear Array and Next-Generation Sequencing in cervical samples from Western Mexico

Abstract

Background

The Linear Array® (LA) genotyping test is one of the most used methodologies for Human papillomavirus (HPV) genotyping, in that it is able to detect 37 HPV genotypes and co-infections in the same sample. However, the assay is limited to a restricted number of HPV, and sequence variations in the detection region of the HPV probes could give false negatives results. Recently, 454 Next-Generation sequencing (NGS) technology has been efficiently used also for HPV genotyping; this methodology is based on massive sequencing of HPV fragments and is expected to be highly specific and sensitive. In this work, we studied HPV prevalence in cervixes of women in Western Mexico by LA and confirmed the genotypes found by NGS.

Methods

Two hundred thirty three cervical samples from women Without cervical lesions (WCL, n = 48), with Cervical intraepithelial neoplasia grade 1 (CIN I, n = 98), or with Cervical cancer (CC, n = 87) were recruited, DNA was extracted, and HPV positivity was determined by PCR amplification using PGMY09/11 primers. All HPV- positive samples were genotyped individually by LA. Additionally, pools of amplicons from the PGMY-PCR products were sequenced using 454 NGS technology. Results obtained by NGS were compared with those of LA for each group of samples.

Results

We identified 35 HPV genotypes, among which 30 were identified by both technologies; in addition, the HPV genotypes 32, 44, 74, 102 and 114 were detected by NGS. These latter genotypes, to our knowledge, have not been previously reported in Mexican population. Furthermore, we found that LA did not detect, in some diagnosis groups, certain HPV genotypes included in the test, such as 6, 11, 16, 26, 35, 51, 58, 68, 73, and 89, which indicates possible variations at the species level.

Conclusions

There are HPV genotypes in Mexican population that cannot be detected by LA, which is, at present, the most complete commercial genotyping test. More studies are necessary to determine the impact of HPV-44, 74, 102 and 114 on the risk of developing CC. A greater number of samples must be analyzed by NGS for the most accurate determination of Mexican HPV variants.

Background

Cervical cancer (CC) is the most frequent female cancer worldwide, occupying fourth place among cancers in women. In Mexico, the incidence of this pathology is high (13,960 new cases each year), reaching second place for cancers in women and causing 4769 deaths each year according to Globocan 2012 [1]. Infection with HPV is the principal etiological factor involved in the development of CC [2]. To date, >202 HPV genotypes have been characterized [3] and >40 have been reported to infect the female genital tract [4, 5]. Globally, the HPV genotypes mainly associated with CC include HPV-16, 18, and 45 [6], but these data vary depending on the geographic area. For example, while HPV-35 and 45 are more prevalent in Africa, the genotypes 52 and 58 are more frequently observed in Asia, whereas 33 and 31 are more prevalent in Europe and the USA [79]. Regarding Cervical intraepithelial neoplasia grade 1 (CIN I), HPV-42 and 84 were reported as more prevalent in Sweden [10], while in Chile, HPV-61 and 89 [11] are most prevalent, and in Asia, HPV-52 [12] and in the USA (in rural American Indian women), HPV-61 [13]. In Mexico, a recent study of our research group, using genotyping by Linear Array, show that the HPV genotypes most frequently found in cervical cancers were 16 (62.8 %), 18 (11.6 %), 45 (8.3 %), 52/58 (6.6 %), and 39 (5.8 %). Notably, a high percentage (58–64 %) of coinfections were found in CIN I—CIN III, suggesting that coinfections could contribute to regression or progression of the cervical lesions [14].

All of these data suggest the need to develop detection and prevention strategies consistent with the characteristics of the geographic population.

There are few commercial methods available for HPV genotyping [15], the majority of these are based on PCR and reverse hybridization. However, this approach entertains the drawback of only recognizing HPV genotypes with sufficient dissimilarity for discrimination by probes used in hybridization, allowing the detection of a restricted number of HPV genotypes and not discriminating among nucleotide variants. Next-Generation Sequencing (NGS), or Massive parallel sequencing (MPS), has emerged as a powerful tool for genotyping, possessing the advantage of analyzing many samples at the same time. Recently, this methodology has been used for the genotyping of HPV [1618]. In Mexico, there are few reports of HPV genotyping by means of LA and none, to our knowledge, employing NGS. In this work, we identified, by LA and NGS, the HPV genotypes infecting the cervixes of woman from Western Mexico.

Results

HPV genotyping by Linear Array test

From the 48 WCL samples, 12 (25 %) were positive for HPV. Considering only HPV-positive samples, 58.3 % were reported as SI and 41.7 % as MI (Fig. 1). In the CIN I group, 31 samples were found to be HPV-negative (31.6 %) and 67 were HPV-positive (68.4 %). Among the HPV-positive samples, 37 samples were SI (55.2 %) and 30 were MI (44.8 %). In the CC sample group, 100 % (n = 87) of the samples were infected with HPV, 67.8 % with SI and 32.2 % with MI.

Fig. 1
figure 1

Schematic representation of samples genotyped by Linear Array (LA). This figure depicts the number of patients analyzed in each diagnosis group (n), the percentage of Human papillomavirus (HPV) positivity, the percentage of HPV present as Single infection (SI) or Multiple infection (MI), and the top five HPV genotypes present in each group. WCL = Without cervical lesion, CIN I = Cervical intraepithelial neoplasia grade I, CC = Cervical cancer. (−) HPV negative with LA, (+) HPV positive with LA. HPV genotypes carcinogenic to humans are marked in red (International Agency for Research on Cancer [IARC] classification)

Calculating the frequency of HPV genotypes detected by LA in the different groups analyzed (Table 1), we observed that the most frequent HPV genotypes found in WCL samples were 16, 51, and 84 (each found in the 25 % of HPV-infected women), followed by 56, 59, and 66 (16.7 % each), and last, HPV genotypes 6, 11, 31, 39, 52, 68, 81, and 89 (8.3 % frequency). Analysis in the CIN I group showed that HPV-16 was the most frequent genotype found in 20 of the 67 HPV-infected patients (29.9 %), followed by genotypes 84 (17.9 %), 89 (13.4 %), 53 (10.4 %), 31 (9 %), and HPVs-39, 56, 58, 61, and 62 (7.5 %). HPV-66 was present in 6 % of the women of this group, 16 additional HPV genotypes were found at frequencies less than 4.5 %, as detailed in Table 1. On the other hand, the CC group exhibited that HPV-16 was also the most frequent genotype found (in 42 of 87 HPV-infected patients, 48.3 %), followed by 18 (18.4 %), 45 (10.3 %), 39 (9.2 %), 31, 52, and 71 (with 6.9 % frequency each), 68 (5.7 %), and 19 additional HPV genotypes with frequencies ranging between 4.6 and 1.1 %. It is noteworthy that HPV-16 was the most frequent genotype found in all analyzed groups. Interestingly, HPV-52, 56, and 68 were found only in patients with MI; and HPV-26, 69, 70, and 72, were detected only in the CC group.

Table 1 Frequencies of each Human papillomavirus (HPV) genotype identified by Linear Array (LA) according to diagnosis group

HPV genotyping by Next-Generation Sequencing

To corroborate the genotypes detected by LA and to determine additional HPV not included in this test, NGS was performed. All samples positive for HPV by LA, as well as negative for LA but positive with PGMY09/11 primers, were included. Samples were pooled according to diagnosis and SI or MI infection. Samples negative by LA but positive with PGMY09/11 primers were designated as WCL (−) or CIN (−). As depicted in Fig. 2a, after comparison of the sequences obtained with the human L1 papillomavirus database downloaded from PaVE, 75.4 % (58,481) of the reads were reported as unmapped or not recognized as belonging to any reported HPV, 7.2 % (5549 reads) were tooshort (<100 pb), and finally, 17.4 % (13,482 reads) were found to be mapped with HPV sequences. Median length of HPV reads was 450 nucleotides (nt), ranging between 404 and 468, while the modal length of the reads was 460 nt.

Fig. 2
figure 2

Next-Generation Sequencing (NGS) data. a Percentage of total reads obtained in the entire sequencing, grouped as unmapped, Human papillomavirus (HPV)-mapped, or TooShort (eliminated). b Diagram of total reads obtained in each pooled study group (gray), including HPV mapped reads (blue) obtained in each group utilizing the HPV reference sequence from the Papillomavirus Episteme (PaVE) database. c Total reads mapped with each HPV genotype independent of the diagnosis group

Comparing the number of reads obtained by group, it was observed that the greatest number of reads was obtained in the WCL (−) group, followed by the CIN I (−), CIN I SI, CC SI, CC MI, WCL (+), and CIN I MI groups (Fig. 2b). However when we analyzed the proportion of reads belonging to HPV by group, we observed that the CC MI group produced the greatest number of reads, followed by the CC SI, CIN SI, and CIN MI groups. Groups CIN (−), WCL (−), and WCL (+) had low numbers of reads belonging to HPV in relation to the numbers of reads obtained in these groups. Of the total reads that mapped to HPV, 228 sequences belonged to the WCL (−) group, 121 to WCL (+), 118 to CIN I (−), 2685 to CIN I SI, 403 to CIN I MI, 4562 to CC SI, and 5365 to the CC MI group (Fig. 2b). Analyzing the total reads obtained for each HPV genotype, the majority of reads were from HPV-16 (4359 reads), followed by HPV-45 (3188 reads), HPV-68 (1244 reads), HPV-18 (1040 reads), and HPV-31 (617 reads) (Fig. 2c).

As depicted in Table 2, reads for each HPV were also determine according to diagnosis group. In the WCL (−) group (samples WCL negative for HPV by LA), we observed the presence of HPV-26, 32, 44, and 51; of these HPV, it is noteworthy that HPV-26 and 51 are included in the LA test. In the WCL (+) group, nearly all genotypes reported by LA were found, but additionally, HPV-52, 81, and 89 were also found (despite that all of these genotypes are included in the LA test). On the other hand, the presence of HPV-6, 11, 31, 39, 56, and 68, which were determined by LA in this group, was not confirmed by NGS. Regarding the CIN I (−) group, we were able to find, by means of NGS, 85 reads mapping to HPV-16 and 33 reads that mapped to HPV-35, despite that the test included both HPV genotypes. With respect to the CIN I SI group, HPV-44 and 74 were additionally found, genotypes that are not included in LA. Similarly, in the CIN I MI group, additional HPV genotypes that are not included in LA were detected by NGS, such as 44, 102, and 114. In terms of the CC SI group, HPV-11, 58, 68, 73, 81, 89, and 114 were identified by NGS, but not by LA. Among these, the LA test possesses probes for nearly all HPV, but not for HPV-102 nor for HPV-114. Finally, in the CC MI group, HPV-58, 74, 102, and 114 were additionally detected (Table 2). For comparison, all results obtained by LA in each sample analyzed are contained in Additional file 1: Table S1.

Table 2 Number of reads obtained for each Human papillomavirus (HPV) genotype according to study group

For a better overview of all NGS results, a summary of HPV genotypes found according to diagnosis group are depicted in Table 3, in which the HPV genotypes found by LA and NGS are schematized in black; HPV genotypes, in blue, found only by NGS despite their being included in the LA test; and in pink, genotypes found only by NGS that are not included in LA. Those HPV genotypes detected by LA and not by NGS were also included (asterisk). Interestingly, by using NGS, the following genotypes not included in the LA test were identified: HPV-32, 44, 74, 102, and 114. If we analyze the presence of each HPV by diagnosis group, we may observe that HPV-16 was detected in six of the seven analyzed groups, HPV-59 was detected in five groups, HPV-31 and 66 in four groups, and the remainder of HPV genotypes was found in three or two groups. Genotypes identified in only one group comprised HPV-11 in CC SI, 32 in WCL (−), 56 in CIN I MI, 67 in CIN I SI, 69 in CC MI, and 70 in CC MI, as can be observed in Table 3.

Table 3 Summary of all Human papillomavirus (HPV) genotypes determined by Next-Generation Sequencing (NGS) in each study group

Finally, consensus sequences (contigs) of each HPV genotype found exclusively among the Mexican samples analyzed were aligned and a phylogenetic tree was constructed (Fig. 3). As expected, HPV contigs were grouped according to the present classification. Exception included HPV-54, found more closely related with HPV-6, 11, 44, and 74. To supplement these data, next to each genotype is displayed the diagnosis group in which the sequences were found (WCL, CIN I, or CC). To contrast these results, we noted that HPV-11, 18, 68, 69, and 70 were only found in CC, whereas HPV-32 and 51 were only found in the WCL group. On the other hand, HPV-31, 33, 35, 39, 45, 54, 58, 62, 73, 74, 102, and 114 were present in CIN I and CC, whereas HPV-44 and 84 were found in the WCL and CIN I groups.

Fig. 3
figure 3

Phylogenetic tree built with the total Human papillomavirus (HPV) consensus sequences identified in Mexican cervical samples analyzed. Analysis was performed by MEGA6 software displayed as Topology. The diagnosis group in which each HPV genotype was detected is also depicted. Blue circle: WCL (Without cervical lesion), Orange circle: CIN I (Cervical intraepithelial neoplasia grade I), and Red circle: CC (Cervical cancer)

Discussion

Cervical cancer remains a public health problem in Mexico, despite the prevention and vaccination campaigns implemented by the government. Vaccination is expected to decrease the incidence of this malignancy in the near future, and this probably would impact differentially, depending on the geographic region, because it has been shown that HPV genotype frequency varies among the different countries [1922]. Thus, there are local efforts to develop population-specific directed HPV vaccines; however, to develop a vaccine specific for a geographic population, it is essential to first know the HPV genotypes circulating in the population. In this study, we detected the HPV genotypes that are infecting women from Western Mexico, not only by LA, which is the commercial test with the major number of genotypes included, but also by NGS, to identify the widest possible genotypes. Results obtained by LA were very similar to those previously reported and discussed by our IMSS Research Network on HPV [14]; however, the use of 454 NGS allowed us to identify previously not reported HPV genotypes in Mexican population. NGS is a technology that has been efficiently used for HPV genotyping; this methodology is based on Massive Sequencing (MS) of HPV fragments from a great number of genotypes, even when these are in low in copy number [16, 18].

Using 454 NGS, we could detect genotypes not included in LA: 32, 44, 74, 102, and 114 (1, 79, 2, 2, and 30 reads, respectively). We found HPV-32 only in the WCL group; this genotype has been associated with oral focal epithelial hyperplasia [23, 24] and oral papillomatous lesions from patients with Human immunodeficiency virus (HIV) infection [25]. In another study, in a search for uncommon and rare HPV genotypes in cervical carcinomas, HPV-32 was not found [26]. On the other hand, characterization of the properties of HPV-32 E6 and E7 oncoproteins have provided evidences of the benign nature of this genotype [27].

With respect to HPV-44, we detected this genotype in WCL or CIN I samples, but not in CC. This HPV was isolated for the first time from a vulvar condyloma, and was found in five cases of 439 normal squamous epithelium samples and in three cases of 112 CIN I samples, but was not found in any of the 56 CC analyzed [28]. However, this genotype has been found in patients with invasive cervical cancer in Africa with a frequency of 1.8 % of co-infection with other HPV genotypes [29]. Additionally, it has been detected in patients with HIV with a frequency of 12.5 % in South Africa, Ghana, and Nigeria [29], and has been also reported in anogenital carcinoma [30]. In other studies, the prevalence of HPV-44 has been reported as ranging from 0.4–5 % [3134].

HPV-74 has been found in low-grade vaginal intraepithelial neoplasia in immunosuppressed women and has been related phylogenetically with HPV-6, 11, 44, and 55 [35, 36]. We found this genotype only in some patients of the CIN I SI and CC MI groups. The number of reads obtained for HPV-74 was very low, suggesting a low prevalence. In this regard, low prevalence of HPV-74 has been detected by LiPA assay in genital warts in a German population [34]; additionally, HPV-74 has also been reported as present in male populations from Brazil, Mexico, and the USA (by direct sequencing of products from nested PCR of the L1 region) [37].

HPV-102 has been also found in this study only by NGS; we detected it in CIN I and CC, both in MI (Table 3). Despite that we found only one read in each group, we consider this genotype as present, because the sequences obtained are sufficient (337 and 411 bp) and identity was 99 and 100 %, respectively. Additionally, phylogenetic tree analysis places this genotype in alpha species 3 (Fig. 3), as also reported [36]. It is important to mention that we used all HPV reference sequences from the L1 alpha species for mapping. It has been reported that “simultaneous mapping” significantly decreases the read numbers obtained [38], but we consider that it is better because more solid data are obtained. In agreement with our findings, this genotype has been detected in cervicovaginal cells from Hispanic females [39].

Other additional genotype detected in our study, only by NGS, was HPV-114; this is a recently characterized papillomavirus isolated from a low-grade cervical lesion, and it has been found in low-grade samples with 1.7 % prevalence [40]. We found this genotype in CIN I MI and CC, both SI and MI; its presence in CC suggest a possible oncogenic risk that needs to be evaluated individually in more samples. HPV-114 has also been detected by NGS in a Swedish cervical sample [17], and has been found in one case of Vulvar intraepithelial neoplasia (VIN) [41].

To complement and summarize our results, it was in our interest to performed a phylogenetic tree exclusively with sequences obtained by NGS from Mexican samples (Fig. 3); reconstruction shows a very similar pattern according to the classification of de Villiers [36] despite using approximately 450 bp of the L1 region amplified by the PGMY09/11 set of primers. Bernard et al. in 1994 demonstrated that a short 291-bp fragment of the L1 gene is sufficient to generate phylogenetic comparisons [42].

It is accepted that in order to recognize a new Papillomavirus (PV) type, its complete genome must be cloned and its L1 gene sequence must be different in 10 % from the closest type known; differences between 2 and 10 % represent a subtype, and differences <2 %, a variant [43]. By analysis with the GS Reference Mapper, a variety of putative high-confidence nucleotide differences were determined, principally in HPV genotypes 31, 39, 45, 62, 66, 68, and 69. These nucleotide differences were found in two or more diagnosis groups (taking into account differences detected in at least 50 reads). Therefore, we continue to not discard the existence of HPV variants present in Mexico.

Regarding the classification of carcinogenicity for humans used by the International Agency for Research on Cancer (IARC) [44], there are some HPV genotypes that are classified as “probably” or “possibly” carcinogenic or not classifiable. Among those, we found HPVs-11, 68, 70, and 69 only in CC samples. In contrast, HPVs-32, 44, 51, and 84 were observed in only the WCL and/or in the CIN I group. With regard to HPV genotypes not previously reported in Mexican population, HPV-74, 102, and 114 were detected in the CIN I and CC groups.

One limitation of this study was that groups of samples (according to the diagnosis), but not individually samples, were sequenced. Additionally, the HPV-genotypes detected were limited to the amplicons obtained by using the PGMY09/11 set of primers.

We continue to recruit and analyze additional cervical samples in order to corroborate and report these observations. Genotyping of HPV by NGS will permit the detection of more possible HPV variants, whose sequences could increase or diminish HPV recognition by the immune system or modify the oncogenic potential of HPV proteins [45].

We continue to need more studies with methodologies such as NGS to determine the HPV genotypes and variants circulating in Mexican population and at the different stages during CC progression.

Conclusions

To our knowledge, this is the first time that data has been reported on HPV genotypes circulating in Mexico by using Massive parallel sequencing (MPS). Our data provide helpful information because we were able to find not previously detected HPV genotypes in Mexican population. The methodological strategy employing the PGMY09/11 primers appears to be a very good system.

Methods

Sample collection

Approximately 250 cervical samples were obtained from the Centro Médico Nacional de Occidente (CMNO, IMSS) located in Guadalajara, Jalisco, Western Mexico. The women whose samples were taken attended a medical examination due to a cervical ailment. Diagnosis was first performed only by colposcopy observation (visual inspection using acetic acid and Lugol’s iodine solution); when cervical lesions/cancer were (was) observed, a biopsy was taken for diagnosis confirmation by histopathology. From the 250 samples obtained, 17 samples were excluded (those samples without sufficient quantity and quality of DNA), therefore, 233 samples were included in this study: 48 WCL; 98 CIN I, and 87 CC. Ages of the participating women ranged from 25–65 years. Samples were obtained by inserting a cytobrush into the endocervical canal, rotating it for 3–5 full turns, and placing it into ThinPrep PreservCyt® transport medium solution (Cat. no. 70097–002; Hologic, Inc., Marlborough, MA, USA).

HPV screening and Linear Array® HPV genotyping test

Total DNA was extracted and purified from each cervical sample using the AmpliLute Liquid Media Extraction Amplicor Kit (Cat. no. 03750540 190; Roche Molecular Systems, Inc., Branchburg, NJ, USA.). A first screening for HPV positivity was performed by conventional single-round PCR utilizing PGMY09/11 primers, which amplified a 450-bp fragment of L1 [46]. HPV-positive samples were genotyped utilizing the Linear Array® HPV Genotyping test (Cat. no. 04391853 190; Roche Molecular Diagnostics, Pleasanton, CA, USA), which is capable of detecting simultaneously up to 37 HPV genotypes in a sample [47]. This assay includes, as control, β-globin amplification in the same PCR reaction. According to the results obtained by LA, samples were subdivided into HPV-negative (−) or HPV-positive (+), and into Single- or Multiple infection (SI or MI, respectively).

HPV sequencing by 454 NGS

All HPV-positive samples identified by conventional PCR utilizing PGMY09/11 primers were included for MS. Amplicons for the NGS Library were generated by taking approximately 500–1000 ng of DNA from each sample and amplification by conventional PCR employing the following specifically designed primers: Forward 5′-[Universal Multiplicom tail A: AAGACTCGGCAGCATCTCCA] – [individually PGMY11 primer]-3′, and Reverse 5′-[Universal Multiplicom tail B: GCGATCGTCACTGTTCTCCA] – [individually PGMY09 primer]-3′. The PGMY sequences used were those reported in Table 2 in Gravitt et al. [46]. Universal Multiplicom tails A and B were the same used in Cat. no. MR-0020.024, Multiplicom NV CFTR; Molecular Diagnostics, Niel, Belgium. PCR conditions included 1 cycle at 94 °C for 3 min, followed by 40 cycles at 94 °C for 45 s, 48 °C for 30 s, and 72 °C for 60 s, and a final cycle at 72 °C for 10 min. The PCR products (~450 pb) were electrophoresed in 1.5 % agarose gels, and cut and purified with the Wizard® SV Gel and PCR Clean-Up System (Cat. no. A9282; Promega, Madison, WI, USA).

Purified amplicons generated from each sample were quantified with the Qubit® dsDNA HS fluorometry assay Kit (Cat. no. Q32851; Life Technologies™, Eugene, OR, USA) to carry out 5-ng amplicon pools per sample. Amplicon pools were generated for the NGS Library, taking into account diagnosis, HPV negativity or positivity observed by LA, and the SI- or MI-HPV infection detected in a same sample. The following groups were generated: WCL (−), samples WCL-negative for HPV by LA; WCL (+), samples WCL-positive for HPV by LA; CIN I (−), CIN I samples negative for HPV by LA; CIN I SI, CIN I samples with single HPV infection; CIN I MI, CIN I samples infected with 2 or more HPVs, and CC SI and CC MI.

A second PCR reaction was performed to add the Multiplex identifiers (MID) barcodes, the key sequence (TCAG), and A or B sequencing adapters by using Multiplicom 454 MID Kits (Cat. no. ML-0008.192 and ML-0016.192, Multiplicom NV CFTR; Molecular Diagnostics) according to the manufacturer’s instructions. Amplicons where then purified with Agentcourt® AMPure® XP beads (Cat. no. A63881; Beckman Coulter Genomics, Danvers, MA, USA) and evaluated with the 2100 Bioanalyzer system (Cat. no. G2940CA; Agilent Technologies, Santa Clara, CA, USA).

The emulsion Polymerase chain reaction (emPCR) was performed with the GS Junior Titanium emPCR kit Lib A (Cat. no. 05996520001; Roche Diagnostics), and amplicons were sequenced in a GS Junior Sequencer utilizing the GS Junior Titanium Sequencing Kit (Cat. no. 05996554001) and the GS Junior Titanium PicoTiterPlate Kit (Cat. no. 0599661900), following the manufacturer’s instruction. All kits were supplied by Roche Diagnostics, Basel, Switzerland.

Data analysis

Quality control for sequencing data was performed, sequences with a Phred quality score of <20 were eliminated, and primer sequences were trimmed. The Roche GS Reference Mapper (ver. 2.9) was used to identify HPV genotypes by comparing fragment sequences with the L1 region of alpha and beta papillomavirus reference sequences, downloaded from Papillomavirus Episteme (PaVE) [48]. Specific reads for each diagnosis group were identified in accordance with their associated MID barcode.

Analysis of phylogenic ancestry

Construction of a phylogenetic tree with consensus sequences of each HPV was performed. First, all consensus sequences were aligned utilizing MEGA6 software (employing the ClustalW ver. 1.6 algorithm), followed by the construction of the phylogenetic tree with the maximum likelihood method. Traditional tree-branch style was chosen and the Topology display was selected.

Ethical approval

This study was approved by the Ethics and Research Committees of the Instituto Mexicano del Seguro Social (IMSS), with registration number R-2014-785-036. All women provided written informed consent before the sample collection.

Abbreviations

HPV:

Human papillomavirus

WCL:

Without cervical lesion

CIN I:

Cervical intraepithelial neoplasia grade I

CC:

Cervical cancer

LA:

Linear Array

SI:

Single infection

MI:

Multiple infection

PaVE:

Papillomavirus Episteme

NGS:

Next-Generation Sequencing

References

  1. Globocan. Estimated Cancer Incidence, Mortality and Prevalence Worldwide in 2012 IARC. 2012. http://globocan.iarc.fr/Default.aspx. Accessed 05-13-2015 2015.

  2. Zur Hausen H. Papillomavirus infections--a major cause of human cancers. Biochim Biophys Acta. 1996;1288(2):F55–78.

    PubMed  Google Scholar 

  3. HPVCenter. International Human Papillomavirus Reference Center. Davit Bzhalava, Karolinska Institutet. http://www.hpvcenter.se/html/refclones.html. Accessed February 11 2015.

  4. 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. J Natl Cancer Inst. 1995;87(11):796–802.

    Article  CAS  PubMed  Google Scholar 

  5. Clifford G, Smith J, Plummer M, Munoz N, Franceschi S. Human papillomavirus types in invasive cervical cancer worldwide: a meta-analysis. Br J Cancer. 2003;88(1):63–73.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, et al. A review of human carcinogens--Part B: biological agents. Lancet Oncol. 2009;10(4):321–2.

    Article  PubMed  Google Scholar 

  7. de Sanjose S, Quint WG, 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(11):1048–56. doi:10.1016/S1470-2045(10)70230-8.

    Article  PubMed  Google Scholar 

  8. Guan P, Howell-Jones R, Li N, Bruni L, de Sanjose S, Franceschi S, et al. Human papillomavirus types in 115,789 HPV-positive women: a meta-analysis from cervical infection to cancer. Int J Cancer. 2012;131(10):2349–59. doi:10.1002/ijc.27485.

    Article  CAS  PubMed  Google Scholar 

  9. Tjalma WA, Fiander A, Reich O, Powell N, Nowakowski AM, Kirschner B, et al. Differences in human papillomavirus type distribution in high-grade cervical intraepithelial neoplasia and invasive cervical cancer in Europe. Int J Cancer. 2013;132(4):854–67. doi:10.1002/ijc.27713.

    Article  CAS  PubMed  Google Scholar 

  10. Froberg M, Johansson B, Hjerpe A, Andersson S. Human papillomavirus 'reflex' testing as a screening method in cases of minor cytological abnormalities. Br J Cancer. 2008;99(4):563–8. doi:10.1038/sj.bjc.6604504.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Ferreccio C, Corvalan A, Margozzini P, Viviani P, Gonzalez C, Aguilera X, et al. Baseline assessment of prevalence and geographical distribution of HPV types in Chile using self-collected vaginal samples. BMC Public Health. 2008;8:78. doi:10.1186/1471-2458-8-78.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Wong OG, Lo CK, Chow JN, Tsun OK, Szeto E, Liu SS, et al. Comparison of the GenoFlow human papillomavirus (HPV) test and the Linear Array assay for HPV screening in an Asian population. J Clin Microbiol. 2012;50(5):1691–7. doi:10.1128/JCM.05933-11.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Schmidt-Grimminger DC, Bell MC, Muller CJ, Maher DM, Chauhan SC, Buchwald DS. HPV infection among rural American Indian women and urban white women in South Dakota: an HPV prevalence study. BMC Infect Dis. 2011;11:252. doi:10.1186/1471-2334-11-252.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Aguilar-Lemarroy A, Vallejo-Ruiz V, Cortes-Gutierrez EI, Salgado-Bernabe ME, Ramos-Gonzalez NP, Ortega-Cervantes L, et al. Human papillomavirus infections in Mexican women with normal cytology, precancerous lesions, and cervical cancer: type-specific prevalence and HPV coinfections. J Med Virol. 2015;87(5):871–84. doi:10.1002/jmv.24099.

    Article  PubMed  Google Scholar 

  15. Poljak M, Kocjan BJ. Commercially available assays for multiplex detection of alpha human papillomaviruses. Expert Rev Anti Infect Ther. 2010;8(10):1139–62. doi:10.1586/eri.10.104.

    Article  CAS  PubMed  Google Scholar 

  16. Barzon L, Militello V, Lavezzo E, Franchin E, Peta E, Squarzon L, et al. Human papillomavirus genotyping by 454 next generation sequencing technology. J Clin Virol. 2011;52(2):93–7. doi:10.1016/j.jcv.2011.07.006.

    Article  CAS  PubMed  Google Scholar 

  17. Arroyo LS, Smelov V, Bzhalava D, Eklund C, Hultin E, Dillner J. Next generation sequencing for human papillomavirus genotyping. J Clin Virol. 2013;58(2):437–42. doi:10.1016/j.jcv.2013.07.013.

    Article  CAS  PubMed  Google Scholar 

  18. Militello V, Lavezzo E, Costanzi G, Franchin E, Di Camillo B, Toppo S, et al. Accurate human papillomavirus genotyping by 454 pyrosequencing. Clin Microbiol Infect. 2013;19(10):E428–34. doi:10.1111/1469-0691.12219.

    Article  CAS  PubMed  Google Scholar 

  19. Schiffman M, Wentzensen N, Wacholder S, Kinney W, Gage JC, Castle PE. Human papillomavirus testing in the prevention of cervical cancer. J Natl Cancer Inst. 2011;103(5):368–83.

  20. 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. doi:10.1002/ijc.22527.

    Article  CAS  PubMed  Google Scholar 

  21. Anderson L, O'Rorke M, Jamison J, Wilson R, Gavin A. members HPVWG. Prevalence of human papillomavirus in women attending cervical screening in the UK and Ireland: new data from northern Ireland and a systematic review and meta-analysis. J Med Virol. 2013;85(2):295–308. doi:10.1002/jmv.23459.

    Article  PubMed  Google Scholar 

  22. Quek SC, Lim BK, Domingo E, Soon R, Park JS, Vu TN, et al. Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical intraepithelial neoplasia across 5 countries in Asia. Int J Gynecol Cancer. 2013;23(1):148–56. doi:10.1097/IGC.0b013e31827670fd.

    Article  PubMed  Google Scholar 

  23. Beaudenon S, Praetorius F, Kremsdorf D, Lutzner M, Worsaae N, Pehau-Arnaudet G, et al. A new type of human papillomavirus associated with oral focal epithelial hyperplasia. J Invest Dermatol. 1987;88(2):130–5.

    Article  CAS  PubMed  Google Scholar 

  24. Falaki F, Amir Chaghmaghi M, Pakfetrat A, Delavarian Z, Mozaffari PM, Pazooki N. Detection of human papilloma virus DNA in seven cases of focal epithelial hyperplasia in Iran. J Oral Pathol Med. 2009;38(10):773–6. doi:10.1111/j.1600-0714.2009.00784.x.

    Article  CAS  PubMed  Google Scholar 

  25. Volter C, He Y, Delius H, Roy-Burman A, Greenspan JS, Greenspan D, et al. Novel HPV types present in oral papillomatous lesions from patients with HIV infection. Int J Cancer. 1996;66(4):453–6. doi:10.1002/(SICI)1097-0215(19960516)66:4<453::AID-IJC7>3.0.CO;2-V.

    Article  CAS  PubMed  Google Scholar 

  26. Kim NR, Kang M, Lee SP, Kim H, An J, Chung DH, et al. Uncommon and rare human papillomavirus genotypes relating to cervical carcinomas. Korean journal of pathology. 2014;48(1):43–9. doi:10.4132/KoreanJPathol.2014.48.1.43.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Caldeira S, Dong W, Tomakidi P, Paradiso A, Tommasino M. Human papillomavirus type 32 does not display in vitro transforming properties. Virology. 2002;301(1):157–64.

    Article  CAS  PubMed  Google Scholar 

  28. Lorincz AT, Quinn AP, Goldsborough MD, Schmidt BJ, Temple GF. Cloning and partial DNA sequencing of two new human papillomavirus types associated with condylomas and low-grade cervical neoplasia. J Virol. 1989;63(6):2829–34.

    PubMed Central  CAS  PubMed  Google Scholar 

  29. 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. Int J Cancer. 2014;134(6):1389–98.

    Article  CAS  PubMed  Google Scholar 

  30. Guimera N, Lloveras B, Lindeman J, Alemany L, van de Sandt M, Alejo M, et al. The occasional role of low-risk human papillomaviruses 6, 11, 42, 44, and 70 in anogenital carcinoma defined by laser capture microdissection/PCR methodology: results from a global study. Am J Surg Pathol. 2013;37(9):1299–310. doi:10.1097/PAS.0b013e31828b6be4.

    Article  PubMed  Google Scholar 

  31. Sandri MT, Riggio D, Salvatici M, Passerini R, Zorzino L, Boveri S, et al. Typing of human papillomavirus in women with cervical lesions: prevalence and distribution of different genotypes. J Med Virol. 2009;81(2):271–7. doi:10.1002/jmv.21382.

    Article  PubMed  Google Scholar 

  32. Stevens MP, Garland SM, Tan JH, Quinn MA, Petersen RW, Tabrizi SN. HPV genotype prevalence in women with abnormal pap smears in Melbourne, Australia. J Med Virol. 2009;81(7):1283–91. doi:10.1002/jmv.21515.

    Article  PubMed  Google Scholar 

  33. Nielsen A, Iftner T, Norgaard M, Munk C, Junge J, Kjaer SK. The importance of low-risk HPV infection for the risk of abnormal cervical cytology/histology in more than 40 000 Danish women. Sex Transm Infect. 2012;88(8):627–32. doi:10.1136/sextrans-2011-050307.

    Article  PubMed  Google Scholar 

  34. Petry KU, Luyten A, Justus A, Iftner A, Strehlke S, Schulze-Rath R, et al. Prevalence of low-risk HPV types and genital warts in women born 1988/89 or 1983/84 -results of WOLVES, a population-based epidemiological study in Wolfsburg, Germany. BMC infectious diseases. 2012;12:367. doi:10.1186/1471-2334-12-367.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Longuet M, Cassonnet P, Orth G. A novel genital human papillomavirus (HPV), HPV type 74, found in immunosuppressed patients. J Clin Microbiol. 1996;34(7):1859–62.

    PubMed Central  CAS  PubMed  Google Scholar 

  36. de Villiers EM. Cross-roads in the classification of papillomaviruses. Virology. 2013;445(1–2):2–10. doi:10.1016/j.virol.2013.04.023.

    Article  PubMed  Google Scholar 

  37. Giuliano AR, Lazcano-Ponce E, Villa LL, Flores R, Salmeron J, Lee JH, et al. The human papillomavirus infection in men study: human papillomavirus prevalence and type distribution among men residing in Brazil, Mexico, and the United States. Cancer Epidemiol Biomarkers Prev. 2008;17(8):2036–43. doi:10.1158/1055-9965.EPI-08-0151.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Meiring TL, Salimo AT, Coetzee B, Maree HJ, Moodley J, Hitzeroth II, et al. Next-generation sequencing of cervical DNA detects human papillomavirus types not detected by commercial kits. Virol J. 2012;9:164. doi:10.1186/1743-422X-9-164.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Chen Z, Schiffman M, Herrero R, Burk RD. Identification and characterization of two novel human papillomaviruses (HPVs) by overlapping PCR: HPV102 and HPV106. J Gen Virol. 2007;88(Pt 11):2952–5. doi:10.1099/vir.0.83178-0.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Ekstrom J, Forslund O, Dillner J. Three novel papillomaviruses (HPV109, HPV112 and HPV114) and their presence in cutaneous and mucosal samples. Virology. 2010;397(2):331–6. doi:10.1016/j.virol.2009.11.027.

    Article  PubMed  Google Scholar 

  41. de Sanjose S, Alemany L, Ordi J, Tous S, Alejo M, Bigby SM, et al. Worldwide human papillomavirus genotype attribution in over 2000 cases of intraepithelial and invasive lesions of the vulva. Eur J Cancer. 2013;49(16):3450–61. doi:10.1016/j.ejca.2013.06.033.

    Article  PubMed  Google Scholar 

  42. Bernard HU. Coevolution of papillomaviruses with human populations. Trends Microbiol. 1994;2(4):140–3.

    Article  CAS  PubMed  Google Scholar 

  43. De Villiers EM, Fauquet C, Broker TR, Bernard HU, Zur Hausen H. Classification of papillomaviruses. Virology. 2004;324(1):17–27. doi:10.1016/j.virol.2004.03.033.

    Article  PubMed  Google Scholar 

  44. IARC. International Agency for Research on Cancer. http://www.iarc.fr/en/about/index.php. 2015. Accessed May 2015.

  45. Lichtig H, Algrisi M, Botzer LE, Abadi T, Verbitzky Y, Jackman A, et al. HPV16 E6 natural variants exhibit different activities in functional assays relevant to the carcinogenic potential of E6. Virology. 2006;350(1):216–27. doi:10.1016/j.virol.2006.01.038.

    Article  CAS  PubMed  Google Scholar 

  46. Gravitt PE, Peyton CL, Alessi TQ, Wheeler CM, Coutlee F, Hildesheim A, et al. Improved amplification of genital human papillomaviruses. J Clin Microbiol. 2000;38(1):357–61.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Coutlee F, Rouleau D, Petignat P, Ghattas G, Kornegay JR, Schlag P, et al. Enhanced detection and typing of human papillomavirus (HPV) DNA in anogenital samples with PGMY primers and the Linear array HPV genotyping test. J Clin Microbiol. 2006;44(6):1998–2006. doi:10.1128/JCM.00104-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. PAVE. Papillomavirus Episteme. National Institute of Allergy and Infectious Diseases and National Institutes of Health. 2015. http://pave.niaid.nih.gov/. 2015.

Download references

Acknowledgments

MGF-M, LAT-R, and LA-R are grateful for the scholarship obtained from CONACYT-México.

Funding

This study was supported by Fondo de Investigación en Salud-IMSS (FIS/IMSS/PROT/PRIO/14/033 granted to AA-L).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Adriana Aguilar-Lemarroy.

Additional information

Competing interest

All of the authors declare no potential conflict of interest.

Authors’ contributions

MGF-M extracted DNA and performed PCR reactions. MGF-M and LA-R performed the LA test, and MGF-M, LAT-R, and VR-R conducted NGS. MGF-M, LAT-R, and SAR-M carried out the bioinformatics analyses; LMB-P was involved in patient interviews and sample recruitment. PP-S, VV-R, and EIC-G also participated in the general study design and critically revised the manuscript. AA-L and LFJ-S conceived of and designed the study, supervised all experiments and analyses, and wrote the manuscript. All authors read and approved the final manuscript.

Luis Felipe Jave-Suárez and Adriana Aguilar-Lemarroy contributed equally to this work.

Additional file

Additional file 1:

Table S1. Human papillomavirus (HPV) genotypes identified by Linear Array (LA). The Table shows the HPV genotypes found in each sample analyzed and grouped according to diagnosis. WCL (Without cervical lesion), CIN I (Cervical intraepithelial neoplasia grade I), CC (Cervical cancer), SI (Single infection), and MI (Multiple infection). (PDF 20 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

Flores-Miramontes, M.G., Torres-Reyes, L.A., Alvarado-Ruíz, L. et al. Human papillomavirus genotyping by Linear Array and Next-Generation Sequencing in cervical samples from Western Mexico. Virol J 12, 161 (2015). https://doi.org/10.1186/s12985-015-0391-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-015-0391-4

Keywords