DNA binding properties of two Arabidopsis MADS domain proteins: binding consensus and dimer formation

Dual specificity and target gene selection by the MADS-domain protein FRUITFULL

How transcription factors attain their target gene specificity and how this specificity may be modulated, acquiring different regulatory functions through the development of plant tissues, is an open question. Here we characterized different regulatory roles of the MADS-domain transcription factor FRUITFULL (FUL) in flower development and mechanisms modulating its activity. We found that the dual role of FUL in regulating floral transition and pistil development is associated with its different in vivo patterns of DNA binding in both tissues. Characterization of FUL protein complexes by liquid chromatography-tandem mass spectrometry and SELEX-seq experiments shows that aspects of tissue-specific target site selection can be predicted by tissue-specific variation in the composition of FUL protein complexes with different DNA binding specificities, without considering the chromatin status of the target region. This suggests a role for dynamic changes in FUL TF complex composition in reshaping the regulatory functions of FUL during flower development.

The identification of type I MADS box genes as the upstream activators of an endosperm-specific invertase inhibitor in Arabidopsis

BACKGROUND

Nuclear endosperm development is a common mechanism among Angiosperms, including Arabidopsis. During nuclear development, the endosperm nuclei divide rapidly after fertilization without cytokinesis to enter the syncytial phase, which is then followed by the cellularized phase. The endosperm can be divided into three spatial domains with distinct functions: the micropylar, peripheral, and chalazal domains. Previously, we identified two putative small invertase inhibitors, InvINH1 and InvINH2, that are specifically expressed in the micropylar region of the syncytial endosperm. In addition, ectopically expressing InvINH1 in the cellularized endosperm led to a reduction in embryo growth rate. However, it is not clear what are the upstream regulators responsible for the specific expression of InvINHs in the syncytial endosperm.

RESULTS

Using protoplast transient expression system, we discovered that a group of type I MADS box transcription factors can form dimers to activate InvINH1 promoter. Promoter deletion assays carried out in the protoplast system revealed the presence of an enhancer region in InvINH1 promoter, which contains several consensus cis-elements for the MADS box proteins. Using promoter deletion assay in planta, we further demonstrated that this enhancer region is required for InvINH1 expression in the syncytial endosperm. One of the MADS box genes, AGL62, is a key transcription factor required for syncytial endosperm development. Using promoter-GFP reporter assay, we demonstrated that InvINH1 and InvINH2 are not expressed in agl62 mutant seeds. Collectively, our data supports the role of AGL62 and other type I MADS box genes as the upstream activators of InvINHs expression in the syncytial endosperm.

CONCLUSIONS

Our findings revealed several type I MADS box genes that are responsible for activating InvINH1 in the syncytial endosperm, which in turn regulates embryo growth rate during early stage of seed development.

Ectopic expression of Triticum polonicum VRT-A2 underlies elongated glumes and grains in hexaploid wheat in a dosage-dependent manner

Flower development is an important determinant of grain yield in crops. In wheat (Triticum spp.), natural variation for the size of spikelet and floral organs is particularly evident in Triticum turgidum ssp. polonicum (also termed Triticum polonicum), a tetraploid subspecies of wheat with long glumes, lemmas, and grains. Using map-based cloning, we identified VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT2), which encodes a MADS-box transcription factor belonging to the SHORT VEGETATIVE PHASE family, as the gene underlying the T. polonicum long-glume (P1) locus. The causal P1 mutation is a sequence rearrangement in intron-1 that results in ectopic expression of the T. polonicum VRT-A2 allele. Based on allelic variation studies, we propose that the intron-1 mutation in VRT-A2 is the unique T. polonicum subspecies-defining polymorphism, which was later introduced into hexaploid wheat via natural hybridizations. Near-isogenic lines differing for the P1 locus revealed a gradient effect of P1 across spikelets and within florets. Transgenic lines of hexaploid wheat carrying the T. polonicum VRT-A2 allele show that expression levels of VRT-A2 are highly correlated with spike, glume, grain, and floral organ length. These results highlight how changes in expression profiles, through variation in cis-regulation, can affect agronomic traits in a dosage-dependent manner in polyploid crops.

Variation for heterodimerization and nuclear localization among known and novel oil palm SHELL alleles

Oil palm breeding involves crossing dura and pisifera palms to produce tenera progeny with greatly improved oil yield. Oil yield is controlled by variant alleles of a type II MADS-box gene, SHELL, that impact the presence and thickness of the endocarp, or shell, surrounding the fruit kernel. We identified six novel SHELL alleles in noncommercial African germplasm populations from the Malaysian Palm Oil Board. These populations provide extensive diversity to harness genetic, mechanistic and phenotypic variation associated with oil yield in a globally critical crop. We investigated phenotypes in heteroallelic combinations, as well as SHELL heterodimerization and subcellular localization by yeast two-hybrid, bimolecular fluorescence complementation and gene expression analyses. Four novel SHELL alleles were associated with fruit form phenotype. Candidate heterodimerization partners were identified, and interactions with EgSEP3 and subcellular localization were SHELL allele-specific. Our findings reveal allele-specific mechanisms by which variant SHELL alleles impact yield, as well as speculative insights into the potential role of SHELL in single-gene oil yield heterosis. Future field trials for combinability and introgression may further optimize yield and improve sustainability.

Overexpression of zmm28 increases maize grain yield in the field

Increasing maize grain yield has been a major focus of both plant breeding and genetic engineering to meet the global demand for food, feed, and industrial uses. We report that increasing and extending expression of a maize MADS-box transcription factor gene, zmm28, under the control of a moderate-constitutive maize promoter, results in maize plants with increased plant growth, photosynthesis capacity, and nitrogen utilization. Molecular and biochemical characterization of zmm28 transgenic plants demonstrated that their enhanced agronomic traits are associated with elevated plant carbon assimilation, nitrogen utilization, and plant growth. Overall, these positive attributes are associated with a significant increase in grain yield relative to wild-type controls that is consistent across years, environments, and elite germplasm backgrounds.

Rice transcription factor OsMADS57 regulates plant height by modulating gibberellin catabolism

BACKGROUND

The MADS-box transcription factors mainly function in floral organ organogenesis and identity specification. Few research on their roles in vegetative growth has been reported.

RESULTS

Here we investigated the functions of OsMADS57 in plant vegetative growth in rice (Oryza sativa). Knockdown of OsMADS57 reduced the plant height, internode elongation and panicle exsertion in rice plants. Further study showed that the cell length was remarkably reduced in the uppermost internode in OsMADS57 knockdown plants at maturity. Moreover, OsMADS57 knockdown plants were more sensitive to gibberellic acid (GA₃), and contained less bioactive GA₃ than wild-type plants, which implied that OsMADS57 is involved in gibberellin (GA) pathway. Expectedly, the transcript levels of OsGA2ox3, encoding GAs deactivated enzyme, were significantly enhanced in OsMADS57 knockdown plants. The level of EUI1 transcripts involved in GA deactivation was also increased in OsMADS57 knockdown plants. More importantly, dual-luciferase reporter assay and electrophoretic mobility shift assay showed that OsMADS57 directly regulates the transcription of OsGA2ox3 as well as EUI1 through binding to the CArG-box motifs in their promoter regions. In addition, OsMADS57 also modulated the expression of multiple genes involved in GA metabolism or GA signaling pathway, indicating the key and complex regulatory role of OsMADS57 in GA pathway in rice.

CONCLUSIONS

These results indicated that OsMADS57 acts as an important transcriptional regulator that regulates stem elongation and panicle exsertion in rice via GA-mediated regulatory pathway.

MADS-Box Genes Are Key Components of Genetic Regulatory Networks Involved in Abiotic Stress and Plastic Developmental Responses in Plants

Plants, as sessile organisms, adapt to different stressful conditions, such as drought, salinity, extreme temperatures, and nutrient deficiency, via plastic developmental and growth responses. Depending on the intensity and the developmental phase in which it is imposed, a stress condition may lead to a broad range of responses at the morphological, physiological, biochemical, and molecular levels. Transcription factors are key components of regulatory networks that integrate environmental cues and concert responses at the cellular level, including those that imply a stressful condition. Despite the fact that several studies have started to identify various members of the MADS-box gene family as important molecular components involved in different types of stress responses, we still lack an integrated view of their role in these processes. In this review, we analyze the function and regulation of MADS-box gene family members in response to drought, salt, cold, heat, and oxidative stress conditions in different developmental processes of several plants. In addition, we suggest that MADS-box genes are key components of gene regulatory networks involved in plant responses to stress and plant developmental plasticity in response to seasonal changes in environmental conditions.

