Classification and phylogeny of the MADS-box multigene family suggest defined roles of MADS-box gene subfamilies in the morphological evolution of eukaryotes

Omics approaches to understand the MADS-box gene family in common bean (Phaseolus vulgaris L.) against drought stress

MADS-box genes are known to play important roles in diverse aspects of growth/devolopment and stress response in several plant species. However, no study has yet examined about MADS-box genes in P. vulgaris. In this study, a total of 79 PvMADS genes were identified and classified as type I and type II according to the phylogenetic analysis. While both type I and type II PvMADS classes were found to contain the MADS domain, the K domain was found to be present only in type II PvMADS proteins, in agreement with the literature. All chromosomes of the common bean were discovered to contain PvMADS genes and 17 paralogous gene pairs were identified. Only two of them were tandemly duplicated gene pairs (PvMADS-19/PvMADS-23 and PvMADS-20/PvMADS-24), and the remaining 15 paralogous gene pairs were segmentally duplicated genes. These duplications were found to play an important role in the expansion of type II PvMADS genes. Moreover, the RNAseq and RT-qPCR analyses showed the importance of PvMADS genes in response to drought stress in P. vulgaris.

Evolutionary analysis of MADS-box genes in buckwheat species and functional study of FdMADS28 in flavonoid metabolism

The MADS-box gene family is a transcription factor family that is widely expressed in plants. It controls secondary metabolic processes in plants and encourages the development of tissues like roots and flowers. However, the phylogenetic analysis and evolutionary model of MADS-box genes in Fagopyrum species has not been reported yet. This study identified the MADS-box genes of three buckwheat species at the whole genome level, and conducted systematic evolution and physicochemical analysis. The results showed that these genes can be divided into four subfamilies, with fragment duplication being the main way for the gene family expansion. During the domestication process from golden buckwheat to tartary buckwheat and the common buckwheat, the Ka/Ks ratio indicated that most members of the family experienced strong purification selection pressure, and with individual gene pairs experiencing positive selection. In addition, we combined the expression profile data of the MADS genes, mGWAS data, and WGCNA data to mine genes FdMADS28/48/50 that may be related to flavonoid metabolism. The results also showed that overexpression of FdMADS28 could increase rutin content by decreasing Kaempferol pathway content in hairy roots, and increase the resistance and growth of hairy roots to PEG and NaCl. This study systematically analyzed the evolutionary relationship of MADS-box genes in the buckwheat species, and elaborated on the expression patterns of MADS genes in different tissues under biotic and abiotic stresses, laying an important theoretical foundation for further elucidating their role in flavonoid metabolism.

MADS-Box Family Genes in Lagerstroemia indica and Their Involvement in Flower Development

MADS-box is a key transcription factor regulating the transition to flowering and flower development. Lagerstroemia indica 'Xiang Yun' is a new cultivar of crape myrtle characterized by its non-fruiting nature. To study the molecular mechanism underlying the non-fruiting characteristics of 'Xiang Yun', 82 MADS-box genes were identified from the genome of L. indica. The physicochemical properties of these genes were examined using bioinformatics methods, and their expression as well as endogenous hormone levels at various stages of flower development were analyzed. The results showed that LiMADS genes were primarily classified into two types: type I and type II, with the majority being type II that contained an abundance of cis-acting elements in their promoters. By screening nine core proteins by predicted protein interactions and performing qRT-PCR analysis as well as in combination with transcriptome data, we found that the expression levels of most MADS genes involved in flower development were significantly lower in 'Xiang Yun' than in the wild type 'Hong Ye'. Hormonal analysis indicated that 'Xiang Yun' had higher levels of iP, IPR, TZR, and zeatin during its early stages of flower development than 'Hong Ye', whereas the MeJA content was substantially lower at the late stage of flower development of 'Hong Ye'. Finally, correlation analysis showed that JA, IAA, SA, and TZR were positively correlated with the expression levels of most type II genes. Based on these analyses, a working model for the non-fruiting 'Xiang Yun' was proposed. During the course of flower development, plant hormone response pathways may affect the expression of MADS genes, resulting in their low expression in flower development, which led to the abnormal development of the stamen and embryo sac and ultimately affected the fruiting process of 'Xiang Yun'.

Two SEPALLATA MADS-Box Genes, SlMBP21 and SlMADS1, Have Cooperative Functions Required for Sepal Development in Tomato

MADS-box transcription factors have crucial functions in numerous physiological and biochemical processes during plant growth and development. Previous studies have reported that two MADS-box genes, SlMBP21 and SlMADS1, play important regulatory roles in the sepal development of tomato, respectively. However, the functional relationships between these two genes are still unknown. In order to investigate this, we simultaneously studied these two genes in tomato. Phylogenetic analysis showed that they were classified into the same branch of the SEPALLATA (SEP) clade. qRT-PCR displayed that both SlMBP21 and SlMADS1 transcripts are preferentially accumulated in sepals, and are increased with flower development. During sepal development, SlMBP21 is increased but SlMADS1 is decreased. Using the RNAi, tomato plants with reduced SlMBP21 mRNA generated enlarged and fused sepals, while simultaneous inhibition of SlMBP21 and SlMADS1 led to larger (longer and wider) and fused sepals than that in SlMBP21-RNAi lines. qRT-PCR results exhibited that the transcripts of genes relating to sepal development, ethylene, auxin and cell expansion were dramatically changed in SlMBP21-RNAi sepals, especially in SlMBP21-SlMADS1-RNAi sepals. Yeast two-hybrid assay displayed that SlMBP21 can interact with SlMBP21, SlAP2a, TAGL1 and RIN, and SlMADS1 can interact with SlAP2a and RIN, respectively. In conclusion, SlMBP21 and SlMADS1 cooperatively regulate sepal development in tomato by impacting the expression or activities of other related regulators or via interactions with other regulatory proteins.

Regulatory frameworks involved in the floral induction, formation and developmental programming of woody horticultural plants: a case study on blueberries

Flowering represents a crucial stage in the life cycles of plants. Ensuring strong and consistent flowering is vital for maintaining crop production amidst the challenges presented by climate change. In this review, we summarized key recent efforts aimed at unraveling the complexities of plant flowering through genetic, genomic, physiological, and biochemical studies in woody species, with a special focus on the genetic control of floral initiation and activation in woody horticultural species. Key topics covered in the review include major flowering pathway genes in deciduous woody plants, regulation of the phase transition from juvenile to adult stage, the roles of CONSTANS (CO) and CO-like gene and FLOWERING LOCUS T genes in flower induction, the floral regulatory role of GA-DELLA pathway, and the multifunctional roles of MADS-box genes in flowering and dormancy release triggered by chilling. Based on our own research work in blueberries, we highlighted the central roles played by two key flowering pathway genes, FLOWERING LOCUS T and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1, which regulate floral initiation and activation (dormancy release), respectively. Collectively, our survey shows both the conserved and diverse aspects of the flowering pathway in annual and woody plants, providing insights into the potential molecular mechanisms governing woody plants. This paves the way for enhancing the resilience and productivity of fruit-bearing crops in the face of changing climatic conditions, all through the perspective of genetic interventions.

Unregulated GmAGL82 due to Phosphorus Deficiency Positively Regulates Root Nodule Growth in Soybean

Nitrogen fixation, occurring through the symbiotic relationship between legumes and rhizobia in root nodules, is crucial in sustainable agriculture. Nodulation and soybean production are influenced by low levels of phosphorus stress. In this study, we discovered a MADS transcription factor, GmAGL82, which is preferentially expressed in nodules and displays significantly increased expression under conditions of phosphate (Pi) deficiency. The overexpression of GmAGL82 in composite transgenic plants resulted in an increased number of nodules, higher fresh weight, and enhanced soluble Pi concentration, which subsequently increased the nitrogen content, phosphorus content, and overall growth of soybean plants. Additionally, transcriptome analysis revealed that the overexpression of GmAGL82 significantly upregulated the expression of genes associated with nodule growth, such as GmENOD100, GmHSP17.1, GmHSP17.9, GmSPX5, and GmPIN9d. Based on these findings, we concluded that GmAGL82 likely participates in the phosphorus signaling pathway and positively regulates nodulation in soybeans. The findings of this research may lay the theoretical groundwork for further studies and candidate gene resources for the genetic improvement of nutrient-efficient soybean varieties in acidic soils.

Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.

