The ABC model and the diversification of floral organ identity

STAMENLESS1 activates SUPERWOMAN 1 and FLORAL ORGAN NUMBER 1 to control floral organ identities and meristem fate in rice

Floral patterns are unique to rice and contribute significantly to its reproductive success. SL1 encodes a C2H2 transcription factor that plays a critical role in flower development in rice, but the molecular mechanism regulated by it remains poorly understood. Here, we describe interactions of the SL1 with floral homeotic genes, SPW1, and DL in specifying floral organ identities and floral meristem fate. First, the sl1 spw1 double mutant exhibited a stamen-to-pistil transition similar to that of sl1, spw1, suggesting that SL1 and SPW1 may located in the same pathway regulating stamen development. Expression analysis revealed that SL1 is located upstream of SPW1 to maintain its high level of expression and that SPW1, in turn, activates the B-class genes OsMADS2 and OsMADS4 to suppress DL expression indirectly. Secondly, sl1 dl displayed a severe loss of floral meristem determinacy and produced amorphous tissues in the third/fourth whorl. Expression analysis revealed that the meristem identity gene OSH1 was ectopically expressed in sl1 dl in the fourth whorl, suggesting that SL1 and DL synergistically terminate the floral meristem fate. Another meristem identity gene, FON1, was significantly decreased in expression in sl1 background mutants, suggesting that SL1 may directly activate its expression to regulate floral meristem fate. Finally, molecular evidence supported the direct genomic binding of SL1 to SPW1 and FON1 and the subsequent activation of their expression. In conclusion, we present a model to illustrate the roles of SL1, SPW1, and DL in floral organ specification and regulation of floral meristem fate in rice.

Exploring the patterns of evolution: Core thoughts and focus on the saltational model

The Modern Synthesis, a pillar in biological thought, united Darwin's species origin concepts with Mendel's laws of character heredity, providing a comprehensive understanding of evolution within species. Highlighting phenotypic variation and natural selection, it elucidated the environment's role as a selective force, shaping populations over time. This framework integrated additional mechanisms, including genetic drift, random mutations, and gene flow, predicting their cumulative effects on microevolution and the emergence of new species. Beyond the Modern Synthesis, the Extended Evolutionary Synthesis expands perspectives by recognizing the role of developmental plasticity, non-genetic inheritance, and epigenetics. We suggest that these aspects coexist in the plant evolutionary process; in this context, we focus on the saltational model, emphasizing how saltation events, such as dichotomous saltation, chromosomal mutations, epigenetic phenomena, and polyploidy, contribute to rapid evolutionary changes. The saltational model proposes that certain evolutionary changes, such as the rise of new species, may result suddenly from single macromutations rather than from gradual changes in DNA sequences and allele frequencies within a species over time. These events, observed in domesticated and wild higher plants, provide well-defined mechanistic bases, revealing their profound impact on plant diversity and rapid evolutionary events. Notably, next-generation sequencing exposes the likely crucial role of allopolyploidy and autopolyploidy (saltational events) in generating new plant species, each characterized by distinct chromosomal complements. In conclusion, through this review, we offer a thorough exploration of the ongoing dissertation on the saltational model, elucidating its implications for our understanding of plant evolutionary processes and paving the way for continued research in this intriguing field.

The regulation of chromatin configuration at AGAMOUS locus by LFR-SYD-containing complex is critical for reproductive organ development in Arabidopsis

Switch defective/sucrose non-fermentable (SWI/SNF) chromatin remodeling complexes are evolutionarily conserved, multi-subunit machinery that play vital roles in the regulation of gene expression by controlling nucleosome positioning and occupancy. However, little is known about the subunit composition of SPLAYED (SYD)-containing SWI/SNF complexes in plants. Here, we show that the Arabidopsis thaliana Leaf and Flower Related (LFR) is a subunit of SYD-containing SWI/SNF complexes. LFR interacts directly with multiple SWI/SNF subunits, including the catalytic ATPase subunit SYD, in vitro and in vivo. Phenotypic analyses of lfr-2 mutant flowers revealed that LFR is important for proper filament and pistil development, resembling the function of SYD. Transcriptome profiling revealed that LFR and SYD shared a subset of co-regulated genes. We further demonstrate that the LFR and SYD interdependently activate the transcription of AGAMOUS (AG), a C-class floral organ identity gene, by regulating the occupation of nucleosome, chromatin loop, histone modification, and Pol II enrichment on the AG locus. Furthermore, the chromosome conformation capture (3C) assay revealed that the gene loop at AG locus is negatively correlated with the AG expression level, and LFR-SYD was functional to demolish the AG chromatin loop to promote its transcription. Collectively, these results provide insight into the molecular mechanism of the Arabidopsis SYD-SWI/SNF complex in the control of higher chromatin conformation of the floral identity gene essential to plant reproductive organ development.

AtSAMS regulates floral organ development by DNA methylation and ethylene signaling pathway

S-adenosylmethionine synthase is the key enzyme involved in the biosynthesis of S-adenosylmethionine, which serves as the universal methyl group donor and a common precursor for the biosynthesis of ethylene and polyamines. However, little is known about how SAMS controls plant development. Here, we report that the abnormal floral organ development in the AtSAMS-overexpressing plants is caused by DNA demethylation and ethylene signaling. The whole-genome DNA methylation level decreased, and ethylene content increased in SAMOE. Wild-type plants treated with DNA methylation inhibitor mimicked the phenotypes and the ethylene levels in SAMOE, suggesting that DNA demethylation enhanced ethylene biosynthesis, which led to abnormal floral organ development. DNA demethylation and elevated ethylene resulted in changes in the expression of ABCE genes, which is essential for floral organ development. Furthermore, the transcript levels of ACE genes were highly correlated to their methylation levels, except for the down-regulation of the B gene, which might have resulted from demethylation-independent ethylene signaling. SAMS-mediated methylation and ethylene signaling might create crosstalk in the process of floral organ development. Together, we provide evidence that AtSAMS regulates floral organ development by DNA methylation and ethylene signaling pathway.

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.