OsMADS25 regulates root system development via auxin signalling in rice

The phytohormone auxin is essential for root development in plants. OsMADS25 is a homologue of the AGL17-clade MADS-box genes in rice. Despite recent progress, the molecular mechanisms underlying the regulation of root development by OsMADS25 are not well known. It is unclear whether OsMADS25 regulates root development via auxin signalling. In this study, we examined the role of OsMADS25 in root development and characterized the signalling pathway through which OsMADS25 regulates root system development in rice. OsMADS25 overexpression significantly increased, but RNAi gene silencing repressed primary root (PR) length and lateral root (LR) density. Moreover, OsMADS25 promoted LR development in response to NO₃^- . Further study showed that OsMADS25 increased auxin accumulation in the root system by enhancing auxin biosynthesis and transport, while also reducing auxin degradation, therefore stimulating root development. More importantly, OsMADS25 was found to regulate OsIAA14 expression directly by binding to the CArG-box in the promoter region of OsIAA14, which encodes an Aux/indole acetic acid (IAA) transcriptional repressor of auxin signalling. Elevated auxin levels and decreased OsIAA14 expression might lead to reduced OsIAA14 protein accumulation, as a mechanism to regulate auxin signalling. Therefore, our findings reveal a molecular mechanism by which OsMADS25 modulates root system growth and development in rice, at least partilly, via Aux/IAA-based auxin signalling.

The floral homeotic protein SEPALLATA3 recognizes target DNA sequences by shape readout involving a conserved arginine residue in the MADS-domain

SEPALLATA3 of Arabidopsis thaliana is a MADS-domain transcription factor (TF) and a key regulator of flower development. MADS-domain proteins bind to sequences termed 'CArG-boxes' [consensus 5'-CC(A/T)₆ GG-3']. Because only a fraction of the CArG-boxes in the Arabidopsis genome are bound by SEPALLATA3, more elaborate principles have to be discovered to better understand which features turn CArG-boxes into genuine recognition sites. Here, we investigate to what extent the shape of the DNA is involved in a 'shape readout' that contributes to the binding of SEPALLATA3. We determined in vitro binding affinities of SEPALLATA3 to DNA probes that all contain the CArG-box motif, but differ in their predicted DNA shape. We found that binding affinity correlates well with a narrow minor groove of the DNA. Substitution of canonical bases with non-standard bases supports the hypothesis of minor groove shape readout by SEPALLATA3. Analysis of mutant SEPALLATA3 proteins further revealed that a highly conserved arginine residue, which is expected to contact the DNA minor groove, contributes significantly to the shape readout. Our studies show that the specific recognition of cis-regulatory elements by a plant MADS-domain TF, and by inference probably also of other TFs of this type, heavily depends on shape readout mechanisms.

MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution

The floral quartet model of floral organ specification poses that different tetramers of MIKC-type MADS-domain transcription factors control gene expression and hence the identity of floral organs during development. Here, we provide a brief history of the floral quartet model and review several lines of recent evidence that support the model. We also describe how the model has been used in contemporary developmental and evolutionary biology to shed light on enigmatic topics such as the origin of land and flowering plants. Finally, we suggest a novel hypothesis describing how floral quartet-like complexes may interact with chromatin during target gene activation and repression.

Differences in DNA Binding Specificity of Floral Homeotic Protein Complexes Predict Organ-Specific Target Genes

Floral organ identities in plants are specified by the combinatorial action of homeotic master regulatory transcription factors. However, how these factors achieve their regulatory specificities is still largely unclear. Genome-wide in vivo DNA binding data show that homeotic MADS domain proteins recognize partly distinct genomic regions, suggesting that DNA binding specificity contributes to functional differences of homeotic protein complexes. We used in vitro systematic evolution of ligands by exponential enrichment followed by high-throughput DNA sequencing (SELEX-seq) on several floral MADS domain protein homo- and heterodimers to measure their DNA binding specificities. We show that specification of reproductive organs is associated with distinct binding preferences of a complex formed by SEPALLATA3 and AGAMOUS. Binding specificity is further modulated by different binding site spacing preferences. Combination of SELEX-seq and genome-wide DNA binding data allows differentiation between targets in specification of reproductive versus perianth organs in the flower. We validate the importance of DNA binding specificity for organ-specific gene regulation by modulating promoter activity through targeted mutagenesis. Our study shows that intrafamily protein interactions affect DNA binding specificity of floral MADS domain proteins. Differential DNA binding of MADS domain protein complexes plays a role in the specificity of target gene regulation.

Gains and Losses of Cis-regulatory Elements Led to Divergence of the Arabidopsis APETALA1 and CAULIFLOWER Duplicate Genes in the Time, Space, and Level of Expression and Regulation of One Paralog by the Other

How genes change their expression patterns over time is still poorly understood. Here, by conducting expression, functional, bioinformatic, and evolutionary analyses, we demonstrate that the differences between the Arabidopsis (Arabidopsis thaliana) APETALA1 (AP1) and CAULIFLOWER (CAL) duplicate genes in the time, space, and level of expression were determined by the presence or absence of functionally important transcription factor-binding sites (TFBSs) in regulatory regions. In particular, a CArG box, which is the autoregulatory site of AP1 that can also be bound by the CAL protein, is a key determinant of the expression differences. Because of the CArG box, AP1 is both autoregulated and cross-regulated (by AP1 and CAL, respectively), and its relatively high-level expression is maintained till to the late stages of sepal and petal development. The observation that the CArG box was gained recently further suggests that the autoregulation and cross-regulation of AP1, as well as its function in sepal and petal development, are derived features. By comparing the evolutionary histories of this and other TFBSs, we further indicate that the divergence of AP1 and CAL in regulatory regions has been markedly asymmetric and can be divided into several stages. Specifically, shortly after duplication, when AP1 happened to be the paralog that maintained the function of the ancestral gene, CAL experienced certain degrees of degenerate evolution, in which several functionally important TFBSs were lost. Later, when functional divergence allowed the survival of both paralogs, CAL remained largely unchanged in expression, whereas the functions of AP1 were gradually reinforced by gains of the CArG box and other TFBSs.

The VviMYB80 Gene is Abnormally Expressed in Vitis vinifera L. cv. 'Zhong Shan Hong' and its Expression in Tobacco Driven by the 35S Promoter Causes Male Sterility

Anther development is a very precise and complicated process. In Arabidopsis, the AtMYB80 transcription factor regulates genes involved in pollen development and controls the timing of tapetal programmed cell death (PCD). In this study, we isolated and characterized cDNA for VviMYB80 expressed in flower buds of male-sterile Vitis vinifera L. cv. 'Zhong Shan Hong', a late-maturing cultivar derived from self-progeny of cv. 'Wink'. VviMYB80 belongs to the MYB80 subfamily and clusters with AtMYB35/TDF1 in a distinct clade. We found that in flower buds, expression of the VviMYB80 gene in cv. 'Zhong Shan Hong' sharply increased at the tetrad stage, resulting in a higher and earlier transcript level than that found in cv. 'Wink'. Expression of the VviMYB80 gene, driven by the 35S promoter, caused pleiotropic effects on the stamens, including smaller and shriveled anthers, delayed dehiscence, fewer seeds, shorter anther filaments, distorted pollen shape and a lack of cytoplasm, with the tapetum exhibiting hypertrophy in transformed tobacco. These results suggest that VviMYB80 may play an important role in stamen development and that expression of VviMYB80 driven by the 35S promoter in tobacco induces male sterility.

SVP-like MADS-box protein from Carya cathayensis forms higher-order complexes

To properly regulate plant flowering time and construct floral pattern, MADS-domain containing transcription factors must form multimers including homo- and hetero-dimers. They are also active in forming hetero-higher-order complexes with three to five different molecules. However, it is not well known if a MADS-box protein can also form homo-higher-order complex. In this study a biochemical approach is utilized to provide insight into the complex formation for an SVP-like MADS-box protein cloned from hickory. The results indicated that the protein is a heterogeneous higher-order complex with the peak population containing over 20 monomers. Y2H verified the protein to form homo-complex in yeast cells. Western blot of the hickory floral bud sample revealed that the protein exists in higher-order polymers in native. Deletion assays indicated that the flexible C-terminal residues are mainly responsible for the higher-order polymer formation and the heterogeneity. Current results provide direct biochemical evidences for an active MADS-box protein to be a high order complex, much higher than a quartermeric polymer. Analysis suggests that a MADS-box subset may be able to self-assemble into large complexes, and thereby differentiate one subfamily from the other in a higher-order structural manner. Present result is a valuable supplement to the action of mechanism for MADS-box proteins in plant development.

