Novel members of the AGAMOUS LIKE 6 subfamily of MIKCC-type MADS-box genes in soybean
© Wong et al.; licensee BioMed Central Ltd. 2013
Received: 1 November 2012
Accepted: 11 July 2013
Published: 20 July 2013
The classical (C) MIKC-type MADS-box transcription factors comprise one gene family that plays diverse roles in the flowering process ranging from floral initiation to the development of floral organs. Despite their importance in regulating developmental processes that impact crop yield, they remain largely unexplored in the major legume oilseed crop, soybean.
We identified 57 MIKCc-type transcription factors from soybean and determined the in silico gene expression profiles of the soybean MIKCc-type genes across different tissues. Our study implicates three MIKCc-type transcription factors as novel members of the AGAMOUS LIKE 6 (AGL6) subfamily of the MIKCC-type MADS-box genes, and we named this sister clade PsMADS3. While similar genes were identified in other legume species, poplar and grape, no such gene is represented in Arabidopsis thaliana or rice. RT-PCR analysis on these three soybean PsMADS3 genes during early floral initiation processes revealed their temporal expression similar to that of APETALA1, a gene known to function as a floral meristem identity gene. However, RNA in situ hybridisation showed that their spatial expression patterns are markedly different from those of APETALA1.
Legume flower development system differs from that in the model plant, Arabidopsis. There is an overlap in the initiation of different floral whorls in soybean, and inflorescent meristems can revert to leaf production depending on the environmental conditions. MIKCC-type MADS-box genes have been shown to play key regulatory roles in different stages of flower development. We identified members of the PsMADS3 sub-clade in legumes that show differential spatial expression during floral initiation, indicating their potential novel roles in the floral initiation process. The results from this study will contribute to a better understanding of legume-specific floral developmental processes.
KeywordsDevelopment Floral meristem MADS-box transcription factor Soybean
Flower development in plants involves tightly regulated processes starting from floral initiation to flower formation. The underlying processes have been extensively investigated, as flower development is an important agronomic trait that determines crop yield. Various transcription factors are essential in regulating these developmental processes, including the family of MADS-box transcription factors.
The MADS-box transcription factors, especially the plant-specific classical (C) MIKC-type MADS-box genes, are known to play key regulatory roles in different stages of flower development. Their roles in coordinating floral developmental processes have been revealed by functional studies largely carried out in the model plant, Arabidopsis thaliana. The MIKCC-type genes are characterised by a conserved structural organisation of the MADS-box, Intervening-, Keratin-like- and C-domains. The highly conserved MADS-domain and the weakly conserved I-domain are required for DNA binding, while the strongly conserved K-domain and the variable C-domain regulate protein interactions .
Genome-wide analyses of the MIKCC-type genes have been carried out in Arabidopsis , rice  and poplar . While Arabidopsis and rice genomes have similar numbers of MIKCC-type genes (39 vs. 38), poplar has 55 of these genes, suggesting a higher birth rate compared to Arabidopsis or rice. These MIKCc-type genes can be divided into 15 distinct gene clades with each clade named after the first member identified [3, 5]. All but two (TM8 and OsMADS32) are found in Arabidopsis [3, 5], and the FLC clade may be absent from the rice genome . It remains unclear whether all clades are present in the poplar genome, as no TM8 genes were used in the phylogenetic analysis .
The SQUA subfamily clade includes four Arabidopsis members, APETALA1 (AP1), CAULIFLOWER (CAL), FRUITFULL (FUL) and AGAMOUS-LIKE 79 (AGL79) . The functions of AP1, CAL and FUL have been characterised, indicating that they play partially redundant roles in determining floral meristem identity . The SEPALLATA (SEP) family belongs to the AGL2 clade, and there are four members documented in Arabidopsis . All four members (SEP1, SEP2, SEP3 and SEP4) play redundant functions in determining floral organ identity and floral meristem determinacy. AP1 has been shown to bind directly to the SEP3 promoter, hence increasing the expression of SEP3 rapidly . The AGL6 subfamily has a relatively small representation (only two genes, AGL6 and its paralog AGL13) in Arabidopsis. While no knockout phenotype has been described for either of these genes in Arabidopsis, studies in rice, maize and Petunia hybrida have largely demonstrated the roles of the AGL6 subfamily in regulating floral organ identity and floral meristem determinacy, indicating redundant roles with closely related genes including SEP[9–11]. A phylogenetic study showed that subfamilies of SQUA, SEP and AGL6 are always rooted together in one superclade, which may be correlated with their overlapping roles in regulating flower development.
A total of 212 MADS-box genes were predicted in the recent genome sequence of soybean . Earlier we reported the diversification of some gene expression and microRNAs in legume SAM [13–16]. However, much remains to be learned about these genes, especially given their potential impact on crop production. Soybean is the largest legume crop in the world and accounts for greater than 50% of the global oilseed production. In this study, we identified all the soybean MIKCc-type MADS-box genes using the current Glyma1.0 gene set and identified potential phylogenetic relationships to their Arabidopsis, rice and poplar counterparts. We examined the expression patterns across different soybean tissues for the entire family. Intriguingly, the results revealed a novel AGL6 sister clade of MIKCc-type genes in soybean, and we focused our subsequent analysis on members of this novel sub-clade.
Results and discussion
Molecular evolutionary analysis of soybean MIKCc-type MADS-box transcription factors
When we searched the soybean predicted gene set available at Phytozome  for sequences containing both a MADS-box and K-domain, we identified a total of 57 sequences. Subsequent inspection revealed three of the sequences were incomplete. We attempted to obtain a full-length sequence for these genes using gene prediction software on the genome sequence surrounding these partial sequences but did not yield any results. Therefore, we omitted these sequences from further analysis. To investigate their phylogenetic relationships with MIKCc-type genes from Arabidopsis, rice and poplar, reported MIKCc group protein sequences [3–5] were retrieved from their respective databases. A total of 159 conceptually translated protein sequences were used in the phylogenetic analysis.
In silico analysis of soybean MIKCcgene expression patterns
All identified soybean MIKCc-genes are expressed in at least one of the three reproductive tissues represented (reproductive SAM, flower and pod), except for members of the AGL12 clade, which are only expressed in the root (Figure 3). Arabidopsis AGL12 is preferentially expressed in the root, and recent loss-of-function analyses have revealed its roles in not only regulating root meristem cell proliferation but also flowering transition [20, 21]. Based on the soybean AGL12-LIKE expression profile, it is tempting to speculate that their functions in floral regulation may have been lost. A similar expression pattern was observed for most MIKCc-genes clustered within a clade. All duplicated genes are transcribed and have comparable expression profiles, especially in the reproductive SAM (Figure 3), suggesting their functional significance.
As for the superclade consisting of AGL2, SQUA, and AGL6, there are some notable differences in the gene expression profiles among the three clades (Figure 3). For example, all AGL2-LIKE genes except one (Glyma01g08130) are absent from the SAM during the early floral transition process but are expressed later in the floral developmental process in the flower and pod. This pattern is expected as these genes are known to be activated following AP1 induction in Arabidopsis . The phylogenetic tree indicates that Glyma01g08130 is the counterpart for Arabidopsis SEP4. In addition to being found in flower and pod like the rest of the AGL2-LIKE genes, it is also expressed in the SAMs during the floral initiation process and very highly in nodules. This pattern implies a likely diverged function of GmSEP4 with additional roles in the early floral initiation process as well as in nodule formation. Glyma01g08150, one of the four soybean counterparts of Arabidopsis AP1, also likely plays a role in nodule formation. Although the expression of Glyma01g08150 is drastically induced on 4SD in the SAM (20 RPKM), the level of expression is 6-fold less than that in the nodule (134 RPKM; Figure 3). Intriguingly, its homeolog Glyma02g13420 is not expressed in the nodule but rather has the highest expression in the reproductive SAM (105 RPKM; Figure 1 & 3), suggesting a functional divergence between this homeolog pair.
Although members of soybean AGL6 genes are expressed in the reproductive SAMs, changes in their transcript levels are not comparable with those of PsMADS3-LIKE and SQUA-LIKE genes during the early floral transition process, suggesting that the latter two clades are likely to play more prominent roles in the developmental transition process. Because there is no information available for PsMADS3-LIKE genes, we focused our study on members of this novel sister clade.
Expression of GmMADS3 during the floral initiation process
The expression of Glyma16g32540 is distinct from that of GmAP1 and GmMADS3. A weak signal associated with its expression was detected in the centre of the inflorescence meristem (Figure 5g); on 6SD, its expression was also observed in the centre of the newly emerged floral meristem (Figure 5h). The expression of Glyma16g32540 in the centre of the meristem indicates its potential regulatory roles in orchestrating events in the inflorescence meristem. The spatial expression pattern of the soybean PsMADS3-LIKE genes supports that these genes are novel, as their expression differs markedly from the spatial expression of closely related family members such as GmAP1 (this study) or Arabidopsis AGL6.
In contrast to Arabidopsis where the initiation timing of floral whorls does not overlap, the legume soybean has a flower development system with overlapping whorls . Furthermore, unlike Arabidopsis that usually cannot undergo flowering reversion , the soybean inflorescent meristem can revert to leaf production when the environmental growth conditions are switched from SD to LD . Because the MIKCC-type MADS-box genes play key regulatory roles in different stages of flower development, it is conceivable that members of the PsMADS3 sub-clade identified in this study could contribute to developmental plasticity in cooperation with key floral regulators such as GmAP1 or GmFLC. Future studies aimed at defining the interacting partners of these genes will aid in our understanding of the floral transition process.
Sequence and phylogenetic analysis
Conceptually translated protein sequences were retrieved from public databases (Phytozome, Rice Genome Annotation Project, TAIR and LjGDB). For the initial identification of the soybean MIKCC-type MADS-box transcription factors, all annotated genes were screened for both the MADS-box domain (PFAM00319) and K-domain (PFAM01486). The results were then manually inspected and filtered for truncated protein sequences, resulting in a total of 54 sequences (Additional file 1: Table S1). Sequences were imported into MEGA version 5 software for subsequent phylogenetic and molecular evolutionary analyses . MUSCLE alignments of protein sequences spanning the MADS-, I- and K-domains were carried out using the default settings in MEGA. After alignment, the evolutionary history was estimated using the Maximum Likelihood method based on the JTT matrix-based model as implemented in Mega 5 with bootstrap analysis set at 200 replicates.
For the expression profile analysis, two separate transcriptome sequencing data were used [19, 29]. The abundance for each gene was normalised within each dataset and expressed in reads per kb per million (RPKM) and are provided in Additional file 1.
Plant growth and RNA extraction
Soybean plants [Glycine max. (L) Merr. Cv. Bragg] were grown in a greenhouse located at the University of Melbourne, Victoria, Australia. To induce flowering, 10-day-old plants were shifted to a growth chamber maintained at a constant temperature of 25°C with a 10-hr day (150 μmol m-2 s-1) and 14-hr night (short-day). Shoot apical meristems (SAMs) were micro-dissected, as previously described . Total RNA was extracted from dissected SAM (approximately 80 SAMs per extraction) using the Qiagen RNeasy Mini Kit (Qiagen, Victoria, Australia) with on-column DNAse digestion.
The Qiagen one-step RT-PCR kit was used according to the manufacturer’s instructions in all RT-PCR analyses. Total RNA (20 ng) isolated from the SAMs of 10-day-old soybean seedlings (0 SD) and from the SAMs of plants subjected to different short-day treatments (2, 4 or 6 SD) was used as the template in a 10-μl reaction volume for 25 amplification cycles. Primers used are:
The soybean actin gene was used as an internal control. The PCR reactions were separated on 1% agarose gels containing 0.1 μg/μl ethidium bromide and visualised under UV light.
RNA in situhybridisation
The soybean shoot apices were dissected and fixed with 4% paraformaldehyde (Sigma, Victoria, Australia) in phosphate-buffered saline overnight at 4°C after vacuum infiltration. Subsequent fixation and hybridisation steps were followed as previously described .
We wish to thank Professor Bernie Carroll for soybean seeds and Dr. Sheh May Tam and Dr. Lim Chee Liew for help with the phylogenetic analyses. Financial support from the Australian Research Council in the form of the ARC Centre of Excellence for Integrative Legume Research (CE0348212) andARC DP0988972 is also gratefully acknowledged.
- Kaufmann K, Melzer R, Theißen G: MIKC-type MADS-domain proteins: structural modularity, protein interactions and network evolution in land plants. Gene. 2005, 347 (2): 183-198. 10.1016/j.gene.2004.12.014.PubMedView ArticleGoogle Scholar
- Parenicova L, de Folter S, Kieffer M, Horner DS, Favalli C, Busscher J, Cook HE, Ingram RM, Kater MM, Davies B, et al: Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: New openings to the MADS world. Plant Cell. 2003, 15 (7): 1538-1551. 10.1105/tpc.011544.PubMedPubMed CentralView ArticleGoogle Scholar
- Arora R, Agarwal P, Ray S, Singh A, Singh V, Tyagi A, Kapoor S: MADS-box gene family in rice: genome-wide identification, organization and expression profiling during reproductive development and stress. BMC Genomics. 2007, 8 (1): 242-10.1186/1471-2164-8-242.PubMedPubMed CentralView ArticleGoogle Scholar
- Leseberg CH, Li AL, Kang H, Duvall M, Mao L: Genome-wide analysis of the MADS-box gene family in Populus trichocarpa. Gene. 2006, 378: 84-94.PubMedView ArticleGoogle Scholar
- Becker A, Theissen G: The major clades of MADS-box genes and their role in the development and evolution of flowering plants. Mol Phylogenet Evol. 2003, 29 (3): 464-489. 10.1016/S1055-7903(03)00207-0.PubMedView ArticleGoogle Scholar
- Ferrandiz C, Gu Q, Martienssen R, Yanofsky MF: Redundant regulation of meristem identity and plant architecture by FRUITFULL, APETALA1 and CAULIFLOWER. Dev. 2000, 127 (4): 725-734.Google Scholar
- Ditta G, Pinyopich A, Robles P, Pelaz S, Yanofsky MF: The SEP4 gene of Arabidopsis thaliana functions in floral organ and meristem identity. Curr Biol. 2004, 14 (21): 1935-1940. 10.1016/j.cub.2004.10.028.PubMedView ArticleGoogle Scholar
- Kaufmann K, Wellmer F, Muino JM, Ferrier T, Wuest SE, Kumar V, Serrano-Mislata A, Madueno F, Krajewski P, Meyerowitz EM, et al: Orchestration of Floral Initiation by APETALA1. Sci. 2010, 328 (5974): 85-89. 10.1126/science.1185244.View ArticleGoogle Scholar
- Ohmori S, Kimizu M, Sugita M, Miyao A, Hirochika H, Uchida E, Nagato Y, Yoshida H: MOSAIC FLORAL ORGANS1, an AGL6-like MADS box gene, regulates floral organ identity and meristem fate in rice. Plant Cell. 2009, 21 (10): 3008-3025. 10.1105/tpc.109.068742.PubMedPubMed CentralView ArticleGoogle Scholar
- Rijpkema AS, Zethof J, Gerats T, Vandenbussche M: The petunia AGL6 gene has a SEPALLATA-like function in floral patterning. Plant J. 2009, 60 (1): 1-9. 10.1111/j.1365-313X.2009.03917.x.PubMedView ArticleGoogle Scholar
- Thompson BE, Bartling L, Whipple C, Hall DH, Sakai H, Schmidt R, Hake S: Bearded-ear encodes a MADS box transcription factor critical for maize floral development. Plant Cell. 2009, 21 (9): 2578-2590. 10.1105/tpc.109.067751.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, et al: Genome sequence of the palaeopolyploid soybean. Nature. 2010, 463 (7278): 178-183. 10.1038/nature08670.PubMedView ArticleGoogle Scholar
- Haerizadeh F, Wong CE, Singh MB, Bhalla PL: Genome-wide analysis of gene expression in soybean shoot apical meristem. Plant Mol Biol. 2009, 69 (6): 711-727. 10.1007/s11103-008-9450-1.PubMedView ArticleGoogle Scholar
- Liang D, Wong CE, Singh MB, Beveridge CA, Phipson B, Smyth GK, Bhalla PL: Molecular Dissection of Pea Shoot Apical Meristem. J Exp Bot. 2009, 60 (14): 4201-4213. 10.1093/jxb/erp254.PubMedPubMed CentralView ArticleGoogle Scholar
- Wong CE, Khor SY, Bhalla PL, Singh MB: Novel spatial expression of soybean WUSCHEL in the incipient floral primordial. Planta. 2011, 233: 553-560. 10.1007/s00425-010-1320-9.PubMedView ArticleGoogle Scholar
- Wong CE, Zhao Y-T, Wang X-J, Croft L, Wang Z-H, Haerizadeh F, Mattick JS, Singh MB, Carroll BJ, Bhalla PL: MicroRNAs in the shoot apical meristem of soybean. J Exp Bot. 2011, 62 (8): 2495-2506. 10.1093/jxb/erq437.PubMedView ArticleGoogle Scholar
- Hecht V, Foucher F, Ferrandiz C, Macknight R, Navarro C, Morin J, Vardy ME, Ellis N, Beltran JP, Rameau C, et al: Conservation of Arabidopsis flowering genes in model legumes. Plant Physiol. 2005, 137 (4): 1420-1434. 10.1104/pp.104.057018.PubMedPubMed CentralView ArticleGoogle Scholar
- Viaene T, Vekemans D, Becker A, Melzer S, Geuten K: Expression divergence of the AGL6 MADS domain transcription factor lineage after a core eudicot duplication suggests functional diversification. BMC Plant Biol. 2010, 10 (1): 148-10.1186/1471-2229-10-148.PubMedPubMed CentralView ArticleGoogle Scholar
- Libault M, Farmer A, Joshi T, Takahashi K, Langley RJ, Franklin LD, He J, Xu D, May G, Stacey G: An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. Plant J. 2010, 63 (1): 86-99.PubMedGoogle Scholar
- Tapia-Lopez R, Garcia-Ponce B, Dubrovsky JG, Garay-Arroyo A, Perez-Ruiz RV, Kim SH, Acevedo F, Pelaz S, Alvarez-Buylla ER: An AGAMOUS-related MADS-box gene, XAL1 (AGL12), regulates root meristem cell proliferation and flowering transition in Arabidopsis. Plant Physiol. 2008, 146 (3): 1182-1192. 10.1104/pp.107.108647.PubMedPubMed CentralView ArticleGoogle Scholar
- Rounsley SD, Ditta GS, Yanofsky MF: Diverse Roles for MADS Box Genes in Arabidopsis Development. Plant Cell. 1995, 7 (8): 1259-1269.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmid M, Uhlenhaut NH, Godard F, Demar M, Bressan R, Weigel D, Lohmann JU: Dissection of floral induction pathways using global expression analysis. Dev. 2003, 130 (24): 6001-6012. 10.1242/dev.00842.View ArticleGoogle Scholar
- Wong CE, Singh MB, Bhalla PL: Molecular processes underlying the floral transition in the soybean shoot apical meristem. Plant J. 2009, 57 (5): 832-845. 10.1111/j.1365-313X.2008.03730.x.PubMedView ArticleGoogle Scholar
- Koo SC, Bracko O, Park MS, Schwab R, Chun HJ, Park KM, Seo JS, Grbic V, Balasubramanian S, Schmid M, et al: Control of lateral organ development and flowering time by the Arabidopsis thaliana MADS-box Gene AGAMOUS-LIKE6. Plant J. 2010, 62 (5): 807-816. 10.1111/j.1365-313X.2010.04192.x.PubMedView ArticleGoogle Scholar
- Tucker SC: Floral Development in Legumes. Plant Physiol. 2003, 131 (3): 911-926. 10.1104/pp.102.017459.PubMedPubMed CentralView ArticleGoogle Scholar
- Tooke F, Ordidge M, Chiurugwi T, Battey N: Mechanisms and function of flower and inflorescence reversion. J Exp Bot. 2005, 56 (420): 2587-2599. 10.1093/jxb/eri254.PubMedView ArticleGoogle Scholar
- Washburn CF, Thomas JF: reversion of flowering in Glycine max (Fabaceae). Am J Bot. 2000, 87 (10): 1425-1438. 10.2307/2656869.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Wong CE, Singh MB, Bhalla PL: The dynamics of soybean leaf and shoot apical meristem transcriptome undergoing floral initiation process. PLoS One. 2013, 8 (6): e65319-10.1371/journal.pone.0065319.PubMedPubMed CentralView ArticleGoogle Scholar