MIKCC-type MADS-box genes in Rosa chinensis: the remarkable expansion of ABCDE model genes and their roles in floral organogenesis

Flower development and a functional analysis of related genes in Impatiens uliginosa

Impatiens uliginosa is a plant of the Impatiens, with peculiar flowers. In this study, we combined morphogenesis and molecular biology to explore its development of flowers. An analysis of basic observational data and paraffin sectioning showed that it took approximately 13 d for the floral organs to differentiate. An analysis of the development of inflorescences and floral organs by scanning electron microscopy showed that the inflorescence of I. uliginosa is a spiral raceme. The floral organs differentiated in the following order: lateral sepals (Ls), posterior sepal (Ps), anterior sepals (As), anterior petal (Ap), lateral petals (Lp), stamens (St) and gynoecium (Gy). I. uliginosa was found to have four sepals, and the connate stamens are caused by the fusion and growth of filament appendages. The results of fluorescence quantification and virus-induced gene silencing showed that I. uliginosa had its own unique model for flower development, and there was functional diversity of IuAP1 and IuDEF. Among them, IuAP1 controls the formation of bract s (Br), regulates the morphogenesis of posterior sepal, controls the anthocyanin precipitation of the anterior petals and the formation of lateral petals. IuDEF regulates the morphogenesis of lateral sepals, the length of development of the spur, and controls the size of yellow flower color plaques of the lateral petals. In this study, the process of flower development and the function of flower development genes of I. uliginosa were preliminarily verified. This study provides basic guidance and new concepts that can be used to study the development of Impatiens flowers.

Characterization of the MIKCC-type MADS-box gene family in blueberry and its possible mechanism for regulating flowering in response to the chilling requirement

MAIN CONCLUSION

The expression peak of VcAP1.4, VcAP1.6, VcAP3.1, VcAP3.2, VcAG3, VcFLC2, and VcSVP9 coincided with the endo-dormancy release of flower buds. Additionally, GA4+7 not only increased the expression of these genes but also promoted flower bud endo-dormancy release. The MIKCC-type MADS-box gene family is involved in the regulation of flower development. A total of 109 members of the MIKCC-type MADS-box gene family were identified in blueberry. According to the phylogenetic tree, these 109 MIKCC-type MADS-box proteins were divided into 13 subfamilies, which were distributed across 40 Scaffolds. The results of the conserved motif analysis showed that among 20 motifs, motifs 1, 3, and 9 formed the MADS-box structural domain, while motifs 2, 4, and 6 formed the K-box structural domain. The presence of 66 pairs of fragment duplication events in blueberry suggested that gene duplication events contributed to gene expansion and functional differentiation. Additionally, the presence of cis-acting elements revealed that VcFLC2, VcAG3, and VcSVP9 might have significant roles in the endo-dormancy release of flower buds. Meanwhile, under chilling conditions, VcAP3.1 and VcAG7 might facilitate flower bud dormancy release. VcSEP11 might promote flowering following the release of endo-dormancy, while the elevated expression of VcAP1.7 (DAM) could impede the endo-dormancy release of flower buds. The effect of gibberellin (GA4+7) treatment on the expression pattern of MIKCC-type MADS-box genes revealed that VcAP1.4, VcAP1.6, VcAP3.1, VcAG3, and VcFLC2 might promote flower bud endo-dormancy release, while VcAP3.2, VcSEP11, and VcSVP9 might inhibit its endo-dormancy release. These results indicated that VcAP1.4, VcAP1.6, VcAP1.7 (DAM), VcAP3.1, VcAG3, VcAG7, VcFLC2, and VcSVP9 could be selected as key regulatory promoting genes for controlling the endo-dormancy of blueberry flower buds.

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

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

RcSPL1-RcTAF15b regulates the flowering time of rose (Rosa chinensis)