In silico characterization and analysis of RTBP1 and NgTRF1 protein through MD simulation and molecular docking - A comparative study

Gaining access to sequence and structure information of telomere binding proteins helps in understanding the essential biological processes involve in conserved sequence specific interaction between DNA and the proteins. Rice telomere binding protein (RTBP1) and Nicotiana glutinosa telomere repeat binding factor (NgTRF1) are helix turn helix motif type of proteins that plays role in telomeric DNA protection and length regulation. Both the proteins share same type of domain but till now there is very less communication on the in silico studies of these complete proteins.Here we intend to do a comparative study between two proteins through modeling of the complete proteins, physiochemical characterization, MD simulation and DNA-protein docking. I-TASSER and CLC protein work bench was performed to find out the protein 3D structure as well as the different parameters to characterize the proteins. MD simulation was completed by GROMOS forcefield of GROMACS for 10 ns of time stretch. The simulated 3D structures were docked with template DNA (3D DNA modeled through 3D-DART) of TTTAGGG conserved sequence motif using HADDOCK web server.Digging up all the facts about the proteins it was reveled that around 120 amino acids in the tail part was showing a good sequence similarity between the proteins. Molecular modeling, sequence characterization and secondary structure prediction also indicates the similarity between the protein's structure and sequence. The result of MD simulation highlights on the RMSD, RMSF, Rg, PCA and Energy plots which also conveys the similar type of motional behavior between them. The best complex formation for both the proteins in docking result also indicates for the first interaction site which is mainly the helix3 region of the DNA binding domain. The overall computational analysis reveals that RTBP1 and NgTRF1 proteins display good amount of similarity in their physicochemical properties, structure, dynamics and binding mode.

Functional conservation and divergence of four ginger AP1/AGL9 MADS-box genes revealed by analysis of their expression and protein-protein interaction, and ectopic expression of AhFUL gene in Arabidopsis

Alpinia genus are known generally as ginger-lilies for showy flowers in the ginger family, Zingiberaceae, and their floral morphology diverges from typical monocotyledon flowers. However, little is known about the functions of ginger MADS-box genes in floral identity. In this study, four AP1/AGL9 MADS-box genes were cloned from Alpinia hainanensis, and protein-protein interactions (PPIs) and roles of the four genes in floral homeotic conversion and in floral evolution are surveyed for the first time. AhFUL is clustered to the AP1 lineage, AhSEP4 and AhSEP3b to the SEP lineage, and AhAGL6-like to the AGL6 lineage. The four genes showed conserved and divergent expression patterns, and their encoded proteins were localized in the nucleus. Seven combinations of PPI (AhFUL-AhSEP4, AhFUL-AhAGL6-like, AhFUL-AhSEP3b, AhSEP4-AhAGL6-like, AhSEP4-AhSEP3b, AhAGL6-like-AhSEP3b, and AhSEP3b-AhSEP3b) were detected, and the PPI patterns in the AP1/AGL9 lineage revealed that five of the 10 possible combinations are conserved and three are variable, while conclusions cannot yet be made regarding the other two. Ectopic expression of AhFUL in Arabidopsis thaliana led to early flowering and floral organ homeotic conversion to sepal-like or leaf-like. Therefore, we conclude that the four A. hainanensis AP1/AGL9 genes show functional conservation and divergence in the floral identity from other MADS-box genes.

ZmSOC1, a MADS-box transcription factor from Zea mays, promotes flowering in Arabidopsis

Zea mays is an economically important crop, but its molecular mechanism of flowering remains largely uncharacterized. The gene, SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1), integrates multiple flowering signals to regulate floral transition in Arabidopsis. In this study, ZmSOC1 was isolated from Zea mays. Sequence alignment and phylogenetic analysis demonstrated that the ZmSOC1 protein contained a highly conserved MADS domain and a typical SOC1 motif. ZmSOC1 protein was localized in the nucleus in protoplasts and showed no transcriptional activation activity in yeast cells. ZmSOC1 was highly expressed in maize reproductive organs, including filaments, ear and endosperm, but expression was very low in embryos; on the other hand, the abiotic stresses could repress ZmSOC1 expression. Overexpression of ZmSOC1 resulted in early flowering in Arabidopsis through increasing the expression of AtLFY and AtAP1. Overall, these results suggest that ZmSOC1 is a flowering promoter in Arabidopsis.

DEF- and GLO-like proteins may have lost most of their interaction partners during angiosperm evolution

BACKGROUND AND AIMS

DEFICIENS (DEF)- and GLOBOSA (GLO)-like proteins constitute two sister clades of floral homeotic transcription factors that were already present in the most recent common ancestor (MRCA) of extant angiosperms. Together they specify the identity of petals and stamens in flowering plants. In core eudicots, DEF- and GLO-like proteins are functional in the cell only as heterodimers with each other. There is evidence that this obligate heterodimerization contributed to the canalization of the flower structure of core eudicots during evolution. It remains unknown as to whether this strict heterodimerization is an ancient feature that can be traced back to the MRCA of extant flowering plants or if it evolved later during the evolution of the crown group angiosperms.

METHODS

The interactions of DEF- and GLO-like proteins of the early-diverging angiosperms Amborella trichopoda and Nuphar advena and of the magnoliid Liriodendron tulipifera were analysed by employing yeast two-hybrid analysis and electrophoretic mobility shift assay (EMSA). Character-state reconstruction, including data from other species as well, was used to infer the ancestral interaction patterns of DEF- and GLO-like proteins.

KEY RESULTS

The yeast two-hybrid and EMSA data suggest that DEF- and GLO-like proteins from early-diverging angiosperms both homo- and heterodimerize. Character-state reconstruction suggests that the ability to form heterodimeric complexes already existed in the MRCA of extant angiosperms and that this property remained highly conserved throughout angiosperm evolution. Homodimerization of DEF- and GLO-like proteins also existed in the MRCA of all extant angiosperms. DEF-like protein homodimerization was probably lost very early in angiosperm evolution and was not present in the MRCA of eudicots and monocots. GLO-like protein homodimerization might have been lost later during evolution, but very probably was not present in the MRCA of eudicots.

CONCLUSIONS

The flexibility of DEF- and GLO-like protein interactions in early-diverging angiosperms may be one reason for the highly diverse flower morphologies observed in these species. The results strengthen the hypothesis that a reduction in the number of interaction partners of DEF- and GLO-like proteins, with DEF-GLO heterodimers remaining the only DNA-binding dimers in core eudicots, contributed to developmental robustness, canalization of flower development and the diversification of angiosperms.

Arabidopsis SEPALLATA proteins differ in cooperative DNA-binding during the formation of floral quartet-like complexes

The SEPALLATA (SEP) genes of Arabidopsis thaliana encode MADS-domain transcription factors that specify the identity of all floral organs. The four Arabidopsis SEP genes function in a largely yet not completely redundant manner. Here, we analysed interactions of the SEP proteins with DNA. All of the proteins were capable of forming tetrameric quartet-like complexes on DNA fragments carrying two sequence elements termed CArG-boxes. Distances between the CArG-boxes for strong cooperative DNA-binding were in the range of 4-6 helical turns. However, SEP1 also bound strongly to CArG-box pairs separated by smaller or larger distances, whereas SEP2 preferred large and SEP4 preferred small inter-site distances for binding. Cooperative binding of SEP3 was comparatively weak for most of the inter-site distances tested. All SEP proteins constituted floral quartet-like complexes together with the floral homeotic proteins APETALA3 (AP3) and PISTILLATA (PI) on the target genes AP3 and SEP3. Our results suggest an important part of an explanation for why the different SEP proteins have largely, but not completely redundant functions in determining floral organ identity: they may bind to largely overlapping, but not identical sets of target genes that differ in the arrangement and spacing of the CArG-boxes in their cis-regulatory regions.

Tetramer formation in Arabidopsis MADS domain proteins: analysis of a protein-protein interaction network