Arabidopsis leaf-expressed AGAMOUS-LIKE 24 mRNA systemically specifies floral meristem differentiation

Plants can record external stimuli in mobile mRNAs and systemically deliver them to distal tissues to adjust development. Despite the identification of thousands of mobile mRNAs, the functional relevance of mobile mRNAs remains limited. Many mobile mRNAs are synthesized in the source cells that perceive environmental stimuli, but specifically exert their functions upon transportation to the recipient cells. However, the translation of mobile mRNA-encoded protein in the source cells could locally activate downstream target genes. How plants avoid ectopic functions of mobile mRNAs in the source cells to achieve tissue specificity remains to be elucidated. Here, we show that Arabidopsis AGAMOUS-LIKE 24 (AGL24) is a mobile mRNA whose movement is necessary and sufficient to specify floral organ identity. Although AGL24 mRNA is expressed in vegetative tissues, AGL24 protein exclusively accumulates in the shoot apex. In leaves, AGL24 proteins are degraded to avoid ectopically activating its downstream target genes. Our results reveal how selective protein degradation in source cells provides a strategy to limit the local effects associated with proteins encoded by mobile mRNAs, which ensures that mobile mRNAs specifically trigger systemic responses only in recipient tissues.

Jujube Witches' Broom Phytoplasma Effector Zaofeng3, a Homologous Effector of SAP54, Induces Abnormal Floral Organ Development and Shoot Proliferation

Plant-pathogenic phytoplasmas secrete specific virulence proteins into a host plant to modulate plant function for their own benefit. Identification of phytoplasmal effectors is a key step toward clarifying the pathogenic mechanisms of phytoplasma. In this study, Zaofeng3, also known as secreted jujube witches' broom phytoplasma protein 3 (SJP3), was a homologous effector of SAP54 and induced a variety of abnormal phenotypes, such as phyllody, malformed floral organs, witches' broom, and dwarfism in Arabidopsis thaliana. Zaofeng3 can also induce small leaves, dwarfism, and witches' broom in Ziziphus jujuba. Further experiments showed that the three complete α-helix domains predicted in Zaofeng3 were essential for induction of disease symptoms in jujube. Yeast two-hybrid library screening showed that Zaofeng3 mainly interacts with proteins involved in flower morphogenesis and shoot proliferation. Bimolecular fluorescence complementation assays confirmed that Zaofeng3 interacted with these proteins in the whole cell. Overexpression of zaofeng3 in jujube shoot significantly altered the expression patterns of ZjMADS19, ZjMADS47, ZjMADS48, ZjMADS77, and ZjTCP7, suggesting that overexpressing zaofeng3 might induce floral organ malformation and witches' broom by altering the expression of the transcriptional factors involved in jujube morphogenesis.

Genome-wide identification and analysis of MIKC-type MADS-box genes expression in Chimonanthus salicifolius

BACKGROUND

MIKC type MADS-box transcription factors are one of the largest gene families and play a pivotal role in flowering time and flower development. Chimonanthus salicifolius belongs to the family Calycanthaceae and has a unique flowering time and flowering morphology compared to other Chimonanthus species, but the research on MIKC type MADS-box gene family of C. salicifolius has not been reported.

OBJECTIVE

Identification, comprehensive bioinformatic analysis, the expression pattern of MIKC-type MADS-box gene family from different tissues of C. salicifolius.

METHODS

Genome-wide investigation and expression pattern under different tissues of the MIKC-type MADS-box gene family in C. salicifolius, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, promoter cis-acting element were performed.

RESULTS

A total of 29 MIKC-type MADS-box genes were identified from the whole genome sequencing. Interspecies synteny analysis revealed more significant collinearity between C. salicifolius and the magnoliids species compared to eudicots and monocots. MIKC-type MADS-box genes from the same subfamily share similar distribution patterns, gene structure, and expression patterns. Compared with Arabidopsis thaliana, Nymphaea colorata, and Chimonanthus praecox, the FLC genes were absent in C. salicifolius, while the AGL6 subfamily was expanded in C. salicifolius. The selectively expanded promoter (AGL6) and lack of repressor (FLC) genes may explain the earlier flowering in C. salicifolius. The loss of the AP3 homologous gene in C. salicifolius is probably the primary cause of the morphological distinction between C. salicifolius and C. praecox. The csAGL6a gene is specifically expressed in the flowering process and indicates the potential function of promoting flowering.

CONCLUSION

This study offers a genome-wide identification and expression profiling of the MIKC-types MADS-box genes in the C. salicifolius, and establishes the foundation for screening flowering development genes and understanding the potential function of the MIKC-types MADS-box genes in the C. salicifolius.

Updated Phylogeny and Protein Structure Predictions Revise the Hypothesis on the Origin of MADS-box Transcription Factors in Land Plants

MADS-box transcription factors (TFs), among the first TFs extensively studied, exhibit a wide distribution across eukaryotes and play diverse functional roles. Varying by domain architecture, MADS-box TFs in land plants are categorized into Type I (M-type) and Type II (MIKC-type). Type I and II genes have been considered orthologous to the SRF and MEF2 genes in animals, respectively, presumably originating from a duplication before the divergence of eukaryotes. Here, we exploited the increasing availability of eukaryotic MADS-box sequences and reassessed their evolution. While supporting the ancient duplication giving rise to SRF- and MEF2-types, we found that Type I and II genes originated from the MEF2-type genes through another duplication in the most recent common ancestor (MRCA) of land plants. Protein structures predicted by AlphaFold2 and OmegaFold support our phylogenetic analyses, with plant Type I and II TFs resembling the MEF2-type structure, rather than SRFs. We hypothesize that the ancestral SRF-type TFs were lost in the MRCA of Archaeplastida (the kingdom Plantae sensu lato). The retained MEF2-type TFs acquired a Keratin-like domain and became MIKC-type before the divergence of Streptophyta. Subsequently in the MRCA of land plants, M-type TFs evolved from a duplicated MIKC-type precursor through loss of the Keratin-like domain, leading to the Type I clade. Both Type I and II TFs expanded and functionally differentiated in concert with the increasing complexity of land plant body architecture. The recruitment of these originally stress-responsive TFs into developmental programs, including those underlying reproduction, may have facilitated the adaptation to the terrestrial environment.

Phylogenomic analysis provides insights into MADS-box and TCP gene diversification and floral development of the Asteraceae, supported by de novo genome and transcriptome sequences from dandelion (Taraxacum officinale)

The Asteraceae is the largest angiosperm family with more than 25,000 species. Individual studies have shown that MADS-box and TCP transcription factors are regulators of the development and symmetry of flowers, contributing to their iconic flower-head (capitulum) and floret. However, a systematic study of MADS-box and TCP genes across the Asteraceae is lacking. We performed a comparative analysis of genome sequences of 33 angiosperm species including our de novo assembly of diploid sexual dandelion (Taraxacum officinale) and 11 other Asteraceae to investigate the lineage-specific evolution of MADS-box and TCP genes in the Asteraceae. We compared the phylogenomic results of MADS-box and TCP genes with their expression in T. officinale floral tissues at different developmental stages to demonstrate the regulation of genes with Asteraceae-specific attributes. Here, we show that MADS-box MIKC^c and TCP-CYCLOIDEA (CYC) genes have expanded in the Asteraceae. The phylogenomic analysis identified AGAMOUS-like (AG-like: SEEDSTICK [STK]-like), SEPALATA-like (SEP3-like), and TCP-PROLIFERATING CELL FACTOR (PCF)-like copies with lineage-specific genomic contexts in the Asteraceae, Cichorioideae, or dandelion. Different expression patterns of some of these gene copies suggest functional divergence. We also confirm the presence and revisit the evolutionary history of previously named “Asteraceae-Specific MADS-box genes (AS-MADS).” Specifically, we identify non-Asteraceae homologs, indicating a more ancient origin of this gene clade. Syntenic relationships support that AS-MADS is paralogous to FLOWERING LOCUS C (FLC) as demonstrated by the shared ancient duplication of FLC and SEP3.

Expression Characterization of ABCDE Class MADS-Box Genes in Brassica rapa with Different Pistil Types

