Mark-specific hazard ratio model with multivariate continuous marks: an application to vaccine efficacy

Quantifying how single dose Ad26.COV2.S vaccine efficacy depends on Spike sequence features

In the ENSEMBLE randomized, placebo-controlled phase 3 trial (NCT04505722), estimated single-dose Ad26.COV2.S vaccine efficacy (VE) was 56% against moderate to severe-critical COVID-19. SARS-CoV-2 Spike sequences were determined from 484 vaccine and 1,067 placebo recipients who acquired COVID-19. In this set of prespecified analyses, we show that in Latin America, VE was significantly lower against Lambda vs. Reference and against Lambda vs. non-Lambda [family-wise error rate (FWER) p < 0.05]. VE differed by residue match vs. mismatch to the vaccine-insert at 16 amino acid positions (4 FWER p < 0.05; 12 q-value ≤ 0.20); significantly decreased with physicochemical-weighted Hamming distance to the vaccine-strain sequence for Spike, receptor-binding domain, N-terminal domain, and S1 (FWER p < 0.001); differed (FWER ≤ 0.05) by distance to the vaccine strain measured by 9 antibody-epitope escape scores and 4 NTD neutralization-impacting features; and decreased (p = 0.011) with neutralization resistance level to vaccinee sera. VE against severe-critical COVID-19 was stable across most sequence features but lower against the most distant viruses.

Prevention efficacy of the broadly neutralizing antibody VRC01 depends on HIV-1 envelope sequence features

In the Antibody Mediated Prevention (AMP) trials (HVTN 704/HPTN 085 and HVTN 703/HPTN 081), prevention efficacy (PE) of the monoclonal broadly neutralizing antibody (bnAb) VRC01 (vs. placebo) against HIV-1 acquisition diagnosis varied according to the HIV-1 Envelope (Env) neutralization sensitivity to VRC01, as measured by 80% inhibitory concentration (IC80). Here, we performed a genotypic sieve analysis, a complementary approach to gaining insight into correlates of protection that assesses how PE varies with HIV-1 sequence features. We analyzed HIV-1 Env amino acid (AA) sequences from the earliest available HIV-1 RNA-positive plasma samples from AMP participants diagnosed with HIV-1 and identified Env sequence features that associated with PE. The strongest Env AA sequence correlate in both trials was VRC01 epitope distance that quantifies the divergence of the VRC01 epitope in an acquired HIV-1 isolate from the VRC01 epitope of reference HIV-1 strains that were most sensitive to VRC01-mediated neutralization. In HVTN 704/HPTN 085, the Env sequence-based predicted probability that VRC01 IC80 against the acquired isolate exceeded 1 µg/mL also significantly associated with PE. In HVTN 703/HPTN 081, a physicochemical-weighted Hamming distance across 50 VRC01 binding-associated Env AA positions of the acquired isolate from the most VRC01-sensitive HIV-1 strain significantly associated with PE. These results suggest that incorporating mutation scoring by BLOSUM62 and weighting by the strength of interactions at AA positions in the epitope:VRC01 interface can optimize performance of an Env sequence-based biomarker of VRC01 prevention efficacy. Future work could determine whether these results extend to other bnAbs and bnAb combinations.

Baseline malaria infection status and RTS,S/AS01E malaria vaccine efficacy

Background

The only licensed malaria vaccine, RTS,S/AS01 E , confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy (VE).

Methods

1,500 children aged 5-17 months were randomized to receive four different RTS,S/AS01 E regimens or a rabies control vaccine in a phase 2b clinical trial in Ghana and Kenya. We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods in over 36K participant specimens. We performed a post hoc analysis of VE based on malaria infection status at first vaccination and force of infection.

Results

We observed significant and comparable VE (25-43%, 95% CI union 9-53%) against first new infection for all four RTS,S/AS01 E regimens across both follow-up periods (12 and 20 months). Each RTS,S/AS01 E regimen significantly reduced the number of new infections in the 20-month follow-up period (control mean 4.1 vs. RTS,S/AS01 E mean 2.6-3.0). VE against first new infection was significantly higher in participants who were malaria-infected (68%; 95% CI, 50 to 80%) versus uninfected (37%; 95% CI, 23 to 48%) at the first vaccination (P=0.0053) and in participants experiencing greater force of infection between dose 1 and 3 (P=0.059).

Conclusions

All tested dosing regimens blocked some infections to a similar degree. Improved VE in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. ( ClinicalTrials.gov number, NCT03276962 ).

