AGAMOUS-LIKE6 (AGL6) genes play essential roles in flower development, but whether and how they work with floral organ identity genes remain less understood. Here, we describe interactions of the rice (Oryza sativa) AGL6 gene MADS6 with other rice floral homeotic genes in flower development. Genetic analyses revealed that MADS6 specifies the identity of the three inner whorls and floral meristem determinacy redundantly with SUPERWOMAN1/MADS16 (B-gene) or MADS3 (C-gene). MADS6 was shown to define carpel/ovule development and floral determinacy by interacting with MADS13 (D-gene) and control the palea and floral meristem identities together with the YABBY gene DROOPING LEAF. Expression analyses revealed that the transcript levels of six B-, C-, and E-class genes were reduced in mads6-1 at the early flower developmental stage, suggesting that MADS6 is a key regulator of early flower development. Moreover, MADS6 can directly bind to a putative regulatory motif on MADS58 (C-gene), and mads6-1 mads58 displayed phenotypes similar to that of mads6-1. These results suggest that MADS6 is a key player in specifying flower development via interacting with other floral homeotic genes in rice, thus providing new insights into the mechanism by which flower development is controlled.

Grass plants develop unique floral patterns that determine grain production. However, the molecular mechanism underlying the specification of floral organ identities and meristem determinacy, including the interaction among floral homeotic genes, remains largely unknown in grasses. Here, we report the interactions of rice (Oryza sativa) floral homeotic genes, OsMADS3 (a C-class gene), OsMADS13 (a D-class gene), and DROOPING LEAF (DL), in specifying floral organ identities and floral meristem determinacy. The interaction among these genes was revealed through the analysis of double mutants. osmads13-3 osmads3-4 displayed a loss of floral meristem determinacy and generated abundant carpelloid structures containing severe defective ovules in the flower center, which were not detectable in the single mutant. In addition, in situ hybridization and yeast two-hybrid analyses revealed that OsMADS13 and OsMADS3 did not regulate each other’s transcription or interact at the protein level. This indicates that OsMADS3 plays a synergistic role with OsMADS13 in both ovule development and floral meristem termination. Strikingly, osmads3-4 dl-sup6 displayed a severe loss of floral meristem determinacy and produced supernumerary whorls of lodicule-like organs at the forth whorl, suggesting that OsMADS3 and DL synergistically terminate the floral meristem. Furthermore, the defects of osmads13-3 dl-sup6 flowers appeared identical to those of dl-sup6, and the OsMADS13 expression was undetectable in dl-sup6 flowers. These observations suggest that DL and OsMADS13 may function in the same pathway specifying the identity of carpel/ovule and floral meristem. Collectively, we propose a model to illustrate the role of OsMADS3, DL, and OsMADS13 in the specification of flower organ identity and meristem determinacy in rice.

As a superfamily of transcription factors (TFs), the basic helix-loop-helix (bHLH) proteins have been characterized functionally in many plants with a vital role in the regulation of diverse biological processes including growth, development, response to various stresses, and so on. However, no systemic analysis of the bHLH TFs has been reported in Brachypodium distachyon, an emerging model plant in Poaceae. A total of 146 bHLH TFs were identified in the Brachypodium distachyon genome and classified into 24 subfamilies. BdbHLHs in the same subfamily share similar protein motifs and gene structures. Gene duplication events showed a close relationship to rice, maize and sorghum, and segment duplications might play a key role in the expansion of this gene family. The amino acid sequence of the bHLH domains were quite conservative, especially Leu-27 and Leu-54. Based on the predicted binding activities, the BdbHLHs were divided into DNA binding and non-DNA binding types. According to the gene ontology (GO) analysis, BdbHLHs were speculated to function in homodimer or heterodimer manner. By integrating the available high throughput data in public database and results of quantitative RT-PCR, we found the expression profiles of BdbHLHs were different, implying their differentiated functions. One hundred fourty-six BdbHLHs were identified and their conserved domains, sequence features, phylogenetic relationship, chromosomal distribution, GO annotations, gene structures, gene duplication and expression profiles were investigated. Our findings lay a foundation for further evolutionary and functional elucidation of BdbHLH genes.

MYB transcription factors are widespread in plants and play key roles in plant development. Although MYB transcription factors have been thoroughly characterized in many plants, genome-wide analysis of the MYB gene family has not yet been undertaken in Brachypodium distachyon. In this study, 122 BdMYB transcription factors were identified, comprising 85 MYB-R2R3, 34 MYB-related and three MYB-R1R2R3. Phylogenetic analysis showed that BdMYBs, OsMYBs and AtMYBs with similar functions were clustered in the same subgroup, and the phylogenetic relationships of BdMYB transcription factors were supported by highly conserved motifs and gene structures. Two cis-elements were found in the promoters of BdMYB genes. One is related to plant growth/development, the other is related to stress responses. Gene Ontology (GO) analysis indicated that most of the BdMYB genes are involved in various biological processes. The chromosome distribution pattern strongly indicated that genome-wide tandem and segment duplication mainly contributed to the expansion of the BdMYB gene family. Synteny analysis showed that 56, 58 and 61 BdMYB genes were orthologous to rice, maize and sorghum, respectively. We further demonstrated that BdMYB genes have evolved under strong purifying selection. The expression profiles indicated that most BdMYB genes might participate in floral development and respond to abiotic stresses. Additionally, 338 pairs of proteins were predicted to interact by constructing the interaction network. This work laid the foundation and provided clues for understanding the biological functions of these transcription factors.

Heat stress greatly affects plant growth/development and influences the output of crops. With the increased occurrence of extreme high temperature, the negative influence on cereal products from heat stress becomes severer and severer. It is urgent to reveal the molecular mechanism in response to heat stress in plants. In this research, we used RNA-seq technology to identify differentially expressed genes (DEGs) in leaves of seedlings, leaves and inflorescences at heading stage of Brachypodium distachyon, one model plant of grasses. Results showed many genes in responding to heat stress. Of them, the expression level of 656 DEGs were altered in three groups of samples treated with high temperature. Gene ontology (GO) analysis showed that the highly enriched DEGs were responsible for heat stress and protein folding. According to KEGG pathway analysis, the DEGs were related mainly to photosynthesis-antenna proteins, the endoplasmic reticulum, and the spliceosome. Additionally, the expression level of 454 transcription factors belonging to 49 gene families was altered, as well as 1,973 splicing events occurred after treatment with high temperature. This research lays a foundation for characterizing the molecular mechanism of heat stress response and identifying key genes for those responses in plants. These findings also clearly show that heat stress regulates the expression of genes not only at transcriptional level, but also at post-transcriptional level.