Abstract

Analysis of transcriptional regulatory interactions and their comparisons across multiple species are crucial for progress in various fields in biology, from functional genomics to the evolution of signal transduction pathways. However, despite the rapidly growing body of data on regulatory interactions in several eukaryotes, no databases exist to provide curated high-quality information on transcription factor–target gene interactions for multiple species. Here, we address this gap by introducing the TFLink gateway, which uniquely provides experimentally explored and highly accurate information on transcription factor–target gene interactions (∼12 million), nucleotide sequences and genomic locations of transcription factor binding sites (∼9 million) for human and six model organisms: mouse, rat, zebrafish, fruit fly, worm and yeast by integrating 10 resources. TFLink provides user-friendly access to data on transcription factor–target gene interactions, interactive network visualizations and transcription factor binding sites, with cross-links to several other databases. Besides containing accurate information on transcription factors, with a clear labelling of the type/volume of the experiments (small-scale or high-throughput), the source database and the original publications, TFLink also provides a wealth of standardized regulatory data available for download in multiple formats. The database offers easy access to high-quality data for wet-lab researchers, supplies data for gene set enrichment analyses and facilitates systems biology and comparative gene regulation studies.

Database URLhttps://tflink.net/

Introduction

Gene regulation in eukaryotes is complex, with many layers of regulation including transcription factors that are the key protein regulators of gene expression. Phenotypic evolution is closely linked to changes of gene regulation. To explore the evolution of gene regulatory networks, first we need to study and understand the individual regulators and the evolution of their interactions (1). Thus, the accurate identification of transcription factor–target gene interactions and transcription factor binding sites in multiple species is of paramount importance for studying gene regulation.

During the last two decades, our knowledge on transcriptional regulation has substantially expanded, owing to the wide utilization of experimental assays, including small- and large-scale approaches. Small-scale experimental methods can be used to verify an interaction between a specific transcription factor and its potential target gene on a case-by-case basis (see the methods in Table 1). Therefore, these methods can provide information on both the location and the nucleotide sequence of the binding site for a particular transcription factor. Large-scale methods are utilized to identify potential transcription factor–target gene interactions and transcription factor binding sites in a high-throughput manner. In vitro methods can identify protein binding events on entire nucleotide libraries, while in vivo methods are exploited to characterize transcription factor binding on the whole genome.

Table 1.

Small- and large-scale experimental methods to explore transcription factor–target gene interactions or transcription factor binding sites

Evidence typeMethod nameReference
Small-scaleDNase-I footprintingGalas DJ & Schmitz, 1978 (2)
EMSA: electrophoretic mobility shift assayGarner MM & Revzin A, 1981 (3)
SELEX: systematic evolution of ligands by exponential enrichmentPollock R & Treisman R, 1990 (4)
Promoter deletion analyses coupled to reporter assaysde Wet JR, et al., 1987 (5)
Large-scale: in vitroChIP-on-chip: chromatin immunoprecipitation with DNA microarrayRen B, et al., 2000 (6)
HT-SELEX: high-throughput systematic evolution of ligands by exponential enrichmentJolma A, et al., 2010 (7)
MITOMI: Mechanically induced trapping of molecular interactionsRockel S, et al., 2012 (8)
Large-scale: in vivoChIP-seq: chromatin immunoprecipitation coupled with high-throughput sequencingJohnson DS, et al., 2007 (9); Robertson G, et al., 2007 (10)
Evidence typeMethod nameReference
Small-scaleDNase-I footprintingGalas DJ & Schmitz, 1978 (2)
EMSA: electrophoretic mobility shift assayGarner MM & Revzin A, 1981 (3)
SELEX: systematic evolution of ligands by exponential enrichmentPollock R & Treisman R, 1990 (4)
Promoter deletion analyses coupled to reporter assaysde Wet JR, et al., 1987 (5)
Large-scale: in vitroChIP-on-chip: chromatin immunoprecipitation with DNA microarrayRen B, et al., 2000 (6)
HT-SELEX: high-throughput systematic evolution of ligands by exponential enrichmentJolma A, et al., 2010 (7)
MITOMI: Mechanically induced trapping of molecular interactionsRockel S, et al., 2012 (8)
Large-scale: in vivoChIP-seq: chromatin immunoprecipitation coupled with high-throughput sequencingJohnson DS, et al., 2007 (9); Robertson G, et al., 2007 (10)
Table 1.

Small- and large-scale experimental methods to explore transcription factor–target gene interactions or transcription factor binding sites

Evidence typeMethod nameReference
Small-scaleDNase-I footprintingGalas DJ & Schmitz, 1978 (2)
EMSA: electrophoretic mobility shift assayGarner MM & Revzin A, 1981 (3)
SELEX: systematic evolution of ligands by exponential enrichmentPollock R & Treisman R, 1990 (4)
Promoter deletion analyses coupled to reporter assaysde Wet JR, et al., 1987 (5)
Large-scale: in vitroChIP-on-chip: chromatin immunoprecipitation with DNA microarrayRen B, et al., 2000 (6)
HT-SELEX: high-throughput systematic evolution of ligands by exponential enrichmentJolma A, et al., 2010 (7)
MITOMI: Mechanically induced trapping of molecular interactionsRockel S, et al., 2012 (8)
Large-scale: in vivoChIP-seq: chromatin immunoprecipitation coupled with high-throughput sequencingJohnson DS, et al., 2007 (9); Robertson G, et al., 2007 (10)
Evidence typeMethod nameReference
Small-scaleDNase-I footprintingGalas DJ & Schmitz, 1978 (2)
EMSA: electrophoretic mobility shift assayGarner MM & Revzin A, 1981 (3)
SELEX: systematic evolution of ligands by exponential enrichmentPollock R & Treisman R, 1990 (4)
Promoter deletion analyses coupled to reporter assaysde Wet JR, et al., 1987 (5)
Large-scale: in vitroChIP-on-chip: chromatin immunoprecipitation with DNA microarrayRen B, et al., 2000 (6)
HT-SELEX: high-throughput systematic evolution of ligands by exponential enrichmentJolma A, et al., 2010 (7)
MITOMI: Mechanically induced trapping of molecular interactionsRockel S, et al., 2012 (8)
Large-scale: in vivoChIP-seq: chromatin immunoprecipitation coupled with high-throughput sequencingJohnson DS, et al., 2007 (9); Robertson G, et al., 2007 (10)

