Prediction of Transcription Factor Regulators and Gene Regulatory Networks in Tomato Using Binding Site Information

Gene regulatory networks (GRNs) represent the regulatory links between transcription factors (TF) and their target genes. In plants, they are essential to understand transcriptional programs that control important agricultural traits such as yield or (a)biotic stress response. Although several high- and low-throughput experimental methods have been developed to map GRNs in plants, these are sometimes expensive, come with laborious protocols, and are not always optimized for tomato, one of the most important horticultural crops worldwide. In this chapter, we present a computational method that covers two protocols: one protocol to map gene identifiers between two different tomato genome assemblies, and another protocol to predict putative regulators and delineate GRNs given a set of functionally related or coregulated genes by exploiting publicly available TF-binding information. As an example, we applied the motif enrichment protocol on tomato using upregulated genes in response to jasmonate, as well as upregulated and downregulated genes in plants with genotypes OENAM1 and nam1, respectively. We found that our protocol accurately infers the expected TFs as top enriched regulators and identifies GRNs functionally enriched in biological processes related with the experimental context under study.