BACKGROUND

MADS domain proteins are transcription factors that coordinate several important developmental processes in plants. These proteins interact with other MADS domain proteins to form dimers, and it has been proposed that they are able to associate as tetrameric complexes that regulate transcription of target genes. Whether the formation of functional tetramers is a widespread property of plant MADS domain proteins, or it is specific to few of these transcriptional regulators remains unclear.

RESULTS

We analyzed the structure of the network of physical interactions among MADS domain proteins in Arabidopsis thaliana. We determined the abundance of subgraphs that represent the connection pattern expected for a MADS domain protein heterotetramer. These subgraphs were significantly more abundant in the MADS domain protein interaction network than in randomized analogous networks. Importantly, these subgraphs are not significantly frequent in a protein interaction network of TCP plant transcription factors, when compared to expectation by chance. In addition, we found that MADS domain proteins in tetramer-like subgraphs are more likely to be expressed jointly than proteins in other subgraphs. This effect is mainly due to proteins in the monophyletic MIKC clade, as there is no association between tetramer-like subgraphs and co-expression for proteins outside this clade.

CONCLUSIONS

Our results support that the tendency to form functional tetramers is widespread in the MADS domain protein-protein interaction network. Our observations also suggest that this trend is prevalent, or perhaps exclusive, for proteins in the MIKC clade. Because it is possible to retrodict several experimental results from our analyses, our work can be an important aid to make new predictions and facilitates experimental research on plant MADS domain proteins.

Identification and characterization of FaSOC1, a homolog of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 from strawberry

A MADS-box gene SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) integrates multiple flowering signals to regulate floral transition in Arabidopsis. Strawberry (Fragaria spp.) is an economically important fruit crop, but its molecular control of flowering is largely unknown. In this study, a SOC1-like gene, FaSOC1, was isolated and characterized from strawberry. The open reading frame of FaSOC1 was 648bp, encoding a protein of 215 amino acids. Sequence alignment and phylogenetic analysis showed that the FaSOC1 protein contained a highly conserved MADS domain and a SOC1 motif, and that it was a member of the SOC1-like genes of dicots. The FaSOC1 protein mainly localized in the cytoplasm of onion epidermal cells and Arabidopsis protoplasts, and showed no transcriptional activation activity in yeast cells. Under the floral induction conditions, the expression of FaSOC1 increased during the first 2weeks of short-day treatment, but declined dramatically during three to 4weeks. FaSOC1 was highly expressed in reproductive organs, including shoot apices, floral buds, flowers, stamens and sepals. Overexpression of FaSOC1 in wild-type Arabidopsis caused early flowering and upregulated the expression of flowering time genes LFY and AP1. In addition, the yeast two-hybrid and BiFC assays confirmed that FaSOC1 could interact with AGL24. In conclusion, these results suggest that FaSOC1 is a flowering promoter in strawberry.

The oil palm SHELL gene controls oil yield and encodes a homologue of SEEDSTICK

A key event in the domestication and breeding of the oil palm, Elaeis guineensis, was loss of the thick coconut-like shell surrounding the kernel. Modern E. guineensis has three fruit forms, dura (thick-shelled), pisifera (shell-less) and tenera (thin-shelled), a hybrid between dura and pisifera^1–4. The pisifera palm is usually female-sterile but the tenera yields far more oil than dura, and is the basis for commercial palm oil production in all of Southeast Asia⁵. Here, we describe the mapping and identification of the Shell gene responsible for the different fruit forms. Using homozygosity mapping by sequencing we found two independent mutations in the DNA binding domain of a homologue of the MADS-box gene SEEDSTICK (STK) which controls ovule identity and seed development in Arabidopsis. The Shell gene is responsible for the tenera phenotype in both cultivated and wild palms from sub-Saharan Africa, and our findings provide a genetic explanation for the single gene heterosis attributed to Shell, via heterodimerization. This gene mutation explains the single most important economic trait in oil palm, and has implications for the competing interests of global edible oil production, biofuels and rainforest conservation⁶.

Tomato FRUITFULL homologues act in fruit ripening via forming MADS-box transcription factor complexes with RIN

The tomato MADS-box transcription factor RIN acts as a master regulator of fruit ripening. Here, we identified MADS-box proteins that interact with RIN; we also provide evidence that these proteins act in the regulation of fruit ripening. We conducted a yeast two-hybrid screen of a cDNA library from ripening fruit, for genes encoding proteins that bind to RIN. The screen identified two MADS-box genes, FUL1 and FUL2 (previously called TDR4 and SlMBP7), both of which have high sequence similarity to Arabidopsis FRUITFULL. Expression analyses revealed that the FUL1 mRNA and FUL1 protein accumulate in a ripening-specific manner in tomato fruits and FUL2 mRNA and protein accumulate at the pre-ripening stage and throughout ripening. Biochemical analyses confirmed that FUL1 and FUL2 form heterodimers with RIN; this interaction required the FUL1 and FUL2 C-terminal domains. Also, the heterodimers bind to a typical target DNA motif for MADS-box proteins. Chromatin immunoprecipitation assays revealed that FUL1 and FUL2 bind to genomic sites that were previously identified as RIN-target sites, such as the promoter regions of ACS2, ACS4 and RIN. These findings suggest that RIN forms complexes with FUL1 and FUL2 and these complexes regulate expression of ripening-related genes. In addition to the functional redundancy between FUL1 and FUL2, we also found they have potentially divergent roles in transcriptional regulation, including a difference in genomic target sites.

The interaction between OsMADS57 and OsTB1 modulates rice tillering via DWARF14

Rice tillering is a multigenic trait that influences grain yield, but its regulation molecular module is poorly understood. Here we report that OsMADS57 interacts with OsTB1 (TEOSINTE BRANCHED1) and targets D14 (Dwarf14) to control the outgrowth of axillary buds in rice. An activation-tagged mutant osmads57-1 and OsMADS57-overexpression lines showed increased tillers, whereas OsMADS57 antisense lines had fewer tillers. OsMIR444a-overexpressing lines exhibited suppressed OsMADS57 expression and tillering. Furthermore, osmads57-1 was insensitive to strigolactone treatment to inhibit axillary bud outgrowth, and OsMADS57’s function in tillering was dependent on D14. D14 expression was downregulated in osmads57-1, but upregulated in antisense and OsMIR444a-overexpressing lines. OsMADS57 bound to the CArG motif [C(A/T)TTAAAAAG] in the promoter and directly suppressed D14 expression. Interaction of OsMADS57 with OsTB1 reduced OsMADS57 inhibition of D14 transcription. Therefore, OsMIR444a-regulated OsMADS57, together with OsTB1, target D14 to control tillering. This regulation mechanism could have important application in rice molecular breeding programs focused on high grain yield.

Tillering is a multigenic complex trait that influences grain yield in cereal; however, the molecular network for its regulation remains unclear. Guo et al. show that OsMADS57, a transcription factor controlled by miR444a, interacts with OsTEOSINTE BRANCHED1 and targets DWARF14 to control tillering in rice.

Regulatory patterns of a large family of defensin-like genes expressed in nodules of Medicago truncatula

Root nodules are the symbiotic organ of legumes that house nitrogen-fixing bacteria. Many genes are specifically induced in nodules during the interactions between the host plant and symbiotic rhizobia. Information regarding the regulation of expression for most of these genes is lacking. One of the largest gene families expressed in the nodules of the model legume Medicago truncatula is the nodule cysteine-rich (NCR) group of defensin-like (DEFL) genes. We used a custom Affymetrix microarray to catalog the expression changes of 566 NCRs at different stages of nodule development. Additionally, bacterial mutants were used to understand the importance of the rhizobial partners in induction of NCRs. Expression of early NCRs was detected during the initial infection of rhizobia in nodules and expression continued as nodules became mature. Late NCRs were induced concomitantly with bacteroid development in the nodules. The induction of early and late NCRs was correlated with the number and morphology of rhizobia in the nodule. Conserved 41 to 50 bp motifs identified in the upstream 1,000 bp promoter regions of NCRs were required for promoter activity. These cis-element motifs were found to be unique to the NCR family among all annotated genes in the M. truncatula genome, although they contain sub-regions with clear similarity to known regulatory motifs involved in nodule-specific expression and temporal gene regulation.

Characterization of an AGAMOUS-like MADS box protein, a probable constituent of flowering and fruit ripening regulatory system in banana