MADS-box is a vital transcription factor family that functions in plant growth and development. Apart from APETALA2, all genes in the ABCDE model that explain the molecular mechanism of floral organ development belong to the MADS-box family. Carpel and ovule numbers in plants are essential agronomic traits that determine seed yield, and multilocular siliques have great potential for the development of high-yield varieties of Brassica. In this study, ABCDE genes in the MADS-box family from Brassica rapa were identified and characterized. Their tissue-specific expression patterns in floral organs and their differential expression in different pistil types of B. rapa were revealed by qRT-PCR. A total of 26 ABCDE genes were found to belong to the MADS-box family. Our proposed ABCDE model of B. rapa is consistent with that of Arabidopsis thaliana, indicating that ABCDE genes are functionally conserved. These results of qRT-PCR showed that the expression levels of class C and D genes were significantly different between the wild-type (wt) and tetracarpel (tetrac) mutant of B. rapa. Interestingly, the expression of the homologs of class E genes was imbalanced. Therefore, it is speculated that class C, D, and E genes are involved in developing the carpel and ovule of B. rapa. Our findings reveal the potential for the selection of candidate genes to improve yield traits in Brassica crops.

Comparative transcriptomic analysis reveals differences in MADS-box genes of different hypericum in Changbai Mountains

To explore the differences between the hypericum in the Changbai Mountains, we carried out a transcriptome analysis of two common hypericums in the area, which was Hypericum attenuatum Choisy and Hypericum longistylum Oliv. We screened the MADS-box genes to analyze divergence time and evolutionary selection expression, and determine their expression levels. The results showed that we detected 9287 differentially expressed genes in the two species, of which shared 6044 genes by the two species. Analysis of the selected MADS genes revealed that the species was in an environment adapted to its natural evolution. The divergence time estimation showed that the segregation of these genes in the two species was related to the changes of external environment and genome replication events. The results of relative expression showed that the later flowering period of Hypericum attenuatum Choisy was related to the higher expression of the SVP (SHORT VEGETATIVE PHASE) and the AGL12 (AGAMOUS LIKE 12), while the lower expression of the FUL (FRUITFULL).

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Genome-wide analysis of the MADS-box gene family involved in salt and waterlogging tolerance in barley (Hordeum vulgare L.)

MADS-box transcription factors are crucial members of regulatory networks underlying multiple developmental pathways and abiotic stress regulatory networks in plants. Studies on stress resistance-related functions of MADS-box genes are very limited in barley. To gain insight into this gene family and elucidate their roles in salt and waterlogging stress resistance, we performed genome-wide identification, characterization and expression analysis of MADS-box genes in barley. A whole-genome survey of barley revealed 83 MADS-box genes, which were categorized into type I (Mα, Mβ and Mγ) and type II (AP1, SEP1, AGL12, STK, AGL16, SVP and MIKC*) lineages based on phylogeny, protein motif structure. Twenty conserved motifs were determined and each HvMADS contained one to six motifs. We also found tandem repeat duplication was the driven force for HvMADS gene family expansion. Additionally, the co-expression regulatory network of 10 and 14 HvMADS genes was predicted in response to salt and waterlogging stress, and we proposed HvMADS11,13 and 35 as candidate genes for further exploration of the functions in abiotic stress. The extensive annotations and transcriptome profiling reported in this study ultimately provides the basis for MADS functional characterization in genetic engineering of barley and other gramineous crops.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)₆GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKC^C-type MTFs (M^C-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. M^C-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which M^C-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of M^C-MTFs). A special focus of our review is on the structural prerequisites of M^C-MTFs that enable the specific recognition of target genes.

Independent recruitment of FRUITFULL-like transcription factors in the convergent origins of vernalization-responsive grass flowering

Flowering in response to low temperatures (vernalization) has evolved multiple times independently across angiosperms as an adaptation to match reproductive development with the short growing season of temperate habitats. Despite the context of a generally conserved flowering time network, evidence suggests that the genes underlying vernalization responsiveness are distinct across major plant clades. Whether different or similar mechanisms underlie vernalization-induced flowering at narrower (e.g., family-level) phylogenetic scales is not well understood. To test the hypothesis that vernalization responsiveness has evolved convergently in temperate species of the grass family (Poaceae), we carried out flowering time experiments with and without vernalization in several representative species from different subfamilies. We then determined the likelihood that vernalization responsiveness evolved through parallel mechanisms by quantifying the response of Pooideae vernalization pathway FRUITFULL (FUL)-like genes to extended periods of cold. Our results demonstrate that vernalization-induced flowering has evolved multiple times independently in at least five grass subfamilies, and that different combinations of FUL-like genes have been recruited to this pathway on several occasions.

Cone-setting in spruce is regulated by conserved elements of the age-dependent flowering pathway

Reproductive phase change is well characterized in angiosperm model species, but less studied in gymnosperms. We utilize the early cone-setting acrocona mutant to study reproductive phase change in the conifer Picea abies (Norway spruce), a gymnosperm. The acrocona mutant frequently initiates cone-like structures, called transition shoots, in positions where wild-type P. abies always produces vegetative shoots. We collect acrocona and wild-type samples, and RNA-sequence their messenger RNA (mRNA) and microRNA (miRNA) fractions. We establish gene expression patterns and then use allele-specific transcript assembly to identify mutations in acrocona. We genotype a segregating population of inbred acrocona trees. A member of the SQUAMOSA BINDING PROTEIN-LIKE (SPL) gene family, PaSPL1, is active in reproductive meristems, whereas two putative negative regulators of PaSPL1, miRNA156 and the conifer specific miRNA529, are upregulated in vegetative and transition shoot meristems. We identify a mutation in a putative miRNA156/529 binding site of the acrocona PaSPL1 allele and show that the mutation renders the acrocona allele tolerant to these miRNAs. We show co-segregation between the early cone-setting phenotype and trees homozygous for the acrocona mutation. In conclusion, we demonstrate evolutionary conservation of the age-dependent flowering pathway and involvement of this pathway in regulating reproductive phase change in the conifer P. abies.

Genome-Wide Identification, Evolution, and Expression Characterization of the Pepper (Capsicum spp.) MADS-box Gene Family

MADS domain transcription factors play roles throughout the whole lifecycle of plants from seeding to flowering and fruit-bearing. However, systematic research into MADS-box genes of the economically important vegetable crop pepper (Capsicum spp.) is still lacking. We identified 174, 207, and 72 MADS-box genes from the genomes of C. annuum, C. baccatum, and C. chinense, respectively. These 453 MADS-box genes were divided into type I (Mα, Mβ, Mγ) and type II (MIKC* and MIKC^C) based on their phylogenetic relationships. Collinearity analysis identified 144 paralogous genes and 195 orthologous genes in the three Capsicum species, and 70, 114, and 10 MADS-box genes specific to C. annuum, C. baccatum, and C. chinense, respectively. Comparative genomic analysis highlighted functional differentiation among homologous MADS-box genes during pepper evolution. Tissue expression analysis revealed three main expression patterns: highly expressed in roots, stems, leaves, and flowers (CaMADS93/CbMADS35/CcMADS58); only expressed in roots; and specifically expressed in flowers (CaMADS26/CbMADS31/CcMADS11). Protein interaction network analysis showed that type II CaMADS mainly interacted with proteins related to flowering pathway and flower organ development. This study provides the basis for an in-depth study of the evolutionary features and biological functions of pepper MADS-box genes.

Genome-wide identification and expression analysis of MIKCC genes in rose provide insight into their effects on flower development

The MIKC^C-type gene family plays important roles in plant growth, development, and tolerance of biotic and abiotic stress, especially during floral organ differentiation. However, there have been no studies of MIKC^C-type genes in rose, and functional differentiation of family members has not been explored. In this study, we identified 42 MIKC^C-type genes in rose, classified the genes into 12 subfamilies, and constructed a phylogenetic tree. We performed expression analysis of these genes, and found that expression patterns correlated with the predicted subfamily, indicating that the features of MIKC^C-type genes were broadly conserved during evolution. Collinear analysis of MIKC^C genes among Rosaceae species confirmed the occurrence of whole genome duplications (WGD) and revealed some species-specific MIKC^C genes. Transcriptome analysis showed that the expression of some MIKC^C-type genes responded to low temperatures (4°C, 24 h) during flower organ differentiation. These conserved, duplicated, and novel expression patterns of MIKC^C-type genes may have facilitated the adaptation of rose to various internal and external environmental changes. The results of this study provide a theoretical basis for future functional analysis of the MIKC^C genes in rose and investigation of the evolutionary pattern of the MIKC^C gene family in the Rosaceae genome.

