- Research article
- Open Access
Prevalence of high-risk human papillomavirus types in Mexican women with cervical intraepithelial neoplasia and invasive carcinoma
Infectious Agents and Cancer volume 3, Article number: 3 (2008)
Prevalence of high risk (HR) human papillomavirus (HPV) types in the states of San Luis Potosí (SLP) and Guanajuato (Gto), Mexico, was determined by restriction fragment length-polymorphism (RFLP) analysis on the E6 ~250 bp (E6-250) HR-HPV products amplified from cervical scrapings of 442 women with cervical intraepithelial neoplasia and invasive carcinoma (280 from SLP and 192 from Gto). Fresh cervical scrapings for HPV detection and typing were obtained from all of them and cytological and/or histological diagnoses were performed on 383.
Low grade intraepithelial squamous lesions (LSIL) were diagnosed in 280 cases (73.1%), high grade intraepithelial squamous lesions (HSIL) in 64 cases (16.7%) and invasive carcinoma in 39 cases (10.2%). In the 437 cervical scrapings containing amplifiable DNA, only four (0.9%) were not infected by HPV, whereas 402 (92.0%) were infected HR-HPV and 31 (7.1%) by low-risk HPV. RFLP analysis of the amplifiable samples identified infections by one HR-HPV type in 71.4%, by two types in 25.9% and by three types in 2.7%. The overall prevalence of HR-HPV types was, in descending order: 16 (53.4%) > 31 (15.6%) > 18 (8.9%) > 35 (5.6) > 52 (5.4%) > 33 (1.2%) > 58 (0.7%) = unidentified types (0.7%); in double infections (type 58 absent in Gto) it was 16 (88.5%) > 31 (57.7%) > 35 (19.2%) > 18 (16.3%) = 52 (16.3%) > 33 (2.8%) = 58 (2.8%) > unidentified types (1.0%); in triple infections (types 33 and 58 absent in both states) it was 16 (100.0%) > 35 (54.5%) > 31 (45.5%) = 52 (45.5%) > 18 (27.3%). Overall frequency of cervical lesions was LSIL (73.1%) > HSIL (16.7%) > invasive cancer (10.2%). The ratio of single to multiple infections was inversely proportional to the severity of the lesions: 2.46 for LSIL, 2.37 for HSIL and 2.15 for invasive cancer. The frequency of HR-HPV types in HSIL and invasive cancer lesions was 16 (55.0%) > 31 (18.6%) > 35 (7.9%) > 52 (7.1%) > 18 (4.3%) > unidentified types (3.6%) > 33 (2.9%) > 58 (0.7%).
Ninety percent of the women included in this study were infected by HR-HPV, with a prevalence 1.14 higher in Gto. All seven HR-HPV types identifiable with the PCR-RFLP method used circulate in SLP and Gto, and were diagnosed in 99.3% of the cases. Seventy-one percent of HR-HPV infections were due to a single type, 25.9% were double and 2.7% were triple. Overall frequency of lesions was LSIL (73.1%) > HSIL (16.7%) > invasive cancer (10.2%), and the ratio of single to multiple infections was inversely proportional to severity of the lesions: 2.46 for LSIL, 2.37 for HSIL and 2.15 for invasive cancer. The frequency of HR-HPV types found in HSIL and invasive cancer was 16 (55.0%) > 31 (18.6%) > 35 (7.9%) > 52 (7.1%) > 18 (4.3%) > unidentified types (3.6%) > 33 (2.9%) > 58 (0.7%). Since the three predominant types (16, 31 and 18) cause 77.9% of the HR-HPV infections and immunization against type 16 prevents type 31 infections, in this region the efficacy of the prophylactic vaccine against types 16 and 18 would be close to 80%.
Cervical cancer (CC) is the second cause of death by cancer among women in the world and the first in most developing countries [1, 2]. In the year 2000, 190,000 deaths and 80% of the 500,000 new CC cases occurred in the developing world, and Latin American countries were among those with the highest incidence rates, together with countries from Sub-Saharan Africa, South and South East Asia . The risk of developing CC increases with early start of sexual activity, number of sexual partners, prolonged use of oral contraceptives and smoking . In Mexico it is associated with poverty related factors such as low schooling, unemployment, residence in rural areas and lack of access to health services [5–7].
zur Hausen  proposed and later demonstrated  that human papillomaviruses (HPV) are the infectious agents responsible of the neoplastic transformation of the cervical epithelium. This hypothesis was validated by finding HPV genome sequences in 99.7% of invasive CC cases . It is currently accepted that HPV is the most common sexually transmitted pathogen  and that infection of the cervix by 'high risk' (HR) HPV types is a necessary factor for CC development [10, 12]. Low risk (LR) HPV types are those usually found in warts and benign lesions whereas HR types are those found in invasive CC .
There are 11 major LR-HPV types (designated as 6, 11, 40, 42, 43, 44, 54, 61, 70, 72, 81) and 15 major HR-HPV types (designated as 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68, 73, 82) . The two most frequently associated to malignancies are HR-HPV types 16 and 18. The first is responsible for nearly 50% and the second for nearly 20% of all invasive CC cases in the world [10, 15]. Several HR-HPV variants with higher oncogenic potential are more prevalent in developing countries, where they appear to contribute to the higher incidence and mortality rates. In Mexico the probability of developing CC by the Asian-American variant of HPV type 16 is several times higher than by the European variant, and around 25% of invasive CC cases are attributed to the AA variant .
During the year 2000 the incidence of CC in Mexico was 40.5 cases per 100,000 women over 25 years of age and the mortality rate was 17.1 , which means 12 deaths per day. In the Mexican states of San Luis Potosí (SLP) and Guanajuato (Gto) the mortality rates for the same year were respectively 15.5 and 21.7 . Lazcano-Ponce et al.  found a HPV prevalence of 16.7% in Mexican women below 25 years of age which decreased to 3.7% in those of 25–44 years and increased to 23% in those above 65 years, and that HR-HPV types are more prevalent in all age groups whereas LR-HPV types are less prevalent in women under 25 years and increase with age.
In the present work we determined the prevalence of the HR-HPV types circulating in the neighboring states of San Luis Potosí and Guanajuato by restriction fragment length polymorphism (RFLP) analysis with the method of Fujinaga et al. , in women with precancerous and cancerous lesions of the cervix that had been diagnosed by cytological analysis; in two thirds of them biopsies or surgical specimens were obtained and the microscopic diagnosis confirmed.
The 442 women enrolled in this study ranged from 16 to 78 years of age. Fresh cervical scrapings for HPV detection and typing were obtained from all of them whereas cytological and/or histological diagnoses were performed on 383. Low grade intraepithelial squamous lesions (LSIL) were diagnosed in 280 cases (73.1%), high grade intraepithelial squamous lesions (HSIL) in 64 cases (16.7%) and invasive carcinoma in 39 cases (10.2%).
Algorithm for HPV DNA amplification of cervical scrapings
Amplification of the HR-HPV E6-250 product by direct and nested PCR was compared. Direct PCR was performed with primer set pU 1M/2R whereas nested PCR was performed in two successive reactions: in the first one ('direct') the E6-650 segment was amplified using primer set LCRS/E7AS, whereas in the second one ('nested') a portion of the first reaction mixture (presumed to contain the E6-650 product) was used as template for the internal primers (pU 1M/2R) to generate E6-250.
Amplification of 35 random samples showed that 11 (31.4%) were positive by direct PCR and 24 (68.6%) by nested PCR. Since E6-250 detection was 2.2 times more sensitive with nested PCR, this modality was chosen to analyze all problem samples.
To determine the presence of HPV DNA in cervical scrapings the algorithm shown in Fig. 1 was used. Each sample is first subjected to nested PCR to successively amplify the E6-650 and E6-250 products. If the second one is visible, the sample is considered positive for HR-HPV (examples in Fig. 2A) and typified by RFLP (examples in Fig. 2B and 3C). If the second product is not observed, the sample could be either positive, inadequate or negative. To define if it is positive it is subjected to PCR with the primer set MY 09/11; if L1-450 is observed, the sample is considered positive for low risk HPV and if L1-450 is not visible, the sample may be negative or of inadequate quality and subjected to PCR with the β-globin primer set. If the β-globin product is observed the sample is considered to be amplifiable but HPV-negative; if it is not visible, the sample is considered inadequate (i.e., not amplifiable) and undefined in relation to HPV infection.
In the analysis of samples from women with cervical dysplasia and neoplasia, as is the case of this study, this algorithm saves labor, time and reagents because it first verifies the presence of HPV through nested PCR (expected to be positive in most samples) and then the amplificability of the HPV-negative samples (expected to be rare) by performing a PCR assay for a single copy housekeeping gene such as that of β-globin.
Using the algorithm for HPV detection (Fig. 1) the 437 cervical scrapings with amplifiable DNA were analyzed. HPV infection was demonstrated in 433 samples (99.1%) of which 402 (92.0%) corresponded to high risk viral types (Table 1).
Algorithm for HR-HPV typing
The sizes of E6-250 products and their corresponding restriction fragments for the seven HR-HPV types that can be identified with the PCR-RFLP method used appear in Table 2.
To identify HR-HPV types the algorithm of Fig. 3 was used. The E6-250 product from each HR-HPV type has a specific size and a single restriction site for each of the five enzymes used, which generate two fragments of specific sizes allowing the identification of seven of the most frequent HR-HPV types. Ava II, the first enzyme used on each HR-HPV positive sample, identifies types 16, 18 and 33 (see examples in Fig. 2A). If Ava II does not digest the E6-250 product or does it only partially, sample aliquots are simultaneously incubated with each of the four remaining enzymes to identify types 31, 52, 35 and 58 (see examples in Fig. 2C).
HR-HPV types circulating in the region
In 402 of the 435 samples with HPV infection, one to three of the seven identifiable HR-HPV types were detected.
The overall prevalence of HR-HPV types was, in descending order: 16 (53.4%) > 31 (15.6%) > 18 (8.9%) > 35 (5.6) > 52 (5.4%) > 33 (1.2%) > 58 (0.7%) = unidentified types (0.7%). Types 16, 31 and 18 therefore accounted for 77.9% of all the infections whereas the other four types accounted for only 12.9% of them. HPV type 16 predominated slightly in SLP (53.3%) over Gto (52.9%), type 31 in Gto (22.3%) over SLP (10.1%), and type 18 was equally prevalent in SLP and Gto (8.8%) (Table 3).
HR-HPV types in single and multiple infections
The overall prevalence of infections by a single HR-HPV type was 71.4%, i.e., 2.5 times higher than that of infections by two and three types (Table 4).
The prevalence of infections by a single HR-HPV type was slightly higher in SLP (75.0%) than in Gto (67.4%), whereas the overall prevalence of infections by two and three HR-HPV types was respectively 25.9% and 2.7%, of which 22.6% and 2.4% corresponded to SLP and 29.5% and 3.2% to Gto (Table 4).
The HR-HPV types of single, double and triple infections as well as the combinations of viral types in multiple infections found in all samples and per state are shown in Table 5.
HPV type 16 was the most prevalent (72.6%) followed by types 18 (10.2%) and 31 (9.9%) and then by types 52 (3.3%), 33 (1.3%) = 35 (1.3%), and 58 (0.3%) (Table 5).
The descending order of frequency of the viral types in single infections was clearly different in each state. In SLP the type 16 (77.4%) predominated, followed by types 18 (9.4%), 31 (8.2%), 52 (2.5%), 33 (1.3%), 35 (0.6%) = 58 (0.6%). In Gto the prevalence of type 16 was lower (63.3%) and type 31 was second (13.3%), followed by types 18 (12.5%), 52 (4.7%), 35 (2.3%) and 33 (1.6%) (Table 5).
The overall prevalence of HPV type 16 increased in double infections and reached its maximum possible value (100%) in triple infections (Table 6). The overall prevalence of HPV type 18 also increased in double infections and was accompanied by increases of types 31, 52 and 35 and the last two types were even higher in triple infections (Table 6). In double and triple infections of SLP the predominant types were 16, 31, 18, 35 and 52 whereas in Gto types 16 and 31 predominated in double infections and types 16 and 31 in the triple ones (Table 6).
Single and multiple infections and HR-HPV types in lesions of increasing severity
The overall frequency of cervical lesions was, in decreasing order, LSIL (73.1%) > HSIL (16.7%) > invasive cancer (10.2%) (Table 7). The ratio of single to multiple infections was inversely proportional to the severity of the lesions: 2.46 for LSIL, 2.37 for HSIL and 2.15 for invasive cancer (Table 8).
The frequency of HR-HPV types found in the cervical HSIL and invasive cancer (n = 140) was, in decreasing order, 16 (55.0%) > 31 (18.6%) > 35 (7.9%) > 52 (7.1%) > 18 (4.3%) > unidentified types (3.6%) > 33 (2.9%) > 58 (0.7%) (Table 9).
Since our aim was to determine the prevalence of the HR-HPV types causing cervical infections in the states of SLP and Gto, the study included women with cervical intraepithelial lesions (73.1% low grade, 16.7% high grade and 10.2% invasive carcinoma), expected to be HPV-positive. This prediction was fulfilled because nearly all women were shown to be infected with HPV, and HR types were identified in most of them.
To determine HPV infection and to identify the high-risk types involved we used the PCR-RFLP method based on amplification of the E6-250 HR-HPV product followed by digestion with endonucleases  and fragment sizing in high resolution agarose gels . Obtaining cervical scrapings with the cytobrush allowed us to obtain more than enough amplifiable DNA in nearly all samples and PCR optimization assured efficient amplification with the primer sets employed.
The E6-250 HR-HPV product was amplified through nested PCR because this modality turned out to be at least twice as sensitive as direct PCR. The universal primer sets for the first ('direct') reaction amplified the E6-650 product [19, 26] and the primer set used in the second ('nested') reaction amplified the E6-250 product . The use of duplex PCR was discarded to amplify E6-250 and L1-450 due to the interference observed when simultaneous amplification of both products was attempted (Fig. 4). The scarce E6-negative samples were tested later with a universal primer set that amplifies L1-450 segments of both low risk- as well as HR-HPV types . From these observations we concluded that to determine HR-HPV infection in women with dysplasia and neoplasia it is convenient first to identify the HPV-positive samples and then to check the amplificability of the negative ones.
Eight HR-HPV types (16, 18, 31, 33, 35, 45, 52 and 58) are responsible for 95% of the invasive CC cases in the world . The PCR-RFLP used in this work identified all of them except type 45.
Considering all cervical samples included in this work, 98.8% contained amplifiable DNA, 97.8% had HPV infections which in 90.2% of the cases corresponded to HR-HPV types whose prevalence was higher in Gto (99.0%) than in SLP (82.6%) (Table 3).
In the 437 samples with amplifiable DNA included in this study, 528 cervical infections were demonstrated by HR-HPV types that were identified in 99.1% of the cases (Table 4). In this way we demonstrated that in our region circulate the seven HR-HPV types identified with the method used and their overall prevalence is, in descending order, 16 (53.4%) > 31 (15.6%) > 18 (8.9%) > 35 (5.6) > 52 (5.4%) > 33 (1.2%) > 58 (0.7%) = unidentified types (0.7%). From these results we recommend that for analysis of cervical samples with dysplasia it is convenient to identify first the types 16, 18 and 33 by restriction with Ava II and then the types 31, 35, 52 and 58 by restriction with Rsa I, Ava I, Bgl II and Acc I.
The prevalence of high and low risk HPV types found by us is similar to that already seen in Mexico and other regions of high CC incidence [20, 21, 29]. The predominant HR-HPV types in our region are 16, 31 and 18 which together represent 77.9% of all cervical infections. The four remaining types cause 12.9%, and unidentified types only 0.7% of the HR-HPV infections. Types 16 and 18 are slightly more prevalent in SLP, whereas type 31 is nearly two times more prevalent in Gto.
The overall prevalence of HPV type 16 in our region (53.4%) is almost identical to the world average (53%) calculated by Munoz (2000). In contrast, the regional prevalence of type 18 (8.9%) is substantially lower than the world average (15%) calculated by Munoz . The prevalence of types 31 (15.6%) and 35 (5.6%) is higher than the average for Latin America , whereas that of type 33 (1.2%) is nearly 10 times lower than the average for Central America (11.8%) recently obtained by Clifford et al. .
Single and triple infections were more frequent in SLP, whereas double infections were more frequent in Gto. Single infections included all seven HR-HPV identifiable types where the three predominant ones were, in descending order, 16 > 18 > 31 with a cumulative frequency of 92.7%. The prevalence of these types in SLP was, in descending order, 16 > 18 > 31 with a cumulative frequency of 95.0%, whereas in Gto it was 16 > 31 > 18 with a cumulative frequency of 89.1%.
The typing method allowed us to identify a considerable proportion of infections by two viral types (overall prevalence of 25.9%, with 22.6% for SLP and 29.5% for Gto) and three viral types (overall prevalence of 2.7%, with 2.4% for SLP and 3.2% for Gto) (Table 4). Type 16 was found in all multiple infections.
Double infections from SLP were caused by the seven identifiable types but type 58 was absent in Gto. The overall predominant types in double infections were 16 > 31 > 35 > 18 = 52. Their prevalence differed in each state since the descending order in SLP was 16 > 31 > 18 > 35 > 52, whereas in Gto it was 16 > 31 > 35 > 52 > 18 (Table 6).
Types 33 and 58 were absent from triple infections in both states and the overall predominant types were 16 > 35 > 31 = 52 > 18. There were also differences in the prevalence by state since in SLP the descending order was 16 > 35 = 52 > 31 > 18 whereas in Gto it was 16 > 35 = 52 > 31 > 18 (Table 6).
Our algorithm for typing, designed to identify first the high-risk viral types and then the low risk ones, probably underestimates the prevalence of low risk types in samples where high risk types are demonstrated. However, it is worth noting that in low grade intraepithelial cervical lesions practically all genital HPV types are found, whereas in high grade lesions only high risk types are found .
The proportion of multiple infections found by us (28.6%) is over three times higher than the world average (8.1%) calculated by Molano et al. . This feature is probably due to the fact that the women in our study are themselves a high risk group since they were included for having cervical dysplasia and cancer. This idea is reinforced by our finding that 10.6% of the double infections are due to the 16/18 type couple, whose prevalence was nearly six times higher than the world average (1.8%) estimated by Clifford et al. .
The five predominant HR-HPV types (16, 31, 18, 35 and 52) found in the 140 HSIL and invasive cancer lesions accounted for 92.9% of the cases and their frequency was similar to their overall prevalence (Table 9).
Infections by two or more HPV types are markers of persistent cervical disease, of multiple cervical lesions and of progression from low to high grade lesions . There is scarce evidence on the interaction of HPV types to enhance CC pathogenesis and it is unknown if some HPV types promote or exclude infection by others. Careful follow up of multiple infections may help to pinpoint the relevance of the interactions of different HPV types in disease progression .
HPV16, the most frequent type in single and multiple infections in SLP and Gto, just as in the rest of the world [14, 21] appears to be the most persistent type [32, 35, 36] increasing the risk of developing severe preneoplastic lesions .
Our work confirms the usefulness of the PCR-RFLP method of Fujinaga et al. [19, 26, 27] for diagnosis and typing cervical infections by HR-HPV, whose sensitivity was at least doubled by amplifying the E6-250 product through nested rather than direct PCR and whose typing precision was improved by using high resolution electrophoresis. In this way we could identify the HR-HPV types circulating in our region and their distribution in the two neighboring states analyzed, which reflect the endemic behavior of the various HPV types as well as local environmental or genetic conditions of the hosts .
Besides being of use to carry out molecular epidemiology studies and to predict the individual risk of each patient, HR-HPV typing is essential to plan the administration of prophylactic vaccines directed against one or more specific HPV types, which is questionable in populations where prevalence of the circulating types is unknown [31, 39, 40]. Since in our region HPV type 16 has the highest overall prevalence (53.4%), type 31 is second (15.6%) and type 18 is third (8.9%), the vaccine against HR-HPV types 16 and 18  whose license was granted first in Mexico and the United States [42, 43] would be expected to cover at least these three types  and would have a 78% efficacy.
Enough DNA is obtained from cervical scrapings of women with dysplasia in order to detect high risk-HPV by PCR amplification followed by typing through restriction fragment length polymorphism analysis.
Sensitivity of HR-HPV detection is nearly doubled by nested rather than direct PCR amplification.
Ninety-two percent of the 437 Mexican women with cervical dysplasia and cancer included in this study were infected by HR-HPV, with a prevalence 1.14 times higher in the state of Guanajuato than in San Luis Potosí.
All seven HR-HPV types identifiable with the method used (16, 18, 31, 33, 35, 52, 58) circulate in SLP and Gto and were diagnosed in nearly all (99.3%) the cases.
Seventy-one percent of the infections were caused by one HR-HPV type, 25.9% by two types and 2.7% by three types.
The overall frequency of cervical lesions was LSIL (73.1%) > HSIL (16.7%) > invasive cancer (10.2%) and the ratio of single to multiple infections was inversely proportional to the severity of the lesions: 2.46 for LSIL, 2.37 for HSIL and 2.15 for invasive cancer.
The frequency of HR-HPV types found in HSIL and invasive cancer was 16 (55.0%) > 31 (18.6%) > 35 (7.9%) > 52 (7.1%) > 18 (4.3%) > unidentified types (3.6%) > 33 (2.9%) > 58 (0.7%).
Since the three predominant HR-HPV types found (16, 31 and 18) cause 77.9% of the HR-HPV infections associated to cervical dysplasia in this region, and immunization against type 16 prevents type 31 infections, in this region the efficacy of the prophylactic vaccine against types 16 and 18 would be close to 80%.
Source of women
The 442 women included in this study were selected for having precancerous or cancerous lesions of the cervical epithelium confirmed by cytologic analysis. Two hundred and eighteen resided in the city of San Luis Potosí, capital of the state of San Luis Potosí (SLP), and 192 in the cities of León, Celaya and Irapuato from the state of Guanajuato (Gto), Mexico. Cervical scrapings from SLP were obtained at the Dysplasia Clinic of the SLP Health Services, the School of Nursing of Universidad Autónoma de San Luis Potosí, and from private patients. Samples from Gto were obtained at the Dysplasia Clinics of the cities of Irapuato, Celaya and León and sent to our laboratory by the Guanajuato Health Secretariat. The study was authorized by the Health Secretariats of SLP and Gto and performed with the informed consent of all participating women.
Cervical scrapings and DNA extraction
Each scraping was taken with an endocervical brush ('cytobrush') that was immediately inserted in a 5 ml polypropylene tube (Nalge Nunc, Rochester, NY) containing 1 ml phosphate buffered saline supplemented with sterile disodium ethtylene-diamino-tetraacetate (PBS-EDTA: 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 2 mM KH2HPO2, 25 mM disodium EDTA, pH 7.4). Once detached from the cytobrush and suspended in the PBS-EDTA vehicle, each sample was fixed by addition of 1.5 ml 96% ethanol and processed to extract DNA on the same day or up to 30 days after being kept at room temperature. Reagents were purchased from J.T. Baker (Xalostoc, México) unless other source is specified.
To extract the DNA each fixed sample was mixed by vortexing and 1 ml transfered to a 1.5 ml tube and spun in a Hettich Mikro 20 microcentrifuge (Cologne, Germany) for 5 min at 13,000 rpm (16,250×g). The supernatant was discarded by decantation and to each pellet were added 500 μl of Tris-EDTA-saline (TES: 10 mM Tris-HCl; 2 mM disodium EDTA, 0.4 M NaCl, pH 8.0 at 25°C), 50 μl of 10% sodium dodecyl sulphate and 20 μl of proteinase K (20 mg/ml). Mixtures were incubated at 55°C for 3 h, at the end of which 150 μl of 5 M NaCl were added and centrifuged again for 15 min. Each supernatant was aspirated and transferred to a tube to which 577 μl cold isopropanol were immediately added and then left stand for 10 min at 4°C to precipitate the nucleic acids. The tubes were centrifuged again for 10 min and supernatants discarded by decantation. Each pellet was washed by vortexing with 1 ml of 70% cold ethanol and centrifuged for 10 min at 10,000 rpm (9,615×g) and room temperature. Supernatants were discarded by aspiration and the pellets dried out by inverting the tubes for 15 min on a paper towel. Each pellet was dissolved with 50 μl TE (10 mM Tris-HCl, 1 mM disodium EDTA, pH 8.0 at 25°C).
The quality of extracted DNA was verified by electrophoretic analysis in 1% agarose gels with TAE (40 mM Tris-acetate, 1 mM disodium EDTA, pH 8.2 at 25°C). Two-μl from each sample were applied to gels which were run at 60 V for 90 min. λ-phage? DNA digested with Hind III (Sigma-Aldrich, Mexico) was used as marker. After staining for 20 min with ethidium bromide (1 μg/ml) gels were illuminated with ultraviolet light and their fluorescence recorded with the Bio-Rad ChemiDoc EQ (Hercules, CA) photo documenter.
DNA was quantitated by fluorometry with the PicoGreen dsDNA Quantitation kit (Molecular Probes; Eugene, OR) by interpolation in a standard curve containing up to 50 ng of λ-phage DNA. To each well of a black FIA 96 well plate (Greiner Bio-One, Frickenhausen, Germany) 198 μl of the assay solution (PicoGreen diluted 1:400 in TE) and 2 μl of standard DNA or problem samples were added, and their fluorescence determined using a 485 nm excitation filter and a 535 emission filter in the GENios Pro fluorometer (Tecan Systems, San Jose, CA) with the Magellan 4 software.
Since amplification of the β-globin gene required = 25 ng DNA per PCR mixture, 100 ng was the minimum amount of DNA required from each cervical scraping in order to prepare at least four PCR mixtures.
The range of pure DNA obtained from the first 370 cervical scrapings was 53–16,590 ng. The amount of purified DNA was sufficient (i.e., = 100 ng) in all samples except five of the initial ones.
The oligonucleotide primer sets and the sizes expected for the HPV amplification products are depicted in Table 10. To optimize PCR with each primer set variable DNA and MgCl2 concentrations and annealing temperatures were tested. The β-globin primer set was added to PCR mixtures containing double-serial dilutions of DNA from a HPV-positive sample (range: 0.78–50 ng). To optimize conditions for HPV primer sets, 1–5 mM MgCl2 concentrations and annealing temperatures from 55°C to 61.4°C were tested.
To maximize HPV detection through PCR amplification, three variables were optimized: DNA content, annealing temperature and magnesium concentration.
With the PC04/GH20 primer set for β-globin the minimum amount of template DNA required in PCR mixtures was determined with double serial dilutions leading to a maximum of 50 ng DNA per mixture. Fluorescence of the E6-250 product increased with DNA content up to 25 ng, which was the amount selected to amplify HPV DNA from each sample. Annealing temperature in the range of 55°C to 62°C was assessed with each primer set. As major products attained maximum fluorescence at 57°C this annealing temperature was used with all primer sets for the rest of the work. MgCl2 concentrations in the 1–5 mM range were tested. Maximum efficiency was attained at 2.5 mM, concentration chosen for the rest of the work.
The primer set LCRS/E7AS amplifies the E6-650 segment and the set pU 1M/2R amplifies the 'nested' E6-250 segment included in the first one. In the HPV16 genome, the set LCRS/E7AS spans nucleotides 26–671 (645 bp) and the set pU 1M/2R spans nt 419–656 (237 bp). The set MY 09/11 amplifies the L1-450 segment which spans nt 6584–7035 (451 bp). To explore the feasibility of using duplex PCR to simultaneously amplify E6-250 with primer set pU 1M/2R and L1-450 with primer set MY 09/11, we tested both sets at 0.6 μM concentration. However, our finding that in duplex mixtures the L1-450 band was nearly lost and the E6-250 band was less intense than in uniplex mixtures (Fig. 4) led us to discard the use of duplex mixtures.
Detection of the HPV DNA amplification products
LCR/E7 ~650 bp (E6-650), E6/E7 ~250 bp (E6-250) and L1 ~450 bp (L1-450) products were amplified with primer sets pU 1M/2R, LCRS/E7AS and MY 09/11, respectively (Table 7). Each 50 μl PCR mixture was prepared in 200 μl Axygen tubes (Union City, CA). PCR mixtures contained 2.5 mM MgCl2, 0.4 mM of each dNTP, 0.6 μM of each oligonucleotide, 1.5 units of recombinant Taq DNA polymerase and 25 ng DNA, 20 mM Tris-HCl, 50 mM KCl, pH 8.4 at 25°C; oligonucleotides were synthesized by Accesolab (Mexico) and other components were purchased from Invitrogen (Mexico).
PCR mixtures were incubated 40 cycles in a Techne Touchgene Gradient thermocycler (Staffordshire, England). Before starting they were incubated at 94°C for 4 min. Each cycle consisted of 1 min periods at 94°C (denaturation), 57°C (annealing) and 72°C (extension), followed by 10 min final extension at 72°C. Each run included a positive control with 25 ng of HeLa cell DNA known to contain 10–50 copies of the HPV type 18 genome per cell , and a negative control with water instead of DNA.
PCR products (10 μl from each mixture) were subjected to electrophoresis in 1.5% agarose gels with 5 mM sodium borate, pH 8.5  that were run at 120 V for 90 min.
High-risk HPV typing
To identify HR-HPV types the method of Fujinaga et al.  was used. The size of the two E6-250 restriction fragments generated with endonucleases for each of the seven HR-HPV types that can be identified with this method are shown in Table 2.
Ava I, Ava II and Bgl II endonucleases and their buffers were purchased from Biolabs (Beverly, MA), whereas Rsa I and Acc I and their buffers were obtained from Invitrogen (Carlsbad, CA). Digestion mixtures containing 8 μl of PCR amplification mixture, 1.5 μl endonuclease solution (10 U/μl), 2 μl buffer and 8.5 μl water were incubated 3 h at 37°C. The restriction fragments present in 10 μl from each digestion mixture were analyzed by electrophoresis in 2% agarose gels with 5 mM sodium borate, pH 8.5 (Brody and Kern 2004) that were run at 120 V for 90 min.
Pisani P, Parkin DM, Bray F, Ferlay J: Estimates of the worldwide mortality from 25 cancers in 1990. Int J Cancer. 1999, 83: 18-29. 10.1002/(SICI)1097-0215(19990924)83:1<18::AID-IJC5>3.0.CO;2-M.
Pisani P, Bray F, Parkin DM: Estimates of the world-wide prevalence of cancer for 25 sites in the adult population. Int J Cancer. 2002, 97: 72-81. 10.1002/ijc.1571.
Arrossi S, Sankaranarayanan R, Parkin DM: Incidence and mortality of cervical cancer in Latin America. Salud Publica Mex. 2003, 45 Suppl 3: S306-14.
Reeves WC, Brinton LA, Garcia M, Brenes MM, Herrero R, Gaitan E, Tenorio F, de Britton RC, Rawls WE: Human papillomavirus infection and cervical cancer in Latin America. N Engl J Med. 1989, 320: 1437-1441.
Soe MM: Screening for cervical cancer in developing countries. Bmj. 1992, 304: 983-984.
Lazcano-Ponce EC, Moss S, Alonso de Ruiz P, Salmeron Castro J, Hernandez Avila M: Cervical cancer screening in developing countries: why is it ineffective? The case of Mexico. Arch Med Res. 1999, 30: 240-250. 10.1016/S0188-0128(99)00006-8.
Palacio-Mejia LS, Rangel-Gomez G, Hernandez-Avila M, Lazcano-Ponce E: Cervical cancer, a disease of poverty: mortality differences between urban and rural areas in Mexico. Salud Publica Mex. 2003, 45 Suppl 3: S315-25.
zur Hausen H: Human papillomaviruses and their possible role in squamous cell carcinomas. Curr Top Microbiol Immunol. 1977, 78: 1-30.
Zur Hausen H: Papillomaviruses in human cancers. Mol Carcinog. 1988, 1: 147-150. 10.1002/mc.2940010302.
Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, Snijders PJ, Peto J, Meijer CJ, Munoz N: Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol. 1999, 189: 12-19. 10.1002/(SICI)1096-9896(199909)189:1<12::AID-PATH431>3.0.CO;2-F.
Volkow P, Rubi S, Lizano M, Carrillo A, Vilar-Compte D, Garcia-Carranca A, Sotelo R, Garcia B, Sierra-Madero J, Mohar A: High prevalence of oncogenic human papillomavirus in the genital tract of women with human immunodeficiency virus. Gynecol Oncol. 2001, 82: 27-31. 10.1006/gyno.2001.6244.
Munoz N, Bosch FX, de Sanjose S, Vergara A, del Moral A, Munoz MT, Tafur L, Gili M, Izarzugaza I, Viladiu P: Risk factors for cervical intraepithelial neoplasia grade III/carcinoma in situ in Spain and Colombia. Cancer Epidemiol Biomarkers Prev. 1993, 2: 423-431.
de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H: Classification of papillomaviruses. Virology. 2004, 324: 17-27. 10.1016/j.virol.2004.03.033.
Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, Snijders PJ, Meijer CJ: Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003, 348: 518-527. 10.1056/NEJMoa021641.
Bosch FX, Manos MM, Munoz N, Sherman M, Jansen AM, Peto J, Schiffman MH, Moreno V, Kurman R, Shah KV: 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. 10.1093/jnci/87.11.796.
Berumen J, Ordonez RM, Lazcano E, Salmeron J, Galvan SC, Estrada RA, Yunes E, Garcia-Carranca A, Gonzalez-Lira G, Madrigal-de la Campa A: Asian-American variants of human papillomavirus 16 and risk for cervical cancer: a case-control study. J Natl Cancer Inst. 2001, 93: 1325-1330. 10.1093/jnci/93.17.1325.
Tasa de mortalidad por cáncer cérvico uterino. Sistema para el Seguimiento de la Situación de la Mujer en México (SISESIM). [http://dgcnesyp.inegi.gob.mx/sisesim/sisesim.html]
Lazcano-Ponce E, Herrero R, Munoz N, Cruz A, Shah KV, Alonso P, Hernandez P, Salmeron J, Hernandez M: Epidemiology of HPV infection among Mexican women with normal cervical cytology. Int J Cancer. 2001, 91: 412-420. 10.1002/1097-0215(20010201)91:3<412::AID-IJC1071>3.0.CO;2-M.
Fujinaga Y, Shimada M, Okazawa K, Fukushima M, Kato I, Fujinaga K: Simultaneous detection and typing of genital human papillomavirus DNA using the polymerase chain reaction. J Gen Virol. 1991, 72 ( Pt 5): 1039-1044. 10.1099/0022-1317-72-5-1039.
Torroella-Kouri M, Morsberger S, Carrillo A, Mohar A, Meneses A, Ibarra M, Daniel RW, Ghaffari AM, Solorza G, Shah KV: HPV prevalence among Mexican women with neoplastic and normal cervixes. Gynecol Oncol. 1998, 70: 115-120. 10.1006/gyno.1998.5055.
Munoz N: Human papillomavirus and cancer: the epidemiological evidence. J Clin Virol. 2000, 19: 1-5. 10.1016/S1386-6532(00)00125-6.
Gonzalez-Losa Mdel R, Rosado-Lopez I, Valdez-Gonzalez N, Puerto-Solis M: High prevalence of human papillomavirus type 58 in Mexican colposcopy patients. J Clin Virol. 2004, 29: 202-205.
Berumen J, Unger ER, Casas L, Figueroa P: Amplification of human papillomavirus types 16 and 18 in invasive cervical cancer. Hum Pathol. 1995, 26: 676-681. 10.1016/0046-8177(95)90175-2.
Hernandez-Avila M, Lazcano-Ponce EC, Berumen-Campos J, Cruz-Valdez A, Alonso de Ruiz PP, Gonzalez-Lira G: Human papilloma virus 16-18 infection and cervical cancer in Mexico: a case-control study. Arch Med Res. 1997, 28: 265-271.
Brody JR, Kern SE: Sodium boric acid: a Tris-free, cooler conductive medium for DNA electrophoresis. Biotechniques. 2004, 36: 214-216.
Hwang T: Detection and typing of human papillomavirus DNA by PCR using consensus primers in various cervical lesions of Korean women. J Korean Med Sci. 1999, 14: 593-599.
Noda T, Sasagawa T, Dong Y, Fuse H, Namiki M, Inoue M: Detection of human papillomavirus (HPV) DNA in archival specimens of benign prostatic hyperplasia and prostatic cancer using a highly sensitive nested PCR method. Urol Res. 1998, 26: 165-169. 10.1007/s002400050041.
Bernard HU, Chan SY, Manos MM, Ong CK, Villa LL, Delius H, Peyton CL, Bauer HM, Wheeler CM: Identification and assessment of known and novel human papillomaviruses by polymerase chain reaction amplification, restriction fragment length polymorphisms, nucleotide sequence, and phylogenetic algorithms. J Infect Dis. 1994, 170: 1077-1085.
Wright TC, Denny L, Kuhn L, Pollack A, Lorincz A: HPV DNA testing of self-collected vaginal samples compared with cytologic screening to detect cervical cancer. Jama. 2000, 283: 81-86. 10.1001/jama.283.1.81.
Clifford GM, Rana RK, Franceschi S, Smith JS, Gough G, Pimenta JM: Human papillomavirus genotype distribution in low-grade cervical lesions: comparison by geographic region and with cervical cancer. Cancer Epidemiol Biomarkers Prev. 2005, 14: 1157-1164. 10.1158/1055-9965.EPI-04-0812.
Garland SM: Human papillomavirus update with a particular focus on cervical disease. Pathology. 2002, 34: 213-224. 10.1080/00313020212469.
Molano M, Van den Brule A, Plummer M, Weiderpass E, Posso H, Arslan A, Meijer CJ, Munoz N, Franceschi S: Determinants of clearance of human papillomavirus infections in Colombian women with normal cytology: a population-based, 5-year follow-up study. Am J Epidemiol. 2003, 158: 486-494. 10.1093/aje/kwg171.
Becker TM, Wheeler CM, McGough NS, Parmenter CA, Jordan SW, Stidley CA, McPherson RS, Dorin MH: Sexually transmitted diseases and other risk factors for cervical dysplasia among southwestern Hispanic and non-Hispanic white women. Jama. 1994, 271: 1181-1188. 10.1001/jama.271.15.1181.
Hubbard RA: Human papillomavirus testing methods. Arch Pathol Lab Med. 2003, 127: 940-945.
Londesborough P, Ho L, Terry G, Cuzick J, Wheeler C, Singer A: Human papillomavirus genotype as a predictor of persistence and development of high-grade lesions in women with minor cervical abnormalities. Int J Cancer. 1996, 69: 364-368. 10.1002/(SICI)1097-0215(19961021)69:5<364::AID-IJC2>3.0.CO;2-3.
Richardson H, Kelsall G, Tellier P, Voyer H, Abrahamowicz M, Ferenczy A, Coutlee F, Franco EL: The natural history of type-specific human papillomavirus infections in female university students. Cancer Epidemiol Biomarkers Prev. 2003, 12: 485-490.
Woodman CB, Collins S, Winter H, Bailey A, Ellis J, Prior P, Yates M, Rollason TP, Young LS: Natural history of cervical human papillomavirus infection in young women: a longitudinal cohort study. Lancet. 2001, 357: 1831-1836. 10.1016/S0140-6736(00)04956-4.
Hildesheim A, Wang SS: Host and viral genetics and risk of cervical cancer: a review. Virus Res. 2002, 89: 229-240. 10.1016/S0168-1702(02)00191-0.
Williams OM, Hart KW, Wang EC, Gelder CM: Analysis of CD4(+) T-cell responses to human papillomavirus (HPV) type 11 L1 in healthy adults reveals a high degree of responsiveness and cross-reactivity with other HPV types. J Virol. 2002, 76: 7418-7429. 10.1128/JVI.76.15.7418-7429.2002.
Corona Gutierrez CM, Tinoco A, Navarro T, Contreras ML, Cortes RR, Calzado P, Reyes L, Posternak R, Morosoli G, Verde ML, Rosales R: Therapeutic vaccination with MVA E2 can eliminate precancerous lesions (CIN 1, CIN 2, and CIN 3) associated with infection by oncogenic human papillomavirus. Hum Gene Ther. 2004, 15: 421-431. 10.1089/10430340460745757.
Villa LL, Costa RL, Petta CA, Andrade RP, Ault KA, Giuliano AR, Wheeler CM, Koutsky LA, Malm C, Lehtinen M, Skjeldestad FE, Olsson SE, Steinwall M, Brown DR, Kurman RJ, Ronnett BM, Stoler MH, Ferenczy A, Harper DM, Tamms GM, Yu J, Lupinacci L, Railkar R, Taddeo FJ, Jansen KU, Esser MT, Sings HL, Saah AJ, Barr E: Prophylactic quadrivalent human papillomavirus (types 6, 11, 16, and 18) L1 virus-like particle vaccine in young women: a randomised double-blind placebo-controlled multicentre phase II efficacy trial. Lancet Oncol. 2005, 6: 271-278. 10.1016/S1470-2045(05)70101-7.
Markowitz LE, Dunne EF, Saraiya M, Lawson HW, Chesson H, Unger ER: Quadrivalent Human Papillomavirus Vaccine: Recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR Recomm Rep. 2007, 56: 1-24.
McLemore MR: Gardasil: Introducing the new human papillomavirus vaccine. Clin J Oncol Nurs. 2006, 10: 559-560. 10.1188/06.CJON.559-560.
Harper DM, Franco EL, Wheeler CM, Moscicki AB, Romanowski B, Roteli-Martins CM, Jenkins D, Schuind A, Costa Clemens SA, Dubin G: Sustained efficacy up to 4.5 years of a bivalent L1 virus-like particle vaccine against human papillomavirus types 16 and 18: follow-up from a randomised control trial. Lancet. 2006, 367: 1247-1255. 10.1016/S0140-6736(06)68439-0.
Chardonnet Y, Lizard G, Chignol MC, Schmitt D: Analytical methods for evaluation on whole cells of human papillomavirus infection. Bull Cancer. 1995, 82: 107-113.
Sotlar K, Diemer D, Dethleffs A, Hack Y, Stubner A, Vollmer N, Menton S, Menton M, Dietz K, Wallwiener D, Kandolf R, Bultmann B: Detection and typing of human papillomavirus by e6 nested multiplex PCR. J Clin Microbiol. 2004, 42: 3176-3184. 10.1128/JCM.42.7.3176-3184.2004.
We wish to thank all the women included in this study, as well as Julio Ortiz, Ma. de Lourdes Martínez, Yolanda Terán and Dalila Álvarez for providing the cervical scrapings from San Luis Potosí, and Raúl Rojas and Silvia Quintana for providing those from Guanajuato. We also thank Aurora Londoño for her efficient collaboration in sample processing. The interest and support of the health authorities and of numerous collaborators from San Luis Potosí and Guanajuato, especially Fernando Toranzo and Raúl Martínez of the San Luis Potosí Health Secretariat, and Éctor Jaime Ramírez Barba and Raúl Rojas of the Guanajuato Health Secretariat, is highly appreciated. The research was partially funded by grants FONINV-Gto-C01-5751 and FMSLP-4441. LAMC received a CONACYT scholarship during the course of the work.
The author(s) declare that they have no competing interests.
LAMC performed most of the molecular studies. MSG suggested and set up the PCR-RFLP method. RLR conceived and designed the study, obtained the funds to carry it out and drafted the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
López-Revilla, R., Martínez-Contreras, L.A. & Sánchez-Garza, M. Prevalence of high-risk human papillomavirus types in Mexican women with cervical intraepithelial neoplasia and invasive carcinoma. Infect Agents Cancer 3, 3 (2008). https://doi.org/10.1186/1750-9378-3-3
- Cervical Cancer
- Cervical Intraepithelial Neoplasia
- Multiple Infection
- High Grade Intraepithelial Squamous Lesion
- Cervical Dysplasia