All's well that ends well: the timing of floral meristem termination

Floral meristem termination (FMT) represents one of the defining features of a floral meristem relative to a vegetative meristem. Timing of FMT is a major determinant of the total number of organs in a flower, and canalization toward relatively rapid FMT is considered to have been a major force in shaping angiosperm evolution. For decades, investigation of FMT has been focused on model systems that only produce four whorls of organs in a flower, while little is known about the molecular basis that underlies nature variation in the timing of FMT. Here, we hypothesize on how known pathways could have been modified to generate variation in FMT and explain how developing new model systems will help to deepen our understanding of the genetic control and evolution of FMT.

Genome-Wide Identification and Analysis of the MADS-Box Transcription Factor Genes in Blueberry (Vaccinium spp.) and Their Expression Pattern during Fruit Ripening

MADS-box is a class of transcriptional regulators that are ubiquitous in plants and plays important roles in the process of plant growth and development. Identification and analysis of blueberry MADS-box genes can lay a foundation for their function investigations. In the present study, 249 putative MADS-box genes were identified in the blueberry genome. Those MADS-box genes were distributed on 47 out of 48 chromosomes. The phylogenetic and evolutionary analyses showed that blueberry MADS-box genes were divided into 131 type I members and 118 type II members. The type I genes contained an average of 1.89 exons and the type II genes contained an average of 7.83 exons. Motif analysis identified 15 conserved motifs, of which 4 were related to the MADS domain and 3 were related to the K-box domain. A variety of cis-acting elements were found in the promoter region of the blueberry MADS-box gene, indicating that the MADS-box gene responded to various hormones and environmental alterations. A total of 243 collinear gene pairs were identified, most of which had a Ka/Ks value of less than 1. Nine genes belonging to SEP, AP3/PI, and AGL6 subfamilies were screened based on transcriptomic data. The expression patterns of those nine genes were also verified using quantitative PCR, suggesting that VcMADS6, VcMADS35, VcMADS44, VcMADS58, VcMADS125, VcMADS188, and VcMADS212 had potential functions in blueberry fruit ripening. The results of this study provide references for an in-depth understanding of the biological function of the blueberry MADS-box genes and the mechanism of blueberry fruit ripening.

Identification and Molecular Characterization of RWP-RK Transcription Factors in Soybean

The RWP-RK is a small family of plant-specific transcription factors that are mainly involved in nitrate starvation responses, gametogenesis, and root nodulation. To date, the molecular mechanisms underpinning nitrate-regulated gene expression in many plant species have been extensively studied. However, the regulation of nodulation-specific NIN proteins during nodulation and rhizobial infection under nitrogen starvation in soybean still remain unclear. Here, we investigated the genome-wide identification of RWP-RK transcription factors and their essential role in nitrate-inducible and stress-responsive gene expression in soybean. In total, 28 RWP-RK genes were identified from the soybean genome, which were unevenly distributed on 20 chromosomes from 5 distinct groups during phylogeny classification. The conserved topology of RWP-RK protein motifs, cis-acting elements, and functional annotation has led to their potential as key regulators during plant growth, development, and diverse stress responses. The RNA-seq data revealed that the up-regulation of GmRWP-RK genes in the nodules indicated that these genes might play crucial roles during root nodulation in soybean. Furthermore, qRT-PCR analysis revealed that most GmRWP-RK genes under Phytophthora sojae infection and diverse environmental conditions (such as heat, nitrogen, and salt) were significantly induced, thus opening a new window of possibilities into their regulatory roles in adaptation mechanisms that allow soybean to tolerate biotic and abiotic stress. In addition, the dual luciferase assay indicated that GmRWP-RK1 and GmRWP-RK2 efficiently bind to the promoters of GmYUC2, GmSPL9, and GmNIN, highlighting their possible involvement in nodule formation. Together, our findings provide novel insights into the functional role of the RWP-RK family during defense responses and root nodulation in soybean.

Genome-Wide Identification of MADS-Box Family Genes in Safflower (Carthamus tinctorius L.) and Functional Analysis of CtMADS24 during Flowering

Safflower is an important economic crop with a plethora of industrial and medicinal applications around the world. The bioactive components of safflower petals are known to have pharmacological activity that promotes blood circulation and reduces blood stasis. However, fine-tuning the genetic mechanism of flower development in safflower is still required. In this study, we report the genome-wide identification of MADS-box transcription factors in safflower and the functional characterization of a putative CtMADS24 during vegetative and reproductive growth. In total, 77 members of MADS-box-encoding genes were identified from the safflower genome. The phylogenetic analysis divided CtMADS genes into two types and 15 subfamilies. Similarly, bioinformatic analysis, such as of conserved protein motifs, gene structures, and cis-regulatory elements, also revealed structural conservation of MADS-box genes in safflower. Furthermore, the differential expression pattern of CtMADS genes by RNA-seq data indicated that type II genes might play important regulatory roles in floral development. Similarly, the qRT-PCR analysis also revealed the transcript abundance of 12 CtMADS genes exhibiting tissue-specific expression in different flower organs. The nucleus-localized CtMADS24 of the AP1 subfamily was validated by transient transformation in tobacco using GFP translational fusion. Moreover, CtMADS24-overexpressed transgenic Arabidopsis exhibited early flowering and an abnormal phenotype, suggesting that CtMADS24 mediated the expression of genes involved in floral organ development. Taken together, these findings provide valuable information on the regulatory role of CtMADS24 during flower development in safflower and for the selection of important genes for future molecular breeding programs.

Transcriptomic heterochrony and completely cleistogamous flower development in the mycoheterotrophic orchid Gastrodia