Adaptive evolvability through direct selection instead of indirect, second-order selection

Can evolvability itself be the product of adaptive evolution? To answer this question is challenging, because any DNA mutation that alters only evolvability is subject to indirect, "second order" selection on the future effects of this mutation. Such indirect selection is weaker than "first-order" selection on mutations that alter fitness, in the sense that it can operate only under restrictive conditions. Here I discuss a route to adaptive evolvability that overcomes this challenge. Specifically, a recent evolution experiment showed that some mutations can enhance both fitness and evolvability through a combination of direct and indirect selection. Unrelated evidence from gene duplication and the evolution of gene regulation suggests that mutations with such dual effects may not be rare. Through such mutations, evolvability may increase at least in part because it provides an adaptive advantage. These observations suggest a research program on the adaptive evolution of evolvability, which aims to identify such mutations and to disentangle their direct fitness effects from their indirect effects on evolvability. If evolvability is itself adaptive, Darwinian evolution may have created more than life's diversity. It may also have helped create the very conditions that made the success of Darwinian evolution possible.

BSISTER transcription factors directly binds to the promoter of IAA19 and IAA29 genes to up-regulate gene expression and promote the root development

Roots play an important role in the growth and development of plants and auxin participates in regulating plant root development. Some studies have shown that BS (BSISTER) gene (the closest gene of class B gene) is involved in plant root development, but whether BS regulates root development via auxin signaling still not clear. To explore VviBS1 and VviBS2 roles in root development, VviBS1 and VviBS2 were overexpressedin Arabidopsis tt16 mutant and we found that they could restore the phenotype of shorter PR (primary roots) and high density of LR (lateral root) of tt16 compared with the wild type Ws Arabidopsis seedlings. However, the addition of exogenous NAA (naphthalene acetic acid) could not significantly promote the PR length of tt16 Arabidopsis, and the auxin signal transduction of tt16 may be blocked. The expression levels of auxin signal transduction pathway genes in Ws, tt16, p35s:VviBS1 in tt16 and p35s:VviBS2 in tt16 seedlings were detected. It was found that the expression of AtARF2, AtARF12, AtARF14, AtARF15, AtARF20, AtGH3, AtGH3-2 and AtSAUR51 genes in tt16 seedlings was higher than that in Ws, while the expression of AtIAA19 and AtIAA29 in Ws seedlings was higher than that of tt16. More importantly, BS may up regulate AtIAA19 and AtIAA29 expression directly by binding to their promoter. In addition, VviBS1 and VviBS2 also affect seed germination and may regulate leaf yellowing by regulating ethylene synthase. Therefore, our findings reveal a molecular mechanism that BS may modulate root system development via Aux/IAA-based auxin signaling, and provide insight into the BS function in regulation of leaf yellowing.

Functional divergence of GhAP1.1 and GhFUL2 associated with flowering regulation in upland cotton (Gossypium hirsutum L.)

The AP1/FUL transcription factors are important for floral development, but the underlying molecular mechanisms remain unclear. In this study, we cloned and identified two AP1/FUL-like genes, GhAP1.1 and GhFUL2, in upland cotton, which is a commonly cultivated economically valuable crop. Sequence alignment and phylogenetic analysis indicated that GhAP1.1 and GhFUL2, which are encoded by genes in the AP1/FUL clade, have conserved N-terminal regions but diverse C-terminal domains. Quantitative real-time PCR analysis revealed that GhAP1.1 and GhFUL2 were expressed in the flower and root, and showed opposite expression patterns during shoot apical meristem development. The upregulated expression of GhAP1.1 in Arabidopsis did not result in significant changes to the flowering time or floral organ development, and the transcript levels of the florigen FT increased and those of LFY decreased. Overexpression of GhFUL2 in Arabidopsis delayed flowering and promoted bolting by decreasing FT and LFY transcript levels. Silencing GhFUL2 in cotton dramatically increased the expression of GhFT and GhAP1.3 and promoted flowering. Yeast two-hybrid and bimolecular fluorescence complementation assays indicated that GhAP1.1 could interact with the SVP homolog GhSVP2.2, whereas GhFUL2 formed heterodimers with GhSEP3/GhSEP4 homologs and GhSVP2.2. The present results demonstrated that the functional divergence of GhAP1.1 and GhFUL2, which involved changes in sequences and expression patterns, influenced the regulation of cotton flower development.

Testing candidate genes linked to corolla shape variation of a pollinator shift in Rhytidophyllum (Gesneriaceae)

Floral adaptations to specific pollinators like corolla shape variation often result in reproductive isolation and thus speciation. But despite their ecological importance, the genetic bases of corolla shape transitions are still poorly understood, especially outside model species. Hence, our goal was to identify candidate genes potentially involved in corolla shape variation between two closely related species of the Rhytidophyllum genus (Gesneriaceae family) from the Antilles with contrasting pollination strategies. Rhytidophyllum rupincola has a tubular corolla and is strictly pollinated by hummingbirds, whereas R. auriculatum has more open flowers and is pollinated by hummingbirds, bats, and insects. We surveyed the literature and used a comparative transcriptome sequence analysis of synonymous and non-synonymous nucleotide substitutions to obtain a list of genes that could explain floral variation between R. auriculatum and R. rupincola. We then tested their association with corolla shape variation using QTL mapping in a F₂ hybrid population. Out of 28 genes tested, three were found to be good candidates because of a strong association with corolla shape: RADIALIS, GLOBOSA, and JAGGED. Although the role of these genes in Rhytidophyllum corolla shape variation remains to be confirmed, these findings are a first step towards identifying the genes that have been under selection by pollinators and thus involved in reproductive isolation and speciation in this genus.

Genome-Wide Analysis of the MADS-Box Gene Family in Holoparasitic Plants (Balanophora subcupularis and Balanophora fungosa var. globosa)

MADS-box is an important transcription factor family that is involved in the regulation of various stages of plant growth and development, especially flowering regulation and flower development. Being a holoparasitic plant, the body structure of Balanophoraceae has changed dramatically over time, and its vegetative and reproductive organs have been extensively modified, with rudimentary flower organs. Meanwhile, extraordinary gene losses have been identified in holoparasitic plants compared with autotrophs. Our study reveals that the MADS-box gene family contracted sharply in Balanophora subcupularis and Balanophora fungosa var. globosa, and some subfamilies were lost, exhibiting reduced redundancy in both. The genes that functioned in the transition from the vegetative to floral production stages suffered a significant loss, but the ABCE model genes remained intact. We further investigated genes related to flowering regulation in B. subcupularis and B. fungosa var. globosa, vernalization and autonomous ways of regulating flowering time remained comparatively integrated, while genes in photoperiod and circadian clock pathways were almost lost. Convergent gene loss in flowering regulation occurred in Balanophora and another holoparasitic plant Sapria himalayana (Rafflesiaceae). The genome-wide analysis of the MADS-box gene family in Balanophora species provides valuable information for understanding the classification, gene loss pattern, and flowering regulation mechanism of MADS-box gene family in parasitic plants.

PIN3 positively regulates the late initiation of ovule primordia in Arabidopsis thaliana

Ovule initiation determines the maximum ovule number and has great impact on seed number and yield. However, the regulation of ovule initiation remains largely elusive. We previously reported that most of the ovule primordia initiate asynchronously at floral stage 9 and PINFORMED1 (PIN1) polarization and auxin distribution contributed to this process. Here, we further demonstrate that a small amount of ovule primordia initiate at floral stage 10 when the existing ovules initiated at floral stage 9 start to differentiate. Genetic analysis revealed that the absence of PIN3 function leads to the reduction in pistil size and the lack of late-initiated ovules, suggesting PIN3 promotes the late ovule initiation process and pistil growth. Physiological analysis illustrated that, unlike picloram, exogenous application of NAA can’t restore these defective phenotypes, implying that PIN3-mediated polar auxin transport is required for the late ovule initiation and pistil length. qRT-PCR results indicated that the expression of SEEDSTICK (STK) is up-regulated under auxin analogues treatment while is down-regulated in pin3 mutants. Meanwhile, overexpressing STK rescues pin3 phenotypes, suggesting STK participates in PIN3-mediated late ovule initiation possibly by promoting pistil growth. Furthermore, brassinosteroid influences the late ovule initiation through positively regulating PIN3 expression. Collectively, this study demonstrates that PIN3 promotes the late ovule initiation and contributes to the extra ovule number. Our results give important clues for increasing seed number and yield of cruciferous and leguminous crops.