Estimation of vaccine efficacy for variants that emerge after the placebo group is vaccinated

SARS‐CoV‐2 continues to evolve and the vaccine efficacy against variants is challenging to estimate. It is now common in phase III vaccine trials to provide vaccine to those randomized to placebo once efficacy has been demonstrated, precluding a direct assessment of placebo controlled vaccine efficacy after placebo vaccination. In this work, we extend methods developed for estimating vaccine efficacy post placebo vaccination to allow variant specific time varying vaccine efficacy, where time is measured since vaccination. The key idea is to infer counterfactual strain specific placebo case counts by using surveillance data that provide the proportions of the different strains. This blending of clinical trial and observational data allows estimation of strain‐specific time varying vaccine efficacy, or sieve effects, including for strains that emerge after placebo vaccination. The key requirements are that the surveillance strain distribution accurately reflects the strain distribution for a placebo group throughout follow‐up after placebo group vaccination, and that at least one strain is present before and after placebo vaccination. For illustration, we develop a Poisson approach for an idealized design under a rare disease assumption and then use a proportional hazards model to address staggered entry, staggered crossover, and smoothly varying strain specific vaccine efficacy. We evaluate these methods by theoretical work and simulations, and demonstrate that useful estimation of the efficacy profile is possible for strains that emerge after vaccination of the placebo group. An important principle is to incorporate sensitivity analyses to guard against misspecification of the strain distribution.

Associations of human leukocyte antigen with neutralizing antibody titers in a tetravalent dengue vaccine phase 2 efficacy trial in Thailand

The recombinant, live, attenuated, tetravalent dengue vaccine CYD-TDV has shown efficacy against all four dengue serotypes. In this exploratory study (CYD59, NCT02827162), we evaluated potential associations of host human leukocyte antigen (HLA) alleles with dengue antibody responses, CYD-TDV vaccine efficacy, and virologically-confirmed dengue (VCD) cases. Children 4-11 years old, who previously completed a phase 2b efficacy study of CYD-TDV in a single center in Thailand, were included in the study. Genotyping of HLA class I and II loci was performed by next-generation sequencing from DNA obtained from 335 saliva samples. Dengue neutralizing antibody titers (NAb) were assessed as a correlate of risk and protection. Regression analyses were used to assess associations between HLA alleles and NAb responses, vaccine efficacy, and dengue outcomes. Month 13 NAb log geometric mean titers (GMTs) were associated with decreased risk of VCD. In the vaccine group, HLA-DRB1*11 was significantly associated with higher NAb log GMT levels (beta: 0.76; p = 0.002, q = 0.13). Additionally, in the absence of vaccination, HLA associations were observed between the presence of DPB1*03:01 and increased NAb log GMT levels (beta: 1.24; p = 0.005, q = 0.17), and between DPB1*05:01 and reduced NAb log GMT levels (beta: -1.1; p = 0.001, q = 0.07). This study suggests associations of HLA alleles with NAb titers in the context of dengue outcomes. This study was registered with clinicaltrials.gov: NCT02827162.

Estimation of Vaccine Efficacy for Variants that Emerge After the Placebo Group Is Vaccinated

SARS-CoV-2 continues to evolve and the vaccine efficacy against variants is challenging to estimate. It is now common in phase III vaccine trials to provide vaccine to those randomized to placebo once efficacy has been demonstrated, precluding a direct assessment of placebo controlled vaccine efficacy after placebo vaccination. In this work we extend methods developed for estimating vaccine efficacy post placebo vaccination to allow variant specific time varying vaccine efficacy, where time is measured since vaccination. The key idea is to infer counterfactual strain specific placebo case counts by using surveillance data that provide the proportions of the different strains. This blending of clinical trial and observational data allows estimation of strain-specific time varying vaccine efficacy, or sieve effects, including for strains that emergent after placebo vaccination. The key requirements are that surveillance strain distribution accurately reflect the strain distribution for a placebo group, throughout follow-up after placebo group vaccination and that at least one strain is present before and after placebo vaccination. For illustration, we develop a Poisson approach for an idealized design under a rare disease assumption and then use a proportional hazards modeling to better reflect the complexities of field trials with staggered entry, crossover, and smoothly varying strain specific vaccine efficacy We evaluate these by theoretical work and simulations, and demonstrate that useful estimation of the efficacy profile is possible for strains that emerge after vaccination of the placebo group. An important principle is to incorporate sensitivity analyses to guard against mis-specfication of the strain distribution. We also provide an approach for use when genotyping of the infecting strains of the trial participants has not been done.

