Genome-wide identification of binding sites for NAC and YABBY transcription factors and co-regulated genes during soybean seedling development by ChIP-Seq and RNA-Seq
BMC Genomics volume 14, Article number: 477 (2013)
Two plant-specific transcription factors, NAC and YABBY, are involved in important plant developmental processes. However their molecular mechanisms, especially DNA binding sites and co-regulated genes, are largely unknown during soybean seedling development.
In order to identify genome-wide binding sites of specific members of the NAC and YABBY transcription factors and co-regulated genes, we performed Chromatin Immunoprecipitation Sequencing (ChIP-Seq) and RNA Sequencing (RNA-Seq) using cotyledons from soybean seedling developmental stages. Our RNA-Seq data revealed that these particular NAC and YABBY transcription factors showed a clear pattern in their expression during soybean seedling development. The highest level of their expression was found in seedling developmental stage 4 when cotyledons undergo a physiological transition from non-photosynthetic storage tissue to a metabolically active photosynthetic tissue. Our ChIP-Seq data identified 72 genes potentially regulated by the NAC and 96 genes by the YABBY transcription factors examined. Our RNA-Seq data revealed highly differentially expressed candidate genes regulated by the NAC transcription factor include lipoxygense, pectin methyl esterase inhibitor, DEAD/DEAH box helicase and homeobox associated proteins. YABBY-regulated genes include AP2 transcription factor, fatty acid desaturase and WRKY transcription factor. Additionally, we have identified DNA binding motifs for the NAC and YABBY transcription factors.
Genome-wide determination of binding sites for NAC and YABBY transcription factors and identification of candidate genes regulated by these transcription factors will advance the understanding of complex gene regulatory networks during soybean seedling development. Our data imply that there is transcriptional reprogramming during the functional transition of cotyledons from non-photosynthetic storage tissue to metabolically active photosynthetic tissue.
Understanding the transcriptional regulatory network during developmental stages has been a focus for many years at the single gene level. Recently genome-wide mapping of protein–DNA interactions enables of complete understanding of transcriptional regulation. A precise map of binding sites for transcription factors and other DNA-binding proteins is vital for deciphering the gene regulatory networks that underlie various biological processes [1–4]. Transcription factors bind to DNA sequences and regulate gene expression both in animals and plants [5–8]. Virtually all biological processes are directly regulated or influenced by transcription factors. More specifically, developmental processes are highly regulated by the appearance and disappearance of particular transcription factors [5, 7–9]. Therefore, the identification of transcription factor binding sites has immense importance for unraveling the gene regulation during developmental stages.
To study the genome-wide profiling of transcription factor binding sites, chromatin immunoprecipitation followed by sequencing (ChIP–Seq) has become one of the principal techniques [1, 2, 10]. Owing to the tremendous progress in next-generation sequencing technology, ChIP–Seq offers higher resolution, less noise and greater coverage than its array-based predecessor ChIP–chip [11, 12]. With the decreasing cost of sequencing, ChIP–Seq has become an indispensable tool for studying gene regulation. It has been widely used to study transcription factor binding, histone modifications and DNA methylation [4, 13]. The combination of ChIP-Seq and quantitative measurements of transcriptomes (RNA-Seq) is increasingly used to decipher the regulation of gene expression by transcription factors [10, 14].
Soybean (Glycine max) is one of the most economically important crops cultivated all over the world. It is an excellent source of vegetable oil and protein. However, the regulation of gene expression during soybean seedling development is still largely unexplored. Thus it is important to investigate how the genes are regulated during different developmental stages of soybean seedlings. The soybean genome was sequenced about two years ago  and this information has widened the molecular research on soybean. Soybean gene annotation revealed a number of transcription factors and many of them play vital roles during seedling development [8, 16]. Seven different stages of soybean seedling development were based on time, size of radicles, hypocotyls, roots, emergence and development of unifoliolate, and appearance of germinating cotyledons (Figure 1). It has been reported that changes in developmental stages shift the level of gene expression in cotyledon of germinating soybean seedlings  and other developmental factors . These changes are accomplished by the changes in the expression level of transcription factors [5, 8, 9]. A number of soybean developmental stage associated transcription factors have already been reported, including NAC , SNF2, zinc finger transcription factor , transcription factor homolog BTF3-like protein and transcription repressor ROM1 [17, 18]. To deeply investigate the transcription factor mediated regulation of gene expression during soybean seedling development, we constructed ChIP-Seq and RNA-Seq libraries using developmental stage specific cotyledons. There is a functional transition between developmental stage 4 and stage 5 when cotyledons undergo a physiological transition from mainly a nutrient and food reserve tissue (yellow) to an active photosynthetic (green) tissue. The physiological transition of the cotyledon is a complex process that must be under strict gene control and regulation. To investigate the expression level of transcription factors and co-regulated genes, we constructed seven different RNA-Seq libraries specific to each developmental stage. Based on our RNA-Seq data, we selected specific NAC and YABBY transcription factors which showed a clear pattern in their expression from lower level to higher level during the physiological transition of soybean seedlings.
These two transcription factors, NAC and YABBY, are involved in numerous biological processes [16, 19, 20]. The NAC family (NAM/ATAF/CUC) constitutes a highly prolific group of plant-specific TFs, with representatives in monocots, dicots, conifers and mosses [21, 22]. Many plants of commercial and scientific interest, such as Glycine max, have more than 100 different NAC proteins [16, 23]. In this study, we focused on four specific genes encoding NAC transcription factors which showed a clear expression pattern during seedling development. NAC transcription factors play important roles in plant growth, development, and stress responses [23–25]. Previously it was shown in Arabidopsis that a consensus DNA sequence (CGT[GA]) to which NAC and other relatively distant NAC TFs bind . But there is no experimental evidence in legumes for their DNA binding sites and co-regulated genes.
Another group of plant specific transcription factors is YABBY, which plays a critical role in determining organ polarity [27–30]. It is involved in the establishment of abaxial-adaxial polarity in lateral organs [19, 20, 31]. YABBY family transcription factors contain a zinc-finger domain in the amino-terminal region and a YABBY domain in the carboxyl-terminal region [19, 20, 32]. There is not much known about the molecular mechanisms of this YABBY transcription factor including the DNA binding sites and co-regulated genes during legumes development.
In order to dissect the developmental stage specific gene regulation by NAC and YABBY transcription factors, we developed polyclonal antibodies against synthetic peptides representing specific members of those transcription factors and performed ChIP-Seq experiment using pooled cotyledons from stage 4 and stage 5. We investigated the differential expression of gene models identified by ChIP-Seq between stage 3 (before the functional transition) and stage 6 (after the functional transition) using our RNA-Seq data. Our ChIP-Seq data identified many candidate genes that are regulated by specific members of NAC and YABBY transcription factors. Motif analysis discovered three separate motifs for NAC and YABBY transcription factors. RNA-Seq analysis revealed the expression in reads per kilobase of gene model per million mapped reads (RPKM) for these motif-containing Glyma models during different developmental stages. The expression analysis efficiently identified differentially expressed genes between stage 3 and stage 6. The identification of NAC and YABBY transcription factor binding sites and the potential genes regulated by these transcription factors will advance our understanding of gene regulation during legume development.
RNA-Seq reveals the differential expression of NAC and YABBY transcription factor
To understand the molecular mechanisms involved in the functional transition during soybean seedling development, we constructed seven different RNA-Seq libraries using cotyledons from each developmental stage separately. High throughput sequencing-by-synthesis (SBS) of these libraries produced 46 million to 76 million reads. Most of these reads mapped to the soybean reference genome and transcriptome of cv. Williams 82 available at the Phytozome database . Analysis of RNA-Seq data from different developmental stages revealed the differential expression of many genes including transcription factors. In this study, we focused on NAC and YABBY transcription factors which showed a clear expression pattern during soybean seedling development (Table 1 and Figure 2). Their expression gradually increases from stage 1 to stage 4 of soybean germinating cotyledons. The highest level of expression was found at stage 4. Then it gradually decreased as the germinating cotyledons develop a mature seedling. There are four specific members of NAC family and two specific members of YABBY family showed this expression pattern. The expression pattern is shown graphically for two specific members of NAC family among four and two specific members of YABBY family (Figure 2). There was a 5 to 10-fold range in differences of RPKM values between stage 1 to stage 4 and stage 5 (Table 1 and Figure 2).
ChIP-Seq libraries and detection of peaks for NAC and YABBY transcription factors
ChIP-Seq libraries were constructed using pooled cotyledons from soybean seedling developmental stage 4 and stage 5. After cross linking of cotyledon samples, chromatin complexes were isolated and sonicated for appropriate fragmentation for sequencing. We performed this experiment using different polyclonal antibodies raised against specific NAC and YABBY transcription factors using synthetic peptides (Table 2). The DNA-antibody complexes were precipitated and DNA was recovered from the complexes. Subsequently, ChIP-Seq libraries for NAC and YABBY transcription factors were constructed and high-throughput sequencing was performed. For NAC ChIP-Seq libraries, we obtained 21 million raw reads for the control library and 34 million raw reads for the antibody treated library (Table 3). Similarly sequencing of YABBY ChIP-Seq libraries generated 95 million raw reads for the control library and 86 million raw reads for the antibody treated library (Table 3).
Millions of raw reads obtained from ChIP-Seq libraries were aligned to the reference soybean genome using the ultrafast Bowtie aligner  to obtain quantitative data for genome matched reads. There are numerous peak detection algorithms available for analyzing ChIP-Seq data sets. In this experiment MACS software  was used to call peaks representing enriched binding sites. The Bowtie alignment outputs for both control and antibody treated libraries were used together as input in the MACS software. For the NAC ChIP-Seq data set, MACS detected 8246 enriched peaks with p-value < 0.05 and for the YABBY ChIP-Seq data set, it detected 18064 (Table 4, Additional file 1 and Additional file 2). The distributions of MACS detected peaks in soybean chromosomes for both NAC and YABBY transcription factors were visualized using Integrative Genomics Viewer (IGV) software  (Figure 3). Additionally, MACS software builds the peak models for NAC and YABBY transcription factors separately using the bimodal distribution of forward (+) and reverse (−) sequence tags [35, 37] (Figure 4). It calculates the estimated DNA fragment size, d, which is the distance between the peak in the forward and reverse strand. Then MACS shifts all the tags by d/2 towards the 3′ ends to get the most likely protein-DNA interaction sites [10, 11, 35, 37]. For the NAC ChIP-Seq dataset, MACS shifted all the tags 55 bp towards the 3′end to get most likely protein-DNA interactions whereas the shift was 52 bp for the YABBY transcription factor.
Locations of detected peaks and discovery of common motifs in the promoter regions
The genomic locations of MACS detected peaks were identified from the soybean gene annotation using a custom made Python programming script. We found that significant numbers of these peaks are located in the promoter region. For the YABBY ChIP-Seq dataset, 1526 peaks are located in the promoter region (Figure 5). Similarly for the NAC ChIP-Seq dataset, 974 peaks are located in the promoter region. Additionally we found that peaks are located in close proximity to the transcription start sites (TSS) (Figure 6). A motif search was performed using the most commonly used Multiple EM for Motif Elicitation (MEME) software . For MEME analysis, we included those Glyma models whose promoter region contained at least one detected peak with a fold enrichment of 3 or more. The motif analysis discovered three separate motifs for the NAC and YABBY transcription factors (Table 5). For the NAC transcription factor, three commonly found motifs were G[AT]G[AG]G[AG]GA, C[AC]C[GA][TC][GA]CC and TGGGCC (Figure 7). The first one matched to a known zinc finger motif and the last two were identified as leucine zippers in the database of plant transcription factor binding motifs, JASPAR CORE plants [39, 40]. Similarly the three most commonly found motifs for YABBY transcription factors are CC [CA][TC]C[TA][CT]C, GA[AG]AGAAA and CCCCAC (Figure 7). The first two motifs matched to a known zinc finger motif and the last one was an AP2 MBD-like motif.
ChIP-Seq coupled with RNA-Seq reveals candidate genes regulated by NAC and YABBY transcription factors
Our ChIP-Seq data identified 72 potential genes which are regulated by specific members of NAC transcription factor (Table 6 and Additional file 3). In similar way, we identified 96 potential candidate genes which are regulated by YABBY transcription factors (Table 7 and Additional file 4). We investigated the differential expression of NAC and YABBY regulated genes between stage 3 (before the functional transition) and stage 6 (after the functional transition) using our RNA-Seq data. RNA-Seq analysis revealed the expression in RPKM for these NAC and YABBY regulated candidate genes during different developmental stages. For differential expression analysis the DESeq package  efficiently identified differentially expressed genes between stage 3 and stage 6 with statistically significant P-values < 0.05. For the NAC transcription factor, we found that 10 candidate genes are up-regulated, 21 candidate genes are down-regulated and 41 candidate genes show no significant difference in their expression. Similarly for the YABBY transcription factor, we found that 19 candidate genes are up-regulated, 27 candidate genes are down-regulated and 50 candidate genes show no significant difference in their expression. The majority of these differentially expressed candidate genes are found to be involved in plant developmental processes.
The regulation of gene expression by transcription factors is a quite complex and coordinated process. Recently the ENCODE (ENCyclopedia of DNA Elements) project has generated chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) data sets for a large number of transcription factors using different human cell lines to identify genome-wide functional and regulatory DNA elements [42, 43]. Although better understood in model plants such as Arabidopsis and rice [5, 44, 45] knowledge is scarce in most other plants. A number of transcription factor binding sites have been identified in Arabidopsis using recently developed techniques such as Chromatin Immunoprecipitation Sequencing (ChIP-Seq) [46–49]. Soybean is a polyploid crop having a complex and large genome . To date, there are no reports of identification of soybean transcription factor binding sites using the high throughput ChIP-Seq technique. In order to identify two developmentally important transcription factors binding sites during soybean seedling development, we used a combination of experimental and bioinformatics approaches. In this study, ChIP-Seq and RNA-Seq were used to dissect the gene regulatory networks for NAC and YABBY transcription factors during soybean seedling development. We constructed seven RNA-Seq libraries using cotyledons from seven different seedling developmental stages separately to see the expression level of transcription factors and their co-regulated genes. Later we constructed separate ChIP-Seq libraries for specific NAC and YABBY transcription factors using pooled cotyledons from soybean seedling developmental stage 4 and stage 5 when the cotyledons undergo a functional transition from non-photosynthetic storage tissues to metabolically active photosynthetic tissues.
The NAC transcription factor is a plant specific transcription factor family which plays important roles in plant growth, development and stress responses [23–25]. Glycine max has more than 100 different NAC proteins [16, 23]. Although NAC transcription factor family is quite large, our RNA-Seq data showed that there are only four specific members of NAC family expressed and showed a clear expression pattern during soybean seedlings development (Table 1 and Figure 2). Additionally, we performed the multiple sequence alignment of these four members of NAC family and found a high homology among their sequences (Additional file 5). These four members of NAC family possess that short peptide sequence used for developing the antibody and they are closely related. For the ChIP-Seq experiment, we used germinating cotyledons from stage 4 and stage 5 which are the transition stages. Thus, our antibody is specific for these four members of the NAC family since they show high homology in their sequences and are the only members expressed during the physiological transition at stage 4 and stage 5. The analysis of ChIP-Seq libraries (control and antibody treated) for the NAC transcription factor using MACS software detected 8246 highly enriched peaks with statistical significance P < 0.05. A significant number of these peaks are associated with soybean gene models. We found that 974 peaks are located in the promoter region of soybean gene models. For MEME analysis, we selected those Glyma models whose promoter region contains at least one detected peak with a fold enrichment of 3 or more over the control. We found three common DNA binding motifs, two of them matched to leucine zipper and one matched to a zinc finger (Figure 7 and Table 5). Previously it had been reported in Arabidopsis that the NAC transcription factor binding site contains the consensus DNA sequence (CGT[GA]) . One of our identified common motifs was C[AC]C[GA][TC][GA]CC which contains the previously identified motif in Arabidopsis, thus corroborating our discovery of DNA binding motifs for the NAC transcription factor in soybean.
To dissect the gene regulatory network of a particular transcription factor, it is important to study the expression of co-regulated genes. In this study, we have identified 72 genes potentially regulated by a NAC transcription factor based on our ChIP-Seq and RNA-Seq data (Table 6 and Additional file 3). Using our developmental stage specific RNA-Seq data, we investigated their expression levels. Our particular interest was on developmental stage 3 which is before the functional transition and developmental stage 6 which is after the functional transition. DESeq analysis showed differential expression of a number of candidate genes at p-value < 0.05. We focused on 10 up-regulated and 21 down-regulated genes to see the level of expression difference in between stage 3 and stage 6. From the RNA expression data, the highest level of expression difference was found with genes annotated as lipoxygense, pectin methylesterase inhibitor (PMEI), DEAD/DEAH box helicase and Homeobox associated proteins. DESeq analysis also showed very low p-values corresponding to these gene models indicating they are significantly differentially expressed.
Among those highly differentially expressed genes, lipoxygenase has been proposed to be involved in reserve lipid mobilization during soybean seed germination [17, 50, 51]. Our RNA-Seq data showed that the lipoxygenase gene is up-regulated. Once germination is triggered, lipids need to be mobilized by the action of lipoxygenase and ultimately triacylglycerols are degraded to act as a carbon and energy source for the developing seedlings. Another candidate gene encodes a cupin domain containing protein, which has been reported to be involved in seed germination and early seedling development [52, 53]. We found that the cupin gene is down-regulated by the specific NAC transcription factor. This might be due to the fact that after the functional transition the seedling tends to shift towards a normal photosynthetic cycle instead of the glyoxylate cycle. In this study, we found pectin methylesterase inhibitors (PMEIs) are regulated by the specific members of NAC transcription factor. Pectin, one of the main components of the plant cell wall, is continually modified and remodeled during plant growth and development [54, 55]. Pectin methylesterases (PMEs) catalyse the demethylesterification of cell wall pectins [55, 56]. In many developmental processes, PMEs are regulated by either differential expression or posttranslational control by pectin methylesterases inhibitors (PMEIs) . These PMEI inhibitors play significant roles in plant growth, cell division, and expansion [56, 57]. We identified PMEI as a NAC regulated potential candidate gene and the expression of PMEI gene is up-regulated by this NAC transcription factor, indicating that PMEI reduces the activity of PMEs during later stages of seedling development.
Another important gene regulated by the NAC transcription factor is the DEAD/DEAH box helicases which are ubiquitous enzymes that catalyze the unwinding of energetically stable duplex DNA (DNA helicases) or duplex RNA secondary structures (RNA helicases) [58–60]. Most helicases are members of the DEAD-box protein superfamily and play essential roles in basic cellular processes such as DNA replication, repair, recombination, transcription, ribosome biogenesis and translation initiation [59, 61]. Therefore, helicases might be playing an important role in regulating plant growth and development. Our ChIP-Seq results showed the potential NAC transcription factor binding sites in the promoter region of DEAD/DEAH box helicase gene (Figure 8a). ChIP-Seq and RNA-Seq analysis together showed that the DEAD/DEAH box helicase gene is down-regulated by the NAC transcription factor. Using our RNA-Seq data, we found that a particular Homeotic (HOX) gene was down- regulated by the NAC transcription factor. The plausible explanation for this down regulation might be that HOX genes are known to be involved in flower development [62, 63] and thus it is expected to be down-regulated during seedling development. Taken together, these results indicate that NAC transcription factors act in multiple pathways to regulate gene expression that facilitate the functional transition of the cotyledons during legume seedling growth.
The YABBY transcription factor family is a group of plant specific transcription factors which play important roles in organ polarity determination [28–30]. Members of the YABBY gene family are characterized by two conserved domains, a C2C2 zinc finger-like domain towards the amino terminus and a helix-loop-helix, which is called the YABBY domain, with sequence similarity to the first two helices of the HMG box towards the carboxyl end of the protein [19, 20, 32]. Glycine max has multiple isoforms of the YABBY transcription factor; however there is not much known about the molecular mechanisms for this transcription factor specifically its DNA binding sites and co-regulated genes. Our RNA-Seq data revealed that there are ten members of the YABBY transcription factor family expressed during the physiological transition. We performed a multiple sequence alignment among the ten members of the YABBY transcription factors. Two of them showed high sequence homology and they are closely related (Additional file 6). Only these two members possess the short peptide sequence used for developing the antibody. Thus, our antibody is specific for those two members of the YABBY family that showed the clear expression pattern during the functional transition. The analysis of ChIP-Seq libraries (control and antibody treated) for the YABBY transcription factor detected 18064 highly enriched peaks with statistical p-value < 0.05. A significant number of these peaks are associated with gene models. We found that 1526 peaks are located in the promoter region of soybean gene models. We found three common DNA binding motifs by MEME analysis; two of them match to the zinc finger motif and one matches to the AP2 MBD-like motif (Figure 7 and Table 5).
In this study, we have identified YABBY regulated genes based on our ChIP-Seq and RNA-Seq data. We found that there are 96 genes potentially regulated by the YABBY transcription factor (Table 7 and Additional file 4). Using our developmental stage specific RNA-Seq data, we have investigated their expression levels before the functional transition (developmental stage 3) and after the functional transition (developmental stage 6). DESeq analysis showed differential expression of a number of candidate genes at p-value < 0.05. We carefully looked at the expression data for the 19 up-regulated and 27 down-regulated genes to see the level of expression difference between stage 3 and stage 6. From the RNA expression data, the highest level of expression difference was found with genes annotated as protein AP2 (APETALA2) transcription factor, fatty acid desaturase and WRKY DNA binding domain protein as well as some other genes with no functional annotation. Among them AP2 is a very well-known transcription factor unique to plants, whose distinguishing characteristic is its AP2 DNA-binding domain [64–66]. It plays a key role in several developmental processes like floral organ identity determination and control of leaf epidermal cell identity and is under strict regulation during developmental processes [65, 67–70]. Our ChIP-Seq results showed the potential YABBY transcription factor binding sites in the promoter region of AP2 gene (Figure 8b). Since it is mainly involved in flower developmental, it is down-regulated at the later stages of seedling development.
The existence of repeated DNA-binding domains not unique to AP2 transcription factors, the WRKY family of transcription factors also may contain a repeated DNA-binding domain . Our RNA-Seq data revealed that there are 53 genes encoding WRKY transcription factors expressed at level (≥ 10 RPKM) in at least one of seven stages of soybean seedling development (Additional file 7). Based on our ChIP-Seq result, a specific WRKY Transcription factor is one among five major candidates regulated by YABBY transcription factor. Previous reports show that the WRKY transcription factor has been involved in the numerous plant developmental processes [72–76]. Specific members of WRKY transcription factor family are playing important role in seed development [72, 73]. However it is quite difficult to pinpoint the regulation by the WRKY transcription factor since there are large numbers of WRKY transcription factors in soybean. Another YABBY regulated candidate gene is Fatty Acid Desaturase (FAD). It encodes the main enzyme responsible for polyunsaturated lipid synthesis in developing seeds of oil crops [77, 78]. Our results showed that fatty acid desaturase was down-regulated by the YABBY transcription factor in agreement with lipid catabolism being more preferable during seedling development than lipid biosynthesis which needs fatty acid desaturase. Our ChIP-Seq results reveal the genome-wide view of binding sites for the YABBY transcription factor and RNA-Seq demonstrates the resultant changes in expression of regulated genes that influence the physiological transition of the soybean cotyledon from a storage tissue to a metabolically active tissue during seedling growth.
ChIP-Seq demonstrates promising potential as a new tool in understanding genome-wide binding sites for transcription factors and transcriptional gene regulatory networks. Our genome-wide identification of NAC and YABBY transcription factor binding sites using antibodies to synthetic peptides representing these rare abundance transcription factors will help to better understand the transcriptional gene regulatory network during the functional transition of cotyledons from a storage tissue to a metabolically active photosynthetic tissue. The discovery of common DNA binding motifs and identification of regulated genes opens a new avenue to pinpoint the molecular mechanisms of these two important transcription factors during seedling growth. Combining ChIP-Seq and RNA-Seq results advances understanding of the underlying genetic mechanisms involved in the functional transition as well as their regulation and control systems throughout the soybean seedling developmental process.
Plant materials and growth conditions
Four soybean (Glycine max cv. Williams) seeds were planted per small pot (4.5 inches) containing Universal SB300 soil mix. A total of 25 pots (100 seeds) were initially used to collect and pool 6 individual cotyledon samples per developmental stage. Plants were grown for approximately 7–8 days with regular watering. A biological replicate was performed with another 25 pots to collect tissues in a similar way. Seven different stages during the development of soybean seedlings were defined based on time, size of radicles, hypocotyls, roots and appearance of germinating cotyledons. Stage 1: Imbibed seed for 24 hours; pre-emerging hypocotyls. Stage 2: Yellow cotyledons; emerging radicle 8–10 mm long. Stage 3: Yellow cotyledons with slightly green edges; hypocotyls15-20 mm long. Stage 4: Yellow-green cotyledons; hypocotyls 30–35 mm long. Stage 5: Yellow-green cotyledons above the ground; primary roots starting to develop. Stage 6: Mostly green cotyledons above the ground; growing straight from the hypocotyl. Stage 7: Fully green cotyledons; plants 6–7 cm long above the ground; the root system fully developed; cotyledons upright; unifoliolate exposed. For the RNA-Seq experiment, cotyledons from each of these developmental stages were collected and then frozen in liquid nitrogen. Subsequently the tissue was freeze dried and stored at −80°C. For ChIP-Seq experiment, fresh cotyledons from stage 4 and stage 5 were used.
RNA-Seq library construction and data analysis
Total RNA was extracted separately for seven different developmental stages from freeze dried cotyledons using a modified McCarty method  using phenol-chloroform extraction and lithium chloride precipitation. A biological replicate was performed to extract RNA in a similar way. Library construction and high-throughput sequencing were carried out by the Illumina HiSeq2000 at the Keck Center, University of Illinois at Urbana-Champaign. The 100 bp RNA-Seq reads were mapped to the 78,773 high and low confidence soybean gene models  using the ultrafast Bowtie aligner  with up to 3 mismatches. RNA-Seq data was normalized in reads per kilobase of gene model per million mapped reads (RPKM) . The DESeq package  was used to determine differential expression between developmental stage 3 and stage 6 and calculate p-values. If the p-value was < 0.05, we considered that gene as significantly differentially expressed gene between two developmental stages. The expression could be up-regulated or down-regulated based on the corresponding RPKM values.
ChIP-Seq library construction
Cotyledons from soybean seedling developmental stage 4 and stage 5 were collected for the ChIP-Seq experiment performed using previously described methods [80, 81]. Briefly, 0.08 g of fresh weight of soybean cotyledons from stage 4 or stage 5 were cross sectioned with a razor blade and then cross linked with 1% formaldehyde under vacuum. Immediately the samples were ground to powder in liquid nitrogen. The chromatin complexes were isolated following previously established protocols [80, 81]. Later, the chromatin was sonicated to shear DNA into 200–600 bp fragments using the 15% power setting and fifteen times for 20 second pulses using a Branson digital probe sonifier. Sample containing tubes were kept on ice while the sonication was performed. Subsequently, the sonicated DNA was incubated with a polyclonal antibody developed against the YABBY or NAC transcription factors. All the antibodies were produced by GenScript Corporation. They used the Jameson and Wolf (JW) prediction algorithm  to design synthetic peptides for the production of antibody against YABBY and NAC transcription factor [Table 2]. Separate controls which were not treated with antibody, but used preimmune sera, were used for each experiment. Then DNA-antibody complexes were precipitated following standard protocol [81, 83] and DNA was recovered by dissociating the complexes. ChIP-Seq library construction and high-throughput sequencing was carried out by the Illumina HiSeq2000 at the Keck Center, University of Illinois at Urbana-Champaign.
ChIP-Seq data analysis
Sequencing of ChIP-Seq libraries produced millions of raw reads which were aligned to the reference soybean genome using the ultrafast Bowtie aligner  to get the number of genome matched reads. The length of our sequence reads was 100 bp and we allowed 3 mismatches for Bowtie alignment. The experiment was conducted in two conditions, the control library and the antibody treated library. MACS software  with specific parameters (bandwidth 300 bp; mfold, 30; P-value of 1.00e-05) was used to call peaks representing enriched binding sites. The Bowtie alignment output for both control and antibody treated libraries was used together as input to the MACS software to detect a number of peaks for the potential binding sites for the YABBY or NAC transcription factors separately. Since ChIP-DNA fragments are equally likely to be sequenced from both ends, the tag density around a true binding site should show a bimodal enrichment pattern, with forward strand tags enriched upstream of binding sites and reverse strand tags enriched downstream of binding sites [10, 11, 35, 37]. MACS software takes advantage of this bimodal pattern to empirically model the shifting size to better locate the precise binding sites. It randomly samples 1,000 of these high-quality peaks, separates their forward (+) and reverse (−) tags, and aligns them by the midpoint between their forward and reverse tag centers [35, 37]. MACS calculated estimated DNA fragment size, d which is the distance between the peak in the forward and reverse strand. Then MACS shifts all the tags by d/2 toward the 3′ ends to get the most likely protein-DNA interaction sites [10, 11, 35, 37]. Then the genomic locations of these peaks were identified from the soybean gene annotation file from the Phytozome database  using a custom made Python programming script. Using that programming script, all binding peaks were sorted based on the following criteria: (1) if a binding site resides in the gene body, it will be further categorized according to its location in the gene body (i.e., 5′-untranslated region, exon, intron, or 3′-untranslated region); (2) if a binding site is localized in the 1000-bp region upstream of the transcription start site of a gene, it is classified as a binding site in the promoter region in our study; the binding sites not selected by the above criteria were defined as the binding sites in the intergenic regions. The outputs of the analysis, specifically the detected peaks were visualized in the Integrative Genomics Viewer genome browser .
A motif search was performed using the most widely used MEME software . For MEME analysis, gene models were selected based on the location of detected peaks and fold enrichment. In this analysis, we included those gene models whose promoter region contains at least one detected peak and a fold enrichment of 3 or more. For promoter associated peaks, 250 bp sequences from both sides of peak summits were retrieved. These 500 bp sequences for associated gene models were given as input in MEME software to identify common motifs. Some of these identified motifs were matched to known motifs in the plant transcription factor binding sites database, JASPAR CORE plants [39, 40].
The high-throughput sequencing data for ChIP-Seq libraries are available under NCBI-GEO  series accession no. GSE42422. In addition, RNA-Seq data for seven developmental stages with two biological replicates are available under NCBI-GEO  series accession no. GSE42550.
Johnson DS, Mortazavi A, Myers RM, Wold B: Genome-wide mapping of in vivo protein-DNA interactions. Sci. 2007, 316: 1497-1502. 10.1126/science.1141319.
Jothi R, Cuddapah S, Barski A, Cui K, Zhao K: Genome-wide identification of in vivo protein–DNA binding sites from ChIP-seq data. Nucleic Acids Res. 2008, 36: 5221-5231. 10.1093/nar/gkn488.
Valouev A, Johnson DS, Sundquist A, Medina C, Anton E, Batzoglou S, Myers RM, Sidow A: Genome-wide analysis of transcription factor binding sites based on ChIP-Seq data. Nat Methods. 2008, 5: 829-834. 10.1038/nmeth.1246.
Barski A, Zhao K: Genomic location analysis by ChIP-seq. J Cell Biochem. 2009, 107: 11-18. 10.1002/jcb.22077.
Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, et al: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Sci. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.
Thibaud-Nissen F, Wu H, Richmond T, Redman JC, Johnson C, Green R, Arias J, Town CD: Development of Arabidopsis whole-genome microarrays and their application to the discovery of binding sites for the TGA2 transcription factor in salicylic acid-treated plants. Plant J. 2006, 47: 152-162. 10.1111/j.1365-313X.2006.02770.x.
MacArthur S, Li XY, Li J, Brown JB, Chu HC, et al: Developmental roles of 21 Drosophila transcription factors are determined by quantitative differences in binding to an overlapping set of thousands of genomic regions. Genome Biol. 2009, 10: R80-10.1186/gb-2009-10-7-r80.
Udvardi MK, Kakar K, Wandrey M, Montanari O, Murray J, Andriankaja A, Zhang J-Y, Benedito V, Hofer JMI, Chueng F, Town CD: Legume transcription factors: global regulators of plant development and response to the environment. Plant Physiol. 2007, 144: 538-549. 10.1104/pp.107.098061.
Singh KB: Transcriptional regulation in plants: the importance of combinatorial control. Plant Physiol. 1998, 118: 1111-1120. 10.1104/pp.118.4.1111.
Pepke S, Wold B, Mortazavi A: Computation for ChIP-seq and RNA-seq studies. Nat Methods. 2009, 6: S22-S32. 10.1038/nmeth.1371.
Park PJ: ChIP-seq: advantages and challenges of a maturing technology. Nat Rev Genet. 2009, 10: 669-680.
Laajala TD, Raghav S, Tuomela S, Lahesmaa R, Aittokallio T, Elo LL: A practical comparison of methods for detecting transcription factor binding sites in ChIP-seq experiments. BMC Genomics. 2009, 10: 618-10.1186/1471-2164-10-618.
Aleksic J, Russell S: ChIPing away at the genome: the new frontier travel guide. Mol Biosyst. 2009, 5: 1421-1428. 10.1039/b906179g.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.
Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, et al: Genome sequence of the palaeopolyploid soybean. Nat. 2010, 463: 178-183. 10.1038/nature08670.
Libault M, Farmer A, Joshi T, Takahashi K, Langley RJ, Franklin LD, He J, Xu D, May G, Stacey G: An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. Plant J. 2010, 63: 86-99.
Gonzalez DO, Vodkin LO: Specific elements of the glyoxylate pathway play a significant role in the functional transition of the soybean cotyledon during seedling development. BMC Genomics. 2007, 8: 468-10.1186/1471-2164-8-468.
Jones SI, Gonzalez DO, Vodkin LO: Flux of transcript patterns during soybean seed development. BMC Genomics. 2010, 11: 136-10.1186/1471-2164-11-136.
Siegfried KR, Eshed Y, Baum SF, Otsuga D, Drews GN, Bowman JL: Members of the YABBY gene family specify abaxial cell fate in Arabidopsis. Dev. 1999, 126: 4117-4128.
Bowman JL: The YABBY gene family and abaxial cell fate. Curr Opin Plant Biol. 2000, 3: 17-22. 10.1016/S1369-5266(99)00035-7.
Ooka H, Satoh K, Doi K, Nagata T, Otomo Y, Murakami K, Matsubara K, Osato N, Kawai J, Carninci P, et al: Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana. DNA Res. 2003, 10: 239-247. 10.1093/dnares/10.6.239.
Hu R, Qi G, Kong Y, Kong D, Gao Q, Zhou G: Comprehensive analysis of NAC domain transcription factor gene family in Populus trichocarpa. BMC Plant Biol. 2010, 10: 145-167. 10.1186/1471-2229-10-145.
Tran LS, Quach TN, Guttikonda SK, Aldrich DL, Kumar R, Neelakandan A, Valliyodan B, Nguyen HT: Molecular characterization of stress-inducible GmNAC genes in soybean. Mol Genet Genomics. 2009, 281: 647-664. 10.1007/s00438-009-0436-8.
Kim YS, Kim SG, Park JE, Park HY, Lim MH, Chua NH, Park CM: A membrane-bound NAC transcription factor regulates cell division in Arabidopsis. Plant Cell. 2006, 18: 3132-3144. 10.1105/tpc.106.043018.
Zhong R, Demura T, Ye ZH: SND1, a NAC domain transcription factor, is a key regulator of secondary wall synthesis in fibers of Arabidopsis. Plant Cell. 2006, 18: 3158-3170. 10.1105/tpc.106.047399.
Olsen AN, Ernst HA, Leggio LL, Skriver K: DNA-binding specificity and molecular functions of NAC transcription factors. Plant Sci. 2005, 169: 785-797. 10.1016/j.plantsci.2005.05.035.
Sieber P, Petrascheck M, Barberis A, Schneitz K: Organ polarity in Arabidopsis. NOZZLE physically interacts with members of the YABBY family. Plant Physiol. 2004, 135: 2171-2185.
Kumaran MK, Bowman JL, Sundaresan V: YABBY polarity genes mediate the repression of KNOX homeobox genes in Arabidopsis. Plant Cell. 2002, 14: 2761-2770. 10.1105/tpc.004911.
Eshed Y, Baum SF, Perea JV, Bowman JL: Establishment of polarity in lateral organs of plants. Curr Biol. 2001, 11: 1251-1260. 10.1016/S0960-9822(01)00392-X.
Bowman JL, Eshed Y, Baum SF: Establishment of polarity in angiosperm lateral organs. Trends Genet. 2002, 18: 134-141. 10.1016/S0168-9525(01)02601-4.
Goldshmidt A, Alvarez JP, Bowman JL, Eshed Y: Signals derived from YABBY gene activities in organ primordia regulate growth and partitioning of Arabidopsis shoot apical meristems. Plant Cell. 2008, 20: 1217-1230. 10.1105/tpc.107.057877.
Bowman JL, Smyth DR: CRABS CLAW, a gene that regulates carpel and nectary development in Arabidopsis, encodes a novel protein with zinc finger and helix-loop-helix domains. Dev. 1999, 126: 2387-2396.
Joint Genome Institute/Phytozome/. http://www.phytozome.net/soybean.php,
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.
Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP: Integrative genomics viewer. Nat Biotechnol. 2011, 29: 24-26. 10.1038/nbt.1754.
Feng J, Liu T, Qin B, Zhang Y, Liu XS: Identifying ChIP-seq enrichment using MACS. Nat Protocols. 2012, 7 (9): 1728-1740. 10.1038/nprot.2012.101.
Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34: W369-W373. 10.1093/nar/gkl198.
Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32: D91-D94. 10.1093/nar/gkh012.
Lenhard B, Wasserman WW: TFBS: computational framework for transcription factor binding site analysis. Bioinformatics. 2002, 18: 1135-1136. 10.1093/bioinformatics/18.8.1135.
Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.
Gerstein MB, Kundaje A, Hariharan M, Landt SG, Yan KK, Cheng C, Mu XJ, Khurana E, Rozowsky J, et al: Architecture of the human regulatory network derived from ENCODE data. Nat. 2012, 489 (7414): 91-100. 10.1038/nature11245.
Bernstein BE, Birney E, Dunham I, Green ED, Gunter C, Snyder M: An integrated encyclopedia of DNA elements in the human genome. Nat. 2012, 489 (7414): 57-74. 10.1038/nature11247.
Vandepoele K, Quimbaya M, Casneuf T, De Veylder L, Peer Y: Unraveling Transcriptional Control in Arabidopsis Using cis-Regulatory Elements and Coexpression Networks. Plant Physiol. 2009, 150 (2): 535-546. 10.1104/pp.109.136028.
Xue LJ, Zhang JJ, Xue HW: Genome-wide analysis of the complex transcriptional networks of rice developing seeds. PLoS One. 2012, 7: e31081-10.1371/journal.pone.0031081.
Brandt R, Salla-Martret M, Bou-Torrent J, Musielak T, Stahl M, Lanz C, Ott F, et al: Genome-wide binding-site analysis of REVOLUTA reveals a link between leaf patterning and light-mediated growth responses. Plant J. 2012, 72: 31-42. 10.1111/j.1365-313X.2012.05049.x.
Immink R, Pose D, Ferrario S, Ott F, Kaufmann K, Leal Valentim F, De Folter S, Van der Wal F, van Dijk ADJ, Schmid M, et al: Characterisation of SOC1’s central role in flowering by the identification of its upstream and downstream regulators. Plant Physiol. 2012, 160: 433-449. 10.1104/pp.112.202614.
Oh E, Kang H, Yamaguchi S, Park J, Lee D, Kamiya Y, Choi G: Genome-wide analysis of genes targeted by PHYTOCHROME INTERACTING FACTOR 3-LIKE5 during seed germination in Arabidopsis. Plant Cell. 2009, 21: 403-419. 10.1105/tpc.108.064691.
Ouyang X, Li J, Li G, Li B, Chen B, Shen H, Huang X, Mo X, Wan X, Lin R, et al: Genome-wide binding site analysis of FAR-RED ELONGATED HYPOCOTYL3 reveals its novel function in Arabidopsis development. Plant Cell. 2011, 23: 2514-2535. 10.1105/tpc.111.085126.
Rylott EL, Hooks MA, Graham IA: Co-ordinate regulation of genes involved in storage lipid mobilization in Arabidopsis thaliana. Biochem Soc Trans. 2001, 29: 283-287. 10.1042/BST0290283.
Wang C, Croft KPC, Jarlfors U, Hildebrand DF: Subcellular localization studies indicate that lipoxygenases 1 to 6 are not involved in lipid mobilization during soybean germination. Plant Physiol. 1999, 120: 227-235. 10.1104/pp.120.1.227.
Lapik YR, Kaufman LS: The Arabidopsis cupin domain protein atpirin1 interacts with the G protein alpha-subunit GPA1 and regulates seed germination and early seedling development. Plant Cell. 2003, 15: 1578-1590. 10.1105/tpc.011890.
Dunwell JM, Purvis A, Khuri S: Cupins: The most functionally diverse protein superfamily?. Phytochemistry. 2004, 65 (1): 7-17. 10.1016/j.phytochem.2003.08.016.
Ridley BL, O’Neill MA, Mohnen D: Pectins: structure, biosynthesis, and oligogalacturonide-related signaling. Phytochemistry. 2001, 57: 929-967. 10.1016/S0031-9422(01)00113-3.
Micheli F: Pectin methylesterases: cell wall enzymes with important roles in plant physiology. Trends Plant Sci. 2001, 6: 414-419. 10.1016/S1360-1385(01)02045-3.
Jolie RP, Duvetter T, Van Loey AM, Hendrickx ME: Pectin methylesterase and its proteinaceous inhibitor: a review. Carbohydr Res. 2010, 345: 2583-2595. 10.1016/j.carres.2010.10.002.
Hong MJ, Kim DJ, Lee TG, Jeon WB, Seo YW: Functional characterization of pectin methylesterase inhibitor (PMEI) in wheat. Genes Genet Syst. 2010, 85: 97-106. 10.1266/ggs.85.97.
Linder P, Lasko PF, Ashburner M, Leroy P, Nielsen PJ, Nishi K, et al: Birth of the D-E-A-D box. Nat. 1989, 337: 121-122. 10.1038/337121a0.
Tuteja N, Tuteja R: Unraveling DNA helicases: motif, structure, mechanism and functions. Physica A. 2006, 372: 70-83. 10.1016/j.physa.2006.05.014.
Aubourg S, Kreis M, Lecharny A: The DEAD box RNA helicase family in Arabidopsis thaliana. Nucleic Acids Res. 1999, 27: 628-636. 10.1093/nar/27.2.628.
Umate P, Tuteja R, Tuteja N: Architectures of the unique domains associated with the DEAD-box helicase motif. Cell Cycle. 2010, 9: 4228-4235. 10.4161/cc.9.20.13635.
Ma H, McMullen MD, Finer JJ: Identification of a homeobox-containing gene with enhanced expression during soybean (Glycine max L.) somatic embryo development. Plant Mol Biol. 1994, 24: 465-473. 10.1007/BF00024114.
Chan RL, Gago GM, Palena CM, Gonzalez DH: Homeoboxes in plant development. Biochim Biophys Acta. 1998, 1442: 1-19. 10.1016/S0167-4781(98)00119-5.
Dietz KJ, Vogel MO, Viehhauser A: AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signaling. Protoplasma. 2010, 245: 3-14. 10.1007/s00709-010-0142-8.
Jofuku KD, den Boer BG, Van Montagu M, Okamuro JK: Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 1994, 6: 1211-1225.
Wang C, Wang H, Zhang J, Chen S: A seed-specific AP2- domain transcription factor from soybean plays a certain role in regulation of seed germination. Sci China Ser C. 2008, 51: 336-345. 10.1007/s11427-008-0044-6.
Maes T, Van de Steene N, Zethof J, Karimi M, D’Hauw M, Mares G, Van Montagu M, Gerats T: Petunia AP2-like genes and their role in flower and seed development. Plant Cell. 2001, 13: 229-244.
Shamimuzzaman M, Vodkin L: Identification of soybean seed developmental stage-specific and tissue-specific miRNA targets by degradome sequencing. BMC Genomics. 2012, 13: 310-10.1186/1471-2164-13-310.
Zabala G, Campos E, Varala KK, Bloomfield S, Jones SI, Win H, Tuteja JH, Calla B, Clough SJ, Hudson M, Vodkin LO: Divergent patterns of endogenous small RNA populations from seed and vegetative tissues of Glycine max. BMC Plant Biol. 2012, 12: 177-10.1186/1471-2229-12-177.
Tsuwamoto R, Yokoi S, Takahata Y: Arabidopsis EMBRYOMAKER encoding an AP2 domain transcription factor plays a key role in developmental change from vegetative to embryonic phase. Plant Mol Biol. 2010, 73: 481-492. 10.1007/s11103-010-9634-3.
Rushton PJ, Macdonald H, Huttly AK, Lazarus CM, Hooley R: Members of a new family of DNA-binding proteins bind to a conserved cis-element in the promoters of α-Amy2 genes. Plant Mol Biol. 1995, 29: 691-702. 10.1007/BF00041160.
Zhang ZL, Xie Z, Zou X, Casaretto J, Ho TH, Shen QJ: A rice WRKY gene encodes a transcriptional repressor of the gibberellin signaling pathway in aleurone cells. Plant Physiol. 2004, 134 (4): 1500-1513. 10.1104/pp.103.034967.
Xie Z, Zhang ZL, Zou X, Yang G, Komatsu S, Shen QJ: Interactions of two abscisic-acid induced WRKY genes in repressing gibberellin signaling in aleurone cells. Plant J. 2006, 46 (2): 231-242. 10.1111/j.1365-313X.2006.02694.x.
Eulgem T: The WRKY superfamily of plant transcription factors. Trends Plant Sci. 2000, 5: 199-206. 10.1016/S1360-1385(00)01600-9.
Ciolkowski I: Studies on DNA-binding selectivity of WRKY transcription factors lend structural clues into WRKY domain function. Plant Mol Biol. 2008, 68: 81-92. 10.1007/s11103-008-9353-1.
Jiang W, Yu D: Arabidopsis WRKY2 transcription factor mediates seed germination and postgermination arrest of development by abscisic acid. BMC Plant Biol. 2009, 9: 96-10.1186/1471-2229-9-96.
Gu K, Yi C, Tian D, Sangha JS, Hong Y, Yin Z: Expression of fatty acid and lipid biosynthetic genes in developing endosperm of Jatropha curcas. Biotechnol Biofuels. 2012, 5: 47-10.1186/1754-6834-5-47.
Zhang J, Liu H, Sun J, Li B, Zhu Q, et al: Arabidopsis Fatty Acid Desaturase FAD2 is required for salt tolerance during seed germination and early seedling growth. PLoS One. 2012, 7 (1): e30355-10.1371/journal.pone.0030355.
McCarty DR: A simple method for extractions of RNA from maize tissue. Maize Genetics Cooperation News Letter. 1986, 60: 61-
Kaufmann K, Muino JM, Osteras M, Farinelli L, Krajewski P, Angenent GC: Chromatin immunoprecipitation (ChIP) of plant transcription factors followed by sequencing (ChIP-SEQ) or hybridization to whole genome arrays (ChIP-CHIP). Nat Protocols. 2010, 5 (3): 457-472. 10.1038/nprot.2009.244.
Ricardi MM, González RM, Iusem ND: Protocol: fine-tuning of a Chromatin Immunoprecipitation (ChIP) protocol in tomato. Plant Methods. 2010, 6: 11-10.1186/1746-4811-6-11.
Jameson BA, Wolf H: The antigenic index: a novel algorithm for predicting antigenic determinants. Comput Appl Biosci. 1988, 4 (1): 181-186.
Haring M, Offermann S, Danker T, Horst I, Peterhansel C, Stam M: Chromatin immunoprecipitation: optimization, quantitative analysis and data normalization. Plant Methods. 2007, 3: 11-10.1186/1746-4811-3-11.
Gene Expression Omnibus, a genomics repository database. http://www.ncbi.nlm.nih.gov/geo/,
We are grateful to Achira Kulasekara and Cameron Lowe for help with data analysis. We would like to extend our thanks to Sarah Jones for reading manuscript and making valuable comments. The research was funded by support from the United Soybean Board, the Illinois Soybean Association, and the USDA.
The authors declare that they have no competing interests.
MS designed experiments, performed RNA extractions, conducted ChIP-Seq experiment, analyzed RNA-Seq and ChIP-Seq data and drafted the manuscript. LOV designed initial approach, obtained funding, led and coordinated the project, and edited the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Complete list of candidate genes potentially regulated by NAC transcription factor during soybean seedling development.(XLSX 15 KB)
Additional file 4: Complete list of candidate genes potentially regulated by YABBY transcription factor during soybean seedling development.(XLSX 17 KB)
Additional file 5: Multiple sequence alignment of four members of NAC family expressed during the functional transition of the cotyledons.(PDF 207 KB)
Additional file 6: Multiple sequence alignment of ten members of YABBY family expressed during the functional transition of the cotyledons.(PDF 226 KB)
Additional file 7: List of 53 gene models encoding WRKY transcription factors expressed at level ( ≥ 10 RPKM) in at least one of seven stages of soybean seedling development.(XLSX 26 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Shamimuzzaman, M., Vodkin, L. Genome-wide identification of binding sites for NAC and YABBY transcription factors and co-regulated genes during soybean seedling development by ChIP-Seq and RNA-Seq. BMC Genomics 14, 477 (2013). https://doi.org/10.1186/1471-2164-14-477
- Transcription Factor Binding Site
- Functional Transition
- WRKY Transcription Factor
- Yabby Family
- Yabby Transcription Factor