Cleistogamy, in which plants can reproduce via self-fertilization within permanently closed flowers, has evolved in > 30 angiosperm lineages; however, consistent with Darwin's doubts about its existence, complete cleistogamy - the production of only cleistogamous flowers - has rarely been recognized. Thus far, the achlorophyllous orchid genus, Gastrodia, is the only known genus with several plausible completely cleistogamous species. Here, we analyzed the floral developmental transcriptomes of two recently evolved, completely cleistogamous Gastrodia species and their chasmogamous sister species to elucidate the possible changes involved in producing common cleistogamous traits. The ABBA-BABA test did not support introgression and protein sequence convergence as evolutionary mechanisms leading to cleistogamy, leaving convergence in gene expression as a plausible mechanism. Regarding transcriptomic differentiation, the two cleistogamous species had common modifications in the expression of developmental regulators, exhibiting a gene family-wide signature of convergent expression changes in MADS-box genes. Our transcriptomic pseudotime analysis revealed a prolonged juvenile state and eventual maturation, a heterochronic pattern consistent with partial neoteny, in cleistogamous flower development. These findings indicate that transcriptomic partial neoteny, arising from changes in the expression of conserved developmental regulators, might have contributed to the rapid and repeated evolution of cleistogamous flowers in Gastrodia.

Overexpression of the AGL42 gene in cotton delayed leaf senescence through downregulation of NAC transcription factors

Premature leaf senescence negatively influences the physiology and yield of cotton plants. The conserved IDLNL sequence in the C-terminal region of AGL42 MADS-box determines its repressor potential for the down regulation of senescence-related genes. To determine the delay in premature leaf senescence, Arabidopsis AGL42 gene was overexpressed in cotton plants. The absolute quantification of transgenic cotton plants revealed higher mRNA expression of AGL42 compared to that of the non-transgenic control. The spatial expression of GUS fused with AGL42 and the mRNA level was highest in the petals, abscission zone (flower and bud), 8 days post anthesis (DPA) fiber, fresh mature leaves, and senescenced leaves. The mRNA levels of different NAC senescence-promoting genes were significantly downregulated in AGL42 transgenic cotton lines than those in the non-transgenic control. The photosynthetic rate and chlorophyll content were higher in AGL42 transgenic cotton lines than those in the non-transgenic control. Fluorescence in situ hybridization of the AG3 transgenic cotton line revealed a fluorescent signal on chromosome 1 in the hemizygous form. Moreover, the average number of bolls in the transgenic cotton lines was significantly higher than that in the non-transgenic control because of the higher retention of floral buds and squares, which has the potential to improve cotton fiber yield.

Genetic architecture underlying variation in floral meristem termination in Aquilegia

Floral organs are produced by floral meristems (FMs), which harbor stem cells in their centers. Since each flower only has a finite number of organs, the stem cell activity of an FM will always terminate at a specific time point, a process termed floral meristem termination (FMT). Variation in the timing of FMT can give rise to floral morphological diversity, but how this process is fine-tuned at a developmental and evolutionary level is poorly understood. Flowers from the genus Aquilegia share identical floral organ arrangement except for stamen whorl number (SWN), making Aquilegia a well-suited system for investigation of this process: differences in SWN between species represent differences in the timing of FMT. By crossing A. canadensis and A. brevistyla, quantitative trait locus (QTL) mapping has revealed a complex genetic architecture with seven QTL. We explored potential candidate genes under each QTL and characterized novel expression patterns of select loci of interest using in situ hybridization. To our knowledge, this is the first attempt to dissect the genetic basis of how natural variation in the timing of FMT is regulated, and our results provide insight into how floral morphological diversity can be generated at the meristematic level.

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 association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

Functional identification of the different regions in B-class floral homeotic MADS-box proteins IiAP3 and IiPI from Isatis indigotica

APETALA3 (AP3) and PISTILLATA (PI) are B-class MADS-box floral homeotic genes of Arabidopsis and are involved in specifying the identity of petals and stamens. In the present work, IiAP3 and IiPI, the respective orthologous genes of AP3 and PI, were cloned from Isatis indigotica. By expressing in ap3-6 and pi-1 homozygous mutant and in wild-type Arabidopsis under the control of AP3 promoter or CaMV 35S promoter, we demonstrated that IiAP3 and IiPI were functionally equivalent to AP3 and PI of Arabidopsis. Referring to previous reports and the research results in the present work, expression patterns of AP3 and PI homologs are not the same in different angiosperms possessing diverse floral structures. It suggests that the alterations in expression may contribute to the changing morphology of flowers. To further determine the relationship between IiAP3 and IiPI, the coding sequences of the different structural regions in these two proteins were swapped with each other, and the data collected from transgenic Arabidopsis plants of the chimeric constructs suggested that MADS domain was irreplaceable for the function of IiAP3, K domain of IiAP3 was involved in specifying the identity of stamens, K domain of IiPI was mainly related to the formation of petals, and C-terminal region of IiPI was involved in characterization of stamens. In addition, a complete KC region of these two proteins was more effective in phenotypic complementation of the mutants.

Conjunctive Analyses of Bulk Segregant Analysis Sequencing and Bulk Segregant RNA Sequencing to Identify Candidate Genes Controlling Spikelet Sterility of Foxtail Millet

Foxtail millet has gradually become a model gramineous C4 crop owing to its short growth period and small genome. Research on the development of its spikelets is not only directly related to the yield and economic value of foxtail millet but also can provide a reference for studying the fertility of other C4 crops. In this study, a hybrid population containing 200 offspring was constructed from the Xinong8852 and An15 parental lines, and two extreme trait populations were constructed from the F2 generation for analysis of the spikelet sterility. The F2 population conformed to a 3:1 Mendelian segregation ratio, and it was thus concluded that this trait is likely controlled by a single recessive gene. Bulk segregant analysis sequencing (BSA-Seq) was used to determine the candidate regions and candidate genes related to the development of foxtail millet spikelets. Additionally, the functional analysis of differentially expressed genes in populations with different traits was conducted by bulk segregant RNA sequencing (BSR-Seq). Finally, conjunctive analysis of BSA-Seq and BSR-Seq results, combined with biological information analysis, revealed six genes on chromosome VII that were ultimately identified as candidate genes controlling foxtail millet spikelet development. This study provides a new reference for research on foxtail millet sterility and lays a solid foundation for the examination of fertility in other gramineous crops.

Reflections on the Triptych of Meristems That Build Flowering Branches in Tomato