Two Randomized Trials of Neutralizing Antibodies to Prevent HIV-1 Acquisition

BACKGROUND

Whether a broadly neutralizing antibody (bnAb) can be used to prevent human immunodeficiency virus type 1 (HIV-1) acquisition is unclear.

METHODS

We enrolled at-risk cisgender men and transgender persons in the Americas and Europe in the HVTN 704/HPTN 085 trial and at-risk women in sub-Saharan Africa in the HVTN 703/HPTN 081 trial. Participants were randomly assigned to receive, every 8 weeks, infusions of a bnAb (VRC01) at a dose of either 10 or 30 mg per kilogram (low-dose group and high-dose group, respectively) or placebo, for 10 infusions in total. HIV-1 testing was performed every 4 weeks. The VRC01 80% inhibitory concentration (IC80) of acquired isolates was measured with the TZM-bl assay.

RESULTS

Adverse events were similar in number and severity among the treatment groups within each trial. Among the 2699 participants in HVTN 704/HPTN 085, HIV-1 infection occurred in 32 in the low-dose group, 28 in the high-dose group, and 38 in the placebo group. Among the 1924 participants in HVTN 703/HPTN 081, infection occurred in 28 in the low-dose group, 19 in the high-dose group, and 29 in the placebo group. The incidence of HIV-1 infection per 100 person-years in HVTN 704/HPTN 085 was 2.35 in the pooled VRC01 groups and 2.98 in the placebo group (estimated prevention efficacy, 26.6%; 95% confidence interval [CI], -11.7 to 51.8; P = 0.15), and the incidence per 100 person-years in HVTN 703/HPTN 081 was 2.49 in the pooled VRC01 groups and 3.10 in the placebo group (estimated prevention efficacy, 8.8%; 95% CI, -45.1 to 42.6; P = 0.70). In prespecified analyses pooling data across the trials, the incidence of infection with VRC01-sensitive isolates (IC80 <1 μg per milliliter) per 100 person-years was 0.20 among VRC01 recipients and 0.86 among placebo recipients (estimated prevention efficacy, 75.4%; 95% CI, 45.5 to 88.9). The prevention efficacy against sensitive isolates was similar for each VRC01 dose and trial; VRC01 did not prevent acquisition of other HIV-1 isolates.

CONCLUSIONS

VRC01 did not prevent overall HIV-1 acquisition more effectively than placebo, but analyses of VRC01-sensitive HIV-1 isolates provided proof-of-concept that bnAb prophylaxis can be effective. (Supported by the National Institute of Allergy and Infectious Diseases; HVTN 704/HPTN 085 and HVTN 703/HPTN 081 ClinicalTrials.gov numbers, NCT02716675 and NCT02568215.).

A Hybrid Approach for the Stratified Mark-Specific Proportional Hazards Model with Missing Covariates and Missing Marks, with Application to Vaccine Efficacy Trials

Deployment of the recently licensed CYD-TDV dengue vaccine requires understanding of how the risk of dengue disease in vaccine recipients depends jointly on a host biomarker measured after vaccination (neutralization titer - NAb) and on a "mark" feature of the dengue disease failure event (the amino acid sequence distance of the dengue virus to the dengue sequence represented in the vaccine). The CYD14 phase 3 trial of CYD-TDV measured NAb via case-cohort sampling and the mark in dengue disease failure events, with about a third missing marks. We addressed the question of interest by developing inferential procedures for the stratified mark-specific proportional hazards model with missing covariates and missing marks. Two hybrid approaches are investigated that leverage both augmented inverse probability weighting and nearest neighborhood hot deck multiple imputation. The two approaches differ in how the imputed marks are pooled in estimation. Our investigation shows that NNHD imputation can lead to biased estimation without properly selected neighborhood. Simulations show that the developed hybrid methods perform well with unbiased NNHD imputations from proper neighborhood selection. The new methods applied to CYD14 show that NAb is strongly inversely associated with risk of dengue disease in vaccine recipients, more strongly against dengue viruses with shorter distances.

Assessing trends in vaccine efficacy by pathogen genetic distance

