Detecting heterogeneity in single-cell RNA-Seq data by non-negative matrix factorization

PeerJ

Volumn 2017, Issue 1, 2017, Pages

  Zhu, Xun ^a,b   Ching, Travers ^a,b   Pan, Xinghua ^c   Weissman, Sherman M ^c   Garmire, Lana ^a  

^a UNIVERSITY OF HAWAII CANCER CENTER   (United States)

^b UNIVERSITY OF HAWAII   (United States)

^c YALE UNIVERSITY SCHOOL OF MEDICINE   (United States)

Author keywords

Clustering; Feature gene; Heterogeneity; Modularity; Non negative matrix factorization; RNA Seq; Single cell; Single cell sequencing; Single cell; Subpopulation

Indexed keywords

ARTICLE; BONE MARROW; CELL DIFFERENTIATION; DOWNSTREAM PROCESSING; EPIGENETIC REPRESSION; GENE EXPRESSION; GENE ONTOLOGY; GENETIC HETEROGENEITY; GLIOBLASTOMA; HEMATOPOIETIC STEM CELL; HUMAN; MOLECULAR GENETICS; PLURIPOTENT STEM CELL; POLYMERASE CHAIN REACTION; RNA SEQUENCE; WORKFLOW;

EID: 85013204492     PISSN: None     EISSN: 21678359     Source Type: Journal    
DOI: 10.7717/peerj.2888     Document Type: Article

References (66)

1. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, Robertson CL, Serova N, Davis S, Soboleva A. 2013. NCBI GEO: archive for functional genomics data sets-update. Nucleic Acids Research 41:D991-D995 DOI 10.1093/nar/gks1193.
2. Biase F, Cao X, Zhong S. 2014. Cell fate inclination within 2-cell and 4-cell mouse em-bryos revealed by single-cell RNA sequencing. Genome Research 24(11):1787-1796 DOI 10.1101/gr.177725.114.
3. Blake-Palmer KG, Su Y, Smith AN, Karet FE. 2007. Molecular cloning and character-ization of a novel form of the human vacuolar HC-ATPase e-subunit: an essential proton pump component. Gene 393:94-100 DOI 10.1016/j.gene.2007.01.020.
4. Brennecke P, Anders S, Kim JK, Kołodziejczyk AA, Zhang X, Proserpio V, Baying B, Benes V, Teichmann SA, Marioni JC. 2013. Accounting for technical noise in single-cell RNA-seq experiments. Nature Methods 10(11):1093-1095 DOI 10.1038/nmeth.2645.
5. Brunet J-P, Tamayo P, Golub TR, Mesirov JP. 2004. Metagenes and molecular pattern discovery using matrix factorization. Proceedings of the National Academy of Sciences of the United States of America 101:4164-4169 DOI 10.1073/pnas.0308531101.
6. Brunskill EW, Park J-S, Chung E, Chen F, Magella B, Potter SS. 2014. Single cell dissection of early kidney development: multilineage priming. Development 141:3093-3101 DOI 10.1242/dev.110601.
7. Buettner F, Natarajan KN, Casale FP, Proserpio V, Scialdone A, Theis FJ, Teichmann SA, Marioni JC, Stegle O. 2015. Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells. Nature Biotechnology 33:155-160 DOI 10.1038/nbt.3102.
8. Bushati N, Smith J, Briscoe J, Watkins C. 2011. An intuitive graphical visualization technique for the interrogation of transcriptome data. Nucleic Acids Research 39:7380-7389 DOI 10.1093/nar/gkr462.
9. Ding B, Zheng L, Zhu Y, Li N, Jia H, Ai R, Wildberg A, Wang W. 2015. Normal-ization and noise reduction for single cell RNA-seq experiments. Bioinformatics 31(13):2225-2227 DOI 10.1093/bioinformatics/btv122.
10. Edgar R, Domrachev M, Lash AE. 2002. Gene expression omnibus: NCBI gene expres-sion and hybridization array data repository. Nucleic Acids Research 30:207-210 DOI 10.1093/nar/30.1.207.
11. Fan J-BJJ-BJ, Salathia N, Liu R, Kaeser GE, Yung YC, Herman JL, Kaper F, Fan JBJJ-BJ, Zhang K, Chun J, Kharchenko PV. 2016. Characterizing transcriptional heterogeneity through pathway and gene set overdispersion analysis. Nature Methods 13(3):241-244 DOI 10.1038/nmeth.3734.
12. Gao Y, Church G. 2005. Improving molecular cancer class discovery through sparse non-negative matrix factorization. Bioinformatics 21:3970-3975 DOI 10.1093/bioinformatics/bti653.
13. Gaujoux R, Seoighe C. 2010. A flexible R package for nonnegative matrix factorization. BMC Bioinformatics 11:367 DOI 10.1186/1471-2105-11-367.
14. George L, Irmler M, Mitra A, Thimraj TA, Beckers J, Upadhyay S, Schulz H, Leikauf GD, Ganguly K. 2016. Mouse lung developmental gene expression profiling study identifies novel candidate genes for pulmonary function. B58. Big and bigger (data): omics and biomarkers of COPD and other chronic lung diseases, American Journal of Respiratory and Critical Care Medicine 2016 193:A4078-A4078.
15. Hicks SC, Teng M, Irizarry RA. 2015. On the widespread and critical impact of systematic bias and batch effects in single-cell RNA-Seq data. bioRxiv Preprint. DOI 10.1101/025528.
16. Huang DW, Sherman BT, Lempicki RA. 2008. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature Protocols 4:44-57 DOI 10.1038/nprot.2008.211.
17. Huang DW, Sherman BT, Lempicki RA. 2009. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Research 37:1-13 DOI 10.1093/nar/gkn923.
18. Huang M, Zhong Z, Lv M, Shu J, Tian Q, Chen J. 2016. Comprehensive analysis of dif-ferentially expressed profiles of lncRNAs and circRNAs with associated co-expression and ceRNA networks in bladder carcinoma. Oncotarget 7(30):47186-47200 DOI 10.18632/oncotarget.9706.
19. Jia Z, Zhang X, Guan N, Bo X, Barnes MR, Luo Z. 2015. Gene ranking of RNA-seq data via discriminant non-negative matrix factorization. PLOS ONE 10:e0137782 DOI 10.1371/journal.pone.0137782.
20. Jiao Y, Widschwendter M, TeschendorffAE. 2014. A systems-level integrative frame-work for genome-wide DNA methylation and gene expression data identifies differential gene expression modules under epigenetic control. Bioinformatics 30(16):2360-2366 DOI 10.1093/bioinformatics/btu316.
21. Junker JP, Noël ES, Guryev V, Peterson KA, Shah G, Huisken J, McMahon AP, Berezikov E, Bakkers J, Van oudenaarden A. 2014. Genome-wide RNA tomography in the zebrafish embryo. Cell 159:662-675 DOI 10.1016/j.cell.2014.09.038.
22. Karolchik D, Barber GP, Casper J, Clawson H, Cline MS, Diekhans M, Dreszer TR, Fujita PA, Guruvadoo L, Haeussler M. 2014. The UCSC genome browser database: 2014 update. Nucleic Acids Research 42:D764-D770 DOI 10.1093/nar/gkt1168.
23. Kharchenko PV, Silberstein L, Scadden DT. 2014. Bayesian approach to single-cell dif-ferential expression analysis. Nature Methods 11:740-742 DOI 10.1038/nmeth.2967.
24. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. 2013. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biology 14(4):Article 1 DOI 10.1186/gb-2013-14-4-r36.
25. Kopp S, Slumstrup L, Corydon TJ, Sahana J, Aleshcheva G, Islam T, Magnusson NE, Wehland M, Bauer J, Infanger M. 2016. Identifications of novel mecha-nisms in breast cancer cells involving duct-like multicellular spheroid formation after exposure to the Random Positioning Machine. Scientific Reports 6:26887 DOI 10.1038/srep26887.
26. Kumar RM, Cahan P, Shalek AK, Satija R, Daleykeyser AJ, Li H, Zhang J, Pardee K, Gennert D, Trombetta JJ, Ferrante TC, Regev A, Daley GQ, Collins JJ. 2014. Deconstructing transcriptional heterogeneity in pluripotent stem cells. Nature 516:56-61 DOI 10.1038/nature13920.
27. Leinonen R, Sugawara H, Shumway M, Collaboration on behalf of the INSD. 2011. The sequence read archive. Nucleic Acids Research 39:D19-D21 DOI 10.1093/nar/gkq1019.
28. Li L, Akashi K. 2003. Unraveling the molecular components and genetic blueprints of stem cells. Biotechniques 35:1233-1239.
29. Liao Y, Smyth GK, Shi W. 2014. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30:923-930 DOI 10.1093/bioinformatics/btt656.
30. Love MI, Huber W, Anders S. 2014. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. bioRxiv Preprint DOI 10.1101/002832.
31. Ma K, Wang H, Li X, Li T, Su G, Yang P, Wu J. 2015b. Long noncoding RNA MALAT1 associates with the malignant status and poor prognosis in glioma. Tumor Biology 36:3355-3359 DOI 10.1007/s13277-014-2969-7.
32. Ma X-Y, Wang J-H, Wang J-L, Ma CX, Wang X-C, Liu F-S. 2015a. Malat1 as an evolutionarily conserved lncRNA, plays a positive role in regulating proliferation and maintaining undifferentiated status of early-stage hematopoietic cells. BMC Genomics 16:676 DOI 10.1186/s12864-015-1881-x.
33. Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, Tirosh I, Bialas AR, Kamitaki N, Martersteck EM. 2015. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell 161:1202-1214 DOI 10.1016/j.cell.2015.05.002.
34. McDavid A, Finak G, Chattopadyay PK, Dominguez M, Lamoreaux L, Ma SS, Roederer M, Gottardo R. 2013. Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments. Bioinformatics 29:461-467 DOI 10.1093/bioinformatics/bts714.
35. Monga V, Mihçak MK. 2007. Robust and secure image hashing via non-negative matrix factorizations. IEEE Transactions on Information Forensics and Security 2:376-390 DOI 10.1109/TIFS.2007.902670.
36. Moore HC, Johnston M, Nicol SM, Bourdon J-C, Thompson AM, Hutvagner G, Fuller-Pace FV. 2011. An evolutionarily conserved, alternatively spliced, intron in the p68/DDX5 DEAD-box RNA helicase gene encodes a novel miRNA. RNA 17:555-562 DOI 10.1261/rna.2591611.
37. Pan X. 2014. Single cell analysis: from technology to biology and medicine. Single Cell Biology 3:Article 1 DOI 10.4172/2168-9431.1000106.
38. Panebianco F, Mazzanti C, Tomei S, Aretini P, Franceschi S, Lessi F, Di Coscio G, Bevilacqua G, Marchetti I. 2015. The combination of four molecular markers improves thyroid cancer cytologic diagnosis and patient management. BMC Cancer 15:1 DOI 10.1186/1471-2407-15-1.
39. Park I, Qian D, Kiel M, Becker MW, Pihalja M, Weissman IL, Morrison SJ, Clarke MF. 2003. Bmi-1 is required for maintenance of adult self-renewing haematopoietic stem cells. Nature 423:302-305 DOI 10.1038/nature01587.
40. Patel AP, Tirosh I, Trombetta JJ, Shalek AK, Gillespie SM, Wakimoto H, Cahill DP, Nahed BV, Curry WT, Martuza RL, Louis DN, Rozenblatt-Rosen O, Suvà ML, Regev A, Bernstein BE. 2014. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 344:1396-1401 DOI 10.1126/science.1254257.
41. Pearson-White S, Deacon D, Crittenden R, Brady G, Iscove N, Quesenberry PJ. 1995. The ski/sno protooncogene family in hematopoietic development. Blood 86:2146-2155.
42. Phillips RL, Ernst RE, Brunk B, Ivanova N, Mahan MA, Deanehan JK, Moore KA, Overton GC, Lemischka IR. 2000. The genetic program of hematopoietic stem cells. Science 288:1635-1640 DOI 10.1126/science.288.5471.1635.
43. Poirion OB, Zhu X, Ching T, Garmire L. 2016. Single-cell transcriptomics bioinformat-ics and computational challenges. Frontiers in Genetics 7:Article 163.
44. Qi Q, Zhao Y, Li M, Simon R. 2009. Non-negative matrix factorization of gene ex-pression profiles: a plug-in for BRB-ArrayTools. Bioinformatics 25:545-547 DOI 10.1093/bioinformatics/btp009.
45. Rajapakse M, Tan J, Rajapakse J. 2004. Color channel encoding with NMF for face recognition. In: 2004 international conference on image processing (ICIP 2004). Piscataway: IEEE, 2007-2010.
46. Rand WM. 1971. Objective criteria for the evaluation of clustering methods. Journal of the American Statistical Association 66:846-850 DOI 10.1080/01621459.1971.10482356.
47. Reichardt J, Bornholdt S. 2006. Statistical mechanics of community detection. Physical Review E 74(1):Article 016110 DOI 10.1103/PhysRevE.74.016110.
48. Risso D, Ngai J, Speed TP, Dudoit S. 2014. Normalization of RNA-seq data using factor analysis of control genes or samples. Nature Biotechnology 32:896-902 DOI 10.1038/nbt.2931.
49. Robinson MD, McCarthy DJ, Smyth GK. 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139-140 DOI 10.1093/bioinformatics/btp616.
50. Schlitzer A, Sivakamasundari V, Chen J, Sumatoh HRB, Schreuder J, Lum J, Malleret B, Zhang S, Larbi A, Zolezzi F, Renia L, Poidinger M, Naik S, Newell EW, Robson P, F Ginhoux. 2015. Identification of cDC1-and cDC2-committed DC progenitors reveals early lineage priming at the common DC progenitor stage in the bone marrow. Nature Immunology 16:718-728 DOI 10.1038/ni.3200.
51. Smaragdis P. 2004. Non-negative matrix factor deconvolution; extraction of multiple sound sources from monophonic inputs. In: Independent component analysis and blind signal separation. Berlin Heidelberg: Springer, 494-499.
52. Summers AR, Fischer MA, Stengel KR, Zhao Y, Kaiser JF, Wells CE, Hunt A, Bhaskara S, Luzwick JW, Sampathi S. 2013. HDAC3 is essential for DNA replication in hematopoietic progenitor cells. The Journal of Clinical Investigation 123:3112-3123 DOI 10.1172/JCI60806.
53. Tamayo P, Scanfeld D, Ebert BL, Gillette MA, Roberts CWM, Mesirov JP. 2007. Metagene projection for cross-platform, cross-species characterization of global transcriptional states. Proceedings of the National Academy of Sciences of the United States of America 104:5959-5964 DOI 10.1073/pnas.0701068104.
54. Tian H, Biehs B, Warming S, Leong KG, Rangell L, Klein OD, De Sauvage FJ. 2011. A reserve stem cell population in small intestine renders Lgr5-positive cells dispensable. Nature 478:255-259 DOI 10.1038/nature10408.
55. Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, Lennon NJ, Livak KJ, Mikkelsen TS, Rinn JL. 2014. Pseudo-temporal ordering of individual cells reveals dynamics and regulators of cell fate decisions. Nature Biotechnology 32:381-386 DOI 10.1038/nbt.2859.
56. Treutlein B, Brownfield DG, Wu AR, NeffNF, Mantalas GL, Espinoza FH, Desai TJ, Krasnow MA, Quake SR. 2014. Reconstructing lineage hierarchies of the distal lung epithelium using single-cell RNA-seq. Nature 509:371-375 DOI 10.1038/nature13173.
57. Tung P-Y, Blischak JD, Hsiao C, Knowles DA, Burnett JE, Pritchard JK, Gilad Y. 2016. Batch effects and the effective design of single-cell gene expression studies. bioRxiv Preprint. DOI 10.1101/025528.
58. Usoskin D, Furlan A, Islam S, Abdo H, Lönnerberg P, Lou D, Hjerling-Leffler J, Haeggström J, Kharchenko O, Kharchenko PV. 2014. Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. Nature Neuroscience 18:145-153 DOI 10.1038/nn.3881.
59. Van der Maaten L. 2013. Barnes-hut-sne. ArXiv preprint. arXiv:1301.3342.
60. Van der Maaten L, Hinton G. 2008. Visualizing data using t-SNE. Journal of Machine Learning Research 9:2579-2605.
61. Yang Z, Zhang H, Yuan Z, Oja E. 2011. Kullback-Leibler divergence for nonnegative matrix factorization. In: Artificial neural networks and machine learning-ICANN 2011. Berlin Heidelberg: Springer, 250-257.
62. Ye X, Liu J, Wu F-X. 2012. Dynamic miRNA-TF-mRNA circuits in mouse lung develop-ment. In: 2012 IEEE 6th international conference on systems biology (ISB). Piscataway: IEEE, 256-262.
63. Yuan Z, Oja E. 2005. Projective nonnegative matrix factorization for image compression and feature extraction. In: Image analysis. Berlin Heidelberg: Springer, 333-342.
64. Zeisel A, Muñoz-Manchado AB, Codeluppi S, Lönnerberg P, La Manno G, Juréus A, Marques S, Munguba H, He L, Betsholtz C. 2015. Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. Science 347:1138-1142 DOI 10.1126/science.aaa1934.
65. Zhang P, Cao H-Y, Bai L-L, Li W-N, Wang Y, Chen S-Y, Zhang L, Yang L-H, Xu HT, Wang E-H. 2015. The high expression of TC1 (C8orf4) was correlated with the expression of β-catenin and cyclin D1 and the progression of squamous cell carcino-mas of the tongue. Tumor Biology 36:7061-7067 DOI 10.1007/s13277-015-3423-1.
66. Zhu P, Wang Y, Du Y, He L, Huang G, Zhang G, Yan X, Fan Z. 2015. C8orf4 negatively regulates self-renewal of liver cancer stem cells via suppression of NOTCH2 signalling. Nature Communications 6:7122 DOI 10.1038/ncomms8122.