Gene Duplication and the Evolution of Plant MADS-box Transcription Factors

Evolutionary analysis of the DHHCs in Saccharinae

The DHHC domain genes are crucial for protein lipid modification, a key post-translational modification influencing membrane targeting, subcellular trafficking, and protein function. Despite their significance, the DHHC gene family in Saccharinae remains understudied. Here, we identified 32 (110 alleles), 28, 53, and 48 DHHC genes in Saccharum spontaneum Np-X, Erianthus rufipilus, Miscanthus sinensis, and Miscanthus lutarioriparius, respectively. Collinearity analysis uncovered the loss of two M. lutarioriparius genes, homologues of EruDHHC1C and EruDHHC3A. Phylogenetic and classification analyses categorized DHHC family members into six subgroups (A-F). Ka/Ks ratio analysis indicated that gene duplication in these species was primarily driven by whole-genome duplication (WGD) and dispersed duplication (DSD), with DHHC genes evolving under strong purifying selection. Gene expression and trait correlation analysis revealed a significant negative correlation between SspDHHC28A expression in S. spontaneum and sucrose content, suggesting a role in photosynthesis product transport during rapid growth. This study deepens our understanding of the DHHC gene family's functional dynamics and evolutionary path in Saccharinae, laying a foundation for future research.

Evolutionary trajectory of transcription factors and selection of targets for metabolic engineering

Transcription factors (TFs) provide potentially powerful tools for plant metabolic engineering as they often control multiple genes in a metabolic pathway. However, selecting the best TF for a particular pathway has been challenging, and the selection often relies significantly on phylogenetic relationships. Here, we offer examples where evolutionary relationships have facilitated the selection of the suitable TFs, alongside situations where such relationships are misleading from the perspective of metabolic engineering. We argue that the evolutionary trajectory of a particular TF might be a better indicator than protein sequence homology alone in helping decide the best targets for plant metabolic engineering efforts. This article is part of the theme issue 'The evolution of plant metabolism'.

A molecular module connects abscisic acid with auxin signals to facilitate seasonal wood formation in Populus

Perennial trees have a recurring annual cycle of wood formation in response to environmental fluctuations. However, the precise molecular mechanisms that regulate the seasonal formation of wood remain poorly understood. Our prior study indicates that VCM1 and VCM2 play a vital role in regulating the activity of the vascular cambium by controlling the auxin homoeostasis of the cambium zone in Populus. This study indicates that abscisic acid (ABA) affects the expression of VCM1 and VCM2, which display seasonal fluctuations in relation to photoperiod changes. ABA-responsive transcription factors AREB4 and AREB13, which are predominantly expressed in stem secondary vascular tissue, bind to VCM1 and VCM2 promoters to induce their expression. Seasonal changes in the photoperiod affect the ABA amount, which is linked to auxin-regulated cambium activity via the functions of VCM1 and VCM2. Thus, the study reveals that AREB4/AREB13-VCM1/VCM2-PIN5b acts as a molecular module connecting ABA and auxin signals to control vascular cambium activity in seasonal wood formation.

Decrease in purifying selection pressures on wheat homoeologous genes: tetraploidization versus hexaploidization

A series of polyploidizations in higher-order polyploids is the main event affecting gene content in a genome. Each polyploidization event can lead to massive functional divergence because of the subsequent decrease in selection pressure on duplicated genes; however, the causal relationship between multiple rounds of polyploidization and the functional divergence of duplicated genes is poorly understood. We focused on the Triticum-Aegilops complex lineage and compared selection pressure before and after tetraploidization and hexaploidization events. Although both events led to decreased selection pressure on homoeologous gene pairs (compared with diploids and tetraploids), the initial tetraploidization had a greater impact on selection pressure on homoeologous gene pairs than did subsequent hexaploidization. Consistent with this, selection pressure on expression patterns for the initial event relaxed more than those for the subsequent event. Surprisingly, the decreased selection pressure on these homoeologous genes was independent of the existence of in-paralogs within the same subgenome. Wheat homoeologous pairs had different evolutionary consequences compared with orthologs related to other mechanisms (ancient allopolyploidization, ancient autopolyploidization, and small-scale duplication). Furthermore, tetraploidization and hexaploidization also seemed to have different evolutionary consequences. This suggests that homoeologous genes retain unique functions, including functions that are unlikely to be preserved in genes generated by the other duplication mechanisms. We found that their unique functions differed between tetraploidization and hexaploidization (e.g., reproductive and chromosome segregation processes). These findings imply that the substantial number of gene pairs resulting from multiple allopolyploidization events, especially initial tetraploidization, may have been a unique source of functional divergence.

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.

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'.

Genome-wide identification and expression analysis of the C2H2-zinc finger transcription factor gene family and screening of candidate genes involved in floral development in Coptis teeta Wall. (Ranunculaceae)

Background: C2H2-zinc finger transcription factors comprise one of the largest and most diverse gene superfamilies and are involved in the transcriptional regulation of flowering. Although a large number of C2H2 zinc-finger proteins (C2H2-ZFPs) have been well characterized in a number of model plant species, little is known about their expression and function in Coptis teeta. C. teeta displays two floral phenotypes (herkogamy phenotypes). It has been proposed that the C2H2-zinc finger transcription factor family may play a crucial role in the formation of floral development and herkogamy observed in C. teeta. As such, we performed a genome-wide analysis of the C2H2-ZFP gene family in C. teeta. Results: The complexity and diversity of C. teeta C2H2 zinc finger proteins were established by evaluation of their physicochemical properties, phylogenetic relationships, exon-intron structure, and conserved motifs. Chromosome localization showed that 95 members of the C2H2 zinc-finger genes were unevenly distributed across the nine chromosomes of C. teeta, and that these genes were replicated in tandem and segmentally and had undergone purifying selection. Analysis of cis-acting regulatory elements revealed a possible involvement of C2H2 zinc-finger proteins in the regulation of phytohormones. Transcriptome data was then used to compare the expression levels of these genes during the growth and development of the two floral phenotypes (F-type and M-type). These data demonstrate that in groups A and B, the expression levels of 23 genes were higher in F-type flowers, while 15 genes showed higher expressions in M-type flowers. qRT-PCR analysis further revealed that the relative expression was highly consistent with the transcriptome data. Conclusion: These data provide a solid basis for further in-depth studies of the C2H2 zinc finger transcription factor gene family in this species and provide preliminary information on which to base further research into the role of the C2H2 ZFPs gene family in floral development in C. teeta.

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

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

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

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

Genome-wide identification of glutamate receptor genes in adzuki bean and the roles of these genes in light and rust fungal response