Preventive vaccines are an effective public health intervention for reducing the burden of infectious diseases, but have yet to be developed for several major infectious diseases. Vaccine sieve analysis studies whether and how the efficacy of a vaccine varies with the genetics of the infectious pathogen, which may help guide future vaccine development and deployment. A standard statistical approach to sieve analysis compares the effect of the vaccine to prevent infection and disease caused by pathogen types defined dichotomously as genetically near or far from a reference pathogen strain inside the vaccine construct. For example, near may be defined by amino acid identity at all amino acid positions considered in a multiple alignment and far defined by at least one amino acid difference. An alternative approach is to study the efficacy of the vaccine as a function of genetic distance from a pathogen to a reference vaccine strain where the distance cumulates over the set of amino acid positions. We propose a nonparametric method for estimating and testing the trend in the effect of a vaccine across genetic distance. We illustrate the operating characteristics of the estimator via simulation and apply the method to a recent preventive malaria vaccine efficacy trial.

A unified evaluation of differential vaccine efficacy

Many infectious diseases are well prevented by proper vaccination. However, when a vaccine is not completely efficacious, there is great interest in determining how the vaccine effect differs in subgroups conditional on measured immune responses postvaccination and also according to the type of infecting agent (eg, strain of a virus). The former is often called correlate of protection (CoP) analysis, while the latter has been called sieve analysis. We propose a unified framework for simultaneously assessing CoP and sieve effects of a vaccine in a large Phase III randomized trial. We use flexible parametric models treating times to infection from different agents as competing risks and estimated maximum likelihood to fit the models. The parametric models under competing risks allow for estimation of both cumulative incidence-based contrasts and instantaneous rates. We outline the assumptions with which we can link the observable data to the causal contrasts of interest, propose hypothesis testing procedures, and evaluate our proposed methods in an extensive simulation study.

Prediction of VRC01 neutralization sensitivity by HIV-1 gp160 sequence features

The broadly neutralizing antibody (bnAb) VRC01 is being evaluated for its efficacy to prevent HIV-1 infection in the Antibody Mediated Prevention (AMP) trials. A secondary objective of AMP utilizes sieve analysis to investigate how VRC01 prevention efficacy (PE) varies with HIV-1 envelope (Env) amino acid (AA) sequence features. An exhaustive analysis that tests how PE depends on every AA feature with sufficient variation would have low statistical power. To design an adequately powered primary sieve analysis for AMP, we modeled VRC01 neutralization as a function of Env AA sequence features of 611 HIV-1 gp160 pseudoviruses from the CATNAP database, with objectives: (1) to develop models that best predict the neutralization readouts; and (2) to rank AA features by their predictive importance with classification and regression methods. The dataset was split in half, and machine learning algorithms were applied to each half, each analyzed separately using cross-validation and hold-out validation. We selected Super Learner, a nonparametric ensemble-based cross-validated learning method, for advancement to the primary sieve analysis. This method predicted the dichotomous resistance outcome of whether the IC50 neutralization titer of VRC01 for a given Env pseudovirus is right-censored (indicating resistance) with an average validated AUC of 0.868 across the two hold-out datasets. Quantitative log IC50 was predicted with an average validated R2 of 0.355. Features predicting neutralization sensitivity or resistance included 26 surface-accessible residues in the VRC01 and CD4 binding footprints, the length of gp120, the length of Env, the number of cysteines in gp120, the number of cysteines in Env, and 4 potential N-linked glycosylation sites; the top features will be advanced to the primary sieve analysis. This modeling framework may also inform the study of VRC01 in the treatment of HIV-infected persons.

Viral genetic diversity and protective efficacy of a tetravalent dengue vaccine in two phase 3 trials

Two phase 3 placebo-controlled trials of the CYD-TDV vaccine, evaluated in children aged 2-14 y (CYD14) and 9-16 y (CYD15), demonstrated vaccine efficacy (VE) of 56.5% and 60.8%, respectively, against symptomatic virologically confirmed dengue (VCD). Sieve analyses were conducted to evaluate whether and how VE varied with amino acid sequence features of dengue viruses (DENVs). DENV premembrane/envelope amino acid sequences from VCD endpoint cases were aligned with the vaccine insert sequences, and extensions of the proportional hazards model were applied to assess variation in VE with amino acid mismatch proportion distances from vaccine strains, individual amino acid residues, and phylogenetic genotypes. In CYD14, VE against VCD of any serotype (DENV-Any) decreased significantly with increasing amino acid distance from the vaccine, whereas in CYD15, VE against DENV-Any was distance-invariant. Restricting to the common age range and amino acid distance range between the trials and accounting for differential VE by serotype, however, showed no evidence of VE variation with distance in either trial. In serotype-specific analyses, VE against DENV4 decreased significantly with increasing amino acid distance from the DENV4 vaccine insert and was significantly greater against residue-matched DENV4 at eight signature positions. These effects were restricted to 2- to 8-y-olds, potentially because greater seropositivity of older children at baseline might facilitate a broader protective immune response. The relevance of an antigenic match between vaccine strains and circulating DENVs was also supported by greater estimated VE against serotypes and genotypes for which the circulating DENVs had shorter amino acid sequence distances from the vaccine.

