Dynamics of dendritic cell maturation are identified through a novel filtering strategy applied to biological time-course microarray replicates

BackgroundDendritic cells (DC) play a central role in primary immune responses and become potent stimulators of the adaptive immune response after undergoing the critical process of maturation. Understanding the dynamics of DC maturation would provide key insights into this important process. Time course microarray experiments can provide unique insights into DC maturation dynamics. Replicate experiments are necessary to address the issues of experimental and biological variability. Statistical methods and averaging are often used to identify significant signals. Here a novel strategy for filtering of replicate time course microarray data, which identifies consistent signals between the replicates, is presented and applied to a DC time course microarray experiment.ResultsThe temporal dynamics of DC maturation were studied by stimulating DC with poly(I:C) and following gene expression at 5 time points from 1 to 24 hours. The novel filtering strategy uses standard statistical and fold change techniques, along with the consistency of replicate temporal profiles, to identify those differentially expressed genes that were consistent in two biological replicate experiments. To address the issue of cluster reproducibility a consensus clustering method, which identifies clusters of genes whose expression varies consistently between replicates, was also developed and applied. Analysis of the resulting clusters revealed many known and novel characteristics of DC maturation, such as the up-regulation of specific immune response pathways. Intriguingly, more genes were down-regulated than up-regulated. Results identify a more comprehensive program of down-regulation, including many genes involved in protein synthesis, metabolism, and housekeeping needed for maintenance of cellular integrity and metabolism.ConclusionsThe new filtering strategy emphasizes the importance of consistent and reproducible results when analyzing microarray data and utilizes consistency between replicate experiments as a criterion in both feature selection and clustering, without averaging or otherwise combining replicate data. Observation of a significant down-regulation program during DC maturation indicates that DC are preparing for cell death and provides a path to better understand the process. This new filtering strategy can be adapted for use in analyzing other large-scale time course data sets with replicates.

[1]  T. Taniguchi,et al.  Distinct and Essential Roles of Transcription Factors IRF-3 and IRF-7 in Response to Viruses for IFN-α/β Gene Induction , 2000 .

[2]  M. Colonna,et al.  Deficit of CD47 Results in a Defect of Marginal Zone Dendritic Cells, Blunted Immune Response to Particulate Antigen and Impairment of Skin Dendritic Cell Migration1 , 2006, The Journal of Immunology.

[3]  R. Steinman,et al.  Dendritic cells and the control of immunity , 1998, Nature.

[4]  Brad T. Sherman,et al.  DAVID: Database for Annotation, Visualization, and Integrated Discovery , 2003, Genome Biology.

[5]  N. Le Meur,et al.  Profiling dendritic cell maturation with dedicated microarrays , 2005, Journal of leukocyte biology.

[6]  J. Dopazo,et al.  Methods and approaches in the analysis of gene expression data. , 2001, Journal of immunological methods.

[7]  J. Lord,et al.  MHC class II‐mediated apoptosis of mature dendritic cells proceeds by activation of the protein kinase C‐δ isoenzyme , 2002 .

[8]  R. Koup,et al.  Toll-Like Receptor Ligands Modulate Dendritic Cells to Augment Cytomegalovirus- and HIV-1-Specific T Cell Responses 1 , 2003, The Journal of Immunology.

[9]  Zachary D. Smith,et al.  Unbiased Reconstruction of a Mammalian Transcriptional Network Mediating Pathogen Responses , 2009 .

[10]  D. Botstein,et al.  Cluster analysis and display of genome-wide expression patterns. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Hong Yan,et al.  Cluster Analysis of Gene Expression Data , 2009, Encyclopedia of Artificial Intelligence.

[12]  M. Rosenzwajg,et al.  Double-stranded RNA stimulation or CD40 ligation of monocyte-derived dendritic cells as models to study their activation and maturation process. , 2004, European cytokine network.

[13]  John Okyere,et al.  How to decide? Different methods of calculating gene expression from short oligonucleotide array data will give different results , 2006, BMC Bioinformatics.

