Cervical cancer is the second most common cancer in women worldwide. Every year, 450,000 women are diagnosed with cervical cancer and 220,000 succumb to this disease (
27). Current approaches to cervical cancer control involve lifelong screening using the Papanicolau (Pap) test (
13). The goal of screening is to detect precancerous lesions so that they can be removed prior to the development of cancer. Despite widespread Pap testing, there were an estimated 10,520 new cases of cervical cancer and nearly 4,000 cervical cancer-related deaths in the United States in 2004 (
1). The national health care burden of current screening systems combined with direct costs of treating precancerous and cancerous lesions is in excess of 3.5 billion U.S. dollars per annum (
7).
Infection with human papillomavirus (HPV) is the first and obligate step in the development of cervical cancer (
3,
4). Infection of the cervical epithelium with HPV results in expression of the E6 and E7 proteins, which have been shown to be potent oncogenes. More than 35 different HPV types are capable of infecting the human genital tract (
2,
4,
28). Of these, four types cause the majority of the HPV-related cervical pathology. HPV 16 and 18 together account for 74.6% of all cervical cancers (
23), whereas HPV6 and -11 cause a significant fraction of precancerous lesions which rarely develop into cervical cancer but morphologically are indistinguishable from lesions from more dangerous HPV types (
37). HVP 6 and 11 are responsible for approximately 90% of all genital wart cases (
37).
The HPV LI capsid protein, when expressed recombinantly, assembles into empty viral capsids or “virus like particles” (VLPs) (
12,
15,
16,
29). Several prophylactic vaccines based on HPV LI VLPs are currently in phases II and III clinical development (
14,
17,
36). The VLPs in the vaccine present the immune system with the conformational, neutralizing epitopes found on the natural virus and prime the immune system to generate antibodies that neutralize the virus and prevent infection upon future exposure. Recently, we have shown that a prototype HPV 16 vaccine was 100% efficacious in preventing acquisition of HPV 16 infection and cervical disease among women who were HPV 16 naïve at baseline (
19). These results have led to phase II and III studies of a quadrivalent vaccine to HPV 6, 11, 16, and 18. Early results from a randomized, placebo-controlled, phase II trial have shown this quadrivalent vaccine to be 90% efficacious against HPV 6-, 11-, 16-, and 18-related infection or disease (
36).
The vaccines under development in Merck's HPV vaccine program are designed to be effective when given to subjects who are naïve for at least one of the vaccine HPV types. However, these studies do not include a screening phase. Therefore, it is expected that a small cohort of subjects who are positive to at least one of the vaccine HPV types at day 1 will be enrolled. Thus, in order to fully assess the vaccines clinical impact, the serology assays must be able to reliably distinguish “positive” from “negative” samples, and the serostatus cutoff must be defined. Furthermore, as several phase II clinical studies have shown that vaccine-induced anti-HPV antibody levels are protective against HPV infection and disease, a reliable measure of the duration of the immune response, and hence the duration of efficacy, requires an assay to measure a broad range of antibody titers within diverse clinical populations. These vaccine-induced immune responses must be accurately measured across an extensive range of HPV types without interference from the immune response to other HPV types. Using the multiplexed, competitive Luminex Immunoassay (cLIA) described by Opalka et al. (
26), antibody titers to neutralizing epitopes on HPV 6, 11, 16, and 18 virions can be simultaneously measured in the serum of vaccinated or naturally infected subjects. While the original version of the assay generated accurate results, we modified and formally validated the assay to ensure it was sensitive, specific, and robust enough to meet the needs of testing samples from our large epidemiology studies and phase II and III vaccine clinical trials. To this end, we evaluated the analytical sensitivity, specificity, accuracy, precision, and reproducibility of the assay. Lastly, we determined the clinical sensitivity and specificity of the assay and set serostatus cutoffs for HPV 6, 11, 16, and 18 using serum samples from individuals of diverse geographical regions and age groups with low or high probability of having HPV antibodies.
MATERIALS AND METHODS
Virus-like particles (VLPs).
The L1 genes of HPV types 6, 11, 16, and 18 were expressed in the yeast
Saccharomyces cerevisiae. The products were then purified from lysates of the yeast as previously described (
15,
16,
24,
29) with modifications (
11). The VLPs used in the development of the assay were from consistency lots used in phase III vaccine clinical trials.
Covalent conjugation of HPV VLPs to Luminex microspheres.
The HPV VLPs for types 6, 11, 16, and 18 were each covalently conjugated to the free carboxyl groups on xMAP Multi-Analyte COOH Microspheres (Luminex Corporation, Austin, TX) 6, 11, 16, and 18, respectively, using an
N-hydroxysulfosuccinimide (sulfo-NHS)-enhanced carbodiimide-mediated conjugation chemistry (
31) as previously described (
26). Carboxylated microspheres were brought to room temperature (RT) and sonicated followed by vortexing for approximately 1 min to obtain a homogeneous distribution of microspheres. An aliquot of 1.25 × 10
7 microspheres from the stock vial was added to a clear 1.5-ml copolymer microcentrifuge tube (USA Scientific, Ocala, FL) and pelleted by centrifugation. Microspheres were then washed with 500 μl of 0.1 M NaH
2PO
4 (phosphate buffer) pH 6.2 resuspended by sonication and then repelleted. To activate the carboxyl sites on the microspheres, 50 μl of a 50-mg/ml solution of
N-hydroxysulfosuccinimide (Pierce, Rockford, IL) and 50 μl of a 50-mg/ml solution of 1-ethyl-3-(3-dimethylaminopropyl)-carbodiimide hydrochloride (EDC; Pierce, Rockford, IL) solution were added to the microspheres and 400 μl of phosphate-buffered saline (PBS). The tubes were sonicated for 1 min and then incubated for 20 min at RT in the dark. Following activation, the microspheres were washed with 500 μl of either 50 mM 2-(
N-morpholino)ethane sulfonic acid buffer, pH 6.0 (MES; Sigma, St. Louis, MO) or PBS, pH 7.4, to remove residual EDC and sulfo-NHS.
The HPV VLPs for types 6, 11, 16, and 18 were diluted in PBS or MES coupling buffer at concentrations ranging from 12.5 to 300 μg/ml. Upon dilution, 500 μl of each VLP solution was added to the corresponding microsphere (i.e., VLP 6, 11, 16, and 18 to microsphere 6, 11, 16, and 18, respectively). Each tube of microspheres was then sonicated and vortexed for approximately 1 min, covered in foil, and placed on a rotator overnight at RT.
Following the overnight incubation, VLP-microspheres were washed twice with 1 ml of wash buffer, in either PBS, pH 7.4, containing 0.05% Tween-20 or PBS, pH 7.4, containing 1% Triton X-100 (TX100). After washing, VLP-microspheres were resuspended in 1.0 ml of histidine buffer (20 mM histidine, 0.5 M NaCl, pH 6.2) containing 1% bovine serum albumin (histidine-BSA) to block any remaining activated carboxyl sites on the microspheres. VLP-microspheres were diluted to 1.0 × 106 beads/ml in histidine-BSA buffer and stored at 4°C in the dark.
Monoclonal antibodies.
The H6.B10.5 (
11), H6.M48 (
11), H11.B2, monoclonal antibody (MAb) 8740 (Chemicon, Temecula, CA) (
10), K11.B2, H16.V5 (
11), and H18.J4 (
11) were used to detect neutralizing epitopes on VLPs 6, 11, 16, and 18. The monoclonal antibodies were directly conjugated to R-phycoerythrin (PE) via succinyl esters containing either a thiol or maleimide reactive group (
21) by Chromaprobe Inc. (Maryland Heights, MO). The use of these heterobifunctional cross-linkers results in a nonreducible, thioether linkage.
HPV competitive Luminex Immunoassay (cLIA).
The basic procedures of the assay have previously been described by Opalka et al. (
26). The MAbs for HPV types 6, 11, 16, and 18 were pooled at varying concentrations (0.25 to 40 μg/ml) and diluted in PBS with 1% BSA and 0.05% sodium azide (PBS-BN) containing varying concentrations of TX100 ranging from 0% to 20%. The final concentrations of the MAbs in each well ranged from 0.0625 to 10 μg/ml. The VLP-microspheres for types 6, 11, 16, and 18 were pooled in PBS-BN or PBS-BN plus 1% TX100 to a concentration of 2.0 × 10
5 VLP-microsphere per ml per VLP type.
A quadriplex reference standard solution was prepared at a concentration of 1,000 milli-Merck units per milliliter (mMU/ml) for each HPV type in a sample matrix. Various sample matrices were examined including normal human serum (NHS), normal goat serum (NGS), PBS, PBS with 1% BSA, PBS with 1% recombinant human albumin (rHA) (Delta Biotech, Nottingham, United Kingdom), defibrinated plasma (DFP) or fibrin-free human plasma (Serological Corporation, Norcross, GA), and 3% (wt/vol) nonfat dried milk (Sigma, St. Louis, MO) in PBS and antibody-depleted human serum (ADHS) (Human Serum Special Stripped; Valley Biomedical, Winchester, VA). The four reference standards were serum samples from African green monkeys (AGM) hyperimmunized with HPV L1 VLP type 6, 11, 16, or 18. Each standard was heat inactivated prior to use for 30 ± 2 min at 56 ± 1°C. The quadriplex reference standard was diluted in a twofold serial dilution in the sample matrix to create a 12-point standard curve in duplicate with final well concentrations ranging from 0.25 to 500 mMU/ml. Fifty microliters of the reference standard was added to 25 μl of the MAb-PE cocktail and 25 μl of the VLP-microspheres.
The assay plate setup contained four controls, each in duplicate. Three of the controls used in the assay were ADHS spiked with specific amounts of each heat-inactivated reference standard. The “high” control contained HPV6, -11, -16, and -18 standards at concentrations of 120, 140, 140, and 140 mMU/ml, respectively. The “low” control contained HPV6, -11, -16, and -18 standards at concentrations of 34, 40, 40, and 40 mMU/ml, respectively. The “negative” control contained HPV 6, 11, 16, and 18 standards at concentrations of 17, 20, 20, and 20 mMU/ml, respectively, and the fourth control was ADHS alone.
Prior to testing, samples were heat inactivated for 30 ± 2 min at 56 ± 1°C. After heat inactivation, 25 μl of control or sample were added to 25 μl of ADHS. All standards, controls, and samples were tested in duplicate. To each well on the assay plate was added 25 μl of the MAb-PE quadriplex solution, followed by 25 μl of the quadriplex VLP-microsphere solution. Final assay conditions consisted of 25 μl of serum, 25 μl of ADHS, 25 μl of MAb-PE quadriplex (0.1 μg/ml for each MAb), and 25 μl of VLP-microspheres (5,000 VLP-microspheres per well per type). The final concentration of TX100 present in the MAb and VLP-microsphere assay buffers was 0.3%. After final addition of the VLP-microspheres, the plate was covered with a foil seal and placed on a shaker (600 to 800 rpm) at RT for 2 to 25 h. The contents of the assay plate were then transferred to a 1.2 μm Low Protein Binding Filter Plate (Millipore, Billerica, MA) prewet with 200 μl of assay wash buffer consisting of PBS with or without 1% TX100. The plates were washed three times with 200 μl of assay wash buffer and resuspended in 125 μl of assay wash buffer for analysis on a Bio-Plex Suspension Array System (Bio-Rad, Hercules, CA).
Differences between the proof-of-concept assay described by Opalka et al. (
26) and the current validated assay being used to support HPV vaccine clinical trials are shown in Table
1.
Assay validation.
Validation of the cLIA for HPV 6, 11, 16, and 18 was performed over 6 days during a 2-week period. Three analysts, using one specific Bio-Plex instrument each, generated two assay plates per day either manually or on a “TECAN Genesis Workstation”. Two MAb-PE lots and three separate preparations of VLP-microspheres were examined. From the standard curves generated during the validation, the limits of detection, and the quantifiable ranges of the HPV cLIA were determined for each of the HPV types. Additionally, specification limits for the slope and 50% effective concentration (EC50) parameters were determined to qualify future assay runs.
Nine human serum sample panels were tested in the validation of the assay (Table
2). Panels were divided into four groups. Group 1 (panels A, B, and C) contained sera likely to be true negatives. These three panels contained samples from male adolescents, female adolescents, and young adult women (16 to 23 years of age) who reported having two or fewer lifetime sexual partners and who repeatedly tested negative in HPV6, -11, -16, and -18 PCR assays. Group 2 (panels D, E, and F) consisted of sera likely to be true positives. These samples were obtained either from individuals that were enrolled in clinical trials and had received an HPV6, -11, -16, and -18 vaccine or from nonvaccinated subjects with four or more lifetime sexual partners. Panel F was also used to evaluate the dilutability of the assay. Group 3 (panel G) consisted of sera randomly selected from American women who had no history of HPV vaccination. Panel G was used to evaluate the performance of the redefined serostatus cutoff and to provide an estimate of the baseline percent positivity for HPV types 6, 11, 16, and 18 in the primary clinical population. Lastly, group 4 (panels H and I) was used to assess assay robustness and precision between analysts, MAb-PE lots, and VLP-microsphere lots. Panel H (samples with low to medium titers) was tested twice per day per analyst over the 2-week period. Panel I (negative, low, and high spiked ADHS-spiked controls) was tested on each of the 42 validation plates.
Data analysis.
The serum antibody titer was obtained by measuring the ability of a test sample to inhibit MAb-PE binding to the VLPs and comparing that value to that of a dilution series of the reference standard. The reference standard dilution series was modeled using the four-parameter logistic function (
25), and the test sample concentration was calibrated from the fitted standard curve. The reference sera used for the standard curve were assigned arbitrary antibody concentrations expressed in milli-Merck units per milliliter. Data, in median fluorescence intensity (MFI) units, were processed using either the Bio-Plex Manager 3.0 software (Bio-Rad, Hercules, CA) or using a Microsoft Excel spreadsheet (Redmond, WA).
Since samples were duplicated within a run, intra-assay precision was estimated using the variability estimates between the two or three replicates of each control sample within each run \(({\hat{{\sigma}}}_{S}^{2})\). To determine inter-assay variability, estimates of analyst to analyst variability \(({\hat{{\sigma}}}_{A}^{2})\), run within analyst variability \([{\hat{{\sigma}}}_{R(A)}^{2}]\), sample by analyst variability \(({\hat{{\sigma}}}_{SxA}^{2})\), and sample by run within analyst variability \([{\hat{{\sigma}}}_{SxR(A)}^{2}]\)were combined. Variability estimates were obtained on the natural log transformed titers using the MIXED procedure in SAS (SAS Institute Inc., Cary, NC). An estimate of the total assay precision (% relative standard deviation [RSD]) was calculated as 100 × \([e\sqrt{{\hat{{\sigma}}}_{S}^{2}\ {+}\ {\hat{{\sigma}}}_{A}^{2}\ {+}\ {\hat{{\sigma}}}_{R(A)}^{2}\ {+}\ {\hat{{\sigma}}}_{SxA}^{2}\ {+}\ {\hat{{\sigma}}}_{SxA}^{2}}\ {-}1]\).
DISCUSSION
An immunoassay that measures HPV type-specific antibodies simultaneously is preferred to running multiple separate tests. Thus, a competitive immunoassay that simultaneously measures type-specific antibodies to neutralizing epitopes on HPV 6, 11, 16, and 18 was developed. We have modified an HPV Luminex immunoassay to be more sensitive and more robust for testing serum samples from individuals enrolled in epidemiology studies and HPV vaccine clinical trials. The increase in sensitivity was necessary to measure the low titers of anti-HPV antibodies in the serum from naturally infected individuals, a requirement for epidemiological and vaccine studies. By analyzing each step of the assay, we have characterized the effects of reagents, timing, and operator variability on the assay. The planned long-term use of the assay required the development of a reference standard to monitor the immune response and to normalize results across runs. In addition, we have developed negative-, low-, medium-, and high-titer quality control panels to monitor the consistency of our results periodically. The assay was formatted for automation so that it could be run on a high-throughput TECAN Genesis Workstation. Following optimization, validation of the assay was performed to (i) assess the reproducibility and precision of the assay, (ii) establish a quantifiable range, (iii) assess the effect of dilution, (iv) establish a serostatus cutoff, and (v) evaluate the performance of the serostatus cutoff.
To perform the competitive multiplexed assay, MAbs specific to a single HPV VLP type were needed. To improve the analytical specificity of the assay we switched from the H6.B10.5 MAb to the H6.M48 MAb and from the H11.B2 MAb to the K11.B2 MAb (Fig.
1B). Another advantage of using H6.M48 over H6.B10.5 is that H6.B10.5 was a difficult antibody to purify and to conjugate to PE. Because the HPV cLIA is used to measure the immune response to four separately manufactured VLPs, it is important that the MAbs recognize type-specific, neutralizing epitopes on each of the VLPs. Ideally, these MAbs would also recognize the immunodominant epitopes that are elicited following natural infection, like has been shown for H16.V5 (
39).
To reduce the number of steps in the assay and reduce sources of error, we chose to directly conjugate R-phycoerythrin to the four MAbs. The concentration of MAb and VLP was also optimized. As expected for a competitive immunoassay, optimization experiments showed that as the MAb concentration (Fig.
2A) or VLP concentration (Fig.
2B) decreased, the sensitivity of the assay increased. However, if the MAb or VLP concentration was too low, the assay became less reproducible. The final concentrations established were 0.1 μg/ml for each of the MAb-PEs and 25, 25, 100, and 25 μg/ml of VLPs for HPV 6, 11, 16, and 18, respectively. These optimized concentrations provided acceptable sensitivity and a broad dynamic range and were robust in routine operation.
Several different sample matrices were evaluated to improve the clinical sensitivity of the assay. Previously, normal human serum was used to limit any differences between the sample and the sample matrix, but identifying large quantities of normal human sera that were negative for HPV 6, 11, 16, and 18 and consistent from donor to donor was extremely difficult. Experiments comparing different sample matrices used to dilute the samples and reference sera were performed and identified ADHS as the best alternative matrix (Fig.
3). Seven other matrices were also evaluated, including PBS, PBS with 1% BSA, PBS with 1% rHA, NGS, DFP, and 3% nonfat milk in PBS, but none demonstrated the reproducibility and similarity to normal human sera as ADHS. In addition, ADHS was similar in composition to human sera and easy to obtain and was consistent from lot to lot. The PBS-based matrices were not used due to significant deviations from the standard curves generated in NHS. Because of the erratic shapes of the standard curves generated in PBS-based matrices and other matrices, such as NGS, ADHS was chosen as the sample diluent to replace NHS.
Several different preanalytical factors can adversely affect results in a clinical assay. These include sample collection and processing variables such as anticoagulants in the collection container, storage and transport temperatures, lipemia, and hemolysis. An important finding during the development of the assay was the observation that a heat-labile factor, present in many individuals' sera, could interfere with binding of the detection MAb, resulting in higher than expected antibody titers. Heat inactivation of freshly thawed samples significantly decreased the number of positive samples (Fig.
4B) to be more in line with expectations of 5 to 20%, as supported by natural history data (
6,
22,
32). The heat-labile factor(s) that caused the false positive results has not yet been identified. Heat inactivation did not impact the titers of samples that had been stored at 4°C for several weeks and did not affect antibodies in the reference standard in comparison to a non-heat-inactivated reference standard (Fig.
4C). These findings led us to include a heat inactivation step of 56 ± 1°C for 30 ± 2 min as part of our sample preparation procedures.
Following assay optimization, a formal validation of the assay for clinical testing was performed. To examine the assay in the target population of females ranging in age from 9 to 45 years and to define the accuracy and precision of the assay, we created nine human serum panels (Table
2). Panels H and I were used to determine the robustness and precision of the assay and also to set specification limits on the controls that are tested in routine operation. Panels A, B, and C contained samples that were most likely to be HPV negative. Panels A and C contained samples from HPV PCR-negative women with two or fewer lifetime sexual partners. Panel B was comprised of samples from both male and female adolescents ages 10 to 12 and was used to determine whether there were differences between adults and adolescents. Titer results from these samples were used to determine the background titers of HPV-negative individuals and to set the serostatus cutoffs. The serostatus cutoffs were established by identifying the lowest or next to lowest dilution titer that all of the samples in panels A, B, and C yielded negative results. Using this methodology, the serostatus cutoffs were set at 20, 16, 20, and 24 mMU/ml for HPV6, -11, -16, and -18, respectively. Of note, one sample in panel B yielded a titer of 31 mMU/ml for HPV16 which was greater than the serostatus cutoff of 20 mMU/ml.
Panels D, E, and F included samples from individuals likely to be positive for HPV antibodies because of prior HPV vaccination or sexual history. These samples were used to assess the assay in the positive region of the standard curve. Panel F (serum from vaccinated individuals) was also used to assess the dilutability of the assay. Results showed that the dilutability of the HPV cLIA ranged from 4.2% to 13.6% per twofold increase in dilution for samples diluted from 1:4 to 1:512 (Table
3). Since the upper limit of quantitation for all four types is 400 mMU/ml, dilutions are necessary to calculate titers greater than 400 mMU/ml and have been set at 1:40, 1:400, and 1:4,000. To be conservative, the lowest dilution titer that fell within the quantifiable range was reported.
Panel G consisted of randomly selected baseline sera from subjects enrolled in Merck HPV vaccine clinical trials and were used to assess the newly defined serostatus cutoff. The percentage of positive samples in panel G was 10.1%, 3.0%, 12.9%, and 3.0% for HPV 6, 11, 16, and 18 which was in line with our expectation of less than 20% for all four HPV types. Previous serologic studies have demonstrated that 20 to 50% of women with HPV DNA do not have detectable type-specific anti-HPV antibodies (
6). In our preliminary analysis, about 50% of the women that had a PCR-positive result also had an antibody-positive result.
It is worth noting that the clinical sensitivity and specificity of the validation was performed using predominantly samples from women. Other studies have shown that serology assays may underestimate past or current HPV infection in males (
34) and have reported HPV seroprevalence to be generally higher in women than in men (
20,
33).
The challenges of developing an HPV serology assay to monitor both natural infection and response to vaccination have also been described by others (
18,
38). The most formidable challenge for developing a serology assay to measure past or current infection is that not all people seroconvert following exposure and that no reference panels exist to determine the clinical sensitivity and specificity of an assay. Furthermore, women with detectable viral DNA in genital epithelial cells may not have detectable HPV antibodies, as the antibody response may occur many months after initial infection (
5,
35). As the majority of HPV infections of the anogenital tract resolve over 3 to 5 years (
30), serologic testing is a better measure of cumulative exposure rather than ongoing infection. Data from ongoing epidemiology studies and clinical trials will provide an opportunity to study the natural history of antibody responses in relationship to HPV DNA detection and clinical outcome, to measure vaccine-induced immune responses over time, and to identify clinically relevant immunologic correlates of protection.
In summary, the HPV 6, 11, 16, and 18 cLIA has proven to be a robust, sensitive, and high-throughput assay and is currently being used for epidemiology studies and evaluating the efficacy of prophylactic HPV vaccines.