Study of regulatory elements in promoter regions of co-expressed genes in tomato when facing low temperatures

Document Type : Research Paper

Authors

1 Assistant Professor, Department of Agronomy and Plant Breeding, Faculty of Agriculture, Shahrood University of Technology, Shahrood, Iran.

2 Postdoc Researcher, INET Institute, Munchen, Germany.

Abstract

Next-generation sequencing (NGS) techniques such as RNA-Seq can provide expression information per gene and can be used to construct co-expression networks. These networks can produce valuable information to understand plant responses when facing environmental changes. In this study, a complete profile of tomato gene expressions under low temperature treatment was analyzed by “Weighted Gene Co-expression Network Analysis”. The clustering results allowed to group 20,267 genes into eight modules. The results of Gene ontology (GO) term analysis highlighted different functions between co-expression groups in terms of biological process and molecular function. Accordingly, the first module overrepresented to heat shock proteins (HSPs). Protein kinases and proteins involved in post-translational modification were in the second and third modules, respectively. The promoter of genes in fifth module was screened according to the different expression patterns in the tomato genotypes. The cis-elements such as ACGTG and TTGAC were observed in promoter site of all genes. The transcription factor families of AP2 / ERF and ZF-HD have the highest binding status in promoter of genes of the fifth module. The results present the key cis-regulatory elements that effect on co-expressed genes network and involve in anti-oxidant process under cold stress.

Keywords

Main Subjects


Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biology. doi: 10.1186/gb-2010-11-10-r106.
Aoki K, Ogata Y, Shibata D (2007) Approaches for extracting practical information from gene co-expression networks in plant biology. Plant cell physiology 48: 381-390.
Bergmann S, Ihmels J, Barkai N (2004) Similarities and differences in genome-wide expression data of six organisms. PLoS Biol. 2(1): 85-93.
Bhattacharjee A, Ghangal R, Garg R, Jain M (2015) Genome-Wide Analysis of Homeobox Gene Family in Legumes: Identification, Gene Duplication and Expression Profiling. PLoS ONE 10(3): e0119198. https://doi.org/10.1371/ journal.pone.0119198.
 
1. Systemic acquired resistance
2. Sulfur-responsive element
Brand LH, Fischer NM, Harter K, Kohlbacher O, Wanke D (2013) Elucidating the evolutionary conserved DNA-binding specificities of WRKY transcription factors by molecular dynamics and in vitro binding assays. Nucleic Acids Research. 41(21): 9764-9778.
 Buchel AS, Brederode FT, Bol JF, Linthorst HJ (1999) Mutation of GT-1 binding sites in the Pr-1A promoter influences the level of inducible gene expression in vivo. Plant. Mol. Biol. 40: 387-396.