The MADS-box family of genes has been shown to play a significant role in the development of reproductive organs, including dry and fleshy fruits. In this study, the molecular properties of an AGAMOUS like MADS box transcription factor in banana cultivar Giant governor (Musa sp, AAA group, subgroup Cavendish) has been elucidated. We have detected a CArG-box sequence binding AGAMOUS MADS-box protein in banana flower and fruit nuclear extracts in DNA-protein interaction assays. The protein fraction in the DNA-protein complex was analyzed by mass spectrometry and using this information we have obtained the full length cDNA of the corresponding protein. The deduced protein sequence showed ∼95% amino acid sequence homology with MA-MADS5, a MADS-box protein described previously from banana. We have characterized the domains of the identified AGAMOUS MADS-box protein involved in DNA binding and homodimer formation in vitro using full-length and truncated versions of affinity purified recombinant proteins. Furthermore, in order to gain insight about how DNA bending is achieved by this MADS-box factor, we performed circular permutation and phasing analysis using the wild type recombinant protein. The AGAMOUS MADS-box protein identified in this study has been found to predominantly accumulate in the climacteric fruit pulp and also in female flower ovary. In vivo and in vitro assays have revealed specific binding of the identified AGAMOUS MADS-box protein to CArG-box sequence in the promoters of major ripening genes in banana fruit. Overall, the expression patterns of this MADS-box protein in banana female flower ovary and during various phases of fruit ripening along with the interaction of the protein to the CArG-box sequence in the promoters of major ripening genes lead to interesting assumption about the possible involvement of this AGAMOUS MADS-box factor in banana fruit ripening and floral reproductive organ development.

Characterization of SOC1's central role in flowering by the identification of its upstream and downstream regulators

The transition from vegetative to reproductive development is one of the most important phase changes in the plant life cycle. This step is controlled by various environmental signals that are integrated at the molecular level by so-called floral integrators. One such floral integrator in Arabidopsis (Arabidopsis thaliana) is the MADS domain transcription factor SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1). Despite extensive genetic studies, little is known about the transcriptional control of SOC1, and we are just starting to explore the network of genes under the direct control of SOC1 transcription factor complexes. Here, we show that several MADS domain proteins, including SOC1 heterodimers, are able to bind SOC1 regulatory sequences. Genome-wide target gene analysis by ChIP-seq confirmed the binding of SOC1 to its own locus and shows that it also binds to a plethora of flowering-time regulatory and floral homeotic genes. In turn, the encoded floral homeotic MADS domain proteins appear to bind SOC1 regulatory sequences. Subsequent in planta analyses revealed SOC1 repression by several floral homeotic MADS domain proteins, and we show that, mechanistically, this depends on the presence of the SOC1 protein. Together, our data show that SOC1 constitutes a major hub in the regulatory networks underlying floral timing and flower development and that these networks are composed of many positive and negative autoregulatory and feedback loops. The latter seems to be crucial for the generation of a robust flower-inducing signal, followed shortly after by repression of the SOC1 floral integrator.

Identification and validation of sex-linked SCAR markers in dioecious Hippophae rhamnoides L. (Elaeagnaceae)

The actinorhizal plant seabuckthorn (Hippophae rhamnoides L., Elaeagnaceae) is a wind pollinated dioecious crop. To distinguish male genotypes from female genotypes early in the vegetative growth phase, we have developed robust PCR-based marker(s). DNA bulk samples from 20 male and 20 female plants each were screened with 60 RAPD primers. Two primers, OPA-04 and OPT-06 consistently amplified female-specific (FS) polymorphic fragments of 1,164 and 868 bp, respectively, that were absent in the male samples. DNA sequence of the two markers did not exhibit significant similarity to previously characterized sequences. A sequence-characterized amplified region marker HrX1 (JQ284019) and HrX2 (JQ284020) designed for the two fragments, continued to amplify the FS allele in 120 female plants but not in 100 male plants tested in the current study. Thus, HrX1 and HrX2 are FS markers that can determine the sex of seabuckthorn plants in an early stage and expedite cultivations for industrial applications.

A cotton gene encoding novel MADS-box protein is preferentially expressed in fibers and functions in cell elongation

Cotton fibers, as natural fibers, are widely used in the textile industry in the world. In order to find genes involved in fiber development, a cDNA (designated as GhMADS11) encoding a novel MADS protein with 151 amino acid residues was isolated from cotton fiber cDNA library. The deduced protein shares high similarity with Arabidopsis AP1 and AGL8 in MADS domain. However, the GhMADS11 protein (being absent of the partial K-domain and normal C-terminus) is shorter than AP1 and AGL8 by the reason of gene frameshift mutation during evolution. The experimental results revealed that GhMADS11 was not a transcriptional activator, and it did not form homodimer. GhMADS11 transcripts were specifically accumulated in elongating fibers, but no or very low signals of its expression were detected in other tissues of cotton. Overexpression of GhMADS11 in fission yeast promotes atypical cell elongation by 1.4-2.0-fold. Furthermore, morphological analysis indicated that the transformed cells expressing GhMADS11m, a MIKC-type derivative of GhMADS11 by the site-directed mutation, displayed the same phenotype as that of the transformed cells with GhMADS11. The concurrence of these data sets suggested that GhMADS11 protein may function in fiber cell elongation, and its MADS domain and partial K-domain are sufficient for this function.

Cell-type specific analysis of translating RNAs in developing flowers reveals new levels of control

Combining translating ribosome affinity purification with RNA-seq for cell-specific profiling of translating RNAs in developing flowers.

Cell type comparisons of cell type-specific hormone responses, promoter motifs, coexpressed cognate binding factor candidates, and splicing isoforms.

Widespread post-transcriptional regulation at both the intron splicing and translational stages.

A new class of noncoding RNAs associated with polysomes.

What constitutes a differentiated cell type? How much do cell types differ in their transcription of genes? The development and functions of tissues rely on constant interactions among distinct and nonequivalent cell types. Answering these questions will require quantitative information on transcriptomes, proteomes, protein–protein interactions, protein–nucleic acid interactions, and metabolomes at cellular resolution. The systems approaches emerging in biology promise to explain properties of biological systems based on genome-wide measurements of expression, interaction, regulation, and metabolism. To facilitate a systems approach, it is essential first to capture such components in a global manner, ideally at cellular resolution.

Recently, microarray analysis of transcriptomes has been extended to a cellular level of resolution by using laser microdissection or fluorescence-activated sorting (for review, see Nelson et al, 2008). These methods have been limited by stresses associated with cellular separation and isolation procedures, and biases associated with mandatory RNA amplification steps. A newly developed method, translating ribosome affinity purification (TRAP; Zanetti et al, 2005; Heiman et al, 2008; Mustroph et al, 2009), circumvents these problems by epitopetagging a ribosomal protein in specific cellular domains to selectively purify polysomes. We combined TRAP with deep sequencing, which we term TRAP-seq, to provide cell-level spatiotemporal maps for Arabidopsis early floral development at single-base resolution.

Flower development in Arabidopsis has been studied extensively and is one of the best understood aspects of plant development (for review, see Krizek and Fletcher, 2005). Genetic analysis of homeotic mutants established the ABC model, in which three classes of regulatory genes, A, B and C, work in a combinatorial manner to confer organ identities of four whorls (Coen and Meyerowitz, 1991). Each class of regulatory gene is expressed in a specific and evolutionarily conserved domain, and the action of the class A, B and C genes is necessary for specification of organ identity (Figure 1A).

Using TRAP-seq, we purified cell-specific translating mRNA populations, which we and others call the translatome, from the A, B and C domains of early developing flowers, in which floral patterning and the specification of floral organs is established. To achieve temporal specificity, we used a floral induction system to facilitate collection of early stage flowers (Wellmer et al, 2006). The combination of TRAP-seq with domain-specific promoters and this floral induction system enabled fine spatiotemporal isolation of translating mRNA in specific cellular domains, and at specific developmental stages.

Multiple lines of evidence confirmed the specificity of this approach, including detecting the expression in expected domains but not in other domains for well-studied flower marker genes and known physiological functions (Figures 1B–D and 2A–C). Furthermore, we provide numerous examples from flower development in which a spatiotemporal map of rigorously comparable cell-specific translatomes makes possible new views of the properties of cell domains not evident in data obtained from whole organs or tissues, including patterns of transcription and cis-regulation, new physiological differences among cell domains and between flower stages, putative hormone-active centers, and splicing events specific for flower domains (Figure 2A–D). Such findings may provide new targets for reverse genetics studies and may aid in the formulation and validation of interaction and pathway networks.