Rose (Rosa chinensis), which is an economically valuable floral species worldwide, has three types, namely once-flowering (OF), occasional or re-blooming (OR), and recurrent or continuous flowering (CF). However, the mechanism underlying the effect of the age pathway on the duration of the CF or OF juvenile phase is largely unknown. In this study, we observed that the RcSPL1 transcript levels were substantially upregulated during the floral development period in CF and OF plants. Additionally, accumulation of RcSPL1 protein was controlled by rch-miR156. The ectopic expression of RcSPL1 in Arabidopsis thaliana accelerated the vegetative phase transition and flowering. Furthermore, the transient overexpression of RcSPL1 in rose plants accelerated flowering, whereas silencing of RcSPL1 had the opposite phenotype. Accordingly, the transcription levels of floral meristem identity genes (APETALA1, FRUITFULL, and LEAFY) were significantly affected by the changes in RcSPL1 expression. RcTAF15b protein, which is an autonomous pathway protein, was revealed to interact with RcSPL1. The silencing and overexpression of RcTAF15b in rose plants led to delayed and accelerated flowering, respectively. Collectively, the study findings imply that RcSPL1-RcTAF15b modulates the flowering time of rose plants.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

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.

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

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

Evolutionary Analysis and Functional Identification of Clock-Associated PSEUDO-RESPONSE REGULATOR (PRRs) Genes in the Flowering Regulation of Roses

Pseudo-response regulators (PRRs) are the important genes for flowering in roses. In this work, clock PRRs were genome-wide identified using Arabidopsis protein sequences as queries, and their evolutionary analyses were deliberated intensively in Rosaceae in correspondence with angiosperms species. To draw a comparative network and flow of clock PRRs in roses, a co-expression network of flowering pathway genes was drawn using a string database, and their functional analysis was studied by silencing using VIGS and protein-to-protein interaction. We revealed that the clock PRRs were significantly expanded in Rosaceae and were divided into three major clades, i.e., PRR5/9 (clade 1), PRR3/7 (clade 2), and TOC1/PRR1 (clade 3), based on their phylogeny. Within the clades, five clock PRRs were identified in Rosa chinensis. Clock PRRs had conserved RR domain and shared similar features, suggesting the duplication occurred during evolution. Divergence analysis indicated the role of duplication events in the expansion of clock PRRs. The diverse cis elements and interaction of clock PRRs with miRNAs suggested their role in plant development. Co-expression network analysis showed that the clock PRRs from Rosa chinensis had a strong association with flowering controlling genes. Further silencing of RcPRR1b and RcPRR5 in Rosa chinensis using VIGS led to earlier flowering, confirming them as negative flowering regulators. The protein-to-protein interactions between RcPRR1a/RcPRR5 and RcCO suggested that RcPRR1a/RcPRR5 may suppress flowering by interfering with the binding of RcCO to the promoter of RcFT. Collectively, these results provided an understanding of the evolutionary profiles as well as the functional role of clock PRRs in controlling flowering in roses.

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

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

MADS-box transcription factors determine the duration of temporary winter dormancy in closely related evergreen and deciduous Iris spp

Winter dormancy (WD) is a crucial strategy for plants coping with potentially deadly environments. In recent decades, this process has been extensively studied in economically important perennial eudicots due to changing climate. However, in evergreen monocots with no chilling requirements, dormancy processes are so far a mystery. In this study, we compared the WD process in closely related evergreen (Iris japonica) and deciduous (I. tectorum) iris species across crucial developmental time points. Both iris species exhibit a 'temporary' WD process with distinct durations, and could easily resume growth under warm conditions. To decipher transcriptional changes, full-length sequencing for evergreen iris and short read RNA sequencing for deciduous iris were applied to generate respective reference transcriptomes. Combining results from a multipronged approach, SHORT VEGETATIVE PHASE and FRUITFULL (FUL) from MADS-box was associated with a dormancy- and a growth-related module, respectively. They were co-expressed with genes involved in phytohormone signaling, carbohydrate metabolism, and environmental adaptation. Also, gene expression patterns and physiological changes in the above pathways highlighted potential abscisic acid and jasmonic acid antagonism in coordinating growth and stress responses, whereas differences in carbohydrate metabolism and reactive oxygen species scavenging might lead to species-specific WD durations. Moreover, a detailed analysis of MIKCCMADS-box in irises revealed common features described in eudicots as well as possible new roles for monocots during temporary WD, such as FLOWERING LOCUS C and FUL. In essence, our results not only provide a portrait of temporary WD in perennial monocots but also offer new insights into the regulatory mechanism underlying WD in plants.

