N-of-1-pathways MixEnrich: advancing precision medicine via single-subject analysis in discovering dynamic changes of transcriptomes

• Published in: BMC medical genomics Vol. 10; no. Suppl 1; pp. 27 - 16
• Main Authors: Li, Qike, Schissler, A. Grant, Gardeux, Vincent, Achour, Ikbel, Kenost, Colleen, Berghout, Joanne, Li, Haiquan, Zhang, Hao Helen, Lussier, Yves A.
• Format: Journal Article
• Language: English
• Published: London BioMed Central 24.05.2017; BioMed Central Ltd; BMC
• Subjects: Bioinformatics; Biomedical and Life Sciences; Biomedicine; Cancer; Competition; Consortia; Data processing; Datasets; Gene Expression; Gene Expression Profiling - methods; Genetic regulation; Genetic research; Genomes; Genomics; Head & neck cancer; Head and neck; Head and neck carcinoma; Head and Neck Neoplasms - genetics; Head and neck squamous cell carcinomas (HNSCCs); Health aspects; Human Genetics; Humans; Microarrays; Mixture Model; N-of-1-pathways; Neoplasms, Squamous Cell - genetics; Ontology; Patients; Precision Medicine; Properties; Ribonucleic acid; RNA; RNA-Seq; ROC Curve; Single-Subject Analysis; Squamous cell carcinoma; Statistical analysis; Transcription factors; Head and neck squamous cell carcinomas (HNSCCs); Single-Subject Analysis; RNA-Seq; Mixture Model; N-of-1-pathways; Precision Medicine
• ISSN: 1755-8794; 1755-8794
• DOI: 10.1186/s12920-017-0263-4

Background Transcriptome analytic tools are commonly used across patient cohorts to develop drugs and predict clinical outcomes. However, as precision medicine pursues more accurate and individualized treatment decisions, these methods are not designed to address single-patient transcriptome analyses. We previously developed and validated the N-of-1-pathways framework using two methods, Wilcoxon and Mahalanobis Distance (MD), for personal transcriptome analysis derived from a pair of samples of a single patient. Although, both methods uncover concordantly dysregulated pathways, they are not designed to detect dysregulated pathways with up- and down-regulated genes (bidirectional dysregulation) that are ubiquitous in biological systems. Results We developed N-of-1-pathways MixEnrich, a mixture model followed by a gene set enrichment test, to uncover bidirectional and concordantly dysregulated pathways one patient at a time. We assess its accuracy in a comprehensive simulation study and in a RNA-Seq data analysis of head and neck squamous cell carcinomas (HNSCCs). In presence of bidirectionally dysregulated genes in the pathway or in presence of high background noise, MixEnrich substantially outperforms previous single-subject transcriptome analysis methods, both in the simulation study and the HNSCCs data analysis (ROC Curves; higher true positive rates; lower false positive rates). Bidirectional and concordant dysregulated pathways uncovered by MixEnrich in each patient largely overlapped with the quasi-gold standard compared to other single-subject and cohort-based transcriptome analyses. Conclusion The greater performance of MixEnrich presents an advantage over previous methods to meet the promise of providing accurate personal transcriptome analysis to support precision medicine at point of care.