Plant glutamate receptor proteins (GLRs) are involved in plant development, biotic stress, and light-signal transduction. Vigna angularis is a traditional crop with important economic value in China, and the identification of functional genes can facilitate the breeding of stress resistant varieties. Here, we identified the members of the GLR gene family in the adzuki bean genome and investigated gene expression under light and rust fungal (Uromyces vignae) stimuli. Sixteen GLR genes were identified in V. angularis (VaGLRs), and these genes clustered in a single clade (clade III) with two groups. Evolutionary analysis showed that three VaGLRs result from tandem duplications and four result from whole genome/segmental duplications. To understand the regulation of expression of VaGLRs, cis-acting elements were analyzed in the promoter regions of the VaGLRs including cis-acting elements associated with light and stress responsiveness. Expression analysis by qRT-PCR revealed transcripts of eight and 10 VaGLRs in response to light stimuli and rust infection, respectively. For light responsiveness, the expression levels of XP_017430569.1 and XP_017425299.1 were higher under light condition than in darkness, while the expression levels of XP_017406996.1, XP_017425763.1, and XP_017423557.1 gradually recovered during dark treatment. Additionally, the relative expression levels of XP_017413816.1, XP_017436268.1, and XP_017425299.1 were significantly elevated during U. vignae infection in a resistant cultivar compared to the expression levels in a susceptible cultivar. XP_017425299.1 expression was induced both by light and rust infection, suggesting this gene may link light and disease resistance signaling pathways. Our results provide insight into how the VaGLRs contribute to adzuki bean response to light stimulus and pathogen attack. These identified VaGLRs also provide important reference to improve adzuki bean germplasm resources.

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 characterization of R2R3-MYB gene family in Santalum album and their expression analysis under cold stress

Sandalwood (Santalum album) is a high-value multifunctional tree species that is rich in aromatic substances and is used in medicine and global cosmetics. Due to the scarcity of land resources in tropical and subtropical regions, land in temperate regions is a potential resource for the development of S. album plantations in order to meet the needs of S. album production and medicine. The R2R3-MYB transcription factor family is one of the largest in plants and plays an important role in the response to various abiotic stresses. However, the R2R3-MYB gene family of S. album has not been studied. In this study, 144 R2R3-MYB genes were successfully identified in the assembly genome sequence, and their characteristics and expression patterns were investigated under various durations of low temperature stress. According to the findings, 31 of the 114 R2R3-MYB genes showed significant differences in expression after cold treatment. Combining transcriptome and weighted gene co-expression network analysis (WGCNA) revealed three key candidate genes (SaMYB098, SaMYB015, and SaMYB068) to be significantly involved in the regulation of cold resistance in S. album. The structural characteristics, evolution, and expression pattern of the R2R3-MYB gene in S. album were systematically examined at the whole genome level for the first time in this study. It will provide important information for future research into the function of the R2R3-MYB genes and the mechanism of cold stress response in S. album.

Gene Structural Specificity and Expression of MADS-Box Gene Family in Camellia chekiangoleosa

MADS-box genes encode transcription factors that affect plant growth and development. Camellia chekiangoleosa is an oil tree species with ornamental value, but there have been few molecular biological studies on the developmental regulation of this species. To explore their possible role in C. chekiangoleosa and lay a foundation for subsequent research, 89 MADS-box genes were identified across the whole genome of C. chekiangoleosa for the first time. These genes were present on all the chromosomes and were found to have expanded by tandem duplication and fragment duplication. Based on the results of a phylogenetic analysis, the 89 MADS-box genes could be divided into either type I (38) or type II (51). Both the number and proportion of the type II genes were significantly greater than those of Camellia sinensis and Arabidopsis thaliana, indicating that C. chekiangoleosa type II genes experienced a higher duplication rate or a lower loss rate. The results of both a sequence alignment and a conserved motif analysis suggest that the type II genes are more conserved, meaning that they may have originated and differentiated earlier than the type I genes did. At the same time, the presence of extra-long amino acid sequences may be an important feature of C. chekiangoleosa. Gene structure analysis revealed the number of introns of MADS-box genes: twenty-one type I genes had no introns, and 13 type I genes contained only 1~2 introns. The type II genes have far more introns and longer introns than the type I genes do. Some MIKC^C genes have super large introns (≥15 kb), which are rare in other species. The super large introns of these MIKC^C genes may indicate richer gene expression. Moreover, the results of a qPCR expression analysis of the roots, flowers, leaves and seeds of C. chekiangoleosa showed that the MADS-box genes were expressed in all those tissues. Overall, compared with that of the type I genes, the expression of the type II genes was significantly higher. The CchMADS31 and CchMADS58 genes (type II) were highly expressed specifically in the flowers, which may in turn regulate the size of the flower meristem and petals. CchMADS55 was expressed specifically in the seeds, which might affect seed development. This study provides additional information for the functional characterization of the MADS-box gene family and lays an important foundation for in-depth study of related genes, such as those involved in the development of the reproductive organs of C. chekiangoleosa.

Genome-wide identification of MADS-box gene family in orchardgrass and the positive role of DgMADS114 and DgMADS115 under different abiotic stress

Abiotic stress, a major factor limit growth and productivity of major crops. Orchardgrass is one of the most important cool-season forage grasses in the world, and it is highly tolerant to abiotic stress. The MADS-box transcription factor family is one of the largest families in plants, and it plays vital roles in multiple biological processes. However, MADS-box transcription factors in orchardgrass, especially those involved in abiotic stress, have not yet been elucidated. Here, 123 DgMADS-box members were identified in orchardgrass and a detailed overview has been presented. Syntenic analysis indicated that the expansion of the DgMADS-box genes in orchardgrass is mainly dependent on tandem duplication events. Some DgMADS-box genes were induced by multiple abiotic stresses, indicating that these genes may play critical regulatory roles in orchardgrass response to various abiotic stresses. Heterologous expression showed that DgMADS114 and DgMADS115 could enhance stress tolerance of transgenic Arabidopsis, as revealed by longer root length or higher survival rates under PEG, NaCl, ABA, and heat stress. The results of this study provide a scientific basis for clarifying the functional characterization of MADS-box genes in orchardgrass in response to environmental stress can be further used to improve forages and crops via breeding programs.

Functional Divergence Analysis of AGL6 Genes in Prunus mume

