- Open Access
Sequence homology in eukaryotes (SHOE): interactive visual tool for promoter analysis
BMC Genomicsvolume 19, Article number: 715 (2018)
Microarray and DNA-sequencing based technologies continue to produce enormous amounts of data on gene expression. This data has great potential to illuminate our understanding of biology and medicine, but the data alone is of limited value without computational tools to allow human investigators to visualize and interpret it in the context of their problem of interest.
We created a web server called SHOE that provides an interactive, visual presentation of the available evidence of transcriptional regulation and gene co-expression to facilitate its exploration and interpretation. SHOE predicts the likely transcription factor binding sites in orthologous promoters of humans, mice, and rats using the combined information of 1) transcription factor binding preferences (position-specific scoring matrix (PSSM) libraries such as Transfac32, Jaspar, HOCOMOCO, ChIP-seq, SELEX, PBM, and iPS-reprogramming factor), 2) evolutionary conservation of putative binding sites in orthologous promoters, and 3) co-expression tendencies of gene pairs based on 1,714 normal human cells selected from the Gene Expression Omnibus Database.
SHOE enables users to explore potential interactions between transcription factors and target genes via multiple data views, discover transcription factor binding motifs on top of gene co-expression, and visualize genes as a network of gene and transcription factors on its native gadget GeneViz, the CellDesigner pathway analyzer, and the Reactome database to search the pathways involved. As we demonstrate here when using the CREB1 and Nf-κB datasets, SHOE can reliably identify experimentally verified interactions and predict plausible novel ones, yielding new biological insights into the gene regulatory mechanisms involved. SHOE comes with a manual describing how to run it on a local PC or via the Garuda platform (www.garuda-alliance.org), where it joins other popular gadgets such as the CellDesigner pathway analyzer and the Reactome database, as part of analysis workflows to meet the growing needs of molecular biologists and medical researchers. SHOE is available from the following URL http://ec2-54-150-223-65.ap-northeast-1.compute.amazonaws.com
A video demonstration of SHOE can be found here: https://www.youtube.com/watch?v=qARinNb9NtE
The analysis of gene regulatory regions is a centrally important problem in biology. Many experimental [1,2,3] and computational methods [4, 5] have been developed to address this problem. Despite these efforts and considerable progress, the analysis of eukaryotic gene regulatory regions remains difficult. One fundamental reason for this is that the binding sites of transcription factors (TFs) are only partially determined by their intrinsic sequence specificity; they are also strongly affected by factors including post-translational modification, interactions with other proteins, and the epigenetic state of the genome. Moreover, TF binding events are not necessarily all functional. Thus, effective promoter analysis is not just a simple matter of reporting motif matches or scores, but rather it demands careful consideration of multiple sources of supporting evidence such as the evolutionary conservation of potential binding sites and the coherence of the set of genes with promoters containing them regarding co-expression and cellular pathways.
Ideally a computer program could automatically and reliably combine all available evidence, but this has not yet been achieved. While many excellent motif discovery tools, such as CONSENSUS , the Gibbs sampler [7, 8], CRMD , and MEME , and alignment-based programs, such as rVista , ConSite [12, 13], Footer , and GPminer , have been developed, the current analytical demands require more than a simple output of putative TF binding sites. One example of an attempt to build an analysis workflow is PAINT  where regulatory analysis is represented in gene networks utilizing public software tools with original analysis.
SHOE takes a different approach; it tries to provide human biology experts with an interactive, visual presentation of the available evidence to facilitate the exploration and interpretation of transcription regulation analysis results. We demonstrate the SHOE analytical workflow by analyzing the CREB1 and Nf-κB datasets.
Input and output
As input, SHOE takes a list of human, mouse, or rat Refseq gene ids. First, SHOE assesses if the received gene has orthologs in the other two species. If so, SHOE extracts the promoter and undertakes the computational steps depicted in Fig. 1 and all steps are descibed in the following Algorithm part. As output, SHOE returns a list of genes with motif hits that are common among the three species with the SHOE criteria and gene co-expression. Users can view this output in several ways, as shown in Fig. 2. A motif hit table and gene network window are provided with options that can help users to select the part of the data of their interest, sort by scores, change cutoffs, and save data in text and graphical formats. In order to obtain information on the transcription factor binding, ArrayExpess  ChIP-seq data on Liver has been added to the results table. By clicking on the motif in the table User can judge how close the peak location to the motif identified. For further addition of tissue-specific data upon users request the workflow is provided in SHOE Manual “Adding ChIP-seq Array Express tissue-specific data” page. The user can change the expression interaction threshold, and these changes will be reflected in the results table and saved. All tables can be reset. SHOE returns a gene co-expression network. Pareto front selection allows the user to stipulate how strongly evolutionary conservation should be weighted vis-à-vis PSSM scores. Emphasizing the former might be suitable for those who are looking for novel motifs in well-conserved alignment blocks, while the dominance of the latter will be more interesting for clinicians working on human data.
A flowchart of the execution of the web-tool SHOE algorithm is shown in Fig. 1.
An ortholog gene list of 6,669 human, mouse, and rat genes was obtained from the DBTSS database  of experimentally verified transcriptional start sites on the level of mRNA expression.
Local pairwise and multiple local alignments
Each human promoter was aligned with his/her orthologous partners from the list of orthologs (mouse and rat) by using the SEARCH local alignment program ; the execution of this program was repeated for promoters defined with length 1000, 2000 and 5000 nt, respectively. In the cases where the same human promoter region was aligned to mouse and rat promoters with a similarity higher than 50%, the three genomes’ respective regions were extracted, realigned with ClustalW , and finally stored in a MySQL database to avoid recomputing.
Multiple alignment (MA) score
To evaluate the degree to which an observed region in an alignment of the three species should contribute to our belief that it is part of a conserved region, we adopted an estimate of the likelihood ratio of observing the region in an alignment of orthologous promoters versus that of observing the region in an alignment of unrelated promoters.
To obtain this estimate, first, we randomly selected 1,000 orthologous triplets (human, mouse, and rat genes) and repeatedly aligned their orthologous promoters of length 5,000 nt. After 835 orthologous three-species alignments (238,800 bp total length) were obtained with ClustalW, the frequency of each alignment column was observed and stored in an orthologous alignment frequency table, i.e., an ortho table. In the same manner, 1,000 non-orthologous triples were randomly selected, and their promoters of length 5,000 nt were repeatedly aligned. After 1,260 non-orthologous three-species alignments (239,600 bp total length) were obtained with ClustalW, the frequency of each possible alignment column was observed and stored as a random alignment frequency in a table.
We refer to the obtained alignment column frequencies as the “ortho” alignment and “random” alignment frequencies, respectively. Using those column frequency tables, we define the multiple alignment (MA) scores using the following formula:
where c is the probability of the observed pattern in each column (calculated using pattern frequencies from orthoalignment and random alignment tables), and m is the length of the alignment region, which for our application is equal to the motif length m.
Position-specific scoring matrix (PSSM) score
After multiple alignment scores have been calculated for the human-mouse-rat regions of length m to evaluate the degree of similarity between the identified conserved sequence and known transcription factor binding sites, SHOE uses its motif-scoring module, which calculates similarity scores with matrices in public databases such as Transfac32 (492 matrices) , Jaspar (138) , HOCOMOCO (1,986) , ChIP-seq (381) , SELEX (791) , PBM (302) , and iPS factor matrices (9) , which was created by the authors and is described in the SHOE manual page.
The PSSM-related score is calculated using the following formula:
where pseudocount = 1, and m is the motif length.
To represent the scores as positive values, the PSSM score is subtracted from 10.
Pareto-optimal front for motif selection
Since SHOE identified putative TF binding sites via two sources of information (multiple alignment scores and PSSM scores), we decided to apply the Pareto-optimal front method to optimize this multi-objective solution . The problem of exploring solutions under multiple objectives has traditionally been tackled in engineering, utilizing a so-called “desirability function”, whose value is 1 when the response takes values considered valid by the analyst and 0 otherwise. The Pareto number evaluates the degree to which both the PSSM score and the MA score are favorable. The application of the Pareto front in SHOE is shown in Additional file 1: Figure S1.
The Pareto number of (xi,yi), where xi is the MA score and yi is the PSSM score, is calculated by
where MA and PSSM are the sets of scores.
Motif enrichment score (MES)
To compare the frequencies of motifs in two user datasets, we acquired the motif enrichment score (MES) , which is applied to the standard deviation of a binomial distribution and is a measure of the evolutionary conservation of motifs. Here, we use that metric in a different context as a sample size-dependent measure of the enrichment of matches to a motif in a gene set. In our application, the MES of a motif m is given by
where K is the number of genes with a match to motif m in the gene set of interest, N is the total number of genes in the gene set, and p0 is the frequency of the same motif in a background set of genes.
Computing gene co-regulation on Eigen cell
Considering that genes with similar transcriptional profiles might be affected by the same transcriptional mechanism, we calculated the Pearson correlation coefficients for all gene pairs in the datasets based on 112 “Eigen cell” synthetic expression profiles available from the Cell Montage web site (http://cellmontage.cbrc.jp/) .
For the reader’s convenience, here we briefly outline whom those Eigen cells were computed for. First, 5000 human cell profiles from the Gene Expression Omnibus database  were classified into 1,714 normal cell profiles consisting of 89 cell or tissue types. Then, a standard principal component analysis was conducted to reduce the number of dimensions. The 1,714 profiles were transformed by eigenvectors and reduced to the top 112 informative “Eigen cells” using the Keiser-Guttman criteria . The reduced “highly informative gene expression data” or “Eigen cell profiles” as well as the raw cell data are freely available on the Cell Montage website and are adopted in SHOE for calculation of co-expression of genes in the gene set received by SHOE as an input.
To test SHOE, from ChIP Atlas (http://chip-atlas.org) we selected 1000 human promoters (1000 nt length) as the target genes of the CREB1 transcription factor and 1000 human promoters (1000 nt length) as the target genes of the Nf-κB1 transcription factor. With predicted target genes, gene-interaction maps have been visualized in CellDesigner (Figs. 3, 4 and 5). The following two paragraphs give a biological discussion of the results of SHOE analysis of CREB1 and Nf-κB1.
For the 1000 target genes tested with CREB1, 37 had orthologs in both mouse and rat genes. Additional file 2: Figure S2 shows the changes in the number of genes in the dataset depending on the similarity threshold with the consensus sequence. For the visualization, we applied a strict selection criteria, selecting only 13 genes with a similarity PSSM consensus score of ≥0.75 and a Pearson correlation score of |r| ≥ 0.6, as shown in Additional file 2: Figure S2. Those 13 genes predicted as CREB1 targets are visualized in the CellDesigner  map through the original SHOE-CellDesigner plugin, as shown in Fig. 3, and compared to the known interactions in GeneMANIA , as shown in Fig. 4. By comparing our CREB map to GeneMANIA data, we found that several SHOE predicted interactions have already been discovered and experimentally verified.
Verified gene/protein experiment-based interactions between SHOE and GeneMANIA
The blue lines on the CellDesigner map in Fig. 3 and GeneMANIA in Fig. 4 point to the interactions based on co-retention frequencies using radiation hybrid genotyping data . Due to the fact that many of the potential interactions identified by SHOE are known, it is highly likely that some of the remaining predicted interactions are genuine novel ones. Another example is the protein interaction between AP2B1 and TXN2 (Figs. 3 and 4, interaction no. 6, respectively) that was found when screening all pairwise combinations of open reading frames from human ORFeome v 5.1 generated by a binary protein-protein interaction map . As AP2B1 is part of intracellular transports and TXN2 plays a role in the control of mitochondrial reactive oxygen species homeostasis and apoptosis regulation according to the DAVID database , the two proteins could interact in multiple pathways.
Verified experiment-based co-localization between SHOE and GeneMANIA
The co-localization of BTG1 and RPL23 (Figs. 3 and 4, interaction No. 8) was identified when microarray data was combined with extensive genome annotations . The co-localization of BTG1 and TPT1 (Figs. 3 and 4, interaction no. 9, respectively) was found when microarrays were used to monitor splicing, providing experimental evidence for alternative splicing events . As TPT1 is a regulator of cellular growth and proliferation, BTG1 is a regulator of the cell cycle in the DAVID database, and they are co-localized, it is likely that these genes will genuinely interact. The prediction of an interaction between TPT1 and RPL23 (Figs. 3 and 4, interaction no. 11, respectively) was found by constructing a functional protein network by extending curated pathways with non-curated sources .
Verified predictions between SHOE and GeneMANIA
Some of the genes identified such as mitogen-activated protein kinase kinase 1 (MAP2K1), ribosomal protein L23 (RPL23), and ribosomal protein S29 (RPS29) are part of the conserved structures or pathways in the DAVID database. RPL23 and RPS29 (Figs. 3 and 4, interaction no. 7, respectively) are ribosomal proteins that take part in nuclear-transcribed mRNA catabolic processes essential to the survival of cells. MAP2K1 is part of multiple signaling pathways and many pathways in cancer. Identifying MAP2K1 as a potential target of CREB is a step closer to uncovering its regulation, network, and what role it might play in different cancer types. CREB was identified as part of the transcriptional regulation of antioxidant enzymes in brain tissue and as a regulator so that the animal responds properly to stressful conditions . Several of the genes SHOE identified such as MAP2K1 and TXN2 are listed on DAVID as part of a stress response. TXN2, which is an oxidoreductase, was identified as being a potential CREB target gene . This further supports our predictions and strengthens the credibility of SHOE as a predictor of interactions between genes.
Out of the 1000 genes tested for Nf-κB, SHOE identified 17 genes containing conserved Nf-κB binding domains with a similarity score in the PSSM consensus of ≥0.5. Thirteen genes depicted in the CellDesigner Nf-κB map, where genes whose ratio of similarity scores to consensus sequence is higher than 0.6, are visualized in Fig. 5. Ten genes are confirmed by the literature to be target genes of Nf-κB (Fig. 5, red rectangles), making the other three the predicted new candidates (Fig. 5, yellow rectangles).
It is known that transcription of tumor necrosis factor (TNF), which is involved in apoptosis, is induced by Nf-κB in mice . HIF1A included in our map demonstrates that Nf-κB is a direct modulator of HIF1A expression . The literature reported that it is possible for NR4A1 (Nur77) to assist Nf-κB binding to promoters of anti-apoptotic genes . Another apoptotic gene in our network is LGALS1 (galectin-1) whose expression is controlled by the Nf-κB signaling pathway, and the LGALS1 gene is a direct target of the Nf-κB p50 subunit .
Immune system and insulin
Many Nf-κB target genes are part of the immune system such as the ICAM1 gene. qPCR analysis identified ICAM1 as a target of Nf-κB , and Nf-κB activation was shown to induce ICAM1 expression . PSMB9 (LMP2) was found to be coordinately regulated by using an Sp1-GC box and Nf-κB site . Sp-1 itself is predicted by SHOE to target PSMB9. Nf-κB is found to be essential for the basal activity of the mouse TAPBP (tapasin) promoter . Furthermore, in neuroblastoma cells, Nf-κB was found to synergize with IRF1 in enhancing tapasin . The IRF1 regulatory region itself was found to bind Nf-κB . Published experimental results demonstrate that in colonic epithelial cells, using the Nf-κB pathway, IL-1β may induce CXCL10 . INPPL1 (SHIP2) is involved in the regulation of insulin and also plays a role in actin remodeling . It is known that palmitate could induce SHIP2 expression in skeletal muscle via the activation of Nf-κB pathways .
Known and new transcription factors predicted by SHOE
SHOE predicted another three genes that have not yet been associated with Nf-κB in mice and rats but have been predicted by SHOE: TRIP10, STRN, and GABBR1. Also, several transcription factor candidates were found to potentially co-regulate these genes together with Nf-κB.
TRIP10 (CIP4), which is involved in insulin signaling and actin reorganization, has been identified as a potential target gene of Nf-κB. TRIP10 (CIP4) is required for the translocation of GLUT4 to the plasma membrane (http://www.uniprot.org).
SHOE predicted that LEF1 and Sp1 among other transcription factors affect TRIP10 expression. Nf-κB was found to regulate the LEF1 transcription factor and interact with the Sp1 transcription factor . The above result supports the idea that Nf-kB co-regulates with TRIP10, LEF1, and Sp1. Since both TRIP10 and INPPL1 are part of insulin signaling and actin remodeling and are predicted to be regulated by Nf-kB, they might be involved in the insulin-signaling pathway together.
STRN, which is a calmodulin-binding protein, and GABBR1, which is a receptor for gamma-aminobutyric acid, are also predicted by SHOE to have Nf-kB binding sites in their promoters. Interestingly, Nf-kB is reported to be a constitutive transcription factor in glutamatergic neurons . Thus, STRN and GABBR1, as part of an insulin pathway, may be involved in learning and memory.
We report that CREB, Egr-3, E2F, and p53 regulate STRN, and accordingly, Nf-κB was found to regulate CREB . Egr-3  and E2F1,  both separately interact with Nf-kB. P53 and Nf-kB were found to synergistically upregulate multiple genes. Sp1 and E2F bind to the GABBR1 promoter, and these two transcription factors interact with Nf-kB.
Bach2 , IRF family members [58,59,60], SP3 , SRF , ETS proteins , SMAD4, and TGF-β  transcription factors are also known to interact with Nf-kB, which give us a chance to extend the present biological insights into Nf-kB activity. Out of 33 predicted transcription factors, 13 had been confirmed in combinatorial interactions with Nf-κB (Fig. 5, orange), supporting the possibility of Nf-κB co-regulating these genes.
Exploration ChIP-seq peaks with SHOE
In order to provide Users with the insight of how the results of SHOE are overlapping with experimentally predicted transcription factor binding sites, we incorporated ChIP-seq Liver enhancer data into SHOE. With this amendment whichever geneset is analyzed the results can be viewed in overlap with ChiP-seq peaks in Liver, or other tissues by User request. On SHOE Manual page we put the detailed protocol of how to add ChIP-seq data from Array Express to SHOE. Thus all dataset we/users analyzed could be investigated on the presence of overlap with ChiP-seq peaks. As an additional example, we brought the dataset of 154 genes overexpressed in type 2 Diabetic Mouse Liver and present in SHOE orthologous list . SHOE 24 identified genes of 154 as having cross-species conserved regions (Liver_ChIP-seq dataset in “Query List”). Several genes such as FIF4B, PKLR, IGFBP4, ATP5D shown overlap with ChIP-seq peaks (Additional file 3: Figure S3a). In the mTOR pathway dataset ChiP-seq locations have been found for two genes TNF and EIF4B, and in the promoter of TNF gene two peak locations were correctly predicted (Additional file 3: Figure S3b). Other genes in the mTOR pathway dataset did not have peak information. The above means that the analysis by SHOE might bring helpful information for the identification of possible transcription factor candidates within ChIP-seq regions of binding peaks.
To capture more Chip-seq peaks several considerations should be included into the methodology: a) since enhancer regions are less conserved comparing to promoters, they can have weaker scores in multiple sequence alignment, thus being left out, as in the case of ATP5D gene, when both sides of peak are very well conserved but peak itself is missed because of weaker multiple alignment similarity; b) motifs in peaks are more corrupted/(having more mismatches) with consensus sequence of known binding site, as it is observed by SHOE (similarity scores to consensus around 0.5, thus the cutoff threshold should be set significantly lower for ChIP-seq peaks analysis); c) some peaks are distantly allocated from the TSS (to more than 10,000 bp), which means that more sensitive, i.e. partial sequence search directly on those regions might be more appropriate. Despite the above considerations are not included in current SHOE methodology, SHOE is successfully finding ChiP-seq locations in a range of cases thus increasing the confidence of the results that will augment the interest of the user to the software.
Through case studies we demonstrate the utility of SHOE in visualizing and exploring potential regulatory interactions involving TFs. The analysis made by SHOE further emphasizes the evolutionary importance of CREB in stress-regulated responses and the regulation of antioxidant enzymes. Moreover, for the example of Nf-kB, 10 out of 13 predicted target genes had evidence in the literature supporting the predictions. Not only known but also novel potential target genes in the insulin pathway and the nervous system were identified, which may give insight into new pathways and how Nf-kB is involved in their regulation. For 11 out of 33 transcription factors, there were studies showing interactions with Nf-kB, portraying SHOE as a strong predictor of not only target genes but also potential transcription factors that might co-regulate target genes. Analysis on the MAPK pathway using the SHOE method have been discussed in a previous study . By identifying the roles of genes and their relations to other genes, SHOE can be used to help create a potential network of gene interactions.
SHOE is connected to such tools as the CellDesigner pathway editor and analyzer , Percellome database , and Reactome database  via the Garuda platform and uses its native visualizer GeneViz to represent/compare several networks at once. As future work, shortly we plan to add other tools to join the SHOE analytic workflow to meet the growing needs of molecular biologists and medical researchers.
We also are considering adding the option to consider tissue-specific data such as ChIP-seq binding data and epigenetic data such as DNA methylation and histone modification.
Availability and requirements
Project name: SHOE: Interactive visual tool for promoter analysis.
Project home page: http://ec2-54-150-223-65.ap-northeast-1.compute.amazonaws.com
Project demo: https://www.youtube.com/watch?v=qARinNb9NtE
Operating system(s): Windows / macOS / Linux.
Other requirements: Vagrant, VirtualBox.
Restrictions for use by non-academics: None.
Djordjevic M. SELEX experiments: new prospects, applications and data analysis in inferring regulatory pathways. Biomol Eng. 2007;24(2):179–89.
Horak CE, Snyder M. ChIP-chip: a genomic approach for identifying transcription factor binding sites. Methods Enzymol. 2002;350:469–83.
Khambata-Ford S, Liu Y, Gleason C, Dickson M, Altman RB, Batzoglou S, Myers RM. Identification of promoter regions in the human genome by using a retroviral plasmid library-based functional reporter gene assay. Genome Res. 2003;13(7):1765–74.
Narlikar L, Ovcharenko I. Identifying regulatory elements in eukaryotic genomes. Brief Funct Genomic Proteomic. 2009;8(4):215–30.
Tompa M, Li N, Bailey TL, Church GM, De Moor B, Eskin E, Favorov AV, Frith MC, Fu Y, Kent WJ, et al. Assessing computational tools for the discovery of transcription factor binding sites. Nat Biotechnol. 2005;23(1):137–44.
Hertz GZ, Stormo GD. Identifying DNA and protein patterns with statistically significant alignments of multiple sequences. Bioinformatics. 1999;15(7–8):563–77.
Lawrence CE, Altschul SF, Boguski MS, Liu JS, Neuwald AF, Wootton JC. Detecting subtle sequence signals: a Gibbs sampling strategy for multiple alignment. Science. 1993;262(5131):208–14.
Liu X, Brutlag DL, Liu JS. BioProspector: discovering conserved DNA motifs in upstream regulatory regions of co-expressed genes. Pac Symp Biocomput. 2001:127–38.
Li G, Chan TM, Leung KS, Lee KH. A cluster refinement algorithm for motif discovery. IEEE/ACM Trans Comput Biol Bioinform. 2010;7(4):654–68.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–8.
Loots GG, Ovcharenko I, Pachter L, Dubchak I, Rubin EM. rVista for comparative sequence-based discovery of functional transcription factor binding sites. Genome Res. 2002;12(5):832–9.
Corcoran DL, Feingold E, Benos PV. FOOTER: a web tool for finding mammalian DNA regulatory regions using phylogenetic footprinting. Nucleic Acids Res. 2005;33(Web Server issue):W442–6.
Lenhard B, Sandelin A, Mendoza L, Engstrom P, Jareborg N, Wasserman WW. Identification of conserved regulatory elements by comparative genome analysis. J Biol. 2003;2(2):13.
Lee TY, Chang WC, Hsu JB, Chang TH, Shien DM. GPMiner: an integrated system for mining combinatorial cis-regulatory elements in mammalian gene group. BMC Genomics. 2012;13(Suppl 1):S3.
Vadigepalli R, Chakravarthula P, Zak DE, Schwaber JS, Gonye GE. PAINT: a promoter analysis and interaction network generation tool for gene regulatory network identification. OMICS. 2003;7(3):235–52.
Ghosh S, Matsuoka Y, Asai Y, Hsin KY, Kitano H. Software for systems biology: from tools to integrated platforms. Nat Rev Genet. 2011;12(12):821–32.
Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, Holloway E, Kolesnykov N, Lilja P, Lukk M, et al. ArrayExpress--a public database of microarray experiments and gene expression profiles. Nucleic Acids Res. 2007;35(Database issue):D747–50.
Suzuki A, Wakaguri H, Yamashita R, Kawano S, Tsuchihara K, Sugano S, Suzuki Y, Nakai K. DBTSS as an integrative platform for transcriptome, epigenome and genome sequence variation data. Nucleic Acids Res. 2015;43(Database issue):D87–91.
Gotoh O. An improved algorithm for matching biological sequences. J Mol Biol. 1982;162(3):705–8.
Thompson JD, Gibson TJ, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc. Bioinformatics. 2002;Chapter 2:Unit 2.3.
Wingender E, Dietze P, Karas H, Knuppel R. TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Res. 1996;24(1):238–41.
Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004;32(Database issue):D91–4.
Kulakovskiy IV, Vorontsov IE, Yevshin IS, Soboleva AV, Kasianov AS, Ashoor H, Ba-Alawi W, Bajic VB, Medvedeva YA, Kolpakov FA, et al. HOCOMOCO: expansion and enhancement of the collection of transcription factor binding sites models. Nucleic Acids Res. 2016;44(D1):D116–25.
Hume MA, Barrera LA, Gisselbrecht SS, Bulyk ML. UniPROBE, update 2015: new tools and content for the online database of protein-binding microarray data on protein-DNA interactions. Nucleic Acids Res. 2015;43(Database issue):D117–22.
Polouliakh N. Reprogramming resistant genes: in-depth comparison of gene expressions among iPS, ES, and somatic cells. Front Physiol. 2013;4:7.
Reguera C, Sanchez MS, Ortiz MC, Sarabia LA. Pareto-optimal front as a tool to study the behaviour of experimental factors in multi-response analytical procedures. Anal Chim Acta. 2008;624(2):210–22.
Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M. Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature. 2005;434(7031):338–45.
Fujibuchi W, Kiseleva L, Taniguchi T, Harada H, Horton P. CellMontage: similar expression profile search server. Bioinformatics. 2007;23(22):3103–4.
Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, Rudnev D, Lash AE, Fujibuchi W, Edgar R. NCBI GEO: mining millions of expression profiles--database and tools. Nucleic Acids Res. 2005;33(Database issue):D562–6.
Guttman L. Some necessary conditions for common-factor analysis. Psychometrika. 1954;19(2):149–61.
Funahashi A, Jouraku A, Matsuoka Y, Kitano H. Integration of CellDesigner and SABIO-RK. In Silico Biol. 2007;7(2 Suppl):S81–90.
Mostafavi S, Ray D, Warde-Farley D, Grouios C, Morris Q. GeneMANIA: a real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 2008;9(Suppl 1):S4.
Lin A, Wang RT, Ahn S, Park CC, Smith DJ. A genome-wide map of human genetic interactions inferred from radiation hybrid genotypes. Genome Res. 2010;20(8):1122–32.
Rolland T, Tasan M, Charloteaux B, Pevzner SJ, Zhong Q, Sahni N, Yi S, Lemmens I, Fontanillo C, Mosca R, et al. A proteome-scale map of the human interactome network. Cell. 2014;159(5):1212–26.
Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, et al. DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35(Web Server issue):W169–75.
Schadt EE, Edwards SW, GuhaThakurta D, Holder D, Ying L, Svetnik V, Leonardson A, Hart KW, Russell A, Li G, et al. A comprehensive transcript index of the human genome generated using microarrays and computational approaches. Genome Biol. 2004;5(10):R73.
Johnson JM, Castle J, Garrett-Engele P, Kan Z, Loerch PM, Armour CD, Santos R, Schadt EE, Stoughton R, Shoemaker DD. Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science. 2003;302(5653):2141–4.
Wu G, Feng X, Stein L. A human functional protein interaction network and its application to cancer data analysis. Genome Biol. 2010;11(5):R53.
Kikuchi D, Tanimoto K, Nakayama K. CREB is activated by ER stress and modulates the unfolded protein response by regulating the expression of IRE1alpha and PERK. Biochem Biophys Res Commun. 2016;469(2):243–50.
Pellegrini M, Cheng JC, Voutila J, Judelson D, Taylor J, Nelson SF, Sakamoto KM. Expression profile of CREB knockdown in myeloid leukemia cells. BMC Cancer. 2008;8:264.
Collart MA, Baeuerle P, Vassalli P. Regulation of tumor necrosis factor alpha transcription in macrophages: involvement of four kappa B-like motifs and of constitutive and inducible forms of NF-kappa B. Mol Cell Biol. 1990;10(4):1498–506.
van Uden P, Kenneth NS, Rocha S. Regulation of hypoxia-inducible factor-1alpha by NF-kappaB. Biochem J. 2008;412(3):477–84.
de Leseleuc L, Denis F. Inhibition of apoptosis by Nur77 through NF-kappaB activity modulation. Cell Death Differ. 2006;13(2):293–300.
Toscano MA, Campagna L, Molinero LL, Cerliani JP, Croci DO, Ilarregui JM, Fuertes MB, Nojek IM, Fededa JP, Zwirner NW, et al. Nuclear factor (NF)-kappaB controls expression of the immunoregulatory glycan-binding protein galectin-1. Mol Immunol. 2011;48(15–16):1940–9.
Bunting K, Rao S, Hardy K, Woltring D, Denyer GS, Wang J, Gerondakis S, Shannon MF. Genome-wide analysis of gene expression in T cells to identify targets of the NF-kappa B transcription factor c-Rel. J Immunol. 2007;178(11):7097–109.
van de Stolpe A, Caldenhoven E, Stade BG, Koenderman L, Raaijmakers JA, Johnson JP, van der Saag PT. 12-O-tetradecanoylphorbol-13-acetate- and tumor necrosis factor alpha-mediated induction of intercellular adhesion molecule-1 is inhibited by dexamethasone. Functional analysis of the human intercellular adhesion molecular-1 promoter. J Biol Chem. 1994;269(8):6185–92.
Wright KL, White LC, Kelly A, Beck S, Trowsdale J, Ting JP. Coordinate regulation of the human TAP1 and LMP2 genes from shared bidirectional promoter. J Exp Med. 1995;181(4):1459–71.
Herrmann F, Trowsdale J, Huber C, Seliger B. Cloning and functional analyses of the mouse tapasin promoter. Immunogenetics. 2003;55(6):379–88.
Lorenzi S, Forloni M, Cifaldi L, Antonucci C, Citti A, Boldrini R, Pezzullo M, Castellano A, Russo V, van der Bruggen P, et al. IRF1 and NF-kB restore MHC class I-restricted tumor antigen processing and presentation to cytotoxic T cells in aggressive neuroblastoma. PLoS One. 2012;7(10):e46928.
Yeruva S, Ramadori G, Raddatz D. NF-kappaB-dependent synergistic regulation of CXCL10 gene expression by IL-1beta and IFN-gamma in human intestinal epithelial cell lines. Int J Color Dis. 2008;23(3):305–17.
Gorgani-Firuzjaee S, Ahmadi S, Meshkani R. Palmitate induces SHIP2 expression via the ceramide-mediated activation of NF-kappaB, and JNK in skeletal muscle cells. Biochem Biophys Res Commun. 2014;450(1):494–9.
Shindo M, Chaudhary PM. The ectodermal dysplasia receptor represses the Lef-1/beta-catenin-dependent transcription independent of NF-kappaB activation. Biochem Biophys Res Commun. 2004;315(1):73–78.
Kaltschmidt B, Kaltschmidt C. NF-kappaB in the nervous system. Cold Spring Harb Perspect Biol. 2009;1(3):a001271.
McKay LI, Cidlowski JA. CBP (CREB binding protein) integrates NF-kappaB (nuclear factor-kappaB) and glucocorticoid receptor physical interactions and antagonism. Mol Endocrinol. 2000;14(8):1222–34.
Wieland GD, Nehmann N, Muller D, Eibel H, Siebenlist U, Suhnel J, Zipfel PF, Skerka C. Early growth response proteins EGR-4 and EGR-3 interact with immune inflammatory mediators NF-kappaB p50 and p65. J Cell Sci. 2005;118(Pt 14):3203–12.
Kundu M, Guermah M, Roeder RG, Amini S, Khalili K. Interaction between cell cycle regulator, E2F-1, and NF-kappaB mediates repression of HIV-1 gene transcription. J Biol Chem. 1997;272(47):29468–74.
Hunter JE, Butterworth JA, Zhao B, Sellier H, Campbell KJ, Thomas HD, Bacon CM, Cockell SJ, Gewurz BE, Perkins ND. The NF-kappaB subunit c-Rel regulates Bach2 tumour suppressor expression in B-cell lymphoma. Oncogene. 2016;35(26):3476–84.
Chae M, Kim K, Park SM, Jang IS, Seo T, Kim DM, Kim IC, Lee JH, Park J. IRF-2 regulates NF-kappaB activity by modulating the subcellular localization of NF-kappaB. Biochem Biophys Res Commun. 2008;370(3):519–24.
Iwanaszko M, Kimmel M. NF-kappaB and IRF pathways: cross-regulation on target genes promoter level. BMC Genomics. 2015;16:307.
Saura M, Zaragoza C, Bao C, McMillan A, Lowenstein CJ. Interaction of interferon regulatory factor-1 and nuclear factor kappaB during activation of inducible nitric oxide synthase transcription. J Mol Biol. 1999;289(3):459–71.
Carver BJ, Plosa EJ, Stinnett AM, Blackwell TS, Prince LS. Interactions between NF-kappaB and SP3 connect inflammatory signaling with reduced FGF-10 expression. J Biol Chem. 2013;288(21):15318–25.
Franzoso G, Carlson L, Brown K, Daucher MB, Bressler P, Siebenlist U. Activation of the serum response factor by p65/NF-kappaB. EMBO J. 1996;15(13):3403–12.
Bassuk AG, Anandappa RT, Leiden JM. Physical interactions between Ets and NF-kappaB/NFAT proteins play an important role in their cooperative activation of the human immunodeficiency virus enhancer in T cells. J Virol. 1997;71(5):3563–73.
Grau AM, Datta PK, Zi J, Halder SK, Beauchamp RD. Role of Smad proteins in the regulation of NF-kappaB by TGF-beta in colon cancer cells. Cell Signal. 2006;18(7):1041–50.
Zhang F, Xu X, Zhang Y, Zhou B, He Z, Zhai Q. Gene expression profile analysis of type 2 diabetic mouse liver. PLoS One. 2013;8(3):e57766.
Polouliakh N, Natsume T, Harada H, Fujibuchi W, Horton P. Comparative genomic analysis of transcription regulation elements involved in human map kinase G-protein coupling pathway. J Bioinforma Comput Biol. 2006;4(2):469–82.
Kanno J, Aisaki K, Igarashi K, Nakatsu N, Ono A, Kodama Y, Nagao T. “Per cell” normalization method for mRNA measurement by quantitative PCR and microarrays. BMC Genomics. 2006;7:64.
Stein LD. Using the Reactome database. Curr Protoc Bioinformatics. 2004;Chapter 8:Unit 8.7.
We thank Prof. Jun Kanno (Japan Bioassay Research Center, Japan Organization of Occupational Health and Safety), Dr. Ken-ichi Aisaki (National Institute of Health Sciences, Japan), Dr. Satoshi Kitajima (National Institute of Health Sciences, Japan) and Prof. Kenta Nakai (Institute of Medical Sciences, The University of Tokyo) for their insightful discussions and comments.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Pareto front score optimization on SHOE. A) Illustration of the conflict between two scores (MA score and PSSM score) and the concept of dominance; B) Visualization of several groups of scores using Pentachlorophenol response dataset (PCP), PCP_raw (red) denotes random score zone, PCP_ma (purple) denotes zone where MA score is high, PCP_pssm (blue) is the zone where the PSSM score is high, and PCPtrimmedLog (green) is the zone where a trade-off between two scores is taken. C) Example of the analysis with 30 top Pareto fronts with each type of motif shown in different shape and color. (PDF 360 kb)
Figure S2. Trade-off on the number of genes in the dataset basing of motif similarity threshold to the consensus when Pearson correlation thresholds |r| ≥ 0.0 and |r| ≥ 0.6 are applied. (PDF 20 kb)
Figure S3. Visualization of ChiP-seq peaks from ArrayExperss database idenyified in SHOE predictions. A) Demosntrates two genes from overexpressed in mouse liver in Diabet 2 condition; B) Demonstrate TNF genes of mTOR human pathway in which promoter two peaks according ChIP-seq analysis have been identified. (ZIP 499 kb)