title: Optimal deconvolution of transcriptional profiling data using quadratic programming with application to complex clinical blood samples.
abstract: Large-scale molecular profiling technologies have assisted the identification of disease biomarkers and facilitated the basic understanding of cellular processes. However, samples collected from human subjects in clinical trials possess a level of complexity, arising from multiple cell types, that can obfuscate the analysis of data derived from them. Failure to identify, quantify, and incorporate sources of heterogeneity into an analysis can have widespread and detrimental effects on subsequent statistical studies.We describe an approach that builds upon a linear latent variable model, in which expression levels from mixed cell populations are modeled as the weighted average of expression from different cell types. We solve these equations using quadratic programming, which efficiently identifies the globally optimal solution while preserving non-negativity of the fraction of the cells. We applied our method to various existing platforms to estimate proportions of different pure cell or tissue types and gene expression profilings of distinct phenotypes, with a focus on complex samples collected in clinical trials. We tested our methods on several well controlled benchmark data sets with known mixing fractions of pure cell or tissue types and mRNA expression profiling data from samples collected in a clinical trial. Accurate agreement between predicted and actual mixing fractions was observed. In addition, our method was able to predict mixing fractions for more than ten species of circulating cells and to provide accurate estimates for relatively rare cell types (<10% total population). Furthermore, accurate changes in leukocyte trafficking associated with Fingolomid (FTY720) treatment were identified that were consistent with previous results generated by both cell counts and flow cytometry. These data suggest that our method can solve one of the open questions regarding the analysis of complex transcriptional data: namely, how to identify the optimal mixing fractions in a given experiment.
authors: Biomarker Development, Novartis Institutes for BioMedical Research, Cambridge, Massachusetts, United States of America. ting.gong@novartis.com
authors: Ting
authors: T
authors: Gong
authors: Nicole
authors: N
authors: Hartmann
authors: Isaac S
authors: IS
authors: Kohane
authors: Volker
authors: V
authors: Brinkmann
authors: Frank
authors: F
authors: Staedtler
authors: Martin
authors: M
authors: Letzkus
authors: Sandrine
authors: S
authors: Bongiovanni
authors: Joseph D
authors: JD
authors: Szustakowski
mesh_terms: D000328:Adult; D000465:Algorithms; D001769:Blood / Q000378:metabolism*; D002460:Cell Line; D019295:Computational Biology / Q000379:methods*; D005260:Female; D020869:Gene Expression Profiling / Q000379:methods*; D006801:Humans; D016018:Least-Squares Analysis; D016014:Linear Models; D012333:RNA, Messenger / Q000097:blood / Q000235:genetics / Q000378:metabolism; D014158:Transcription, Genetic / Q000235:genetics*
chemical_list: D012333:RNA, Messenger