A valid strategy for precise identifications of transcription factor binding sites in combinatorial regulation using bioinformatic and experimental approaches
© Wang et al.; licensee BioMed Central Ltd. 2013
Received: 1 April 2013
Accepted: 13 August 2013
Published: 24 August 2013
Transcription factor (TF) binding sites (cis element) play a central role in gene regulation, and eukaryotic organisms frequently adapt a combinatorial regulation to render sophisticated local gene expression patterns. Knowing the precise cis element on a distal promoter is a prerequisite for studying a typical transcription process; however, identifications of cis elements have lagged behind those of their associated trans acting TFs due to technical difficulties. Consequently, gene regulations via combinatorial TFs, as widely observed across biological processes, have remained vague in many cases.
We present here a valid strategy for identifying cis elements in combinatorial TF regulations. It consists of bioinformatic searches of available databases to generate candidate cis elements and tests of the candidates using improved experimental assays. Taking the MYB and the bHLH that collaboratively regulate the anthocyanin pathway genes as examples, we demonstrate how candidate cis motifs for the TFs are found on multi-specific promoters of chalcone synthase (CHS) genes, and how to experimentally test the candidate sites by designing DNA fragments hosting the candidate motifs based on a known promoter (us1 allele of Ipomoea purpurea CHS-D in our case) and applying site-mutagenesis at the motifs. It was shown that TF-DNA interactions could be unambiguously analyzed by assays of electrophoretic mobility shift (EMSA) and dual-luciferase transient expressions, and the resulting evidence precisely delineated a cis element. The cis element for R2R3 MYBs including Ipomoea MYB1 and Magnolia MYB1, for instance, was found to be ANCNACC, and that for bHLHs (exemplified by Ipomoea bHLH2 and petunia AN1) was CACNNG. A re-analysis was conducted on previously reported promoter segments recognized by maize C1 and apple MYB10, which indicated that cis elements similar to ANCNACC were indeed present on these segments, and tested positive for their bindings to Ipomoea MYB1.
Identification of cis elements in combinatorial regulation is now feasible with the strategy outlined. The working pipeline integrates the existing databases with experimental techniques, providing an open framework for precisely identifying cis elements. This strategy is widely applicable to various biological systems, and may enhance future analyses on gene regulation.
Keywordscis element MYB bHLH EMSA Dual-luciferase transient expression assay
The transcription process can be highly dynamic and sophisticated in eukaryotic cells, and its initiation typically involves a recognition between a transcription factor (TF) and a cis element at the upstream region on a gene. Knowing the precise sequence of a cis element for its associated TF is therefore critical for understanding the transcription process of a given gene. The coupling process, however, has not been understood in a balanced way, as cis elements have not been clarified for many well characterized trans acting TFs.
Both technical and inherent difficulties appear to contribute to the situation. In the previous efforts devoted to cis identification, a widely taken strategy was segment-dissection. Typically, the 5’ promoter region of a gene was cut by segments and examined how such manipulations affected the expression of the gene, then an inference was made on which segment was likely responsible for the transcription initiation. The protocol can be laborious and frequently result in only rough estimates of cis elements in simple cases or inconclusiveness in cases of multiple TFs (such as those in the combinatorial regulation). When the correct segments are identified, their lengths can be inconsistent between reports due to uncertainties about the exact binding sites.
The situation may become inherently more chaotic in some cases. The DNA binding of a TF per se may not require a full participation of all nucleotides within the binding region (previously referred to as gapped or degenerate cis element) on a promoter, causing more or less varied sequence content and consequently an increased difficulty for identifying the relevant cis element. The bordering nucleotide sites around a cis element may also influence TF-DNA binding without being part of the TF-DNA complex. More subtly, different species’ versions of the same member of a TF family may recognize somewhat varied cis elements due to evolution of the TF-DNA interaction.
It seems lack of precise identifications of cis elements cannot be circumvented by the existing approaches. A recent survey on MYB cis elements shows that rich experimental data has been collected on MYBs across kingdoms, and many cis elements have been reported, but MYB-DNA interactions remain vague. In plants, most MYBs have R2 and R3 domains, and the cis elements for R2R3 MYBs have been reported in various lengths from 4 (AACA) in the case of rice OSMYB5 by DNase I footprinting analyses to 14 nucleotides (TAT AAC GGT TTT TT) in that of soybean GmMYBs by yeast one hybrid. While these reports do not show coherence in the length of R2R3 binding region, there is no reason to believe that the binding sites for plant R2R3 MYBs should vary greatly in length. In contrast, the R2R3 domains of mouse c-MYB had been shown to bind to AACNG via the nuclear magnetic resonance, and the crystallization of protozoan tv MYB1-DNA suggested the binding sequence to be a/gACGAT. In other words, the range of R2R3 MYB binding sites has been poorly defined in planta, causing problems for further analysis. Apparently, a more focused identification method is desired for pinpointing the cis element for a given TF.
In an analysis of the combinatorial regulation on the flavonoid network, we have developed a strategy that is highly effective in identifications of cis elements, particularly for plant systems. The strategy features bioinformatically generated cis candidates and their effective validations via experimental means based on some of the existing protocols. The initial step places a bioinformatic analysis in the forefront, which mines candidate cis sites species- and locus-wide (assuming a tractable cis element). It takes advantage of partially known cis information for one of the TFs involved in combinatorial regulation in order to predict the likely cis element for another TF in interest. For instance, the binding sites for bHLHs have been known for some members of the TF family. Both mouse c-myc and Brassica bHLH recognize CACGTG, which was initially known as the core of the G-box (−TCTTACACGTGGCAYY-) on the promoter of a small subunit of ribulose 1, 5-bisphosphate carboxylase gene. Since the regulation of the anthocyanin pathway requires both MYB and bHLH (there is evidence for the bHLH’s binding to CACGTG), this attribute of bHLH may serve as an anchor for obtaining the cis sites for MYB via our bioinformatic approaches.
Regulation by combinatorial TFs features many biological processes. Here we focused on better known MYB, bHLH and WD-repeat protein (WDR) families in plants. These three kinds of TFs may form a complex to influence trichome formation and proanthocyanidin synthesis, in additional to anthocyanin synthesis. As part of the flavonoid network, the components of anthocyanin pathway are relatively well-defined, mainly consisting of chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3- dioxygenase (F3H), flavonoid 3’-monooxygenase (F3’H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), and UDP-glucose flavonoid 3-O-glucosyltransferase (3GT). These supposed regulons and their TFs have been identified in several species. The MYB-bHLH-WDR complex have been known to include C1- B/R- PAC1 in maize[20, 21], AN2 - AN1 - AN11 in petunia[22, 23], and MYB1 - bHLH2 - WDR1 in Ipomoea[24–26], respectively. Meanwhile, MYB and bHLH (but not WDR) have been known to interact with promoters in order to fulfill their roles in gene regulations. The anthocyanin pathway system is thus an ideal system for examining TF-DNA interactions.
We take CHS as an example here, showing that once candidate cis motifs are generated from the bioinformatic analysis, their validations may be effectively tested through site-directed mutagenesis and experiments targeting specific TF-DNA interactions. Numerous tests have suggested that electrophoretic mobility shift assays (EMSAs)[27, 28] and transient expression assays using living cells[29, 30] are highly effective in TF-DNA interaction analysis. For EMSA, we have obtained the best resolution with commercialized fluorescent dyes including SYBR® Green and SYPRO® Ruby (Molecular probes/Life technologies), which may bind to nucleotide and protein, respectively. Sequential applications of these dyes to the same gel and exposures of the gel under different light conditions lead to detections of unambiguous signals of DNA-protein interactions. For dual-luciferase transient expression assays, we have made promoter constructs with desired site-mutations, and engaged particle bombardment and transient gene expressions in living-cells for analyzing candidate motifs. The complete working pipeline is detailed below for precise cis identifications on genes under combinatorial regulation.
Results and discussion
Bioinformatic searches for candidate cis motifs
Results of promoter searches using the GenBank data on plants
Number of promoter sequences
bHLH2 _bh2b: gi|97974090|dbj|BAE94388.1
It was no coincidence that the CHS promoters in our search presented the Motifs most readily, as available CHS promoters greatly outnumbered those of other anthocyanin genes in our dataset (Table 1). Clearly, a large locus-wide promoter sample size may increase the odds of detecting candidates for cis elements. The strategy of taking bHLH binding sites as bait is not restricted to the anthocyanin pathway. Many members of the bHLH family are known to interact with other TFs, and the bioinformatic search protocol (Additional file1) is also applicable to cases of TFs with other known cis elements. Compared to earlier bioinformatic treatments on cis elements enriched via various approaches such as the case regarding Arabidopsis anthocyanin regulator PAP1, our analysis is much easier to interpret since a clear biological function is associated with the locus-based searches. More importantly, the candidate cis motifs by our protocol do not depend on any prior knowledge of the cis to be searched, which distinguishes itself from other bioinformatic methods.
Validation experiments on the cis element for anthocyanin bHLH
Since Motif 1 appeared on target for the MYB1-bHLH2-WDR1 combinatorial regulation in living cells, we further examined Motif 1 in EMSAs to make sure that it was bHLH2 that specifically bound to the cis sites. To this end, we designed series of probes from the CHS-D us1 promoter region hosting Motif 1. Besides the local sequence, probes also included mutated sequences at Motif 1 and those at the 5’ and 3’ bordering sites (Figure 2C). We prepared the bHLH2 protein with the pMAL system (NewEngland BioLabs), which used maltose-binding protein (MBP) to tag a protein (Additional file4). The combined protein was then separated through an amylose resin-based column. Since a whole bHLH2 could not be properly expressed and folded in the E. coli system, its partial expression was adopted to include the intact bHLH binding domain and the immediately adjacent 3’ conserved region. Including the 3’ conserved region has been shown effective previously and in our tests. To increase the binding signals, target mutations were duplicated twice on the probes. The duplications had been examined in our prior tests to make sure it enhanced signal intensities only without compromising the binding patterns (data not shown). After each probe was exposed to the bHLH2, we examined the TF-DNA interaction on a nondenaturing polyacrylamide gel. If the probe binds to the protein, the TF-DNA complex will move slower than the free probe and protein inside the gel, causing a slower band visible under light for viewing a DNA or protein staining. Since a positive binding can be detected by independent fluorescent signals on the same gel, the evidence for a DNA-protein interaction is more convincing than one of the competition design in previous EMSAs. Here, we observed a positive binding of the MBP-bHLH2 to a probe containing quadruple CACGTG (Figure 2D). Mutations at the 5’ bordering site and the first three and last nucleotides of CACGTG significantly reduced the binding activity while mutations at the fourth and fifth nucleotides of CACGTG and 3’ bordering site did not significantly affect the binding pattern. To make sure of the results, we also tested petunia AN1 (AF260919) and witnessed the same pattern (Figure 2E). The G-box sequence (CACGTG) was hence considered the bHLH recognition element (BRE) for the tested TFs.
In summary, both EMSA and transient expression experiments agree that BREs for the anthocyanin pathway bHLH take the form of CACNNG(T). For the cis structure, the fourth (allowing G → A but not G → C) and fifth positions may tolerate certain changes, and the last site prefers G or T. Here, we did not examine all likely substitutions at the 5’ adjacent site to CACGTG since it may be specific to an interacting TF, as we will see later with the case of MYB-DNA interactions.
Verification and delineation of cis motif for anthocyanin MYB
Transient expression experiments were also examined for their effectiveness in delineating and verifying cis elements. We constructed reporters with the promoters mutated around the MRE and performed the living cell tests using pCHSD-337 (Figure 3C). The results confirmed that mutation at the 5’ bordering site did not affect the promoter activity, agreeing with the previous EMSAs on the 5’ border of the cis element. Our delineated cis element (AGCTACC) resembles ACCTACCAACC in maize at the 5’ side, and agrees with the middle portion of ACTGGTAGCTATT in apple in the reverse complement (GGTAGCT). These previously inconsistent results now converge with our MRE due to the precise identifications presented here.
The IpMYB1 based analysis was examined again in a different species to know whether or not IpMYB1 homologues can recognize the cis element identified in Ipomoea. We cloned MlMYB1 (KC794950) from Mulan magnolia (Magnolia liliflora), whose amino acid sequence is 49% identical to that of IpMYB1. When the two MYBs were independently incubated with the same probes based on the promoter of IpMYB1 (KC794943), we found essentially the same binding pattern with only varied binding intensities (Figure 3D). MlMYB1prefered the first probe more than the third one, though the two probes differed only at the nucleotide adjacent to the MRE, while IpMYB1 seemed to favor the third probe. The results suggest that the bordering site of a cis element may cause different DNA-binding affinities between TFs.
MYB-DNA interactions on previous anthocyanin promoters
By taking the combinatorial regulation on the anthocyanin pathway as an example, we have demonstrated the validity of the whole strategy for precisely identifying cis elements on a relevant promoter. It is evidently not species- or tissues- specific. Promoters of anthocyanin genes (CHS, MYB1, DFR, 3GT) across species (I. purpurea, Zea mays, Malus domestica, Gerbera hybrida) have been successfully scrutinized for their cis elements. TFs across species (I. purpurea, M. liliflora, P. hybrida) have been tested for their binding capabilities. Among the previous tools such as crystallization, DNA footprinting, yeast one-hybrid, EMSA, and transient expressions in cis research, we found the combination of site-directed mutagenesis with enhanced EMSA protocol to be the most cost-effective and time-efficient in defining a cis element. As the evidence gathered from the EMSA was for physical binding outside the cell, we devised the second system with mutagenesis/particle bombardment/transient expression to collect relevant evidence inside living cells. Obviously, further species-specific in vivo approaches including transgenic techniques can be integrated to strengthen the cis analysis. The strategy outlined here is versatile enough to accommodate additional techniques of the future.
With ever increasing genomic data and promoter sequences, bioinformatic searches will offer an indispensable summary of the available data and prove effective in establishing candidate motifs for tractable cis elements, about which no prior information is required. The hypothetical motifs can be experimentally tested and delineated with properly designed site-directed mutagenesis and in vitro and in vivo approaches. This logic is broadly applicable to biological systems. The application will hopefully facilitate both gene functional analysis and gene regulation at large.
We performed a blastp search on the Genbank database (http://www.ftp.ncbi.nih.gov/genbank/, up to 10 April 2009) to assemble protein homologs for each of the ten anthocyanin gene templates, using the criteria of sequence identity > 40% and length coverage > 80%. The searches led to a collection of the DNA sequences (n = 571) that were at least 200 bp from the ATG start codon at the 5’ region. We then applied clustalw (http://www.ebi.ac.uk/Tools/msa/clustalo/) to the sequences and only picked the longer sequence for each of the clusters that scored higher than 80% in identity. This process resulted in a non-redundant sample (n = 271) for the subsequent data analysis. Assuming that the binding site of bHLH2 took the form of CACGTG and allowed one mutation at any of the sites on the collected promoters, we searched on both directions for motifs of 4-nucleotides (XXXX) that satisfied CACGTG(…)XXXX under the conditions that (1) the occurrence frequency of each combination > 4 and (2) the distance between CACGTG and XXXX ≤ 50 bp. Since a given combination of CACGTG and XXXX may have slightly varied separating distances, we classified them using the nearest neighbour method to reduce scattering of a particular combination. For each of the classified combinations, we dissected out 100 bp promoter sequences flanking the sites (Additional file1), and prepared the dataset in fasta format for MEME (http://www.meme.nbcr.net/meme/, version 4.3.0). The parameters for MEME were set as: occurrences of a single motif among the sequences: any number; minimum width of each motif: 6 bp; maximum width of each motif: 10 bp; Maximum number of motifs to find: 5; Search: given strand only. Results for each combination of CACGTG(…)XXXX or XXXX(…)CACGTG were visually checked for likely patterns of motifs, which led to patterns of CACGTG(6–20)CTAC found on multiple promoters.
The common morning glory (Ipomoea purpurea) seeds were planted in the fields of the Botanical Garden of the Institute of Botany, Chinese Academy of Sciences, in Beijing, to provide petals and young leaves for TF and promoter amplifications and cloning. Mulan magnolia (Magnolia liliflora) inside the botanical garden provided floral tissue for MYB cloning. Petals of petunia R27 (kind gift from Dr. Quattrocchio) in a growth chamber were used for AN1 cloning. Plants from the ivory seeds of Japanese morning glory (Ipomoea nil, strain baichou) purchased locally provided leaf materials for the transient expression assays. For consistency, the assays were conducted on only the second young leaf (about 3 cm in diameter) of each plant stem in a growth chamber.
Promoter sequencing and cloning
The promoter of CHS-D (AF358654-9) was obtained from leaf genomic DNA via standard PCRs. The promoter of MYB1 was isolated via thermal asymmetric interlaced - PCR by specific and random primers. Three specific primers were used for obtaining the IpMYB1 promoter, including ltMyb + 131 (gctctaaagggaactaggtgccat), ltMyb + 56 (acgatggactccagtccgaccaagcaccctttctcact), and ltMyb + 30 (cggagaccaccttgcagaagaattaacc).
The coding sequences of MYB1, bHLH2, WDR1, firefly luciferase gene and renilla luciferase gene were integrated into pJIT163 at appropriate restriction sites (Additional file3). Plant transformations by particle bombardment (BioRad PDS-1000/He equipment) were performed on young leaves of the WDR1 mutant with the abaxial side facing up on the solidified water medium (0.7% agar). Effector and reporter plasmid DNAs each in 1 mg were mixed with 0.2 mg of reference plasmid pJIT163-35S: RUC. Microparticles (1.0 μm, Bio-Rad) were prepared at the concentration of 50 mg/ml, and mixed with prepared plasmids at 2.0 μl microparticles per 1.0 μg plasmid DNA. The particle-DNA complex was then mixed with 2.5 M CaCl2 and 0.1 M spermidine, and centrifuged to remove the supernatant. After washing in 70% and 100% ethanol successively, the particle-DNA complex was re-suspended in appropriate amount of 100% ethanol (10 μl for one shot). The leaves were subjected to bombardment at 1100 psi and 26 mm Hg vacuum, and kept in dark at 24 C for about 24 hrs. The luciferase activities were detected by a Glomax 20/20 luminometer following the manufactory’s instruction on Dual-Luciferase® Reporter (DLR™) Assay System (Promega).
Protein expression and purification
Constructs were built from pMAL-c2G (Additional file4), as instructed for the pMAL™ protein fusion and purification system (New England Biolabs). The host strain E. coli BL21 (DE3) was transformed by the target vector and grown at 37 C in LB solution with ampicillin (50 μg/ml) and chloramphenicol (50 μg/ml), and added with IPTG to 0.5 mM when the solution reached 0.5 at OD600, then cultured under 25 C for six hours before protein extraction. One liter of the cultured cell solution was cleaned twice with the loading buffer then ultrasonicated at 400 W in pulses of 5 second for 5 min in ice and centrifuged at 4 C at 12000 rpm for 15 min to collect the supernatant. Affinity-based purification of the proteins was performed on an amylose resin column following the standard protocol, and the separated samples were examined with 12% SDS-PAGE electrophoresis to pick ones with the desired protein expression. Qualified samples were quantified via Bradford method using standard bovine gamma globulin following the manufactory’s instructions on Bio-Rad protein assays.
Electrophoretic mobility shift assays (EMSA)
Two complementary oligonucleotides synthesized by Sangon (Sangon Biotech) were annealed at 95°C for 10 min in 0.2 × SSC buffer (3 mM sodium citrate, 30 mM NaCl, pH 7.0) then slowly cooled in three hrs to room temperature to make a probe. Each probe was quantified with known quantity of λDNA (Life Technologies) using Picogreen in the quantification module of Rotor Gene 3000 (Corbett Research).
A 10 μl binding reaction mixture, containing 2 μl 5× binding buffer [25% glycerol, 40 mM MgCl2, 5 mM EDTA, 5 mM DTT, 250 mM NaCl, and 50 mM Tris–HCl (pH 7.5)], 20 pmol probe, and 2 μg MBP- target proteins, was incubated at 25 C for 10 min. Samples loaded to a 12% (wt/vol) nondenaturing polyacrylamide gel were subject to electrophoresis in 0.25X TBE buffer at 130 V for 1 h. The gel was covered in 1XSYBR Green EMSA stain (Life Technologies) for 20 min at room temperature in dark and washed twice with dd H2O, then visualized for the presence of DNA using a G-Box ChemiXL Fluorescent Imager (Syngene) under 300 nm excitation and 500 - 600 nm emission. The gel was further treated in SYPRO Ruby EMSA stain (Life Technologies) and destained in 10% methanol and 7% acetic acid overnight. Before acquiring its image, the gel was washed twice, then detected for the presence of protein under 300 nm excitation and 550 to 640 nm (EtBr/UV Mid pass filters) emission settings.
We thank the reviewers for constructive comments, the groups of Dr. Kang Chong and Dr. Hongya Gu for equipment access. The study was supported by the National Science Foundation of China (31070263, 30770176), the Chinese Academy of Sciences (KSCX2-YW-N-043), and the Ministry of Science and Technology.
- Lelli KM, Slattery M, Mann RS: Disentangling the many layers of eukaryotic transcriptional regulation. Annu Rev Genet. 2012, 46: 43-68. 10.1146/annurev-genet-110711-155437.PubMed CentralView ArticlePubMedGoogle Scholar
- Walhout AJM: Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping. Genome Res. 2006, 16: 1445-1454. 10.1101/gr.5321506.View ArticlePubMedGoogle Scholar
- Elnitski L, Jin VX, Farnham PJ, Jones SJM: Locating mammalian transcription factor binding sites: a survey of computational and experimental techniques. Genome Res. 2006, 16: 1455-1464. 10.1101/gr.4140006.View ArticlePubMedGoogle Scholar
- Hannenhalli S: Eukaryotic transcription factor binding sites - modeling and integrative search methods. Bioinformatics. 2008, 24: 1325-1331. 10.1093/bioinformatics/btn198.View ArticlePubMedGoogle Scholar
- Prouse MB, Campbell MM: The interaction between MYB proteins and their target DNA binding sites. Biochimica et biophysica acta. 1819, 2012: 67-77.Google Scholar
- Stracke R, Werber M, Weisshaar B: The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2001, 4: 447-456. 10.1016/S1369-5266(00)00199-0.View ArticlePubMedGoogle Scholar
- Suzuki A, Wu CY, Washida H, Takaiwa F: Rice MYB protein OSMYB5 specifically binds to the AACA motif conserved among promoters of genes for storage protein glutelin. Plant and Cell Physiology. 1998, 39: 555-559. 10.1093/oxfordjournals.pcp.a029404.View ArticlePubMedGoogle Scholar
- Liao Y, Zou HF, Wang HW, Zhang WK, Ma B, Zhang JS, Chen SY: Soybean GmMYB76, GmMYB92, and GmMYB177 genes confer stress tolerance in transgenic Arabidopsis plants. Cell Res. 2008, 18: 1047-1060. 10.1038/cr.2008.280.View ArticlePubMedGoogle Scholar
- Ogata K, Morikawa S, Nakamura H, Sekikawa A, Inoue T, Kanai H, Sarai A, Ishii S, Nishimura Y: Solution structure of a specific DNA complex of the myb DNA-binding domain with cooperative recognition helices. Cell. 1994, 79: 639-648. 10.1016/0092-8674(94)90549-5.View ArticlePubMedGoogle Scholar
- Jiang I, Tsai CK, Chen SC, Wang SH, Amiraslanov I, Chang CF, Wu WJ, Tai JH, Liaw YC, Huang TH: Molecular basis of the recognition of the ap65-1 gene transcription promoter elements by a Myb protein from the protozoan parasite Trichomonas vaginalis. Nucleic Acids Res. 2011, 39: 8992-9008. 10.1093/nar/gkr558.PubMed CentralView ArticlePubMedGoogle Scholar
- Blackwell TK, Kretzner L, Blackwood EM, Eisenman RN, Weintraub H: Sequence-specific dna-binding by the c-myc protein. Science. 1990, 250: 1149-1151. 10.1126/science.2251503.View ArticlePubMedGoogle Scholar
- Waldmuller S, Link G: Sequence and expression characteristics of three g-Box-binding factor cDNAs from Brassica napus. Plant Mol Biol. 1995, 29: 507-517. 10.1007/BF00020981.View ArticlePubMedGoogle Scholar
- Giuliano G, Hoffman NE, Ko K, Scolnik PA, Cashmore AR: A light-entrained circadian clock controls transcription of several plant genes. Embo Journal. 1988, 7: 3635-3642.PubMed CentralPubMedGoogle Scholar
- Kong Q, Pattanaik S, Feller A, Werkman JR, Chai CL, Wang YQ, Grotewold E, Yuan L: Regulatory switch enforced by basic helix-loop-helix and ACT-domain mediated dimerizations of the maize transcription factor R. Proc Natl Acad Sci U S A. 2012, 109: E2091-E2097. 10.1073/pnas.1205513109.PubMed CentralView ArticlePubMedGoogle Scholar
- Feller A, Machemer K, Braun EL, Grotewold E: Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant Journal. 2011, 66: 94-116. 10.1111/j.1365-313X.2010.04459.x.View ArticlePubMedGoogle Scholar
- Payne CT, Zhang F, Lloyd AM: GL3 encodes a bHLH protein that regulates trichome development in Arabidopsis through interaction with GL1 and TTG1. Genetics. 2000, 156: 1349-1362.PubMed CentralPubMedGoogle Scholar
- Baudry A, Heim MA, Dubreucq B, Caboche M, Weisshaar B, Lepiniec L: TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis thaliana. Plant Journal. 2004, 39: 366-380. 10.1111/j.1365-313X.2004.02138.x.View ArticlePubMedGoogle Scholar
- Koes R, Verweij W, Quattrocchio F: Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005, 10: 236-242. 10.1016/j.tplants.2005.03.002.View ArticlePubMedGoogle Scholar
- Winkel-Shirley B: It takes a garden. How work on diverse plant species has contributed to an understanding of flavonoid metabolism. Plant Physiol. 2001, 127: 1399-1404. 10.1104/pp.010675.PubMed CentralView ArticlePubMedGoogle Scholar
- Goff SA, Cone KC, Chandler VL: Functional-analysis of the transcriptional activator encoded by the maize-B gene - evidence for a direct functional interaction between 2 classes of regulatory proteins. Genes Dev. 1992, 6: 864-875. 10.1101/gad.6.5.864.View ArticlePubMedGoogle Scholar
- Carey CC, Strahle JT, Selinger DA, Chandler VL: Mutations in the pale aleurone color1 regulatory gene of the Zea mays anthocyanin pathway have distinct phenotypes relative to the functionally similar TRANSPARENT TESTA GLABRA1 gene in Arabidopsis thaliana. Plant Cell. 2004, 16: 450-464. 10.1105/tpc.018796.PubMed CentralView ArticlePubMedGoogle Scholar
- Quattrocchio F, Wing JF, Leppen HTC, Mol JNM, Koes RE: Regulatory genes-controlling anthocyanin pigmentation are functionally conserved among plant-species and have distinct sets of target genes. Plant Cell. 1993, 5: 1497-1512.PubMed CentralView ArticlePubMedGoogle Scholar
- deVetten N, Quattrocchio F, Mol J, Koes R: The an11 locus controlling flower pigmentation in petunia encodes a novel WD-repeat protein conserved in yeast, plants, and animals. Genes Dev. 1997, 11: 1422-1434. 10.1101/gad.11.11.1422.View ArticleGoogle Scholar
- Chang SM, Lu YQ, Rausher MD: Neutral evolution of the nonbinding region of the anthocyanin regulatory gene Ipmyb1 in Ipomoea. Genetics. 2005, 170: 1967-1978. 10.1534/genetics.104.034975.PubMed CentralView ArticlePubMedGoogle Scholar
- Park KI, Ishikawa N, Morita Y, Choi JD, Hoshino A, Iida S: A bHLH regulatory gene in the common morning glory, Ipomoea purpurea, controls anthocyanin biosynthesis in flowers, proanthocyanidin and phytomelanin pigmentation in seeds, and seed trichome formation. Plant Journal. 2007, 49: 641-654. 10.1111/j.1365-313X.2006.02988.x.View ArticlePubMedGoogle Scholar
- Morita Y, Saitoh M, Hoshino A, Nitasaka E, Iida S: Isolation of cDNAs for R2R3-MYB, bHLH and WDR transcriptional regulators and identification of c and ca mutations conferring white flowers in the Japanese morning glory. Plant and Cell Physiology. 2006, 47: 457-470. 10.1093/pcp/pcj012.View ArticlePubMedGoogle Scholar
- Fried M, Crothers DM: Equilibria and kinetics of lac repressor-operator interactions by polyacrylamide-gel electrophoresis. Nucleic Acids Res. 1981, 9: 6505-6525. 10.1093/nar/9.23.6505.PubMed CentralView ArticlePubMedGoogle Scholar
- Hellman LM, Fried MG: Electrophoretic mobility shift assay (EMSA) for detecting protein–nucleic acid interactions. Nat Protoc. 2007, 2: 1849-1861. 10.1038/nprot.2007.249.PubMed CentralView ArticlePubMedGoogle Scholar
- Ow DW, Wood KV, Deluca M, Dewet JR, Helinski DR, Howell SH: Transient and stable expression of the firefly luciferase gene in plant-cells and transgenic plants. Science. 1986, 234: 856-859. 10.1126/science.234.4778.856.View ArticlePubMedGoogle Scholar
- McNabb DS, Reed R, Marciniak RA: Dual luciferase assay system for rapid assessment of gene expression in Saccharomyces cerevisiae. Eukaryot Cell. 2005, 4: 1539-1549. 10.1128/EC.4.9.1539-1549.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren JY, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, 37: W202-W208. 10.1093/nar/gkp335.PubMed CentralView ArticlePubMedGoogle Scholar
- Toledo-Ortiz G, Huq E, Quail PH: The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003, 15: 1749-1770. 10.1105/tpc.013839.PubMed CentralView ArticlePubMedGoogle Scholar
- Dare AP, Schaffer RJ, Lin-Wang K, Allan AC, Hellens RP: Identification of a cis-regulatory element by transient analysis of co-ordinately regulated genes. Plant Methods. 2008, 4: 17-10.1186/1746-4811-4-17.PubMed CentralView ArticlePubMedGoogle Scholar
- Roth BA, Goff SA, Klein TM, Fromm ME: C1-dependent and r-dependent expression of the maize bz1 gene requires sequences with homology to mammalian myb and myc binding-sites. Plant Cell. 1991, 3: 317-325.PubMed CentralView ArticlePubMedGoogle Scholar
- Sainz MB, Grotewold E, Chandler VL: Evidence for direct activation of an anthocyanin promoter by the maize C1 protein and comparison of DNA binding by related Myb domain proteins. Plant Cell. 1997, 9: 611-625.PubMed CentralView ArticlePubMedGoogle Scholar
- Lesnick ML, Chandler VL: Activation of the maize anthocyanin gene a2 is mediated by an element conserved in many anthocyanin promoters. Plant Physiol. 1998, 117: 437-445. 10.1104/pp.117.2.437.PubMed CentralView ArticlePubMedGoogle Scholar
- Espley RV, Brendolise C, Chagne D, Kutty-Amma S, Green S, Volz R, Putterill J, Schouten HJ, Gardiner SE, Hellens RP, Allan AC: Multiple repeats of a promoter segment causes transcription factor autoregulation in red apples. Plant Cell. 2009, 21: 168-183. 10.1105/tpc.108.059329.PubMed CentralView ArticlePubMedGoogle Scholar
- Elomaa P, Uimari A, Mehto M, Albert VA, Laitinen RAE, Teeri TH: Activation of anthocyanin biosynthesis in Gerbera hybrida (Asteraceae) suggests conserved protein-protein and protein-promoter interactions between the anciently diverged monocots and eudicots. Plant Physiol. 2003, 133: 1831-1842. 10.1104/pp.103.026039.PubMed CentralView ArticlePubMedGoogle Scholar
- Cover T, Hart P: Nearest neighbor pattern classification. IEEE Transactions on Information Theory. 1967, 13: 21-27.View ArticleGoogle Scholar
- Liu YG, Whittier RF: Thermal asymmetric interlaced PCR - automatable amplification and sequencing of insert end fragments from P1 and YAC clones for chromosome walking. Genomics. 1995, 25: 674-681. 10.1016/0888-7543(95)80010-J.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://www.creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.