Adaptive evolution in the Arabidopsis MADS-box gene family inferred from its complete resolved phylogeny

Structural and Functional Analysis of the MADS-Box Genes Reveals Their Functions in Cold Stress Responses and Flower Development in Tea Plant (Camellia sinensis)

MADS-box genes comprise a large family of transcription factors that play crucial roles in all aspects of plant growth and development. However, no detailed information on the evolutionary relationship and functional characterization of MADS-box genes is currently available for some representative lineages, such as the Camellia plant. In this study, 136 MADS-box genes were detected from a reference genome of the tea plant (Camellia sinensis) by employing a 569 bp HMM (Hidden Markov Model) developed using nucleotide sequencing including 73 type I and 63 type II genes. An additional twenty-seven genes were identified, with five MIKC-type genes. Truncated and/or inaccurate gene models were manually verified and curated to improve their functional characterization. Subsequently, phylogenetic relationships, chromosome locations, conserved motifs, gene structures, and gene expression profiles were systematically investigated. Tea plant MIKC genes were divided into all 14 major eudicot subfamilies, and no gene was found in Mβ. The expansion of MADS-box genes in the tea plant was mainly contributed by WGD/fragment and tandem duplications. The expression profiles of tea plant MADS-box genes in different tissues and seasons were analyzed, revealing widespread evolutionary conservation and genetic redundancy. The expression profiles linked to cold stress treatments suggested the wide involvement of MADS-box genes from the tea plant in response to low temperatures. Moreover, a floral 'ABCE' model was proposed in the tea plant and proved to be both conserved and ancient. Our analyses offer a detailed overview of MADS-box genes in the tea plant, allowing us to hypothesize the potential functions of unknown genes and providing a foundation for further functional characterizations.

Genome-Wide Identification of MADS-Box Genes in Taraxacum kok-saghyz and Taraxacum mongolicum: Evolutionary Mechanisms, Conserved Functions and New Functions Related to Natural Rubber Yield Formation

MADS-box transcription regulators play important roles in plant growth and development. However, very few MADS-box genes have been isolated in the genus Taraxacum, which consists of more than 3000 species. To explore their functions in the promising natural rubber (NR)-producing plant Taraxacum kok-saghyz (TKS), MADS-box genes were identified in the genome of TKS and the related species Taraxacum mongolicum (TM; non-NR-producing) via genome-wide screening. In total, 66 TkMADSs and 59 TmMADSs were identified in the TKS and TM genomes, respectively. From diploid TKS to triploid TM, the total number of MADS-box genes did not increase, but expansion occurred in specific subfamilies. Between the two genomes, a total of 11 duplications, which promoted the expansion of MADS-box genes, were identified in the two species. TkMADS and TmMADS were highly conserved, and showed good collinearity. Furthermore, most TkMADS genes exhibiting tissue-specific expression patterns, especially genes associated with the ABCDE model, were preferentially expressed in the flowers, suggesting their conserved and dominant functions in flower development in TKS. Moreover, by comparing the transcriptomes of different TKS lines, we identified 25 TkMADSs related to biomass formation and 4 TkMADSs related to NR content, which represented new targets for improving the NR yield of TKS.

Comparative transcriptome analysis on the mangrove Acanthus ilicifolius and its two terrestrial relatives provides insights into adaptation to intertidal habitats

Acanthus is a unique genus covering both mangroves and terrestrial species, and thus is an ideal system to comparatively analyze the mechanisms of mangrove adaptation to intertidal habitats. We performed RNA sequencing of the mangrove plant Acanthus ilicifolius and its two terrestrial relatives, Acanthus leucostachyus and Acanthus mollis. A total of 91,125, 118,290, and 141,640 unigenes were obtained. Simple sequence repeats (SSR) analysis showed that A. ilicifolius had more SSRs, the highest frequency of distribution, and higher in polymorphism potential compared to the two terrestrial relatives. Phylogenetic analyses suggested a relatively recent split between A. ilicifolius and A. leucostachyus, i.e., about 16.76 million years ago (Mya), after their ancestor divergence with A. mollis (32.11 Mya), indicating that speciation of three Acanthus species occurred in the Early to Middle Miocene. Gene Ontology (GO) enrichment revealed that the unique unigenes in A. ilicifolius are predominantly related to rhythmic process, reproductive process and response to stimuli. The accelerated evolution and positive selection analyses indicated that the genus Acanthus migrated from terrestrial to intertidal habitats, where 311 pairs may be under positive selection. Functional enrichment analysis revealed that these genes associated with essential metabolism and biosynthetic pathways such as oxidative phosphorylation, plant hormone signal transduction, photosynthetic carbon fixation and arginine and proline metabolism, are related to the adaptation of A. ilicifolius to intertidal habitats, which are characterized by high salinity and hypoxia. Our results indicate the evolutionary processes and the mechanisms underlying the adaptability of Acanthus to various harsh environments from the arid terrestrial to intertidal habitats.

The MADS transcription factor GhFYF is involved in abiotic stress responses in upland cotton (Gossypium hirsutum L.)

Cotton is an important textile industry raw material crops, which plays a critical role in the development of society. MADS transcription factors (TFs) play a key role about the flowering time, flower development, and abiotic stress responses in plants, but little is known about their functions on abiotic stress in cotton. In this study, a MIKC^C subfamily gene from cotton, GhFYF (FOREVER YOUNG FLOWER), was isolated and characterized. Our data showed that GhFYF localized to the nucleus. A β-glucuronidase (GUS) activity assay revealed that the promoter of GhFYF was mainly expressed in the flower and seed of ProGhFYF::GUS transgenic A. thaliana plants. The GUS staining of flowers and seeds was deepened after drought, salt treatment, and the expression level of the GUS gene and corresponding stress genes AtERD10, AtAnnexin1 are up-regulated in the inflorescence. Overexpression GhFYF in A. thaliana could promote the seed germination and growth under different salt concentrations, and determin the proline content. Yeast two-hybrid (Y2H) assays showed that GhFYF interacted with the HAD-like protein GhGPP2, which has responds to abiotic stress. Our findings indicate that GhFYF is involved in abiotic stress responses, especially for salt stress. This work establishes a solid foundation for further functional analysis of the GhFYF gene in cotton.

Functional Conservation and Divergence of Five AP1/FUL-like Genes in Marigold (Tagetes erecta L.)

Members of AP1/FUL subfamily genes play an essential role in the regulation of floral meristem transition, floral organ identity, and fruit ripping. At present, there have been insufficient studies to explain the function of the AP1/FUL-like subfamily genes in Asteraceae. Here, we cloned two euAP1 clade genes TeAP1-1 and TeAP1-2, and three euFUL clade genes TeFUL1, TeFUL2, and TeFUL3 from marigold (Tagetes erecta L.). Expression profile analysis demonstrated that TeAP1-1 and TeAP1-2 were mainly expressed in receptacles, sepals, petals, and ovules. TeFUL1 and TeFUL3 were expressed in flower buds, stems, and leaves, as well as reproductive tissues, while TeFUL2 was mainly expressed in flower buds and vegetative tissues. Overexpression of TeAP1-2 or TeFUL2 in Arabidopsis resulted in early flowering, implying that these two genes might regulate the floral transition. Yeast two-hybrid analysis indicated that TeAP1/FUL proteins only interacted with TeSEP proteins to form heterodimers and that TeFUL2 could also form a homodimer. In general, TeAP1-1 and TeAP1-2 might play a conserved role in regulating sepal and petal identity, similar to the functions of MADS-box class A genes, while TeFUL genes might display divergent functions. This study provides a theoretical basis for the study of AP1/FUL-like genes in Asteraceae species.