While the number of such experiments continue to grow, the results of the available studies are stored in scattered resources. Existing transcription factor databases either (i) capture only a few species: e.g. TRRUST (human and mouse) (11), REDfly (fruit fly) (12) or YEASTRACT (yeast) (13) (see more examples in Supplementary Table); (ii) contain a mixture of high-quality small-scale and less accurate large-scale experimental data without clear labelling: e.g. HTRIdb (14), ORegAnno (15) and YEASTRACT or (iii) fail to associate the interactions with the binding sites of transcription factors: e.g. TRRUST, YEASTRACT or JASPAR (16). The lack of a curated, large, multispecies database of transcription factor–target gene interactions and transcription factor binding sites makes comparative studies of transcription factors difficult and labour-intensive. Here, we introduce TFLink, a gateway that integrates experimental transcription factor data from a benchmarked list of existing databases (Table 2). TFLink uniquely provides highly accurate information on transcription factor–target gene interactions, genomic locations and nucleotide sequences of transcription factor binding sites for human and six eukaryotic model organisms: mouse (Mus musculus), rat (Rattus norvegicus), zebrafish (Danio rerio), fruit fly (Drosophila melanogaster), nematode (Caenorhabditis elegans) and yeast (Saccharomyces cerevisiae).

Table 2.

Source databases of TFLink

Source databaseOriginal URLVersionDownloading dateType of dataaNr. of integrated entriesOrganismsb
DoRothEAchttps://saezlab.github.io/dorothea/219/06/2020SS interactions3453Hs
GTRDhttp://gtrd.biouml.org/20.0602/07/2020LS interactions10 685 122Hs, Mm, Rn, Dr, Dm, Ce, Sc
HTRIdbhttp://www.lbbc.ibb.unesp.br/htri/d129/04/2017SS interactions2020Hs
LS interactions47 140
JASPARhttp://jaspar.genereg.net/202022/07/2020SS binding sites3048Hs, Mm, Rn, Dm, Ce
LS binding sites8 567 469
ORegAnnohttp://www.oreganno.org/e324/05/2017SS interactions1979Hs, Mm, Rn, Dm, Ce, Sc
LS interactions160 096
SS binding sites47 304
LS binding sites705 121
REDflyhttp://redfly.ccr.buffalo.edu/6.0.216/06/2020SS interactions683Dm
LS interactions90
SS binding sites2240
LS binding sites27
ReMaphttp://remap.univ-amu.fr/1.216/07/2018LS interactions2 933 177Hs
TREDhttp://rulai.cshl.edu/cgi-bin/TREDf08/06/2018SS interactions8693Hs, Mm
TRRUSThttps://www.grnpedia.org/trrust/230/07/2018SS interactions16 570Hs, Mm
Yeastracthttp://www.yeastract.com/202020/07/2020SS interactions5349Sc
LS interactions188 072
Source databaseOriginal URLVersionDownloading dateType of dataaNr. of integrated entriesOrganismsb
DoRothEAchttps://saezlab.github.io/dorothea/219/06/2020SS interactions3453Hs
GTRDhttp://gtrd.biouml.org/20.0602/07/2020LS interactions10 685 122Hs, Mm, Rn, Dr, Dm, Ce, Sc
HTRIdbhttp://www.lbbc.ibb.unesp.br/htri/d129/04/2017SS interactions2020Hs
LS interactions47 140
JASPARhttp://jaspar.genereg.net/202022/07/2020SS binding sites3048Hs, Mm, Rn, Dm, Ce
LS binding sites8 567 469
ORegAnnohttp://www.oreganno.org/e324/05/2017SS interactions1979Hs, Mm, Rn, Dm, Ce, Sc
LS interactions160 096
SS binding sites47 304
LS binding sites705 121
REDflyhttp://redfly.ccr.buffalo.edu/6.0.216/06/2020SS interactions683Dm
LS interactions90
SS binding sites2240
LS binding sites27
ReMaphttp://remap.univ-amu.fr/1.216/07/2018LS interactions2 933 177Hs
TREDhttp://rulai.cshl.edu/cgi-bin/TREDf08/06/2018SS interactions8693Hs, Mm
TRRUSThttps://www.grnpedia.org/trrust/230/07/2018SS interactions16 570Hs, Mm
Yeastracthttp://www.yeastract.com/202020/07/2020SS interactions5349Sc
LS interactions188 072
a

Abbreviations: SS: small-scale experiments; LS: large-scale experiments.

b

Abbreviations: Hs: Homo sapiens, Mm: Mus musculus, Rn: Rattus norvegicus, Dr: Danio rerio, Dm: D. melanogaster, Ce: Caenorhabditis elegans and Sc: Saccharomyces cerevisiae.