Comparative Analysis of the MADS-Box Genes Revealed Their Potential Functions for Flower and Fruit Development in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is an important economic crop widely planted in tropical and subtropical regions, and flower and fruit development play decisive effects on the longan yield and fruit quality formation. MCM1, AGAMOUS, DEFICIENS, Serum Response Factor (MADS)-box transcription factor family plays important roles for the flowering time, floral organ identity, and fruit development in plants. However, there is no systematic information of MADS-box family in longan. In this study, 114 MADS-box genes were identified from the longan genome, phylogenetic analysis divided them into type I (Mα, Mβ, Mγ) and type II (MIKC*, MIKC ^C ) groups, and MIKC ^C genes were further clustered into 12 subfamilies. Comparative genomic analysis of 12 representative plant species revealed the conservation of type II in Sapindaceae and analysis of cis-elements revealed that Dof transcription factors might directly regulate the MIKC ^C genes. An ABCDE model was proposed for longan based on the phylogenetic analysis and expression patterns of MADS-box genes. Transcriptome analysis revealed that MIKC ^C genes showed wide expression spectrums, particularly in reproductive organs. From 35 days after KClO₃ treatment, 11 MIKC genes were up-regulated, suggesting a crucial role in off-season flower induction, while DlFLC, DlSOC1, DlSVP, and DlSVP-LIKE may act as the inhibitors. The gene expression patterns of longan fruit development indicated that DlSTK, DlSEP1/2, and DlMADS53 could be involved in fruit growth and ripening. This paper carried out the whole genome identification and analysis of the longan MADS-box family for the first time, which provides new insights for further understanding its function in flowers and fruit.

Circlize package in R and Analytic Hierarchy Process (AHP): Contribution values of ABCDE and AGL6 genes in the context of floral organ development

The morphological diversity of floral organs can largely be attributed to functional divergence in the MADS-box gene family. Nonetheless, research based on the ABCDE model has yet to conclusively determine whether the AGAMOUS-LIKE 6 (AGL6) subgroup has a direct influence on floral organ development. In the current study, the ABCDE model was used to quantify the contributions of ABCDE and AGL6 genes in the emergence of floral organs. We determined that the flower formation contribution values of the ABCDE and AGL6 genes were as follows: A gene, 0.192; B gene, 0.231; CD gene, 0.192; E gene, 0.385; and AGL6, 0.077. As AGL6 does not directly influence floral structure formation, the contribution value of AGL6 to flower formation was low. Furthermore, the gradient values of the floral organs were as follows: sepals, 0.572; petals, 1.606; stamens, 2.409; and carpels, 2.288. We also performed detailed analysis of the ABCDE and AGL6 genes using the Circlize package in R. Our results suggest that these genes likely emerged in one of two orders: 1) B genes→CD genes→AGL6→E genes→A genes; or 2) B genes→CD genes→AGL6/E genes→A genes. We use the analytic hierarchy process (AHP) to prove the contribution values and gradient values of floral organs. This is the first study to understand the contribution values of ABCDE and AGL6 genes using the AHP and the Circlize package in R.

Genome-Wide Analysis and the Expression Pattern of the MADS-Box Gene Family in Bletilla striata

Bletilla striata (Thunb. ex A. Murray) Rchb. f., a species of the perennial herb Orchidaceae, has potent anti-inflammatory and antiviral biological activities. MADS-box transcription factors play critical roles in the various developmental processes of plants. Although this gene family has been extensively investigated in many species, it has not been analyzed for B. striata. In total, 45 MADS-box genes were identified from B. striata in this study, which were classified into five subfamilies (Mδ, MIKC, Mα, Mβ, and Mγ). Meanwhile, the highly correlated protein domains, motif compositions, and exon-intron structures of BsMADSs were investigated according to local B. striata databases. Chromosome distribution and synteny analyses revealed that segmental duplication and homologous exchange were the main BsMADSs expansion mechanisms. Further, RT-qPCR analysis revealed that BsMADSs had different expression patterns in response to various stress treatments. Our results provide a potential theoretical basis for further investigation of the functions of MADS genes during the growth of B. striata.

Utilizing MIKC-type MADS-box protein SOC1 for yield potential enhancement in maize

Overexpression of Zea mays SOC gene promotes flowering, reduces plant height, and leads to no reduction in grain production per plant, suggesting enhanced yield potential, at least, through increasing planting density. MIKC-type MADS-box gene SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) is an integrator conserved in the plant flowering pathway. In this study, the maize SOC1 (ZmSOC1) gene was cloned and overexpressed in transgenic maize Hi-II genotype. The T₀ plants were backcrossed with nontransgenic inbred B73 to produce first generation backcross (BC₁) seeds. Phenotyping of both transgenic and null segregant (NT) BC₁ plants was conducted in three independent experiments. The BC₁ transgenic plants showed new attributes such as increased vegetative growth, accelerated flowering time, reduced overall plant height, and increased grain weight. Second generation backcross (BC₂) plants were evaluated in the field using two planting densities. Compared to BC₂ NT plants, BC₂ transgenic plants, were 12-18% shorter, flowered 5 days earlier, and showed no reduction in grain production per plant and an increase in fat, starch, and simple sugars in the grain. Transcriptome comparison in young leaves of 56-day-old BC₁ plants revealed that the overexpressed ZmSOC1 resulted in 107 differentially expressed genes. The upregulated transcription factor DNA BINDING WITH ONE FINGER 5.4 (DOF5.4) was among the genes responsible for the reduced plant height. Modulating expression of SOC1 opens a new and effective approach to promote flowering and reduce plant height, which may have potential to enhance crop yield and improve grain quality.

Genome-wide identification of the MADS-box transcription factor family in autotetraploid cultivated alfalfa (Medicago sativa L.) and expression analysis under abiotic stress

Background

Alfalfa, the “queen of forage”, is the most extensively cultivated forage legume in the world. The development and yield of alfalfa are seriously limited by abiotic stress. MADS-box transcription factors are one of the largest gene families and play a pivotal role in plant development and abiotic stress. However, little is known regarding the MADS-box transcription factors in autotetraploid cultivated alfalfa.

Results

In the present study, we identified 120 MsMADS-box genes in the alfalfa genome. Phylogenetic analysis indicated that 75 type-I MsMADS-box genes were classified into the Mα, Mβ, and Mγ subgroups, and 45 type-II MsMADS-box genes were classified into 11 subgroups. The promoter region of MsMADS-box genes containing several hormone and stress related elements. Chromosomal location analysis revealed that 117 MsMADS-box genes were unevenly distributed on 32 chromosomes, and the remaining three genes were located on unmapped scaffolds. A total of nine pairs of segmental duplications and four groups of tandem duplications were found. Expression analysis showed that MsMADS-box genes were differentially expressed in various tissues and under abiotic stresses. qRT-PCR analysis revealed that the expression profiles of eight selected MsMADS-box genes were distinct under various stresses.

Conclusions

In this study, MsMADS-box genes were identified in the cultivated alfalfa genome based on autotetraploid level, and further confirmed by Gene Ontology (GO) analysis, phylogenetic analysis, sequence features and expression analysis. Taken together, these findings will provide clues for further study of MsMADS-box functions and alfalfa molecular breeding.

Our study is the first to systematically identify and characterize the MADS-box transcription factors in autotetraploid cultivated alfalfa (Medicago sativa L.), and eight MsMADS-box genes were significantly involved in response to various stresses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07911-9.

Evolutionary divergence of motifs in B-class MADS-box proteins of seed plants

Background

MADS-box transcription factors function as homo- or heterodimers and regulate many aspects of plant development; moreover, MADS-box genes have undergone extensive duplication and divergence. For example, the morphological diversity of floral organs is closely related to the functional divergence of the MADS-box gene family. B-class genes (such as Arabidopsis thaliana APETALA3 [AP3] and PISTILLATA [PI]) belong to a subgroup of MADS-box genes. Here, we collected 97 MADS-box B protein sequences from 21 seed plant species and examined their motifs to better understand the functional evolution of B proteins.