Branching is an important component determining crop yield. In tomato, the sympodial pattern of shoot and inflorescence branching is initiated at floral transition and involves the precise regulation of three very close meristems: (i) the shoot apical meristem (SAM) that undergoes the first transition to flower meristem (FM) fate, (ii) the inflorescence sympodial meristem (SIM) that emerges on its flank and remains transiently indeterminate to continue flower initiation, and (iii) the shoot sympodial meristem (SYM), which is initiated at the axil of the youngest leaf primordium and takes over shoot growth before forming itself the next inflorescence. The proper fate of each type of meristems involves the spatiotemporal regulation of FM genes, since they all eventually terminate in a flower, but also the transient repression of other fates since conversions are observed in different mutants. In this paper, we summarize the current knowledge about the genetic determinants of meristem fate in tomato and share the reflections that led us to identify sepal and flower abscission zone initiation as a critical stage of FM development that affects the branching of the inflorescence.

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

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

The underlying molecular conservation and diversification of dioecious flower and leaf buds provide insights into the development, dormancy breaking, flowering, and sex association of willows

As harbingers of bursting growth, flower buds and leaf buds generally show similar surface morphologies but different structural and functional changes. Dioecious plants further generate four types of Female/Male Flower/Leaf Buds (FFB, FLB, MFB, and MLB), showing a complex regulation. However, little is known about their underlying molecular mechanisms. Here, we exemplify the woody dioecious Salix linearistipularis to investigate their morphological characteristics and potential molecular mechanisms by combining cytological, physiological, phenological, and transcriptomic datasets. First, FFB and MFB have simultaneous development dynamics and so do FLB and MLB. Interestingly, FLB and MLB show very similar expression profiles preparing for photosynthesis and stress-tolerance, whereas FFB and MFB show great similarities but also striking sexual differences. Comparing flower buds and leaf buds after their revival from dormancy shows different cold- and vernalization-responsive genes (e.g. SliVRN1, SliAGL19, and SliAGL24), implying different programming processes for dormancy breaking between the buds. Moreover, except SliAP3, the expression of ABCDE model genes is consistent with their roles in the buds, suggesting a conserved mechanism of flower development between dioecious Salix and hermaphrodite Arabidopsis. Finally, considering sex-associated genes (e.g. SliCLE25, SliTPS21, and SliARR9) on Salix chromosomes and other reports, we hypothesize a dynamic model of sex determination on chromosomes 15 and 19 in the last ancestor of Salix and Populus but evolutionarily on 15 in Salix after their divergence. Together, our study provides new insights into the molecular mechanisms of dioecious four-type buds by showing the genes involved in their development, dormancy breaking, flowering, and sexual association.

An Optimized Transformation System and Functional Test of CYC-Like TCP Gene CpCYC in Chirita pumila (Gesneriaceae)

The development of an ideal model plant located at a key phylogenetic node is critically important to advance functional and regulatory studies of key regulatory genes in the evolutionary developmental (evo-devo) biology field. In this study, we selected Chirita pumila in the family Gesneriaceae, a basal group in Lamiales, as a model plant to optimize its genetic transformation system established previously by us through investigating a series of factors and further conduct functional test of the CYC-like floral symmetry gene CpCYC. By transforming a RNAi:CpCYC vector, we successfully achieved the desired phenotypes of upright actinomorphic flowers, which suggest that CpCYC actually determines the establishment of floral zygomorphy and the horizontal orientation of flowers in C. pumila. We also confirmed the activities of CpCYC promoter in dorsal petals, dorsal/lateral staminodes, as well as the pedicel by transferring a CpCYC promoter:GUS vector into C. pumila. Furthermore, we testified the availability of a transient gene expression system using C. pumila mesophyll protoplasts. The improved transformation system together with the inherent biological features would make C. pumila an attractive new model in functional and regulatory studies for a broad range of evo-devo issues.

Distance-based measurement determines the coexistence of B protein hetero- and homodimers in lily tepal and stamen tetrameric complexes

The floral quartet model proposes that plant MADS box proteins function as higher order tetrameric complexes. However, in planta evidence for MADS box tetramers remains scarce. Here, we applied a strategy using in vivo fluorescence resonance energy transfer (FRET) based on the distance change and distance symmetry of stable tetrameric complexes in tobacco (Nicotiana benthamiana) leaf cells to improve the accuracy of the estimation of heterotetrameric complex formation. This measuring system precisely verified the stable state of Arabidopsis petal (AP3/PI/SEP3/AP1) and stamen (AP3/PI/SEP3/AG) complexes and showed that the lily (Lilium longiflorum) PI co-orthologs LMADS8 and LMADS9 likely formed heterotetrameric petal complexes with Arabidopsis AP3/SEP3/AP1, which rescued petal defects of pi mutants. However, L8/L9 did not form heterotetrameric stamen complexes with Arabidopsis AP3/SEP3/AG to rescue the stamen defects of the pi mutants. Importantly, this system was applied successfully to find complicated tepal and stamen heterotetrameric complexes in lily. We found that heterodimers of B function AP3/PI orthologs (L1/L8) likely coexist with the homodimers of PI orthologs (L8/L8, L9/L9) to form five (two most stable and three stable) tepal- and four (one most stable and three stable) stamen-related heterotetrameric complexes with A/E and C/E function proteins in lily. Among these combinations, L1 preferentially interacted with L8 to form the most stable heterotetrameric complexes, and the importance of the L8/L8 and L9/L9 homodimers in tepal/stamen formation in lily likely decreased to a minor part during evolution. The system provides substantial improvements for successfully estimating the existence of unknown tetrameric complexes in plants.

Multifunctional evolution of B and AGL6 MADS box genes in orchids