The AGAMOUS-LIKE6 (AGL6) lineage is an important clade of MADS-box transcription factors that play essential roles in floral organ development. The genome of Prunus mume contains two homoeologous AGL6 genes that are replicated as gene fragments. In this study, two AGL6 homologs, PmAGL6-1 and PmAGL6-2, were cloned from P. mume and then functionally characterized. Sequence alignment and phylogenetic analyses grouped both genes into the AGL6 lineage. The expression patterns and protein-protein interaction patterns showed significant differences between the two genes. However, the ectopic expression of the two genes in Arabidopsis thaliana resulted in similar phenotypes, including the promotion of flowering, alteration of floral organ structure, participation in the formation of the floral meristem and promotion of pod bending. Therefore, gene duplication has led to some functional divergence of PmAGL6-1 and PmAGL6-2 but their functions are similar. We thus speculated that AGL6 genes play a crucial role in flower development in P. mume.

Genome-wide analysis of MADS-box gene family in kiwifruit (Actinidia chinensis var. chinensis) and their potential role in floral sex differentiation

Kiwifruit (Actinidia chinensis Planch.) is a functionally dioecious plant, which displays diverse morphology in male and female flowers. MADS-box is an ancient and huge gene family that plays a key role in plant floral organ differentiation. In this study, we have identified 89 MADS-box genes from A. chinensis Red 5 genome. These genes are distributed on 26 chromosomes and are classified into type I (21 genes) and type II (68 genes). Overall, type II AcMADS-box genes have more complex structures than type I with more exons, protein domains, and motifs, indicating that type II genes may have more diverse functions. Gene duplication analysis showed that most collinearity occurred in type II AcMADS-box genes, which was consistent with a large number of type II genes. Analysis of cis-acting elements in promoters showed that AcMADS-box genes are mainly associated with light and phytohormone responsiveness. The expression profile of AcMADS-box genes in different tissues showed that most genes were highly expressed in flowers. Further, the qRT-PCR analysis of the floral organ ABCDE model-related genes in male and female flowers revealed that AcMADS4, AcMADS56, and AcMADS70 were significantly expressed in female flowers. It indicated that those genes may play an important role in the sex differentiation of kiwifruit. This work provided a comprehensive analysis of the AcMADS-box genes and may help facilitate our understanding of the sex differentiation regulatory mechanism in kiwifruit.

Tracing the Evolution of the SEPALLATA Subfamily across Angiosperms Associated with Neo- and Sub-Functionalization for Reproductive and Agronomically Relevant Traits

SEPALLATA transcription factors (SEP TFs) have been extensively studied in angiosperms as pivotal components of virtually all the MADS-box tetrameric complex master regulators of floral organ identities. However, there are published reports that suggest that some SEP members also regulate earlier reproductive events, such as inflorescence meristem determinacy and inflorescence architecture, with potential for application in breeding programs in crops. The SEP subfamily underwent a quite complex pattern of duplications during the radiation of the angiosperms. Taking advantage of the many whole genomic sequences now available, we present a revised and expanded SEP phylogeny and link it to the known functions of previously characterized genes. This snapshot supports the evidence that the major SEP3 clade is highly specialized for the specification of the three innermost floral whorls, while its sister LOFSEP clade is functionally more versatile and has been recruited for diverse roles, such as the regulation of extra-floral bract formation and inflorescence determinacy and shape. This larger pool of angiosperm SEP genes confirms previous evidence that their evolution was driven by whole-genome duplications rather than small-scale duplication events. Our work may help to identify those SEP lineages that are the best candidates for the improvement of inflorescence traits, even in far distantly related crops.

Genome-wide analysis of MADS-box families and their expressions in flower organs development of pineapple (Ananas comosus (L.) Merr.)

MADS-box genes play crucial roles in plant vegetative and reproductive growth, better development of inflorescences, flower, and fruit. Pineapple is a typical collective fruit, and a comprehensive analysis of the MADS-box gene family in the development of floral organs of pineapple is still lacking. In this study, the whole-genome survey and expression profiling of the MADS-box family in pineapple were introduced. Forty-four AcMADS genes were identified in pineapple, 39 of them were located on 18 chromosomes and five genes were distributed in five scaffolds. Twenty-two AcMADS genes were defined as 15 pairs of segmental duplication events. Most members of the type II subfamily of AcMADS genes had higher expression levels in floral organs compared with type I subfamily, thereby suggesting that AcMADS of type II may play more crucial roles in the development of floral organs of pineapple. Six AcMADS genes have significant tissue-specificity expression, thereby suggesting that they may participate in the formation of one or more floral organs. This study provides valuable insights into the role of MADS-box gene family in the floral organ development of pineapple.

Genome-wide identification, structural analysis and expression profiles of short internodes related sequence gene family in quinoa

Based on the whole genome data information of Chenopodium quinoa Willd, the CqSRS gene family members were systematically identified and analyzed by bioinformatics methods, and the responses of CqSRS genes to NaCl (100 mmol/L), salicylic acid (200 umol/L) and low temperature (4°C) were detected by qRT-PCR. The results showed that a total of 10 SHI related sequence genes were identified in quinoa, and they were distributed on 9 chromosomes, and there were four pairs of duplicated genes. The number of amino acids encoded ranged from 143 aa to 370 aa, and the isoelectric point ranged from 4.81 to 8.90. The secondary structure was mainly composed of random coil (Cc). Most of the SRS gene encoding proteins were located in the cytoplasm (5 CqSRS). Phylogenetic analysis showed that the CqSRS genes were divided into three groups, and the gene structure showed that the number of exons of CqSRS was between two-five. Promoter analysis revealed that there are a total of 44 elements related to plant hormone response elements, light response elements, stress response elements and tissue-specific expression in the upstream regin of the gene. Protein interaction showed that all 10 CqSRS proteins appeared in the known protein interaction network diagram in Arabidopsis. Expression profile analysis showed that CqSRS genes had different expression patterns, and some genes had tissue-specific expression. qRT-PCR showed that all SRS family genes responded to ABA、NaCl、drought and low-temperature treatments, but the expression levels of different CqSRS genes were significantly different under various stresses. This study lays a foundation for further analyzed the function of CqSRS genes.

Deep evaluation of the evolutionary history of the Heat Shock Factor (HSF) gene family and its expansion pattern in seed plants

Heat shock factor (HSF) genes are essential in some of the basic developmental pathways in plants. Despite extensive studies on the structure, functional diversification, and evolution of HSF genes, their divergence history and gene duplication pattern remain unknown. To further illustrate the probable divergence patterns in these subfamilies, we analyzed the evolutionary history of HSF genes using phylogenetic reconstruction and genomic syntenic analyses, taking advantage of the increased sampling of genomic data from pteridophytes, gymnosperms and basal angiosperms. We identified a novel clade that includes HSFA2, HSFA6, HSFA7, and HSFA9 with a complex relationship, which is very likely due to orthologous or paralogous genes retained after frequent gene duplication events. We hypothesized that HSFA9 derives from HSFA2 through gene duplication in eudicots at the ancestral state, and then expanded in a lineage-specific way. Our findings indicate that HSFB3 and HSFB5 emerged before the divergence of ancestral angiosperms, but were lost in the most recent common ancestors of monocots. We also presumed that HSFC2 derives from HSFC1 in ancestral monocots. This work proposes that during the radiation of flowering plants, an era during which there was a differentiation of angiosperms, the size of the HSF gene family was also being adjusted with considerable sub- or neo-functionalization. The independent evolution of HSFs in eudicots and monocots, including lineage-specific gene duplication, gave rise to a new gene in ancestral eudicots and monocots, and lineage-specific gene loss in ancestral monocots. Our analyses provide essential insights for studying the evolutionary history of this multigene family.

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.

