Flexible variable selection in the presence of missing data.

In many applications, it is of interest to identify a parsimonious set of features, or panel, from multiple candidates that achieves a desired level of performance in predicting a response. This task is often complicated in practice by missing data arising from the sampling design or other random mechanisms. Most recent work on variable selection in missing data contexts relies in some part on a finite-dimensional statistical model, e.g., a generalized or penalized linear model. In cases where this model is misspecified, the selected variables may not all be truly scientifically relevant and can result in panels with suboptimal classification performance. To address this limitation, we propose a nonparametric variable selection algorithm combined with multiple imputation to develop flexible panels in the presence of missing-at-random data. We outline strategies based on the proposed algorithm that achieve control of commonly used error rates. Through simulations, we show that our proposal has good operating characteristics and results in panels with higher classification and variable selection performance compared to several existing penalized regression approaches in cases where a generalized linear model is misspecified. Finally, we use the proposed method to develop biomarker panels for separating pancreatic cysts with differing malignancy potential in a setting where complicated missingness in the biomarkers arose due to limited specimen volumes.

CMV viral load kinetics as surrogate endpoints after allogeneic transplantation.

BACKGROUNDViral load (VL) surrogate endpoints transformed development of HIV and hepatitis C therapeutics. Surrogate endpoints for CMV-related morbidity and mortality could advance development of antiviral treatments. Although observational data support using CMV VL as a trial endpoint, randomized controlled trials (RCTs) demonstrating direct associations between virological markers and clinical endpoints are lacking.METHODSWe performed CMV DNA PCR on frozen serum samples from the only placebo-controlled RCT of ganciclovir for early treatment of CMV after hematopoietic cell transplantation (HCT). We used established criteria to assess VL kinetics as surrogates for CMV disease or death by weeks 8, 24, and 48 after randomization and quantified antiviral effects captured by each marker. We used ensemble-based machine learning to assess the predictive ability of VL kinetics and performed this analysis on a ganciclovir prophylaxis RCT for validation.RESULTSVL suppression with ganciclovir reduced cumulative incidence of CMV disease and death for 20 years after HCT. Mean VL, peak VL, and change in VL during the first 5 weeks of treatment fulfilled the Prentice definition for surrogacy, capturing more than 95% of ganciclovir's effect, and yielded highly sensitive and specific predictions by week 48. In the prophylaxis trial, the viral shedding rate satisfied the Prentice definition for CMV disease by week 24.CONCLUSIONSOur results support using CMV VL kinetics as surrogates for CMV disease, provide a framework for developing CMV preventative and therapeutic agents, and support reductions in VL as the mechanism through which antivirals reduce CMV disease.FUNDINGMerck Sharp & Dohme Corp., a subsidiary of Merck & Co., Inc.

Antibody Fc effector functions and IgG3 associate with decreased HIV-1 risk.

HVTN 505 is a preventative vaccine efficacy trial testing DNA followed by recombinant adenovirus serotype 5 (rAd5) in circumcised, Ad5-seronegative men and transgendered persons who have sex with men in the United States. Identified immune correlates of lower HIV-1 risk and a virus sieve analysis revealed that, despite lacking overall efficacy, vaccine-elicited responses exerted pressure on infecting HIV-1 viruses. To interrogate the mechanism of the antibody correlate of HIV-1 risk, we examined antigen-specific antibody recruitment of Fcγ receptors (FcγRs), antibody-dependent cellular phagocytosis (ADCP), and the role of anti-envelope (anti-Env) IgG3. In a prespecified immune correlates analysis, antibody-dependent monocyte phagocytosis and antibody binding to FcγRIIa correlated with decreased HIV-1 risk. Follow-up analyses revealed that anti-Env IgG3 breadth correlated with reduced HIV-1 risk, anti-Env IgA negatively modified infection risk by Fc effector functions, and that vaccine recipients with a specific FcγRIIa single-nucleotide polymorphism locus had a stronger correlation with decreased HIV-1 risk when ADCP, Env-FcγRIIa, and IgG3 binding were high. Additionally, FcγRIIa engagement correlated with decreased viral load setpoint in vaccine recipients who acquired HIV-1. These data support a role for vaccine-elicited anti-HIV-1 Env IgG3, antibody engagement of FcRs, and phagocytosis as potential mechanisms for HIV-1 prevention.

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.

DRUG-NEM: Optimizing drug combinations using single-cell perturbation response to account for intratumoral heterogeneity.

An individual malignant tumor is composed of a heterogeneous collection of single cells with distinct molecular and phenotypic features, a phenomenon termed intratumoral heterogeneity. Intratumoral heterogeneity poses challenges for cancer treatment, motivating the need for combination therapies. Single-cell technologies are now available to guide effective drug combinations by accounting for intratumoral heterogeneity through the analysis of the signaling perturbations of an individual tumor sample screened by a drug panel. In particular, Mass Cytometry Time-of-Flight (CyTOF) is a high-throughput single-cell technology that enables the simultaneous measurements of multiple ([Formula: see text]40) intracellular and surface markers at the level of single cells for hundreds of thousands of cells in a sample. We developed a computational framework, entitled Drug Nested Effects Models (DRUG-NEM), to analyze CyTOF single-drug perturbation data for the purpose of individualizing drug combinations. DRUG-NEM optimizes drug combinations by choosing the minimum number of drugs that produce the maximal desired intracellular effects based on nested effects modeling. We demonstrate the performance of DRUG-NEM using single-cell drug perturbation data from tumor cell lines and primary leukemia samples.