[14]  R. Steinman,et al.  Dendritic cell development and maturation. , 1997, Advances in experimental medicine and biology.

[15]  A. Bowie,et al.  Schlafen-1 Causes a Cell Cycle Arrest by Inhibiting Induction of Cyclin D1* , 2005, Journal of Biological Chemistry.

[16]  Jill P. Mesirov,et al.  Consensus Clustering: A Resampling-Based Method for Class Discovery and Visualization of Gene Expression Microarray Data , 2003, Machine Learning.

[17]  R. Tibshirani,et al.  Significance analysis of microarrays applied to the ionizing radiation response , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Joshua M. Korn,et al.  The plasticity of dendritic cell responses to pathogens and their components. , 2001, Science.

[19]  Baolin Wu,et al.  Differential gene expression detection using penalized linear regression models: the improved SAM statistics , 2005, Bioinform..

[20]  Ka Yee Yeung,et al.  Validating clustering for gene expression data , 2001, Bioinform..

[21]  J. Hiscott,et al.  Regulation of RANTES Chemokine Gene Expression Requires Cooperativity Between NF-κB and IFN-Regulatory Factor Transcription Factors1 , 2000, The Journal of Immunology.

[22]  Henry C. Chang,et al.  Differential Lysosomal Proteolysis in Antigen-Presenting Cells Determines Antigen Fate , 2005, Science.

[23]  Jacquelyn S. Fetrow,et al.  Additional limitations of the clustering validation method figure of merit , 2007, ACM-SE 45.

[24]  Yeung Sam Hung,et al.  Genome-scale cluster analysis of replicated microarrays using shrinkage correlation coefficient , 2008, BMC Bioinformatics.

[25]  Xiaohui Liu,et al.  Consensus clustering and functional interpretation of gene-expression data , 2004, Genome Biology.

[26]  E. Coccia,et al.  Differential responsiveness to IFN-alpha and IFN-beta of human mature DC through modulation of IFNAR expression. , 2006, Journal of leukocyte biology.

[27]  Maria Persico,et al.  Inducible IL-2 production by dendritic cells revealed by global gene expression analysis , 2001, Nature Immunology.

[28]  David Bryant,et al.  DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists , 2007, Nucleic Acids Res..

[29]  J. Chen,et al.  ‘Survival gene’ Bcl-xl potentiates DNA-raised antitumor immunity , 2005, Gene Therapy.

[30]  John D. Storey,et al.  Significance analysis of time course microarray experiments. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[31]  O. Beretta,et al.  Effects of dexamethazone on LPS‐induced activationand migration of mouse dendritic cells revealed by a genome‐wide transcriptional analysis , 2006, European journal of immunology.

[32]  Raffaele Giancarlo,et al.  Computational cluster validation for microarray data analysis: experimental assessment of Clest, Consensus Clustering, Figure of Merit, Gap Statistics and Model Explorer , 2008, BMC Bioinformatics.

[33]  B. Williams,et al.  ISG20, a New Interferon-induced RNase Specific for Single-stranded RNA, Defines an Alternative Antiviral Pathway against RNA Genomic Viruses* , 2003, The Journal of Biological Chemistry.

[34]  E. Coccia,et al.  Differential responsiveness to IFN‐α and IFN‐β of human mature DC through modulation of IFNAR expression , 2006 .

[35]  Roger E Bumgarner,et al.  Clustering gene-expression data with repeated measurements , 2003, Genome Biology.

[36]  M. Knittler,et al.  Molecular Architecture of the TAP-Associated MHC Class I Peptide-Loading Complex1 , 2007, The Journal of Immunology.

[37]  Kevin R. Coombes,et al.  Identifying and Quantifying Sources of Variation in Microarray Data Using High-Density cDNA Membrane Arrays , 2002, J. Comput. Biol..

[38]  David Liu,et al.  DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis , 2007, BMC Bioinformatics.

[39]  Jihoon Kim,et al.  Difference-based clustering of short time-course microarray data with replicates , 2007, BMC Bioinformatics.

[40]  J. Trowsdale,et al.  The transporters associated with antigen presentation. , 1993, Seminars in cell biology.

