Journal of Clinical Microbiology, February 2009, p. 287-293, Vol. 47, No. 2
0095-1137/09/$08.00+0 doi:10.1128/JCM.01531-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Departments of Microbiology,1 Obstetrics and Gynaecology,2 Anatomical and Cellular Pathology,3 Community and Family Medicine, Faculty of Medicine, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, New Territories, Hong Kong Special Administrative Region, People's Republic of China4
Received 8 August 2008/ Returned for modification 10 October 2008/ Accepted 15 November 2008
|
|
|---|
|
|
|---|
Women who have tested positive for high-risk HPV tend to harbor abnormal cervical cytology, and some studies have reported a correlation between viral load and disease severity (19, 20, 24, 28, 31). However, the consistency of these findings has been questioned, and it is now realized that the relationship is a lot more complex (5-7, 35). As the severity of disease increases, HPV is often found integrated into the host genome, making the interpretation of viral load even more complex. Furthermore, the available viral load data have been generated using different methodologies that may not be directly comparable (13, 17, 28, 29).
The HPV genome can exist in two physical forms: a closed circular episomal form or linearized and integrated into the host genome. Integration of the viral genome is often associated with the disruption of the E2 gene, which has regulatory control of the oncogenes E6 and E7. Nevertheless, this concept has been questioned by a recent study where the integration of the HPV16 genome was not necessarily associated with the overexpression of E6 and E7 (14). HPV18 has been found to be strongly associated with adenocarcinoma of the cervix. Studies carried out with HPV18 have shown that the viral genome is almost invariably found to be integrated in high-grade cervical intraepithelial lesions and invasive cancers (2, 9, 22, 27, 32), and these findings potentially allow the use of integration status as a marker to triage HPV18 infections with different clinical implications. However, such applications should be supported by data generated from a large series of clinical samples. In this study, we have used refined methods to determine the viral load, the physical status of the viral genome, and the type of HPV18 infection (whether single infection or coinfection) in a cohort of Chinese women with different grades of cervical lesions in an attempt to provide some information about the usefulness of these viral parameters in clinical setting.
|
|
|---|
HPV detection and genotyping. Total DNA was extracted from cervical scrapes using a commercially available extraction kit (QIAamp DNA Mini Kit; Qiagen GmbH, Hilden, Germany). HPV genotypes were determined using a Linear Array HPV Genotyping Test (Roche, Molecular Systems, Inc., CA) which can detect 37 HPV types. Samples found to be positive for HPV18 were subjected to viral load and integration status determination.
Viral load, E2 physical status, and integration. Real-time quantitative PCR (qPCR) targeting the HPV18 E7 gene was used to determine the total crude (episomal and integrated) viral load. This region was selected as it is retained in both the episomal and integrated forms. To account for variation in the number of cells collected in each cervical scrape sample, levels of the housekeeping gene, beta-actin, were also determined using qPCR. Crude viral load measurements were normalized with the beta-actin gene levels using the following formula: (E7 gene copy number/beta-actin gene copy number) x 2. Values are expressed as the number of viral copies/cell equivalent.
Based on the assumption that HPV integration disrupts the E2 gene, the integrity of E2 was taken as a marker that HPV had integrated into the host genome. To account for the possibility that disruption might have involved only part of the E2 gene as has been shown in other HPV types (5-7), four nonoverlapping PCRs were designed to span the entire E2 gene, covering the amino (N)-terminal, hinge, and carboxyl (C)-terminal regions. Primer sets used for qPCR are listed in Table 1.
|
View this table: [in a new window] |
TABLE 1. Primer sets for real-time PCRs
|
In total, six sets of real-time PCRs were performed for each sample using a Power SYBR Green PCR master mix (Applied Biosystems, Foster City, CA) with a 7900 HT Sequence Detection System (Applied Biosystems, Foster City, CA). A standard curve for E7 and for each of the E2 gene targets was generated by plotting threshold cycle values against 10-fold serial dilutions of plasmids containing the full genome of HPV18 (ATCC, Rockville, MD). The standards used for quantitation of the beta-actin gene were purchased commercially (Applied Biosystems, Foster City, CA).
A 200-µl aliquot of cervical sample was used for DNA extraction. Five microliters of the extracted preparation was amplified in a 25-µl reaction mixture containing 0.125 to 0.25 µmol of primers (Table 1). The cycling conditions were 95°C for 10 min, followed by 40 cycles at 95°C for 15 s, 58°C for 15 s (except for beta-actin, which required 56°C for 15 s), and 72°C for 30 s. All real-time PCR assays showed a wide linear range that covered at least 10 to 10,000,000 copies/µl, and all showed high amplification efficiency. The specificity of each amplification was confirmed by checking the dissociation curve against the expected melting temperature of the amplification product.
Statistical analysis. The differences in viral load levels between groups were compared by the Mann-Whitney U test or the Kruskal-Wallis test using SPSS software (version 15.0; SPSS Inc., IL). The distribution of categorical variables including integration status and location of disruption among groups was assessed by the chi-square test or the Fisher's test, as appropriate. This statistical analysis was performed using the Statcalc program (Epi Info; Centres for Disease Control and Prevention, Atlanta, GA). Two-sided P values less than or equal to 0.05 were regarded as significant.
|
|
|---|
Viral load. The crude E7 gene copy numbers were distributed over a wide range from 4.4 to 21,358,178.0 (median, 6,622.4; interquartile range, 245.0 to 185,546.7) gene copies/µl. The gene copy numbers for beta-actin ranged from 0.05 to 44,125.1 (median, 1,237.2; interquartile range, 106.3 to 8,123.4) gene copies/µl. The viral load after normalization with beta-actin ranged from <0.01 to 37,850.3 (median, 12.8; interquartile range, 3.3 to 96.2) gene copies/cell equivalent. The distributions of crude and normalized viral load and beta-actin gene levels among the different grades of cervical lesion are shown in Fig. 1. The crude total viral load (episomal and integrated) appeared to increase with disease severity. However, when normalized with the beta-actin gene levels, this association was no longer observed (P = 0.830 by a Kruskal-Wallis test). Study subjects were arbitrarily divided into three age groups (20 to 34 years, 35 to 54 years, and 55 to 76 years) to analyze for a possible effect between age and viral load. No significant association between viral load and age was observed (P = 0.591 by a Kruskal-Wallis test). The lack of association was further confirmed by dividing subjects into two age groups, i.e., 20 to 40 years and 41 to 76 years (P = 0.402 by Mann-Whitney U test).
![]() View larger version (16K): [in a new window] |
FIG. 1. Distribution of viral load among different grades of cervical lesions. The middle line indicates the median; the box represents the interquartile range. Lines extending from each box are the upper and lower limits. Viral load was determined as described in Materials and Methods.
|
![]() View larger version (15K): [in a new window] |
FIG. 2. Distribution of viral load according to the physical status of the genome (a) and type of HPV18 infection (single or coinfection) (b). In both panels, the middle line indicates median, the box represents interquartile range, and lines extending from each box are the upper and lower limits. Viral load was determined as described in Materials and Methods.
|
E2 disruption status. Altogether, 31/78 (39.8%) samples were found to have E2 disruptions. The most common disruption occurred at the N1 (nucleotides 2821 to 3043) region in 22/31 (71.0%) samples, followed by the hinge (67.7%), the C-terminal (64.5%), and then the N2 (nucleotides 3109 to 3265) (61.3%) regions. The preferential site of disruption did not change with the severity of lesion (Fig. 3). The proportion of samples with E2 disruption (any site combined) did not show a significant association with the degree of cervical lesion (E2 disruption proportion in normal/LGSIL/CIN 1 cases, 41.7%; in CIN 2/CIN 3 cases, 33.3%; and in ICC cases, 42.9%; PTrend = 0.737, by the chi-square test for trend) (Fig. 3).
![]() View larger version (32K): [in a new window] |
FIG. 3. Proportion and location of E2 disruption according to the type of cervical lesions. Specimens with disruption involving more than one E2 region are counted more than once. H, hinge region.
|
Infection form. Of the 78 subjects, 46 (59.0%) had single infections of HPV18 alone. The other 32 (41.0%) subjects had coinfections of HPV18 with other HPV types (Table 2). The most common coinfecting HPV type was HPV16 (50.0%), followed by HPV52 (15.6%); HPV types 31, 59, 61, and 81 (6.3% each); and HPV types 51, 53, and 84 (3.1% each). Figure 2b shows the distribution of viral load among samples with HPV18 single infections and coinfections. The viral loads obtained for HPV18 single infections were significantly higher than the load in HPV18 coinfections (median [interquartile range] = 21.8 [6.75 to 104.8] versus 4.5 [0.05 to 72.25] viral copies/cell equivalent; P = 0.046 by the Mann-Whitney U test). When the form of the infection (single or coinfection) was analyzed with respect to the severity of cervical lesion, a significant trend of increase in the proportion of single infections was found (normal/LGSIL/CIN 1, 33.3%; CIN 2/3, 58.3%; ICC, 66.7%; PTrend = 0.050 by the chi-square test for trend).
|
View this table: [in a new window] |
TABLE 2. Viral load and viral genome physical status of HPV18 infection according to the degree of cervical lesion
|
Adeno/adenosquamous carcinoma versus squamous cell carcinoma. When the viral load results were analyzed with respect to the histological group of adeno/adenosquamous carcinoma (n = 21) and squamous cell carcinoma (n = 21), no significant difference was found (P = 0.400 by the Mann-Whitney U test) though a wider range of viral load was seen for the squamous cell carcinoma (median, 8.5 copies/cell equivalent; interquartile range, 0.05 to 53.4 copies/cell equivalent) than with adeno/adenosquamous carcinoma (median, 16.3 copies/cell equivalent; interquartile range, 7.0 to 36.4 copies/cell equivalent). The physical status of the viral genome also did not differ between the two histological groups; the pure episomal form was found in 10 (47.6%) cases of squamous cell carcinoma and in 14 (66.7%) cases of adeno/adenosquamous carcinoma (P = 0.212 by the chi-square test).
When the infection form (single or coinfection) of the two histological groups was analyzed, it was found that the proportion of HPV18 single infections was significantly higher for adeno/adenosquamous carcinoma cases (85.7% versus 14.3%; P = 0.009 by the chi-square test).
|
|
|---|
In this study, we found that normalization with a human gene to adjust for variation in the degree of cellularity of the specimen is important for an accurate interpretation of the association between viral load and disease severity. Our results showed that as the disease worsened, the crude viral load appeared to increase. However, when the viral load was normalized, such an association was no longer seen. The observation of a higher crude viral load in high-grade lesions may be explained by the fact that abnormal cells express fewer intracellular adhesion molecules than the normal cells and thus are more readily exfoliated (28). Although we were able to normalize the viral load to minimize biases due to variation in the cellular content among different specimens, it was not possible to distinguish the HPV-infected and noninfected portion of cells collected in the cervical scrape sample. While cervical scrape is the most convenient and feasible sample type, it is inevitable that the sample will contain a variable mixture of relevant (HPV infected) and irrelevant (noninfected) cells.
We found that the interpretation of viral load can be complicated by its association with the physical state of the viral genome and with the form of HPV infection (single or coinfection). Though not reaching statistical significance, specimens harboring the episomal viral genome form seemed to have a higher viral load. This observation is in line with the fact that viral replication can occur only in the episomal form. While coinfection with multiple HPV types is common, the biological interaction between HPV types in such circumstances is unknown. Our finding of a significantly lower viral load in coinfections is in line with the hypothesis that the coinfecting types are competing with each other for the limited microelements required for viral replication. Overall, our results indicated that measuring viral load to assess HPV18 infection is premature for clinical use at the present stage.
Many previous studies have found that almost all HPV18 infections exist in an integrated form in women with invasive cancer (2, 9, 15, 22, 32). In the current study, we found a rather high proportion of ICC (57.1%) harboring only the episomal form of the virus. We believe one explanation for this discrepancy could be the sensitivity of the method used. Methods such as Southern blotting or two-dimensional gel electrophoresis were often employed to document the physical state of the HPV genome in previous studies (2, 9, 10, 22, 27). One disadvantage of these methods is that they often require a large amount of DNA to produce a positive signal. In contrast, the method employed in the current study, real-time PCR, is able to detect trace amounts of DNA. A recent study by Huang et al. (17) also reported that the episomal form of HPV18 was seen in all lesion grades ranging from CIN 1 to cervical cancer, with the highest proportion of the episomal form seen in cervical cancer stages II to IV (66.9%). We also obtained a similar proportion of the episomal form (57.1%) for ICC cases. A well-established notion is that disruption of the E2 gene results in the overexpression of E6 and E7, which is the principal oncogenic mechanism conferred by high-risk HPV. Our findings together with those from Huang et al. (17), suggest that there might be alternative oncogenic pathways that bypass E2 disruption.
To date, only a few studies have examined the integration state of the HPV18 genome in normal or low-grade lesions (17). We found that integration was not uncommon (41.7%) in the normal/LGSIL/CIN 1 group. This indicates that integration occurs early in the disease progression process, and probably integration alone cannot drive the full progress to high-grade or invasive cancer.
We observed that the most common site of disruption, regardless of lesion grade, was the N1 region. The available information on the location of disruption is mainly derived from HPV16, where the hinge region was identified as the most common site of disruption (1, 21). However, if the hinge region had been selected for our study, we could potentially have missed 32.3% of the E2 disruptions. These data emphasize the importance of selecting the right target of the E2 gene as a surrogate marker for determining the state of genome integration.
One should note that integration of the HPV genome could also lead to disruption of gene regions other than the E2. The current study might have missed integrations that disrupt only viral gene regions outside the E2. Nevertheless, the oncogenic consequence of such non-E2-disrupting integrations is unknown, and there is no evidence to suggest that these integrations could lead to an overexpression of E6 and E7. Thus, the approach of this study was to perform a more comprehensive analysis on four nonoverlapping regions of the most important gene, E2.
Adenocarcinoma of the cervix has a strong association with HPV18 (35). Here, we have made some subgroup analyses according to the histological type. Although no significant differences in viral load or the physical state of the genome between adeno/adenosquamous carcinomas and squamous cell carcinoma were seen, we found that the proportion of single infections was significantly higher among adeno/adenosquamous carcinomas than in squamous cell carcinomas. Further studies are required to explore the implication of this finding.
Our study has provided further information on the complex relationship between viral load, integration of the viral genome, and disease severity. These data, which can be generated by a widely available technique, real-time PCR, on a routinely feasible specimen, cervical scrape, should inform researchers about the potential advantages and drawbacks of applying viral load and integration testing in a clinical setting.
Published ahead of print on 26 November 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»