Beside cellular heterogeneity, the transcriptome is regulated at several steps through the life of mRNA molecules, which are not directly available through traditional transcriptome profiling of total mRNA abundance. By comparing the translatome and transcriptome, we integratively profiled two key posttranscriptional control points, intron splicing and translation state. From our translatome-wide profiling, we (i) confirmed that both posttranscriptional regulation control points were used by a large portion of the transcriptome; (ii) identified a number of cis-acting features within the coding or noncoding sequences that correlate with splicing or translation state; and (iii) revealed correlation between each regulation mechanism and gene function. Our transcriptome-wide surveys have highlighted target genes transcripts of which are probably under extensive posttranscriptional regulation during flower development.

Finally, we reported the finding of a large number of polysome-associated ncRNAs. About one-third of all annotated ncRNA in the Arabidopsis genome were observed co-purified with polysomes. Coding capacity analysis confirmed that most of them are real ncRNA without conserved ORFs. The group of polysome-associated ncRNA reported in this study is a potential new addition to the expanding riboregulator catalog; they could have roles in translational regulation during early flower development.

Determining both the expression levels of mRNA and the regulation of its translation is important in understanding specialized cell functions. In this study, we describe both the expression profiles of cells within spatiotemporal domains of the Arabidopsis thaliana flower and the post-transcriptional regulation of these mRNAs, at nucleotide resolution. We express a tagged ribosomal protein under the promoters of three master regulators of flower development. By precipitating tagged polysomes, we isolated cell type-specific mRNAs that are probably translating, and quantified those mRNAs through deep sequencing. Cell type comparisons identified known cell-specific transcripts and uncovered many new ones, from which we inferred cell type-specific hormone responses, promoter motifs and coexpressed cognate binding factor candidates, and splicing isoforms. By comparing translating mRNAs with steady-state overall transcripts, we found evidence for widespread post-transcriptional regulation at both the intron splicing and translational stages. Sequence analyses identified structural features associated with each step. Finally, we identified a new class of noncoding RNAs associated with polysomes. Findings from our profiling lead to new hypotheses in the understanding of flower development.

Molecular interactions of orthologues of floral homeotic proteins from the gymnosperm Gnetum gnemon provide a clue to the evolutionary origin of 'floral quartets'

Several lines of evidence suggest that the identity of floral organs in angiosperms is specified by multimeric transcription factor complexes composed of MADS-domain proteins. These bind to specific cis-regulatory elements ('CArG-boxes') of their target genes involving DNA-loop formation, thus constituting 'floral quartets'. Gymnosperms, angiosperms' closest relatives, contain orthologues of floral homeotic genes, but when and how the interactions constituting floral quartets were established during evolution has remained unknown. We have comprehensively studied the dimerization and DNA-binding of several classes of MADS-domain proteins from the gymnosperm Gnetum gnemon. Determination of protein-protein and protein-DNA interactions by yeast two-hybrid, in vitro pull-down and electrophoretic mobility shift assays revealed complex patterns of homo- and heterodimerization among orthologues of floral homeotic class B, class C and class E proteins and B(sister) proteins. Using DNase I footprint assays we demonstrate that both orthologues of class B with C proteins, and orthologues of class C proteins alone, but not orthologues of class B proteins alone can loop DNA in floral quartet-like complexes. This is in contrast to class B and class C proteins from angiosperms, which require other factors such as class E floral homeotic proteins to 'glue' them together in multimeric complexes. Our findings suggest that the evolutionary origin of floral quartet formation is based on the interaction of different DNA-bound homodimers, does not depend on class E proteins, and predates the origin of angiosperms.

Determination of sexual organ development

Plant sexual organ development is initiated from the floral meristem. At early stages, the activation of a set of genes that encode transcription factors determines the identity of the floral organs. These transcription factors are known as organ identity genes, and they form multimeric complexes that bind to target genes to control their expression. The transcriptional regulation of target genes triggers the formation of an organ by activating pathways required for its development initiating a cascade of events that leads to sexual plant reproduction. Here, I review the complex mechanisms involved in transcriptional regulation of organ identity genes and how they determine sexual organ development. Their expression is the result of complex interactions between repressors and activators that are often coexpressed. After the production of floral identity proteins, the formation of multimeric complexes defines target specificity and exerts a transcriptional regulatory effect on the target. Thanks to an increasing knowledge of the molecular control of sexual organ development in multiple species, we are beginning to understand how these genes evolved and how reproductive organ development occurs in different groups of plants. Comparative studies will, in future, provide a new insight into mechanisms of sexual organ development.

Systematic identification of X1-homologous genes reveals a family involved in stress responses in rice

X1-homologous genes (XHS) encode plant-specific proteins containing three major domains (XH, XS, zf-XS), but their functions are largely unknown. We report the systematic identification and characterization of XHS genes in the rice genome. Eleven putative XHS protein sequences (OXHS1-11) were identified in the sequenced genome of Oryza sativa japonica cv. Nipponbare, and these sequences, along with other plant XHS homologues, were classified into five subgroups based on phylogenetic analysis. Distinct diversification of the XHS proteins occurred between monocotyledon and dicotyledon plants. The OXHS family has diverse exon-intron structures and organizations of putative domains and motifs. The OXHS proteins showed no transactivation activity, and no interaction between the XH domain and the XS domain in yeast. Four representative OXHS proteins were targeted to cytoplasm, which contradicts the previous speculation that XHS proteins are putative transcription factors. All the OXHS genes are predominantly expressed in floral organs, and some are expressed in a wide range of tissues or organs in indica rice Minghui 63. Nine OXHS genes are responsive to at least one of the abiotic stresses including drought, salt, cold, and abscisic acid treatment. Over-expression of one stress-responsive gene OXHS2 in rice resulted in reduced tolerance to salt and drought stresses. These results suggest that the OXHS family may be functionally diversified and some members of this family may play important roles in regulating stress tolerance in rice.

Evolution of plant MADS box transcription factors: evidence for shifts in selection associated with early angiosperm diversification and concerted gene duplications

Phylogenomic analyses show that gene and genome duplication events have led to the diversification of transcription factor gene families throughout the evolutionary history of land plants and that gene duplications have played an important role in shaping regulatory networks influencing key phenotypic characters including floral development and flowering time. A molecular evolutionary investigation of the mode and tempo of selection acting on the angiosperm MADS box AP1/SQUA, AP3/PI, AG/AGL11, and SEP gene subfamilies revealed site-specific patterns of shifting evolutionary constraint throughout angiosperm history. Specific positions in the four canonical MADS box gene regions, especially K domains and C-terminal regions of all four of these MADS box gene subfamilies exhibited clade-specific shifts in selective constraint following concerted duplication events. Moreover, the frequency of site-specific shifts in constraint was correlated with gene duplications and early angiosperm diversification. We hypothesize that coevolution among interacting MADS box proteins may be responsible for simultaneous increases in the ratio of nonsynonymous to synonymous substitutions (d(N)/d(S) = omega) early in angiosperm history and following concerted duplication events.

Target genes of the MADS transcription factor SEPALLATA3: integration of developmental and hormonal pathways in the Arabidopsis flower

The molecular mechanisms by which floral homeotic genes act as major developmental switches to specify the identity of floral organs are still largely unknown. Floral homeotic genes encode transcription factors of the MADS-box family, which are supposed to assemble in a combinatorial fashion into organ-specific multimeric protein complexes. Major mediators of protein interactions are MADS-domain proteins of the SEPALLATA subfamily, which play a crucial role in the development of all types of floral organs. In order to characterize the roles of the SEPALLATA3 transcription factor complexes at the molecular level, we analyzed genome-wide the direct targets of SEPALLATA3. We used chromatin immunoprecipitation followed by ultrahigh-throughput sequencing or hybridization to whole-genome tiling arrays to obtain genome-wide DNA-binding patterns of SEPALLATA3. The results demonstrate that SEPALLATA3 binds to thousands of sites in the genome. Most potential target sites that were strongly bound in wild-type inflorescences are also bound in the floral homeotic agamous mutant, which displays only the perianth organs, sepals, and petals. Characterization of the target genes shows that SEPALLATA3 integrates and modulates different growth-related and hormonal pathways in a combinatorial fashion with other MADS-box proteins and possibly with non-MADS transcription factors. In particular, the results suggest multiple links between SEPALLATA3 and auxin signaling pathways. Our gene expression analyses link the genomic binding site data with the phenotype of plants expressing a dominant repressor version of SEPALLATA3, suggesting that it modulates auxin response to facilitate floral organ outgrowth and morphogenesis. Furthermore, the binding of the SEPALLATA3 protein to cis-regulatory elements of other MADS-box genes and expression analyses reveal that this protein is a key component in the regulatory transcriptional network underlying the formation of floral organs.