Cai B, Li CH, Huang J (2014) Systematic identification of cell-wall related genes in populous based on analysis of functional modules in co-expression network. PloS One. 9: e95176. Doi:10.1371/journal.pone.0095176.
Chen H, Chen X, Chen D, Li J, Zhang Y, Wang A (2015) A comparison of the low temperature transcriptomes of two tomato genotypes that differ in freezing tolerance: Solanum lycopersicum and Solanum habrochaites. BMC Plant Biology. 15:132. doi:10.1186/s12870-015-0521-6.
Chen C, Chen Z (2002) Potentiation of Developmentally Regulated Plant Defense Response by AtWRKY18, a Pathogen-Induced Arabidopsis Transcription Factor. Plant Physiology. 129(2): 706-716.
Chinnusamy V, Ohta M, Kanrar S, Lee B, Hong X (2007) ICE1: a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis. Genes & Development 17: 1043-1054.
Chow CN, Zheng H.Q, Wu NY, Chien C.H, Huang HD, Lee TY, Chang W.C (2016) PlantPAN 2.0: an update of plant promoter analysis navigator for reconstructing transcriptional regulatory networks in plants. Nucleic Acids Research. 44 (Database issue): D1154-1160.
de Silva K, Laska B, Brown C, Sederoff H, Khodakovskaya M (2011) Arabidopsis thaliana calcium-dependent lipid-binding protein (AtCLB): a novel repressor of abiotic stress response. Journal of Experimental Botany 62(8): 2679-2689.
Falcon S, Gentleman R (2006) Using GOstats to Test Gene Lists for GO Term Association. Bioinformatics 23(2): 257-258.
Fernandez-Pozo N, Menda N, Edwards JD, Saha S, Tecle IY, Strickler SR, Mueller L. A (2015) The Sol Genomics Network (SGN)-from genotype to phenotype to breeding. Nucleic Acids Research. 43(Database issue): D1036-D1041.
Haralampidis K, Milioni D, Rigas S, Hatzopoulos P (2002) Combinatorial Interaction of Cis Elements Specifies the Expression of the Arabidopsis AtHsp90-1 Gene. Plant Physiology. 129(3):1138-1149.
Heidari P, Ahmadizadeh M, Najafi-Zarrini H (2015) In Silico Analysis of Cis-Regulatory Elements on Co-Expressed Genes. J. Biol. Environ. Sci. 9(25): 1-9.
Hudson ME, Quail PH (2003) Identification of Promoter Motifs Involved in the Network of Phytochrome A-Regulated Gene Expression by Combined Analysis of Genomic Sequence and Microarray Data. Plant Physiology 133(4): 1605-1616.
Imamura A, Hanaki N, Umeda H, Nakamura A, Suzuki T, Ueguchi C, Mizuno T (1998) Response regulators implicated in His-to-Asp phosphotransfer signaling in Arabidopsis. Proceedings of the National Academy of Sciences of the United States of America. 95(5): 2691-2696.
Langfelder P, Horvath S (2008) WGCNA: An R package for weighted correlation network analysis. BMC bioinformatics. 9(1): 559. https://doi.org/10.1186/1471-2105-9-559.
Licausi F, Ohme-Takagi M, Perata P (2013) APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytologist. 199 (3): 639-649.
Lindemose S, O’Shea C, Jensen MK, Skriver K (2013) Structure, function and networks of transcription factors involved in abiotic stress responses. International Journal of Molecular Sciences. 14: 5842-5878. 56. 
Liu H, Ouyang B, Zhang J, Wang T, Li H, et al. (2012) Differential Modulation of Photosynthesis, Signaling, and Transcriptional Regulation between Tolerant and Sensitive Tomato Genotypes under Cold Stress. PLoS ONE. 7(11): e50785. https://doi.org/ 10.1371/journal.pone.0050785.
Mailloux RJ, Bériault R, Lemire J, Singh R, Chénier DR, Hamel RD, et al. (2007) The Tricarboxylic Acid Cycle, an Ancient Metabolic Network with a Novel Twist. PLoS ONE 2(8): e690. https://doi.org/10.1371/journal.pone.0000690
Maruyama-Nakashita A, Nakamura Y, Watanabe-Takahashi A, Inoue E, Yamaya T, Takahashi H (2005) Identification of a novel cis-acting element conferring sulfur deficiency response in Arabidopsis roots. Plant Journal 42(3): 305-314.
Nakashima K, Fujita Y, Katsura K, Maruyama, K, Narusaka Y, Seki M, Shinozaki K, Yamaguchi-Shinozaki K (2006) Plant Mol Biol. 60: 51-68.
Planchais S, Perennes C, Glab N, Mironov V, Inzé D, Bergounioux C (2002) Characterization of cis-acting element involved in cell cycle phase-independent activation of Arath; CycB1;1 transcription and identification of putative regulatory proteins. Plant Mol Biol 50: 109-125.
Ross E, Stone JM, Elowsky CG, Arredondo-Peter R, Klucas R, Sarath G (2004) Activation of the Oryza sativa non-symbiotic haemoglobin-2 promoter by the cytokinin-regulated transcription factor, ARR1. Journal of Experimental Botany. 55 (403): 1721-1731.
Saibo NJ, Lourenco T, Oliveira MM (2009) Transcription factors and regulation of photosynthetic and related metabolism under environmental stresses. Annals of Botany. 103: 609-623.
Sornaraj P, Luang S, Lopato S, Hrmova M (2016) Basic leucine zipper (bZIP) transcription factors involved in abiotic stresses: A molecular model of a wheat bZIP factor and implications of its structure in function. Biochimica et Biophysica Acta (BBA). 1860 (1): 46-56.
Supek F, Bošnjak M, Škunca N, Šmuc T (2011) REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS ONE 6(7): e21800. https://doi.org/10.1371/journal.pone.0021800.
Verma G, Dhar YV, Srivastava D, Kidwai M, Chauhan PS, Bag SK, et al. (2017) Genome-wide analysis of rice dehydrin gene family: Its evolutionary conservedness and expression pattern in response to PEG induced dehydration stress. PLoS ONE 12(5): e0176399. https://doi.org/10.1371/journal.pone.0176399.
Yadav V, Mallappa C, Gangappa SN, Bhatia S, Chattopadhyay S (2005) A Basic Helix-Loop-Helix Transcription Factor in Arabidopsis, MYC2, Acts as a Repressor of Blue Light–Mediated Photomorphogenic Growth. The Plant Cell. 17(7): 1953-1966.
Yu D, Chen C, Chen Z (2001) Evidence for an important role of WRKY DNA binding proteins in the regulation of NPR1 gene expression. Plant Cell. 13: 1527-1540.
Zhou j, Li C, Gao F, Luo X, Li Q, Zhao H, Yao H, Chen H, Wang A, Wu Q (2016) Characterization of Three Glucosyltransferase Genes in Tartary Buckwheat and Their Expression after Cold Stress. Journal of Agricultural and Food Chemistry 64 (37): 6930-6938.