Shikha Mittal1
Pooja Banduni1
Mallana G. Mallikarjuna1
Atmakuri R. Rao2
Prashant A. Jain3
Prasanta K. Dash4
Nepolean Thirunavukkarasu1*- 1Division of Genetics, ICAR-Indian Agricultural Research Institute, New Delhi, India
- 2Centre for Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi, India
- 3Department of Computational Biology & Bioinformatics, J.I.B.B., Sam Higginbottom University of Agriculture, Technology and Sciences, Allahabad, India
- 4National Research Centre on Plant Biotechnology, New Delhi, India
Drought is one of the major threats to the maize yield especially in subtropical production systems. Understanding the genes and regulatory mechanisms of drought tolerance is important to sustain the yield. Transcription factors (TFs) play a major role in gene regulation under drought stress. In the present study, a set of 15 major TF families comprising 1,436 genes was structurally and functionally characterized. The functional annotation indicated that the genes were involved in ABA signaling, ROS scavenging, photosynthesis, stomatal regulation, and sucrose metabolism. Duplication was identified as the primary force in divergence and expansion of TF families. Phylogenetic relationship was developed for individual TF and combined TF families. Phylogenetic analysis clustered the genes into specific and mixed groups. Gene structure analysis revealed that more number of genes were intron-rich as compared to intron-less. Drought-responsive cis-regulatory elements such as ABREA, ABREB, DRE1, and DRECRTCOREAT have been identified. Expression and interaction analyses identified leaf-specific bZIP TF, GRMZM2G140355, as a potential contributor toward drought tolerance in maize. Protein-protein interaction network of 269 drought-responsive genes belonging to different TFs has been provided. The information generated on structural and functional characteristics, expression, and interaction of the drought-related TF families will be useful to decipher the drought tolerance mechanisms and to breed drought-tolerant genotypes in maize.
Introduction
Maize serves as a model system for genetic studies from simple inheritance, physical linkage, cytological crossing over, and recombination analyses to modern epigenetic silencing, genome imprinting, and transposition (Bennetzen and Hake, 2009). Although maize production is considerably affected by many biotic and abiotic stresses, drought is the major one (Leach et al., 2011). Being a water-intensive field crop, shortage of water severely limits the productivity of maize. The stress manifests through several symptoms, including increased leaf senescence and smaller leaves. Cellular damage and impaired carbon fixation during drought lead to imbalanced carbon assimilation and yield reduction (Hall et al., 1981). Drought is thus a serious threat to food security in all maize production systems.
During drought, plant TFs promote multiple responses at the molecular level by regulating gene expression and develop interactions with basal TFs at target gene promoters in response to developmental and intracellular signals, and control the target gene function either by activation or repression (Lu et al., 2012). Through drought-TF interfaces, the plants can tightly control gene expression in response to drought or water stress signals. Many drought-responsive TFs, including AP2, bZIP, NAC, HD-Zip, MYB, and NF-YA, play an active role in drought regulation, and the function of these genes is important for drought tolerance (Mittal et al., 2017; Thirunavukkarasu et al., 2017).
The structural and functional characterization of transcription factors have been studied in different crops. Physical mapping of WRKY in Setaria italica, NAC in Brachypodium distachyon, and MYB in Arabidopsis and rice revealed an uneven distribution of genes due to distinct duplication patterns (Katiyar et al., 2012; Muthamilarasan et al., 2015). Grouping of TFs into different families could be done based on 3D structural similarity in DNA-binding and acidic domains (Yilmaz et al., 2009). Phylogenetic studies on HD-Zip in maize, and bZIP in cucumber and rice also showed the formation of different clades on the basis of similarities and divergences (Baloglu et al., 2014; Mao et al., 2016). TF families have specific roles in plant response to abiotic stress; for example, the basic leucine zipper (bZIP) family consists of a common bZIP domain for DNA binding at the N terminus and a leucine-rich motif for dimerization at the C terminus and functions in an ABA-dependent manner by associating with ABRE (ABA-responsive cis-regulatory element) in its promoter regions (Zou et al., 2008). WRKY family has a conserved WRKY domain that binds to the W-box (a conserved stretch of 60 amino acids) to control the expression of genes (Rushton et al., 2010). All WRKY genes contain a conserved amino acid sequence (WRKYGQK) at the N terminus and a zinc finger motif at the C terminus. Several studies have shown that the role of WRKY in drought response is through modulation of an ABA-dependent pathway.
In maize, different TF families have been identified as responsible for modulating the gene regulation in response to drought stress. TFs such as MYB, NAC, WRKY, bZIP, bHLH, HD-zip via ABA-dependent or ABA-independent pathways play a significant role in drought tolerance. ZmNF-YB2 in maize is shown to have an equal role as AtNF-YB1 in Arabidopsis to confer improved crop performance under drought conditions (Nelson et al., 2007). bZIP plays a vital role in ABA signaling along with other functions in plant growth and abiotic stresses. ZmbZIP72, a bZIP TF in maize was over-expressed in various organs in response to drought, salinity, and ABA stress at seedling stage. AP2, ERF, dehydration-responsive element-binding protein (DREB), Cys2His2 Zinc Finger (C2H2 ZF), MYB, bHLH were identified as important TFs for drought tolerance (Ying et al., 2012). A set of 10 drought-responsive TF families (AP2, ERF, MYB, bZIP, bHLH, GRAS, WRKY, NAC, NF-YA, and NF-YB) have been identified in maize to regulate different physiological and molecular functions namely, stomatal regulation, hormone signaling, root development, and osmoregulation (Mittal et al., 2017).
The role of TFs in various stress tolerances in different crop plants including maize is well-documented. However, there is no substantial report on comprehensive structural and functional analysis of the major drought TF families in maize. Therefore, the present investigation has been designed to identify and characterize major TF families related to drought tolerance. These families were classified based on their putative DNA binding domain and the transcription regulatory domain (TRD) or transcriptional activation domain (TAD). The study covered the following aspects: detailed analysis of TF-specific domains, physical mapping of the TF families, prediction of gene functions, phylogeny, duplication events, paralogs, gene structure, cis-regulatory elements (CREs) and motifs analysis, expression assay, and protein-protein interaction analysis. A hypothesis on how different TFs by binding to cis-regulatory elements regulate drought tolerance in maize has been developed. The results provide insights into the role of TF families in relation to drought. Further a functional analysis has identified important TF genes which can be used as candidate for drought breeding.
Materials and Methods
Identification of TF Gene Families and Domain Analysis in Maize
A set of 15 major drought-related TF families namely, AP2, ARF, bHLH, bZIP, C2H2, ERF, GATA, GRAS, HD-Zip, HSF, MYB, NAC, NF-YA, NF-YB, and WRKY were selected for structural and functional characterization. These TF families were considered as major drought-related TF families (Thirunavukkarasu et al., 2014; Mittal et al., 2017; Van Gioi et al., 2017) and a total of 1,436 genes representing those families were selected from the Plant Transcription Factor Database 3.0 (http://planttfdb.cbi.pku.edu.cn) and in-house data (Thirunavukkarasu et al., 2014; Mittal et al., 2017). We noticed that although some member proteins lacked a specific TF family domain, they were nevertheless annotated as integral members of that TF family. Therefore, we excluded such redundant member proteins from consideration. To identify the number of domains present in a TF family protein, we searched the Pfam database (http://pfam.xfam.org) with an expectation cut-off (E-value) of 1.0, and only the significant domains found in the maize protein sequences were considered as valid domains.
Chromosomal Distribution
The chromosomal position of each gene encoding a maize protein was determined by Phytozome database (https://phytozome.jgi.doe.gov/pz/portal.html). Each TF family was then mapped onto the maize genome (n = 10). Genes were plotted separately onto 10 chromosomes of maize by using R package “synbreed.”
Gene Ontology
Nucleotide and protein sequences of the 1,436 TF genes and primary information about them, including the chromosome on which they are located, were retrieved from the Ensembl plant's website (http://asia.ensembl.org/index.html). GO annotations of all maize TF family proteins were executed using the Blast2GO server (https://www.blast2go.com). The Blast2GO annotation links the genes with the GO terms using hierarchical vocabularies. The GO terms have been classified into three categories: biological process, cellular components, and molecular function. The visual display of GO annotation has been drawn using WEGO tool.
Phylogenetic Analysis
Phylogenetic analysis was performed using full-length protein sequences to understand the evolutionary relationship among and within the TFs. To generate the phylogenetic tree for TFs families, amino acid sequences were imported into MEGA 6.06 to perform multiple sequence alignment with gap open and gap extension penalties of 10 and 0.1, respectively. Then the alignment file obtained using MSA was used for the construction of unrooted phylogenetic tree using neighbor-joining method through MEGA ver. 6.06 (http://www.megasoftware.net) and was rendered using iTOL (http://itol.embl.de).The reliability of the constructed unrooted tree was established using 1,000 rapid bootstrap replicates.
Gene Duplication and Paralog Identification
Information on gene duplication events, such as tandem duplication and block duplication, was retrieved using PLAZA, a tool (http://bioinformatics.psb.ugent.be/plaza/) available at Plant Genome Duplication Database (PGDD) (http://chibba.agtec.uga.edu/duplication). Each gene was aligned against the 1,436 genes, and paralogs with a minimum of 75% identity were searched using BLAST (https://blast.ncbi.nlm.nih.gov). Adjacent genes of the same sub-family located within 10 predicted genes apart or within 30 Kbp of each other were characterized as tandem duplications (Shiu and Bleecker, 2003; Du et al., 2013).
Gene Structure Analysis
Gene structure or the exon-intron structure of the genes was identified using the Gene Structure Display Server 2.0 (http://gsds.cbi.pku.edu.cn). The genomic sequences were aligned with their corresponding coding sequences for identifying the exon-intron structure in GSDS 2.0. The genomic and coding sequences of the genes were extracted from Ensembl genomes BioMart (http://ensemblgenomes.org).
Motif Prediction
The conserved protein and nucleotide motifs for all 15 TF families were discovered using MEME 4.10.2, the multiple maximizations for motif elicitation analysis tool (http://meme-suite.org/doc/download.html). The MEME suite was used to search 5 motifs of each TF family with a base width of a minimum of 12 and maximum of 60 Kb and an E-value smaller than 0.01.
Identification of cis-Regulatory Elements
Upstream promoter sequences (1 Kb) of all the major drought-related TF families have been retrieved using Phytozome Biomart (https://phytozome.jgi.doe.gov/pz/portal.html). These sequences were then used to predict and locate the cis-elements using PlantPAN (http://plantpan2.itps.ncku.edu.tw).
Expression Analysis of Drought-Responsive TF Genes in Zea mays
Expression analysis of TF genes was performed using maize RNA-Seq data generated under drought stress. A total of two independent RNA-Seq experiments (Liu et al., 2015; Thatcher et al., 2016) available in the NCBI's SRA database for Zea mays were considered for in-silico expression analysis (Accession numbers: SRP061631 and SRP062027; https://www.ncbi.nlm.nih.gov). In the first experiment, a drought-tolerant maize genotype ZD619 was used to generate RNA-Seq data under drought stress and controlled condition using root, stem, and leaves. In the second experiment, RNA-Seq data was generated using B73 under drought stress in different growth stages such as V12, V14, V18, and R1. Expression data from the above-mentioned experiments were considered for the validation of genes.
Protein Interaction Network Analysis of Drought-Responsive TFs
Based on the expression data TF interactions were studied for the selected 269 drought-responsive genes. The STRING 10.0 (Search Tool for the Retrieval of Interacting Genes/Proteins) was employed to analyze and visualize the TF interaction network (http://string-db.org). The network was constructed using the confidence score of 0.5 to assure a high reliability of the interaction network and all other parameters were the default.
Results and Discussion
In the present study, major drought-related TF families have been selected on the basis of their role in drought tolerance in maize and in different species revealed by literature and in-house experiments (Thirunavukkarasu et al., 2014; Mittal et al., 2017; Van Gioi et al., 2017). For example, SNAC1 gene confers drought tolerance in rice by regulation of downstream genes (You et al., 2014). NAC has been characterized in different species including Chinese cabbage, soybean, maize, rice, Arabidopsis and grape (Shao et al., 2015). The over-expression of ANACO19, ANACO55, and ANACO72 elevated the drought expression in Arabidopsis (Joshi et al., 2016). In Arabidopsis, rice, and soybean, AREB1 over-expression improves the drought tolerance (Oh et al., 2005; Barbosa et al., 2013; Yoshida et al., 2015). These major 15 drought-related TF families comprise of 1,436 genes of which, MYB family had the highest number of genes (281), followed by ERF (181), whereas NF-YA had the lowest (16) (Supplementary Table S1).
Chromosomal Distribution
To investigate the genome organization of drought-related TF genes, all 1436 genes were mapped onto different chromosomes of maize. Uneven distribution of genes was found in the genome where chromosome 1 was gene-rich, containing 15% of total genes. Among all TFs, MYB family identified on chromosome 3 showed the highest gene density, whereas NF-YA showed the lowest gene density. AP2 and NF-YB is absent in chromosome 4 and 10, respectively. Chromosome 7 was devoid of the genes from ARF and GATA, and chromosomes 4, 6, 8, and 9 showed no genes from NF-YA (Figure 1A).
Figure 1. Chromosomal mapping & Blast2GO annotation of 1436 genes. (A) Genes associated with 15 drought-related TFs mapped on 10 chromosomes of maize according to their position. Different colors correspond to different TF families. (B) Genes of TF families representing various functional categories.
Genes from ARF, bHLH, bZIP, C2H2, ERF, GRAS, HD-Zip, NAC, and WRKY were discretely clustered on the chromosomes, whereas genes from AP2, GATA, HSF, NF-YA, and NF-YB, were more widely and sparsely distributed. Earlier research based on the physical and genetic framework for the reference genome (B73) has also shown such discrete distribution, which accounts for the genomic variability of the 10 chromosomes of maize and the underlying diverse changes associated with chromatin composition (Schnable et al., 2009). Additionally, movement of transposable elements during maize evolution might have contributed to the distribution of TF genes on the maize genome (Tian et al., 2009).
Annotation of TF Families
The GO annotation of TF genes was carried out to understand the involvement of the genes in different biological and molecular functions. The total GO terms identified in 1,436 genes belonging to 15 drought-related TFs were 22,506 (Supplementary Table S2). Of which, 63.30% of genes participated in biological process (BP), 40.53% in cellular components (CC), and 88.51% in molecular function (MF) (Figure 1B). Approximately 90% of the genes were predictive of their molecular function, and the most common sub-function in the molecular category was DNA binding. The GO annotations also indicated that 33% of the genes were drought-related and most of them were involved in ABA signaling, ROS scavenging, photosynthesis, stomatal development, and sucrose metabolism. The major drought-related function in ERF, WRKY, and bZIP was ABA signaling; similarly, defense and detoxification response were the major drought-related functions for MYB and C2H2. Among all genes, DREB, ABSCISIC ACID-INSENSITIVE 5-LIKE PROTEIN, ZINC FINGER PROTEIN 2, CUP-SHAPED COTYLEDON, V-TYPE PROTON ATPase SUBUNIT H, GLUTATHIONE S-TRANSFERASE, and SPEECHLESS were the important genes for stress response. It was proved that over-expression of DREB2ca in Arabidopsis and soybean improved drought tolerance without any growth impairment (Engels et al., 2013). Even in mature dry seeds, the over-expression of ABI5 regulates LEA genes (Delseny et al., 2001), thus had a role in dehydration tolerance. In Arabidopsis, various ZFPs such as AZF1, AZF2, AZF3, ZAT6, ZAT7, and ZAT10 were identified for drought-responsiveness (Sakamoto et al., 2000; Ciftci-Yilmaz and Mittler, 2008). ZAT12 played a significant role in ROS scavenging (Davletova, 2005) thus helps in stress tolerance. CUP-SHAPED COTYLEDON was involved in abiotic stress tolerant mechanisms (Jogaiah et al., 2013; Wang et al., 2016). In wheat, a V-H+-PPase gene, TVP1, was identified for improving the drought tolerance (Brini et al., 2007). The role of drought-induced expression of GLUTATHIONE S-TRANSFERASE was well documented in rice (Sharma et al., 2014). SPEECHLESS and MUTE, belonging to the bHLH family, responded to drought stress through regulating stomatal development, meristemoid differentiation, and guard cell morphogenesis (Pillitteri et al., 2007). All drought-related genes identified by their GO annotation can be considered as candidates and could be further validated in maize genotypes.
Phylogenetic Relationships Among and Within Drought-Related TF Families
Phylogenetic analysis was performed to understand the evolutionary relationship among and within drought-related TF families of maize. Evolutionary association among TFs indicated that the TF genes fell in the mixed clusters category, including those from the MYB family, in which about 99% were clustered into three groups. Of which two were well-differentiated (Figure 2A) while the third one was mixed with genes distributed in different clusters. The remaining 1% was present in association with their homologous genes from other TF families, including ERF, NAC, GRAS, bZIP, and C2H2. Nearly 10–15% of genes in the latter TFs were scattered outside their specific clusters, which suggested that the irregular distribution of few genes outside the TF-family could be due to their homologous association with other TF families (Figure 2A). The mixed grouping of genes belonging to various TFs suggested that these TFs could have evolved from common ancestral genes during the process of evolution. Studies have reported that mixed clustering of gene families based on sequence homology correlates with their function (Mochida et al., 2010). Several rounds of duplication at the genome level are the major contributors toward the expansion of TF families in plant evolution, and many TF families expanded in parallel with the evolution of multi-cellularity (Pires and Dolan, 2010). The lineage-specific parallel expansion of TF genes was reported in Arabidopsis-rice TF orthologous groups (Shiu et al., 2005). Evidences suggested that many of the eukaryotic genomes, including maize, have undergone multiple large-scale duplication events in their evolutionary history (Edger and Pires, 2009). Furthermore, joint contribution of trans-regulatory changes coupled with cis-regulatory changes accelerated the functional diversification of TFs (Schmitz et al., 2016). The phylogenetic association within the TF genes in maize (Figure 2B) showed that the majority of the TF trees were bifurcated into five clades (AP2, ARF, bZIP, and ERF), followed by three (bHLH, NF-YA, and NF-YB) and six clades (GRAS, HSF, and MYB). The maximum number of clades was observed in NAC TFs (10) followed by GATA and WRKY (7 each). However, C2H2 TF tree showed a significant number of polytomies and the presence of polytomies in C2H2 tree (Figure 2B: C2H2) indicated the rapid evolution of C2H2 genes. Schmitz et al. reported that repeats of sequence motifs such as C2H2 zinc fingers are responsible for the accelerated expansions of TF families (Schmitz et al., 2016).
Figure 2. An unrooted phylogenetic tree of drought-related TF gene families. (A) An unrooted phylogenetic tree of the maize drought related TF families comprising of 1436 TF genes from 15 TF families. (B) An unrooted phylogenetic tree of TF gene families. The phylogenetic relationship of the TF genes was analyzed and depicted within the family viz., AP2, ARF, bHLH, bZIP, C2H2, ERF, GRAS, GATA, HD-ZIP, HSF, MYB, NAC, NF-YA, NF-YB, and WRKY.
Same TF-families does have multiple functions since they regulate different mechanisms which improve drought tolerance in maize as well as in other species. MYB family has been involved in different functions including, plant development, hormonal signal transduction, and abiotic stress tolerance (Allan et al., 2008; Cominelli and Tonelli, 2009). In maize, rice, Arabidopsis, and soybean, MYB members play significant role in variety of cell processes such as cell cycle and cell morphogenesis (Jin and Martin, 1999; Stracke et al., 2001; Feller et al., 2011). Different members of WRKY family play vital role in abiotic stress response, germination, plant development, and senescence (Rushton et al., 2010).
According to our results, we found that the most of the drought-related genes identified using gene expression assay were found in specific sub-clades. For example—genes belonging to HD-Zip, HSF, ERF, WRKY, and MYB were present in the cluster in different specific sub-clades of the evolutionary tree. Most of the TF genes were present in specific clades while very few genes were scattered across clades. The result indicated that the genes belonging to a specific TF group play specific function under given condition. On the other hand, the scattered ones could have shared the common functional relationship along with their respective TF group.
Gene Duplication and Paralog Identification
Gene duplication is considered as one of the crucial events in evolution. To elucidate the expanded mechanism of drought-related TF genes in maize during evolution process, the duplication events including tandem and segmental were studied. Duplication events are the main factor responsible for expansion of gene families. Duplication events, such as tandem and block duplications, have paved the way for new gene functions and, in turn, evolutionary freshness. The maize genome has already undergone many rounds of genome duplications, including a paleopolyploid event ~70 million years ago (Blanc and Wolfe, 2014) and a subsequent whole genome duplication event ~11 million years ago (Blanc and Wolfe, 2014). Further studies based on gene duplications in maize have confirmed that a second large duplication event is responsible for gene duplication in most of the gene families (Schnable et al., 2009). For ease of analysis, we classified all the duplication events into three types: Type I, which corresponds only with tandem duplications; Type II, only with block duplications; and Type III, with both tandem and block duplications. In the present study, the largest number of Type I, Type II, and Type III events were present in NAC (31%), HD-Zip (64%), and WRKY (10%), respectively. Type II events were more frequent than Type I events, and only 11 TFs fell into the Type III category. The range of Type I, II, and III events was 0–16.27%, 6.54–16.57%, and 6.10–15.85%, respectively. Furthermore, Type I and Type II events were prominent on chromosome 1, whereas chromosome 8 was rich in Type III events (Supplementary Figure S1). These events signify that tandem duplication was a common one for rapidly evolving families; block duplication was a feature of the slowly evolving gene families, and both played a significant role in the expansion of gene families (Liu et al., 2011).
Subsequent identification and analysis of paralogs for each TF confirmed that only 981 (68.31%) genes had paralogs out of 1,436 genes. The distribution of paralogs across all chromosomes varied from 6.62 to 16.30%: the highest number of paralogs (160 genes) was found on chromosome 1 and the lowest on chromosome 10 (65 genes).
Gene Structure and Intron Distribution
To reveal the structural features of drought-related TFs, gene structure was analyzed using GSDS server. The study of gene structure is crucial in biology and would provide important clues on gene evolution. The structural orientation of intron-exon helps to determine the evolutionary path of a gene family and provides an opportunity to investigate the evolutionary mechanisms of the origin of gene families. In multigene families, variation in gene structure forms the basis of evolution (Cao and Shi, 2012). To provide insight into this process, we determined the number of exons and introns with respect to each gene in each TF family. Further, genes were divided into intron-rich and intronless. It was found that 72% were intron-rich and 28% were intronless. The ERF family contributed to the major share of intronless genes (77%), whereas NF-YA showed only 6.25% of genes as intronless. Among intron-rich genes, 21% had a single intron, 23% had two, 8% had three, and 20% had four or more (Figure 3).
Figure 3. Details of gene structure: Colored bars correspond to the number of introns in each gene family; different colors are assigned on the basis of the number of introns within each gene family.
The intronless genes in eukaryotes could have evolved either through horizontal gene transfer from ancient prokaryotes or duplication of existing intronless genes and retrotransposition of intron-containing genes (He et al., 2011). The results indicated that the gain or loss of exons during evolution could be responsible for the functional diversity of these different gene families including drought tolerance. These loss and gain of introns are essential for structural diversity and complexity of genes (Zhang and Kishino, 2004).
Conserved motifs
Motifs play an important role in the transcriptional regulation. To understand the similarity and diversity of gene motifs in different TF genes, five conserved motifs were identified. The importance of motifs generally lies in their association with the structure and function of proteins because protein networks are regulated by specific peptide motifs. Domain reorganization and sequence differentiation in proteins are the essential prerequisites for novel protein functions (Björklund et al., 2005). In the present study, most of the identified peptide motifs were extremely TF-specific, such as a 28-amino-acid-long sequence that was 100% conserved in all ARF genes and a 41-amino-acid-long sequence with 100% conservation in AP2 genes (Supplementary Table S3). In contrast, a 39-amino acid-long sequence conserved only in 81% genes was the least conserved sequence in HSF (Figure 4). Similarly, consistent results were obtained in nucleotide motifs: dissimilar types of 60-nucleotide-long motifs were 100% conserved in ARF, HD-Zip, bHLH, GATA, and NF-YB whereas a 28-nucleotide-long sequence conserved in 94% genes, was the least conserved motif in NF-YA (Figure 5). To gain deeper insight into the structures of conserved motifs, their composition matches were analyzed. The analysis revealed the significant features of unique and highly conserved motifs to each particular TF: for example, arginine was found at high frequencies in motif-2 (AP2), motif-2 (ARF), motif-1 (bZIP), motif 1 (HD-Zip), and motif-1 (NF-YA). In nucleotide motifs, guanine was the most conserved nucleotide across all conserved nucleotide motifs, such as motif-1 (bHLH, C2H2, GATA, GRAS, HD-Zip, NF-YA, NF-YB), motif-3 (ARF, MYB), and motif-5 (HSF). Similarly, cytosine was the second most conserved sequence in motif-1 (AP2, ERF, and NF-YB), motif-2 (NAC), motif-4 (bZIP), and motif-5 (WRKY). The conservation and differences in motif sequences between proteins represent functional similarity and diversity in terms of various aspects of biological functions (Puranik et al., 2012). Conserved motifs can provide evidence for further classification into different subgroups as it is possible that proteins within a subgroup, which share identical motifs, are likely to exhibit similar functions. Motifs information can be used to develop resistance genes and markers, perform clustering (Broin et al., 2015), assess the gene expression (Jensen et al., 2005; Huber and Bulyk, 2006), discover the homology relations, classify the families (Blekas et al., 2005; Jensen et al., 2005), discover the sub-families in large protein families (Leonardi and Galves, 2005) and identify new signaling pathways (Ma et al., 2013).
Figure 4. Conserved motif sequences in each of the 15 gene families of maize. The protein sequence logo for highly conserved motifs in each transcription factor is shown above.
Figure 5. Conserved motif sequences in each of the 15 gene families of maize. The nucleotide sequence logo for highly conserved motifs in each transcription factor is shown above.
cis-Regulatory Elements
The CREs in promoter region are the short motifs on which the TF binds to regulate the gene expression. To understand the gene regulation, the knowledge of cis-regulatory elements is indispensable (Mochida et al., 2011). To further investigate the transcriptional regulation and the potential functions of these drought-related TFs in maize, cis-regulatory motifs, a 1 kb upstream promoter region was selected in all TFs using Plant Promoter Analysis Navigator (PlantPAN). A total of 42 types of CREs were identified, of which 28 were common to all TFs (Table 1). These CREs help in determining transcriptional activity by specifying its time and location. The most abundant CRE was DOFCOREZM, with 1,168 duplications, followed by POLASIG3 andDRE1COREZMRAB17, with 1,106 and 1,059 duplications, respectively, whereas ANAEROBICCISZMGAPC4 was the least abundant, with only 134 duplications (Supplementary Table S4). The functions of CREs include responding to dehydration, sugar signaling, responding to light, fermentation and response to phyto-hormones. Of 28 CREs, only 4 were specific to drought stress. ABREAZMRAB28 (ABRE A), found at −148 to −139, and ABREBZMRAB28 (ABRE B) found at −105 to −96, were involved in ABA-responsive function in maize embryos via the rab28 gene (Busk and Pagès, 1997). Other drought-related CREs, such as DRE1COREZMRAB17 (DRE1), DRE/CRT (dehydration responsive element/C-repeat), and DRECRTCOREAT were also found (Busk et al., 1997; Dubouzet et al., 2003). The identification of these CREs could help in understanding the diversity, distribution and evolutionary relationships of drought-related genes in maize (Wittkopp and Kalay, 2011; Hernandez-Garcia and Finer, 2014; Lemmon et al., 2014).
Table 1. Typical cis-regulatory elements found in maize with their consensus sequences and functions.
Interaction Among cis-Regulatory Elements and TF Assigning Drought Tolerance in Plants
Plants respond to drought stress at physiological and molecular levels through complex processes, mediated by TF interactions (Figure 6). The drought-tolerant mechanisms are regulated through TFs by binding to different cis-acting elements in maize (Figure 6). DREB, AREB/ABF, MYB, WRKY, and ERF bind to cis-regulatory elements such as DRE/CRT, ABRE, MYBRS, W-boxes, and GCC-boxes, respectively that are present in the promoter region of drought-responsive genes to regulate the gene expression (Nakashima et al., 2007). For example- MYB by binding to MYBR regulates different drought tolerant mechanisms including stomatal closure, ROS detoxification, and scavenging. ARF by binding to AuxRE regulates plant growth and development and reduced lateral branching. NAC was involved in maintaining cell turgor pressure, ROS homeostasis and plant development. NF-YB was involved in controlling higher photosynthetic rate by binding with CBF cis-element. While AP2, ARF, GRAS, GATA, C2H2, WRKY, MYB, bHLH, and bZIP were found to regulate plant growth and development via ABA/GA signaling.
Figure 6. The genome-TF interaction in assigning drought tolerance in plants. The TFs binds to cis-acting elements to modulate the physiology of plants under drought stress. A light green color box indicates TFs, and the sky-blue color indicates their cis-acting elements.
The interaction of TFs in ABA-dependent and ABA-independent pathways are known to regulate the plant transcriptional response under drought stress by activating stress-responsive genes, resulting in stomatal closure to check excessive water loss by transpiration (Nakashima et al., 2014). The involvement of AtMYB60 and AtMYB61 in stomatal regulation was clearly demonstrated under drought stress (Cominelli et al., 2005; Liang et al., 2005). In rice, SNAC1 and DST (a C2H2-type Zinc finger-containing protein) were involved in stomatal opening/closure (Hu et al., 2006; Huang et al., 2009). The interaction between AREB/ABF and NAC was previously reported under drought stress. ATAF1, a SNAC TF, binds to the promoter of NCED3, which directly controls ABA hormone levels, suggested the role of SNAC in regulating ABA-dependent gene expression (Jensen et al., 2013). In response to drought and osmotic stress, it was reported that ANAC096 cooperated with ABF2 and ABF4 factors (Xu et al., 2013).
Drought stress increases the water-use-efficiency (WUE) through the promotion of deeper root system and inhibition of lateral root development to maximize the water absorption from deeper soil layers (Overvoorde et al., 2010). During drought, several NAC genes regulate the root morphology to enhance the WUE in plants (Jeong et al., 2010). Detoxification of ROS formed during drought stress is important to adjust the physiology of the plants to cope up with the stress. ROS accumulation is dependent on the equilibrium between ROS scavenging and ROS production (Mittler et al., 2004). Antioxidants, such as ascorbic acid, glutathione, and ROS-scavenging enzymes such as glutathione peroxidase, superoxide dismutase, and ascorbate peroxidase, are required for ROS detoxification (Mittler et al., 2004). In rice, several genes, including NAC (Hu et al., 2006), DST (Huang et al., 2009), MYB (Dai et al., 2007), and ZIP (Xiang et al., 2008), play significant roles in drought tolerance through ROS scavenging. C2H2-type zinc finger proteins play vital role in the synthesis of osmoprotectants and detoxification of ROS (Davletova, 2005). WRKY genes were suggested to be involved in ROS detoxification and ROS scavenging. Crosstalk between ROS signaling, ABA-signaling, and other stress-induced redox and metabolite signaling has also been observed. In wheat, to prevent oxidative damage in plant cells, WRKY44 decreases ROS levels (Wang et al., 2015). Cross-talks between AREB with DREB led to the complex interaction of TF families (Narusaka et al., 2003; Lee et al., 2010; Kim et al., 2011). In Arabidopsis, ANAC096 interacts with AREB to provide the drought tolerance (Xu et al., 2013) and interaction between DREB and AP2/ERF has also been reported (Cheng et al., 2013). OsMyb4 integrated with ABA and bZIP that directly regulates the ABA-responsive genes (Nakashima et al., 2009). It was also reported that NF-YB and NF-YC form a dimer and interacts with NF-YA in the nucleus and plays prominent role in root development (He et al., 2016). Understanding the inter-connection among the drought tolerant mechanisms through TFs regulation would help to breed drought-tolerant maize in the changing climatic scenario in a more efficient way. The information could be also useful in short-listing the best combination of genes that work in tandem under drought stress condition for further exploitation.
Expression Analysis of Drought-Responsive TFs Genes Under Drought
To understand the expression pattern of TF families in response to drought stress, RNA-Seq data available in the public domain were used. Several reports have discussed the role of drought-responsive TFs genes in the modulation of drought tolerance through differential expression patterns in tissue-specific and generic manners (Gahlaut et al., 2016). To investigate the expression pattern of these TFs families under drought stress, the TF genes were considered from the experiments of Liu et al. (2015) and Thatcher et al. (2016). A total of 269 genes belonging to 14 TF families with fold change of >2 was selected for further analysis (Figure 7A). Results explained that the genes of the ARF family did not show any differential expression under drought stress conditions. The maximum number of differentially expressed genes was in ERF and MYB families (56 each) and the least was in AP2, GATA, and NF-YA (2 each). Furthermore, among tissues, roots showed maximum number of differentially expressed genes (DEGs) (109), followed by leaves (92), shoot (72), ear (22) and the lowest was in the tassel (18). No common TF genes were differentially expressed in the tissues, indicating the tissue specificity in TF gene expression. Genes were over-expressed in response to drought stress in the majority of the TFs families. It was noted that the roots regulate drought stress by plant hormones and ABA pathways. An early response to drought stress induces ABA in roots, which further transmits signals from roots to aerial parts through a cascade of reactions (Janiak et al., 2016). In roots, drought stress induces up-regulation of TF genes was reported in several plant species, including rice, maize, and wheat (Janiak et al., 2016; Opitz et al., 2016; Zhang et al., 2016). Shiriga et al. (2014) reported an up-regulation of NAC genes (GRMZM2G336533 and GRMZM2G347043) in the tolerant genotype (HKI577) under drought stress compared to the sensitive genotype (HKIPC3-3). It was demonstrated that the down-regulation of ATHB7, a member of HD-Zip, was associated with drought-induced ABA signal transduction (Gao et al., 2015). Furthermore, several TF genes did not show any differential expression pattern in the selected tissues due to an environmental condition, genetic background effect and/or growth stages of the plant (Wang et al., 2017). Among bZIP genes, GRMZM2G140355 showed leaf-specific down-regulation under drought stress which was also confirmed by Jiang et al. (2012).
Figure 7. Heatmap and interaction network of 269 drought-responsive genes. (A) Overview of differentially expressed TF genes in response to drought stress in various tissues with log2-fold changes of 2 or more and with a p-value of <0.05. A total of 269 DEGs from 14 of 15 TF families were represented in different tissues. (B) Drought-responsive TF network of 269 genes. A total of seven networks were identified. The color of the nodes indicates the category of TFs.
The role of TF families under drought stress has been documented in several crops. For example- the over-expression of DREB1/CBF improved drought tolerance in tomato, chrysanthemum and potato (Iwaki et al., 2013), soybean (de Paiva Rolla et al., 2014), rice (Nakashima et al., 2014; Paul et al., 2015), tobacco (Phuong et al., 2015), sugarcane (Augustine et al., 2014), groundnut (Tiwari et al., 2015), peanut (Bhatnagar-Mathur et al., 2014), and wheat (Shavrukov et al., 2016). Majority of the DREB genes were involved in ABA-independent stress responses whereas ABA-dependent responses was also reported (Sakuma et al., 2006; Nakashima et al., 2009). In rice, regulation of downstream gene OsPP18 by SNAC1 conferred drought tolerance via ABA-independent pathway (You et al., 2014). Over-expression of ANACO19, ANACO55, and ANACO72 also explained drought tolerance in Arabidopsis. In barley, HvSNAC1 not only imparted drought tolerance but also provided resistance against leaf spot (Al Abdallat et al., 2014; McGrann et al., 2015). AtMYB60 and AtMYB96 were involved in drought response and regulate stomata and disease resistance via ABA-dependent pathways (Cominelli et al., 2005; Seo et al., 2009). In rice, OsWRKY45 was up-regulated under drought stress treatment and involved in ABA synthesis that induces a signaling cascade resulting in enhanced drought tolerance (Qiu and Yu, 2009; Supplementary Table S5).
Interaction of Drought-Responsive TFs in Assigning Drought Tolerance
Protein-protein interaction is an important component of gene regulation and is known to mediate most of the cellular processes (Ding et al., 2014). The interactions among different TFs are the major sources of diversity in the gene expression pattern. The physical interactions among multiple TFs may mediate the specific set of conditions under which gene expression will be regulated (Reményi et al., 2004; Hussain et al., 2011). Therefore, understanding the interactions of various drought-responsive TF proteins provides a deeper understanding of molecular mechanisms and pathways operating in drought tolerance. A set of 269 genes showing variable expression levels under drought stress in different tissues of maize was selected for interaction analysis. The results revealed the presence of seven networks with 269 nodes and 212 edges with a clustering coefficient of 0.79 (Figure 7B). Network I showed maximum number of interfamily interactions among the TF families, including bZIP, GATA, HSF, MYB, NF-YB, and NF-YA.
In Network I, bZIP (GRMZM2G140355) showed maximum interactions (54), including 11 HSF and 43 MYB. From the phylogenetic tree (Figure 2), the MYB genes are grouped together with bZIP genes. Similarly, the HSF genes were grouped in a single cluster. The grouping of genes of TF families suggested that they could be involved in similar functions under drought stress. MYB, HSF, and bZIP explained the higher level of structural relationship since more than 95% of the cis-elements were same. From an evolutionary perspective, among different TF families, it has been noticed that some TF families have higher rate of expansion including bZIP, BHLH and MYB over others (Shiu et al., 2005; Feller et al., 2006; Riechmann et al., 2016). The interactions among the bZIP and MYB in Arabidopsis are known to provide light responsiveness in plants and induce stress signaling pathways (Hartmann et al., 2005). bZIP and MYB, along with tobacco bZIP transcription activator (TAF), play direct roles in drought tolerance through ABA signaling pathway (Golldack et al., 2011). OsMyb4 was integrated with ABA and other signaling pathways through bZIP that directly regulates the ABA-responsive genes (Nakashima et al., 2009). The GATA genes, AC194965.4_FG004, and GRMZM2G113098 interacted with MYB and NF-YB families (IN WHICH CROP OR WHAT PURPOSE/ROLE). During seed germination, a repressor of gal-3 (RGA) and RGA-LIKE2 (RGL2) interacts with NF-YB proteins to regulate GA, ABA, and GA signaling pathways (Hou et al., 2014; Liu et al., 2016). NF-YA (GRMZM2G026157, GRMZM2G038303) interacted with six NF-YB proteins (GRMZM2G180947, GRMZM2G444073, GRMZM2G473152, GRMZM2G180947, GRMZM2G444073, and GRMZM2G473152) (IN WHICH CROP OR WHAT PURPOSE/ROLE). The interaction of NF-YA and NF-YB families could be attributed to the high level of similarity of cis-elements (98%) between them. Several NF-YA and NF-YB have been found to play significant role in drought tolerance (Nelson et al., 2007; Li et al., 2008; Hackenberg et al., 2012). The NF-YA (GRMZM2G026157) and its interaction with other co-factors play a prominent role in the development of roots under stress (He et al., 2016). It has also been reported that NF-YB and NF-YC form a dimer and further by interacting with NF-YA, a trimer is formed in the nucleus. In the promoter of its target genes, this trimer binds to cis-elements called CCAAT-box and further regulates the downstream stress-related genes (Frontini et al., 2004; Kahle et al., 2005). It has been reported that TaNF-YA10–1 by interacting with NF-YA/NF-YB activates target drought-responsive genes (Ma et al., 2015).
Utility of Drought-Related TFs in Stress Breeding
TFs are the potential regulators of various genes that help plants to effectively cope up with drought stress. To unravel the molecular mechanisms in response to drought stress, it is important to mine candidate genes, pathways, and the interaction of genes in different pathways. Genetic manipulation of combination of drought stress-responsive TF genes stands to be a powerful approach for improving stress tolerance than addressing each functional gene individually. DREB, MYB, and bZIP are the important TFs which have been already used in several molecular breeding programmes (Foolad et al., 2003; Fleury et al., 2010). Genes involved in ABA biosynthesis, catabolism, and signaling network represent interesting candidate genes for breeding crops with improved drought tolerance. Different candidate genes that are involved in the detoxification process have been reported to benefit plant growth as well as survival upon abiotic stress and therefore provide an excellent resource for improving drought tolerance in crop breeding.
The major emphasize of the experiment is laid on detailed characterization of drought-related TFs and their structural and functional importance in assigning the drought tolerance. It should be noted that only a low level of phenotypic variation (for example 10–30%) could be achieved through transfer of a major gene or QTL since drought tolerance is a complex trait controlled by many genes (Thirunavukkarasu et al., 2014, 2017). On the other hand, by selecting a combination of genes that are involved in related pathways would provide enhanced stress tolerance. Our experiment provided comprehensive information of TF genes that are related to stress tolerance and their interactions. Presently, the breeding strategies for quantitative traits such as drought are re-oriented from a few genes (through marker-based selection) to genome-wide approaches (through genomic selection) to capture the total phenotypic variation. The identified set of TFs could be targeted as candidates to breed drought tolerant maize genotypes through different selection approaches.
Conclusion
TFs play a key role in the regulation of drought tolerance in crop plants. An effort has been made to characterize the major drought-related TF families encoded by 1,436 genes in maize. The emphasis of the study was to understand the structure, function, and evolution of the drought-related TF families and their utility in maize drought breeding. An extensive in-silico analysis was performed to understand the distribution, structure, expression pattern, interactions, and function in maize. Our results linked the structural characteristics of different TFs with the functional aspects in drought-stress tolerance. The working hypothesis revealed the interaction of TFs with various molecular functions was crucial in assigning drought tolerance in maize. The identified set of TFs could be targeted to breed drought-tolerant maize genotypes through selection approaches, such as genomic selection, to achieve maximum phenotypic expression. The present study thus attempts to provide a novel foundation for such comparative analysis of TF families and the findings would accelerate the understanding of drought tolerance mechanisms in maize in addition to their practical utility in drought breeding programs.
Author Contributions
NT and SM conceived and designed the experiments. SM, PB, MM, AR, and PJ analyzed the data. NT, SM, MM, PB, and PD drafted the manuscript. All authors read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer, SK, and handling Editor declared their shared affiliation.
Acknowledgments
We thank Computational Biology and Agricultural Bioinformatics [Agril.Edn.14 (44)/2014- A&P] and the ICAR Network Projects on Transgenics in Crop Plants (Maize Functional Genomics Component: 21-22) for funding the study. The funders had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00177/full#supplementary-material
References
Al Abdallat, A. M., Ayad, J. Y., Abu Elenein, J. M., Al Ajlouni, Z., and Harwood, W. A. (2014). Overexpression of the transcription factor HvSNAC1 improves drought tolerance in barley (Hordeum vulgare L.). Mol. Breed. 33, 401–414. doi: 10.1007/s11032-013-9958-1
Allan, A. C., Hellens, R. P., and Laing, W. A. (2008). MYB transcription factors that colour our fruit. Trends Plant Sci. 13, 99–102. doi: 10.1016/j.tplants.2007.11.012
Augustine, S. M., Ashwin Narayan, J., Syamaladevi, D. P., Appunu, C., Chakravarthi, M., Ravichandran, V., et al. (2014). Overexpression of EaDREB2 and pyramiding of EaDREB2 with the pea DNA helicase gene (PDH45) enhance drought and salinity tolerance in sugarcane (Saccharum spp. hybrid). Plant Cell Rep. 34, 247–263. doi: 10.1007/s00299-014-1704-6
Baloglu, M. C., Eldem, V., Hajyzadeh, M., and Unver, T. (2014). Genome-wide analysis of the bZIP transcription factors in cucumber. PLoS ONE 9:e96014. doi: 10.1371/journal.pone.0096014
Barbosa, E. G. G., Leite, J. P., Marin, S. R. R., Marinho, J. P., de Fátima Corrêa Carvalho, J., Fuganti-Pagliarini, R., et al. (2013). Overexpression of the ABA-dependent AREB1 transcription factor from Arabidopsis thaliana improves soybean tolerance to water deficit. Plant Mol. Biol. Rep. 31, 719–730. doi: 10.1007/s11105-012-0541-4
Bennetzen, J. L., and Hake, S. (2009). Handbook of Maize: Genetics and Genomics. New York, NY: Springer-Verlag.
Bhatnagar-Mathur, P., Rao, J. S., Vadez, V., Dumbala, S. R., Rathore, A., Yamaguchi-Shinozaki, K., et al. (2014). Transgenic peanut overexpressing the DREB1A transcription factor has higher yields under drought stress. Mol. Breed. 33, 327–340. doi: 10.1007/s11032-013-9952-7
Björklund, A. K., Ekman, D., Light, S., Frey-Skött, J., and Elofsson, A. (2005). Domain rearrangements in protein evolution. J. Mol. Biol. 353, 911–923. doi: 10.1016/j.jmb.2005.08.067
Blanc, G., and Wolfe, K. (2014). Functional divergence of duplicated genes formed by polyploidy during arabidopsis evolution. Plant Cell 16, 1679–1691. doi: 10.1105/tpc.021410
Blekas, K., Fotiadis, D. I., and Likas, A. (2005). Motif-based protein sequence classification using neural networks. J. Comput. Biol. 12, 64–82. doi: 10.1089/cmb.2005.12.64
Brini, F., Hanin, M., Mezghani, I., Berkowitz, G. A., and Masmoudi, K. (2007). Overexpression of wheat Na+/H+ antiporter TNHX1 and H+-pyrophosphatase TVP1 improve salt- and drought-stress tolerance in Arabidopsis thaliana plants. J. Exp. Bot. 58, 301–308. doi: 10.1093/jxb/erl251
Broin, P. Ó., Smith, T. J., and Golden, A. A. J. (2015). Alignment-free clustering of transcription factor binding motifs using a genetic-k-medoids approach. BMC Bioinformatics 16:2. doi: 10.1186/s12859-015-0450-2
Busk, P. K., Jensen, A. B., and Pagès, M. (1997). Regulatory elements in vivo in the promoter of the abscisic acid responsive gene rab17 from maize. Plant J. 11, 1285–1295. doi: 10.1046/j.1365-313X.1997.11061285.x
Busk, P. K., and Pagès, M. (1997). Protein binding to the abscisic acid-responsive element is independent of VIVIPAROUS1 in vivo. Plant Cell 9, 2261–2270. doi: 10.1105/tpc.9.12.2261
Cao, J., and Shi, F. (2012). Evolution of the RALF gene family in plants: gene duplication and selection patterns. Evol. Bioinform. 2012, 271–292. doi: 10.4137/EBO.S9652
Cheng, M.-C., Liao, P.-M., Kuo, W.-W., and Lin, T.-P. (2013). The Arabidopsis ETHYLENE RESPONSE FACTOR1 regulates abiotic stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals. Plant Physiol. 162, 1566–1582. doi: 10.1104/pp.113.221911
Ciftci-Yilmaz, S., and Mittler, R. (2008). The zinc finger network of plants. Cell. Mol. Life Sci. 65, 1150–1160. doi: 10.1007/s00018-007-7473-4
Cominelli, E., Galbiati, M., Vavasseur, A., Conti, L., Sala, T., Vuylsteke, M., et al. (2005). A guard-cell-specific MYB transcription factor regulates stomatal movements and plant drought tolerance. Curr. Biol. 15, 1196–1200. doi: 10.1016/j.cub.2005.05.048
Cominelli, E., and Tonelli, C. (2009). A new role for plant R2R3-MYB transcription factors in cell cycle regulation. Cell Res. 19, 1231–1232. doi: 10.1038/cr.2009.123
Dai, X., Xu, Y., Ma, Q., Xu, W., Wang, T., Xue, Y., et al. (2007). Overexpression of an R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic Arabidopsis. Plant Physiol. 143, 1739–1751. doi: 10.1104/pp.106.094532
Davletova, S. (2005). Cytosolic ascorbate peroxidase 1 is a central component of the reactive oxygen gene network of Arabidopsis. Plant Cell 17, 268–281. doi: 10.1105/tpc.104.026971
de Paiva Rolla, A. A., de Fátima Corrêa Carvalho, J., Fuganti-Pagliarini, R., Engels, C., do Rio, A., Marin, S. R. R., et al. (2014). Phenotyping soybean plants transformed with rd29A:AtDREB1A for drought tolerance in the greenhouse and field. Transgenic Res. 23, 75–87. doi: 10.1007/s11248-013-9723-6
Delseny, M., Bies-etheve, N., Carles, C., Hull, G., Vicient, C., Raynal, M., et al. (2001). Late Embryogenesis Abundant ( LEA ) protein gene regulation during Arabidopsis seed maturation. Plant Physiol. 158, 419–427. doi: 10.1078/0176-1617-00353
Ding, Y. D., Cheng, J. W., Guo, J., Chen, D., Li, S., Xu, Q., et al. (2014). Prediction and functional analysis of the sweet orange protein-protein interaction network. BMC Plant Biol. 14:213. doi: 10.1186/s12870-014-0213-7
Du, H., Wang, Y., Bin Xie, Y., Liang, Z., Jiang, S. J., Zhang, S. S., et al. (2013). Genome-wide identification and evolutionary and expression analyses of MYB-related genes in land plants. DNA Res. 20, 437–448. doi: 10.1093/dnares/dst021
Dubouzet, J. G., Sakuma, Y., Ito, Y., Kasuga, M., Dubouzet, E. G., Miura, S., et al. (2003). OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought, high salt and cold responsive gene expression. Plant J. 33, 751–763. doi: 10.1046/j.1365-313X.2003.01661.x
Edger, P. P., and Pires, J. C. (2009). Gene and genome duplications: the impact of dosage-sensitivity on the fate of nuclear genes. Chromosom. Res. 17, 699–717. doi: 10.1007/s10577-009-9055-9
Engels, C., Fuganti-Pagliarini, R., Marin, S. R. R., Marcelino-Guimarães, F. C., Oliveira, M. C. N., Kanamori, N., et al. (2013). Introduction of the rd29A:AtDREB2A ca gene into soybean (Glycine max L. Merril) and its molecular characterization in leaves and roots during dehydration. Genet. Mol. Biol. 36, 556–565. doi: 10.1590/S1415-47572013000400015
Feller, A., Hernandez, J. M., and Grotewold, E. (2006). An ACT-like domain participates in the dimerization of several plant basic-helix-loop-helix transcription factors. J. Biol. Chem. 281, 28964–28974. doi: 10.1074/jbc.M603262200
Feller, A., MacHemer, K., Braun, E. L., and Grotewold, E. (2011). Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant J. 66, 94–116. doi: 10.1111/j.1365-313X.2010.04459.x
Fleury, D., Jefferies, S., Kuchel, H., and Langridge, P. (2010). Genetic and genomic tools to improve drought tolerance in wheat. J. Exp. Bot. 61, 3211–3222. doi: 10.1093/jxb/erq152
Foolad, M. R., Zhang, L. P., and Subbiah, P. (2003). Genetics of drought tolerance during seed germination in tomato: inheritance and QTL mapping. Genome 46, 536–545. doi: 10.1139/g03-035
Frontini, M., Imbriano, C., Manni, I., and Mantovani, R. (2004). Cell cycle regulation of NF-YC nuclear localization. Cell CycleCell Cycle 217, 217–222. doi: 10.4161/cc.3.2.654
Gahlaut, V., Jaiswal, V., Kumar, A., and Gupta, P. K. (2016). Transcription factors involved in drought tolerance and their possible role in developing drought tolerant cultivars with emphasis on wheat (Triticum aestivum L.). Theor. Appl. Genet. 129, 2019–2042. doi: 10.1007/s00122-016-2794-z
Gao, F., Wang, J., Wei, S., Li, Z., Wang, N., Li, H., et al. (2015). Transcriptomic analysis of drought stress responses in Ammopiptanthus mongolicus leaves using the RNA-seq technique. PLoS ONE 10:e124382. doi: 10.1371/journal.pone.0124382
Golldack, D., Lüking, I., and Yang, O. (2011). Plant tolerance to drought and salinity: stress regulating transcription factors and their functional significance in the cellular transcriptional network. Plant Cell Rep. 30, 1383–1391. doi: 10.1007/s00299-011-1068-0
Hackenberg, D., Wu, Y., Voigt, A., Adams, R., Schramm, P., and Grimm, B. (2012). Studies on differential nuclear translocation mechanism and assembly of the three subunits of the Arabidopsis thaliana transcription factor NF-Y. Mol. Plant 5, 876–888. doi: 10.1093/mp/ssr107
Hall, A. J., Lemcoff, J. H., and Trapani, N. (1981). Water stress before and during flowering in maize and its effects on yield, its components, and their determinants. Maydica 26, 19–38.
Hartmann, U., Sagasser, M., Mehrtens, F., Stracke, R., and Weisshaar, B. (2005). Differential combinatorial interactions of cis-acting elements recognized by R2R3-MYB, BZIP, and BHLH factors control light-responsive and tissue-specific activation of phenylpropanoid biosynthesis genes. Plant Mol. Biol. 57, 155–171. doi: 10.1007/s11103-004-6910-0
He, S., Zou, M., and Guo, B. (2011). The roles and evolutionary patterns of intronless genes in deuterostomes. Comp. Funct. Genomics 2011:680673. doi: 10.1155/2011/680673
He, X., Ma, H., Zhao, X., Nie, S., Li, Y., Zhang, Z., et al. (2016). Comparative RNA-Seq analysis reveals that regulatory network of maize root development controls the expression of Genes in response to N stress. PLoS ONE 11:e0151697. doi: 10.1371/journal.pone.0151697
Hernandez-Garcia, C. M., and Finer, J. J. (2014). Identification and validation of promoters and cis-acting regulatory elements. Plant Sci. 217–218, 109–119. doi: 10.1016/j.plantsci.2013.12.007
Hou, X., Zhou, J., Liu, C., Liu, L., Shen, L., and Yu, H. (2014). Nuclear factor Y-mediated H3K27me3 demethylation of the SOC1 locus orchestrates flowering responses of Arabidopsis. Nat. Commun. 5, 1–14. doi: 10.1038/ncomms5601
Hu, H., Dai, M., Yao, J., Xiao, B., Li, X., Zhang, Q., et al. (2006). Over expressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc. Natl. Acad. Sci. U.S.A. 103, 12987–12992. doi: 10.1073/pnas.0604882103
Huang, X., Chao, D., Gao, J., Zhu, M., and Shi, M. (2009). A previously unknown zinc finger protein, DST, regulates drought and salt tolerance in rice. Genes Dev. 23, 1805–1817. doi: 10.1101/gad.1812409
Huber, B. R., and Bulyk, M. L. (2006). Meta-analysis discovery of tissue-specific DNA sequence motifs from mammalian gene expression data. BMC Bioinformatics 7:229. doi: 10.1186/1471-2105-7-229
Hussain, S. S., Kayani, M. A., and Amjad, M. (2011). Transcription factors as tools to engineer enhanced drought stress tolerance in plants. Biotechnol. Prog. 27, 297–306. doi: 10.1002/btpr.514
Iwaki, T., Guo, L., Ryals, J. A., Yasuda, S., Shimazaki, T., Kikuchi, A., et al. (2013). Metabolic profiling of transgenic potato tubers expressing Arabidopsis dehydration response element-binding protein 1A (DREB1A). J. Agric. Food Chem. 61, 893–900. doi: 10.1021/jf304071n
Janiak, A., Kwaśniewski, M., and Szarejko, I. (2016). Gene expression regulation in roots under drought. J. Exp. Bot. 67, 1003–1014. doi: 10.1093/jxb/erv512
Jensen, M. K., Lindemose, S., de Masi, F., Reimer, J. J., Nielsen, M., Perera, V., et al. (2013). ATAF1 transcription factor directly regulates abscisic acid biosynthetic gene NCED3 in Arabidopsis thaliana. FEBS Open Biol. 3, 321–327. doi: 10.1016/j.fob.2013.07.006
Jensen, S. T., Shen, L., and Liu, J. S. (2005). Combining phylogenetic motif discovery and motif clustering to predict co-regulated genes. Bioinformatics 21, 3832–3839. doi: 10.1093/bioinformatics/bti628
Jeong, J. S., Kim, Y. S., Baek, K. H., Jung, H., Ha, S.-H., Do Choi, Y., et al. (2010). Root-specific expression of OsNAC10 improves drought tolerance and grain yield in rice under field drought conditions. Plant Physiol. 153, 185–197. doi: 10.1104/pp.110.154773
Jiang, Y., Zeng, B., Zhao, H., Zhang, M., Xie, S., and Lai, J. (2012). Genome-wide transcription factor gene prediction and their expressional tissue-specificities in maize. J. Integr. Plant Biol. 54, 616–630. doi: 10.1111/j.1744-7909.2012.01149.x
Jin, H., and Martin, C. (1999). Multifunctionality and diversity within the plant MYB -gene family. Plant Mol. Bio. 41, 577–585.
Jogaiah, S., Govind, S. R., and Tran, L.-S. P. (2013). Systems biology-based approaches toward understanding drought tolerance in food crops. Crit. Rev. Biotechnol. 33, 23–39. doi: 10.3109/07388551.2012.659174
Joshi, R., Wani, S. H., Singh, B., Bohra, A., Dar, Z. A., Lone, A. A., et al. (2016). Transcription factors and plants response to drought stress: current understanding and future directions. Front. Plant Sci. 7:1029. doi: 10.3389/fpls.2016.01029
Kahle, J., Baake, M., Doenecke, D., and Albig, W. (2005). Subunits of the heterotrimeric transcription factor NF-Y are imported into the nucleus by distinct pathways involving importin and importin 13. Mol. Cell. Biol. 25, 5339–5354. doi: 10.1128/MCB.25.13.5339-5354.2005
Katiyar, A., Smita, S., Lenka, S. K., Rajwanshi, R., Chinnusamy, V., and Bansal, K. C. (2012). Genome-wide classification and expression analysis of MYB transcription factor families in rice and Arabidopsis. BMC Genomics 13:544. doi: 10.1186/1471-2164-13-544
Kim, J. S., Mizoi, J., Yoshida, T., Fujita, Y., Nakajima, J., Ohori, T., et al. (2011). An ABRE promoter sequence is involved in osmotic stress-responsive expression of the DREB2A gene, which encodes a transcription factor regulating drought-inducible genes in Arabidopsis. Plant Cell Physiol. 52, 2136–2146. doi: 10.1093/pcp/pcr143
Leach, K. A., Hejlek, L. G., Hearne, L. B., Nguyen, H. T., Sharp, R. E., and Davis, G. L. (2011). Primary root elongation rate and abscisic acid levels of maize in response to water stress. Crop Sci. 51, 157–172. doi: 10.2135/cropsci2009.12.0708
Lee, S. J., Kang, J. Y., Park, H.-J., Kim, M. D., Bae, M. S., Choi, H. I., et al. (2010). DREB2C interacts with ABF2, a bZIP protein regulating abscisic acid-responsive gene expression, and its overexpression affects abscisic acid sensitivity. Plant Physiol. 153, 716–727. doi: 10.1104/pp.110.154617
Lemmon, Z. H., Bukowski, R., Sun, Q., and Doebley, J. F. (2014). The role of cis regulatory evolution in maize domestication. PLoS Genet. 10:e004745. doi: 10.1371/journal.pgen.1004745
Leonardi, F., and Galves, A. (2005). “Sequence motif identification and protein family classification using probabilistic trees,” in Advances in Bioinformatics and Computational Biology. BSB 2005. Lecture Notes in Computer Science, eds J. C. Setubal and S. Verjovski-Almeida. (Berlin; Heidelberg: Springer), 190–193.
Li, W.-X., Oono, Y., Zhu, J., He, X.-J., Wu, J.-M., Iida, K., et al. (2008). The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell 20, 2238–2251. doi: 10.1105/tpc.108.059444
Liang, Y. K., Dubos, C., Dodd, I. C., Holroyd, G. H., Hetherington, A. M., and Campbell, M. M. (2005). AtMYB61, an R2R3-MYB transcription factor controlling stomatal aperture in Arabidopsis thaliana. Curr. Biol. 15, 1201–1206. doi: 10.1016/j.cub.2005.06.041
Liu, X., Hu, P., Huang, M., Tang, Y., Li, Y., Li, L., et al. (2016). The NF-YC–RGL2 module integrates GA and ABA signalling to regulate seed germination in Arabidopsis. Nat. Commun. 7:12768. doi: 10.1038/ncomms12768
Liu, Y., Jiang, H. Y., Chen, W. J., Qian, Y. X., Ma, Q., Cheng, B. J., et al. (2011). Genome-wide analysis of the auxin response factor (ARF) gene family in maize (Zea mays). Plant Growth Regul. 63, 225–234. doi: 10.1007/s10725-010-9519-0
Liu, Y., Zhou, M., Gao, Z., Ren, W., Yang, F., He, H., et al. (2015). RNA-seq analysis reveals MAPKKK family members related to drought tolerance in maize. PLoS ONE 10:e2818302. doi: 10.1145/2818302
Lu, T., Yang, Y., Yao, B., Liu, S., Zhou, Y., and Zhang, C. (2012). Template-based structure prediction and classification of transcription factors in Arabidopsis thaliana. Protein Sci. 21, 828–838. doi: 10.1002/pro.2066
Ma, S., Shah, S., Bohnert, H. J., Snyder, M., and Dinesh-Kumar, S. P. (2013). Incorporating motif analysis into gene co-expression networks reveals novel modular expression pattern and new signaling pathways. PLoS Genet. 9:e1003840. doi: 10.1371/journal.pgen.1003840
Ma, X., Li, C., and Wang, M. (2015). Wheat NF-YA10 functions independently in salinity and drought stress. Bioengineered 6, 245–247. doi: 10.1080/21655979.2015.1054085
Mao, H., Yu, L., Li, Z., Liu, H., and Han, R. (2016). Molecular evolution and gene expression differences within the HD-Zip transcription factor family of Zea mays L. Genetica 144, 243–257. doi: 10.1007/s10709-016-9896-z
McGrann, G. R., Steed, A., Burt, C., Goddard, R., Lachaux, C., Bansal, A., et al. (2015). Contribution of the drought tolerance-related Stress-responsive NAC1 transcription factor to resistance of barley to Ramularia leaf spot. Mol. Plant Pathol. 16, 201–209. doi: 10.1111/mpp.12173
Mittal, S., Arora, K., Rao, A. R., Mallikarjuna, M. G., Gupta, H. S., and Nepolean, T. (2017). Genomic selection for drought tolerance using genome-wide SNPs in maize. Front. Plant Sci. 8:550. doi: 10.3389/fpls.2017.00550
Mittler, R., Vanderauwera, S., Gollery, M., and Van Breusegem, F. (2004). Reactive oxygen gene network of plants. Trends Plant Sci. 9, 490–498. doi: 10.1016/j.tplants.2004.08.009
Mochida, K., Yoshida, T., Sakurai, T., Yamaguchi-Shinozaki, K., Shinozaki, K., and Tran, L. S. (2010). Genome-wide analysis of two-component systems and prediction of stress-responsive two-component system members in soybean. DNA Res. 17, 303–324. doi: 10.1093/dnares/dsq021
Mochida, K., Yoshida, T., Sakurai, T., Yamaguchi-Shinozaki, K., Shinozaki, K., and Tran, L. S. P. (2011). In silico analysis of transcription factor repertoires and prediction of stress-responsive transcription factors from six major gramineae plants. DNA Res. 18, 321–332. doi: 10.1093/dnares/dsr019
Muthamilarasan, M., Bonthala, V. S., Khandelwal, R., Jaishankar, J., Shweta, S., Nawaz, K., et al. (2015). Global analysis of WRKY transcription factor superfamily in Setaria identifies potential candidates involved in abiotic stress signaling. Front. Plant Sci. 6:910. doi: 10.3389/fpls.2015.00910
Nakashima, K., Ito, Y., and Yamaguchi-Shinozaki, K. (2009). Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. Plant Physiol. 149, 88–95. doi: 10.1104/pp.108.129791
Nakashima, K., Jan, A., Todaka, D., Maruyama, K., Goto, S., Shinozaki, K., et al. (2014). Comparative functional analysis of six drought-responsive promoters in transgenic rice. Planta 239, 47–60. doi: 10.1007/s00425-013-1960-7
Nakashima, K., Tran, L. S., Van Nguyen, D., Fujita, M., Maruyama, K., Todaka, D., et al. (2007). Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice. Plant J. 51, 617–630. doi: 10.1111/j.1365-313X.2007.03168.x
Narusaka, Y., Nakashima, K., Shinwari, Z. K., Sakuma, Y., Furihata, T., Abe, H., et al. (2003). Interaction between two cis-acting elements, ABRE and DRE, in ABA-dependent expression of Arabidopsis rd29A gene in response to dehydration and high-salinity stresses. Plant J. 34, 137–148. doi: 10.1046/j.1365-313X.2003.01708.x
Nelson, D. E., Repetti, P. P., Adams, T. R., Creelman, R. A., Wu, J., Warner, D. C., et al. (2007). Plant nuclear factor Y (NF-Y) B subunits confer drought tolerance and lead to improved corn yields on water-limited acres. Proc. Natl. Acad. Sci. U.S.A. 104, 16450–16455. doi: 10.1073/pnas.0707193104
Oh, S.-J., Song, S. I., Kim, Y. S., Jang, H.-J., Kim, S. Y., Kim, M., et al. (2005). Arabidopsis CBF3/DREB1A and ABF3 in transgenic rice increased tolerance to abiotic stress without stunting growth. Plant Physiol. 138, 341–351. doi: 10.1104/pp.104.059147
Opitz, N., Marcon, C., Paschold, A., Malik, W. A., Lithio, A., Brandt, R., et al. (2016). Extensive tissue-specific transcriptomic plasticity in maize primary roots upon water deficit. J. Exp. Bot. 67, 1095–1107. doi: 10.1093/jxb/erv453
Overvoorde, P., Fukaki, H., and Beeckman, T. (2010). Auxin control of root development. Cold Spring Harb. Perspect. Biol. 2:a001537. doi: 10.1101/cshperspect.a001537
Paul, S., Gayen, D., Datta, S. K., and Datta, K. (2015). Dissecting root proteome of transgenic rice cultivars unravels metabolic alterations and accumulation of novel stress responsive proteins under drought stress. Plant Sci. 234, 133–143. doi: 10.1016/j.plantsci.2015.02.006
Phuong, N. D., Tuteja, N., Nghia, P. T., and Hoi, P. X. (2015). Identification and characterization of a stress-inducible gene OsNLI-IF enhancing drought tolerance in transgenic tobacco. Curr. Sci. 109, 541–551.
Pillitteri, L. J., Sloan, D. B., Bogenschutz, N. L., and Torii, K. U. (2007). Termination of asymmetric cell division and differentiation of stomata. Nature 445, 501–505. doi: 10.1038/nature05467
Pires, N., and Dolan, L. (2010). Origin and diversification of basic-helix-loop-helix proteins in plants. Mol. Biol. Evol. 27, 862–874. doi: 10.1093/molbev/msp288
Puranik, S., Sahu, P. P., Srivastava, P. S., and Prasad, M. (2012). NAC proteins: regulation and role in stress tolerance. Trends Plant Sci. 17, 369–381. doi: 10.1016/j.tplants.2012.02.004
Qiu, Y., and Yu, D. (2009). Over-expression of the stress-induced OsWRKY45 enhances disease resistance and drought tolerance in Arabidopsis. Environ. Exp. Bot. 65, 35–47. doi: 10.1016/j.envexpbot.2008.07.002
Reményi, A., Schöler, H. R., and Wilmanns, M. (2004). Combinatorial control of gene expression. Nat. Struct. Mol. Biol. 11, 812–815. doi: 10.1038/nsmb820
Riechmann, A. J. L., Heard, J., Martin, G., Reuber, L., Jiang, C., Keddie, J., et al. (2016). Arabidopsis transcription factors : genome-wide comparative analysis among eukaryotes. Science 290, 2105–2110. doi: 10.1126/science.290.5499.2105
Rushton, P. J., Somssich, I. E., Ringler, P., and Shen, Q. J. (2010). WRKY transcription factors. Trends Plant Sci. 15, 247–258. doi: 10.1016/j.tplants.2010.02.006
Sakamoto, H., Araki, T., Meshi, T., and Iwabuchi, M. (2000). Expression of a subset of the Arabidopsis Cys2/His2-type zinc-finger protein gene family under water stress. Gene 248, 23–32. doi: 10.1016/S0378-1119(00)00133-5
Sakuma, Y., Maruyama, K., Osakabe, Y., Qin, F., Seki, M., Shinozaki, K., et al. (2006). Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell 18, 1292–1309. doi: 10.1105/tpc.105.035881
Schmitz, J. F., Zimmer, F., and Bornberg-Bauer, E. (2016). Mechanisms of transcription factor evolution in Metazoa. Nucleic Acids Res. 44, 6287–6297. doi: 10.1093/nar/gkw492
Schnable, P. S., Ware, D., Fulton, R. S., Stein, J. C., Wei, F., Pasternak, S., et al. (2009). The B73 maize genome: complexity, diversity, and dynamics. Science 326, 1112–1115. doi: 10.1126/science.1178534
Seo, P. J., Xiang, F., Qiao, M., Park, J.-Y., Lee, Y. N., Kim, S.-G., et al. (2009). The MYB96 transcription factor mediates abscisic acid signaling during drought stress response in Arabidopsis. Plant Physiol. 151, 275–289. doi: 10.1104/pp.109.144220
Shao, H., Wang, H., and Tang, X. (2015). NAC transcription factors in plant multiple abiotic stress responses: progress and prospects. Front. Plant Sci. 6:902. doi: 10.3389/fpls.2015.00902
Sharma, R., Sahoo, A., Devendran, R., and Jain, M. (2014). Over-expression of a rice tau class glutathione S-transferase gene improves tolerance to salinity and oxidative stresses in arabidopsis. PLoS ONE 9:e92900. doi: 10.1371/journal.pone.0092900
Shavrukov, Y., Baho, M., Lopato, S., and Langridge, P. (2016). The TaDREB3 transgene transferred by conventional crossings to different genetic backgrounds of bread wheat improves drought tolerance. Plant Biotechnol. J. 14, 313–322. doi: 10.1111/pbi.12385
Shiriga, K., Sharma, R., Kumar, K., Yadav, S. K., Hossain, F., and Nepolean, T. (2014). Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize. Meta Gene 2, 407–417. doi: 10.1016/j.mgene.2014.05.001
Shiu, S.-H., Shih, M.-C., and Li, W.-H. (2005). Transcription factor families have much higher expansion rates in plants than in animals. Plant Physiol. 139, 18–26. doi: 10.1104/pp.105.065110
Shiu, S., and Bleecker, A. B. (2003). Expansion of the receptor-like kinase/pelle gene family and receptor-like proteins in Arabidopsis. Plant Physiol. 132, 530–543. doi: 10.1104/pp.103.021964
Stracke, R., Werber, M., and Weisshaar, B. (2001). The R2R3-MYB gene family in Arabidopsis thaliana. Curr. Opin. Plant Biol. 4, 447–456. doi: 10.1016/S1369-5266(00)00199-0
Thatcher, S. R., Danilevskaya, O. N., Meng, X., Beatty, M., Zastrow-Hayes, G., Harris, C., et al. (2016). Genome-wide analysis of alternative splicing during development and drought stress in maize. Plant Physiol. 170, 586–599. doi: 10.1104/pp.15.01267
Thirunavukkarasu, K., Hossain, F., Arora, K., Sharma, R., Shiriga, K., Mittal, S., et al. (2014). Functional mechanisms of drought tolerance in subtropical maize (Zea mays L.) identified using genome-wide association mapping. BMC Genomics 15:1182. doi: 10.1186/1471-2164-15-1182
Thirunavukkarasu, N., Sharma, R., Singh, N., Shiriga, K., Mohan, S., Mittal, S., et al. (2017). Genomewide expression and functional interactions of genes under drought stress in maize. Int. J. Genomics 2017:2568706. doi: 10.1155/2017/2568706
Tian, Z., Rizzon, C., Du, J., Zhu, L., Jl, B., Sa, J., et al. (2009). Do genetic recombination and gene density shape the pattern of DNA elimination in rice long terminal repeat retrotransposons ? Genome Res. 19, 2221–2230. doi: 10.1101/gr.083899.108
Tiwari, V., Chaturvedi, A. K., Mishra, A., and Jha, B. (2015). Introgression of the SbASR-1 gene cloned from a halophyte Salicornia brachiata enhances salinity and drought endurance in transgenic groundnut (Arachis hypogaea) and acts as a transcription factor. PLoS ONE 10:e0131567. doi: 10.1371/journal.pone.0131567
Van Gioi, H., Mallikarjauna, M., Mittal, S., Banduni, P., Jha, S., Dash, P., et al. (2017). Variable level of dominance of candidate genes controlling drought functional traits in maize hybrids. Front. Plant Sci. 8:940. doi: 10.3389/fpls.2017.00940
Wang, F., Chen, H. W., Li, Q. T., Wei, W., Li, W., Zhang, W. K., et al. (2015). GmWRKY27 interacts with GmMYB174 to reduce expression of GmNAC29 for stress tolerance in soybean plants. Plant J. 83, 224–236. doi: 10.1111/tpj.12879
Wang, G., Weng, L., Li, M., and Xiao, H. (2017). Response of gene expression and alternative splicing to distinct growth environments in tomato. Int. J. Mol. Sci. 18:e475. doi: 10.3390/ijms18030475
Wang, W., Qin, Q., Sun, F., Wang, Y., Xu, D., Li, Z., et al. (2016). Genome-wide differences in DNA methylation changes in two contrasting rice genotypes in response to drought conditions. Front. Plant Sci. 7:1675. doi: 10.3389/fpls.2016.01675
Wittkopp, P. J., and Kalay, G. (2011). Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence. Nat. Rev. Genet. 13, 59–69. doi: 10.1038/nrg3095
Xiang, Y., Tang, N., Du, H., Ye, H., and Xiong, L. (2008). Characterization of OsbZIP23 as a key player of the basic leucine zipper transcription factor family for conferring abscisic acid sensitivity and salinity and drought tolerance in rice. Plant Physiol. 148, 1938–1952. doi: 10.1104/pp.108.128199
Xu, Z.-Y., Kim, S. Y., Hyeon, D. Y., Kim, D. H., Dong, T., Park, Y., et al. (2013). The Arabidopsis NAC transcription factor ANAC096 cooperates with bZIP-type transcription factors in dehydration and osmotic stress responses. Plant Cell 25, 4708–4724. doi: 10.1105/tpc.113.119099
Yilmaz, A., Nishiyama, M. Y., Fuentes, B. G., Souza, G. M., Janies, D., Gray, J., et al. (2009). GRASSIUS: a platform for comparative regulatory genomics across the grasses. Plant Physiol. 149, 171–180. doi: 10.1104/pp.108.128579
Ying, S., Zhang, D.-F., Fu, J., Shi, Y.-S., Song, Y.-C., Wang, T.-Y., et al. (2012). Cloning and characterization of a maize bZIP transcription factor, ZmbZIP72, confers drought and salt tolerance in transgenic Arabidopsis. Planta 235, 253–266. doi: 10.1007/s00425-011-1496-7
Yoshida, T., Fujita, Y., Maruyama, K., Mogami, J., Todaka, D., Shinozaki, K., et al. (2015). Four Arabidopsis AREB/ABF transcription factors function predominantly in gene expression downstream of SnRK2 kinases in abscisic acid signalling in response to osmotic stress. Plant Cell Environ. 38, 35–49. doi: 10.1111/pce.12351
You, J., Zong, W., Hu, H., Li, X., Xiao, J., and Xiong, L. (2014). A STRESS-RESPONSIVE NAC1-regulated protein phosphatase gene rice protein phosphatase18 modulates drought and oxidative stress tolerance through abscisic acid-independent reactive oxygen species scavenging in rice. Plant Physiol. 166, 2100–2114. doi: 10.1104/pp.114.251116
Zhang, C., Liu, J., Zhao, T., Gomez, A., Li, C., Yu, C., et al. (2016). A drought-inducible bZIP transcription factor OsABF1 delays reproductive timing in Rice. Plant Physiol. 171, 334–343. doi: 10.1104/pp.16.01691
Zhang, Z., and Kishino, H. (2004). Genomic background predicts the fate of duplicated genes: evidence from the yeast genome. Genetics 166, 1995–1999. doi: 10.1534/genetics.166.4.1995
Keywords: drought, gene expression, gene interaction, maize, transcription factors
Citation: Mittal S, Banduni P, Mallikarjuna MG, Rao AR, Jain PA, Dash PK and Thirunavukkarasu N (2018) Structural, Functional, and Evolutionary Characterization of Major Drought Transcription Factors Families in Maize. Front. Chem. 6:177. doi: 10.3389/fchem.2018.00177
Received: 31 August 2017; Accepted: 03 May 2018;
Published: 23 May 2018.
Edited by:
Raju Datla, National Research Council Canada (NRC-CNRC), CanadaReviewed by:
Manosh Kumar Biswas, University of Leicester, United KingdomSateesh Kagale, National Research Council Canada (NRC-CNRC), Canada
Copyright © 2018 Mittal, Banduni, Mallikarjuna, Rao, Jain, Dash and Thirunavukkarasu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nepolean Thirunavukkarasu, dG5lcG9sZWFuQGdtYWlsLmNvbQ==; dG5lcG9sZWFuQHlhaG9vLmNvbQ==