Sage Journals: Discover world-class research

Abstract

An immune cell's phenotype expresses through its high-dimensional marker signature. Cluster analyses of data from high-throughput mass and flow cytometry marker panels permit discovery of previously undescribed immune cell phenotypes. Impactful reporting of new phenotypes demands low-dimensional visualization tools that preserve with integrity phenotypes' original high-dimensional structure. For this purpose, we introduce penalized supervised star plots. As designed and as we demonstrate, penalized supervised star plots are two-dimensional projections that tend to preserve separation of clusters as well as information on the relative contributions of various markers in differentiating phenotypes. The new method is robust to markers that do not differentiate phenotypes at all, as shown in a challenge data set. Results include comparison with other popular procedures. Penalized supervised star plots incorporate cross-validation to permit portability of estimated optimal projections to new samples. Supervised star plots are further illustrated with a featured influenza-specific T cell data set as well as a peripheral blood mononuclear cell phenotyping data set.

Get full access to this article

View all access options for this article.

References

Abraham

. Radviz: A R Package for Multi-Dimensional Data Visualization. The Comprehensive R Archive Network (CRAN), 2016. Available at https://cran.cnr.berkeley.edu/ (Last accessed January 23, 2019 ).

Abraham

. 2016. Visualizing multivariate data with Radviz. Retrieved from: https://cran.r-project.org/web/packages/Radviz/vignettes/single_cell_projections.html (Accessed on December 5, 2018 ).

Amir

, Davis

, Tadmor

, et al. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat Biotechnol, 2013; 31:545–552.

Baker

. 1,500 scientists lift the lid on reproducibility. Nature, 2016; 533:452–454.

Bruggner

, Bodenmiller

, Dill

, et al. Automated identification of stratifying signatures in cellular subpopulations. Proc Natl Acad Sci USA, 2014; 111:E2770-7.

Chiu

, Talhouk

, and Liu

. diceR: Diverse cluster ensemble in R. Comprehensive R Archive Network (CRAN), 2018. (Accessed on July 11, 2017 ).

Duchon

Splines minimizing rotation-invariant semi-norms in Sobolev spaces. In: Schempp

, Zeller

, eds. Constructive Theory of Functions of Several Variables. Berlin, Heidelberg: Springer-Verlag (Dold A, Eckmann B, editors. Lecture Notes in Mathematics), 1977.

Hansen

. The risk of James–Stein and lasso shrinkage. Econom Rev, 2016; 35:1456–1470.

Hastie

, Tibshirani

, and Friedman

. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. 2nd ed. New York: Springer, 2009:84–95.

10.

Hastie

, Tibshirani

, and Friedman

. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. 2nd ed. New York: Springer, 2009:501–528.

11.

Hoffman

, Grinstein

, Marx

, et al. DNA visual and analytic data mining. In: Proceedings Visualization'97 (Cat No 97CB36155). Phoenix, AZ: IEEE, 1997:437–441.

12.

Hoffman

, Grinstein

, and Pinkney

. Dimensional anchors: a graphic primitive for multidimensional multivariate information visualizations. In: Proceedings of the 1999 Workshop on New Paradigms in Information Visualization and Manipulation in Conjunction with the Eighth ACM International Conference on Information and Knowledge Management—NPIVM'99. New York, NY: ACM Press, 1999:9–16.

13.

Jombart

. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics, 2008; 24:1403–1405.

14.

Jombart

, Devillard

, and Balloux

. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet, 2010; 11:94.

15.

Kandogan

. Visualizing multi-dimensional clusters, trends, and outliers using star coordinates. In: Proceedings of the Seventh ACM SIGKDD International Conference on Knowledge Discovery and Data Mining—KDD'01. New York, NY: ACM Press, 2001:107–116.

16.

Kim

J-H

. Estimating classification error rate: repeated cross-validation, repeated hold-out and bootstrap. Comput Stat Data Anal, 2009; 53:3735–3745.

17.

Kohonen

. The self-organizing map. Proc IEEE, 1990; 78:1464–1480.

18.

Lemon

. plotrix: a package in the red light district of R. R-News, 2006; 6:8–12.

19.

Lin

, Gupta

, and Maecker

. Intracellular cytokine staining on PBMCs using CyTOF™ mass cytometry. Bio Protoc, 2015; 5:e1370.

20.

van der Maaten

, and Hinton

. Visualizing data using t-SNE. J Machine Learn Res, 2009; 9:2579–2605.

21.

Maurya

. JPEN: Covariance and Inverse Covariance Matrix Estimation Using Joint Penalty. Comprehensive R Archive Network (CRAN), 2015. Available at https://cran.cnr.berkeley.edu/ (Last accessed January 23, 2019 ).

22.

Novomestky

. matrixcalc: Collection of Functions for Matrix Calculations. Comprehensive R Archive Network (CRAN), 2012. Available at https://cran.cnr.berkeley.edu/ (Last accessed January 23, 2019 ).

23.

Orloci

. An agglomerative method for classification of plant communities. J Ecol, 1967; 55:193–206.

24.

Pearson

. L III. On lines and planes of closest fit to systems of points in space. Philos Mag Ser, 1901; 6 2:559–572.

25.

Peterson

, and Ford

. Random matrix theory and covariance matrix filtering for cancer gene expression. In: Peterson

, Masulli

, Russo

, eds. Computational Intelligence Methods for Bioinformatics and Biostatistics. Berlin, Heidelberg: Springer Berlin Heidelberg ( Hutchison

, Kanade

, Kittler

, et al., eds. Lecture Notes in Computer Science), 2013:173–184.

26.

Qiu

, and Joe

. clusterGeneration: Random Cluster Generation (with Specified Degree of Separation). Comprehensive R Archive Network (CRAN), 2015. Available at https://cran.cnr.berkeley.edu/ (Last accessed January 23, 2019 ).

27.

Rousseeuw

. Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. J Comput Appl Math, 1987; 20:53–65.

28.

Schuetzenmeister

, and Dufey

. VCA: Variance Component Analysis. Comprehensive R Archive Network (CRAN), 2018. Available at https://cran.cnr.berkeley.edu/ (Last accessed January 23, 2019 ).

29.

Sharko

, Grinstein

, and Marx

. Vectorized Radviz and its application to multiple cluster datasets. IEEE Trans Vis Comput Graph, 2008; 14:1444–1451.

30.

Stanford Medicine ME/CFS Initiative. R Code for a High-Dimensional Cell-Phenotype Visualization Tool. http://med.stanford.edu/chronicfatiguesyndrome/research/Rcode-Visualization-Tool.html

31.

Suni

, Ghanekar

, Houck

, et al. CD4(+)CD8(dim) T lymphocytes exhibit enhanced cytokine expression, proliferation and cytotoxic activity in response to HCMV and HIV-1 antigens. Eur J Immunol, 2001; 31:2512–2520.

32.

Tibshirani

, Saunders

, Rosset

, et al. Sparsity and smoothness via the fused lasso. J R Stat Soc Ser B (Stat Method), 2005; 67:91–108.

33.

Van Long

, and Linsen

. Visualizing high density clusters in multidimensional data using optimized star coordinates. Comput Stat, 2011; 26:655–678.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.69 MB

0.12 MB

0.19 MB

0.10 MB

0.22 MB

0.13 MB

0.02 MB

0.00 MB

Penalized Supervised Star Plots: Example Application in Influenza-Specific CD4+ T Cells

Abstract

Get full access to this article

References

Supplementary Material