B-cell Ligand Processing Pathways Detected by Large-scale Comparative Analysis

The initiation of B-cell ligand recognition is a critical step for the generation of an immune response against foreign bodies. We sought to identify the biochemical pathways involved in the B-cell ligand recognition cascade and sets of ligands that trigger similar immunological responses. We utilized several comparative approaches to analyze the gene coexpression networks generated from a set of microarray experiments spanning 33 different ligands. First, we compared the degree distributions of the generated networks. Second, we utilized a pairwise network alignment algorithm, BiNA, to align the networks based on the hubs in the networks. Third, we aligned the networks based on a set of KEGG pathways. We summarized our results by constructing a consensus hierarchy of pathways that are involved in B cell ligand recognition. The resulting pathways were further validated through literature for their common physiological responses. Collectively, the results based on our comparative analyses of degree distributions, alignment of hubs, and alignment based on KEGG pathways provide a basis for molecular characterization of the immune response states of B-cells and demonstrate the power of comparative approaches (e.g., gene coexpression network alignment algorithms) in elucidating biochemical pathways involved in complex signaling events in cells.

[1]  B. Efron.  The jackknife, the bootstrap, and other resampling plans , 1987 .

[2]  Bonnie Berger,et al.  IsoRankN: spectral methods for global alignment of multiple protein networks , 2009, Bioinform..

[3]  Shizuo Akira,et al.  Regulation of innate immune responses by autophagy-related proteins , 2010 .

[4]  Hans-Peter Kriegel,et al.  Shortest-path kernels on graphs , 2005, Fifth IEEE International Conference on Data Mining (ICDM'05).

[5]  Roderic D. M. Page,et al.  TreeView: an application to display phylogenetic trees on personal computers , 1996, Comput. Appl. Biosci..

[6]  Vasant Honavar,et al.  Aligning Biomolecular Networks Using Modular Graph Kernels , 2009, WABI.

[7]  Luigi Aloe,et al.  Nerve Growth Factor as a Signaling Molecule for Nerve Cells and also for the Neuroendocrine-Immune Systems , 2009, Reviews in the neurosciences.

[8]  Wojciech Szpankowski,et al.  Pairwise Alignment of Protein Interaction Networks , 2006, J. Comput. Biol..

[9]  Pablo Villoslada,et al.  A computational analysis of protein-protein interaction networks in neurodegenerative diseases , 2008, BMC Systems Biology.

[10]  R. Scheuermann,et al.  Tyrosine kinase activation in the decision between growth, differentiation, and death responses initiated from the B cell antigen receptor. , 2000, Advances in immunology.

[11]  Facundo D Batista,et al.  Early events in B cell activation. , 2010, Annual review of immunology.

[12]  Richard H. Scheuermann,et al.  Analysis of the Major Patterns of B Cell Gene Expression Changes in Response to Short-Term Stimulation with 33 Single Ligands1 , 2004, The Journal of Immunology.

[13]  T. Ideker,et al.  Modeling cellular machinery through biological network comparison , 2006, Nature Biotechnology.

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

[15]  Peng Yang,et al.  Components of the antigen processing and presentation pathway revealed by gene expression microarray analysis following B cell antigen receptor (BCR) stimulation , 2006, BMC Bioinformatics.

[16]  J. Felsenstein.  CONFIDENCE LIMITS ON PHYLOGENIES: AN APPROACH USING THE BOOTSTRAP , 1985, Evolution; international journal of organic evolution.

[17]  Jason Wittenberg,et al.  Clarify: Software for Interpreting and Presenting Statistical Results , 2003 .

[18]  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.

[19]  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.

[20]  Laurin A. J. Mueller,et al.  Integrative Network Biology: Graph Prototyping for Co-Expression Cancer Networks , 2011, PloS one.

[21]  Vincent Frouin,et al.  A Myc-regulated transcriptional network controls B-cell fate in response to BCR triggering , 2009, BMC Genomics.

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

[23]  Anita B. Roberts,et al.  REGULATION OF IMMUNE RESPONSES BY TGF-β* , 1998 .

[24]  Antal F. Novak,et al.  networks Græmlin : General and robust alignment of multiple large interaction data , 2006 .

[25]  Jeet Sukumaran,et al.  A justification for reporting the majority-rule consensus tree in Bayesian phylogenetics. , 2008, Systematic biology.

[26]  G. Altay,et al.  Structural influence of gene networks on their inference: analysis of C3NET. , 2011 .

[27]  Roded Sharan,et al.  Efficient Algorithms for Detecting Signaling Pathways in Protein Interaction Networks , 2006, J. Comput. Biol..

[28]  G. Marsaglia,et al.  Evaluating Kolmogorov's distribution , 2003 .

[29]  M. Bogo,et al.  Adenosine A(2A) receptor agonist (CGS-21680) prevents endotoxin-induced effects on nucleotidase activities in mouse lymphocytes. , 2011, European journal of pharmacology.

[30]  Yoshihiro Yamanishi,et al.  KEGG for linking genomes to life and the environment , 2007, Nucleic Acids Res..

[31]  Vasant Honavar,et al.  Detection of gene orthology from gene co-expression and protein interaction networks , 2010, BMC Bioinformatics.

[32]  Roded Sharan,et al.  NetworkBLAST: comparative analysis of protein networks , 2008 .

[33]  Shankar Subramaniam,et al.  Identification of signaling components required for the prediction of cytokine release in RAW 264.7 macrophages , 2006, Genome Biology.

[34]  J. D. Dal Porto,et al.  B cell antigen receptor signaling 101. , 2004, Molecular immunology.

[35]  Frank Reimann,et al.  TAC3 and TACR3 mutations in familial hypogonadotropic hypogonadism reveal a key role for Neurokinin B in the central control of reproduction , 2009, Nature Genetics.

[36]  John Quackenbush.  Microarray data normalization and transformation , 2002, Nature Genetics.

[37]  D. Chaplin.  Overview of the immune response. , 2003, The Journal of allergy and clinical immunology.

[38]  A. DeFranco,et al.  Molecular aspects of B-lymphocyte activation. , 1987, Annual review of cell biology.

[39]  Ashok Reddy Dinasarapu,et al.  Signaling gateway molecule pages - a data model perspective , 2011, Bioinform..

[40]  Vasant Honavar,et al.  Detection of gene orthology from gene co-expression and protein interaction networks , 2010 .

[41]  Roded Sharan,et al.  PathBLAST: a tool for alignment of protein interaction networks , 2004, Nucleic Acids Res..

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

[43]  Chitta Baral,et al.  Pairwise Alignment of Interaction Networks by Fast Identification of Maximal Conserved Patterns , 2008, Pacific Symposium on Biocomputing.