Plant protein-coding gene families: Their origin and evolution

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses.

Genome-Wide Identification and Analysis of the MADS-Box Gene Family in Theobroma cacao

The MADS-box family gene is a class of transcription factors that have been extensively studied and involved in several plant growth and development processes, especially in floral organ specificity, flowering time and initiation and fruit development. In this study, we identified 69 candidate MADS-box genes and clustered these genes into five subgroups (Mα: 11; Mβ: 2; Mγ: 14; Mδ: 9; MIKC: 32) based on their phylogenetical relationships with Arabidopsis. Most TcMADS genes within the same subgroup showed a similar gene structure and highly conserved motifs. Chromosomal distribution analysis revealed that all the TcMADS genes were evenly distributed in 10 chromosomes. Additionally, the cis-acting elements of promoter, physicochemical properties and subcellular localization were also analyzed. This study provides a comprehensive analysis of MADS-box genes in Theobroma cacao and lays the foundation for further functional research.

Isolation and Functional Analysis of a PISTILLATA-like MADS-Box Gene from Argan Tree (Argania spinosa)

Argan trees (Argania spinosa) belong to a species native to southwestern Morocco, playing an important role in the environment and local economy. Argan oil extracted from kernels has a unique composition and properties. Argan trees were introduced in Tunisia, where hundreds of trees can be found nowadays. In this study, we examined reproductive development in Argan trees from four sites in Tunisia and carried out the functional characterization of a floral homeotic gene in this non-model species. Despite the importance of reproductive development, nothing is known about the genetic network controlling flower development in Argania spinosa. Results obtained in several plant species established that floral organ development is mostly controlled by MADS-box genes and, in particular, APETALA3 (AP3) and PISTILLATA (PI) homologs are required for proper petal and stamen identity. Here, we describe the isolation and functional characterization of a MADS-box gene from Argania spinosa. Phylogenetic analyses showed strong homology with PI-like proteins, and the expression of the gene was found to be restricted to the second and third whorls. Functional homology with Arabidopsis PI was demonstrated by the ability of AsPI to confer petal and stamen identity when overexpressed in a pi-1 mutant background. The identification and characterization of this gene support the strong conservation of PI homologs among distant angiosperm plants.

Genome-Wide Analysis of MADS-Box Genes in Foxtail Millet (Setaria italica L.) and Functional Assessment of the Role of SiMADS51 in the Drought Stress Response

MADS-box transcription factors play vital roles in multiple biological processes in plants. At present, a comprehensive investigation into the genome-wide identification and classification of MADS-box genes in foxtail millet (Setaria italica L.) has not been reported. In this study, we identified 72 MADS-box genes in the foxtail millet genome and give an overview of the phylogeny, chromosomal location, gene structures, and potential functions of the proteins encoded by these genes. We also found that the expression of 10 MIKC-type MADS-box genes was induced by abiotic stresses (PEG-6000 and NaCl) and exogenous hormones (ABA and GA), which suggests that these genes may play important regulatory roles in response to different stresses. Further studies showed that transgenic Arabidopsis and rice (Oryza sativa L.) plants overexpressing SiMADS51 had reduced drought stress tolerance as revealed by lower survival rates and poorer growth performance under drought stress conditions, which demonstrated that SiMADS51 is a negative regulator of drought stress tolerance in plants. Moreover, expression of some stress-related genes were down-regulated in the SiMADS51-overexpressing plants. The results of our study provide an overall picture of the MADS-box gene family in foxtail millet and establish a foundation for further research on the mechanisms of action of MADS-box proteins with respect to abiotic stresses.

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.

Draft genome of Korthalsia laciniosa (Griff.) Mart., a climbing rattan elucidates its phylogenetic position

Korthalsia laciniosa (Griff.) Mart. is a climbing rattan used as a source of durable and flexible cane. In the present study, the draft genome of K. laciniosa was sequenced, de novo assembled and annotated. Genome-wide identification of MADS-Box transcription factors revealed loss of Mβ, and Mγ genes belonging to Type I subclass in the rattan lineage. Mining of the genome revealed presence of 13 families of lignin biosynthetic pathway genes and expression profiling of nine major genes documented relatively lower level of expression in cirrus when compared to leaflet and petiole. The chloroplast genome was re-constructed and analysis revealed the phylogenetic relatedness of this genus to Eugeissona, in contrast with its present taxonomic position. The genomic resource generated in the present study will accelerate population structure analysis, genetic resource conservation, phylogenomics and facilitate understanding the unique developmental processes like gender expression at molecular level.

Genome-Wide Identification and Functional Characterization of the Cation Proton Antiporter (CPA) Family Related to Salt Stress Response in Radish (Raphanus sativus L.)

The CPA (cation proton antiporter) family plays an essential role during plant stress tolerance by regulating ionic and pH homeostasis of the cell. Radish fleshy roots are susceptible to abiotic stress during growth and development, especially salt stress. To date, CPA family genes have not yet been identified in radish and the biological functions remain unclear. In this study, 60 CPA candidate genes in radish were identified on the whole genome level, which were divided into three subfamilies including the Na+/H+ exchanger (NHX), K+ efflux antiporter (KEA), and cation/H+ exchanger (CHX) families. In total, 58 of the 60 RsCPA genes were localized to the nine chromosomes. RNA-seq. data showed that 60 RsCPA genes had various expression levels in the leaves, roots, cortex, cambium, and xylem at different development stages, as well as under different abiotic stresses. RT-qPCR analysis indicated that all nine RsNHXs genes showed up regulated trends after 250 mM NaCl exposure at 3, 6, 12, and 24h. The RsCPA31 (RsNHX1) gene, which might be the most important members of the RsNHX subfamily, exhibited obvious increased expression levels during 24h salt stress treatment. Heterologous over-and inhibited-expression of RsNHX1 in Arabidopsis showed that RsNHX1 had a positive function in salt tolerance. Furthermore, a turnip yellow mosaic virus (TYMV)-induced gene silence (VIGS) system was firstly used to functionally characterize the candidate gene in radish, which showed that plant with the silence of endogenous RsNHX1 was more susceptible to the salt stress. According to our results we provide insights into the complexity of the RsCPA gene family and a valuable resource to explore the potential functions of RsCPA genes in radish.