Genome-wide analysis reveals widespread roles for RcREM genes in floral organ development in Rosa chinensis

Members of the REM (Reproductive Meristem) gene family are expressed primarily in reproductive meristems and floral organs. However, their evolution and their functional profiles in flower development remain poorly understood. Here, we performed genome-wide identification and evolutionary analysis of the REM gene family in Rosaceae. This family has been greatly expanded in rose (Rosa chinensis) compared to other species, primarily through tandem duplication. Expression analysis revealed that most RcREM genes are specifically expressed in reproductive organs and that their specific expression patterns are dramatically altered in rose plants with mutations affecting floral organs. Protein-protein interaction analysis indicated that RcREM14 interact with RcAP1 (one of the homology of A class genes in ABCDE model), highlighting the roles of RcREM genes in floral organ identity. Finally, co-expression network analysis indicated that RcREM genes are co-expressed with a high proportion of key genes that regulate flowering time, floral organ development, and cell proliferation and expansion in R. chinensis.

Genome-Wide Identification MIKC-Type MADS-Box Gene Family and Their Roles during Development of Floral Buds in Wheel Wingnut (Cyclocarya paliurus)

MADS-box transcription factors (TFs) have fundamental roles in regulating floral organ formation and flowering time in flowering plants. In order to understand the function of MIKC-type MADS-box family genes in Cyclocarya paliurus (Batal.) Iljinskaja, we first implemented a genome-wide analysis of MIKC-type MADS-box genes in C. paliurus. Here, the phylogenetic relationships, chromosome location, conserved motif, gene structure, promoter region, and gene expression profile were analyzed. The results showed that 45 MIKC-type MADS-box were divided into 14 subfamilies: BS (3), AGL12 (1), AP3-PI (3), MIKC* (3), AGL15 (3), SVP (5), AGL17 (2), AG (3), TM8 (1), AGL6 (2), SEP (5), AP1-FUL (6), SOC1 (7), and FLC (1). The 43 MIKC-type MADS-box genes were distributed unevenly in 14 chromosomes, but two members were mapped on unanchored scaffolds. Gene structures were varied in the same gene family or subfamily, but conserved motifs shared similar distributions and sequences. The element analysis in promoters’ regions revealed that MIKC-type MADS-box family genes were associated with light, phytohormone, and temperature responsiveness, which may play important roles in floral development and differentiation. The expression profile showed that most MIKC-type MADS-box genes were differentially expressed in six tissues (specifically expressed in floral buds), and the expression patterns were also visibly varied in the same subfamily. CpaF1st24796 and CpaF1st23405, belonging to AP3-PI and SEP subfamilies, exhibited the high expression levels in PA-M and PG-F, respectively, indicating their functions in presenting heterodichogamy. We further verified the MIKC-type MADS-box gene expression levels on the basis of transcriptome and qRT-PCR analysis. This study would provide a theoretical basis for classification, cloning, and regulation of flowering mechanism of MIKC-type MADS-box genes in C. paliurus.

Identification and Characterization of MIKCc-Type MADS-Box Genes in the Flower Organs of Adonis amurensis

Adonis amurensis is a perennial herbaceous flower that blooms in early spring in northeast China, where the night temperature can drop to −15 °C. To understand flowering time regulation and floral organogenesis of A. amurensis, the MIKC^c-type MADS (Mcm1/Agamous/ Deficiens/Srf)-box genes were identified and characterized from the transcriptomes of the flower organs. In this study, 43 non-redundant MADS-box genes (38 MIKC^c, 3 MIKC*, and 2 Mα) were identified. Phylogenetic and conserved motif analysis divided the 38 MIKC^c-type genes into three major classes: ABCDE model (including AP1/FUL, AP3/PI, AG, STK, and SEPs/AGL6), suppressor of overexpression of constans1 (SOC1), and short vegetative phase (SVP). qPCR analysis showed that the ABCDE model genes were highly expressed mainly in flowers and differentially expressed in the different tissues of flower organs, suggesting that they may be involved in the flower organ identity of A. amurensis. Subcellular localization revealed that 17 full-length MADSs were mainly localized in the nucleus: in Arabidopsis, the heterologous expression of three full-length SOC1-type genes caused early flowering and altered the expression of endogenous flowering time genes. Our analyses provide an overall insight into MIKC^c genes in A. amurensis and their potential roles in floral organogenesis and flowering time regulation.

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

