Ritornello: high fidelity control-free chromatin immunoprecipitation peak calling.

With the advent of next generation high-throughput DNA sequencing technologies, omics experiments have become the mainstay for studying diverse biological effects on a genome wide scale. Chromatin immunoprecipitation (ChIP-seq) is the omics technique that enables genome wide localization of transcription factor (TF) binding or epigenetic modification events. Since the inception of ChIP-seq in 2007, many methods have been developed to infer ChIP-target binding loci from the resultant reads after mapping them to a reference genome. However, interpreting these data has proven challenging, and as such these algorithms have several shortcomings, including susceptibility to false positives due to artifactual peaks, poor localization of binding sites and the requirement for a total DNA input control which increases the cost of performing these experiments. We present Ritornello, a new approach for finding TF-binding sites in ChIP-seq, with roots in digital signal processing that addresses all of these problems. We show that Ritornello generally performs equally or better than the peak callers tested and recommended by the ENCODE consortium, but in contrast, Ritornello does not require a matched total DNA input control to avoid false positives, effectively decreasing the sequencing cost to perform ChIP-seq. Ritornello is freely available at https://github.com/KlugerLab/Ritornello.

Gating mass cytometry data by deep learning.

MOTIVATION: Mass cytometry or CyTOF is an emerging technology for high-dimensional multiparameter single cell analysis that overcomes many limitations of fluorescence-based flow cytometry. New methods for analyzing CyTOF data attempt to improve automation, scalability, performance and interpretation of data generated in large studies. Assigning individual cells into discrete groups of cell types (gating) involves time-consuming sequential manual steps, untenable for larger studies. RESULTS: We introduce DeepCyTOF, a standardization approach for gating, based on deep learning techniques. DeepCyTOF requires labeled cells from only a single sample. It is based on domain adaptation principles and is a generalization of previous work that allows us to calibrate between a target distribution and a source distribution in an unsupervised manner. We show that DeepCyTOF is highly concordant (98%) with cell classification obtained by individual manual gating of each sample when applied to a collection of 16 biological replicates of primary immune blood cells, even when measured across several instruments. Further, DeepCyTOF achieves very high accuracy on the semi-automated gating challenge of the FlowCAP-I competition as well as two CyTOF datasets generated from primary immune blood cells: (i) 14 subjects with a history of infection with West Nile virus (WNV), (ii) 34 healthy subjects of different ages. We conclude that deep learning in general, and DeepCyTOF specifically, offers a powerful computational approach for semi-automated gating of CyTOF and flow cytometry data. AVAILABILITY AND IMPLEMENTATION: Our codes and data are publicly available at https://github.com/KlugerLab/deepcytof.git. CONTACT: yuval.kluger@yale.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Removal of batch effects using distribution-matching residual networks.

MOTIVATION: Sources of variability in experimentally derived data include measurement error in addition to the physical phenomena of interest. This measurement error is a combination of systematic components, originating from the measuring instrument and random measurement errors. Several novel biological technologies, such as mass cytometry and single-cell RNA-seq (scRNA-seq), are plagued with systematic errors that may severely affect statistical analysis if the data are not properly calibrated. RESULTS: We propose a novel deep learning approach for removing systematic batch effects. Our method is based on a residual neural network, trained to minimize the Maximum Mean Discrepancy between the multivariate distributions of two replicates, measured in different batches. We apply our method to mass cytometry and scRNA-seq datasets, and demonstrate that it effectively attenuates batch effects. AVAILABILITY AND IMPLEMENTATION: our codes and data are publicly available at https://github.com/ushaham/BatchEffectRemoval.git. CONTACT: yuval.kluger@yale.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Algorithm 971: An Implementation of a Randomized Algorithm for Principal Component Analysis.

Recent years have witnessed intense development of randomized methods for low-rank approximation. These methods target principal component analysis and the calculation of truncated singular value decompositions. The present article presents an essentially black-box, foolproof implementation for Mathworks' MATLAB, a popular software platform for numerical computation. As illustrated via several tests, the randomized algorithms for low-rank approximation outperform or at least match the classical deterministic techniques (such as Lanczos iterations run to convergence) in basically all respects: accuracy, computational efficiency (both speed and memory usage), ease-of-use, parallelizability, and reliability. However, the classical procedures remain the methods of choice for estimating spectral norms and are far superior for calculating the least singular values and corresponding singular vectors (or singular subspaces).

Multiplexed (18-Plex) Measurement of Signaling Targets and Cytotoxic T Cells in Trastuzumab-Treated Patients using Imaging Mass Cytometry.

PURPOSE: Imaging mass cytometry (IMC) uses metal-conjugated antibodies to provide multidimensional, objective measurement of protein targets. We used this high-throughput platform to perform an 18-plex assessment of HER2 ICD/ECD, cytotoxic T-cell infiltration and other structural and signaling proteins in a cohort of patients treated with trastuzumab to discover associations with trastuzumab benefit. EXPERIMENTAL DESIGN: An antibody panel for detection of 18 targets (pan-cytokeratin, HER2 ICD, HER2 ECD, CD8, vimentin, cytokeratin 7, ß-catenin, HER3, MET, EGFR, ERK 1-2, MEK 1-2, PTEN, PI3K p110 α, Akt, mTOR, Ki67, and Histone H3) was used with a selection of trastuzumab-treated patients from the Hellenic Cooperative Oncology Group 10/05 trial (n = 180), and identified a case-control series. RESULTS: Patients that recurred after adjuvant treatment with trastuzumab trended toward a decreased fraction of HER2 ECD pixels over threshold compared with cases without recurrence (P = 0.057). After exclusion of the lowest HER2 expressers, 5-year recurrence events were associated with reduced total extracellular domain (ECD)/intracellular domain (ICD) ratio intensity in tumor (P = 0.044). These observations are consistent with our previous work using quantitative immunofluorescence, but represent the proof on identical cell content. We also describe the association of the ECD of HER2 with CD8 T-cell infiltration on the same slide. CONCLUSIONS: The proximity of CD8 cells as a function of the expression of the ECD of HER2 provides further evidence for the role of the immune system in the mechanism of action of trastuzumab.