A MADS-box gene, EgMADS21, negatively regulates EgDGAT2 expression and decreases polyunsaturated fatty acid accumulation in oil palm (Elaeis guineensis Jacq.)

EgMADS21 regulates PUFA accumulation in oil palm. Oil palm (Elaeis guineensis Jacq.) is the most productive world oil crop, accounting for 36% of world plant oil production. However, the molecular mechanism of the transcriptional regulation of fatty acid accumulation and lipid synthesis in the mesocarp of oil palm by up- or downregulating the expression of genes involved in related pathways remains largely unknown. Here, an oil palm MADS-box gene, EgMADS21, was screened in a yeast one-hybrid assay using the EgDGAT2 promoter sequence as bait. EgMADS21 is preferentially expressed in early mesocarp developmental stages in oil palm fruit and presents a negative correlation with EgDGAT2 expression. The direct binding of EgMADS21 to the EgDGAT2 promoter was confirmed by electrophoretic mobility shift assay. Subsequently, transient expression of EgMADS21 in oil palm protoplasts revealed that EgMADS21 not only binds to the EgDGAT2 promoter but also negatively regulates the expression of EgDGAT2. Furthermore, EgMADS21 was stably overexpressed in transgenic oil palm embryoids by Agrobacterium-mediated transformation. In three independent transgenic lines, EgDGAT2 expression was significantly suppressed by the expression of EgMADS21. The content of linoleic acid (C18:2) in the three transgenic embryoids was significantly decreased, while that of oleic acid (C18:1) was significantly increased. Combined with the substrate preference of EgDGAT2 identified in previous research, the results demonstrate the molecular mechanism by which EgMADS21 regulates EgDGAT2 expression and ultimately affects fatty acid accumulation in the mesocarp of oil palm.

Evolution of SHORT VEGETATIVE PHASE (SVP) genes in Rosaceae: Implications of lineage-specific gene duplication events and function diversifications with respect to their roles in processes other than bud dormancy

MADS-box genes that are homologous to Arabidopsis SHORT VEGETATIVE PHASE (SVP) have been shown to play key roles in the regulation of bud dormancy in perennial species, particularly in the deciduous fruit trees of Rosaceae. However, their evolutionary profiles in Rosaceae have not yet been analyzed systematically. Here, The SVP genes were found to be significantly expanded in Rosaceae when compared with annual species from Brassicaceae. Phylogenetic analysis showed that Rosaceae SVP genes could be classified into five clades, namely, SVP1, SVP2-R1, SVP2-R2, SVP2-R3 and SVP3. The SVP1 clade genes were retained in most of the species, whereas the SVP2-R2 and SVP2-R3 clades were found to be Maleae- and Amygdaleae-specific (Both of the lineages belong to Amygdaloideae), respectively, and SVP2-R1 was Rosoideae-specific in Rosaceae. Furthermore, 10 lineage-specific gene duplication (GD) events (GD1-10) were proposed for the expansion of SVP genes, suggesting that the expansion and divergence of Rosaceae SVP genes were mainly derived by lineage-specific manner during evolution. Moreover, tandem and segmental duplications were the major reasons for the expansion of SVP genes, and interestingly, tandem duplications, a well-known evolutionary feature of SVP genes, were found to be mainly Amygdaloideae-specific. Sequence alignment, selection pressure, and cis-acting element analysis suggested large functional innovations and diversification of SVP genes in different lineages of Rosaceae. Finally, the different growth cycle of Rosa multiflora and their novel expression patterns of RmSVP genes provided new insights into the functional diversification of SVP genes in terms of their roles in processes other than bud dormancy.

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.

Comparative analysis of the P-type ATPase gene family in seven Rosaceae species and an expression analysis in pear (Pyrus bretschneideri Rehd.)

P-type ATPases are integral membrane transporters that play important roles in transmembrane transport in plants. However, a comprehensive analysis of the P-type ATPase gene family has not been conducted in Chinese white pear (Pyrus bretschneideri) or other Rosaceae species. Here, we identified 419 P-type ATPase genes from seven Rosaceae species (Pyrus bretschneideri, Malus domestica, Prunus persica, Fragaria vesca, Prunus mume, Pyrus communis and Pyrus betulifolia). Structural and phylogenetic analyses revealed that P-type ATPase genes can be divided into five subfamilies. Different subfamilies have different conserved motifs and cis-acting elements, which may lead to functional divergence within one gene family. Dispersed duplication and whole-genome duplication may play critical roles in the expansion of the P-type ATPase family. Purifying selection was the primary force driving the evolution of P-type ATPase family genes. Based on the dynamic transcriptome analysis and transient transformation of Chinese white pear fruit, Pbr029767.1 in the P_3A subfamily were found to be associated with malate accumulation during pear fruit development. Using a co-expression network, we identified several transcription factors that may have regulatory relationships with the P-type ATPase gene family. Overall, this study lays a solid foundation for understanding the evolution and functions of P-type ATPase genes in Chinese white pear and six other Rosaceae species.

Identification, Characterization and Functional Analysis of C-Class Genes Associated with Double Flower Trait in Carnation (Dianthus caryphyllus L.)

Flowers with more petals are of more ornamental value. It is well known that AGAMOUS (AG) is the core member of the C-class gene which plays an essential role in double flower formation and identification of stamens and carpels in Arabidopsisthaliana. We searched C-class genes in the genome of the carnation, and found two AG orthologs (DcaAGa, DcaAGb). Phylogenetic analysis showed that the two genes were closely related to the euAG subclade. Then we searched the genomes of other Caryophyllales plants (Beta vulgaris, Spinacia oleracea, Chenopodium quinoa) for C-class genes, and found that their C-class genes all belonged to the euAG subclade. Semi-quantitative PCR (sq-PCR) analysis indicated that the expression of DcaAG genes in the single flower phenotype was higher than that in the double flower phenotype. Quantitative real-time RT-PCR (qRT-PCR) analysis showed that the expressions of DcaAG genes in the flower bud were significantly different from those in the root, stem, and leaf between the single and double flower phenotype carnations, and that DcaAG genes were specifically expressed in the stamen and carpel of carnation. Moreover, the expression of other floral organ identity genes (AP1 and AP2, PI and AP3, SEP1 and SEP3 corresponding to the A-, B-, and E-class of genes, respectively) showed no significant difference in all floral organs between the single and double flower phenotype carnations, suggesting that C-class (DcaAG) genes might play an important role in the double flower phenotype in carnation. Petal loss or decrease, precocious flowering, silique shortening, and seed sterility were observed in 35S::DcaAGa and 35S::DcaAGb transgenic Arabidopsis plants. All these results show that DcaAG genes might affect the petal number negatively and have a specific function in stamen and carpel development in carnation.