Chromosome-Scale Reference Genome of Amphicarpaea edgeworthii: A New Resource for Amphicarpic Plants Research and Complex Flowering Pattern

Amphicarpaea edgeworthii, an annual twining herb, is a widely distributed species and an attractive model for studying complex flowering types and evolutionary mechanisms of species. Herein, we have generated a high-quality assembly of A. edgeworthii by using a combination of PacBio, 10× Genomics libraries, and Hi-C mapping technologies. The final 11 chromosome-level scaffolds covered 90.61% of the estimated genome (343.78Mb), which is a chromosome-scale assembled genome of an amphicarpic plant. Subsequently, we characterized the genetic diversity and population structure of A. edgeworthii species by resequencing individuals collected from their natural area of distribution. Using transcriptome profiling, we observed that specific phenotypes are regulated by a complex network of light, hormones, and MADS-box gene families. These data are beneficial for the discovery of genes that control major agronomic traits and spur genetic improvement of and functional genetic studies in legumes, as well as supply comparative genetic resources for other amphicarpic plants.

A high-quality genome assembly of Jasminum sambac provides insight into floral trait formation and Oleaceae genome evolution

As one of the most economically significant Oleaceae family members, Jasminum sambac is renowned for its distinct sweet, heady fragrance. Using Illumina reads, Nanopore long reads, and HiC-sequencing, we efficiently assembled and annotated the J. sambac genome. The high-quality genome assembly consisted of a total of 507 Mb sequence (contig N50 = 17.6 Mb) with 13 pseudomolecules. A total of 21,143 protein-coding genes and 303 Mb repeat sequences were predicted. An ancient whole-genome triplication event at the base of Oleaceae (~66 million years ago [Ma], Late Cretaceous) was identified and this may have contributed to the diversification of the Oleaceae ancestor and its divergence from the Lamiales. Stress-related (e.g., WRKY) and flowering-related (e.g., MADS-box) genes were located in the triplicated regions, suggesting that the polyploidy event might have contributed adaptive potential. Genes related to terpenoid biosynthesis, for example, FTA and TPS, were observed to be duplicated to a great extent in the J. sambac genome, perhaps explaining the strong fragrance of the flowers. Copy number changes in distinct phylogenetic clades of the MADS-box family were observed in J. sambac genome, for example, AGL6- and Mα- were lost and SOC- expanded, features that might underlie the long flowering period of J. sambac. The structural genes implicated in anthocyanin biosynthesis were depleted and this may explain the absence of vivid colours in jasmine. Collectively, assembling the J. sambac genome provides new insights into the genome evolution of the Oleaceae family and provides mechanistic insights into floral properties.

The intervening domain is required for DNA-binding and functional identity of plant MADS transcription factors

The MADS transcription factors (TF) are an ancient eukaryotic protein family. In plants, the family is divided into two main lineages. Here, we demonstrate that DNA binding in both lineages absolutely requires a short amino acid sequence C-terminal to the MADS domain (M domain) called the Intervening domain (I domain) that was previously defined only in type II lineage MADS. Structural elucidation of the MI domains from the floral regulator, SEPALLATA3 (SEP3), shows a conserved fold with the I domain acting to stabilise the M domain. Using the floral organ identity MADS TFs, SEP3, APETALA1 (AP1) and AGAMOUS (AG), domain swapping demonstrate that the I domain alters genome-wide DNA-binding specificity and dimerisation specificity. Introducing AG carrying the I domain of AP1 in the Arabidopsis ap1 mutant resulted in strong complementation and restoration of first and second whorl organs. Taken together, these data demonstrate that the I domain acts as an integral part of the DNA-binding domain and significantly contributes to the functional identity of the MADS TF.

MADS transcription factors regulate multiple aspects of plant development. Here the authors show that the intervening I domain is conserved in both type I and type II plant MADS lineages and contributes to the functional identity of the protein by influencing both DNA binding activity and dimerisation specificity.

Genome-wide (ChIP-seq) identification of target genes regulated by WRKY33 during submergence stress in Arabidopsis

Background

Hypoxia induced by flooding causes significant losses to crop production almost every year. However, the molecular network of submergence signaling pathway is still poorly understood. According to previous studies, transgenic plants overexpressing the WRKY33 gene showed enhanced resistance to submergence stress. Thus, this transcription factor may regulate a series of target genes in response to submergence. Here, to determine putative downstream targets of WRKY33 at a genome-wide scale in Arabidopsis thaliana, we performed the chromatin immunoprecipitation sequencing (ChIP-seq) using 35S:FLAG-WRKY33 overexpression transgenic lines (WRKY33-OE) after 24 h of submergence treatment.

Results

Using ChIP-seq data, we identified a total of 104 WRKY33-binding genes under submergence stress (WRKY33BGSs). Most WRKY33BGSs are involved in the oxidation-reduction process, programmed cell death in response to reactive oxygen species, lipid biosynthesis process, and other processes related to stress responses. Moreover, the major motif identified in the WRKY33BGSs promoters is a new cis-element, TCTCTC (named here as “TC box”). This cis-element differs from the previously known W box for WRKY33. Further qPCR experiments verified that genes carrying this motif in their promoters could be regulated by WRKY33 upon submergence treatment.

Conclusions

Our study has identified a new putative binding motif of WRKY33 and recovered numerous previously unknown target genes of WRKY33 during submergence stress. The WRKY33 gene positively participates in flooding response probably by transcriptional regulation of the downstream submergence-related target genes via a “TC box”.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00972-5.

MADS-box genes underground becoming mainstream: plant root developmental mechanisms

Plant growth is largely post-embryonic and depends on meristems that are active throughout the lifespan of an individual. Developmental patterns rely on the coordinated spatio-temporal expression of different genes, and the activity of transcription factors is particularly important during most morphogenetic processes. MADS-box genes constitute a transcription factor family in eukaryotes. In Arabidopsis, their proteins participate in all major aspects of shoot development, but their role in root development is still not well characterized. In this review we synthetize current knowledge pertaining to the function of MADS-box genes highly expressed in roots: XAL1, XAL2, ANR1 and AGL21, as well as available data for other MADS-box genes expressed in this organ. The role of Trithorax group and Polycomb group complexes on MADS-box genes' epigenetic regulation is also discussed. We argue that understanding the role of MADS-box genes in root development of species with contrasting architectures is still a challenge. Finally, we propose that MADS-box genes are key components of the gene regulatory networks that underlie various gene expression patterns, each one associated with the distinct developmental fates observed in the root. In the case of XAL1 and XAL2, their role within these networks could be mediated by regulatory feedbacks with auxin.

Evolution and Diversification of FRUITFULL Genes in Solanaceae

Ecologically and economically important fleshy edible fruits have evolved from dry fruit numerous times during angiosperm diversification. However, the molecular mechanisms that underlie these shifts are unknown. In the Solanaceae there has been a major shift to fleshy fruits in the subfamily Solanoideae. Evidence suggests that an ortholog of FRUITFULL (FUL), a transcription factor that regulates cell proliferation and limits the dehiscence zone in the silique of Arabidopsis, plays a similar role in dry-fruited Solanaceae. However, studies have shown that FUL orthologs have taken on new functions in fleshy fruit development, including regulating elements of tomato ripening such as pigment accumulation. FUL belongs to the core eudicot euFUL clade of the angiosperm AP1/FUL gene lineage. The euFUL genes fall into two paralogous clades, euFULI and euFULII. While most core eudicots have one gene in each clade, Solanaceae have two: FUL1 and FUL2 in the former, and MBP10 and MBP20 in the latter. We characterized the evolution of the euFUL genes to identify changes that might be correlated with the origin of fleshy fruit in Solanaceae. Our analyses revealed that the Solanaceae FUL1 and FUL2 clades probably originated through an early whole genome multiplication event. By contrast, the data suggest that the MBP10 and MBP20 clades are the result of a later tandem duplication event. MBP10 is expressed at weak to moderate levels, and its atypical short first intron lacks putative transcription factor binding sites, indicating possible pseudogenization. Consistent with this, our analyses show that MBP10 is evolving at a faster rate compared to MBP20. Our analyses found that Solanaceae euFUL gene duplications, evolutionary rates, and changes in protein residues and expression patterns are not correlated with the shift in fruit type. This suggests deeper analyses are needed to identify the mechanism underlying the change in FUL ortholog function.

