Towards a PBMC "virogram assay" for precision medicine: Concordance between ex vivo and in vivo viral infection transcriptomes

BACKGROUND Understanding individual patient host-response to viruses is key to designing optimal personalized therapy. Unsurprisingly, in vivo human experimentation to understand individualized dynamic response of the transcriptome to viruses are rarely studied because of the obvious limitations stemming from ethical considerations of the clinical risk. OBJECTIVE In this rhinovirus study, we first hypothesized that ex vivo human cells response to virus can serve as a proxy for otherwise controversial in vivo human experimentation. We further hypothesized that the N-of-1-pathways framework, previously validated in cancer, can be effective in understanding the more subtle individual transcriptomic response to viral infection. METHOD N-of-1-pathways computes a significance score for a given list of gene sets at the patient level, using merely the 'omics profiles of two paired samples as input. We extracted the peripheral blood mononuclear cells (PBMC) of four human subjects, aliquoted in two paired samples, one subjected to ex vivo rhinovirus infection. Their dysregulated genes and pathways were then compared to those of 9 human subjects prior and after intranasal inoculation in vivo with rhinovirus. Additionally, we developed the Similarity Venn Diagram, a novel visualization method that goes beyond conventional overlap to show the similarity between two sets of qualitative measures. RESULTS We evaluated the individual N-of-1-pathways results using two established cohort-based methods: GSEA and enrichment of differentially expressed genes. Similarity Venn Diagrams and individual patient ROC curves illustrate and quantify that the in vivo dysregulation is recapitulated ex vivo both at the gene and pathway level (p-values⩽0.004). CONCLUSION We established the first evidence that an interpretable dynamic transcriptome metric, conducted as an ex vivo assays for a single subject, has the potential to predict individualized response to infectious disease without the clinical risks otherwise associated to in vivo challenges. These results serve as a foundational work for personalized "virograms".

[1]  中尾 光輝,et al.  KEGG(Kyoto Encyclopedia of Genes and Genomes)〔和文〕 (特集 ゲノム医学の現在と未来--基礎と臨床) -- (データベース) , 2000 .

[2]  P. Brown,et al.  Stunned Silence: Gene Expression Programs in Human Cells Infected with Monkeypox or Vaccinia Virus , 2011, PloS one.

[3]  Yong Huang,et al.  Single Sample Expression-Anchored Mechanisms Predict Survival in Head and Neck Cancer , 2012, PLoS Comput. Biol..

[4]  Y. Lussier,et al.  Curation-free biomodules mechanisms in prostate cancer predict recurrent disease , 2013, BMC Medical Genomics.

[5]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[6]  E. Holmes,et al.  Intra- and Interhost Evolutionary Dynamics of Equine Influenza Virus , 2010, Journal of Virology.

[7]  Ben S. Wittner,et al.  Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1 , 2009, Nature.

[8]  Yves A. Lussier,et al.  Biomolecular Systems of Disease Buried Across Multiple GWAS Unveiled by Information Theory and Ontology , 2010, Summit on translational bioinformatics.

[9]  J. Jiang,et al.  Multi-word complex concept retrieval via lexical semantic similarity , 1999, Proceedings 1999 International Conference on Information Intelligence and Systems (Cat. No.PR00446).

[10]  Susumu Goto,et al.  Data, information, knowledge and principle: back to metabolism in KEGG , 2013, Nucleic Acids Res..

[11]  L. Carin,et al.  Gene expression signatures diagnose influenza and other symptomatic respiratory viral infections in humans. , 2009, Cell host & microbe.

[12]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[13]  Yves A. Lussier,et al.  Complex-disease networks of trait-associated single-nucleotide polymorphisms (SNPs) unveiled by information theory , 2012, J. Am. Medical Informatics Assoc..

[14]  Ian T. Foster,et al.  ‘N-of-1-pathways’ unveils personal deregulated mechanisms from a single pair of RNA-Seq samples: towards precision medicine , 2014, J. Am. Medical Informatics Assoc..

[15]  J. Banchereau,et al.  Whole Blood Gene Expression Profiles to Assess Pathogenesis and Disease Severity in Infants with Respiratory Syncytial Virus Infection , 2013, PLoS medicine.

[16]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[17]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[18]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[19]  E. Holmes,et al.  Evolution of an Eurasian Avian-like Influenza Virus in Naïve and Vaccinated Pigs , 2012, PLoS pathogens.

[20]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[21]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Jean YH Yang,et al.  Bioconductor: open software development for computational biology and bioinformatics , 2004, Genome Biology.

[23]  G. Cheng,et al.  The host type I interferon response to viral and bacterial infections , 2005, Cell Research.

[24]  Yutein Chung,et al.  Macrophage activation state determines the response to rhinovirus infection in a mouse model of allergic asthma , 2014, Respiratory Research.

[25]  Christopher B Wilson,et al.  Regulation of interferon-gamma during innate and adaptive immune responses. , 2007, Advances in immunology.

[26]  Yong Huang,et al.  Towards Mechanism Classifiers: Expression-anchored Gene Ontology Signature Predicts Clinical Outcome in Lung Adenocarcinoma Patients , 2012, AMIA.

[27]  Yves A. Lussier,et al.  GO-Module: functional synthesis and improved interpretation of Gene Ontology patterns , 2011, Bioinform..

[28]  Yves A Lussier,et al.  Concordance of deregulated mechanisms unveiled in underpowered experiments: PTBP1 knockdown case study , 2014, BMC Medical Genomics.

[29]  Giorgio Valle,et al.  The Gene Ontology in 2010: extensions and refinements , 2009, Nucleic Acids Res..

[30]  Yves A. Lussier,et al.  Translating Mendelian and complex inheritance of Alzheimer's disease genes for predicting unique personal genome variants , 2012, J. Am. Medical Informatics Assoc..

[31]  Carol Friedman,et al.  Information theory applied to the sparse gene ontology annotation network to predict novel gene function , 2007, ISMB/ECCB.

[32]  Yves A. Lussier,et al.  Interpreting Personal Transcriptomes: Personalized Mechanism-Scale Profiling of RNA-seq Data , 2012, Pacific Symposium on Biocomputing.

[33]  Antoine M. van Oijen,et al.  Real-time single-molecule observation of rolling-circle DNA replication , 2009, Nucleic acids research.