We previously found that B and AGL6 proteins form L (OAP3-2/OAGL6-2/OPI) and SP (OAP3-1/OAGL6-1/OPI) complexes to determine lip/sepal/petal identities in orchids. Here, we show that the functional L' (OAP3-1/OAGL6-2/OPI) and SP' (OAP3-2/OAGL6-1/OPI) complexes likely exist and AP3/PI/AGL6 genes have acquired additional functions during evolution. We demonstrate that the presumed L' complex changes the structure of the lower lateral sepals and helps the lips fit properly in the center of the flower. In addition, we find that OAP3-1/OAGL6-1/OPI in SP along with presumed SP' complexes regulate anthocyanin accumulation and pigmentation, whereas presumed L' along with OAP3-2/OAGL6-2/OPI in L complexes promotes red spot formation in the perianth. Furthermore, the B functional proteins OAP3-1/OPI and OAGL6-1 in the SP complex could function separately to suppress sepal/petal senescence and promote pedicel abscission, respectively. These findings expand the current knowledge behind the multifunctional evolution of the B and AGL6 genes in plants.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Comparative RNA-Seq analysis reveals genes associated with masculinization in female Cannabis sativa

Using RNA profiling, we identified several silver thiosulfate-induced genes that potentially control the masculinization of female Cannabis sativa plants. Genetically female Cannabis sativa plants normally bear female flowers, but can develop male flowers in response to environmental and developmental cues. In an attempt to elucidate the molecular elements responsible for sex expression in C. sativa plants, we developed genetically female lines producing both female and chemically-induced male flowers. Furthermore, we carried out RNA-Seq assays aimed at identifying differentially expressed genes responsible for male flower development in female plants. The results revealed over 10,500 differentially expressed genes, of which around 200 potentially control masculinization of female cannabis plants. These genes include transcription factors and other genes involved in male organ (i.e., anther and pollen) development, as well as genes involved in phytohormone signalling and male-biased phenotypes. The expressions of 15 of these genes were further validated by qPCR assay confirming similar expression patterns to that of RNA-Seq data. These genes would be useful for understanding predisposed plants producing flowers of both sex types in the same plant, and help breeders to regulate the masculinization of female plants through targeted breeding and plant biotechnology.

Comparative Transcriptome Reveals the Genes' Adaption to Herkogamy of Lumnitzera littorea (Jack) Voigt

Lumnitzera littorea (Jack) Voigt is among the most endangered mangrove species in China. The morphology and evolution of L. littorea flowers have received substantial attention for their crucial reproductive functions. However, little is known about the genomic regulation of flower development in L. littorea. In this study, we characterized the morphology of two kinds of L. littorea flowers and performed comparative analyses of transcriptome profiles of the two different flowers. Morphological observation showed that some flowers have a column embedded in the petals while others produce a stretched flower style during petal unfolding in flowering. By using RNA-seq, we obtained 138,857 transcripts that were assembled into 82,833 unigenes with a mean length of 1055.48 bp. 82,834 and 34,997 unigenes were assigned to 52 gene ontology (GO) functional groups and 364 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. A total of 4,267 differentially expressed genes (DEGs), including 1,794 transcription factors (TFs), were identified between two types of flowers. These TFs are mainly involved in bHLH, B3, bZIP, MYB-related, and NAC family members. We further validated that 12 MADS-box genes, including 4 MIKC-type and 8 M-type TFs, were associated with the pollinate of L. littorea by herkogamy. Our current results provide valuable information for genetic analysis of L. littorea flowering and may be useful for illuminating its adaptive evolutionary mechanisms.

Evolutionary Variation in MADS Box Dimerization Affects Floral Development and Protein Abundance in Maize

Interactions between MADS box transcription factors are critical in the regulation of floral development, and shifting MADS box protein-protein interactions are predicted to have influenced floral evolution. However, precisely how evolutionary variation in protein-protein interactions affects MADS box protein function remains unknown. To assess the impact of changing MADS box protein-protein interactions on transcription factor function, we turned to the grasses, where interactions between B-class MADS box proteins vary. We tested the functional consequences of this evolutionary variability using maize (Zea mays) as an experimental system. We found that differential B-class dimerization was associated with subtle, quantitative differences in stamen shape. In contrast, differential dimerization resulted in large-scale changes to downstream gene expression. Differential dimerization also affected B-class complex composition and abundance, independent of transcript levels. This indicates that differential B-class dimerization affects protein degradation, revealing an important consequence for evolutionary variability in MADS box interactions. Our results highlight complexity in the evolution of developmental gene networks: changing protein-protein interactions could affect not only the composition of transcription factor complexes but also their degradation and persistence in developing flowers. Our results also show how coding change in a pleiotropic master regulator could have small, quantitative effects on development.

Comparative transcriptome analysis of flower bud transition and functional characterization of EjAGL17 involved in regulating floral initiation in loquat

Floral initiation plays a critical role for reproductive success in plants, especially fruit trees. However, little information is known on the mechanism of the initiation in loquat (Eriobotrya japonica Lindl.). Here, we used transcriptomic, expression and functional analysis to investigate the candidate genes in floral initiation in loquat. Comparative transcriptome analysis showed differentially expressed genes (DEGs) were mainly enriched in the metabolic pathways of plant hormone signal transduction. The DEGs were mainly involved in the gibberellin, auxin, cytokinin, abscisic acid, salicylic acid and ethylene signaling pathways. Meanwhile, some transcription factors, including MADS-box (MCM1, AGAMOUS, DEFICIENS and SRF), MYB (Myeloblastosis), TCP (TEOSINTE BRANCHED 1, CYCLOIDEA and PCF1), WOX (WUSCHEL-related homeobox) and WRKY (WRKY DNA-binding protein), were significantly differentially expressed. Among these key DEGs, we confirmed that an AGL17 ortholog EjAGL17 was significantly upregulated at the flower bud transition stage. Phylogenetic tree analysis revealed that EjAGL17 was grouped into an AGL17 clade of MADS-box transcription factors. Protein sequence alignment showed that EjAGL17 included a distinctive C-terminal domain. Subcellular localization of EjAGL17 was found only in the nucleus. Expression levels of EjAGL17 reached the highest at the development stage of flower bud transition. Moreover, ectopic expression of EjAGL17 in Arabidopsis significantly exhibited early flowering. Our study provides abundant resources of candidate genes for studying the mechanisms underlying the floral initiation in loquat and other Rosaceae species.

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.