Genome-wide survey of potato MADS-box genes reveals that StMADS1 and StMADS13 are putative downstream targets of tuberigen StSP6A

Background

MADS-box genes encode transcription factors that are known to be involved in several aspects of plant growth and development, especially in floral organ specification. To date, the comprehensive analysis of potato MADS-box gene family is still lacking after the completion of potato genome sequencing. A genome-wide characterization, classification, and expression analysis of MADS-box transcription factor gene family was performed in this study.

Results

A total of 153 MADS-box genes were identified and categorized into MIKC subfamily (MIKC^C and MIKC^*) and M-type subfamily (Mα, Mβ, and Mγ) based on their phylogenetic relationships to the Arabidopsis and rice MADS-box genes. The potato M-type subfamily had 114 members, which is almost three times of the MIKC members (39), indicating that M-type MADS-box genes have a higher duplication rate and/or a lower loss rate during potato genome evolution. Potato MADS-box genes were present on all 12 potato chromosomes with substantial clustering that mainly contributed by the M-type members. Chromosomal localization of potato MADS-box genes revealed that MADS-box genes, mostly MIKC, were located on the duplicated segments of the potato genome whereas tandem duplications mainly contributed to the M-type gene expansion. The potato MIKC subfamily could be further classified into 11 subgroups and the TT16-like, AGL17-like, and FLC-like subgroups found in Arabidopsis were absent in potato. Moreover, the expressions of potato MADS-box genes in various tissues were analyzed by using RNA-seq data and verified by quantitative real-time PCR, revealing that the MIKC^C genes were mainly expressed in flower organs and several of them were highly expressed in stolon and tubers. StMADS1 and StMADS13 were up-regulated in the StSP6A-overexpression plants and down-regulated in the StSP6A-RNAi plant, and their expression in leaves and/or young tubers were associated with high level expression of StSP6A.

Conclusion

Our study identifies the family members of potato MADS-box genes and investigate the evolution history and functional divergence of MADS-box gene family. Moreover, we analyze the MIKC^C expression patterns and screen for genes involved in tuberization. Finally, the StMADS1 and StMADS13 are most likely to be downstream target of StSP6A and involved in tuber development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5113-z) contains supplementary material, which is available to authorized users.

Genome-wide (ChIP-seq) identification of target genes regulated by BdbZIP10 during paraquat-induced oxidative stress

BACKGROUND

bZIP transcription factors play a significant role in many aspects of plant growth and development and also play critical regulatory roles during plant responses to various stresses. Overexpression of the Brachypodium bZIP10 (Bradi1g30140) transcription factor conferred enhanced oxidative stress tolerance and increased viability when plants or cells were exposed to the herbicide paraquat. To gain a better understanding of genes involved in bZIP10 conferred oxidative stress tolerance, chromatin immunoprecipitation followed by high throughput sequencing (ChIP-Seq) was performed on BdbZIP10 overexpressing plants in the presence of oxidative stress.

RESULTS

We identified a transcription factor binding motif, TGDCGACA, different from most known bZIP TF motifs but with strong homology to the Arabidopsis zinc deficiency response element. Analysis of the immunoprecipitated sequences revealed an enrichment of gene ontology groups with metal ion transmembrane transporter, transferase, catalytic and binding activities. Functional categories including kinases and phosphotransferases, cation/ion transmembrane transporters, transferases (phosphorus-containing and glycosyl groups), and some nucleoside/nucleotide binding activities were also enriched.

CONCLUSIONS

Brachypodium bZIP10 is involved in zinc homeostasis, as it relates to oxidative stress.

Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation

Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.

Protein interaction evolution from promiscuity to specificity with reduced flexibility in an increasingly complex network

A key question regarding protein evolution is how proteins adapt to the dynamic environment in which they function and how in turn their evolution shapes the protein interaction network. We used extant and resurrected ancestral plant MADS-domain transcription factors to understand how SEPALLATA3, a protein with hub and glue properties, evolved and takes part in network organization. Although the density of dimeric interactions was saturated in the network, many new interactions became mediated by SEPALLATA3 after a whole genome triplication event. By swapping SEPALLATA3 and its ancestors between dimeric networks of different ages, we found that the protein lost the capacity of promiscuous interaction and acquired specificity in evolution. This was accompanied with constraints on conformations through proline residue accumulation, which made the protein less flexible. SHORT VEGETATIVE PHASE on the other hand (non-hub) was able to gain protein-protein interactions due to a C-terminal domain insertion, allowing for a larger interaction interface. These findings illustrate that protein interaction evolution occurs at the level of conformational dynamics, when the binding mechanism concerns an induced fit or conformational selection. Proteins can evolve towards increased specificity with reduced flexibility when the complexity of the protein interaction network requires specificity.

Natural haplotypes of FLM non-coding sequences fine-tune flowering time in ambient spring temperatures in Arabidopsis

Cool ambient temperatures are major cues determining flowering time in spring. The mechanisms promoting or delaying flowering in response to ambient temperature changes are only beginning to be understood. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) regulates flowering in the ambient temperature range and FLM is transcribed and alternatively spliced in a temperature-dependent manner. We identify polymorphic promoter and intronic sequences required for FLM expression and splicing. In transgenic experiments covering 69% of the available sequence variation in two distinct sites, we show that variation in the abundance of the FLM-ß splice form strictly correlate (R² = 0.94) with flowering time over an extended vegetative period. The FLM polymorphisms lead to changes in FLM expression (PRO2+) but may also affect FLM intron 1 splicing (INT6+). This information could serve to buffer the anticipated negative effects on agricultural systems and flowering that may occur during climate change.

DOI:http://dx.doi.org/10.7554/eLife.22114.001

Analysis of the floral MADS-box genes from monocotyledonous Trilliaceae species indicates the involvement of SEPALLATA3-like genes in sepal-petal differentiation

The evolution of greenish sepals from petaloid outer tepals has occurred repeatedly in various lineages of non-grass monocots. Studies in distinct monocot species showed that the evolution of sepals could be explained by the ABC model; for example, the defect of B-class function in the outermost whorl was linked to the evolution of sepals. Here, floral MADS-box genes from three sepal-bearing monocotyledonous Trilliaceae species, Trillium camschatcense, Paris verticillata, and Kinugasa japonica were examined. Unexpectedly, expression of not only A- but also B-class genes was detected in the sepals of all three species. Although the E-class gene is generally expressed across all floral whorls, no expression was detected in sepals in the three species examined here. Overexpression of the E-class SEPALLATA3-like gene from T. camschatcense (TcamSEP) in Arabidopsis thaliana produced phenotypes identical to those reported for orthologs in other monocots. Additionally, yeast hybrid experiments indicated that TcamSEP could form a higher-order complex with an endogenous heterodimer of B-class APETALA3/DEFICIENS-like (TcamDEF) and PISTILLATA/GLOBOSA-like (TcamGLO) proteins. These results suggest a conserved role for Trilliaceae SEPALLATA3-like genes in functionalization of the B-class genes, and that a lack of SEPALLATA3-like gene expression in the outermost whorl may be related to the formation of greenish sepals.