Results

We used the MEME tool to identify conserved sequence motifs in these B proteins; unique motif arrangements and sequences were identified in these B proteins. The keratin-like domains of Malus domestica and Populus trichocarpa B proteins differed from those in other angiosperms, suggesting that a novel regulatory network might have evolved in these species. The MADS domains of Nelumbo nucifera, Glycine max, and Amborella trichopoda B-proteins contained motif 9; in contrast, those of other plants contained motif 1. Protein modelling analyses revealed that MADS domains with motif 9 may lack amino acid sites required for DNA-binding. These results suggested that the three species might share an alternative mechanism controlling floral development.

Conclusions

Amborella trichopoda has B proteins with either motif 1 or motif 9 MADS domains, suggesting that these two types of MADS domains evolved from the ancestral domain into two groups, those with motif 9 (N. nucifera and G. max), and those with motif 1. Moreover, our results suggest that the homodimer/heterodimer intermediate transition structure first appeared in A. trichopoda. Therefore, our systematic analysis of the motifs in B proteins sheds light on the evolution of these important transcription factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40709-021-00144-7.

K-Domain Technology: Constitutive Expression of a Blueberry Keratin-Like Domain Mimics Expression of Multiple MADS-Box Genes in Enhancing Maize Grain Yield

MADS-box genes are considered as the foundation of all agronomic traits because they play essential roles in almost every aspect of plant reproductive development. Keratin-like (K) domain is a conserved protein domain of tens of MIKC-type MADS-box genes in plants. K-domain technology constitutively expresses a K-domain to mimic expression of the K-domains of other MADS-box genes simultaneously and thus to generate new opportunities for yield enhancement, because the increased K-domains can likely prevent MADS-domain proteins from binding to target DNA. In this study, we evaluated utilizing the K-domain technology to increase maize yield. The K-domain of a blueberry's SUPPRESSOR of CONSTITUTIVE EXPRESSION OF CONSTANS 1 (VcSOC1K) has similarities to five MADS-box genes in maize. Transgenic maize plants expressing the VcSOC1K showed 13-100% of more grain per plant than the nontransgenic plants in all five experiments conducted under different experimental conditions. Transcriptome comparisons revealed 982 differentially expressed genes (DEGs) in the leaves from 83-day old plants, supporting that the K-domain technology were powerful and multiple functional. The results demonstrated that constitutive expression of the VcSOC1K was very effective to enhance maize grain production. With the potential of mimicking the K-domains of multiple MADS-box genes, the K-domain technology opens a new approach to increase crop yield.

Identification and characterization of MADS box gene family in pigeonpea for their role during floral transition

MADS box genes are class of transcription factors involved in various physiological and developmental processes in plants. To understand their role in floral transition-related pathways, a genome-wide identification was done in Cajanus cajan, identifying 102 members which were classified into two different groups based on their gene structure. The status of all these genes was further analyzed in three wild species i.e. C. scarabaeoides, C. platycarpus and C. cajanifolius which revealed absence of 31-34 MADS box genes in them hinting towards their role in domestication and evolution. We could locate only a single copy of both FLOWERING LOCUS C (FLC) and SHORT VEGETATIVE PHASE (SVP) genes, while three paralogs of SUPPRESSOR OF ACTIVATION OF CONSTANS 1 (SOC1) were found in C. cajan genome. One of those SOC1 paralogs i.e. CcMADS1.5 was found to be missing in all three wild relatives, also forming separate clade in phylogeny. This SOC1 gene was also lacking the characteristic MADS box domain in it. Expression profiling of major MADS box genes involved in flowering was done in different tissues viz shoot apical meristem, vegetative leaf, reproductive meristem, and reproductive bud. Gene-based time tree of FLC and SOC1 gene dictates their divergence from Arabidopsis before 71 and 23 million year ago (mya), respectively. This study provides valuable insights into the functional characteristics, expression pattern, and evolution of MADS box proteins in grain legumes with emphasis on C. cajan, which may help in further characterizing these genes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02605-7.

Orchid Bsister gene PeMADS28 displays conserved function in ovule integument development

The ovules and egg cells are well developed to be fertilized at anthesis in many flowering plants. However, ovule development is triggered by pollination in most orchids. In this study, we characterized the function of a B_sister gene, named PeMADS28, isolated from Phalaenopsis equestris, the genome-sequenced orchid. Spatial and temporal expression analysis showed PeMADS28 predominantly expressed in ovules between 32 and 48 days after pollination, which synchronizes with integument development. Subcellular localization and protein–protein interaction analyses revealed that PeMADS28 could form a homodimer as well as heterodimers with D-class and E-class MADS-box proteins. In addition, ectopic expression of PeMADS28 in Arabidopsis thaliana induced small curled rosette leaves, short silique length and few seeds, similar to that with overexpression of other species’ B_sister genes in Arabidopsis. Furthermore, complementation test revealed that PeMADS28 could rescue the phenotype of the ABS/TT16 mutant. Together, these results indicate the conserved function of B_sisterPeMADS28 associated with ovule integument development in orchid.

Asynchrony of ovule primordia initiation in Arabidopsis

Plant ovule initiation determines the maximum of ovule number and has a great impact on the seed number per fruit. The detailed processes of ovule initiation have not been accurately described, although two connected processes, gynoecium and ovule development, have been investigated. Here, we report that ovules initiate asynchronously. The first group of ovule primordia grows out, the placenta elongates, the boundaries of existing ovules enlarge and a new group of primordia initiates from the boundaries. The expression pattern of different marker genes during ovule development illustrates that this asynchronicity continues throughout whole ovule development. PIN-FORMED1 polar distribution and auxin response maxima correlate with ovule primordia asynchronous initiation. We have established computational modeling to show how auxin dynamics influence ovule primordia initiation. Brassinosteroid signaling positively regulates ovule number by promoting placentae size and ovule primordia initiation through strengthening auxin response. Transcriptomic analysis demonstrates numerous known regulators of ovule development and hormone signaling, and many new genes are identified that are involved in ovule development. Taken together, our results illustrate that the ovule primordia initiate asynchronously and the hormone signals are involved in the asynchrony.

Summary: Ovule primordia initiation, which determines the maximum ovule number and subsequent seed number in Arabidopsis, is asynchronous and is regulated by PIN1 polar distribution and the auxin response.

Genome-wide study of flowering-related MADS-box genes family in Cardamine hirsuta

MADS-box genes take part in diverse biological functions especially in development of reproductive structures and control of flowering time. Recently, Cardamine hirsuta has emerged as an exclusively powerful genetic system in comparative studies of development. Although the C. hirsuta genome sequence is available but a comprehensive analysis of its MADS-box family genes is still lacking. Here, we determined 50 Cardamine MADS-box genes through bioinformatics tools and classified them into 2 Mβ, 6 Mα and 2 Mγ and 40 MIKC-type (35 MIKCc and 5MIKC*) genes based on a phylogenetic analysis. The C. hirsuta MIKC subfamily could be further classified into 14 subgroups as Arabidopsis. However the number of MADS-box proteins was not equal among these subgroups. Based on the structural diversity among 50 MADS-box genes, 2 lineages were obtained, type I and type II. The lowest number of introns (0 or 1) was found in the Mα, Mβ, and Mγ groups of the type I genes. The most Cardamine MADS-box genes were randomly distributed on only three chromosomes. C. hirsuta had a relatively lower number of flowering MADS-box genes than A. thaliana and probably tandem duplication event resulted in the expansion of FLC, SQUA and TM3 family members in Arabidopsis. Moreover among the conserved motifs, ChMADS5 of SQUA, ChMADS34 of TM3 and ChMADS51 of AGL15 families had no K-domain. This study provides a basis for further functional investigation of MADS-box genes in C. hirsuta.

Molecular cloning and potential role of DiSOC1s in flowering regulation in Davidia involucrata Baill

