NMF-RI: blind spectral unmixing of highly mixed multispectral flow and image cytometry data

MOTIVATION Recent advances in multiplex immunostaining and multispectral cytometry have opened the door to simultaneously visualizing an unprecedented number of biomarkers both in liquid and solid samples. Properly unmixing fluorescent emissions is a challenging task, which normally requires the characterization of the individual fluorochromes from control samples. As the number of fluorochromes increases, the cost in time and use of reagents becomes prohibitively high. Here we present a fully-unsupervised blind spectral unmixing method for the separation of fluorescent emissions in highly mixed spectral data, without the need for control samples. To this end, we extend an existing method based on Non-negative Matrix Factorization, and introduce several critical improvements: initialization based on the theoretical spectra, automated selection of 'sparse' data and use of a re-initialized multi-layer optimizer. RESULTS Our algorithm is exhaustively tested using synthetic data to study its robustness against different levels of colocalization, signal to noise ratio, spectral resolution, and the effect of errors in the initialization of the algorithm. Then we compare the performance of our method to that of traditional spectral unmixing algorithms using novel multispectral flow and image cytometry systems. In all cases, we show that our blind unmixing algorithm performs robust unmixing of highly spatially and spectrally mixed data with an unprecedently low computational cost. In summary, we present the first use of a blind unmixing method in multispectral flow and image cytometry, opening the door to the widespread use of our method to efficiently pre-process multiplex immunostaining samples without the need of experimental controls. AVAILABILITY https://drive.google.com/open?id=11Bb5CwKshAYNS91AFy_mRPB0b3gUhQ6- contains the source code and all datasets used in this manuscript. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Ewa M. Goldys,et al.  Statistically strong label-free quantitative identification of native fluorophores in a biological sample , 2017, Scientific Reports.

[2]  Naoto Yokoya,et al.  Coupled Nonnegative Matrix Factorization Unmixing for Hyperspectral and Multispectral Data Fusion , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[3]  M. Atkins,et al.  Predictive biomarkers for checkpoint inhibitor-based immunotherapy. , 2016, The Lancet. Oncology.

[4]  Carlos Ortiz-de-Solorzano,et al.  Efficient Blind Spectral Unmixing of Fluorescently Labeled Samples Using Multi-Layer Non-Negative Matrix Factorization , 2013, PloS one.

[5]  Yuan Yan Tang,et al.  Matrix-Vector Nonnegative Tensor Factorization for Blind Unmixing of Hyperspectral Imagery , 2017, IEEE Transactions on Geoscience and Remote Sensing.

[6]  Yannick Deville,et al.  Inertia-Constrained Pixel-by-Pixel Nonnegative Matrix Factorisation: a Hyperspectral Unmixing Method Dealing with Intra-class Variability , 2017, Remote. Sens..

[7]  Fabian J Theis,et al.  Blind source separation techniques for the decomposition of multiply labeled fluorescence images. , 2009, Biophysical journal.

[8]  Binjie Qin,et al.  Two-hierarchical nonnegative matrix factorization distinguishing the fluorescent targets from autofluorescence for fluorescence imaging , 2015, BioMedical Engineering OnLine.

[9]  Hairong Qi,et al.  Endmember Extraction From Highly Mixed Data Using Minimum Volume Constrained Nonnegative Matrix Factorization , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[10]  Patrik O. Hoyer,et al.  Non-negative Matrix Factorization with Sparseness Constraints , 2004, J. Mach. Learn. Res..

[11]  Mei Jiang,et al.  Validation of multiplex immunofluorescence panels using multispectral microscopy for immune-profiling of formalin-fixed and paraffin-embedded human tumor tissues , 2017, Scientific Reports.

[12]  M Roederer,et al.  Spectral compensation for flow cytometry: visualization artifacts, limitations, and caveats. , 2001, Cytometry.

[13]  Andrzej Cichocki,et al.  New Algorithms for Non-Negative Matrix Factorization in Applications to Blind Source Separation , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[14]  Antonio J. Plaza,et al.  On Endmember Identification in Hyperspectral Images Without Pure Pixels: A Comparison of Algorithms , 2011, Journal of Mathematical Imaging and Vision.

[15]  Hassan Ghassemian,et al.  Spectral Unmixing of Hyperspectral Imagery Using Multilayer NMF , 2014, IEEE Geoscience and Remote Sensing Letters.

[16]  Yuan Yan Tang,et al.  Hypergraph-Regularized Sparse NMF for Hyperspectral Unmixing , 2016, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[17]  Bartek Rajwa,et al.  Generalized unmixing model for multispectral flow cytometry utilizing nonsquare compensation matrices , 2013, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[18]  C. Le Roy,et al.  Flow cytometry APC‐tandem dyes are degraded through a cell‐dependent mechanism , 2009, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[19]  David J. Field,et al.  What Is the Goal of Sensory Coding? , 1994, Neural Computation.

[20]  Y. Saeys,et al.  Computational flow cytometry: helping to make sense of high-dimensional immunology data , 2016, Nature Reviews Immunology.

[21]  T. Zimmermann Spectral imaging and linear unmixing in light microscopy. , 2005, Advances in biochemical engineering/biotechnology.

[22]  Yuval Garini,et al.  Spectral imaging: Principles and applications , 2006, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[23]  Xiaodong Tao,et al.  Autofluorescence Removal by Non-Negative Matrix Factorization , 2011, IEEE Transactions on Image Processing.

[24]  Chen Hu,et al.  Target/Background Classification Regularized Nonnegative Matrix Factorization for Fluorescence Unmixing , 2016, IEEE Transactions on Instrumentation and Measurement.

[25]  H. Sebastian Seung,et al.  Algorithms for Non-negative Matrix Factorization , 2000, NIPS.