Identification of the target genes of AqAPETALA3-3 (AqAP3-3) in Aquilegia coerulea (Ranunculaceae) helps understand the molecular bases of the conserved and nonconserved features of petals

Identification and comparison of the conserved and variable downstream genes of floral organ identity regulators are critical to understanding the mechanisms underlying the commonalities and peculiarities of floral organs. Yet, because of the lack of studies in nonmodel species, a general picture of the regulatory evolution between floral organ identity genes and their targets is still lacking. Here, by conducting extensive chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), electrophoretic mobility shift assay and bioinformatic analyses, we identify and predict the target genes of a petal identity gene, AqAPETALA3-3 (AqAP3-3), in Aquilegia coerulea (Ranunculaceae) and compare them with those of its counterpart in Arabidopsis thaliana, AP3. In total, 7049 direct target genes are identified for AqAP3-3, of which 2394 are highly confident and 1085 are shared with AP3. Gene Ontology enrichment analyses further indicate that conserved targets are largely involved in the formation of identity-related features, whereas nonconserved targets are mostly required for the formation of species-specific features. These results not only help understand the molecular bases of the conserved and nonconserved features of petals, but also pave the way to studying the regulatory evolution between floral organ identity genes and their targets.

Developmental and molecular characterization of novel staminodes in Aquilegia

BACKGROUND AND AIMS

The ranunculid model system Aquilegia is notable for the presence of a fifth type of floral organ, the staminode, which appears to be the result of sterilization and modification of the two innermost whorls of stamens. Previous studies have found that the genetic basis for the identity of this new organ is the result of sub- and neofunctionalization of floral organ identity gene paralogues; however, we do not know the extent of developmental and molecular divergence between stamens and staminodes.

METHODS

We used histological techniques to describe the development of the Aquilegia coerulea 'Origami' staminode relative to the stamen filament. These results have been compared with four other Aquilegia species and the closely related genera Urophysa and Semiaquilegia. As a complement, RNA sequencing has been conducted at two developmental stages to investigate the molecular divergence of the stamen filaments and staminodes in A. coerulea 'Origami'.

KEY RESULTS

Our developmental study has revealed novel features of staminode development, most notably a physical interaction along the lateral margin of adjacent organs that appears to mediate their adhesion. In addition, patterns of abaxial/adaxial differentiation are observed in staminodes but not stamen filaments, including asymmetric lignification of the adaxial epidermis in the staminodes. The comparative transcriptomics are consistent with the observed lignification of staminodes and indicate that stamen filaments are radialized due to overexpression of adaxial identity, while the staminodes are expanded due to the balanced presence of abaxial identity.

CONCLUSIONS

These findings suggest a model in which the novel staminode identity programme interacts with the abaxial/adaxial identity pathways to produce two whorls of laterally expanded organs that are highly differentiated along their abaxial/adaxial axis. While the ecological function of Aquilegia staminodes remains to be determined, these data are consistent with a role in protecting the early carpels from herbivory and/or pathogens.

Draft genome of the famous ornamental plant Paeonia suffruticosa

Tree peony (Paeonia Sect. Moutan) is a famous ornamental plant, with huge historical, cultural, and economic significance worldwide. In this study, we reported the ~13.79 Gb draft genome of a wide-grown Paeonia suffruticosa cultivar "Luo shen xiao chun," representing the largest sequenced genome in dicots to date. Phylogenetic analyses based on genome sequences demonstrated that P. suffruticosa was placed as sister to Vitales, and they together formed a clade that was sister to Rosids, weakly supporting a relationship of ((Saxifragales and Vitales) and Rosids). The identification and expression analysis of MADS-box genes based on the genome assembly and de novo transcriptome assembly of P. suffruticosa revealed that the function of C class genes was restricted in flower development, which might be responsible for the stamen petalody in tree peony cultivars. Overall, the first sequenced genome in the family Paeoniaceae provides an important resource for the origin, domestication, and evolutionary study as well as cultivar breeding in tree peony.

Homeotic transformation from stamen to petal in Eriobotrya japonica is associated with hormone signal transduction and reduction of the transcriptional activity of EjAG

Double-flower loquat (Eriobotrya japonica) is a new germplasm with homeotic transformation of stamen into petal in whorl 3. However, little information is available on the molecular mechanism of this transformation. Herein, we analyzed the transcriptome, candidate genes and endogenous hormones to investigate the mechanisms underlying this homeotic transformation. Some transcription factors, such as MADS-box, TCP and MYB, were significantly differentially expressed. Importantly, we confirmed that one of these (DN39625_c0_g1), which encoded a C-class floral homeotic protein referred to as AGAMOUS ortholog (EjAG), was significantly downregulated. Subcellular localization of EjAG was found to be in the nucleus. Ectopic expression of EjAG rescued the development of stamens and carpels from the double-flower phenotype in an Arabidopsis ag mutant, suggesting that EjAG expression is associated with double-flower formation. Meanwhile, enrichment analyses showed that the differentially expressed genes (DEGs) were mainly involved in the metabolic pathways of hormone signal transduction. The DEGs of auxin, gibberellin A (GA) and cytokinin signaling pathways were mainly upregulated. However, the DEGs of abscisic acid (ABA) and the ethylene signaling pathway were mainly downregulated. Accordingly, the concentrations of indoleacetic acid, kinetin and GA₃ were high at the petaloid stamen stage, but the ABA concentration remained low. The identified genes and pathways provide abundant sequence resources for studying the mechanisms underlying the homeotic transformation in loquat and other Rosaceae species.

Deep Sequencing and Analysis of Transcriptomes of Pinus koraiensis Sieb. & Zucc