[41]  I. Mellman,et al.  Activation of Lysosomal Function During Dendritic Cell Maturation , 2003, Science.

[42]  Ian B. Jeffery,et al.  Comparison and evaluation of methods for generating differentially expressed gene lists from microarray data , 2006, BMC Bioinformatics.

[43]  Takuma Tsukahara,et al.  Differential gene induction by type I and type II interferons and their combination. , 2006, Journal of interferon & cytokine research : the official journal of the International Society for Interferon and Cytokine Research.

[44]  G. Church,et al.  Systematic determination of genetic network architecture , 1999, Nature Genetics.

[45]  E. Meyer-Bernstein,et al.  Restoration of circadian behavioural rhythms in a period null Drosophila mutant (per01) by mammalian period homologues mPer1 and mPer2 , 2002, Genes to cells : devoted to molecular & cellular mechanisms.

[46]  Kathleen F. Kerr,et al.  Evaluation of methods for oligonucleotide array data via quantitative real-time PCR , 2006, BMC Bioinformatics.

[47]  S. Akira,et al.  Selective contribution of IFN-α/β signaling to the maturation of dendritic cells induced by double-stranded RNA or viral infection , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[48]  R. F. Woolson.  Wilcoxon Signed-Rank Test , 2008 .

[49]  Michael J. Zilliox,et al.  Gene expression patterns in dendritic cells infected with measles virus compared with other pathogens. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[50]  Laurence Zitvogel,et al.  Antigen presentation and T cell stimulation by dendritic cells. , 2002, Annual review of immunology.

[51]  L. Espert,et al.  Interferon-induced exonuclease ISG20 exhibits an antiviral activity against human immunodeficiency virus type 1. , 2005, The Journal of general virology.

[52]  Aidong Zhang,et al.  Cluster analysis for gene expression data: a survey , 2004, IEEE Transactions on Knowledge and Data Engineering.

[53]  D. Toomre,et al.  Dendritic cell maturation triggers retrograde MHC class II transport from lysosomes to the plasma membrane , 2002, Nature.

[54]  Timothy Wilkes,et al.  Microarray data quality - review of current developments. , 2007, Omics : a journal of integrative biology.

[55]  R. Steinman,et al.  Generation of large numbers of dendritic cells from mouse bone marrow cultures supplemented with granulocyte/macrophage colony-stimulating factor , 1992, The Journal of experimental medicine.

[56]  Gary A. Churchill,et al.  Analysis of Variance for Gene Expression Microarray Data , 2000, J. Comput. Biol..

[57]  S. Dudoit,et al.  A prediction-based resampling method for estimating the number of clusters in a dataset , 2002, Genome Biology.

[58]  I. Mellman,et al.  Developmental Control of Endocytosis in Dendritic Cells by Cdc42 , 2000, Cell.

[59]  R. Steinman,et al.  Transport of peptide-MHC class II complexes in developing dendritic cells. , 2000, Science.

[60]  T. Taniguchi,et al.  Distinct and essential roles of transcription factors IRF-3 and IRF-7 in response to viruses for IFN-alpha/beta gene induction. , 2000, Immunity.

[61]  Charles Kung,et al.  Chemical genomic profiling to identify intracellular targets of a multiplex kinase inhibitor. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[62]  D. Tang,et al.  Functional analysis of a novel KRAB/C2H2 zinc finger protein Mipu1. , 2007, Biochemical and biophysical research communications.

[63]  Genechip ® Expression Analysis Data Analysis Fundamentals Genechip ® Expression Analysis Chapter 4 First-order Data Analysis and Data Quality Assessment..................................... 27 Genechip ® Expression Analysis Affymetrix ® Netaffx Analysis of Promoter Sequences of Regulated Transcripts , 2022 .

[64]  M. Rubio,et al.  CD47 Engagement Inhibits Cytokine Production and Maturation of Human Dendritic Cells1 , 2000, The Journal of Immunology.