Exploring potential new floral organ morphogenesis genes of Arabidopsis thaliana using systems biology approach

Flowering is one of the important defining features of angiosperms. The initiation of flower development and the formation of different floral organs are the results of the interplays among numerous genes. But until now, just fewer genes have been found linked with flower development. And the functions of lots of genes of Arabidopsis thaliana are still unknown. Although, the quartet model successfully simplified the ABCDE model to elaborate the molecular mechanism by introducing protein-protein interactions (PPIs). We still don't know much about several important aspects of flower development. So we need to discriminate even more genes involving in the flower development. In this study, we identified seven differentially modules through integrating the weighted gene co-expression network analysis (WGCNA) and Support Vector Machine (SVM) method to analyze co-expression network and PPIs using the public floral and non-floral expression profiles data of Arabidopsis thaliana. Gene set enrichment analysis was used for the functional annotation of the related genes, and some of the hub genes were identified in each module. The potential floral organ morphogenesis genes of two significant modules were integrated with PPI information in order to detail the inherent regulation mechanisms. Finally, the functions of the floral patterning genes were elucidated by combining the PPI and evolutionary information. It was indicated that the sub-networks or complexes, rather than the genes, were the regulation unit of flower development. We found that the most possible potential new genes underlining the floral pattern formation in A. thaliana were FY, CBL2, ZFN3, and AT1G77370; among them, FY, CBL2 acted as an upstream regulator of AP2; ZFN3 activated the flower primordial determining gene AP1 and AP2 by HY5/HYH gene via photo induction possibly. And AT1G77370 exhibited similar function in floral morphogenesis, same as ELF3. It possibly formed a complex between RFC3 and RPS15 in cytoplasm, which regulated TSO1 and CPSF160 in the nucleus, to control the floral organ morphogenesis. This process might also be fine tuning by AT5G53360 in the nucleus.

Plant development regulation: Overview and perspectives

Plant development, as occur in other eukaryotes, is conducted through a complex network of hormones, transcription factors, enzymes and micro RNAs, among other cellular components. They control developmental processes such as embryo, apical root and shoot meristem, leaf, flower, or seed formation, among others. The research in these topics has been very active in last decades. Recently, an explosion of new data concerning regulation mechanisms as well as the response of these processes to environmental changes has emerged. Initially, most of investigations were carried out in the model eudicot Arabidopsis but currently data from other plant species are available in the literature, although they are still limited. The aim of this review is focused on summarize the main molecular actors involved in plant development regulation in diverse plant species. A special attention will be given to the major families of genes and proteins participating in these regulatory mechanisms. The information on the regulatory pathways where they participate will be briefly cited. Additionally, the importance of certain structural features of such proteins that confer ductility and flexibility to these mechanisms will also be reported and discussed.

XAANTAL2 (AGL14) Is an Important Component of the Complex Gene Regulatory Network that Underlies Arabidopsis Shoot Apical Meristem Transitions

In Arabidopsis thaliana, multiple genes involved in shoot apical meristem (SAM) transitions have been characterized, but the mechanisms required for the dynamic attainment of vegetative, inflorescence, and floral meristem (VM, IM, FM) cell fates during SAM transitions are not well understood. Here we show that a MADS-box gene, XAANTAL2 (XAL2/AGL14), is necessary and sufficient to induce flowering, and its regulation is important in FM maintenance and determinacy. xal2 mutants are late flowering, particularly under short-day (SD) condition, while XAL2 overexpressing plants are early flowering, but their flowers have vegetative traits. Interestingly, inflorescences of the latter plants have higher expression levels of LFY, AP1, and TFL1 than wild-type plants. In addition we found that XAL2 is able to bind the TFL1 regulatory regions. On the other hand, the basipetal carpels of the 35S::XAL2 lines lose determinacy and maintain high levels of WUS expression under SD condition. To provide a mechanistic explanation for the complex roles of XAL2 in SAM transitions and the apparently paradoxical phenotypes of XAL2 and other MADS-box (SOC1, AGL24) overexpressors, we conducted dynamic gene regulatory network (GRN) and epigenetic landscape modeling. We uncovered a GRN module that underlies VM, IM, and FM gene configurations and transition patterns in wild-type plants as well as loss and gain of function lines characterized here and previously. Our approach thus provides a novel mechanistic framework for understanding the complex basis of SAM development.

Role of transcriptional regulation in the evolution of plant phenotype: A dynamic systems approach

A growing body of evidence suggests that alterations in transcriptional regulation of genes involved in modulating development are an important part of phenotypic evolution, and this can be documented among species and within populations. While the effects of differential transcriptional regulation in organismal development have been preferentially studied in animal systems, this phenomenon has also been addressed in plants. In this review, we summarize evidence for cis-regulatory mutations, trans-regulatory changes and epigenetic modifications as molecular events underlying important phenotypic alterations, and thus shaping the evolution of plant development. We postulate that a mechanistic understanding of why such molecular alterations have a key role in development, morphology and evolution will have to rely on dynamic models of complex regulatory networks that consider the concerted action of genetic and nongenetic components, and that also incorporate the restrictions underlying the genotype to phenotype mapping process. Developmental Dynamics 244:1074-1095, 2015. © 2015 Wiley Periodicals, Inc.

Evolution of the Plant Reproduction Master Regulators LFY and the MADS Transcription Factors: The Role of Protein Structure in the Evolutionary Development of the Flower

Understanding the evolutionary leap from non-flowering (gymnosperms) to flowering (angiosperms) plants and the origin and vast diversification of the floral form has been one of the focuses of plant evolutionary developmental biology. The evolving diversity and increasing complexity of organisms is often due to relatively small changes in genes that direct development. These "developmental control genes" and the transcription factors (TFs) they encode, are at the origin of most morphological changes. TFs such as LEAFY (LFY) and the MADS-domain TFs act as central regulators in key developmental processes of plant reproduction including the floral transition in angiosperms and the specification of the male and female organs in both gymnosperms and angiosperms. In addition to advances in genome wide profiling and forward and reverse genetic screening, structural techniques are becoming important tools in unraveling TF function by providing atomic and molecular level information that was lacking in purely genetic approaches. Here, we summarize previous structural work and present additional biophysical and biochemical studies of the key master regulators of plant reproduction - LEAFY and the MADS-domain TFs SEPALLATA3 and AGAMOUS. We discuss the impact of structural biology on our understanding of the complex evolutionary process leading to the development of the bisexual flower.

Survival and innovation: The role of mutational robustness in evolution

Biological systems are resistant to perturbations caused by the environment and by the intrinsic noise of the system. Robustness to mutations is a particular aspect of robustness in which the phenotype is resistant to genotypic variation. Mutational robustness has been linked to the ability of the system to generate heritable genetic variation (a property known as evolvability). It is known that greater robustness leads to increased evolvability. Therefore, mechanisms that increase mutational robustness fuel evolvability. Two such mechanisms, molecular chaperones and gene duplication, have been credited with enormous importance in generating functional diversity through the increase of system's robustness to mutational insults. However, the way in which such mechanisms regulate robustness remains largely uncharacterized. In this review, I provide evidence in support of the role of molecular chaperones and gene duplication in innovation. Specifically, I present evidence that these mechanisms regulate robustness allowing unstable systems to survive long periods of time, and thus they provide opportunity for other mutations to compensate the destabilizing effects of functionally innovative mutations. The findings reported in this study set new questions with regards to the synergy between robustness mechanisms and how this synergy can alter the adaptive landscape of proteins. The ideas proposed in this article set the ground for future research in the understanding of the role of robustness in evolution.

Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean

Similar to Arabidopsis thaliana, the wild soybeans (Glycine soja) and many cultivars exhibit indeterminate stem growth specified by the shoot identity gene Dt1, the functional counterpart of Arabidopsis TERMINAL FLOWER1 (TFL1). Mutations in TFL1 and Dt1 both result in the shoot apical meristem (SAM) switching from vegetative to reproductive state to initiate terminal flowering and thus produce determinate stems. A second soybean gene (Dt2) regulating stem growth was identified, which, in the presence of Dt1, produces semideterminate plants with terminal racemes similar to those observed in determinate plants. Here, we report positional cloning and characterization of Dt2, a dominant MADS domain factor gene classified into the APETALA1/SQUAMOSA (AP1/SQUA) subfamily that includes floral meristem (FM) identity genes AP1, FUL, and CAL in Arabidopsis. Unlike AP1, whose expression is limited to FMs in which the expression of TFL1 is repressed, Dt2 appears to repress the expression of Dt1 in the SAMs to promote early conversion of the SAMs into reproductive inflorescences. Given that Dt2 is not the gene most closely related to AP1 and that semideterminacy is rarely seen in wild soybeans, Dt2 appears to be a recent gain-of-function mutation, which has modified the genetic pathways determining the stem growth habit in soybean.

Impact of recurrent gene duplication on adaptation of plant genomes

BACKGROUND

Recurrent gene duplication and retention played an important role in angiosperm genome evolution. It has been hypothesized that these processes contribute significantly to plant adaptation but so far this hypothesis has not been tested at the genome scale.

RESULTS

We studied available sequenced angiosperm genomes to assess the frequency of positive selection footprints in lineage specific expanded (LSE) gene families compared to single-copy genes using a dN/dS-based test in a phylogenetic framework. We found 5.38% of alignments in LSE genes with codons under positive selection. In contrast, we found no evidence for codons under positive selection in the single-copy reference set. An analysis at the branch level shows that purifying selection acted more strongly on single-copy genes than on LSE gene clusters. Moreover we detect significantly more branches indicating evolution under positive selection and/or relaxed constraint in LSE genes than in single-copy genes.

CONCLUSIONS

In this - to our knowledge -first genome-scale study we provide strong empirical support for the hypothesis that LSE genes fuel adaptation in angiosperms. Our conservative approach for detecting selection footprints as well as our results can be of interest for further studies on (plant) gene family evolution.

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.

Expression of paralogous SEP-, FUL-, AG- and STK-like MADS-box genes in wild-type and peloric Phalaenopsis flowers

The diverse flowers of Orchidaceae are the result of several major morphological transitions, among them the most studied is the differentiation of the inner median tepal into the labellum, a perianth organ key in pollinator attraction. Type A peloria lacking stamens and with ectopic labella in place of inner lateral tepals are useful for testing models on the genes specifying these organs by comparing their patterns of expression between wild-type and peloric flowers. Previous studies focused on DEFICIENS- and GLOBOSA-like MADS-box genes because of their conserved role in perianth and stamen development. The "orchid code" model summarizes this work and shows in Orchidaceae there are four paralogous lineages of DEFICIENS/AP3-like genes differentially expressed in each floral whorl. Experimental tests of this model showed the conserved, higher expression of genes from two specific DEF-like gene lineages is associated with labellum development. The present study tests whether eight MADS-box candidate SEP-, FUL-, AG-, and STK-like genes have been specifically duplicated in the Orchidaceae and are also differentially expressed in association with the distinct flower organs of Phalaenopsis hyb. "Athens." The gene trees indicate orchid-specific duplications. In a way analogous to what is observed in labellum-specific DEF-like genes, a two-fold increase in the expression of SEP3-like gene PhaMADS7 was measured in the labellum-like inner lateral tepals of peloric flowers. The overlap between SEP3-like and DEF-like genes suggests both are associated with labellum specification and similar positional cues determine their domains of expression. In contrast, the uniform messenger levels of FUL-like genes suggest they are involved in the development of all organs and their expression in the ovary suggests cell differentiation starts before pollination. As previously reported AG-like and STK-like genes are exclusively expressed in gynostemium and ovary, however no evidence for transcriptional divergence was found in the stage investigated. Gene expression suggests a developmental regulatory system based on the combined activity of duplicate MADS-box genes. We discuss its feasibility based on documented protein interactions and patterns of expression.

Adaptive functional divergence of the warm temperature acclimation-related protein (WAP65) in fishes and the ortholog hemopexin (HPX) in mammals

Gene duplication is an important mechanism that leads to genetic novelty. Different, nonexclusive processes are likely involved, and many adaptive and nonadaptive events may contribute to the maintenance of duplicated genes. In some teleosts, a duplicate copy of the mammalian ortholog Hemopexin (HPX) is present, known as the warm temperature acclimation-related protein (WAP65). Both WAP65 and HPX have been associated with iron homeostasis due to the affinity to bind the toxic-free heme circulating in the blood stream. We have assessed the evolutionary dynamics of WAP65 and HPX genes to understand the adaptive role of positive selection at both nucleotide and amino acid level. Our results showed an asymmetrical evolution between the paralogs WAP65-1 and WAP65-2 after duplication with a slight acceleration of the evolutionary rate in WAP65-1, but not in WAP65-2, and few sites contributing to the functional distinction between the paralogs, whereas the majority of the protein remained under negative selection or relaxed negative selection. WAP65-1 is functionally more distinct from the ancestral protein function than WAP65-2. HPX is phylogenetically closer to WAP65-2 but even so functional divergence was detected between both proteins. In addition, HPX showed a fast rate of evolution when compared with both WAP65-1 and WAP65-2 genes. The assessed 3-dimensional (3-D) structure of WAP65-1 and WAP65-2 suggests that the functional differences detected are not causing noticeable structural changes in these proteins. However, such subtle changes between WAP65 paralogs may be important to understand the differential gene retention of both copies in 20 out of 30 teleosts species studied.

Molecular evolution constraints in the floral organ specification gene regulatory network module across 18 angiosperm genomes

The gene regulatory network of floral organ cell fate specification of Arabidopsis thaliana is a robust developmental regulatory module. Although such finding was proposed to explain the overall conservation of floral organ types and organization among angiosperms, it has not been confirmed that the network components are conserved at the molecular level among flowering plants. Using the genomic data that have accumulated, we address the conservation of the genes involved in this network and the forces that have shaped its evolution during the divergence of angiosperms. We recovered the network gene homologs for 18 species of flowering plants spanning nine families. We found that all the genes are highly conserved with no evidence of positive selection. We studied the sequence conservation features of the genes in the context of their known biological function and the strength of the purifying selection acting upon them in relation to their placement within the network. Our results suggest an association between protein length and sequence conservation, evolutionary rates, and functional category. On the other hand, we found no significant correlation between the strength of purifying selection and gene placement. Our results confirm that the studied robust developmental regulatory module has been subjected to strong functional constraints. However, unlike previous studies, our results do not support the notion that network topology plays a major role in constraining evolutionary rates. We speculate that the dynamical functional role of genes within the network and not just its connectivity could play an important role in constraining evolution.

Characterization of the sequence and expression pattern of LFY homologues from dogwood species (Cornus) with divergent inflorescence architectures

BACKGROUND AND AIMS

LFY homologues encode transcription factors that regulate the transition from vegetative to reproductive growth in flowering plants and have been shown to control inflorescence patterning in model species. This study investigated the expression patterns of LFY homologues within the diverse inflorescence types (head-like, umbel-like and inflorescences with elongated internodes) in closely related lineages in the dogwood genus (Cornus s.l.). The study sought to determine whether LFY homologues in Cornus species are expressed during floral and inflorescence development and if the pattern of expression is consistent with a function in regulating floral development and inflorescence architectures in the genus.

METHODS