Sieve analysis using the number of infecting pathogens

Assessment of vaccine efficacy as a function of the similarity of the infecting pathogen to the vaccine is an important scientific goal. Characterization of pathogen strains for which vaccine efficacy is low can increase understanding of the vaccine's mechanism of action and offer targets for vaccine improvement. Traditional sieve analysis estimates differential vaccine efficacy using a single identifiable pathogen for each subject. The similarity between this single entity and the vaccine immunogen is quantified, for example, by exact match or number of mismatched amino acids. With new technology, we can now obtain the actual count of genetically distinct pathogens that infect an individual. Let F be the number of distinct features of a species of pathogen. We assume a log-linear model for the expected number of infecting pathogens with feature "f," <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>f</mml:mi> <mml:mo>=</mml:mo> <mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mo>…</mml:mo> <mml:mo>,</mml:mo> <mml:mi>F</mml:mi></mml:math> . The model can be used directly in studies with passive surveillance of infections where the count of each type of pathogen is recorded at the end of some interval, or active surveillance where the time of infection is known. For active surveillance, we additionally assume that a proportional intensity model applies to the time of potentially infectious exposures and derive product and weighted estimating equation (WEE) estimators for the regression parameters in the log-linear model. The WEE estimator explicitly allows for waning vaccine efficacy and time-varying distributions of pathogens. We give conditions where sieve parameters have a per-exposure interpretation under passive surveillance. We evaluate the methods by simulation and analyze a phase III trial of a malaria vaccine.

Sieve analysis of breakthrough HIV-1 sequences in HVTN 505 identifies vaccine pressure targeting the CD4 binding site of Env-gp120

Although the HVTN 505 DNA/recombinant adenovirus type 5 vector HIV-1 vaccine trial showed no overall efficacy, analysis of breakthrough HIV-1 sequences in participants can help determine whether vaccine-induced immune responses impacted viruses that caused infection. We analyzed 480 HIV-1 genomes sampled from 27 vaccine and 20 placebo recipients and found that intra-host HIV-1 diversity was significantly lower in vaccine recipients (P ≤ 0.04, Q-values ≤ 0.09) in Gag, Pol, Vif and envelope glycoprotein gp120 (Env-gp120). Furthermore, Env-gp120 sequences from vaccine recipients were significantly more distant from the subtype B vaccine insert than sequences from placebo recipients (P = 0.01, Q-value = 0.12). These vaccine effects were associated with signatures mapping to CD4 binding site and CD4-induced monoclonal antibody footprints. These results suggest either (i) no vaccine efficacy to block acquisition of any viral genotype but vaccine-accelerated Env evolution post-acquisition; or (ii) vaccine efficacy against HIV-1s with Env sequences closest to the vaccine insert combined with increased acquisition due to other factors, potentially including the vaccine vector.

Estimation of Stratified Mark-Specific Proportional Hazards Models under Two-Phase Sampling with Application to HIV Vaccine Efficacy Trials

An objective of preventive HIV vaccine efficacy trials is to understand how vaccine-induced immune responses to specific protein sequences of HIV-1 associate with subsequent infection with specific sequences of HIV, where the immune response biomarkers are measured in vaccine recipients via a two-phase sampling design. Motivated by this objective, we investigate the stratified mark-specific proportional hazards model under two-phase biomarker sampling, where the mark is the genetic distance of an infecting HIV-1 sequence to an HIV-1 sequence represented inside the vaccine. Estimation and inference procedures based on inverse probability weighting of complete-cases and on augmented inverse probability weighting of complete-cases are developed. Asymptotic properties of the estimators are derived and their finite-sample performances are examined in simulation studies. The methods are shown to have satisfactory performance, and are applied to the RV144 vaccine trial to assess whether immune response correlates of HIV-1 infection are stronger for HIV-1 infecting sequences similar to the vaccine than for sequences distant from the vaccine. Augmented inverse probability weighting; auxiliary variables; case-cohort design; censored failure time; competing risks; HIV vaccine efficacy trial.

