Comprehensive Analysis and Selection of Anthrax Vaccine Adsorbed Immune Correlates of Protection in Rhesus Macaques

ABSTRACT Humoral and cell-mediated immune correlates of protection (COP) for inhalation anthrax in a rhesus macaque (Macaca mulatta) model were determined. The immunological and survival data were from 114 vaccinated and 23 control animals exposed to Bacillus anthracis spores at 12, 30, or 52 months after the first vaccination. The vaccinated animals received a 3-dose intramuscular priming series (3-i.m.) of anthrax vaccine adsorbed (AVA) (BioThrax) at 0, 1, and 6 months. The immune responses were modulated by administering a range of vaccine dilutions. Together with the vaccine dilution dose and interval between the first vaccination and challenge, each of 80 immune response variables to anthrax toxin protective antigen (PA) at every available study time point was analyzed as a potential COP by logistic regression penalized by least absolute shrinkage and selection operator (LASSO) or elastic net. The anti-PA IgG level at the last available time point before challenge (last) and lymphocyte stimulation index (SI) at months 2 and 6 were identified consistently as a COP. Anti-PA IgG levels and lethal toxin neutralization activity (TNA) at months 6 and 7 (peak) and the frequency of gamma interferon (IFN-γ)-secreting cells at month 6 also had statistically significant positive correlations with survival. The ratio of interleukin 4 (IL-4) mRNA to IFN-γ mRNA at month 6 also had a statistically significant negative correlation with survival. TNA had lower accuracy as a COP than did anti-PA IgG response. Following the 3-i.m. priming with AVA, the anti-PA IgG responses at the time of exposure or at month 7 were practicable and accurate metrics for correlating vaccine-induced immunity with protection against inhalation anthrax.