Total RNAs were extracted using the CTAB method and the first-strand cDNA was synthesized using the SuperScript III first-strand synthesis system kit (Invitrogen). Expression of CorLFY was investigated by RT-PCR and RNA in situ hybridization. Phylogenetic analyses were conducted using the maximum likelihood methods implemented in RAxML-HPC v7.2.8.

KEY RESULTS

cDNA clones of LFY homologues (designated CorLFY) were isolated from six Cornus species bearing different types of inflorescence. CorLFY cDNAs were predicted to encode proteins of approximately 375 amino acids. The detection of CorLFY expression patterns using in situ RNA hybridization demonstrated the expression of CorLFY within the inflorescence meristems, inflorescence branch meristems, floral meristems and developing floral organ primordia. PCR analyses for cDNA libraries derived from reverse transcription of total RNAs showed that CorLFY was also expressed during the late-stage development of flowers and inflorescences, as well as in bracts and developing leaves. Consistent differences in the CorLFY expression patterns were not detected among the distinct inflorescence types.

CONCLUSIONS

The results suggest a role for CorLFY genes during floral and inflorescence development in dogwoods. However, the failure to detect expression differences between the inflorescence types in the Cornus species analysed suggests that the evolutionary shift between major inflorescence types in the genus is not controlled by dramatic alterations in the levels of CorLFY gene transcript accumulation. However, due to spatial, temporal and quantitative limitations of the expression data, it cannot be ruled out that subtle differences in the level or location of CorLFY transcripts may underlie the different inflorescence architectures that are observed across these species. Alternatively, differences in CorLFY protein function or the expression or function of other regulators (e.g. TFL1 and UFO homologues) may support the divergent developmental trajectories.

The MADS transcription factor XAL2/AGL14 modulates auxin transport during Arabidopsis root development by regulating PIN expression

Elucidating molecular links between cell-fate regulatory networks and dynamic patterning modules is a key for understanding development. Auxin is important for plant patterning, particularly in roots, where it establishes positional information for cell-fate decisions. PIN genes encode plasma membrane proteins that serve as auxin efflux transporters; mutations in members of this gene family exhibit smaller roots with altered root meristems and stem-cell patterning. Direct regulators of PIN transcription have remained elusive. Here, we establish that a MADS-box gene (XAANTAL2, XAL2/AGL14) controls auxin transport via PIN transcriptional regulation during Arabidopsis root development; mutations in this gene exhibit altered stem-cell patterning, root meristem size, and root growth. XAL2 is necessary for normal shootward and rootward auxin transport, as well as for maintaining normal auxin distribution within the root. Furthermore, this MADS-domain transcription factor upregulates PIN1 and PIN4 by direct binding to regulatory regions and it is required for PIN4-dependent auxin response. In turn, XAL2 expression is regulated by auxin levels thus establishing a positive feedback loop between auxin levels and PIN regulation that is likely to be important for robust root patterning.

Structural and functional evolution of positively selected sites in pine glutathione S-transferase enzyme family

Phylogenetic analyses have identified positive selection as an important driver of protein evolution, both structural and functional. However, the lack of appropriate combined functional and structural assays has generally hindered attempts to elucidate patterns of positively selected sites and their effects on enzyme activity and substrate specificity. In this study we investigated the evolutionary divergence of the glutathione S-transferase (GST) family in Pinus tabuliformis, a pine that is widely distributed from northern to central China, including cold temperate and drought-stressed regions. GSTs play important roles in plant stress tolerance and detoxification. We cloned 44 GST genes from P. tabuliformis and found that 26 of the 44 belong to the largest (Tau) class of GSTs and are differentially expressed across tissues and developmental stages. Substitution models identified five positively selected sites in the Tau GSTs. To examine the functional significance of these positively selected sites, we applied protein structural modeling and site-directed mutagenesis. We found that four of the five positively selected sites significantly affect the enzyme activity and specificity; thus their variation broadens the GST family substrate spectrum. In addition, positive selection has mainly acted on secondary substrate binding sites or sites close to (but not directly at) the primary substrate binding site; thus their variation enables the acquisition of new catalytic functions without compromising the protein primary biochemical properties. Our study sheds light on selective aspects of the functional and structural divergence of the GST family in pine and other organisms.

Congratulations, you have been carefully chosen to represent an important developmental regulator!

BACKGROUND

Studying a process in a new species often relies on focusing our attention to a candidate gene, encoding a protein similar to one with a known function. Not all the choices seem to be prudent.

SCOPE

This Viewpoint includes an overview of issues that are encountered during research of candidate genes. Defining a match for a gene of interest, deciding whether variation in ESTs or RNAseq data for a certain transcript, represent more than one gene. The problem of incorrect annotation of genes due to incorrect in-silico splicing, is also mentioned. The author's humble opinion on how to deal with these issues is provided.

CONCLUSIONS

The vast amount of new sequence data provides us with great possibilities for giant leaps in our understanding. Still, we cannot afford to skip over the tedious steps required to confirm that we are indeed studying the correct gene, and try to be sure that the complex expression pattern we observe is not a composite of several genes.

Molecular mapping of vernalization requirement and fertility restoration genes in carrot

Carrot (Daucus carota L.) is a cool-season vegetable normally classified as a biennial species, requiring vernalization to induce flowering. Nevertheless, some cultivars adapted to warmer climates require less vernalization and can be classified as annual. Most modern carrot cultivars are hybrids which rely upon cytoplasmic male-sterility for commercial production. One major gene controlling floral initiation and several genes restoring male fertility have been reported but none have been mapped. The objective of the present work was to develop the first linkage map of carrot locating the genomic regions that control vernalization response and fertility restoration. Using an F(2) progeny, derived from the intercross between the annual cultivar 'Criolla INTA' and a petaloid male sterile biennial carrot evaluated over 2 years, both early flowering habit, which we name Vrn1, and restoration of petaloid cytoplasmic male sterility, which we name Rf1, were found to be dominant traits conditioned by single genes. On a map of 355 markers covering all 9 chromosomes with a total map length of 669 cM and an average marker-to-marker distance of 1.88 cM, Vrn1 mapped to chromosome 2 with flanking markers at 0.70 and 0.46 cM, and Rf1 mapped to chromosome 9 with flanking markers at 4.38 and 1.12 cM. These are the first two reproductive traits mapped in the carrot genome, and their map location and flanking markers provide valuable tools for studying traits important for carrot domestication and reproductive biology, as well as facilitating carrot breeding.

Selaginella Genome Analysis - Entering the "Homoplasy Heaven" of the MADS World

In flowering plants, arguably the most significant transcription factors regulating development are MADS-domain proteins, encoded by Type I and Type II MADS-box genes. Type II genes are divided into the MIKC(C) and MIKC* groups. In angiosperms, these types and groups play distinct roles in the development of female gametophytes, embryos, and seeds (Type I); vegetative and floral tissues in sporophytes (MIKC(C)); and male gametophytes (MIKC*), but their functions in other plants are largely unknown. The complete set of MADS-box genes has been described for several angiosperms and a moss, Physcomitrella patens. Our examination of the complete genome sequence of a lycophyte, Selaginella moellendorffii, revealed 19 putative MADS-box genes (13 Type I, 3 MIKC(C), and 3 MIKC*). Our results suggest that the most recent common ancestor of vascular plants possessed at least two Type I and two Type II genes. None of the S. moellendorffii MIKC(C) genes were identified as orthologs of any floral organ identity genes. This strongly corroborates the view that the clades of floral organ identity genes originated in a common ancestor of seed plants after the lineage that led to lycophytes had branched off, and that expansion of MIKC(C) genes in the lineage leading to seed plants facilitated the evolution of their unique reproductive organs. The number of MIKC* genes and the ratio of MIKC* to MIKC(C) genes is lower in S. moellendorffii and angiosperms than in P. patens, correlated with reduction of the gametophyte in vascular plants. Our data indicate that Type I genes duplicated and diversified independently within lycophytes and seed plants. Our observations on MADS-box gene evolution echo morphological evolution since the two lineages of vascular plants appear to have arrived independently at similar body plans. Our annotation of MADS-box genes in S. moellendorffii provides the basis for functional studies to reveal the roles of this crucial gene family in basal vascular plants.