Davidia involucrata Baill. (dove tree) is unique Tertiary relic plant in China, also known as 'living fossil' and 'giant panda'. The MADS-box family gene SOC1 is involved in the regulatory pathway that integrates flowering signals to promote flowering at the optimal time. In this study, we isolated and identified two dove tree SOC1 homologues, named DiSOC1-a and DiSOC1-b. These two sequences possess highly conserved domains MADS-box and SOC1-motif, as well as the semi-conserved region K-box. DiSOC1-a and DiSOC1-b were expressed at varying levels in all tested tissues of dove tree and shared high levels of expression in the flower buds. The expression tendencies of both genes in bract were initially upward and then downward and were highest in young bracts. Neither DiSOC1-a nor DiSOC1-b was expressed in immature leaves. Proteins encoded by DiSOC1-a and DiSOC1-b were located in the nucleus. In addition, ectopic overexpression of both genes in WT Arabidopsis promoted early flowering and the growth of the main bolt. Taken together, these results suggest that DiSOC1-a and DiSOC1-b are involved in the flowering initiation and the main bolt growth process of dove tree. Our results provide a foundation for horticultural breeding to control flowering time of dove tree.

The MADS transcription factor GhAP1.7 coordinates the flowering regulatory pathway in upland cotton (Gossypium hirsutum L.)

MADS-box gene family plays an important role in the molecular regulatory network of flower development. APETALA1 (AP1), a MADS-box gene, plays an important role in the development of flower organs. Although many studies about MADS-box family genes have been reported, the function of AP1 is still not clear in cotton. In this study, GhAP1.7 (Gh_D03G0922), a candidate gene for cotton flower time and plant height obtained from our previous studies, was cloned from CCRI50 cotton variety and functionally characterized. Subcellular localization demonstrated that GhAP1.7 was located in nucleus. Infection test of Arabidopsis revealed that GhAP1.7 could cause precocious flowering and virus-induced gene silence (VIGS) assay demonstrated that GhAP1.7 could lead to delayed flowering of cotton plants. Yeast one-hybrid assays and transient dual-luciferase assays suggested that floral meristem identity control gene LEAFY (LFY) can bind the promoter of GhAP1.7 and negatively regulate it. Our research indicated that GhAP1.7 might work as a positive regulator in plant flowering. Moreover, GhAP1.7 may negatively regulated by GhLFY in the regulatory pathways. This work laid the foundation for subsequent functional studies of GhAP1.7.

MawuAP1 promotes flowering and fruit development in the basal angiosperm Magnolia wufengensis (Magnoliaceae)

The APETALA1/SQUAMOSA (AP1/SQUA)-like genes of flowering plants play crucial roles in the development processes of floral meristems, sepals, petals and fruits. Although many of the AP1/SQUA-like genes have been characterized in angiosperms, few have been identified in basal angiosperm taxa. Therefore, the functional evolution of the AP1/SQUA subfamily is still unclear. We characterized an AP1 homolog, MawuAP1, from Magnolia wufengensis that is an ornamental woody plant belonging to the basal angiosperms. Gene sequence and phylogenetic analyses suggested that MawuAP1 was clustered with the FUL-like homologous genes of basal angiosperms and had FUL motif and paleoAP1 motif domain, but it did not have the euAP1 motif domain of core eudicots. Expression pattern analysis showed that MawuAP1 was highly expressed in vegetative and floral organs, particularly in the early stage of flower bud development and pre-anthesis. Protein-protein interaction pattern analysis revealed that MawuAP1 has interaction with an A-class gene (MawuAP1), C-class gene (MawuAG-1) and E-class gene (MawuAGL9) of the MADS-box family genes. Ectopic expression in Arabidopsis thaliana indicated that MawuAP1 could significantly promote flowering and fruit development, but it could not restore the sepal and petal formation of ap1 mutants. These results demonstrated that there are functional differences in the specification of sepal and petal floral organs and development of fruits among the AP1/SQUA-like genes, and functional conservation in the regulation of floral meristem. These findings provide strong evidence for the important functions of MawuAP1 in floral meristem determination, promoting flowering and fruit development, and further highlight the importance of AP1/SQUA subfamily in biological evolution and diversity.

Overexpression of MADS-box Gene AGAMOUS-LIKE 12 Activates Root Development in Juglans sp. and Arabidopsis thaliana

Until recently, the roles of plant MADS-box genes have mainly been characterized during inflorescence and flower differentiation. In order to precise the roles of AGAMOUS-LIKE 12, one of the few MADS-box genes preferentially expressed in roots, we placed its cDNA under the control of the double 35S CaMV promoter to produce transgenic walnut tree and Arabidopsis plants. In Juglans sp., transgenic somatic embryos showed significantly higher germination rates but abnormal development of their shoot apex prevented their conversion into plants. In addition, a wide range of developmental abnormalities corresponding to ectopic root-like structures affected the transgenic lines suggesting partial reorientations of the embryonic program toward root differentiation. In Arabidopsis, AtAGL12 overexpression lead to the production of faster growing plants presenting dramatically wider and shorter root phenotypes linked to increased meristematic cell numbers within the root apex. In the upper part of the roots, abnormal cell divisions patterns within the pericycle layer generated large ectopic cell masses that did not prevent plants to grow. Taken together, our results confirm in both species that AGL12 positively regulates root meristem cell division and promotes overall root vascular tissue formation. Genetic engineering of AGL12 expression levels could be useful to modulate root architecture and development.

GsMAS1 Encoding a MADS-box Transcription Factor Enhances the Tolerance to Aluminum Stress in Arabidopsis thaliana

The MADS-box transcription factors (TFs) are essential in regulating plant growth and development, and conferring abiotic and metal stress resistance. This study aims to investigate GsMAS1 function in conferring tolerance to aluminum stress in Arabidopsis. The GsMAS1 from the wild soybean BW69 line encodes a MADS-box transcription factor in Glycine soja by bioinformatics analysis. The putative GsMAS1 protein was localized in the nucleus. The GsMAS1 gene was rich in soybean roots presenting a constitutive expression pattern and induced by aluminum stress with a concentration-time specific pattern. The analysis of phenotypic observation demonstrated that overexpression of GsMAS1 enhanced the tolerance of Arabidopsis plants to aluminum (Al) stress with larger values of relative root length and higher proline accumulation compared to those of wild type at the AlCl3 treatments. The genes and/or pathways regulated by GsMAS1 were further investigated under Al stress by qRT-PCR. The results indicated that six genes resistant to Al stress were upregulated, whereas AtALMT1 and STOP2 were significantly activated by Al stress and GsMAS1 overexpression. After treatment of 50 μM AlCl3, the RNA abundance of AtALMT1 and STOP2 went up to 17-fold and 37-fold than those in wild type, respectively. Whereas the RNA transcripts of AtALMT1 and STOP2 were much higher than those in wild type with over 82% and 67% of relative expression in GsMAS1 transgenic plants, respectively. In short, the results suggest that GsMAS1 may increase resistance to Al toxicity through certain pathways related to Al stress in Arabidopsis.

Expression Pattern and Functional Characterization of PISTILLATA Ortholog Associated With the Formation of Petaloid Sepals in Double-Flower Eriobotrya japonica (Rosaceae)

Double-flower Eriobotrya japonica, of which one phenotype is homeotic transformation of sepals into petals, is a new germplasm for revealing the molecular mechanisms underlying the floral organ transformation. Herein, we analyzed the sequence, expression pattern and functional characterization of EjPI, which encoded a B-class floral homeotic protein referred to as PISTILLATA ortholog, from genetically cognate single-flower and double-flower E. japonica. Phylogenetic analysis suggested that the EjPI gene was assigned to the rosids PI/GLO lineage. Analysis of protein sequence alignments showed that EjPI has typical domains of M, I, K, and C, and includes a distinctive PI motif at the C-terminal region. Compared with asterids PI/GLO lineage, the K1 and K3 subdomains of EjPI both contain a single amino acid difference. Subcellular localization of EjPI was determined to be in the nucleus. Expression pattern analysis revealed that EjPI expressed not only in petals, filament, and anther in single-flower E. japonica, but also in petaloid sepals in double-flower E. japonica. Meanwhile, there were high correlation between EjPI transcript level and petaloid area within a sepal. Furthermore, 35S::EjPI transgenic wild-type Arabidopsis caused the homeotic transformation of the first whorl sepals into petaloid sepals. Ectopic expression of EjPI in transgenic pi-1 mutant Arabidopsis rescued normal petals and stamens. These results suggest expression pattern of EjPI is associated with the formation of petaloid sepal. Our study provides the potential application of EjPI for biotechnical engineering to create petaloid sepals or regulate floral organ identity in angiosperms.