Analysis of MADS-box genes revealed modified flowering gene network and diurnal expression in pineapple

Background

Pineapple is the most important crop with CAM photosynthesis, but its molecular biology is underexplored. MADS-box genes are crucial transcription factors involving in plant development and several biological processes. However, there is no systematic analysis of MADS-box family genes in pineapple (Ananas comosus).

Results

Forty-eight MADS-box genes were identified in the pineapple genome. Based on the phylogenetic studies, pineapple MADS-box genes can be divided into type I and type II MADS-box genes. Thirty-four pineapple genes were classified as type II MADS-box genes including 32 MIKC-type and 2 Mδ-type, while 14 type I MADS-box genes were further divided into Mα, Mβ and Mγ subgroups. A majority of pineapple MADS-box genes were randomly distributed across 19 chromosomes. RNA-seq expression patterns of MADS-box genes in four different tissues revealed that more genes were highly expressed in flowers, which was confirmed by our quantitative RT-PCR results. There is no FLC and CO orthologs in pineapple. The loss of FLC and CO orthologs in pineapple indicated that modified flowering genes network in this tropical plant compared with Arabidopsis. The expression patterns of MADS-box genes in photosynthetic and non-photosynthetic leaf tissues indicated the potential roles of some MADS-box genes in pineapple CAM photosynthesis. The 23% of pineapple MADS-box genes showed diurnal rhythm, indicating that these MADS-box genes are regulated by circadian clock.

Conclusions

MADS-box genes identified in pineapple are closely related to flowering development. Some MADS-box genes are involved in CAM photosynthesis and regulated by the circadian clock. These findings will facilitate research on the development of unusual spiral inflorescences on pineapple fruit and CAM photosynthesis.

Genetic variation and temperature affects hybrid barriers during interspecific hybridization

Genomic imprinting regulates parent-specific transcript dosage during seed development and is mainly confined to the endosperm. Elucidation of the function of many imprinted genes has been hampered by the lack of corresponding mutant phenotypes, and the role of imprinting is mainly associated with genome dosage regulation or allocation of resources. Disruption of imprinted genes has also been suggested to mediate endosperm-based post-zygotic hybrid barriers depending on genetic variation and gene dosage. Here, we have analyzed the conservation of a clade from the MADS-box type I class transcription factors in the closely related species Arabidopsis arenosa, A. lyrata, and A. thaliana, and show that AGL36-like genes are imprinted and maternally expressed in seeds of Arabidopsis species and in hybrid seeds between outbreeding species. In hybridizations between outbreeding and inbreeding species the paternally silenced allele of the AGL36-like gene is reactivated in the hybrid, demonstrating that also maternally expressed imprinted genes are perturbed during hybridization and that such effects on imprinted genes are specific to the species combination. Furthermore, we also demonstrate a quantitative effect of genetic diversity and temperature on the strength of the post-zygotic hybridization barrier. Markedly, a small decrease in temperature during seed development increases the survival of hybrid F1 seeds, suggesting that abiotic and genetic parameters play important roles in post-zygotic species barriers, pointing at evolutionary scenarios favoring such effects. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA562212. All sequences generated in this study have been deposited in the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra/) with project number PRJNA562212.

Conservation and divergence of ancestral AGAMOUS/SEEDSTICK subfamily genes from the basal angiosperm Magnolia wufengensis

AGAMOUS/SEEDSTICK (AG/STK) subfamily genes play crucial roles in the reproductive development of plants. However, most of our current knowledge of AG/STK subfamily genes is restricted to core eudicots and grasses, and the knowledge of ancestral exon-intron structures, expression patterns, protein-protein interaction patterns and functions of AG/STK subfamily genes remains unclear. To determine these, we isolated AG/STK subfamily genes (MawuAG1, MawuAG2 and MawuSTK) from a woody basal angiosperm Magnolia wufengensis (Magnoliaceae). MawuSTK arose from the gene duplication event occurring before the diversification of extant angiosperms, and MawuAG1 and MawuAG2 may result from a gene duplication event occurring before the divergence of Magnoliaceae and Lauraceae. Gene duplication led to apparent diversification in their expression and interaction patterns. It revealed that expression in both stamens and carpels likely represents the ancestral expression profiles of AG lineage genes, and expression of STK-like genes in stamens may have been lost soon after the appearance of the STK lineage. Moreover, AG/STK subfamily proteins may have immediately established interactions with the SEPALLATA (SEP) subfamily proteins following the emergence of the SEP subfamily; however, their interactions with the APETALA1/FRUITFULL subfamily proteins or themselves differ from those found in monocots and basal and core eudicots. MawuAG1 plays highly conserved roles in the determinacy of stamen, carpel and ovule identity, while gene duplication contributed to the functional diversification of MawuAG2 and MawuSTK. In addition, we investigated the evolutionary history of exon-intron structural changes of the AG/STK subfamily, and a novel splice-acceptor mode (GUU-AU) and the convergent evolution of N-terminal extension in the euAG and PLE subclades were revealed for the first time. These results further advance our understanding of ancestral AG/STK subfamily genes in terms of phylogeny, exon-intron structures, expression and interaction patterns, and functions, and provide strong evidence for the significance of gene duplication in the expansion and evolution of the AG/STK subfamily.

Developmental Analysis of Mimulus Seed Transcriptomes Reveals Functional Gene Expression Clusters and Four Imprinted, Endosperm-Expressed Genes

The double fertilization of the female gametophyte initiates embryogenesis and endosperm development in seeds via the activation of genes involved in cell differentiation, organ patterning, and growth. A subset of genes expressed in endosperm exhibit imprinted expression, and the correct balance of gene expression between parental alleles is critical for proper endosperm and seed development. We use a transcriptional time series analysis to identify genes that are associated with key shifts in seed development, including genes associated with secondary cell wall synthesis, mitotic cell cycle, chromatin organization, auxin synthesis, fatty acid metabolism, and seed maturation. We relate these genes to morphological changes in Mimulus seeds. We also identify four endosperm-expressed transcripts that display imprinted (paternal) expression bias. The imprinted status of these four genes is conserved in other flowering plants, suggesting that they are functionally important in endosperm development. Our study explores gene regulatory dynamics in a species with ab initio cellular endosperm development, broadening the taxonomic focus of the literature on gene expression in seeds. Moreover, it is the first to validate genes with imprinted endosperm expression in Mimulus guttatus, and will inform future studies on the genetic causes of seed failure in this model system.