[65]  Z. Berneman,et al.  Gene expression signatures in CD34+-progenitor-derived dendritic cells exposed to the chemical contact allergen nickel sulfate. , 2006, Toxicology and applied pharmacology.

[66]  Robert Tibshirani,et al.  Estimating the number of clusters in a data set via the gap statistic , 2000 .

[67]  I. Mellman,et al.  Invariant Chain Controls H2-M Proteolysis in Mouse Splenocytes and Dendritic Cells , 2000, The Journal of experimental medicine.

[68]  Wenxuan Zhong,et al.  A data-driven clustering method for time course gene expression data , 2006, Nucleic acids research.

[69]  Crispin J. Miller,et al.  Cell Culture , 2010, Cell.

[70]  Isabelle Guyon,et al.  A Stability Based Method for Discovering Structure in Clustered Data , 2001, Pacific Symposium on Biocomputing.

[71]  S. Fujieda,et al.  Double-stranded RNA induces production of RANTES and IL-8 by human nasal fibroblasts. , 2006, Clinical immunology.

[72]  S. Hedrick,et al.  Schlafen, a new family of growth regulatory genes that affect thymocyte development. , 1998, Immunity.

[73]  Hongzhe Li,et al.  Clustering of time-course gene expression data using a mixed-effects model with B-splines , 2003, Bioinform..

[74]  Haruhiko Soma,et al.  Transcriptional oscillation of canonical clock genes in mouse peripheral tissues , 2004, BMC Molecular Biology.

[75]  Kathleen A. Kennedy,et al.  Systems biology approaches identify ATF3 as a negative regulator of Toll-like receptor 4 , 2006, Nature.

[76]  S. Akira,et al.  Selective contribution of IFN-alpha/beta signaling to the maturation of dendritic cells induced by double-stranded RNA or viral infection. , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[77]  T. Speed,et al.  Summaries of Affymetrix GeneChip probe level data. , 2003, Nucleic acids research.

[78]  Takayoshi Mamine,et al.  Molecular mechanism of cell-autonomous circadian gene expression of Period2, a crucial regulator of the mammalian circadian clock. , 2005, Molecular biology of the cell.

[79]  S. Kaufmann,et al.  Modulation of T cell development and activation by novel members of the Schlafen (slfn) gene family harbouring an RNA helicase-like motif. , 2004, International immunology.

[80]  Arpad Kelemen,et al.  Differential and trajectory methods for time course gene expression data , 2005, Bioinform..

[81]  Soheil Shams,et al.  Noise Sampling Method: An ANOVA Approach Allowing Robust Selection of Differentially Regulated Genes Measured by DNA Microarrays , 2003, Bioinform..

[82]  Douglas B. Kell,et al.  Computational cluster validation in post-genomic data analysis , 2005, Bioinform..

[83]  G. Pruneri,et al.  Comparative genome‐wide profiling of post‐transplant lymphoproliferative disorders and diffuse large B‐cell lymphomas , 2006, British journal of haematology.

[84]  P. Ricciardi-Castagnoli,et al.  Brief Definitive Report Dendritic Cell Survival and Maturation Are Regulated by Different Signaling Pathways , 2022 .

[85]  P. Moore,et al.  Dynamics of interferon-specific gene expression in peripheral blood of interferon alfa-naïve patients with genotype 1 chronic hepatitis C infection treated with albumin-interferon alfa. , 2006, Hepatology research : the official journal of the Japan Society of Hepatology.

[86]  J. Lord,et al.  MHC class II-mediated apoptosis of mature dendritic cells proceeds by activation of the protein kinase C-delta isoenzyme. , 2002, International immunology.

[87]  Eric J. Brown,et al.  Bidirectional Negative Regulation of Human T and Dendritic Cells by CD47 and Its Cognate Receptor Signal-Regulator Protein-α: Down-Regulation of IL-12 Responsiveness and Inhibition of Dendritic Cell Activation1 , 2001, The Journal of Immunology.

[88]  Kai Wang,et al.  Comparative analysis of microarray normalization procedures: effects on reverse engineering gene networks , 2007, ISMB/ECCB.