c

TFLink also indicates the original source of the interactions downloaded from DoRothEA (when available), see chapter 2 in Supplementary Notes for the details.

d

HTRIdb is no longer available at its original website, nor in other databases, making the TFLink the only source for this data. We downloaded the content of HTRIdb in 2017, when it was still available.

e

There is no working website, only data can be downloaded in the form of TSV files.

f

TRED is no longer available at its original website; we were able to download its data from the third party via RegNetwork (24) in 2018.

Table 2.

Source databases of TFLink

Source databaseOriginal URLVersionDownloading dateType of dataaNr. of integrated entriesOrganismsb
DoRothEAchttps://saezlab.github.io/dorothea/219/06/2020SS interactions3453Hs
GTRDhttp://gtrd.biouml.org/20.0602/07/2020LS interactions10 685 122Hs, Mm, Rn, Dr, Dm, Ce, Sc
HTRIdbhttp://www.lbbc.ibb.unesp.br/htri/d129/04/2017SS interactions2020Hs
LS interactions47 140
JASPARhttp://jaspar.genereg.net/202022/07/2020SS binding sites3048Hs, Mm, Rn, Dm, Ce
LS binding sites8 567 469
ORegAnnohttp://www.oreganno.org/e324/05/2017SS interactions1979Hs, Mm, Rn, Dm, Ce, Sc
LS interactions160 096
SS binding sites47 304
LS binding sites705 121
REDflyhttp://redfly.ccr.buffalo.edu/6.0.216/06/2020SS interactions683Dm
LS interactions90
SS binding sites2240
LS binding sites27
ReMaphttp://remap.univ-amu.fr/1.216/07/2018LS interactions2 933 177Hs
TREDhttp://rulai.cshl.edu/cgi-bin/TREDf08/06/2018SS interactions8693Hs, Mm
TRRUSThttps://www.grnpedia.org/trrust/230/07/2018SS interactions16 570Hs, Mm
Yeastracthttp://www.yeastract.com/202020/07/2020SS interactions5349Sc
LS interactions188 072
Source databaseOriginal URLVersionDownloading dateType of dataaNr. of integrated entriesOrganismsb
DoRothEAchttps://saezlab.github.io/dorothea/219/06/2020SS interactions3453Hs
GTRDhttp://gtrd.biouml.org/20.0602/07/2020LS interactions10 685 122Hs, Mm, Rn, Dr, Dm, Ce, Sc
HTRIdbhttp://www.lbbc.ibb.unesp.br/htri/d129/04/2017SS interactions2020Hs
LS interactions47 140
JASPARhttp://jaspar.genereg.net/202022/07/2020SS binding sites3048Hs, Mm, Rn, Dm, Ce
LS binding sites8 567 469
ORegAnnohttp://www.oreganno.org/e324/05/2017SS interactions1979Hs, Mm, Rn, Dm, Ce, Sc
LS interactions160 096
SS binding sites47 304
LS binding sites705 121
REDflyhttp://redfly.ccr.buffalo.edu/6.0.216/06/2020SS interactions683Dm
LS interactions90
SS binding sites2240
LS binding sites27
ReMaphttp://remap.univ-amu.fr/1.216/07/2018LS interactions2 933 177Hs
TREDhttp://rulai.cshl.edu/cgi-bin/TREDf08/06/2018SS interactions8693Hs, Mm
TRRUSThttps://www.grnpedia.org/trrust/230/07/2018SS interactions16 570Hs, Mm
Yeastracthttp://www.yeastract.com/202020/07/2020SS interactions5349Sc
LS interactions188 072
a

Abbreviations: SS: small-scale experiments; LS: large-scale experiments.

b

Abbreviations: Hs: Homo sapiens, Mm: Mus musculus, Rn: Rattus norvegicus, Dr: Danio rerio, Dm: D. melanogaster, Ce: Caenorhabditis elegans and Sc: Saccharomyces cerevisiae.

c

TFLink also indicates the original source of the interactions downloaded from DoRothEA (when available), see chapter 2 in Supplementary Notes for the details.

d

HTRIdb is no longer available at its original website, nor in other databases, making the TFLink the only source for this data. We downloaded the content of HTRIdb in 2017, when it was still available.

e

There is no working website, only data can be downloaded in the form of TSV files.

f

TRED is no longer available at its original website; we were able to download its data from the third party via RegNetwork (24) in 2018.

Materials and methods

Data sources

We examined the available transcription factor databases (altogether 66 databases, see Supplementary Table), and integrated data from those that comply with all of the following criteria: (i) include clearly labelled experimental data of transcription factor–target gene interactions and/or transcription factor binding sites, (ii) provide information about the experimental methods or the corresponding publications, (iii) are—at least partially—the primary source of the data and (iv) are freely accessible and can be redistributed. Based on these criteria 10 databases were selected for integration: DoRothEA (17), GTRD (18), HTRIdb (14), JASPAR (16), ORegAnno (15), REDfly (12), ReMap (19), TRED (20), TRRUST (11) and Yeastract (13) (see Table 2 for details). By exploiting these database sources, we integrated accurate, small-scale experimental data and the results of large-scale experiments. We did not synthetize all the available data from the source databases, but only those entries that fulfil the rigorous criteria listed above.

Data processing, curating and database construction