Background

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

Results

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

Conclusions

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

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

Supplementary Information

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

Phytochrome-interacting factors interact with transcription factor CONSTANS to suppress flowering in rose

As day-neutral (DN) woody perennial plants, the flowering time of roses (Rosa spp.) is assumed to be independent of the photoperiodic conditions; however, light responses of rose plants are not well understood. Chinese rose (Rosa chinensis) plants were grown under two light intensities (low light [LL], 92 μmol·m-2·s-1; or high light [HL], 278 μmol·m-2·s-1), and either with or without an end-of-day far-red (EOD-FR) treatment. Flowering was significantly delayed in the LL condition compared with the HL, but was not affected by EOD-FR treatment. The time until flowering positively corresponded with the mRNA and protein levels of phytochrome-interacting factors (PIFs; RcPIFs). The heterologous expression of RcPIF1, RcPIF3, or RcPIF4 in the Arabidopsis (Arabidopsis thaliana) pifq quadruple mutant partially rescued the mutant's shorter hypocotyl length. Simultaneous silencing of three RcPIFs in R. chinensis accelerated flowering under both LL and HL, with a more robust effect in LL, establishing RcPIFs as flowering suppressors in response to light intensity. The RcPIFs interacted with the transcription factor CONSTANS (RcCO) to form a RcPIFs-RcCO complex, which interfered with the binding of RcCO to the promoter of FLOWERING LOCUS T (RcFT), thereby inhibiting its expression. Furthermore, this inhibition was enhanced when RcPIFs were stabilized by LL, leading to delayed flowering under LL compared with HL. Our results not only revealed another layer of PIF functioning in the flowering of woody perennial plants, but also established a mechanism of light response in DN plants.

The SOC1-like gene BoMADS50 is associated with the flowering of Bambusa oldhamii

Bamboo is known for its edible shoots and beautiful texture and has considerable economic and ornamental value. Unique among traditional flowering plants, many bamboo plants undergo extensive synchronized flowering followed by large-scale death, seriously affecting the productivity and application of bamboo forests. To date, the molecular mechanism of bamboo flowering characteristics has remained unknown. In this study, a SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1)-like gene, BoMADS50, was identified from Bambusa oldhamii. BoMADS50 was highly expressed in mature leaves and the floral primordium formation period during B. oldhamii flowering and overexpression of BoMADS50 caused early flowering in transgenic rice. Moreover, BoMADS50 could interact with APETALA1/FRUITFULL (AP1/FUL)-like proteins (BoMADS14-1/2, BoMADS15-1/2) in vivo, and the expression of BoMADS50 was significantly promoted by BoMADS14-1, further indicating a synergistic effect between BoMADS50 and BoAP1/FUL-like proteins in regulating B. oldhamii flowering. We also identified four additional transcripts of BoMADS50 (BoMADS50-1/2/3/4) with different nucleotide variations. Although the protein-CDS were polymorphic, they had flowering activation functions similar to those of BoMADS50. Yeast one-hybrid and transient expression assays subsequently showed that both BoMADS50 and BoMADS50-1 bind to the promoter fragment of itself and the SHORT VEGETATIVE PHASE (SVP)-like gene BoSVP, but only BoMADS50-1 can positively induce their transcription. Therefore, nucleotide variations likely endow BoMADS50-1 with strong regulatory activity. Thus, BoMADS50 and BoMADS50-1/2/3/4 are probably important positive flowering regulators in B. oldhamii. Moreover, the functional conservatism and specificity of BoMADS50 and BoMADS50-1 might be related to the synchronized and sporadic flowering characteristics of B. oldhamii.