Transcriptional Structure of Petunia Clock in Leaves and Petals

The plant circadian clock coordinates environmental signals with internal processes including secondary metabolism, growth, flowering, and volatile emission. Plant tissues are specialized in different functions, and petals conceal the sexual organs while attracting pollinators. Here we analyzed the transcriptional structure of the petunia (Petunia x hybrida) circadian clock in leaves and petals. We recorded the expression of 13 clock genes in petunia under light:dark (LD) and constant darkness (DD). Under light:dark conditions, clock genes reached maximum expression during the light phase in leaves and the dark period in petals. Under free running conditions of constant darkness, maximum expression was delayed, especially in petals. Interestingly, the rhythmic expression pattern of PhLHY persisted in leaves and petals in LD and DD. Gene expression variability differed among leaves and petals, time of day and photoperiod. The transcriptional noise was higher especially in leaves under constant darkness. We found that PhPRR7, PhPRR5, and PhGI paralogs showed changes in gene structure including exon number and deletions of CCT domain of the PRR family. Our results revealed that petunia petals presented a specialized clock.

Genome-Wide Analysis of the MADS-Box Transcription Factor Family in Solanum lycopersicum

MADS-box family genes encode transcription factors that are involved in multiple developmental processes in plants, especially in floral organ specification, fruit development, and ripening. However, a comprehensive analysis of tomato MADS-box family genes, which is an important model plant to study flower fruit development and ripening, remains obscure. To gain insight into the MADS-box genes in tomato, 131 tomato MADS-box genes were identified. These genes could be divided into five groups (Mα, Mβ, Mγ, Mδ, and MIKC) and were found to be located on all 12 chromosomes. We further analyzed the phylogenetic relationships among Arabidopsis and tomato, as well as the protein motif structure and exon–intron organization, to better understand the tomato MADS-box gene family. Additionally, owing to the role of MADS-box genes in floral organ identification and fruit development, the constitutive expression patterns of MADS-box genes at different stages in tomato development were identified. We analyzed 15 tomato MADS-box genes involved in floral organ identification and five tomato MADS-box genes related to fruit development by qRT-PCR. Collectively, our study provides a comprehensive and systematic analysis of the tomato MADS-box genes and would be valuable for the further functional characterization of some important members of the MADS-box gene family.

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.

Dating Whole Genome Duplication in Ceratopteris thalictroides and Potential Adaptive Values of Retained Gene Duplicates

Whole-genome duplications (WGDs) are widespread in plants and frequently coincide with global climatic change events, such as the Cretaceous–Tertiary (KT) extinction event approximately 65 million years ago (mya). Ferns have larger genomes and higher chromosome numbers than seed plants, which likely resulted from multiple rounds of polyploidy. Here, we use diploid and triploid material from a model fern species, Ceratopteris thalictroides, for the detection of WGDs. High-quality RNA-seq data was used to infer the number of synonymous substitutions per synonymous site (Ks) between paralogs; Ks age distribution and absolute dating approach were used to determine the age of WGD events. Evidence of an ancient WGD event with a Ks peak value of approximately 1.2 was obtained for both samples; however, the Ks frequency distributions varied significantly. Importantly, we dated the WGD event at 51–53 mya, which coincides with the Paleocene-Eocene Thermal Maximum (PETM), when the Earth became warmer and wetter than any other period during the Cenozoic. Duplicate genes were preferentially retained for specific functions, such as environment response, further support that the duplicates may have promoted quick adaption to environmental changes and potentially resulted in evolutionary success, especially for pantropical species, such as C. thalictroides, which exhibits higher temperature tolerance.

Genome wide analysis of MADS-box gene family in Brassica oleracea reveals conservation and variation in flower development

BACKGROUND

MADS-box genes play important roles in vegetative growth and reproductive development and are essential for the correct development of plants (particularly inflorescences, flowers, and fruits). However, this gene family has not been identified nor their functions analyzed in Brassica oleracea.

RESULTS

In this study, we performed a whole-genome survey of the complete set of MADS-box genes in B. oleracea. In total, 91 MADS-box transcription factors (TFs) were identified and categorized as type I (Mα, Mβ, Mγ) and type II (MIKCC, MIKC*) groups according to the phylogeny and gene structure analysis. Among these genes, 59 were randomly distributed on 9 chromosomes, while the other 23 were assigned to 19 scaffolds and 9 genes from NCBI had no location information. Both RNA-sequencing and quantitative real-time-PCR analysis suggested that MIKC genes had more active and complex expression patterns than M type genes and most type II genes showed high flowering-related expression profiles. Additional quantitative real-time-PCR analysis of pedicel and four flower whorls revealed that the structure of the B.oleracea MIKC genes was conserved, but their homologues showed variable expression patterns compared to those in Arabidopsis thaliana.

CONCLUSION

This paper gives a detailed overview of the BolMADS genes and their expression patterns. The results obtained in this study provide useful information for understanding the molecular regulation of flower development and further functional characterization of MADS-box genes in B. oleracea.

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

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

MADS-box genes and crop domestication: the jack of all traits

MADS-box genes are key regulators of virtually every aspect of plant reproductive development. They play especially prominent roles in flowering time control, inflorescence architecture, floral organ identity determination, and seed development. The developmental and evolutionary importance of MADS-box genes is widely acknowledged. However, their role during flowering plant domestication is less well recognized. Here, we provide an overview illustrating that MADS-box genes have been important targets of selection during crop domestication and improvement. Numerous examples from a diversity of crop plants show that various developmental processes have been shaped by allelic variations in MADS-box genes. We propose that new genomic and genome editing resources provide an excellent starting point for further harnessing the potential of MADS-box genes to improve a variety of reproductive traits in crops. We also suggest that the biophysics of MADS-domain protein-protein and protein-DNA interactions, which is becoming increasingly well characterized, makes them especially suited to exploit coding sequence variations for targeted breeding approaches.

Gene Duplication and Transference of Function in the paleoAP3 Lineage of Floral Organ Identity Genes

