Improved Analysis of Phage ImmunoPrecipitation Sequencing (PhIP-Seq) Data Using a Z-score Algorithm

Phage ImmunoPrecipitation Sequencing (PhIP-Seq) is a massively multiplexed, phage-display based methodology for analyzing antibody binding specificities, with several advantages over existing techniques, including the uniformity and completeness of proteomic libraries, as well as high sample throughput and low cost. Data generated by the PhIP-Seq assay are unique in many ways. The only published analytical approach for these data suffers from important limitations. Here, we propose a new statistical framework with several improvements. Using a set of replicate mock immunoprecipitations (negative controls lacking antibody input) to generate background binding distributions, we establish a statistical model to quantify antibody-dependent changes in phage clone abundance. Our approach incorporates robust regression of experimental samples against the mock IPs as a means to calculate the expected phage clone abundance, and provides a generalized model for calculating each clone9s expected abundance-associated standard deviation. In terms of bias removal and detection sensitivity, we demonstrate that this z-score algorithm outperforms the previous approach. Further, in a large cohort of autoantibody-defined Sjogren9s Syndrome (SS) patient sera, PhIP-Seq robustly identified Ro52, Ro60, and SSB/La as known autoantigens associated with SS. In an effort to identify novel SS-specific binding specificities, SS z-scores were compared with z-scores obtained by screening Ro-positive sera from patients with systemic lupus erythematosus (SLE). This analysis did not yield any commonly targeted SS-specific autoantigens, suggesting that if they exist at all, their epitopes are likely to be discontinuous or post-translationally modified. In summary, we have developed an improved algorithm for PhIP-Seq data analysis, which was validated using a large set of sera with clinically characterized autoantibodies. This z-score approach will substantially improve the ability of PhIP-Seq to detect and interpret antibody binding specificities. The associated Python code is freely available for download here: https://github.com/LarmanLab/PhIP-Seq-Analyzer.