Genes and genome editing tools for breeding desirable phenotypes in ornamentals

We review the main genes underlying commercial traits in cut flower species and critically discuss the possibility to apply genome editing approaches to produce novel variation and phenotypes. Promoting flowering and flower longevity as well as creating novelty in flower structure, colour range and fragrances are major objectives of ornamental plant breeding. The novel genome editing techniques add new possibilities to study gene function and breed new varieties. The implementation of such techniques, however, relies on detailed information about structure and function of genomes and genes. Moreover, improved protocols for efficient delivery of editing reagents are required. Recent results of the application of genome editing techniques to elite ornamental crops are discussed in this review. Enabling technologies and genomic resources are reviewed in relation to the implementation of such approaches. Availability of the main gene sequences, underlying commercial traits and in vitro transformation protocols are provided for the world's best-selling cut flowers, namely rose, lily, chrysanthemum, lisianthus, tulip, gerbera, freesia, alstroemeria, carnation and hydrangea. Results obtained so far are described and their implications for the improvement of flowering, flower architecture, colour, scent and shelf-life are discussed.

The pineapple MADS-box gene family and the evolution of early monocot flower

Unlike the flower of the model monocot rice, which has diverged greatly from the ancestral monocot flower, the pineapple (Ananas comosus) flower is more typical of monocot flowers. Here, we identified 43 pineapple genes containing MADS-box domains, including 11 type I and 32 type II genes. RNA-seq expression data generated from five pineapple floral organs (sepals, petals, stamens, pistils, and ovules) and quantitative real-time PCR revealed tissue-specific expression patterns for some genes. We found that AcAGL6 and AcFUL1 were mainly expressed in sepals and petals, suggesting their involvement in the regulation of these floral organs. A pineapple ‘ABCDE’ model was proposed based on the phylogenetic analysis and expression patterns of MADS-box genes. Unlike rice and orchid with frequent species-specific gene duplication and subsequent expression divergence, the composition and expression of the ABCDE genes were conserved in pineapple. We also found that AcSEP1/3, AcAG, AcAGL11a/b/c, and AcFUL1 were highly expressed at different stages of fruit development and have similar expression profiles, implicating these genes’ role in fruit development and ripening processes. We propose that the pineapple flower can be used as a model for studying the ancestral form of monocot flowers to investigate their development and evolutionary history.

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.

Genome-wide identification of MADS-box gene family in sacred lotus (Nelumbo nucifera) identifies a SEPALLATA homolog gene involved in floral development

BACKGROUND

Sacred lotus (Nelumbo nucifera) is a vital perennial aquatic ornamental plant. Its flower shape determines the horticultural and ornamental values. However, the mechanisms underlying lotus flower development are still elusive. MADS-box transcription factors are crucial in various features of plant development, especially in floral organogenesis and specification. It is still unknown how the MADS-box transcription factors regulate the floral organogenesis in lotus.

RESULTS

To obtain a comprehensive insight into the functions of MADS-box genes in sacred lotus flower development, we systematically characterized members of this gene family based on the available genome information. A total of 44 MADS-box genes were identified, of which 16 type I and 28 type II genes were categorized based on the phylogenetic analysis. Furthermore, the structure of MADS-box genes and their expressional patterns were also systematically analyzed. Additionally, subcellular localization analysis showed that they are mainly localized in the nucleus, of which a SEPALLATA3 (SEP3) homolog NnMADS14 was proven to be involved in the floral organogenesis.

CONCLUSION

These results provide some fundamental information about the MADS-box gene family and their functions, which might be helpful in not only understanding the mechanisms of floral organogenesis but also breeding of high ornamental value cultivars in lotus.

Genome-wide identification of WD40 genes reveals a functional diversification of COP1-like genes in Rosaceae