T o date, there has not been a systematic evaluation of the relationship between anthrax vaccine-stimulated humoral and cell-mediated immune responses, their relative contributions to protection, or their comparative importance when used singly or in combination to predict the probability of survival in animal models or in humans.
Anthrax toxin protective antigen (PA) is the primary immunogen in licensed anthrax vaccines in the United States and the European Union, as well as in many of the second-generation anthrax vaccines in current development (1). Consequently, the quantitative analysis of anti-PA IgG antibody levels and lethal toxin neutralization activity (TNA) in serum are generally accepted as immunological correlates of protection (COP) for vaccine efficacy in animal models (2). Anti-PA IgG levels and TNA are also considered pivotal for cross-species predictions of anthrax vaccine efficacy in humans, for whom clinical efficacy studies are either impractical or ethically infeasible (3,4) (http://www.fda .gov/AdvisoryCommittees/CommitteesMeetingMaterials/Blood VaccinesandOtherBiologics/VaccinesandRelatedBiologicalProducts AdvisoryCommittee/ucm239733.htm). Anti-PA IgG and TNA levels, however, are but one part of the spectrum of humoral and cell-mediated immune responses that may contribute to protection. The COP for anthrax may differ depending on vaccine formulations, schedules, and routes of administration (5)(6)(7)(8)(9)(10).
The U.S.-licensed anthrax vaccine adsorbed (AVA) (BioThrax) was approved in 1970 for the prevention of anthrax in humans (11)(12)(13)(14). The original regimen for AVA was a subcutaneous (s.c.) six-dose primary schedule at 0, 0.5, 1, 6, 12, and 18 months, with subsequent annual boosters. In May 2012, the U.S. Food and Drug Administration (FDA) approved the AVA regimen as an intramuscular (i.m.) three-dose priming schedule at 0, 1, and 6 months, with boosters at 12 and 18 months and annually thereafter (http://www .fda.gov/BiologicsBloodVaccines/Vaccines/ApprovedProducts /ucm304758.htm). These recent changes in the schedule and ad-ministration route were based on data from the Centers for Disease Control and Prevention Anthrax Vaccine Research Program (AVRP) (12,13). The goals of the AVRP were to improve the AVA safety profile and ensure efficacy while minimizing the number of doses required. The study objectives included determining immunological correlates of protection, documenting vaccine efficacy, and optimizing the vaccination schedule and route of administration (14). Due to the low prevalence of inhalation anthrax and the ethical concerns of conducting an efficacy trial in humans, vaccine efficacy and duration of protection were evaluated in rhesus macaques (Macaca mulatta) (15).
The AVRP nonhuman primate (NHP) study used the 0-, 1-, and 6-month intramuscular priming series (3-i.m.) with a full human dose or saline dilutions of AVA to modulate the immune response. The NHP were challenged with high-dose (median, 504ϫ the 50% lethal dose [LD 50 ]) aerosolized Bacillus anthracis spores at month 12, 30, or 52. The PA-specific humoral and PAstimulated cellular immune response variables were examined during and after the 3-i.m. schedule. In an earlier analysis, the vaccine-induced immune responses were characterized by analysis of variance (ANOVA) and logistic regression. These models were individually fitted to each immunological variable to deter-mine if a variable could predict survival at a specific time point subsequent to the completion of the 3-i.m. priming (15).
In the present analysis, we adopted an alternative strategy to more comprehensively interrogate the AVRP rhesus macaque immunological data to select the best available COP variables. A representative set of 18 immunological responses to PA and 3 response ratios representing Th1/Th2 bias (16,17), interleukin-4 (IL-4) protein to gamma interferon (IFN-␥) protein, IL-4 mRNA to IFN-␥ mRNA, and IL-4-secreting cells to IFN-␥-secreting cells was used to generate a data set of 80 response variables, each of which was considered individually at every available study time point. Together with the AVA dose and time interval between scheduled first vaccination and challenge, we performed variable selection using penalized logistic regressions by two complementing statistical approaches, the stringency of least absolute shrinkage and selection operator (LASSO) and the correlation tolerance of elastic net. Three R software packages, Glmnet (18), Elasticnet (19), and Pensim (20), as well as the Cϩϩ software Bayesian binary regression (BBR) (21), which differ in their optimization algorithms and penalty parameter tuning, were used to avoid having to exclude important predictors and to ensure the selection of a reliable set of COP. The simplest plausible sets of variables were derived (22). The selected variables were evaluated for their correlation with survival, adjusting for variable multicollinearity. The receiver operator characteristic (ROC) area under the curve (AUC) was used as the criterion for testing the accuracy of the prediction model. This comprehensive approach utilized all the available information from the NHP study to select the most practicable and accurate COP for AVA-vaccinated NHP.

MATERIALS AND METHODS
NHP study design, vaccination schedule, and challenge. The rhesus macaque (Macaca mulatta) study design, vaccination, sample schedules, and challenge outcomes have been described in detail (15). Due to the number of vaccinated NHP (n ϭ 114), long duration, and multisite nature of the study, the NHP were organized into 12 groups of 8 to 10 animals plus 2 controls per group, for a total of 137 animals. The control animals were given saline injections instead of AVA but were otherwise treated identically to the vaccinated animals (Table 1). Different groups were vaccinated with the full human dose of AVA (HuAVA), saline-diluted AVA at a 1:5, 1:10, 1:20, or 1:40 dilution, or a saline placebo. The injections were 0.5-ml intramuscular (i.m.) injections administered at 0, 1, and 6 months, followed by aerosol challenge with spores of B. anthracis Ames strain at month 12, 30, or 52.
Humoral and cellular immune responses to AVA. The immunological variables used in the analysis are listed in Table 2. The methods for determining the total anti-PA IgG concentrations (g/ml), TNA titers, anti-PA IgG avidity index (AI), lymphocyte (peripheral blood mononuclear cell [PBMC]) stimulation index (SI), and frequencies of PA IgGspecific B cells, IFN-␥-secreting cells, and IL-4-secreting cells are described in detail elsewhere (15) and are reproduced in brief in the supplemental material. The methods for evaluating the PA-specific induction of gene transcription (mRNA) and translation (protein secretion) for cytokines associated with Th1 (IFN-␥ and IL-2), Th2 (IL-4 and IL-6), and acute-phase (IL-1␤ and tumor necrosis alpha [TNF-␣]) immune responses in PBMC are provided in Tables S1 to S3 and the supplemental methods in the supplemental material).  Data set construction. The construction of data sets is detailed in the supplemental material. Briefly, the data were from control and vaccinated animals that completed the study (Table 1). Except for vaccine dose and the interval between first vaccination and aerosol challenge, the primary data set was constructed, with each variable corresponding to an assay, with measurements taken immediately before each injection, 4 weeks after each injection, and the last available time point prior to aerosol challenge. For an assessment of the relative contributions of humoral and cellular immune responses, ratio variables were generated by dividing the Th2 response-related variables by the Th1 response-related variables. The ratio variables were the ratio of IL-4 mRNA to IFN-␥ mRNA (R_IL4IFNm), the ratio of secreted IL-4 protein to secreted IFN-␥ protein (R_IL4IFNe), and the ratio of the frequency of IL-4-secreting cells to that of IFN-␥-secreting cells (R_IL4IFNeLi). The measurement at each time point was then converted to an individual variable (e.g., anti-PA IgG at month 6 is one variable [IgG_6], and anti-PA IgG at month 7 is a separate variable [IgG_7]). The month-7 time point, which is 1 month after the priming series, was designated "peak," and the last available time point prior to challenge was designated "last" for all NHP. There was a broad range of missing rates for different variables (see Table S4 in the supplemental material). Multiple imputations were used to minimize the bias from missing data (see Table S5 and methods in the supplemental material). Some variables were not imputed due to the high frequency of missing data and therefore were removed from the data set. The final assay variables (n ϭ 80) used in the analysis are listed in Table 2, in addition to the 2 study variables used for vaccine dilution and the time of challenge.
Variable selection by logistic regression penalized by LASSO and elastic net. Variable selection was performed on 20 imputed data sets. Variables that were selected Ն10 times out of the 20 imputed data sets were included in the final variable set for each selection method (see the supplemental material) (23,24). To have the highest confidence that the best correlates were identified, we selected software packages that employed two statistical approaches and that differed in their optimization algorithms and penalty parameter tuning. LASSO and elastic net variable selections were performed in three R packages: Glmnet (18), Elasticnet (19), and Pensim (20), and LASSO variable selection was performed in the Cϩϩ software package BBR (21) (see Table S6 in the supplemental material). LASSO may undergo too-stringent shrinkage and thus ignore important predictors, while elastic net has a grouping effect, selecting important predictors even if they are correlated with each other. Elastic net may consequently select too many predictors, resulting in overfitting in the prediction model (18,19). Repeated (60 times) 10-fold cross-validations were performed to select the best penalty parameters. The optimal sets of variables were selected, with the penalty parameters being where the cross-validation error is minimal or the cross-validated likelihood is maximal. The simplest plausible (parsimonious) sets of variables were selected by applying the "1-standard error rule" (24), choosing the variables with the penalty parameters for which the cross-validation error reached the sum of the minimum cross-validation error and one standard error. The correlation between each selected variable and survival was examined by simple logistic regression in SAS version 9.3 (SAS Institute, Cary, NC) with the original unimputed data set.
Collinearity or multicollinearity diagnosis and evaluation of model performance. Although variables with high rates of missing data are not useful for building practicable models, these variables are still worthy of examination. If the performances of the models that include these variables are compared with those of the models that exclude these variables, the importance of these variables will be known. If these variables do not add much to the performances of the models, these variables can be safely excluded from the models. Otherwise, these variables should be taken into consideration in the future generation of prediction models. Therefore, the model evaluation was performed with the imputed data sets including variables that were heavily imputed.
The emphasis on PA as the primary immunogen and the interrelatedness of the many facets of acquired immunity may result in data from two or more immunoassays being highly correlated with each other. In addition, variables from the same immunoassay but that are measured at different time points may also be highly correlated. Therefore, collinearity or multicollinearity may occur if these variables are present in the same regression models. Logistic regression does not work properly on models containing collinear or multicollinear variables. Therefore, before an evaluation of each variable set for the potential to be a prediction model with good performance, a collinearity or multicollinearity diagnosis was completed (see the supplemental material). To overcome collinearity or multicollinearity, principal component logistic regression (PCLR) was performed to evaluate the performances of the models built with each selected variable set that had collinearity or multicollinearity. For performing PCLR, principal component analysis (PCA) was done using Prcomp in R to generate an eigenvector and score matrix, which were used for logistic regression using glm.fit in R. For variable sets that did not have collinearity or multicollinearity, logistic regression was done using glm.fit in R. The PCLR and logistic regression models are described in detail in the supplemental material. To compare the AUCs between the models, paired permutation tests were performed (25,26), with a Bonferronicorrected significance level of 0.0025 for multiple comparisons.
Practicable and accurate correlates of protection. For further performance evaluation as a COP, PCLR models with multiple variables that had multicollinearity, and logistic models with one variable or two variables that did not have collinearity, were built on the original unimputed data set. PCLR and logistic regressions were performed with the original unimputed data set by glm.fit in R or proc logistic in SAS. To compare the prediction accuracy of the logistic regression models that were built with different variables but on the same original data set, the Roccontrast statement in proc logistic in SAS was used to perform dependent AUC comparisons. To compare the prediction accuracy of the PCLR models with that of the logistic regression models, with both built on the same original data set, paired permutation tests (25,26) were performed. Table 3. Additionally, optimal variable selections were performed with Glmnet (18), Elasticnet (19), the R package Pensim (20), and BBR (21) by choosing the minimum cross-validation error or the maximum cross-validated likelihood, thus minimizing prediction error (see Table S7 in the supplemental material). Additional details for parsimonious and optimal variable selections are included in the supplemental material (see Tables S7 to S13 in the supplemental material). The three variables last anti-PA IgG and SI at months 2 (SI_2) and 6 (SI_6) were selected by all the methods. Among these packages, all the parsimonious variable selections by LASSO (Par_LASSO) chose only these three variables. Parsimonious selection by elastic net with the Elasticnet package (Par_ elastic_Elasticnet) additionally selected both anti-PA IgG and TNA at month 7 (peak). Parsimonious selection by elastic net with the Glmnet package (Par_elastic_Glmnet) further selected IL4e and TNFe at month 1, anti-PA IgG, TNA, and IFNeLi at month 6, and R_IL4IFNeLi and R_IL4IFNm at month 7. The optimal selections chose additional variables IL4eLi at month 1, IL1Be at month 2, IL1Be, TNFe, and R_IL4IFNeLi at month 6, IL1Bm and IL4eLi at month 7, and IL4eLi and R_IL4IFNeLi at the last time point (see Table S7).

Variable selection by logistic regression penalized by LASSO or elastic net. Parsimonious variable selections are summarized in
Simple logistic regression with each of the selected variables showed statistically significant (P Ͻ 0.05) positive correlations between survival and anti-PA IgG at months 6 and 7 (peak), as well as at the last time point, TNA at months 6 and 7, SI at months 2 and 6, and IFN-␥-secreting cell frequency at month 6 ( Table 3). The ratio of the IL-4 to IFN-␥ mRNA levels (R_IL4IFNm) at month 6, which was chosen only in the optimal selections, also had a statistically significant correlation with survival (see Table  S7 in the supplemental material).
Evaluation of survival prediction models with selected sets of variables. The selected variable sets were used to build regression models to test their survival prediction performance. The PCLR models were built with variable sets that had multicollinearity. Logistic regression was applied to the three-variable set (Par_LASSO) that did not have multicollinearity. When the models with the three parsimonious variable sets Par_LASSO, Par_elastic_Elasticnet, and Par_elastic_Glmnet were applied to the imputed data sets, the mean AUCs were 0.8492, 0.8494, and 0.9022, respectively (Table 4; see also Table S14 in the supplemental material). Therefore, the model with variable set Par_LASSO was the most parsimonious among the three. However, this model was limited by the fact that in the unimputed original data, there were 71 out of 137 animals (51.82%) that did not have an SI measurement at month 6. Using the variable SI at month 6 to predict survival with the unimputed data will therefore have low statistical power and may create bias. To test the impact of deleting SI at month 6 on the model performance, a model with SI at month 6 excluded from Par_LASSO was applied to the 20 imputed data sets, and its performance was compared with that of other models. This model had an AUC of 0.8409, which is not statistically significantly different from that of the model with variable set Par_LASSO, suggesting that SI at month 6 did not play a significant role in model performance (Table 4; see also Table S14 in the supplemental material). Therefore, SI at month 6 was excluded in further evaluations.
All together, the model with last anti-PA IgG and SI at month 2 was the most parsimonious. Together with their low missing rates in the unimputed data set, these two variables are worthy of consideration in the generation of practicable prediction models.
Practicable and accurate correlates of protection. The last anti-PA IgG measurement concurrent with aerosol exposure was consistently selected by all the variable selection methods. For situations where a preexposure sample was not available, however, the peak measurement was found to be the best practicable alternative for predicting protection against future exposure. The peak response variables were therefore included in the regression model evaluations.
TNA was statistically significantly correlated with survival but did not provide the same good level of accuracy as anti-PA IgG.  Par_elastic_Glmnet a P value of Ͻ0.05 from Wald chi-square test of the parameter estimate. b •, variables that were selected by all five selection methods. X, variables that were selected by 1 or 2 methods. c Parsimonious selections by LASSO in all the packages selected the same variables that were considered one variable set, Par_LASSO, for subsequent analyses. However, TNA is considered to be species neutral since it does not rely on a species-specific conjugate antibody, allowing a direct comparison in the same assay of serum antibody functional activity across animal species and genera. TNA features prominently in the literature as an immunological bridge to extrapolate data in animals in order to predict anthrax vaccine effectiveness in humans (4,27). Furthermore, anti-PA IgG and TNA levels in NHP were highly correlated (15), and corresponding data in humans are available (12,13,15). Similar to peak anti-PA IgG, peak TNA was included in the regression model evaluations as a potential alternative to last anti-PA IgG. Consequently, the parsimonious variable selection was expanded to include the following five variables: last anti-PA IgG, last TNA, peak anti-PA IgG, peak TNA, and SI at month 2. These variables were evaluated both singly and in combinations in regression models of predicted survival probability in NHP. A collinearity or multicollinearity diagnosis was used to detect col-linearity or multicollinearity and generate the combinations of these variables for the models. Multicollinearity was present in all the variable combinations composed of three or more variables. Among the two-variable combinations, collinearity was not present between SI at month 2 and any of the other four variables, nor between last TNA and peak anti-PA IgG or peak TNA (data not shown). Last anti-PA IgG had only mild collinearity with peak TNA (variance inflation factor [VIF], 2.57; tolerance, 0.39 [the cutoff value is 2.5 for VIF and 0.4 for tolerance]; details are in the supplemental methods and references 12 and 13 in the supplemental material) and peak anti-PA IgG (VIF, 2.68; tolerance, 0.37) (28,29). Because PCLR can be done only with more than three variables, the two-variable logistic regression models containing either last anti-PA IgG and peak TNA or last anti-PA IgG and peak anti-PA IgG were fitted. A PCLR model containing all five variables or the Par_elastic_Elasticnet variable set except SI at month 6 was fitted for comparison with the single-variable and two- Ϫ0.0977 (09300) SI_2 ϩ peak anti-PA IgG ϩ last TNA ϩ peak TNA ϩ last anti-PA IgG variable logistic regression models ( Table 5). The variable set Par_elastic_Glmnet that contains 12 variables was not examined at this step, due to the high frequency of missing data within some of these variables in the original unimputed data set.
Of the single-variable models (Table 5), last anti-PA IgG had the highest AUC (0.8214; 95% confidence interval [CI], 0.7514 to 0.8914), peak anti-PA IgG, peak TNA, and SI at month 2 had lower but not statistically significantly lower AUCs (0.7956, 95% CI, 0.7208 to 0.8703; 0.7918, 95% CI, 0.7158 to 0.8678; and 0.7860, 95% CI, 0.7086 to 0.8633, respectively), and last TNA had statistically significantly lower AUC (0.7270; 95% CI, 0.6526 to 0.8014). Logistic regression using last anti-PA IgG at specific challenge times generated different parameter estimates but similar significant P values compared with those for overall last anti-PA IgG ( Table 5). The different parameter estimates were due to differences in the magnitude of the anti-PA IgG responses at these time points. Last anti-PA IgG had a parameter estimate of 2.1628 ( Table 5), indicating that a 1-log 10 unit (10-fold) increase in last anti-PA IgG can make an 8.70-fold increase in the odds of survival (probability of survival/probability of death). The logistic regression model for last anti-PA IgG is shown in Fig. 1.
In each two-variable model, there was always one variable with a nonsignificant parameter estimate ( Table 5). The 2 two-variable models that contain last anti-PA IgG had AUCs of 0.8247 (95% CI, 0.7546 to 0.8947) and 0.8266 (95% CI, 0.7570 to 0.8963), which were very close to the AUC of the model with last anti-PA IgG only. The model with SI_2 and last anti-PA IgG had an AUC of 0.8424 (95% CI, 0.7758 to 0.9090), the highest accuracy of all the one-variable and two-variable models ( Table 5). The PCLR model containing all five variables had an AUC of 0.8432 (95% CI, 0.7768 to 0.9097), and the PCLR model containing the Par_elas-tic_Elasticnet variable set not including SI at month 6 had an AUC of 0.8403 (95% CI, 0.7744 to 0.9062). Permutation tests showed that compared with the single-variable model containing last anti-PA IgG, only the single-variable model containing last TNA had a statistically significantly lower AUC (0.7270; 95% CI, 0.6526 to 0.8014; P ϭ 0.0009) ( Table 5). The ROC curves for the two-variable models that had a higher AUC than that of the last anti-PA IgG model roughly overlap the ROC curve of the last anti-PA IgG model (Fig. 2).
Overall, anti-PA IgG at the last time point and SI at months 2 and 6 were the three best correlates of protection in rhesus macaques. Of these, the single-variable last anti-PA IgG model had the highest prediction accuracy in this genus and thus is the most appropriate single-variable model for survival predictions (Table 5).

DISCUSSION
Anti-PA antibody levels have been an accepted serological indicator for protection against anthrax since acellular vaccines for human use were first developed (1,(30)(31)(32)(33)(34). More recently, serum TNA has been the focal point for correlates of protection in animal models and estimates of survival probability in humans (4,27,35). Despite this focus, there has been no systematic evaluation of the relationship between PA-stimulated humoral and cell-mediated immune responses, their relative contributions to protection, or their comparative importance when used singly or in combination to predict the probability of survival.
We applied a comprehensive strategy using logistic regression penalized by LASSO or elastic net in four different software applications with dissimilar algorithms, methods, and criteria for penalty parameter tuning. The strategy selected the optimal and par- The two statistical approaches of LASSO and elastic net were applied because LASSO may perform shrinkage too stringently, and the resulting prediction model may exclude important predictors. Alternatively, elastic net may be too lenient in its variable selection, with subsequent overfitting in the prediction model, and consequently, variations within the sample may be disproportionately interpreted as variations in the study population. In order to avoid overfitting of the prediction model, parsimonious sets of variables were derived. The selected variables were then evaluated for correlation with survival. Models without collinearity or multicollinearity issues were evaluated using logistic regression, and models containing variables with multicollinearity were evaluated by PCLR. The AUC was used as the criterion for testing the accuracy of the prediction model. The AUC measures the ability of the model to correctly identify the survival status. AUCs between 0.90 and 1, 0.80 and 0.90, 0.70 and 0.80, and 0.50 and 0.70 represent high, good, moderate, and low discrimination performance of the model, respectively (36). This comprehensive approach utilized all the available information from the study to select the most plausible and practicable COP for AVA-vaccinated rhesus macaques.
Statistically significant positive correlations with survival were identified for anti-PA IgG at months 6 and 7 and the last time point before infectious challenge (last), for TNA at months 6 and 7 (peak), SI at months 2 and 6, and the frequency of IFN-␥-secreting cells at month 6. The ratio of IL-4 mRNA to IFN-␥ mRNA transcription levels at month 6 had a statistically significant negative correlation with survival. The different software applications consistently selected anti-PA IgG responses at challenge (last) and SI at months 2 and 6, thus confirming their importance as immune correlates of protection in rhesus macaques. The frequency of missing data for SI at month 6 was 71 out of 137 (51.82%), causing potential issues of bias and low statistical power. Deleting SI at month 6 from Par_LASSO did not lower the prediction accuracy much. The suitable variables were then used in logistic regression models for survival predictions. Although their performance in the analysis was lower, the peak and last TNA responses were also included in the models, primarily due to their positive correlation with survival and the cross-species utility of TNA for predicting survival probabilities in humans. Therefore, the final selected single-correlate variables were peak anti-PA IgG, peak TNA, last anti-PA IgG, last TNA, and SI at month 2. The final dual-correlate variables were combinations of peak anti-PA IgG, last anti-PA IgG, peak TNA, last TNA, and SI at month 2.
The single-correlate logistic regression models using peak TNA and peak anti-PA IgG responses described the immune status on completion of the 3-i.m. priming series and were valid predictors of survival at later time points. The model with only last anti-PA IgG had good predictive accuracy (AUC, 0.8214; 95% CI, 0.7514 to 0.8914). The predictive accuracy of last TNA was moderate, with an AUC of 0.7270 (95% CI, 0.6526 to 0.8014), which was statistically significantly less than the last anti-PA IgG AUC (P ϭ 0.0009). This difference is likely due to the lower analytical sensitivity of the TNA assay. At the last time point, there were 104 values that were above the limit of detection (LOD) of the anti-PA IgG assay but only 69 that were above the LOD of the TNA assay. Peak anti-PA IgG, peak TNA, and SI at month 2 were moderately accurate predictors of survival, with an AUC of Ͼ0.78.
In the dual-correlate logistic regression models, there was always one parameter estimate that was not statistically significant, indicating that it was not making a significant contribution to the predictive power of the model. Compared with the single-correlate model with last anti-PA IgG only, the predictive accuracy was slightly improved but not statistically significantly different using the dual-correlate models containing last anti-PA IgG or PCLR models containing all five variables or the Par_elastic_Elasticnet variable set except SI at month 6, and it was slightly decreased but not statistically significantly different using all other dual-correlate models. Given the lower technical complexity, higher sample stability, and higher throughput of an enzyme-linked immunosorbent assay (ELISA) than those of a lymphocyte proliferation assay, together with the ability to accurately calibrate ELISA standards and quantify the IgG analyte for each species (15,37), these data confirm that last anti-PA IgG provides an appropriate correlate of protection for cross-species survival predictions (4).
Fay et al. (4) noted in their cross-species meta-analysis that vaccine formulation plays a statistically significant role in the quantitative assessment of correlates of protection. In the AVRP study, both the NHP and human cohorts were treated with AVA in experiments designed to match the schedule and time of measurement as closely as possible. The quantitative model generated from the AVRP NHP data therefore should be used only to bridge to data from human cohorts tested with the same vaccine schedule. The anti-PA IgG concentration at the time of challenge should be a suitable correlate for other PA-based vaccines, although any cross-species bridging should be based on data from clinical and nonclinical studies with matching vaccine formulations, schedules, and times of measurement. Anthrax vaccines that do not rely