The floral organ identity gene APETALA3 (AP3) is a MADS-box transcription factor involved in stamen and petal identity that belongs to the B-class of the ABC model of flower development. Thalictrum (Ranunculaceae), an emerging model in the non-core eudicots, has AP3 homologs derived from both ancient and recent gene duplications. Prior work has shown that petals have been lost repeatedly and independently in Ranunculaceae in correlation with the loss of a specific AP3 paralog, and Thalictrum represents one of these instances. The main goal of this study was to conduct a functional analysis of the three AP3 orthologs present in Thalictrum thalictroides, representing the paleoAP3 gene lineage, to determine the degree of redundancy versus divergence after gene duplication. Because Thalictrum lacks petals, and has lost the petal-specific AP3, we also asked whether heterotopic expression of the remaining AP3 genes contributes to the partial transference of petal function to the first whorl found in insect-pollinated species. To address these questions, we undertook functional characterization by virus-induced gene silencing (VIGS), protein-protein interaction and binding site analyses. Our results illustrate partial redundancy among Thalictrum AP3s, with deep conservation of B-class function in stamen identity and a novel role in ectopic petaloidy of sepals. Certain aspects of petal function of the lost AP3 locus have apparently been transferred to the other paralogs. A novel result is that the protein products interact not only with each other, but also as homodimers. Evidence presented here also suggests that expression of the different ThtAP3 paralogs is tightly integrated, with an apparent disruption of B function homeostasis upon silencing of one of the paralogs that codes for a truncated protein. To explain this result, we propose two testable alternative scenarios: that the truncated protein is a dominant negative mutant or that there is a compensational response as part of a back-up circuit. The evidence for promiscuous protein-protein interactions via yeast two-hybrid combined with the detection of AP3 specific binding motifs in all B-class gene promoters provide partial support for these hypotheses.

Genome-wide analysis of banana MADS-box family closely related to fruit development and ripening

Proteins encoded by MADS-box genes are important transcription factors involved in the regulation of flowering plant growth and development. Currently, no systematic information exists regarding the MADS-box family in the important tropical fruit banana. Ninety-six MADS-box genes were identified from the banana (Pahang) A genome. Phylogenetic analysis indicated that Musa acuminata MCM1-AGAMOUS- DEFICIENS-SRF (MaMADS) could be divided into MIKC^c, MIKC*, Mα/β and Mγ groups. MIKC^c could be further divided into 11 subfamilies, which was further supported by conserved motif and gene structure analyses. Transcriptome analysis on the Feng Jiao (FJ) and BaXi Jiao (BX) banana cultivars revealed that MaMADS genes are differentially expressed in various organs, at different fruit development and ripening stages, indicating the involvement of these genes in fruit development and ripening processes. Interactive network analysis indicated that MaMADS24 and 49 not only interacted with MaMADS proteins themselves, but also interacted with hormone-response proteins, ethylene signal transduction and biosynthesis-related proteins, starch biosynthesis proteins and metabolism-related proteins. This systematic analysis identified candidate MaMADS genes related to fruit development and ripening for further functional characterization in plants, and also provided new insights into the transcriptional regulation of MaMADS genes, facilitating the future genetic manipulation of MADS-mediated fruit development and ripening.

The phytochelatin synthase gene in date palm (Phoenix dactylifera L.): Phylogeny, evolution and expression

We studied date palm phytochelatin synthase type I (PdPCS1), which catalyzes the cytosolic synthesis of phytochelatins (PCs), a heavy metal binding protein, in plant cells. The gene encoding PdPCS1 (Pdpcs) consists of 8 exons and 7 introns and encodes a protein of 528 amino acids. PCs gene history was studied using Notung phylogeny. During evolution, gene loss from several lineages was predicted including Proteobacteria, Bilateria and Brassicaceae. In addition, eleven gene duplication events appeared toward interior nodes of the reconciled tree and four gene duplication events appeared toward the external nodes. These latter sequences belong to species with a second copy of PCs suggesting that this gene evolved through subfunctionalization. Pdpcs1 gene expression was measured in seedling hypocotyls exposed to Cd, Cu and Cr using quantitative real-time polymerase chain reaction (qPCR). A Pdpcs1 overexpression was evidenced in P. dactylifera seedlings exposed to metals suggesting that 1-the Pdpcs1 gene is functional, 2-there is an implication of the enzyme in metal detoxification mechanisms. Additionally, the structure of PdPCS1 was predicted using its homologue from Nostoc (cyanobacterium, NsPCS) as a template in Discovery studio and PyMol software. These analyses allowed us to identify the phytochelatin synthase type I enzyme in date palm (PdPCS1) via recognition of key consensus amino acids involved in the catalytic mechanism, and to propose a hypothetical binding and catalytic site for an additional substrate binding cavity.

Novel Insights into Tree Biology and Genome Evolution as Revealed Through Genomics

Reference genome sequences are the key to the discovery of genes and gene families that determine traits of interest. Recent progress in sequencing technologies has enabled a rapid increase in genome sequencing of tree species, allowing the dissection of complex characters of economic importance, such as fruit and wood quality and resistance to biotic and abiotic stresses. Although the number of reference genome sequences for trees lags behind those for other plant species, it is not too early to gain insight into the unique features that distinguish trees from nontree plants. Our review of the published data suggests that, although many gene families are conserved among herbaceous and tree species, some gene families, such as those involved in resistance to biotic and abiotic stresses and in the synthesis and transport of sugars, are often expanded in tree genomes. As the genomes of more tree species are sequenced, comparative genomics will further elucidate the complexity of tree genomes and how this relates to traits unique to trees.

Increased rates of protein evolution and asymmetric deceleration after the whole-genome duplication in yeasts

Background

Whole-genome duplication (WGD) events have shaped the genomes of eukaryotic organisms. Relaxed selection after duplication along with inherent functional constraints are thought to determine the fate of the paralogs and, ultimately, the evolution of gene function. Here, we investigated the rate of protein evolution (as measured by dN/dS ratios) before and after the WGD in the hemiascomycete yeasts, and the way in which changes in such rates relate to molecular and biological function.

Results

For most groups of orthologous genes (81%) we observed a change in the rates of evolution after genome duplication. Genes with atypically-low dN/dS ratio before the WGD were prone to increase their rates of evolution after duplication. Importantly, the paralogs were often different in their rates of evolution after the WGD (50% cases), however, this was more consistent with an asymmetric deceleration in the protein-evolution rates, rather than an asymmetric increase of the initial rates. Functional-category analysis showed that regulatory proteins such as protein kinases and transcription factors were enriched in genes that increase their rates of evolution after the WGD. While changes in the rate of protein-sequence evolution were associated to protein abundance, content of disordered regions, and contribution to fitness, these features were an attribute of specific functional classes.

Conclusions

Our results indicate that strong purifying selection in ancestral pre-duplication sequences is a strong predictor of increased rates after the duplication in yeasts and that asymmetry in evolution rate is established during the deceleration phase. In addition, changes in the rates at which paralogous sequences evolve before and after WGD are different for specific protein functions; increased rates of protein evolution after duplication occur preferentially in specific protein functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0895-1) contains supplementary material, which is available to authorized users.