Genome-wide identification of WD40-like genes reveals a duplication of COP1-like genes, one of the key players involved in regulation of flowering time and photomorphogenesis, with strong functional diversification in Rosaceae. WD40 proteins play crucial roles in a broad spectrum of developmental and physiological processes. Here, we conducted a systematic characterization of this family of genes in Rosa chinensis 'Old Blush' (OB), a founder genotype for modern rose domestication. We identified 187 rose WD40 genes and classified them into 5 clusters and 15 subfamilies with 11 of RcWD40s presumably generated via tandem duplication. We found RcWD40 genes were expressed differentially following stages of vegetative and reproductive development. We detected a duplication of CONSTITUTIVE PHOTOMORPHOGENIC1-like genes in rose (RcCOP1 and RcCOP1L) and other Rosaceae plants. Featuring a distinct expression pattern and a different profile of cis-regulatory-elements in the transcriptional regulatory regions, RcCOP1 seemed being evolutionarily conserved while RcCOP1L did not dimerize with RcHY5 and RcSPA4. Our data thus reveals a functional diversification of COP1-like genes in Rosacaeae plants, and provides a valuable resource to explore the potential function and evolution of WD40-like genes in Rosaceae plants.

RcMYB84 and RcMYB123 mediate jasmonate-induced defense responses against Botrytis cinerea in rose (Rosa chinensis)

Jasmonates (JAs) are important for pathogen resistance in many plants, but the role of these phytohormones in fungal pathogen resistance in rose is unclear. Here, we determined that exogenous application of methyl jasmonate increased resistance to the important fungal pathogen Botrytis cinerea in Rosa chinensis 'Old blush', whereas silencing the JA biosynthetic pathway gene Allene Oxide Synthase (AOS) and JA co-receptor gene CORONATINE INSENSITIVE 1 (COI1) suppressed this response. Transcriptome profiling identified various MYB transcription factor genes that responded to both JA and B. cinerea treatment. Silencing Ri-RcMYB84/Ri-RcMYB123 increased the susceptibility of rose plants to B. cinerea and inhibited the protective effects of JA treatment, confirming the crucial roles of these genes in JA-induced responses to B. cinerea. JAZ1, a key repressor of JA signaling, directly interacts with RcMYB84 and RcMYB123 to deplete their free pools. The JAZ1-RcMYB84 complex binds to the RcMYB123 promoter via the CAACTG motifs to block its transcription. Upon JA treatment, the expression of RcMYB123 is de-repressed, and free forms of RcMYB84 and RcMYB123 are released due to JAZ1 degradation, thereby activating the defense responses of plants to B. cinerea. These findings shed light on the molecular mechanisms underlying JA-induced pathogen resistance in roses.

Alternate expression of CONSTANS-LIKE 4 in short days and CONSTANS in long days facilitates day-neutral response in Rosa chinensis

Photoperiodic flowering responses are classified into three major types: long day (LD), short day (SD), and day neutral (DN). The inverse responses to daylength of LD and SD plants have been partly characterized in Arabidopsis and rice; however, the molecular mechanism underlying the DN response is largely unknown. Modern roses are economically important ornamental plants with continuous flowering (CF) features, and are generally regarded as DN plants. Here, RcCO and RcCOL4 were identified as floral activators up-regulated under LD and SD conditions, respectively, in the CF cultivar Rosa chinensis 'Old-Blush'. Diminishing the expression of RcCO or/and RcCOL4 by virus-induced gene silencing (VIGS) delayed flowering time under both SDs and LDs. Interestingly, in contrast to RcCO-silenced plants, the flowering time of RcCOL4-silenced plants was more delayed under SD than under LD conditions, indicating perturbed plant responses to day neutrality. Further analyses revealed that physical interaction between RcCOL4 and RcCO facilitated binding of RcCO to the CORE motif in the promoter of RcFT and induction of RcFT. Taken together, the complementary expression of RcCO in LDs and of RcCOL4 in SDs guaranteed flowering under favorable growth conditions regardless of the photoperiod. This finding established the molecular foundation of CF in roses and further shed light on the underlying mechanisms of DN responses.

Genome-wide identification and functional analysis of JmjC domain-containing genes in flower development of Rosa chinensis