To establish a database of well curated data, only the information matching with the above criteria were integrated from the source databases. We also filtered out incomplete or potentially incorrect data: for example, transcription factor binding sites with a nucleotide content of over 1500 bp (21) were excluded from the TFLink database. Data were converted to uniform formats, and the different gene or protein identifiers and names (applied by the source databases) were mapped to UniProt IDs. Besides the UniProt IDs which were applied as the main identifiers of genes or proteins in the database, we established unique TFLink IDs (starting with ‘TFLink’, followed by SS or LS indicating the small- and large-scale experiments correspondingly, and a unique number) for each transcription factor binding site, which makes the identification of binding sites unambiguous. In the integrated dataset, each interaction is present only once, while the original annotation data from all sources are preserved.

We updated the coordinates of transcription factor binding sites to the newest genome version available in any of the source databases (these were hg38 for human, mm10 for M. musculus, rn6 for R. norvegicus, dm6 for D. melanogaster, ce10 for C. elegans and saccer1 for S. cerevisiae) using the online tool, LiftOver of UCSC Genome Browser (https://genome.ucsc.edu/cgi-bin/hgLiftOver). This step ensures that the user finds reliable data on binding sites, confirmed by multiple experiments and databases (Figure 1). We identified the orthologs of the transcription factors and target genes with the Sherlock platform (22), using the ortholog table downloaded from the OMA orthology resource (23).

Strength of evidence (number of supporting experiments) for each transcription factor binding site, referring to the number of binding sites (for a particular transcription factor) that overlap (in at least one nucleotide length) with the investigated binding site. The original data are included as binding site tables in the download parts of the TFLink database.
Figure 1.

Strength of evidence (number of supporting experiments) for each transcription factor binding site, referring to the number of binding sites (for a particular transcription factor) that overlap (in at least one nucleotide length) with the investigated binding site. The original data are included as binding site tables in the download parts of the TFLink database.

TFLink was built by using static pages hosted on github (https://github.com/), on a dedicated server of the research group. The source code of the website is freely available at https://github.com/korcsmarosgroup/TFlink. We are planning to update the content of TFLink, biennially, to include the latest versions of the source databases, or even include more species and databases, as well as to translate the transcription factor binding site locations using the newest genome versions.

Results

Content

TFLink is a curated, comprehensive, non-redundant database that contains data on 3984 transcription factors, 110 808 target genes, 31 486 regulatory interactions identified by small-scale experiments and 11 826 489 interactions from large-scale experiments. It also contains 9 290 526 binding site locations and 9 325 209 corresponding binding site sequences (for species-specific statistics and explanation, see Table 1 in Supplementary Notes).

Web interface and usage

TFLink is publicly available at https://tflink.net/ without registration. The user can ‘Browse’, search and ‘Download’ all the data, find a short description and summary statistics on the main page and get the answers for ‘frequently asked questions’ at the detailed ‘FAQ’ part.

The ‘Browse’ page

After selecting the organism, the user can browse and search within the dataset. The results can be filtered by gene name, UniProt ID, NCBI Gene ID, function (e.g. ‘transcription factor’, ‘target gene’ or ‘transcription factor and target gene’) and according to evidence type (small- or large-scale experiments). TFLink differentiates between ‘transcription factor’ and ‘transcription factor and target gene’ functions based on whether the transcription factor protein regulating the gene for the particular transcription factor is known (present in the TFLink database) or not. Information on the number of interactions a particular gene or protein is involved in is also provided. After selecting an entry (gene or protein) from the Browsing table, an ‘entry page’ is opened. (Links to example entry pages are provided in chapter 3.2 of the Supplementary Notes and in the FAQ part of the TFLink gateway).

Entry pages

Each entry page contains basic information about the transcription factor protein or target gene: gene name, UniProt ID (linked to the corresponding UniProt protein page (25)), NCBI Gene ID (linked to the corresponding NCBI Gene site (26)), organism (the scientific name of the species), its function (transcription factor, target gene or both), the number of its interactions and its orthologs (species name and UniProt ID)—when there are any (Figure 2A). In case the ortholog is also available in the TFLink database, a link is provided to the related entry page. Binding site nucleotide composition frequency matrices (27) and sequence logos (28) of the transcription factors are also available through the JASPAR website to facilitate the prediction of more binding sites.

Sample content of a TFLink entry page. (A) Basic information about the transcription factor or target gene. (B) Visualization of the interaction network of the transcription factor (red) and its target genes (green). (C) Table containing information about the target genes. (D) Table containing genomic location of the binding sites. (E) Transcription factor binding site sequences in FASTA format, based on small-scale evidences. (F) Downloadable FASTA file containing transcription factor binding site sequences based on large-scale evidences.
Figure 2.

Sample content of a TFLink entry page. (A) Basic information about the transcription factor or target gene. (B) Visualization of the interaction network of the transcription factor (red) and its target genes (green). (C) Table containing information about the target genes. (D) Table containing genomic location of the binding sites. (E) Transcription factor binding site sequences in FASTA format, based on small-scale evidences. (F) Downloadable FASTA file containing transcription factor binding site sequences based on large-scale evidences.

Below the basic information section, the user may visualize three layers of information (if available) about the selected transcription factor: (i) target genes of the transcription factor (Figure 2C) and/or (ii) transcription factors for the target gene and (iii) binding sites of the transcription factor. In the target gene and transcription factor tables, the user finds details on gene names (linked to corresponding TFLink entries), UniProt IDs (linked to corresponding UniProt protein pages), NCBI Gene IDs (linked to corresponding NCBI Gene sites), name of the source database(s), method(s) of detection, cross-links to the original publications (when available) and to the publications of the databases at NCBI PubMed (29), and indications of the evidence type (small- or large-scale experiments). Along with these tables, interactive network visualizations are presented, demonstrating the interactions between the transcription factor(s) and target gene(s) (indicated by green and red colours, respectively; Figure 2B) to facilitate the visual inspection of the interactions.

Besides the TFLink ID, the name of the source database(s), the method(s) of detection, the link to the original publications and the indication to clarify whether the evidence is based on a small- or a large-scale experiment, the binding site table also presents information about the genomic location: genome assembly version, chromosome, the coordinates of the start and end points of transcription factor binding sites and the number of overlapping binding sites for the particular transcription factor (Figure 2D). To make the visual exploration of the genomic context easier, each binding site is linked to its particular genomic location at the UCSC genome browser website (30).

In case there are >100 interactions or binding sites available for a particular entry in the TFLink gateway, we only show the first 100 targets/transcription factors/binding sites in the tables on the website, and make the full information available in the form of downloadable table (and in case of binding sites: GFF3 annotation) files.

The sequences of binding sites based on small-scale evidence are shown below the tables in FASTA format. The header of the sequences contains the TFLink and the UniProt IDs, gene name, genome assembly version, chromosome name and the start and end point coordinates of the binding sites (Figure 2E). Some data downloaded from the JASPAR database refer to binding sequences without exact localization, for example, in cases when random sequences were investigated with SELEX. The binding sequences revealed by large-scale experiments are available from the entry pages as downloadable FASTA files (Figure 2F).

The ‘Download’ page

After selecting the organism of interest, the user can download transcription factor–target gene interaction files in various formats, the transcription factor binding site tables and binding site sequence files combined, and also the data based on small- and large-scale experiments separately. The interaction table is a TSV (tab separated value) file that can be opened by any spreadsheet editors as well as by the Cytoscape software (31) and the igraph R package (32) to visualize the interactions. The interaction MITAB files contain transcription factor–target gene interactions in a standardized MITAB 2.8 format (as defined by the Human Proteome Organization—Proteomics Standards Initiative, HUPO-PSI (33)), which can be used as an input to various software tools, including Cytoscape. Interaction GMT (Gene Matrix Transposed) is a tab delimited file format, which describes gene sets (target genes of a transcription factor) in each row, allowing the user to calculate gene set enrichment (by applying e.g. the GSEA software (34)) or to execute overrepresentation analysis (by applying e.g. the MulEA R package, https://github.com/koralgooll/MulEA) to answer questions like ‘which transcription factors regulate the differentially expressed genes in this experiment?’. The information on transcription factor binding sites is available in TSV format and the sequences of binding sites are provided in FASTA file format. We have also created GFF3 (General Feature Format, version 3) annotation files from the binding site tables, which can be opened by various software tools for NGS analysis, e.g. with the IGV genome viewer (35). For a detailed description of these downloadable file formats, please see chapter 3.3 in the Supplementary Notes. Files over the size limit of 100 Mb are compressed by the very effective 7z algorithm.

Summary of the results

On the TFLink website one can search and browse among the species-specific datasets of transcription factors and target genes. Each of these proteins and genes has its own entry page with information about interactions, binding sites, orthologs and hyperlinks to other TFLink pages and external sources: the UniProt, NCBI Gene, PubMed, JASPAR databases and the UCSC genome browser website. Known binding sites are also indicated here and are available for download as FASTA files as well. Furthermore, the interaction network (containing a maximum of 100 nodes) is also visualized on each entry page. Species-specific datasets are downloadable in TSV, MITAB 2.8, GMT, GFF3 and FASTA formats (see chapter 3.3 in the Supplementary Notes for details), allowing direct analysis in widely-used software tools.

TFLink is a FAIR (36) compliant database, since it provides easy access to transcription factor, target gene and binding site data by offering sufficient metadata, user-friendly browsing and file download in standard formats. TFLink is the only source for the data previously available on the HTRIdb. Moreover, TFLink is the only resource with a working graphical interface for data downloaded from the ORegAnno database. Therefore, our gateway facilitates the distribution of plenty of biological data.

Discussion

Here, we have introduced the TFLink database, which uniquely provides comprehensive and highly accurate information on transcription factor–target gene interactions, nucleotide sequences and genomic locations of transcription factor binding sites for human and six eukaryotic model organisms. To establish the TFLink web resource (https://TFlink.net), we have integrated data from 10, mostly unconnected, resources, translated the names and identifiers, created standard downloadable files, visualized the regulation networks and cross-linked several other webpages containing related or supplementary information about the interactions, binding sites, sources or genomic contents.

Comparison with other databases

We have established a gateway, which is considered gap-filling for several reasons: (i) unlike numerous other transcription factor databases (e.g. Cistrome DB (37), DoRothEA (17), hTFtarget (38), HTRIdb (14), REDfly (12), ReMap (19) and 18 others we have investigated during our search for databases worth being integrated into the TFLink, see Supplementary Table), TFLink is a multispecies gateway. Therefore, TFLink facilitates the investigations into the evolutionary changes in transcription factors, their binding sites and target gene repertoire. (ii) Unlike AnimalTFDB (39) and TRANSFAC (40) all data in TFLink are available for download in easily accessible formats, thereby allowing even less experienced users to analyze large amounts of gene regulation data. (iii) While AnimalTFDB, ScerTF (41), YeTFaSCo (42) and Yeastract only provide the consensus sequences or matrices of transcription factor binding sites, TFLink includes the actual binding sequences, thus it allows the user to identify the binding sites in resequenced genomes, compare the binding sites across different transcription factors and design specific experiments to enhance or suppress gene expression. (iv) While AnimalTFDB and PlantTFDB (43) focus more on the sequence and structural properties of the transcription factors themselves, TFLink provides functional information on transcription factors, besides being cross-linked to the corresponding UniProt protein sites.

Limitations and updating plans

We note that the actual version of TFLink has some important limitations: it covers only seven taxa and it only contains data that were already available from other scattered resources. We also note that in some cases TFLink does not contain all the known target genes of a specific transcription factor due to limited availability of small-scale experiments. In contrast, TFLink may overestimate the set of target genes when considering large-scale data due to the presence of false-positives. Combining TFLink data with tissue-specific gene expression measurements from databases, such as Encode (44), GEO (45) or the SRA (46), can ameliorate this issue. Therefore, besides updating TFLink regularly by integrating new data posted in the original resources, in the next few years we plan to collect transcriptional regulation information not only from transcription factor databases, but also from other resources, like research papers and ExTRI text mining transcription regulation database (47).

General applications

TFLink is a useful resource for wet-lab researchers, since it provides easy access to high-quality transcription factor–target gene interaction and transcription factor binding site data, with cross-links to several other databases. TFLink is also a long-awaited resource for bioinformaticians, as it contains large quantities of standardized, downloadable regulatory data in multiple formats. The provided interaction tables can be used as input data for the Cytoscape software or to the igraph package to perform systems and network biology studies. The GMT files are useful for gene set enrichment and overrepresentation analyses. Binding site tables allow the user to investigate the genomic location of binding sites. Users can apply the GFF3 binding site annotation files in various NGS analyses, for example when investigating the mapped RNA-seq reads with IGV genome viewer. The binding site sequences can be applied in binding site predictions, binding site matrix calculations, as well as for investigations of the rate of evolution of transcription factor binding sites. Therefore, TFLink will facilitate benchmarking experiments in several fields of gene regulation research.

Use cases

To facilitate the application of TFLink, we provide some examples on how to use and process the data available at the gateway in form of descriptions, R scripts and unix shell commands:

  • Here we want to check and visualize the common target genes of two transcription factors. We describe how to find transcription factors that share common target genes. We cluster transcription factors based on their common target genes. We create a transcription factor–target gene interaction graph of the STAT5A and STAT5B transcription factors using the igraph R package. We also show how to create the same transcription factor–target gene interaction graph by using the Cytoscape software (Supplementary Use Cases, TFLink use case 1).

  • Here we investigate the functional diversity of target genes of a nuclear hormone receptor transcription factor, the unc-55 in human and a nematode species. We perform Gene Ontology overrepresentation analyses of the target genes in the two species in order to identify shared functional roles that likely represent the ancestral function of unc-55. Furthermore, this comparison will yield insights into the potentially divergent roles unc-55 play in these two distant animal groups (Supplementary Use Cases, TFLink use case 2).

  • Here, we investigate the binding sites of the EGR1 transcription factor. After converting the TFLink binding site table to BED and BAM files, we calculate the ‘coverage’ to reveal the strength of evidence (number of supporting experiments) for each binding site. Then, we plot the binding sites on the human chromosomes, indicating the number of supporting evidences each binding site has. Finally, we investigate specific binding sites using the IGV genome viewer tool (Supplementary Use Cases, TFLink use case 3).

Conclusions

TFLink is a gap-filling, novel, curated, multispecies gateway on eukaryotic transcription regulatory networks and transcription factor binding sites. It provides experimentally determined, highly accurate information collected from scattered third party databases. On the TFLink website (https://tflink.net/) the user can browse and download tables and files containing detailed information on transcription factor–target gene interactions and transcription factor binding sites. All the identifiers, names and files available on the TFLink gateway are provided in a standardized way or format to allow bioinformaticians to process these data easily. Besides being a practical resource for computational biologists, TFLink also targets wet-lab researchers with interactive network visualizations, cross-links to the original publications and to the UniProt database and allows the visual exploration of the genomic content around each binding site via links to the UCSC genome browser. We have established TFLink to facilitate the reuse of high-quality gene expression regulation data in an effective way, to provide an easy to use web interface and downloadable files, as well as to let the researchers launch comparative gene regulation studies.

Supplementary data

Supplementary data are available at Database Online.

Acknowledgements

We are grateful to Tibor Vellai for his valuable comments on the TFLink gateway. We thank Dóra Bokor, PharmD, for proofreading the manuscript.

Funding

This work was supported by the National Research, Development and Innovation Office, Hungary (NKFIH) grant PD [grant number 131839 to E.A.]; an NKFIH KKP [grant number 129814 to B.P.]; GINOP iChamber [grant number 2.3.2-15-2016-00026 to B.P.]; and The European Union’s Horizon 2020 research and innovation programme under grant agreement [grant number 739593 to B.P.]; a Biotechnology and Biological Sciences Research Council (BBSRC) Core Strategic Programme Grant [grant number BB/CSP17270/1 to T.K.]; and a BBSRC ISP grant [grant number BB/R012490/1 to T.K.].

Conflict of interest

None declared.

References

1.

Chen
K.
and
Rajewsky
N.
(
2007
)
The evolution of gene regulation by transcription factors and microRNAs
.
Nat. Rev. Genet.
,
8
,
93
103
.doi: .

2.

Galas
D.J.
and
Schmitz
A.
(
1978
)
DNAse footprinting: a simple method for the detection of protein–DNA binding specificity
.
Nucleic Acids Res.
,
5
,
3157
3170
.doi: .

3.

Garner
M.M.
and
Revzin
A.
(
1981
)
A gel electrophoresis method for quantifying the binding of proteins to specific DNA regions: application to components of the Escherichia coli lactose operon regulatory system+
.
Nucleic Acids Res.
,
9
,
3047
3060
.doi: .

4.

Pollock
R.
and
Treisman
R.
(
1990
)
A sensitive method for the determination of protein-DNA binding specificities
.
Nucleic Acids Res.
,
18
,
6197
6204
.doi: .

5.

de Wet
J.R.
,
Wood
K.V.
,
DeLuca
M.
et al.  (
1987
)
Firefly luciferase gene: structure and expression in mammalian cells
.
Mol. Cell. Biol.
,
7
,
725
737
.doi: .

6.

Ren
B.
,
Robert
F.
,
Wyrick
J.J.
et al.  (
2000
)
Genome-wide location and function of DNA binding proteins
.
Science
,
290
,
2306
2309
.doi: .

7.

Jolma
A.
,
Kivioja
T.
,
Toivonen
J.
et al.  (
2010
)
Multiplexed massively parallel SELEX for characterization of human transcription factor binding specificities
.
Genome Res.
,
20
,
861
873
.doi: .

8.

Rockel
S.
Geertz
M.
Maerkl
S.J.
(
2012
) MITOMI: a microfluidic platform for in vitro characterization of transcription factor–DNA interaction. In:
Deplancke
B
,
Gheldof
N
(eds).
Gene Regulatory Networks: Methods and Protocols, Methods in Molecular Biology
.
Humana Press
,
Totowa, NJ
, pp.
97
114
.

9.

Johnson
D.S.
,
Mortazavi
A.
,
Myers
R.M.
et al.  (
2007
)
Genome-wide mapping of in vivo protein–DNA interactions
.
Science
,
316
,
1497
1502
.doi: .

10.

Robertson
G.
,
Hirst
M.
,
Bainbridge
M.
et al.  (
2007
)
Genome-wide profiles of STAT1 DNA association using chromatin immunoprecipitation and massively parallel sequencing
.
Nat. Methods
,
4
,
651
657
.doi: .

11.

Han
H.
,
Cho
J.-W.
,
Lee
S.
et al.  (
2018
)
TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions
.
Nucleic Acids Res.
,
46
,
D380
D386
.doi: .

12.

Rivera
J.
,
Keränen
S.V.E.
,
Gallo
S.M.
et al.  (
2019
)
REDfly: the transcriptional regulatory element database for Drosophila
.
Nucleic Acids Res.
,
47
,
D828
D834
.doi: .

13.

Monteiro
P.T.
,
Oliveira
J.
,
Pais
P.
et al.  (
2020
)
YEASTRACT+: a portal for cross-species comparative genomics of transcription regulation in yeasts
.
Nucleic Acids Res.
,
48
,
D642
D649
.doi: .

14.

Bovolenta
L.A.
,
Acencio
M.L.
and
Lemke
N.
(
2012
)
HTRIdb: an open-access database for experimentally verified human transcriptional regulation interactions
.
BMC Genomics
,
13
, 405.doi: .

15.

Lesurf
R.
,
Cotto
K.C.
,
Wang
G.
et al.  (
2016
)
ORegAnno 3.0: a community-driven resource for curated regulatory annotation
.
Nucleic Acids Res.
,
44
,
D126
D132
.doi: .

16.

Fornes
O.
,
Castro-Mondragon
J.A.
,
Khan
A.
et al.  (
2020
)
JASPAR 2020: update of the open-access database of transcription factor binding profiles
.
Nucleic Acids Res.
,
48
,
D87
D92
.doi: .

17.

Garcia-Alonso
L.
,
Holland
C.H.
,
Ibrahim
M.M.
et al.  (
2019
)
Benchmark and integration of resources for the estimation of human transcription factor activities
.
Genome Res.
,
29
,
1363
1375
.doi: .

18.

Yevshin
I.
,
Sharipov
R.
,
Kolmykov
S.
et al.  (
2019
)
GTRD: a database on gene transcription regulation—2019 update
.
Nucleic Acids Res.
,
47
,
D100
D105
.doi: .

19.

Chèneby
J.
,
Gheorghe
M.
,
Artufel
M.
et al.  (
2018
)
ReMap 2018: an updated atlas of regulatory regions from an integrative analysis of DNA-binding ChIP-seq experiments
.
Nucleic Acids Res.
,
46
,
D267
D275
.doi: .

20.

Jiang
C.
,
Xuan
Z.
,
Zhao
F.
et al.  (
2007
)
TRED: a transcriptional regulatory element database, new entries and other development
.
Nucleic Acids Res.
,
35
,
D137
D140
.doi: .

21.

Hammal
F.
,
de Langen
P.
,
Bergon
A.
et al.  (
2021
)
ReMap 2022: a database of human, mouse, Drosophila and Arabidopsis regulatory regions from an integrative analysis of DNA-binding sequencing experiments
.
Nucleic Acids Res.
,
50
,
D316
D325
.doi: .

22.

Bohár
B.
,
Fazekas
D.
,
Madgwick
M.
et al.  (
2021
)
Sherlock: an open-source data platform to store, analyze and integrate big data for biology
.
F1000Research
,
10
, 409.doi: .

23.

Altenhoff
A.M.
,
Train
C.-M.
,
Gilbert
K.J.
et al.  (
2021
)
OMA orthology in 2021: website overhaul, conserved isoforms, ancestral gene order and more
.
Nucleic Acids Res.
,
49
,
D373
D379
.doi: .

24.

Liu
Z.-P.
,
Wu
C.
,
Miao
H.
et al.  (
2015
)
RegNetwork: an integrated database of transcriptional and post-transcriptional regulatory networks in human and mouse
.
Database
,
2015
, bav095.doi: .

25.

The UniProt Consortium
. (
2021
)
UniProt: the universal protein knowledgebase in 2021
.
Nucleic Acids Res.
,
49
,
D480
D489
.doi: .

26.

Bethesda (MD)
.
National Library of Medicine (US), National Center for Biotechnology Information (NCBI)
.
2021
.
Gene
, https://www.ncbi.nlm.nih.gov/gene/ (20 December 2021, date last accessed).

27.

Stormo
G.D.
,
Schneider
T.D.
,
Gold
L.
et al.  (
1982
)
Use of the ‘Perceptron’ algorithm to distinguish translational initiation sites in E. coli
.
Nucleic Acids Res.
,
10
,
2997
3011
.doi: .

28.

Schneider
T.D.
and
Stephens
R.M.
(
1990
)
Sequence logos: a new way to display consensus sequences
.
Nucleic Acids Res.
,
18
,
6097
6100
.doi: .

29.

Bethesda (MD)
.
National Library of Medicine (US), National Center for Biotechnology Information (NCBI)
.
2021
.
PubMed
, https://pubmed.ncbi.nlm.nih.gov/ (20 December 2021, date last accessed).

30.

Kent
W.J.
,
Sugnet
C.W.
,
Furey
T.S.
et al.  (
2002
)
The human genome browser at UCSC
.
Genome Res.
,
12
,
996
1006
.doi: .

31.

Shannon
P.
,
Markiel
A.
,
Ozier
O.
et al.  (
2003
)
Cytoscape: a software environment for integrated models of biomolecular interaction networks
.
Genome Res.
,
13
,
2498
2504
.doi: .

32.

Csardi
G.
and
Nepusz
T.
(
2006
)
The igraph software package for complex network research
.
Int. J. Complex Syst.
,
1695
.

33.

Sivade Dumousseau
M.
,
Alonso-López
D.
,
Ammari
M.
et al.  (
2018
)
Encompassing new use cases—level 3.0 of the HUPO-PSI format for molecular interactions
.
BMC Bioinform.
,
19
, 134.doi: .

34.

Subramanian
A.
,
Tamayo
P.
,
Mootha
V.K.
et al.  (
2005
)
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc. Natl. Acad. Sci. U.S.A.
,
102
,
15545
15550
.doi: .

35.

Thorvaldsdóttir
H.
,
Robinson
J.T.
and
Mesirov
J.P.
(
2013
)
Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration
.
Brief. Bioinform.
,
14
,
178
192
.doi: .

36.

Wilkinson
M.D.
,
Dumontier
M.
,
Aalbersberg
I.J.
et al.  (
2016
)
The FAIR guiding principles for scientific data management and stewardship
.
Sci. Data.
,
3
, 160018.doi: .

37.

Zheng
R.
,
Wan
C.
,
Mei
S.
et al.  (
2019
)
Cistrome data browser: expanded datasets and new tools for gene regulatory analysis
.
Nucleic Acids Res.
,
47
,
D729
D735
.doi: .

38.

Zhang
Q.
,
Liu
W.
,
Zhang
H.-M.
et al.  (
2020
)
hTFtarget: a comprehensive database for regulations of human transcription factors and their targets
.
Genom. Proteom. Bioinform.
,
18
,
120
128
.doi: .

39.

Hu
H.
,
Miao
Y.-R.
,
Jia
L.-H.
et al.  (
2019
)
AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors
.
Nucleic Acids Res.
,
47
,
D33
D38
.doi: .

40.

Wingender
E.
,
Chen
X.
,
Hehl
R.
et al.  (
2000
)
TRANSFAC: an integrated system for gene expression regulation
.
Nucleic Acids Res.
,
28
,
316
319
.doi: .

41.

Spivak
A.T.
and
Stormo
G.D.
(
2012
)
ScerTF: a comprehensive database of benchmarked position weight matrices for Saccharomyces species
.
Nucleic Acids Res.
,
40
,
D162
D168
.doi: .

42.

de Boer
C.G.
and
Hughes
T.R.
(
2012
)
YeTFaSCo: a database of evaluated yeast transcription factor sequence specificities
.
Nucleic Acids Res.
,
40
,
D169
D179
.doi: .

43.

Jin
J.
,
Tian
F.
,
Yang
D.-C.
et al.  (
2017
)
PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants
.
Nucleic Acids Res.
,
45
,
D1040
D1045
.doi: .

44.

Luo
Y.
,
Hitz
B.C.
,
Gabdank
I.
et al.  (
2020
)
New developments on the Encyclopedia of DNA elements (ENCODE) data portal
.
Nucleic Acids Res.
,
48
,
D882
D889
.doi: .

45.

Edgar
R.
,
Domrachev
M.
and
Lash
A.E.
(
2002
)
Gene expression omnibus: NCBI gene expression and hybridization array data repository
.
Nucleic Acids Res.
,
30
,
207
210
.doi: .

46.

Kodama
Y.
,
Shumway
M.
,
Leinonen
R.
et al.  (
2012
)
The sequence read archive: explosive growth of sequencing data
.
Nucleic Acids Res.
,
40
,
D54
D56
.doi: .

47.

Vazquez
M.
,
Krallinger
M.
,
Leitner
F.
et al.  (
2022
)
ExTRI: extraction of transcription regulation interactions from literature
.
Biochim. Biophys. Acta Gene Regul. Mech.
,
1865
, 194778.doi: .

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial License (https://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com

Supplementary data