Cloning and functional characterization of a floral repressor gene from Lavandula angustifolia

Using RNA-Seq, we identified genes involved in floral development in lavenders and functionally characterized the floral repressor LaSVP. The molecular aspects of flower initiation and development have not been adequately investigated in lavender (Lavandula). In order to identify genes that control these processes, we employed RNA-Seq to obtain sequence information for transcripts originating from the vegetative shoot apical meristem (SAM) and developing inflorescence tissues of Lavandula angustifolia and Lavandula × intermedia plants, and assemble a comprehensive transcriptome of 105,294 contigs. Homology-based annotation provided gene ontology terms for the majority of transcripts, including over 100 genes homologous to those that control flower initiation and organ identity in Arabidopsis thaliana. Expression analysis revealed that most of these genes are differentially expressed during flower development. For example, LaSVP, a homolog of the floral repressor SHORT VEGETATIVE PHASE (SVP), was strongly expressed in vegetative SAM compared to developing flowers, implicating its potential involvement in flowering repression in lavender. To investigate LaSVP further, we constitutively expressed the gene in transformed A. thaliana plants, evaluating its effects on flower initiation and morphology. Expression of the LaSVP in A. thaliana delayed flowering and affected flower organ identity in a dosage-dependent manner. Two of the highest expressing lines produced sepals instead of petals and were sterile as they failed to develop proper seed pods. This study provides the foundation for future investigations aimed at elucidating flower initiation and development in lavender.

The Role of EjSVPs in Flower Initiation in Eriobotrya japonica

Flowering plants have evolved different flowering habits to sustain long-term reproduction. Most woody trees experience dormancy and then bloom in the warm spring, but loquat blooms in the cold autumn and winter. To explore its mechanism of flowering regulation, we cloned two SHORT VEGETATIVE PHASE (SVP) homologous genes from 'Jiefanzhong' loquat (Eriobotrya japonica Lindl.), namely, EjSVP1 and EjSVP2. Sequence analysis revealed that the EjSVPs were typical MADS-box transcription factors and exhibited a close genetic relationship with other plant SVP/DORMANCY-ASSOCIATED MADS-BOX (DAM) proteins. The temporal and spatial expression patterns showed that EjSVP1 and EjSVP2 were mainly expressed in the shoot apical meristem (SAM) after the initiation of flowering; after reaching their highest level, they gradually decreased with the development of the flower until they could not be detected. EjSVP1 expression levels were relatively high in young tissues, and EjSVP2 expression levels were relatively high in young to mature transformed tissues. Interestingly, EjSVP2 showed relatively high expression levels in various flower tissues. We analyzed the EjSVP promoter regions and found that they did not contain the C-repeat/dehydration-responsive element. Finally, we overexpressed the EjSVPs in wild-type Arabidopsis thaliana Col-0 and found no significant changes in the number of rosette leaves of Arabidopsis thaliana; however, overexpression of EjSVP2 affected the formation of Arabidopsis thaliana flower organs. In conclusion, EjSVPs were found to play an active role in the development of loquat flowering. These findings may provide a reference for exploring the regulation mechanisms of loquat flowering and the dormancy mechanisms of other plants.

Genome-wide Analysis of the MADS-Box Gene Family in Watermelon

MADS-box genes comprise a family of transcription factors that function in the growth and development of plants. To obtain insights into their evolution in watermelon (Citrullus lanatus), we carried out a genome-wide analysis and identified 39 MADS-box genes. These genes were classified into MIKC^c (25), MIKC^*(3), Mα (5), Mβ (3), and Mγ (3) clades according to their phylogenetic relationship with Arabidopsis thaliana and Cucumis sativus; moreover, these 25 genes in the MIKC clade could be classified into 13 subfamilies, and the Flowering Locus C (FLC) subfamily is absent in watermelon. Analysis of the conserved gene motifs showed similar motifs among clades. Continuing chromosomal localizations analysis indicated that MADS-box genes were distributed across 11 chromosomes in watermelon, and these genes were conditioned to be differentially expressed during plant growth and development. This research provides information that will aid further investigations into the evolution of the MADS-box gene family in plants.

Origination and selection of ABCDE and AGL6 subfamily MADS-box genes in gymnosperms and angiosperms

BACKGROUND

The morphological diversity of flower organs is closely related to functional divergence within the MADS-box gene family. Bryophytes and seedless vascular plants have MADS-box genes but do not have ABCDE or AGAMOUS-LIKE6 (AGL6) genes. ABCDE and AGL6 genes belong to the subgroup of MADS-box genes. Previous works suggest that the B gene was the first ABCDE and AGL6 genes to emerge in plant but there are no mentions about the probable origin time of ACDE and AGL6 genes. Here, we collected ABCDE and AGL6 gene 381 protein sequences and 361 coding sequences from gymnosperms and angiosperms and reconstructed a complete Bayesian phylogeny of these genes. In this study, we want to clarify the probable origin time of ABCDE and AGL6 genes is a great help for understanding the role of the formation of the flower, which can decipher the forming order of MADS-box genes in the future.

RESULTS

These genes appeared to have been under purifying selection and their evolutionary rates are not significantly different from each other. Using the Bayesian evolutionary analysis by sampling trees (BEAST) tool, we estimated that: the mutation rate of the ABCDE and AGL6 genes was 2.617 × 10-3 substitutions/site/million years, and that B genes originated 339 million years ago (MYA), CD genes originated 322 MYA, and A genes shared the most recent common ancestor with E/AGL6 296 MYA, respectively.

CONCLUSIONS

The phylogeny of ABCDE and AGL6 genes subfamilies differed. The APETALA1 (AP1 or A gene) subfamily clustered into one group. The APETALA3/PISTILLATA (AP3/PI or B genes) subfamily clustered into two groups: the AP3 and PI clades. The AGAMOUS/SHATTERPROOF/SEEDSTICK (AG/SHP/STK or CD genes) subfamily clustered into a single group. The SEPALLATA (SEP or E gene) subfamily in angiosperms clustered into two groups: the SEP1/2/4 and SEP3 clades. The AGL6 subfamily clustered into a single group. Moreover, ABCDE and AGL6 genes appeared in the following order: AP3/PI → AG/SHP/STK → AGL6/SEP/AP1. In this study, we collected candidate sequences from gymnosperms and angiosperms. This study highlights important events in the evolutionary history of the ABCDE and AGL6 gene families and clarifies their evolutionary path.

Asymmetric birth and death of type I and type II MADS-box gene subfamilies in the rubber tree facilitating laticifer development

The rubber tree (Hevea brasiliensis Muell. Arg.) is a rubber producing crop and contains specialized laticifers. MADS-box genes are a family of transcription factor genes that regulate plant development, especially floral organ and gametophyte development. 97 MADS-box genes were identified in the rubber tree through transcriptomes and genome mining. 93.8% of the genes were mapped onto the genome scaffolds in correspondence to the coverage (93.8%) of current version of sequenced genome. Phylogenetic analysis indicates that type II MADS-box genes have been more actively duplicated than their orthologous genes in Arabidopsis and rice, so that most (70, 72.2%) of the MADS-box genes in the rubber tree belong to type II subfamily. This is a high percentage compared to those in Arabidopsis (43.7%) and rice (56.8%). Moreover, 69 out of 70 type II genes in the rubber tree are transcribed, and they are mostly predominantly expressed in flowers, but some genes are predominantly expressed in laticifers, suggesting their roles in both flower and laticifer development. The number of type I genes in the rubber tree is only 27 (27.8%), a much smaller number compared to their orthologous genes in Arabidopsis (56.3%) and rice (43.2%). At the same time, most of the type I genes (55.6%, 15) in the rubber tree are silent and are probably pseudogenes. The high birth rate and low death rate of type II genes and low birth rate and high death rate of type I genes may corresponds to special developmental requirements in the rubber tree, e.g. the development of laticifer system for biosynthesis of cis-polyisoprene, the rubber. Moreover, atypical MIKC* factors (e.g. HbMADS1 in S-clade, and HbMADS20 in P-clade) are identified. These genes are diverged to typical MIKC* genes in sequences and facilitate functions required in laticifer development and rubber biosynthesis, which is not necessary in Arabidopsis and rice.