We genome-wide identified 28 JmjC domain-containing genes, further spatio-temporal expression profiling and genetic analysis defined them as epigenetic regulators in flowering initiation of Rosa chinensis. The JmjC domain-containing histone demethylases play critical roles in maintaining homeostasis of histone methylations, thus are vital for plant growth and development. Genome-wide identification of the JmjC domain-containing genes have been reported in several species, however, no systematic study has been performed in rose plants. In this paper, we identified 28 JmjC domain-containing genes from the newly published genome database of Rosa chinensis. The JmjC domain-containing proteins in R. chinensis were divided into seven groups, KDM3 was the largest group with 13 members, and JmjC domain-only A and KDM5B were the smallest clades both with only one member. Although all the JmjC domain proteins having a conserved JmjC domain, the gene and protein structure experienced differentiation and specification during the evolution, especially in KDM3 clade, one gene (RcJMJ40) was found carrying site deletions for cofactors Fe (II) and α-KG binding which were crucial for demethylase activities, three genes (RcJMJ41, RcJMJ43 and RcJMJ44) had no intron while two of them had tandem JmjC domains. Spatial expression pattern analysis of these JmjC domain-containing genes in different tissues showed most of them were highly expressed in reproductive tissues such as floral meristem and closed flowers than vegetative tissues, demonstrating their important functions in developmental switch from vegetative to reproductive growth of roses. Temporal expression profiling indicated majority of JmjC domain-containing genes from R. chinensis fluctuated along with floral bud differentiation and development, further proving their essential roles in flower organogenesis. VIGS induced silencing of RcJMJ12 led to delayed flowering time, and decreased the expression levels of flowering integrator such as RcFT, RcSOC1, RcFUL, RcLFY and RcAP1, therefore providing the genetic evidence of RcJMJ12 in flowering initiation. Collectively, spatio-temporal expression profiling and genetic analysis defined the JmjC domain-containing genes as important epigenetic regulators in flower development of R. chinensis.

Differentially Expressed Genes between Carrot Petaloid Cytoplasmic Male Sterile and Maintainer during Floral Development

Petaloid cytoplasmic male sterility (CMS) is a maternally inherited loss of male fertility due to the complete conversion of stamens into petal-like organs, and CMS lines have been widely utilized in carrot breeding. Petaloid CMS is an ideal model not only for studying the mitochondrial-nuclear interaction but also for discovering genes that are essential for floral organ development. To investigate the comprehensive mechanism of CMS and homeotic organ alternation during carrot flower development, we conducted transcriptome analysis between the petaloid CMS line (P2S) and its maintainer line (P2M) at four flower developmental stages (T1-T4). A total of 2838 genes were found to be differentially expressed, among which 1495 genes were significantly downregulated and 1343 genes were significantly upregulated in the CMS line. Functional analysis showed that most of the differentially expressed genes (DEGs) were involved in protein processing in the endoplasmic reticulum, plant hormone signal transduction, and biosynthesis. A total of 16 MADS-box genes were grouped into class A, B, C, and E, but not class D, genes. Several key genes associated with oxidative phosphorylation showed continuously low expression from stage T2 in P2S, and the expression of DcPI and DcAG-like genes also greatly decreased at stage T2 in P2S. This indicated that energy deficiency might inhibit the expression of B- and C-class MADS-box genes resulting in the conversion of stamens into petals. Stamen petaloidy may act as an intrinsic stress, upregulating the expression of heat shock protein (HSP) genes and MADS-box genes at stages T3 and T4 in P2S, which results in some fertile revertants. This study will provide a better understanding of carrot petaloid CMS and floral development as a basis for further research.

Morphological Characterization of Flower Buds Development and Related Gene Expression Profiling at Bud Break Stage in Heterodichogamous Cyclocarya paliurus (Batal.) lljinskaja