Most regulatory genes encode transcription factors, which modulate gene expression by binding to regulatory sequences of their target genes. In plants in particular, which genes are directly controlled by these transcription factors, and the molecular mechanisms of target gene recognition in vivo, are still largely unexplored. One of the best-understood developmental processes in plants is flower development. In different combinations, transcription factors of the MADS-box family control the identities of the different types of floral organs: sepals, petals, stamens, and carpels. Here, we present the first genome-wide analysis of binding sites of a MADS-box transcription factor in plants. We show that the MADS-domain protein SEPALLATA3 (SEP3) binds to the regulatory regions of thousands of potential target genes, many of which are also transcription factors. We provide insight into mechanisms of DNA recognition by SEP3, and suggest roles for other transcription factor families in SEP3 target gene regulation. In addition to effects on genes involved in floral organ identity, our data suggest that SEP3 binds to, and modulates, the transcription of target genes involved in hormonal signaling pathways.

The key floral regulator SEPALLATA3 binds to the promoters of a large number of potential direct target genes to integrate different growth-related and hormonal pathways in flower development.

Evolution of MIR168 paralogs in Brassicaceae

Background

In plants, expression of ARGONAUTE1 (AGO1), the catalytic subunit of the RNA-Induced Silencing Complex responsible for post-transcriptional gene silencing, is controlled through a feedback loop involving the miR168 microRNA. This complex auto-regulatory loop, composed of miR168-guided AGO1-catalyzed cleavage of AGO1 mRNA and AGO1-mediated stabilization of miR168, was shown to ensure the maintenance of AGO1 homeostasis that is pivotal for the correct functioning of the miRNA pathway.

Results

We applied different approaches to studying the genomic organization and the structural and functional evolution of MIR168 homologs in Brassicaeae. A whole genome comparison of Arabidopsis and poplar, phylogenetic footprinting and phylogenetic reconstruction were used to date the duplication events originating MIR168 homologs in these genomes. While orthology was lacking between Arabidopsis and poplar MIR168 genes, we successfully isolated orthologs of both loci present in Arabidopsis (MIR168a and MIR168b) from all the Brassicaceae species analyzed, including the basal species Aethionema grandiflora, thus indicating that (1) independent duplication events took place in Arabidopsis and poplar lineages and (2) the origin of MIR168 paralogs predates both the Brassicaceae radiation and the Arabidopsis alpha polyploidization. Different phylogenetic footprints, corresponding to known functionally relevant regions (transcription starting site and double-stranded structures responsible for microRNA biogenesis and function) or for which functions could be proposed, were found to be highly conserved among MIR168 homologs. Comparative predictions of the identified microRNAs also indicate extreme conservation of secondary structure and thermodynamic stability.

Conclusion

We used a comparative phylogenetic footprinting approach to identify the structural and functional constraints that shaped MIR168 evolution in Brassicaceae. Although their duplication happened at least 40 million years ago, we found evidence that both MIR168 paralogs have been maintained throughout the evolution of Brassicaceae, most likely functionally as indicated by the extremely high conservation of functionally relevant regions, predicted secondary structure and thermodynamic profile. Interestingly, the expression patterns observed in Arabidopsis indicate that MIR168b underwent partial subfunctionalization as determined by the experimental characterization of its expression pattern provided in this study. We found further evolutionary evidence that pre-miR168 lower stem (the RNA-duplex structure adjacent to the miR-miR* stem) is significantly longer than animal lower stems and probably plays a relevant role in multi-step miR168 biogenesis.

Reconstitution of 'floral quartets' in vitro involving class B and class E floral homeotic proteins

Homeotic MADS box genes encoding transcription factors specify the identity of floral organs by interacting in a combinatorial way. The ‘floral quartet model’, published several years ago, pulled together several lines of evidence suggesting that floral homeotic proteins bind as tetramers to two separated DNA sequence elements termed ‘CArG boxes’ by looping the intervening DNA. However, experimental support for ‘floral quartet’ formation remains scarce. Recently, we have shown that the class E floral homeotic protein SEPALLATA3 (SEP3) is sufficient to loop DNA in floral-quartet-like complexes in vitro. Here, we demonstrate that the class B floral homeotic proteins APETALA3 (AP3) and PISTILLATA (PI) do only weakly, at best, form floral-quartet-like structures on their own. However, they can be incorporated into such complexes together with SEP3. The subdomain K3 of SEP3 is of critical importance for the DNA-bound heterotetramers to be formed and is capable to mediate floral quartet formation even in the sequence context of AP3 and PI. Evidence is presented suggesting that complexes composed of SEP3, AP3 and PI form preferentially over other possible complexes. Based on these findings we propose a mechanism of how target gene specificity might be achieved at the level of floral quartet stability.

The class E floral homeotic protein SEPALLATA3 is sufficient to loop DNA in 'floral quartet'-like complexes in vitro

The organs of a eudicot flower are specified by four functional classes, termed class A, B, C and E, of MADS domain transcription factors. The combinatorial formation of tetrameric complexes, so called 'floral quartets', between these classes is widely believed to represent the molecular basis of floral organ identity specification. As constituents of all complexes, the class E floral homeotic proteins are thought to be of critical relevance for the formation of floral quartets. However, experimental support for tetrameric complex formation remains scarce. Here we provide physico-chemical evidence that in vitro homotetramers of the class E floral homeotic protein SEPALLATA3 from Arabidopsis thaliana bind cooperatively to two sequence elements termed 'CArG boxes' in a phase-dependent manner involving DNA looping. We further show that the N-terminal part of SEPALLATA3 lacking K3, a subdomain of the protein-protein interactions mediating K domain, and the C-terminal domain, is sufficient for protein dimerization, but not for tetramer formation and cooperative DNA binding. We hypothesize that the capacity of class E MADS domain proteins to form tetrameric complexes contributes significantly to the formation of floral quartets. Our findings further suggest that the spacing and phasing of CArG boxes are important parameters in the molecular mechanism by which floral homeotic proteins achieve target gene specificity.

DNA-binding specificity, transcriptional activation potential, and the rin mutation effect for the tomato fruit-ripening regulator RIN

The RIN gene encodes a putative MADS box transcription factor that controls tomato fruit ripening, and its ripening inhibitor (rin) mutation yields non-ripening fruit. In this study, the molecular properties of RIN and the rin mutant protein were clarified. The results revealed that the RIN protein accumulates in ripening fruit specifically and is localized in the nucleus of the cell. In vitro studies revealed that RIN forms a stable homodimer that binds to MADS domain-specific DNA sites. Analysis of binding site selection experiments revealed that the consensus binding sites of RIN highly resemble those of the SEPALLATA (SEP) proteins, which are Arabidopsis MADS box proteins that control the identity of floral organs. RIN exhibited a transcription-activating function similar to that exhibited by the SEP proteins. These results indicate that RIN exhibits similar molecular functions to SEP proteins although they play distinctly different biological roles. In vivo assays revealed that RIN binds to the cis-element of LeACS2. Our results also revealed that the rin mutant protein accumulates in the mutant fruit and exhibits a DNA-binding activity similar to that exhibited by the wild-type protein, but has lost its transcription-activating function, which in turn would inhibit ripening in mutant fruit.

Differential gene expression in Arabidopsis wild-type and mutant anthers: insights into anther cell differentiation and regulatory networks

