Inferring activity changes of transcription factors by binding association with sorted expression profiles

BackgroundThe identification of transcription factors (TFs) associated with a biological process is fundamental to understanding its regulatory mechanisms. From microarray data, however, the activity changes of TFs often cannot be directly observed due to their relatively low expression levels, post-transcriptional modifications, and other complications. Several approaches have been proposed to infer TF activity changes from microarray data. In some models, a linear relationship between gene expression and TF-gene binding strength is assumed. In some other models, the target genes of a TF are first determined by a significance cutoff to binding affinity scores, and then expression differentiation is checked between the target and other genes.ResultsWe propose a novel method, referred to as BASE (binding association with sorted expression), to infer TF activity changes from microarray expression profiles with the help of binding affinity data. It searches the maximum association between bind affinity profile of a TF and expression change profile along the direction of sorted differentiation. The method does not make hard target gene selection, rather, the significances of TF activity changes are evaluated by permutation tests of binding association at the end. To show the effectiveness of this method, we apply it to three typical examples using different kinds of binding affinity data, namely, ChIP-chip data, motif discovery data, and positional weighted matrix scanning data, respectively. The implications obtained from all three examples are consistent with established biological results. Moreover, the inferences suggest new and biological meaningful hypotheses for further investigation.ConclusionThe proposed method makes transcription inference from profiles of expression and binding affinity. The same machinery can be used to deal with various kinds of binding affinity data. The method does not require a linear assumption, and has the desirable property of scale-invariance with respect to TF-specific binding affinity. This method is easy to implement and can be routinely applied for transcriptional inferences in microarray studies.

[1]  Henry Horng-Shing Lu,et al.  Statistical methods for identifying yeast cell cycle transcription factors. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[2]  B. Vogelstein,et al.  p53 mutations in human cancers. , 1991, Science.

[3]  H. Bussemaker,et al.  Regulatory element detection using correlation with expression , 2001, Nature Genetics.

[4]  E. Wingender,et al.  MATCH: A tool for searching transcription factor binding sites in DNA sequences. , 2003, Nucleic acids research.

[5]  B. Rayet,et al.  Aberrant rel/nfkb genes and activity in human cancer , 1999, Oncogene.

[6]  Huanying Ge,et al.  Inference of transcription modification in long-live yeast strains from their expression profiles , 2007, BMC Genomics.

[7]  John D. Storey,et al.  Statistical significance for genomewide studies , 2003, Proceedings of the National Academy of Sciences of the United States of America.

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

[9]  S. Hanes,et al.  The ESS1 prolyl isomerase and its suppressor BYE1 interact with RNA pol II to inhibit transcription elongation in Saccharomyces cerevisiae. , 2003, Genetics.

[10]  Michael A. Beer,et al.  Predicting Gene Expression from Sequence , 2004, Cell.

[11]  David G Johnson,et al.  Distinct and Overlapping Roles for E2F Family Members in Transcription, Proliferation and Apoptosis. , 2006, Current molecular medicine.

[12]  Inder M Verma,et al.  NF-kappaB regulation in the immune system. , 2002, Nature reviews. Immunology.

[13]  Patrick J. Killion,et al.  Genetic reconstruction of a functional transcriptional regulatory network , 2007, Nature Genetics.

[14]  Kathleen Marchal,et al.  PKA and Sch9 control a molecular switch important for the proper adaptation to nutrient availability , 2004, Molecular microbiology.

[15]  B. Bauvois,et al.  Interferons Inhibit Tumor Necrosis Factor-α-mediated Matrix Metalloproteinase-9 Activation via Interferon Regulatory Factor-1 Binding Competition with NF-κB* , 2002, The Journal of Biological Chemistry.

[16]  Nicola J. Rinaldi,et al.  Transcriptional regulatory code of a eukaryotic genome , 2004, Nature.

[17]  Charles Boone,et al.  Identifying transcription factor functions and targets by phenotypic activation , 2006, Proceedings of the National Academy of Sciences.

[18]  A. Boulesteix,et al.  Predicting transcription factor activities from combined analysis of microarray and ChIP data: a partial least squares approach , 2005, Theoretical Biology and Medical Modelling.

[19]  Alexander E. Kel,et al.  MATCHTM: a tool for searching transcription factor binding sites in DNA sequences , 2003, Nucleic Acids Res..

[20]  P. Venugopal,et al.  Alterations in IRF1/IRF2 expression in acute myelogenous leukemia , 2001, American journal of hematology.

[21]  G. Church,et al.  Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae. , 2000, Journal of molecular biology.

[22]  Jun S. Liu,et al.  Integrating regulatory motif discovery and genome-wide expression analysis , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Yudong D. He,et al.  Functional Discovery via a Compendium of Expression Profiles , 2000, Cell.