Cyclocarya paliurus (Batal.) Iljinskaja, a unique species growing in southern China, is a multi-function tree species with medicinal, healthcare, material, and ornamental values. So far, sexual reproduction is the main method for extensive cultivation of C. paliurus plantations, but this is limited by low seed plumpness resulted from the character of heterodichogamy. Phenological observations have revealed the asynchronism of flower development in this species. However, its molecular mechanism remains largely unknown. To reveal molecular mechanism of heterodichogamy in C. paliurus, transcriptome of female (F) and male (M) buds from two mating types (protandry, PA; protogyny, PG) at bud break stage were sequenced using Illumina Hiseq 4000 platform. The expression patterns of both 32 genes related to flowering and 58 differentially expressed transcription factors (DETFs) selected from 6 families were divided four groups (PG-F, PG-M, PA-F, and PA-M) into two categories: first flowers (PG-F and PA-M) and later flowers (PA-F and PG-M). The results indicated that genes related to plant hormones (IAA, ABA, and GA) synthesis and response, glucose metabolism, and transcription factors (especially in MIKC family) played significant roles in regulating asynchronism of male and female flowers in the same mating type. The expression of DETFs showed two patterns. One contained DETFs up-regulated in first flowers in comparison to later flowers, and the other was the reverse. Nine genes related to flowering were selected for qRT-PCR to confirm the accuracy of RNA-seq, and generally, the RPKM values of these genes were consistent with the result of qRT-PCR. The results of this work could improve our understanding in asynchronism of floral development within one mating type in C. paliurus at transcriptional level, as well as lay a foundation for further study in heterodichogamous plants.

Targeted deletion of floral development genes in Arabidopsis with CRISPR/Cas9 using the RNA endoribonuclease Csy4 processing system

The formation of flowers in higher plants is controlled by complex gene regulatory networks. The study of floral development in Arabidopsis is promoted and maintained by transposon-tagged mutant lines. In this study, we report a CRISPR/Cas9 genome-editing system based on RNA endoribonuclease Csy4 processing to induce high-efficiency and inheritable targeted deletion of transcription factors involved in floral development in Arabidopsis. Using AP1, SVP, and TFL1 as the target genes, multisite and multiple-gene mutations were achieved with a tandemly arrayed Csy4-sgRNA architecture to express multiplexed sgRNAs from a single transcript driven by the Pol II promoter in transgenic lines. Targeted deletions of chromosomal fragments between the first exon and second exon in either one or three genes were generated by using a single binary vector. Interestingly, the efficiency of site-targeted deletion was comparable to that of indel mutation with the multiplexed sgRNAs. DNA sequencing analysis of RT-PCR products showed that targeted deletions of AP1 and TFL1 could lead to frameshift mutations and introduce premature stop codons to disrupt the open-reading frames of the target genes. In addition, no RT-PCR amplified product was acquired after SVP-targeted deletion. Furthermore, the targeted deletions resulted in abnormal floral development in the mutant lines compared to that of wild-type plants. AP1 and SVP mutations increased plant branching significantly, while TFL1 mutant plants displayed a change from indeterminate to determinate inflorescences. Thus, our results demonstrate that CRISPR/Cas9 with the RNA endoribonuclease Csy4 processing system is an efficient tool to study floral development and improve floral traits rapidly and simply.

Identification of Transcription Factors Involved in the Regulation of Flowering in Adonis Amurensis Through Combined RNA-seq Transcriptomics and iTRAQ Proteomics

Temperature is one of the most important environmental factors affecting flowering in plants. Adonis amurensis, a perennial herbaceous flower that blooms in early spring in northeast China where the temperature can drop to −15 °C, is an ideal model for studying the molecular mechanisms of flowering at extremely low temperatures. This study first investigated global gene expression profiles at different developmental stages of flowering in A. amurensis by RNA-seq transcriptome and iTRAQ proteomics. Finally, 123 transcription factors (TFs) were detected in both the transcriptome and the proteome. Of these, 66 TFs belonging to 14 families may play a key role in multiple signaling pathways of flowering in A. amurensis. The TFs FAR1, PHD, and B3 may be involved in responses to light and temperature, while SCL, SWI/SNF, ARF, and ERF may be involved in the regulation of hormone balance. SPL may regulate the age pathway. Some members of the TCP, ZFP, MYB, WRKY, and bHLH families may be involved in the transcriptional regulation of flowering genes. The MADS-box TFs are the key regulators of flowering in A. amurensis. Our results provide a direction for understanding the molecular mechanisms of flowering in A. amurensis at low temperatures.

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

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