In flowering plants, the anther contains highly specialized reproductive and somatic cells that are required for male fertility. Genetic studies have uncovered several genes that are important for anther development. However, little information is available regarding most genes active during anther development, including possible relationships between these genes and genetically defined regulators. In Arabidopsis, two previously isolated male-sterile mutants display dramatically altered anther cell differentiation patterns. The sporocyteless (spl)/nozzle (nzz) mutant is defective in the differentiation of primary sporogenous cells into microsporocytes, and does not properly form the anther wall. The excess microsporocytes1 (ems1)/extrasporogenous cells (exs) mutants produce excess microsporocytes at the expense of the tapetum. To gain additional insights into microsporocyte and tapetum differentiation and to uncover potential genetic interactions, expression profiles were compared between wild-type anthers (stage 4-6) and those of the spl or ems1 mutants. A total of 1954 genes were found to be differentially expressed in the ems1 and/or spl anthers, and these were grouped into 14 co-expression clusters. The presence of genes with known and predicted functions in specific clusters suggests potential functions for other genes in the same cluster. To obtain clues about possible co-regulation within co-expression clusters, we searched for shared cis-regulatory motifs in putative promoter regions. Our analyses were combined with data from previous studies to develop a model of the anther gene regulatory network. This model includes hypotheses that can be tested experimentally to gain further understanding of the mechanisms controlling anther development.

Identification of a sex-specific SCAR marker in dioecious Pandanus fascicularis L. (Pandanaceae)

Pandanus fascicularis L. is a dioecious plant native to South Asia with significant numbers in coastal areas. With the aim of distinguishing male genotypes from female genotypes early in the vegetative growth phase, the current study was initiated using molecular markers. Based on the principle of bulked segregant analysis, the sampled plants were separated into 2 bulks depending on their sex. Of the 89 random amplified polymorphic DNA and inter-simple sequence repeat markers used, one decamer (OPO-08) consistently amplified a 1263 bp band in the males that was absent in the females. Its DNA sequence did not exhibit significant similarity to previously characterized sequences, but the presence of mononucleotide and dinucleotide repeats suggested that it was a repeat-rich region. A sequence-characterized amplified region marker (MSSRF-01) designed for this fragment continued to amplify the specific allele in all the male plants. Southern hybridization performed using the sex-specific fragment as a probe yielded results consistent with those previously obtained by polymerase chain reaction. These results strongly suggest that MSSRF-01 is a male-specific molecular marker. With no information available on the presence of sex chromosomes in Pandanus , this marker can be used to differentiate the sexes.

AthaMap web tools for the analysis and identification of co-regulated genes

The AthaMap database generates a map of cis-regulatory elements for the whole Arabidopsis thaliana genome. This database has been extended by new tools to identify common cis-regulatory elements in specific regions of user-provided gene sets. A resulting table displays all cis-regulatory elements annotated in AthaMap including positional information relative to the respective gene. Further tables show overviews with the number of individual transcription factor binding sites (TFBS) present and TFBS common to the whole set of genes. Over represented cis-elements are easily identified. These features were used to detect specific enrichment of drought-responsive elements in cold-induced genes. For identification of co-regulated genes, the output table of the colocalization function was extended to show the closest genes and their relative distances to the colocalizing TFBS. Gene sets determined by this function can be used for a co-regulation analysis in microarray gene expression databases such as Genevestigator or PathoPlant. Additional improvements of AthaMap include display of the gene structure in the sequence window and a significant data increase. AthaMap is freely available at .

trans meets cis in MADS science

The interaction between a transcription factor and its binding site at the DNA is an integral part of transcriptional regulatory networks, which is fundamental for an understanding of biological processes. An example is the family of MADS domain transcription factors, which represent key regulators of processes in yeast, animals and plants. However, despite our extensive knowledge of these transcription factors, limited information is available on the cis-elements to which these proteins bind or how these elements are defined. Here, we discuss the current understanding of MADS protein binding sites and compare data from various organisms. This information can help us in developing algorithms to predict binding sites for MADS domain transcription factors, which would be a significant step forward in the identification of "down-stream" target genes and the elucidation of transcriptional networks.

Mutant analysis, protein-protein interactions and subcellular localization of the Arabidopsis B sister (ABS) protein

Recently, close relatives of class B floral homeotic genes, termed B(sister) genes, have been identified in both angiosperms and gymnosperms. In contrast to the B genes themselves, B(sister) genes are exclusively expressed in female reproductive organs, especially in the envelopes or integuments surrounding the ovules. This suggests an important ancient function in ovule or seed development for B(sister) genes, which has been conserved for about 300 million years. However, investigation of the first loss-of-function mutant for a B(sister) gene (ABS/TT16 from Arabidopsis) revealed only a weak phenotype affecting endothelium formation. Here, we present an analysis of two additional mutant alleles, which corroborates this weak phenotype. Transgenic plants that ectopically express ABS show changes in the growth and identity of floral organs, suggesting that ABS can interact with floral homeotic proteins. Yeast-two-hybrid and three-hybrid analyses indicated that ABS can form dimers with SEPALLATA (SEP) floral homeotic proteins and multimeric complexes that also include the AGAMOUS-like proteins SEEDSTICK (STK) or SHATTERPROOF1/2 (SHP1, SHP2). These data suggest that the formation of multimeric transcription factor complexes might be a general phenomenon among MIKC-type MADS-domain proteins in angiosperms. Heterodimerization of ABS with SEP3 was confirmed by gel retardation assays. Fusion proteins tagged with CFP (Cyan Fluorescent Protein) and YFP (Yellow Fluorescent Protein) in Arabidopsis protoplasts showed that ABS is localized in the nucleus. Phylogenetic analysis revealed the presence of a structurally deviant, but closely related, paralogue of ABS in the Arabidopsis genome. Thus the evolutionary developmental genetics of B(sister) genes can probably only be understood as part of a complex and redundant gene network that may govern ovule formation in a conserved manner, which has yet to be fully explored.

Molecular cloning and characterization of a CBF gene from Capsella bursa-pastoris

A new CBF gene was cloned from Capsella bursa-pastoris(shepherd's purse) by rapid amplification of cDNA ends (RACE). The full-length cDNA of C. bursa-pastoris CBF gene (designated as Cbcbf) was 1034 bp long and contained a 657 bp open reading frame (ORF) encoding a putative DRE/CRT (LTRE)-binding protein of 219 amino acids. The predicted CbCBF protein was found to have a potential nuclear localization signal (NLS) in its N-terminal region followed by an AP2 DNA-binding motif and an acidic C-terminal half that might act as an activator domain. Bioinformatic analysis revealed that Cbcbf strongly resembled other CBF genes from Arabidopsis thaliana (cbf1, cbf2, cbf3) and Brassica napus (Bncbf5, Bncbf 7, Bncbf16 and Bncbf17). Subsequent cold acclimation assay showed that Cbcbf was relevant to cold acclimation. Our study implies that Cbcbf might have similar functions possessed by other CBF genes such as inducing the expression of some cold-regulated genes and increasing plants' freezing tolerance.

MIKC-type MADS-domain proteins: structural modularity, protein interactions and network evolution in land plants

MIKC-type proteins represent a class of MADS-domain transcription factors and are defined by a unique domain structure: in addition to the highly conserved DNA-binding MADS-domain, they have three other domains ('I', 'K' and 'C'), with the keratin-like K-domain being the most highly conserved and characteristic one. The number and functional diversity of MIKC-type proteins increased considerably during land plant evolution, culminating in higher flowering plants, where they dominate the control of reproductive development from early to late stages. We wonder how one special class of proteins became important in the control of essentially all stages of a morphogenetic process. All MADS-domain proteins appear to bind to DNA as homo- or heterodimers and may function as part of ternary transcription factor complexes involving non-MADS-domain proteins. Only MIKC-type proteins, however, generate complex intrafamily interaction networks. These are based on the special potential of MIKC-type proteins to form complexes involving more than two homologous proteins constituting transcriptional regulators. We speculate that the potential to form heteromultimers of homologous proteins was achieved by the acquisition of the K-domain during evolution. There is emerging evidence that organismal complexity arises from progressively more elaborate regulation of gene expression. We hypothesize that combinatorial multimer formation of MIKC-type MADS-domain proteins facilitated an unusually efficient and rapid functional diversification based on gene duplication, sequence divergence and fixation. This 'networking' may have enabled a more sophisticated transcriptional control of target genes which was recruited for controlling increasingly complex and diverse developmental pathways during the rapid origin and diversification of plant reproductive structures. Therefore, MIKC-type proteins may owe their evolutionary 'success' and present-day developmental importance in part to their modular domain structure. Investigating the evolution of MIKC-type genes may thus help to better understand origin and diversification of gene regulatory networks.