Mark-specific hazard ratio model with missing multivariate marks

An objective of randomized placebo-controlled preventive HIV vaccine efficacy (VE) trials is to assess the relationship between vaccine effects to prevent HIV acquisition and continuous genetic distances of the exposing HIVs to multiple HIV strains represented in the vaccine. The set of genetic distances, only observed in failures, is collectively termed the 'mark.' The objective has motivated a recent study of a multivariate mark-specific hazard ratio model in the competing risks failure time analysis framework. Marks of interest, however, are commonly subject to substantial missingness, largely due to rapid post-acquisition viral evolution. In this article, we investigate the mark-specific hazard ratio model with missing multivariate marks and develop two inferential procedures based on (i) inverse probability weighting (IPW) of the complete cases, and (ii) augmentation of the IPW estimating functions by leveraging auxiliary data predictive of the mark. Asymptotic properties and finite-sample performance of the inferential procedures are presented. This research also provides general inferential methods for semiparametric density ratio/biased sampling models with missing data. We apply the developed procedures to data from the HVTN 502 'Step' HIV VE trial.

Predicting Overall Vaccine Efficacy in a New Setting by Re-Calibrating Baseline Covariate and Intermediate Response Endpoint Effect Modifiers of Type-Specific Vaccine Efficacy

We develop a transport formula for predicting overall cumulative vaccine efficacy through time t (VE(t)) to prevent clinically significant infection with a genetically diverse pathogen (e.g., HIV infection) in a new setting for which a Phase III preventive vaccine efficacy trial that would directly estimate VE(t) has not yet been conducted. The formula integrates data from (1) a previous Phase III trial, (2) a Phase I/II immune response biomarker endpoint trial in the new setting where a follow-up Phase III trial is planned, (3) epidemiological data on background HIV infection incidence in the new setting; and (4) genomic epidemiological data on HIV sequence distributions in the previous and new settings. For (1), the randomized vaccine versus placebo Phase III trial yields estimates of vaccine efficacy to prevent particular genotypes of HIV in participant subgroups defined by baseline covariates X and immune responses to vaccination S(1) measured at a fixed time point τ (potential outcomes if assigned vaccine); often one or more immune responses to vaccination are available that modify genotype-specific vaccine efficacy. The formula focuses on subgroups defined by X and S(1) and being at-risk for HIV infection at τ under both the vaccine and placebo treatment assignments. For (2), the Phase I/II trial tests the same vaccine in a new setting, or a refined new vaccine in the same or new setting, and measures the same baseline covariates and immune responses as the original Phase III trial. For (3), epidemiological data in the new setting are used to project overall background HIV infection rates in the baseline covariate subgroups in the planned Phase III trial, hence re-calibrating for HIV incidence differences in the two settings; whereas for (4), data bases of HIV sequences measured from HIV infected individuals are used to re-calibrate for differences in the distributions of the circulating HIV genotypes in the two settings. The transport formula incorporates a user-specified bridging assumption function that measures differences in HIV genotype-specific conditional biological-susceptibility vaccine efficacies in the two settings, facilitating a sensitivity analysis. We illustrate the transport formula with application to HIV Vaccine Trials Network (HVTN) research. One application of the transport formula is to use predicted VE(t) as a rational criterion for ranking a set of candidate vaccines being studied in Phase I/II trials for their priority for down-selection into the follow-up Phase III trial.

Leaky vaccines protect highly exposed recipients at a lower rate: implications for vaccine efficacy estimation and sieve analysis

“Leaky” vaccines are those for which vaccine-induced protection reduces infection rates on a per-exposure basis, as opposed to “all-or-none” vaccines, which reduce infection rates to zero for some fraction of subjects, independent of the number of exposures. Leaky vaccines therefore protect subjects with fewer exposures at a higher effective rate than subjects with more exposures. This simple observation has serious implications for analysis methodologies that rely on the assumption that the vaccine effect is homogeneous across subjects. We argue and show through examples that this heterogeneous vaccine effect leads to a violation of the proportional hazards assumption, to incomparability of infected cases across treatment groups, and to nonindependence of the distributions of the competing failure processes in a competing risks setting. We discuss implications for vaccine efficacy estimation, correlates of protection analysis, and mark-specific efficacy analysis (also known as sieve analysis).