[24]  Ivo Pedruzzi,et al.  TOR and PKA signaling pathways converge on the protein kinase Rim15 to control entry into G0. , 2003, Molecular cell.

[25]  T. Hughes,et al.  Exploration of Essential Gene Functions via Titratable Promoter Alleles , 2004, Cell.

[26]  M. Spector,et al.  Hir1p and Hir2p function as transcriptional corepressors to regulate histone gene transcription in the Saccharomyces cerevisiae cell cycle , 1997, Molecular and cellular biology.

[27]  M. Tyers,et al.  Feedback‐regulated degradation of the transcriptional activator Met4 is triggered by the SCFMet30 complex , 2000, The EMBO journal.

[28]  K. Helin,et al.  The E2F transcription factors: key regulators of cell proliferation. , 2000, Biochimica et biophysica acta.

[29]  D. Botstein,et al.  Genomic expression programs in the response of yeast cells to environmental changes. , 2000, Molecular biology of the cell.

[30]  B Hamilton,et al.  Nuclear localization of the C2H2 zinc finger protein Msn2p is regulated by stress and protein kinase A activity. , 1998, Genes & development.

[31]  J. Vaupel,et al.  Chronological aging‐independent replicative life span regulation by Msn2/Msn4 and Sod2 in Saccharomyces cerevisiae , 2004, FEBS letters.

[32]  Claudio De Virgilio,et al.  Short Article TOR and PKA Signaling Pathways Converge on the Protein Kinase Rim15 to Control Entry into G 0 , 2003 .

[33]  J. Nevins,et al.  The Rb/E2F pathway and cancer. , 2001, Human molecular genetics.

[34]  Feng Gao,et al.  Defining transcriptional networks through integrative modeling of mRNA expression and transcription factor binding data , 2004, BMC Bioinformatics.

[35]  J. Deckert,et al.  Multiple elements and auto-repression regulate Rox1, a repressor of hypoxic genes in Saccharomyces cerevisiae. , 1995, Genetics.

[36]  J. Jauniaux,et al.  Gzf3p, a fourth GATA factor involved in nitrogen‐regulated transcription in Saccharomyces cerevisiae , 1997, Molecular microbiology.

[37]  Robert Clarke,et al.  Interferon regulatory factor-1 (IRF-1) exhibits tumor suppressor activities in breast cancer associated with caspase activation and induction of apoptosis. , 2005, Carcinogenesis.

[38]  S. Reed,et al.  Regulation of Transcription by Ubiquitination without Proteolysis Cdc34/SCFMet30-Mediated Inactivation of the Transcription Factor Met4 , 2000, Cell.

[39]  Chiara Sabatti,et al.  Network component analysis: Reconstruction of regulatory signals in biological systems , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[40]  V. Longo,et al.  Regulation of Longevity and Stress Resistance by Sch9 in Yeast , 2001, Science.

[41]  B. Bauvois,et al.  Interferons inhibit tumor necrosis factor-alpha-mediated matrix metalloproteinase-9 activation via interferon regulatory factor-1 binding competition with NF-kappa B. , 2002, The Journal of biological chemistry.

[42]  K. Macleod pRb and E2f-1 in mouse development and tumorigenesis. , 1999, Current opinion in genetics & development.

[43]  Frederick R. Cross,et al.  High Functional Overlap Between MluI Cell-Cycle Box Binding Factor and Swi4/6 Cell-Cycle Box Binding Factor in the G1/S Transcriptional Program in Saccharomyces cerevisiae , 2005, Genetics.

[44]  T. Gilmore,et al.  Introduction to NF-kappaB: players, pathways, perspectives. , 2006, Oncogene.

[45]  Shin Namiki,et al.  IRF‐1 mediates upregulation of LMP7 by IFN‐γ and concerted expression of immunosubunits of the proteasome , 2005, FEBS letters.

[46]  I. Fingerman,et al.  p53-Mediated Transcriptional Activation From Test Tube to Cell , 2004, Cell.

[47]  D. Lane,et al.  The p53 tumour suppressor gene , 1998, The British journal of surgery.

[48]  T. Heinemeyer,et al.  Expanding the TRANSFAC database towards an expert system of regulatory molecular mechanisms , 1999, Nucleic Acids Res..

[49]  Gene H Golub,et al.  Integrative analysis of genome-scale data by using pseudoinverse projection predicts novel correlation between DNA replication and RNA transcription. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[50]  I. Pedruzzi,et al.  Saccharomyces cerevisiae Ras/cAMP pathway controls post‐diauxic shift element‐dependent transcription through the zinc finger protein Gis1 , 2000, The EMBO journal.