Overexpressing the ANR1 MADS-box gene in transgenic plants provides new insights into its role in the nitrate regulation of root development

The expression of the ANR1 MADS-box gene was manipulated in transgenic plants to investigate its role in the NO(3)(-)-dependent regulation of root development in Arabidopsis thaliana. Constitutive overexpression of ANR1 in roots, achieved using GAL4 enhancer trap lines, resulted in more rapid early seedling development, increased lengths and numbers of lateral roots and increased shoot fresh weight. Based on results obtained with five different enhancer trap lines, the overexpression of ANR1 in the lateral root tips appears to be more important for this phenotype than its level of expression in the developing lateral root primordia. Dexamethasone-mediated induction of ANR1 in lines expressing an ANR1-GR (glucocorticoid receptor) fusion protein stimulated lateral root growth but not primary root growth. Short-term (24 h) dexamethasone treatments led to prolonged stimulation of lateral root growth, whether the lateral roots were already mature or still unemerged at the time of treatment. In split-root experiments, localized application of dexamethasone to half of the root system of an ANR1-GR line elicited a localized increase in both the length and numbers of lateral roots, mimicking the effect of a localized NO(3)(-) treatment. In both types of transgenic line, the root phenotype was strongly dependent on the presence of NO(3)(-), indicating that there are additional components involved in ANR1 function that are NO(3)(-) regulated. The implications of these results for our understanding of ANR1's mode of action in the root response to localized NO(3)(-) are discussed.

Emerging knowledge from genome sequencing of crop species

Extensive insights into the genome composition, organization, and evolution have been gained from the plant genome sequencing and annotation ongoing projects. The analysis of crop genomes provided surprising evidences with important implications in plant origin and evolution: genome duplication, ancestral re-arrangements and unexpected polyploidization events opened new doors to address fundamental questions related to species proliferation, adaptation, and functional modulations. Detailed paleogenomic analysis led to many speculation on how chromosomes have been shaped over time in terms of gene content and order. The completion of the genome sequences of several major crops, prompted to a detailed identification and annotation of transposable elements: new hypothesis related to their composition, chromosomal distribution, insertion models, amplification rate, and evolution patterns are coming up. Availability of full genome sequence of several crop species as well as from many accessions within species is providing new keys for biodiversity exploitation and interpretation. Re-sequencing is enabling high-throughput genotyping to identify a wealth of SNP and afterward to produce haplotype maps necessary to accurately associate molecular variation to phenotype. Conservation genomics is emerging as a powerful tool to explain adaptation, genetic drift, natural selection, hybridization and to estimate genetic variation, fitness and population's viability.

When ABC becomes ACB

Understanding how the information contained in genes is mapped onto the phenotypes, and deriving formal frameworks to search for generic aspects of developmental constraints and evolution remains one of the main challenges of contemporary biological research. The Mexican endemic triurid Lacandonia schismatica (Lacandoniaceae), a mycoheterotrophic monocotyledonous plant with hermaphroditic reproductive axes is alone among 250,000 species of angiosperms, as it has central stamens surrounded by a peripheral gynoecium, representing a natural instance of a homeotic mutant. Based on the classical ABC model of flower development, it has recently been shown that the B-function gene APETALA3 (AP3), essential for stamen identity, was displaced toward the flower centre in L. schismatica (ABC to ACB) from the early stages of flower development. A functional conservation of B-function genes from L. schismatica through the rescue of B-gene mutants in Arabidopsis thaliana, as well as conserved protein interactions, has also been demonstrated. Thus, it has been shown that relatively simple genetic alterations may underlie large morphological shifts fixed in extant natural populations. Nevertheless, critical questions remain in order to have a full and sufficient explanation of the molecular genetic mechanisms underlying L. schismatica's unique floral arrangement. Evolutionary approaches to developmental mechanisms and systems biology, including high-throughput functional genomic studies and models of complex developmental gene regulatory networks, constitute two main approaches to meet such a challenge. In this review, the aim is to address some of the pending questions with the ultimate goal of investigating further the mechanisms of L. schismatica's unique homeotic flower arrangement and its evolution.

The plant proteome folding project: structure and positive selection in plant protein families