The objective of this research was to study the differences in endogenous hormone levels and the genes related to reproductive development in Chinese pinenut (Pinus koraiensis) trees of different ages. The apical buds of P. koraiensis were collected from 2-, 5-, 10-, 15-, and 30-year-old plants and also from grafted plants. There were three replicates from each group used for transcriptome sequencing. After assembly and annotation, we identified the differentially expressed genes (DEGs) and performed enrichment analysis, pathway analysis, and expression analysis of the DEGs in each sample. The results showed that unigenes related to reproductive development, such as c64070.graph_c0 and c68641.graph_c0, were expressed at relatively low levels at young ages, and that the relative expression gradually increased with increasing plant age. In addition the highest expression levels were reached around 10 and 15 years of age, after which they gradually decreased. Moreover, some unigenes, such as c61855.graph_c0, were annotated as abscisic acid hydroxylase genes, and the expression of c61855.graph_c0 gradually declined with increasing age in P. koraiensis.

Mass-spectrometry-based draft of the Arabidopsis proteome

Plants are essential for life and are extremely diverse organisms with unique molecular capabilities¹. Here we present a quantitative atlas of the transcriptomes, proteomes and phosphoproteomes of 30 tissues of the model plant Arabidopsis thaliana. Our analysis provides initial answers to how many genes exist as proteins (more than 18,000), where they are expressed, in which approximate quantities (a dynamic range of more than six orders of magnitude) and to what extent they are phosphorylated (over 43,000 sites). We present examples of how the data may be used, such as to discover proteins that are translated from short open-reading frames, to uncover sequence motifs that are involved in the regulation of protein production, and to identify tissue-specific protein complexes or phosphorylation-mediated signalling events. Interactive access to this resource for the plant community is provided by the ProteomicsDB and ATHENA databases, which include powerful bioinformatics tools to explore and characterize Arabidopsis proteins, their modifications and interactions.

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.

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.

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.

Default Sex and Single Gene Sex Determination in Dioecious Plants

A well-established hypothesis for the evolution of dioecy involves two genes linked at a sex-determining region (SDR). Recently there has been increased interest in possible single gene sex determination. Work in Populus has finally provided direct experimental evidence for single gene sex determination in plants using CRISPR-Cas9 to knock out a single gene and convert individuals from female to male. In poplar, the feminizing factor popARR17 acts as a "master regulator", analogous to the mammalian masculinizing factor SRY. The production of fully functional males from females by a simple single gene knockout is experimental evidence that an antagonistic male-determining factor does not exist in Populus. Mammals have a "default sex" (female), as do poplar trees (Populus), although the default sex in poplars is male. The occurrence of single gene sex determination with a default sex may be much commoner in plants than hitherto expected, especially when dioecy evolved via monoecy. The master regulator does not even need to be at the SDR (although it may be). In most poplars the feminizing factor popARR17 is not at the SDR, but instead a negative regulator of it. So far there is little information on how high-level regulators are connected to floral phenotype. A model is presented of how sex-determining genes could lead to different floral morphologies via MADS-box floral developmental genes.

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.

Ectopic expression of an Eriobotrya japonica APETALA3 ortholog rescues the petal and stamen identities in Arabidopsis ap3-3 mutant

APETALA3: (AP3) encodes a floral homeotic class B-function MADS-box protein and plays crucial roles in petal and stamen development. To better understand the functional roles of AP3 orthologs in Eriobotrya, we isolated and identified an AP3 ortholog, referred to as EjAP3, from Eriobotrya japonica. Analyses of protein sequence and phylogenetic tree showed that the EjAP3 was assigned to the rosids euAP3 lineage and included a distinctive PI-derived and euAP3 motifs at the C-terminal domain. Subcellular localization of EjAP3 was determined to be in the nucleus. Expression analysis suggested that EjAP3 expression was restricted only in petals and stamens, but not in sepals and carpels. Importantly, during the floral development, EjAP3 expression level was the highest at the stage of visible floral bud. Furthermore, ectopic expression of EjAP3 in Arabidopsis ap3-3 mutant rescued the second whorl petals and the third whorl stamens. The expression pattern and function characterization of EjAP3 contribute to better understand the roles of AP3 orthologs in Eriobotrya.

OsLUGL is involved in the regulating auxin level and OsARFs expression in rice (Oryza sativa L.)

Specification of floral organ identity is critical for floral morphology and inflorescence architecture. Floral organ identity in plants is controlled by floral homeotic A/B/C/D/E-class genes. Although multiple genes regulate floral organogenesis, our understanding of the regulatory network remains fragmentary. Here, we characterized a rice floral organ gene KAIKOUXIAO (KKX), mutation of which produces an uncharacteristic open hull, abnormal seed and semi-sterility. KKX encodes a putative LEUNIG-like (LUGL) transcriptional regulator OsLUGL. OsLUGL is preferentially expressed in young panicles and its protein can interact with OsSEU, which functions were reported as an adaptor for LEUNIG. OsLUGL-OsSEU functions together as a transcriptional co-regulatory complex to control organ identity. SEP3 (such as OsMADS8) and AP1 (such as OsMADS18) serve as the DNA-binding partner of OsLUGL-OsSEU complex. Further studies indicated that OsMADS8 and OsMADS18 could bind to the promoter of OsGH3-8. The altered expression of OsGH3-8 might cause the increased auxin level and the decreased expression of OsARFs. Overall, our results demonstrate a possible pathway whereby OsLUGL-OsSEU-OsAP1-OsSEP3 complex as a transcriptional co-regulator by targeting the promoter of OsGH3-8, then affecting auxin level, OsARFs expression and thereby influencing floral development. These findings provide a valuable insight into the molecular functions of OsLUGL in rice floral development.

Homologs of LEAFY and UNUSUAL FLORAL ORGANS Promote the Transition From Inflorescence to Floral Meristem Identity in the Cymose Aquilegia coerulea

Homologs of the transcription factor LEAFY (LFY) and the F-box family member UNUSUAL FLORAL ORGANS (UFO) have been found to promote floral meristem identity across diverse dicot model systems. The lower eudicot model Aquilegia produces cymose inflorescences that are independently evolved from the well-studied cymose models Petunia and tomato. We have previously characterized the expression pattern of the Aquilegia homolog AqLFY but in the current study, we add expression data on the two UFO homologs, AqUFO1 and 2, and conduct virus-induced gene silencing of all the loci. Down-regulation of AqLFY or AqUFO1 and 2 does not eliminate floral meristem identity but, instead, causes the transition from inflorescence to floral identity to become gradual rather than discrete. Inflorescences in down-regulated plants generate several nodes of bract/sepal chimeras and, once floral development does commence, flowers initiate several whorls of sepals before finally producing the wildtype floral whorls. In addition, silencing of AqUFO1/2 appears to specifically impact petal identity and/or the initiation of petal and stamen whorls. In general, however, there is no evidence for an essential role of AqLFY or AqUFO1/2 in transcriptional activation of the B or C gene homologs. These findings highlight differences between deeply divergent dicot lineages in the functional conservation of the floral meristem identity program.