Despite its importance, relatively little is known about the relationship between the structure, function, and evolution of proteins, particularly in land plant species. We have developed a database with predicted protein domains for five plant proteomes (http://pfp.bio.nyu.edu) and used both protein structural fold recognition and de novo Rosetta-based protein structure prediction to predict protein structure for Arabidopsis and rice proteins. Based on sequence similarity, we have identified ∼15,000 orthologous/paralogous protein family clusters among these species and used codon-based models to predict positive selection in protein evolution within 175 of these sequence clusters. Our results show that codons that display positive selection appear to be less frequent in helical and strand regions and are overrepresented in amino acid residues that are associated with a change in protein secondary structure. Like in other organisms, disordered protein regions also appear to have more selected sites. Structural information provides new functional insights into specific plant proteins and allows us to map positively selected amino acid sites onto protein structures and view these sites in a structural and functional context.

Agave tequilana MADS genes show novel expression patterns in meristems, developing bulbils and floral organs

Agave tequilana is a monocarpic perennial species that flowers after 5-8 years of vegetative growth signaling the end of the plant's life cycle. When fertilization is unsuccessful, vegetative bulbils are induced on the umbels of the inflorescence near the bracteoles from newly formed meristems. Although the regulation of inflorescence and flower development has been described in detail for monocarpic annuals and polycarpic species, little is known at the molecular level for these processes in monocarpic perennials, and few studies have been carried out on bulbils. Histological samples revealed the early induction of umbel meristems soon after the initiation of the vegetative to inflorescence transition in A. tequilana. To identify candidate genes involved in the regulation of floral induction, a search for MADS-box transcription factor ESTs was conducted using an A. tequilana transcriptome database. Seven different MIKC MADS genes classified into 6 different types were identified based on previously characterized A. thaliana and O. sativa MADS genes and sequences from non-grass monocotyledons. Quantitative real-time PCR analysis of the seven candidate MADS genes in vegetative, inflorescence, bulbil and floral tissues uncovered novel patterns of expression for some of the genes in comparison with orthologous genes characterized in other species. In situ hybridization studies using two different genes showed expression in specific tissues of vegetative meristems and floral buds. Distinct MADS gene regulatory patterns in A. tequilana may be related to the specific reproductive strategies employed by this species.

Robustness and evolvability in the B-system of flower development

BACKGROUND

Gene duplication has often been invoked as a key mechanism responsible for evolution of new morphologies. The floral homeotic B-group gene family has undergone a number of gene duplication events, and yet the functions of these genes appear to be largely conserved. However, detailed comparative analysis has indicated that such duplicate genes have considerable cryptic variability in their functions. In the Solanaceae, two duplicate B-group gene lineages have been retained in three subfamilies. Comparisons of orthologous genes across members of the Solanaceae have demonstrated that the combined function of all four B-gene members is to establish petal and stamen identity, but that this function was partitioned differently in each species. These observations emphasize both the robustness and the evolvability of the B-system.

SCOPE

We provide an overview of how the B-function genes can robustly specify petal and stamen identity and at the same time evolve through changes in protein-protein interaction, gene expression patterns, copy number variation or alterations in the downstream genes they control. By using mathematical models we explore regulatory differences between species and how these impose constraints on downstream gene regulation.

CONCLUSIONS

Evolvability of the B-genes can be understood through the multiple ways in which the B-system can be robust. Quantitative approaches should allow for the incorporation of more biological realism in the representations of these regulatory systems and this should contribute to understanding the constraints under which different B-systems can function and evolve. This, in turn, can provide a better understanding of the ways in which B-genes have contributed to flower diversity.

Vertebrate paralogous MEF2 genes: origin, conservation, and evolution

Background

The myocyte enhancer factor 2 (MEF2) gene family is broadly expressed during the development and maintenance of muscle cells. Although a great deal has been elucidated concerning MEF2 transcription factors' regulation of specific gene expression in diverse programs and adaptive responses, little is known about the origin and evolution of the four members of the MEF2 gene family in vertebrates.

Methodology/Principal Findings

By phylogenetic analyses, we investigated the origin, conservation, and evolution of the four MEF2 genes. First, among the four MEF2 paralogous branches, MEF2B is clearly distant from the other three branches in vertebrates, mainly because it lacks the HJURP_C (Holliday junction recognition protein C-terminal) region. Second, three duplication events might have occurred to produce the four MEF2 paralogous genes and the latest duplication event occurred near the origin of vertebrates producing MEF2A and MEF2C. Third, the ratio (K_a/K_s) of non-synonymous to synonymous nucleotide substitution rates showed that MEF2B evolves faster than the other three MEF2 proteins despite purifying selection on all of the four MEF2 branches. Moreover, a pair model of M0 versus M3 showed that variable selection exists among MEF2 proteins, and branch-site analysis presented that sites 53 and 64 along the MEF2B branch are under positive selection. Finally, and interestingly, substitution rates showed that type II MADS genes (i.e., MEF2-like genes) evolve as slowly as type I MADS genes (i.e., SRF-like genes) in animals, which is inconsistent with the fact that type II MADS genes evolve much slower than type I MADS genes in plants.

Conclusion

Our findings shed light on the relationship of MEF2A, B, C, and D with functional conservation and evolution in vertebrates. This study provides a rationale for future experimental design to investigate distinct but overlapping regulatory roles of the four MEF2 genes in various tissues.

The emerging importance of type I MADS box transcription factors for plant reproduction

Based on their evolutionary origin, MADS box transcription factor genes have been divided into two classes, namely, type I and II. The plant-specific type II MIKC MADS box genes have been most intensively studied and shown to be key regulators of developmental processes, such as meristem identity, flowering time, and fruit and seed development. By contrast, very little is known about type I MADS domain transcription factors, and they have not attracted interest for a long time. A number of recent studies have now indicated a key regulatory role for type I MADS box factors in plant reproduction, in particular in specifying female gametophyte, embryo, and endosperm development. These analyses have also suggested that type I MADS box factors are decisive for setting reproductive boundaries between species.

Conserved and variable correlated mutations in the plant MADS protein network

BACKGROUND

Plant MADS domain proteins are involved in a variety of developmental processes for which their ability to form various interactions is a key requisite. However, not much is known about the structure of these proteins or their complexes, whereas such knowledge would be valuable for a better understanding of their function. Here, we analyze those proteins and the complexes they form using a correlated mutation approach in combination with available structural, bioinformatics and experimental data.

RESULTS

Correlated mutations are affected by several types of noise, which is difficult to disentangle from the real signal. In our analysis of the MADS domain proteins, we apply for the first time a correlated mutation analysis to a family of interacting proteins. This provides a unique way to investigate the amount of signal that is present in correlated mutations because it allows direct comparison of mutations in various family members and assessing their conservation. We show that correlated mutations in general are conserved within the various family members, and if not, the variability at the respective positions is less in the proteins in which the correlated mutation does not occur. Also, intermolecular correlated mutation signals for interacting pairs of proteins display clear overlap with other bioinformatics data, which is not the case for non-interacting protein pairs, an observation which validates the intermolecular correlated mutations. Having validated the correlated mutation results, we apply them to infer the structural organization of the MADS domain proteins.

CONCLUSION

Our analysis enables understanding of the structural organization of the MADS domain proteins, including support for predicted helices based on correlated mutation patterns, and evidence for a specific interaction site in those proteins.

Tracing the evolution of the floral homeotic B- and C-function genes through genome synteny

The evolution of the floral homeotic genes has been characterized using phylogenetic and functional studies. It is possible to enhance these studies by comparing gene content and order between species to determine the evolutionary history of the regulatory genes. Here, we use a synteny-based approach to trace the evolution of the floral B- and C-function genes that are required for specification of the reproductive organs. Consistent with previous phylogenetic studies, we show that the euAP3–TM6 split occurred after the monocots and dicots diverged. The Arabidopsis TM6 and papaya euAP3 genes are absent from the respective genomes, and we have detected loci from which these genes were lost. These data indicate that either the TM6 or the euAP3 lineage genes can be lost without detriment to flower development. In contrast, PI is essential for male reproductive organ development; yet, contrary to predictions, complex genomic rearrangements have resulted in almost complete breakdown of synteny at the PI locus. In addition to showing the evolution of B-function genes through the prediction of ancestral loci, similar reconstructions reveal the origins of the C-function AG and PLE lineages in dicots, and show the shared ancestry with the monocot C-function genes. During our studies, we found that transposable elements (TEs) present in sequenced Antirrhinum genomic clones limited comparative studies. A pilot survey of the Antirrhinum data revealed that gene-rich regions contain an unusually high degree of TEs of very varied types, which will be an important consideration for future genome sequencing efforts.

The Arabidopsis B-sister MADS-box protein, GORDITA, represses fruit growth and contributes to integument development

The MADS-box family of transcription factors have diverse developmental roles in flower pattern formation, gametophyte cell division and fruit differentiation. The B-sister MADS-box proteins are most similar to the B-class floral homeotic proteins, and are expressed in female reproductive organs. The Arabidopsis B-sister MADS-box protein, TT16, is necessary for inner integument differentiation. We have functionally characterized the only other B-sister MADS-box gene in Arabidopsis, AGL63, renamed here as GORDITA (GOA). A loss-of-function mutation in goa or reduction of endogenous GOA expression results in larger fruits, illustrating its novel function in regulating fruit growth. Consistent with its function, GOA expression is detected in the walls of the valves and throughout the replum of the fruit. Our phenotypic and molecular analyses of 35S::GOA and goa plants show that GOA controls organ size via cell expansion. Further, functional studies of goa tt16 double mutants have shown their additive role in controlling seed coat development, and have revealed the importance of GOA expression in the outer integument. Together, our studies provide evidence of a new regulatory role for a B-sister MADS-box gene in the control of organ growth.

Flower development

Flowers are the most complex structures of plants. Studies of Arabidopsis thaliana, which has typical eudicot flowers, have been fundamental in advancing the structural and molecular understanding of flower development. The main processes and stages of Arabidopsis flower development are summarized to provide a framework in which to interpret the detailed molecular genetic studies of genes assigned functions during flower development and is extended to recent genomics studies uncovering the key regulatory modules involved. Computational models have been used to study the concerted action and dynamics of the gene regulatory module that underlies patterning of the Arabidopsis inflorescence meristem and specification of the primordial cell types during early stages of flower development. This includes the gene combinations that specify sepal, petal, stamen and carpel identity, and genes that interact with them. As a dynamic gene regulatory network this module has been shown to converge to stable multigenic profiles that depend upon the overall network topology and are thus robust, which can explain the canalization of flower organ determination and the overall conservation of the basic flower plan among eudicots. Comparative and evolutionary approaches derived from Arabidopsis studies pave the way to studying the molecular basis of diverse floral morphologies.