Structural Basis for Plant MADS Transcription Factor Oligomerization

MADS transcription factors (TFs) are DNA binding proteins found in almost all eukaryotes that play essential roles in diverse biological processes. While present in animals and fungi as a small TF family, the family has dramatically expanded in plants over the course of evolution, with the model flowering plant, Arabidopsis thaliana, possessing over 100 type I and type II MADS TFs. All MADS TFs contain a core and highly conserved DNA binding domain called the MADS or M domain. Plant MADS TFs have diversified this domain with plant-specific auxiliary domains. Plant type I MADS TFs have a highly diverse and largely unstructured Carboxy-terminal (C domain), whereas type II MADS have added oligomerization domains, called Intervening (I domain) and Keratin-like (K domain), in addition to the C domain. In this mini review, we describe the overall structure of the type II "MIKC" type MADS TFs in plants, with a focus on the K domain, a critical oligomerization module. We summarize the determining factors for oligomerization and provide mechanistic insights on how secondary structural elements are required for oligomerization capability and specificity. Using MADS TFs that are involved in flower organ specification as an example, we provide case studies and homology modeling of MADS TFs complex formation. Finally, we highlight outstanding questions in the field.

RNA interference suppression of AGAMOUS and SEEDSTICK alters floral organ identity and impairs floral organ determinacy, ovule differentiation, and seed-hair development in Populus

The role of the floral homeotic gene AGAMOUS (AG) and its close homologues in development of anemophilous, unisexual catkins has not previously been studied. We transformed two RNA interference (RNAi) constructs, PTG and its matrix-attachment-region flanked version MPG, into the early-flowering female poplar clone 6K10 (Populus alba) to suppress the expression of its two duplicate AG orthologues. By early 2018, six out of 22 flowering PTG events and 11 out of 12 flowering MPG events showed modified floral phenotypes in a field trial in Oregon, USA. Flowers in catkins from modified events had 'carpel-inside-carpel' phenotypes. Complete disruption of seed production was observed in seven events, and sterile anther-like organs in 10 events. Events with strong co-suppression of both the two AG and two SEEDSTICK (STK) paralogues lacked both seeds and associated seed hairs. Alterations in all of the modified floral phenotypes were stable over 4 yr of study. Trees from floral-modified events did not differ significantly (P < 0.05) from nonmodified transgenic or nontransgenic controls in biomass growth or leaf morphology. AG and STK genes show strong conservation of gene function during poplar catkin development and are promising targets for genetic containment of exotic or genetically engineered trees.

Comparative de novo flower transcriptome analysis of polygamodioecious tree Garcinia indica

To extend our understanding of molecular mechanism of sex determination in agro-economically important, polygamodioecious tree Garcinia indica (Kokum), high-throughput, next-generation flower transcriptome sequencing (NGS), and comparative analyses were performed to investigate differentially expressed gene in bisexual, female, and male flowers. A total of 49414 unigenes in BS, 45944 unigenes in FL, and 49028 unigenes in ML flowers were annotated. KO annotations revealed that 25 functional categories were large number of genes which were annotated to 'signal transduction'. We identified 33 genes for 'auxin response factor' and 50 for 'ethylene-responsive factor' whose expression changed significantly in all the three paired library combinations. Furthermore, key regulators of floral development such as FLC, SVP, AP1, AP2, AP3, AG, AGL2, AGL4, AGL9, and PI were identified. A total of 327 differentially expressed MADS-box genes were identified in G. indica transcriptome. Analysis of MADS-box genes identified five genes such as MADS AGL11, CRS2-associated factor chloroplastic, conserved hypothetical protein, uncharacterized protein LOC104422218, and MADS-box JOINTLESS-like isoform X3 significantly expressed in only FL flower. In addition, number of DEGs like dynamin 2A, auxin response factor, and spermidine synthase involved in sex expression and reproduction were discovered. The expression patterns of selected genes matched well with the expression levels of unigenes by transcriptome sequencing. Our large-scale comparative analyses may provide valuable hints for the next insights into the molecular mechanism of sex determination in G. indica.

Floral MADS-box protein interactions in the early diverging angiosperm Aristolochia fimbriata Cham. (Aristolochiaceae: Piperales)

Floral identity MADS-box A, B, C, D, E, and AGL6 class genes are predominantly single copy in Magnoliids, and predate the whole genome duplication (WGD) events in monocots and eudicots. By comparison with the model species Arabidopsis thaliana, the expression patterns of B-, C-, and D-class genes in stamen, carpel, and ovules are conserved in Aristolochia fimbriata, whereas A-, E-class, and AGL6 genes have different expression patterns. Nevertheless, the interactions of these proteins that act through multimeric complexes remain poorly known in early divergent angiosperms. This study evaluates protein interactions among all floral MADS-box A. fimbriata proteins using the Yeast Two Hybrid System (Y2H). We found no homodimers and less heterodimers formed by AfimFUL when compared to AfimAGL6, which allowed us to suggest AGL6 homodimers in combination with AfimSEP2 as the most likely tetramer in sepal identity. We found AfimAP3-AfimPI obligate heterodimers and AfimAG-AfimSEP2 protein interactions intact suggesting conserved stamen and carpel tetrameric complexes in A. fimbriata. We observed a broader interaction partner set for AfimSEP2 than for its paralog AfimSEP1. We show conserved and exclusive MADS-box protein interactions in A. fimbriata in comparison with other eudicot and monocot model species in order to establish plesiomorphic MADS-box protein floral networks in angiosperms.

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.