MADS-box genes and floral development: the dark side

Transcriptome wide analysis of MADS box genes in Crocus sativus and interplay of CstMADS19-CstMADS26 in orchestrating apocarotenoid biosynthesis

Flowers of Crocus sativus L. are immensely important not only for arrangement of floral whorls but more because each floral organ is dominated by a different class of specialized compounds. Dried stigmas of C. sativus flowers form commercial saffron, and are known to accumulate unique apocarotenoids like crocin, picrocrocin and safranal. Inspite of being a high value crop, the molecular mechanism regulating flower development in Crocus remains largely unknown. Moreover, it would be very interesting to explore any co-regulatory mechanism which controls floral architecture and secondary metabolic pathways which exist in specific floral organs. Here we report transcriptome wide identification of MADS box genes in Crocus. A total of 39 full length MADS box genes were identified among which three belonged to type I and 36 to type II class. Phylogeny classified them into 11 sub-clusters. Expression pattern revealed some stigma up-regulated genes among which CstMADS19 encoding an AGAMOUS gene showed high expression. Transient over-expression of CstMADS19 in stigmas of Crocus resulted in increased crocin by enhancing expression of pathway genes. Yeast one hybrid assay demonstrated that CstMADS19 binds to promoters of phytoene synthase and carotenoid cleavage dioxygenase 2 genes. Yeast two hybrid and BiFC assays confirmed interaction of CstMADS19 with CstMADS26 which codes for a SEPALATA gene. Co-overexpression of CstMADS19 and CstMADS26 in Crocus stigmas enhanced crocin content more than was observed when genes were expressed individually. Collectively, these findings indicate that CstMADS19 functions as a positive regulator of stigma based apocarotenoid biosynthesis in Crocus.

Progressive meristem and single-cell transcriptomes reveal the regulatory mechanisms underlying maize inflorescence development and sex differentiation

Maize develops separate ear and tassel inflorescences with initially similar morphology but ultimately different architecture and sexuality. The detailed regulatory mechanisms underlying these changes still remain largely unclear. In this study, through analyzing the time-course meristem transcriptomes and floret single-cell transcriptomes of ear and tassel, we revealed the regulatory dynamics and pathways underlying inflorescence development and sex differentiation. We identified 16 diverse gene clusters with differential spatiotemporal expression patterns and revealed biased regulation of redox, programmed cell death, and hormone signals during meristem differentiation between ear and tassel. Notably, based on their dynamic expression patterns, we revealed the roles of two RNA-binding proteins in regulating inflorescence meristem activity and axillary meristem formation. Moreover, using the transcriptional profiles of 53 910 single cells, we uncovered the cellular heterogeneity between ear and tassel florets. We found that multiple signals associated with either enhanced cell death or reduced growth are responsible for tassel pistil suppression, while part of the gibberellic acid signal may act non-cell-autonomously to regulate ear stamen arrest during sex differentiation. We further showed that the pistil-protection gene SILKLESS 1 (SK1) functions antagonistically to the known pistil-suppression genes through regulating common molecular pathways, and constructed a regulatory network for pistil-fate determination. Collectively, our study provides a deep understanding of the regulatory mechanisms underlying inflorescence development and sex differentiation in maize, laying the foundation for identifying new regulators and pathways for maize hybrid breeding and improvement.

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

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

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

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

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

Genome-wide analysis of the MADS-box gene family in Lonicera japonica and a proposed floral organ identity model

BACKGROUND

Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.

RESULTS

In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.

CONCLUSION

This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

Chromosome-scale genome assembly of Prunus pusilliflora provides novel insights into genome evolution, disease resistance, and dormancy release in Cerasus L

Prunus pusilliflora is a wild cherry germplasm resource distributed mainly in Southwest China. Despite its ornamental and economic value, a high-quality assembled P. pusilliflora genome is unavailable, hindering our understanding of its genetic background, population diversity, and evolutionary processes. Here, we de novo assembled a chromosome-scale P. pusilliflora genome using Oxford Nanopore, Illumina, and chromosome conformation capture sequencing. The assembled genome size was 309.62 Mb, with 76 scaffolds anchored to eight pseudochromosomes. We predicted 33 035 protein-coding genes, functionally annotated 98.27% of them, and identified repetitive sequences covering 49.08% of the genome. We found that P. pusilliflora is closely related to Prunus serrulata and Prunus yedoensis, having diverged from them ~41.8 million years ago. A comparative genomic analysis revealed that P. pusilliflora has 643 expanded and 1128 contracted gene families. Furthermore, we found that P. pusilliflora is more resistant to Colletotrichum viniferum, Phytophthora capsici, and Pseudomonas syringae pv. tomato (Pst) DC3000 infections than cultivated Prunus avium. P. pusilliflora also has considerably more nucleotide-binding site-type resistance gene analogs than P. avium, which explains its stronger disease resistance. The cytochrome P450 and WRKY families of 263 and 61 proteins were divided into 42 and 8 subfamilies respectively in P. pusilliflora. Furthermore, 81 MADS-box genes were identified in P. pusilliflora, accompanying expansions of the SVP and AGL15 subfamilies and loss of the TM3 subfamily. Our assembly of a high-quality P. pusilliflora genome will be valuable for further research on cherries and molecular breeding.

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.

Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii

Tree peony (Paeonia ostii) is an economically important ornamental plant native to China. It is also notable for its seed oil, which is abundant in unsaturated fatty acids such as α-linolenic acid (ALA). Here, we report chromosome-level genome assembly (12.28 Gb) of P. ostii. In contrast to monocots with giant genomes, tree peony does not appear to have undergone lineage-specific whole-genome duplication. Instead, explosive LTR expansion in the intergenic regions within a short period (~ two million years) may have contributed to the formation of its giga-genome. In addition, expansion of five types of histone encoding genes may have helped maintain the giga-chromosomes. Further, we conduct genome-wide association studies (GWAS) on 448 accessions and show expansion and high expression of several genes in the key nodes of fatty acid biosynthetic pathway, including SAD, FAD2 and FAD3, may function in high level of ALAs synthesis in tree peony seeds. Moreover, by comparing with cultivated tree peony (P. suffruticosa), we show that ectopic expression of class A gene AP1 and reduced expression of class C gene AG may contribute to the formation of petaloid stamens. Genomic resources reported in this study will be valuable for studying chromosome/genome evolution and tree peony breeding.

NBLAST: a graphical user interface-based two-way BLAST software with a dot plot viewer

BLAST, a basic bioinformatics tool for searching local sequence similarity, has been one of the most widely used bioinformatics programs since its introduction in 1990. Users generally use the web-based NCBI-BLAST program for BLAST analysis. However, users with large sequence data are often faced with a problem of upload size limitation while using the web-based BLAST program. This proves inconvenient as scientists often want to run BLAST on their own data, such as transcriptome or whole genome sequences. To overcome this issue, we developed NBLAST, a graphical user interface-based BLAST program that employs a two-way system, allowing the use of input sequences either as "query" or "target" in the BLAST analysis. NBLAST is also equipped with a dot plot viewer, thus allowing researchers to create custom database for BLAST and run a dot plot similarity analysis within a single program. It is available to access to the NBLAST with http://nbitglobal.com/nblast.

Comparative transcriptomic analysis provides insight into carpel petaloidy in lotus (Nelumbo nucifera)

Lotus (Nelumbo nucifera) is a highly recognized flower with high ornamental value. Flower color and flower morphology are two main factors for flower lotus breeding. Petaloidy is a universal phenomenon in lotus flowers. However, the genetic regulation of floral organ petaloidy in lotus remains elusive. In this study, the transcriptomic analysis was performed among three organs, including petal, carpel petaloidy, and carpel in lotus. A total of 1,568 DEGs related to carpel petaloidy were identified. Our study identified one floral homeotic gene encoded by the MADS-box transcription factor, AGAMOUS (AG) as the candidate gene for petaloid in lotus. Meanwhile, a predicted labile boundary in floral organs of N. nucifera was hypothesized. In summary, our results explored the candidate genes related to carpel petaloidy, setting a theoretical basis for the molecular regulation of petaloid phenotype.

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

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

De novo transcriptome characterization of Iris atropurpurea (the Royal Iris) and phylogenetic analysis of MADS-box and R2R3-MYB gene families

The Royal Irises (section Oncocyclus) are a Middle-Eastern group of irises, characterized by extremely large flowers with a huge range of flower colors and a unique pollination system. The Royal Irises are considered to be in the course of speciation and serve as a model for evolutionary processes of speciation and pollination ecology. However, no transcriptomic and genomic data are available for these plants. Transcriptome sequencing is a valuable resource for determining the genetic basis of ecological-meaningful traits, especially in non-model organisms. Here we describe the de novo transcriptome assembly of Iris atropurpurea, an endangered species endemic to Israel's coastal plain. We sequenced and analyzed the transcriptomes of roots, leaves, and three stages of developing flower buds. To identify genes involved in developmental processes we generated phylogenetic gene trees for two major gene families, the MADS-box and MYB transcription factors, which play an important role in plant development. In addition, we identified 1503 short sequence repeats that can be developed for molecular markers for population genetics in irises. This first reported transcriptome for the Royal Irises, and the data generated, provide a valuable resource for this non-model plant that will facilitate gene discovery, functional genomic studies, and development of molecular markers in irises, to complete the intensive eco-evolutionary studies of this group.

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

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

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.

Genetic and comparative mapping of Lupinus luteus L. highlight syntenic regions with major orthologous genes controlling anthracnose resistance and flowering time

Anthracnose susceptibility and ill-adapted flowering time severely affect Lupinus luteus yield, which has high seed protein content, is excellent for sustainable agriculture, but requires genetic improvement to fulfil its potential. This study aimed to (1) develop a genetic map; (2) define collinearity and regions of synteny with Lupinus angustifolius; and (3) map QTLs/candidate genes for anthracnose resistant and flowering time. A few linkage groups/genomic regions tended to be associated with segregation distortion, but did not affect the map. The developed map showed collinearity, and syntenic regions with L. angustifolius. Major QTLs were mapped in syntenic regions. Alleles from the wild parent and cultivar, explained 75% of the phenotypic variance for anthracnose resistance and 83% for early flowering, respectively. Marker sequences flanking the QTLs showed high homology with the Lanr1 gene and Flowering-locus-T of L. angustifolius. This suggests orthologous genes for both traits in the L. luteus genome. The findings are remarkable, revealing the potential to combine early flowering/anthracnose resistant in fulfilling yield capacity in L. luteus, and can be a major strategy in the genetic improvement and usage of this species for sustainable protein production. Allele sequences and PCR-marker tagging of these genes are being applied in marker assisted selection.

Ectopic expression of an AGAMOUS homologue gene in Jatropha curcas causes early flowering and heterostylous phenotypes

Jatropha curcasseeds are abundant in biodiesel, and low seed yields are linked to poor quality female flowers, which creates a bottleneck for Jatropha seed utilization. Therefore, identifying the genes associated with flowering is crucial for the genetic enrichment of seed yields. Here, we identified an AGAMOUS homologue gene (JcAG) from J. curcas. We found that reproductive organs had higher JcAG expression than vegetative organs, particularly the carpel. Rosette leaves were small and misshapen in 35S:JcAG transgenic lines in comparison with those in wild-type plants. JcAG overexpression caused an extremely early flowering, delayed perianth and stamen filament development, small flowers, and significantly shorter Arabidopsis plants with little fruit. In the JcAG-overexpressing line, the homeotic transformation of sepals into pistillate organs was observed, and floral meristem and organ identity genes were regulated. This study provides insights into the JcAG's function and benefits to our knowledge of the underlying the genetic mechanisms related to floral sex differentiation in Jatropha.

Dissection of the Arabidopsis HUA-PEP gene activity reveals that ovule fate specification requires restriction of the floral A-function

Ovules are essential for sexual plant reproduction and seed formation, and are fundamental for agriculture. However, our understanding of the molecular mechanisms governing ovule development is far from complete. In Arabidopsis, ovule identity is determined by homeotic MADS-domain proteins that define the floral C- (AG) and D- (SHP1/SHP2, STK) functions. Pre-mRNA processing of these genes is critical and mediated by HUA-PEP activity, composed of genes encoding RNA-binding proteins. In strong hua-pep mutants, functional transcripts for C- and D-function genes are reduced, resulting in homeotic transformation of ovules. Thus, hua-pep mutants provide an unique sensitized background to study ovule morphogenesis when C- and D-functions are simultaneously compromised. We found that hua-pep ovules are morphologically sepaloid and show ectopic expression of the homeotic class-A gene AP1. Inactivation of AP1 or AP2 (A-function genes) in hua-pep mutants reduced homeotic conversions, rescuing ovule identity while promoting carpelloid traits in transformed ovules. Interestingly, increased AG dosage led to similar results. Our findings strongly suggest that HUA-PEP activity is required for correct C and D floral functions, which in turn prevents ectopic expression of class-A genes in ovules for their proper morphogenesis, evoking the classic A-C antagonism of the ABC model for floral organ development.

Comparative transcriptomic analysis reveals novel roles of transcription factors and hormones during the flowering induction and floral bud differentiation in sweet cherry trees (Prunus avium L. cv. Bing)

In sweet cherry trees, flowering is commercially important because the flowers, after fertilization, will generate the fruits. In P. avium, the flowering induction and flower organogensis are the first developmental steps towards flower formation and they occur within specialized organs known as floral buds during the summer, nine months before blooming. During this period the number of floral buds per tree and the bud fruitfulness (number of flowers per bud) are stablished affecting the potential yield of orchards and the plant architecture. The floral bud development is sensitive to any type of stress and the hotter and drier summers will interfere with this process and are calling for new adapted cultivars. A better understanding of the underlying molecular and hormonal mechanisms would be of help, but unlike the model plant Arabidopsis, very little is known about floral induction in sweet cherry. To explore the molecular mechanism of floral bud differentiation, high-throughput RNA sequencing was used to detect differences in the gene expression of P. avium floral buds at five differentiation stages. We found 2,982 differentially expressed genes during floral bud development. We identified genes associated with floral initiation or floral organ identity that appear to be useful biomarkers of floral development and several transcription factor families (ERF, MYB, bHLH, MADS-box and NAC gene family) with novel potential roles during floral transition in this species. We analyzed in deep the MADS-box gene family and we shed light about their key role during floral bud and organs development in P. avium. Furthermore, the hormonal-related signatures in the gene regulatory networks and the dynamic changes of absicic acid, zeatin and indolacetic acid contents in buds suggest an important role for these hormones during floral bud differentiation in sweet cherry. These data provide a rich source of novel informacion for functional and evolutionary studies about floral bud development in sweet cherry and new tools for biotechnology and breeding.

JCDB: a comprehensive knowledge base for Jatropha curcas, an emerging model for woody energy plants

Background

Jatropha curcas is an oil-bearing plant, and has seeds with high oil content (~ 40%). Several advantages, such as easy genetic transformation and short generation duration, have led to the emergence of J. curcas as a model for woody energy plants. With the development of high-throughput sequencing, the genome of Jatropha curcas has been sequenced by different groups and a mass of transcriptome data was released. How to integrate and analyze these omics data is crucial for functional genomics research on J. curcas.

Results

By establishing pipelines for processing novel gene identification, gene function annotation, and gene network construction, we systematically integrated and analyzed a series of J. curcas transcriptome data. Based on these data, we constructed a J. curcas database (JCDB), which not only includes general gene information, gene functional annotation, gene interaction networks, and gene expression matrices but also provides tools for browsing, searching, and downloading data, as well as online BLAST, the JBrowse genome browser, ID conversion, heatmaps, and gene network analysis tools.

Conclusions

JCDB is the most comprehensive and well annotated knowledge base for J. curcas. We believe it will make a valuable contribution to the functional genomics study of J. curcas. The database is accessible at http://jcdb.xtbg.ac.cn.

Phenotypic variation of floral organs in Malus using frequency distribution functions

BACKGROUND

Phenotypic diversity of floral organs plays an important role in plant systematic taxonomy and genetic variation studies. Previous research have focused on the direction of variation but disregarded its degree. Phenotypic variation (including directions and degrees) of 17 floral traits from wild to cultivated crabapples were explored by comparing their distributions and deviations in three different dimensions: floral organ number, size, and the shape.

RESULTS

Except for petal number, petal length / petal width, and sepal length / sepal width, the analyzed floral traits of cultivated crabapples all showed downward distributed box bodies in box plot analysis and left deviations of fitted curves in frequency distribution function analysis when compared to the wild, which revealed consistent variation directions of petaloid conversion (pistils or stamens → petals), size miniaturization (large → small), and shape narrowness (petal shape: circular → elliptic; sepal shape: triangular → lanceolate). However, only seven floral traits exhibited significant differences in box plot analysis, while all of the traits in frequency distribution function analysis were obviously offset. The variation degrees were quantitatively characterized by sizing traits > shaping traits > numbering traits and by horizontal dimensions > radial dimensions.

CONCLUSIONS

Frequency distribution function analysis was more sensitive than the box plot analysis, which constructed a theoretical basis for Malus flower type breeding and would provide a new quantitative method for future evaluation of floral variation among different groups of angiosperms at large.

Identification of candidate flowering and sex genes in white Guinea yam (D. rotundata Poir.) by SuperSAGE transcriptome profiling

Dioecy (distinct male and female individuals) and scarce to non-flowering are common features of cultivated yam (Dioscorea spp.). However, the molecular mechanisms underlying flowering and sex determination in Dioscorea are largely unknown. We conducted SuperSAGE transcriptome profiling of male, female and monoecious individuals to identify flowering and sex-related genes in white Guinea yam (D. rotundata), generating 20,236 unique tags. Of these, 13,901 were represented by a minimum of 10 tags. A total 88 tags were significantly differentially expressed in male, female and monoecious plants, of which 18 corresponded to genes previously implicated in flower development and sex determination in multiple plant species. We validated the SuperSAGE data with quantitative real-time PCR (qRT-PCR)-based analysis of the expression of three candidate genes.

We further investigated the flowering patterns of 1938 D. rotundata accessions representing diverse geographical origins over two consecutive years. Over 85% of accessions were either male or non-flowering, less than 15% were female, while monoecious plants were rare. Intensity of flowering varied between male and female plants, with the former flowering more abundantly than the latter. Candidate genes identified in this study can be targeted for further validation and to induce regular flowering in poor to non-flowering cultivars. Findings of the study provide important inputs for further studies aiming to overcome the challenge of flowering in yams and to improve efficiency of yam breeding.

Genome-wide analysis of spatiotemporal gene expression patterns during floral organ development in Brassica rapa

Flowering is a key agronomic trait that directly influences crop yield and quality and serves as a model system for elucidating the molecular basis that controls successful reproduction, adaptation, and diversification of flowering plants. Adequate knowledge of continuous series of expression data from the floral transition to maturation is lacking in Brassica rapa. To unravel the genome expression associated with the development of early small floral buds (< 2 mm; FB2), early large floral buds (2-4 mm; FB4), stamens (STs) and carpels (CPs), transcriptome profiling was carried out with a Br300K oligo microarray. The results showed that at least 6848 known nonredundant genes (30% of the genes of the Br300K) were differentially expressed during the floral transition from vegetative tissues to maturation. Functional annotation of the differentially expressed genes (DEGs) (fold change ≥ 5) by comparison with a close relative, Arabidopsis thaliana, revealed 6552 unigenes (4579 upregulated; 1973 downregulated), including 131 Brassica-specific and 116 functionally known floral Arabidopsis homologs. Additionally, 1723, 236 and 232 DEGs were preferentially expressed in the tissues of STs, FB2, and CPs. These DEGs also included 43 transcription factors, mainly AP2/ERF-ERF, NAC, MADS-MIKC, C2H2, bHLH, and WRKY members. The differential gene expression during flower development induced dramatic changes in activities related to metabolic processes (23.7%), cellular (22.7%) processes, responses to the stimuli (7.5%) and reproduction (1%). A relatively large number of DEGs were observed in STs and were overrepresented by photosynthesis-related activities. Subsequent analysis via semiquantitative RT-PCR, histological analysis performed with in situ hybridization of BrLTP1 and transgenic reporter lines (BrLTP promoter::GUS) of B. rapa ssp. pekinensis supported the spatiotemporal expression patterns. Together, these results suggest that a temporally and spatially regulated process of the selective expression of distinct fractions of the same genome leads to the development of floral organs. Interestingly, most of the differentially expressed floral transcripts were located on chromosomes 3 and 9. This study generated a genome expression atlas of the early floral transition to maturation that represented the flowering regulatory elements of Brassica rapa.

Evolution and functional divergence of MADS-box genes in Pyrus

MADS-box transcription factors widely regulate all aspects of plant growth including development and reproduction. Although the MADS-box gene family genes have been extensively characterized in many plants, they have not been studied in closely related species. In this study, 73 and 74 MADS-box genes were identified in European pear (Pyrus communis) and Chinese pear (Pyrus bretschneideri), respectively. Based on the phylogenetic relationship, these genes could be clustered into five groups (Mα, Mβ, Mr, MIKC^C, MIKC*) and the MIKC^C group was further categorized into 10 subfamilies. The distribution of MADS-box genes on each chromosome was significantly nonrandom. Thirty-seven orthologs, twenty-five PcpMADS (P. communis MADS-box) paralogs and nineteen PbrMADS (P. bretschneideri MADS-box) paralogs were predicted. Among these paralogous genes, two pairs arose from tandem duplications (TD), nineteen from segmental duplication (SD) events and twenty-three from whole genome duplication (WGD) events, indicating SD/WGD events led to the expansion of MADS-box gene family. The MADS-box genes expression profiles in pear fruits indicated functional divergence and neo-functionalization or sub-functionalization of some orthologous genes originated from a common ancestor. This study provided a useful reference for further analysis the mechanisms of species differentiation and biodiversity formation among closely related species.

Evolution and Diversification of FRUITFULL Genes in Solanaceae

Ecologically and economically important fleshy edible fruits have evolved from dry fruit numerous times during angiosperm diversification. However, the molecular mechanisms that underlie these shifts are unknown. In the Solanaceae there has been a major shift to fleshy fruits in the subfamily Solanoideae. Evidence suggests that an ortholog of FRUITFULL (FUL), a transcription factor that regulates cell proliferation and limits the dehiscence zone in the silique of Arabidopsis, plays a similar role in dry-fruited Solanaceae. However, studies have shown that FUL orthologs have taken on new functions in fleshy fruit development, including regulating elements of tomato ripening such as pigment accumulation. FUL belongs to the core eudicot euFUL clade of the angiosperm AP1/FUL gene lineage. The euFUL genes fall into two paralogous clades, euFULI and euFULII. While most core eudicots have one gene in each clade, Solanaceae have two: FUL1 and FUL2 in the former, and MBP10 and MBP20 in the latter. We characterized the evolution of the euFUL genes to identify changes that might be correlated with the origin of fleshy fruit in Solanaceae. Our analyses revealed that the Solanaceae FUL1 and FUL2 clades probably originated through an early whole genome multiplication event. By contrast, the data suggest that the MBP10 and MBP20 clades are the result of a later tandem duplication event. MBP10 is expressed at weak to moderate levels, and its atypical short first intron lacks putative transcription factor binding sites, indicating possible pseudogenization. Consistent with this, our analyses show that MBP10 is evolving at a faster rate compared to MBP20. Our analyses found that Solanaceae euFUL gene duplications, evolutionary rates, and changes in protein residues and expression patterns are not correlated with the shift in fruit type. This suggests deeper analyses are needed to identify the mechanism underlying the change in FUL ortholog function.

Transcriptome analysis provides insight into prickle development and its link to defense and secondary metabolism in Solanum viarum Dunal

Prickles are epidermal outgrowth found on the aerial surface of several terrestrial plants. Microscopic studies on prickles of S. viarum Dunal indicated a crucial role of glandular trichomes (GTs) in their development. A spontaneously obtained prickleless mutant showed normal epidermal GTs, but its downstream developmental process to prickle was perturbed. Thus, prickleless mutant offers an ideal opportunity to unveil molecular regulators working downstream to GTs in the prickle formation. Differential transcriptome analysis of epidermis of prickly and prickleless mutant revealed that expression of several defense regulators like ethylene, salicylic acid, PR-proteins, etc. were significantly down-regulated in prickleless mutant, provide an important link between defense and prickle development. It was also noteworthy that the expression of few essential development related TFs like MADS-box, R2R3-MYB, REM, DRL1, were also down-regulated in the stem, petioles, and leaves of prickleless mutant indicating their potential role in prickle development. Interestingly, the gene expression of terpenoid, steroid, flavonoid, glucosinolate, and lignin biosynthesis pathways were up-regulated in prickleless mutant. The biochemical and qRT-PCR analysis also confirmed metabolite elevation. These results indicated that the loss of prickle was compensated by elevated secondary metabolism in the prickleless mutant which played important role in the biotic and abiotic stress management.

Overexpression of the MADS-box gene K-domain increases the yield potential of blueberry

MADS-box genes play a significant role for plant flowering. Keratin-like (K) domains are involved in protein-to-protein interactions in the formation of the MIKC-type MADS-box domain proteins. In this study, the potential of utilizing the K domain of a Vaccinium corymbosum SOC1-like gene (VcSOC1K) was investigated to modulate expression of other blueberry MADS-box genes for increasing blueberry productivity. Chilled transgenic blueberry plants overexpressing the VcSOC1K showed a significant increase in the number of canes, floral buds, and flower and fruit clusters compared to chilled non-transgenic plants. Additionally, nonchilled transgenic plants flowered whereas nonchilled non-transgenic plants did not. Transgenic plants showed an increase in tolerance to high soil pH. Comparative transcriptome analysis of transgenic and non-transgenic leaves showed differential expression of 17% of the MADS-box genes identified in blueberry. These differentially expressed (DE) MADS-box genes were associated with genes related to plant flowering, phytohormones, and response to various biotic and abiotic stimuli. The phenotypic changes and the DE MADS-box genes caused by the overexpression of VcSOC1K not only reveal that the MADS-box genes are involved in chilling/vernalization-mediated flowering in blueberry but also demonstrated that the overexpression of the K domain can effectively modulate plant reproductive processes.

Isolation and Characterization of AGAMOUS-Like Genes Associated With Double-Flower Morphogenesis in Kerria japonica (Rosaceae)

Double-flower phenotype is more popular and attractive in garden and ornamental plants. There is great interest in exploring the molecular mechanisms underlying the double-flower formation for further breeding and selection. Kerria japonica, a commercial ornamental shrub of the Rosaceae family, is considered an excellent system to determine the mechanisms of morphological alterations, because it naturally has a single-flower form and double-flower variant with homeotic conversion of stamens into petals and carpels into leaf-like carpels. In this study, Sf-KjAG (AGAMOUS homolog of single-flower K. japonica) and Df-KjAG (AGAMOUS homolog of double-flower K. japonica) were isolated and characterized as two AGAMOUS (AG) homologs that occur strictly in single- and double-flower K. japonica, respectively. Our sequence comparison showed that Df-KjAG is derived from ectopic splicing with the insertion of a 2411 bp transposon-like fragment, which might disrupt mRNA accumulation and protein function, into intron 1. Ectopic expression analysis in Arabidopsis revealed that Sf-KjAG is highly conserved in specifying carpel and stamen identities. However, Df-KjAG did not show any putative C-class function in floral development. Moreover, yeast-two-hybrid assays showed that Sf-KjAG can interact with KjAGL2, KjAGL9, and KjAP1, whereas Df-KjAG has lost interactions with these floral identity genes. In addition, loss-of-function of Df-KjAG affected not only its own expression, but also that of other putative floral organ identity genes such as KjAGL2, KjAGL9, KjAP1, KjAP2, KjAP3, and KjPI. In conclusion, our findings suggest that double-flower formation in K. japonica can be attributed to Df-KjAG, which appears to be a mutant produced by the insertion of a transposon-like fragment in the normal AG homolog (Sf-KjAG) of single-flower K. japonica. Highlights:Sf-KjAG and Df-KjAG are different variations only distinguished by a transposon-like fragment insertion which lead to the evolutionary transformation from single-flower to double-flowers morphogenesis in Kerria japonica.

Deregulation of MADS-box transcription factor genes in a mutant defective in the WUSCHEL-LIKE HOMEOBOX gene EVERGREEN of Petunia hybrida

Angiosperm inflorescences develop in two fundamentally different ways. In monopodial plants, for example in Arabidopsis thaliana, the flowers are initiated as lateral appendages of a central indeterminate inflorescence meristem. In sympodial plants, flowers arise by terminal differentiation of the inflorescence meristem, while further inflorescence development proceeds from new sympodial meristems that are generated at the flank of the terminal flower. We have used the sympodial model species Petunia hybrida to investigate inflorescence development. Here, we describe a mutant, bonsai (bns), which is defective in flower formation, inflorescence branching, and control of meristem size. Detailed microscopic analysis revealed that bns meristems retain vegetative charateristics including spiral phyllotaxis. Consistent with a block in flower formation, bns mutants exhibit a deregulated expression of various MADS-box genes. Molecular analysis revealed that the bns mutant carries a transposon insertion in the previously described EVERGREEN (EVG) gene, which belongs to the WUSCHEL-LIKE HOMEOBOX (WOX) transcription factor gene family. EVG falls in the WOX9 subfamily, which has diverse developmental functions in angiosperms. The comparison of WOX9 orthologues in five model species for flowering shows that these genes play functionally divergent roles in monopodial and sympodial plants, indicating that the WOX9 regulatory node may have played an important role in the evolution of shoot architecture.

TM8 represses developmental timing in Nicotiana benthamiana and has functionally diversified in angiosperms

BACKGROUND

MADS-box genes are key regulators of plant reproductive development and members of most lineages of this gene family have been extensively studied. However, the function and diversification of the ancient TM8 lineage remains elusive to date. The available data suggest a possible function in flower development in tomato and fast evolution through numerous gene loss events in flowering plants.

RESULTS

We show the broad conservation of TM8 within angiosperms and find that in contrast to other MADS-box gene lineages, no gene duplicates have been retained after major whole genome duplication events. Through knock-down of NbTM8 by virus induced gene silencing in Nicotiana benthamiana, we show that NbTM8 represses miR172 together with another MADS-box gene, SHORT VEGETATIVE PHASE (NbSVP). In the closely related species Petunia hybrida, PhTM8 is not expressed under the conditions we investigated and consistent with this, a knock-out mutant did not show a phenotype. Finally, we generated transgenic tomato plants in which TM8 was silenced or ectopically expressed, but these plants did not display a clear phenotype. Therefore, no clear function could be confirmed for Solanum lycopersium.

CONCLUSIONS

While the presence of TM8 is generally conserved, it remains difficult to propose a general function in angiosperms. Based on all the available data to date, supplemented with our own results, TM8 function seems to have diversified quickly throughout angiosperms and acts as repressor of miR172 in Nicotiana benthamiana, together with NbSVP.

Molecular cloning, characterization and expression analysis of MADS-box genes associated with reproductive development in Momordica dioica Roxb

The repertoire and functions of MADS-box family transcription factors (TFs) largely remains unexplored with respect to floral organogenesis of Momordica dioica Roxb. Degenerative PCR followed by rapid amplification of cDNA ends was employed in the present study to clone and characterize 17 MADS-box genes (designated as MdMADS01 to MdMADS17) from the floral buds of M. dioica. The cloned genes were clustered into three subgroups (11 MIKC^C, 4 MIKC* and 2 Mα) based on phylogenetic relationships with the MADS-box genes from Cucumis sativus, Cucumis melo and Arabidopsis thaliana. Southern hybridization showed that all the isolated genes were represented by single copy locus in the M. dioica genome. Gene structure analysis revealed 1-8 exons in MdMADS-box genes with the number of exons in MIKC greatly exceeding from that in M-type genes. Motif elicitation of the MdMADS-box genes indicated the presence of additional domains with MIKC type, suggesting that they had more complex structures. Expression analysis of MdMADS genes in six M. dioica transcriptome suggested that, 11 MIKC^C-type genes are associated with floral homeotic functions, 4 MIKC*-type genes (MdMADS12 to MdMADS15) controlled the growth of male gametophyte, while the two M-type genes (MdMADS16 and MdMADS17) played significant role in female gametogenesis and seed development. Overall, these are the first set of MADS-box genes from M. dioica exhibiting a differential expression pattern during floral development. The results from this study will provide valuable information for further functional studies of candidate MADS-box genes in the sexual dimorphism of this economically important dioecious cucurbit.

Analysis of Transcriptional Responses of the Inflorescence Meristems in Jatropha curcas Following Gibberellin Treatment

Jatropha curcas L. seeds an oilseed plant with great potential for biodiesel production. However, low seed yield, which was limited by its lower female flowers, was a major drawback for its utilization. Our previous study found that the flower number and female-to-male ratio were increased by gibberellin treatment. Here, we compared the transcriptomic profiles of inflorescence meristem at different time points after gibberellic acid A3 (GA₃) treatment. The present study showed that 951 differentially expressed genes were obtained in response to gibberellin treatment, compared with control samples. The 6-h time point was an important phase in the response to exogenous gibberellin. Furthermore, the plant endogenous gibberellin, auxin, ethylene, abscisic acid, and brassinolide-signaling transduction pathways were repressed, whereas the genes associated with cytokinin and jasmonic acid signaling were upregulated for 24-h time point following GA₃ treatment. In addition, the floral meristem determinacy genes (JcLFY, JcSOC1) and floral organ identity genes (JcAP3, JcPI, JcSEP1-3) were significantly upregulated, but their negative regulator (JcSVP) was downregulated after GA₃ treatment. Moreover, the effects of phytohormone, which was induced by exogenous plant growth regulator, mainly acted on the female floral differentiation process. To the best of our knowledge, this data is the first comprehensive analysis of the underlying transcriptional response mechanism of floral differentiation following GA₃ treatment in J. curcas, which helps in engineering high-yielding varieties of Jatropha.

How to Evolve a Perianth: A Review of Cadastral Mechanisms for Perianth Identity

The flower of angiosperms is considered to be a major evolutionary innovation that impacted the whole biome. In particular, two properties of the flower are classically linked to its ecological success: bisexuality and a differentiated perianth with sepals and petals. Although the molecular basis for floral organ identity is well understood in extant species and summarized in the famous ABC model, how perianth identity appeared during evolution is still unknown. Here we propose that cadastral mechanisms that maintain reproductive organ identities to the center of the flower could have supported perianth evolution. In particular, repressing B- and C-class genes expression toward the inner whorls of the flower, is a key process to isolate domains with sepal and petal identity in the outer whorls. We review from the literature in model species the diverse regulators that repress B- and C-class genes expression to the center of the flower. This review highlights the existence of both unique and conserved repressors between species, and possible candidates to investigate further in order to shed light on perianth evolution.

Transcriptome profile analysis reveals the regulation mechanism of floral sex differentiation in Jatropha curcas L

The seeds of Jatropha curcas contain a high percentage of biodiesel. However, low seed yield which was limited by its poor female flowers was a bottleneck for its utilization. Here, we compared the transcriptomic profiles of five different samples during floral sex differentiation stages using Illumina Hiseq 4000. Our results showed that hundreds of differentially expressed genes (DEGs) were detected in floral sex initiation period, but thousands of DEGs were involved in the stamens and ovules development process. Moreover, the DEGs were mainly shown up-regulation in male floral initiation, but mainly down-regulation in female floral initiation. Male floral initiation was associated with the flavonoid biosynthesis pathway while female floral initiation was related to the phytohormone signal transduction pathway. Cytokinin (CTK) signaling triggered the initiation of female floral primordium, thereafter other phytohormones co-promoted the female floral development. In addition, the floral organ identity genes played important roles in floral sex differentiation process and displayed a general conservation of the ABCDE model in J. curcas. To the best of our knowledge, this data is the first comprehensive analysis of the underlying regulatory mechanism and the related genes during floral sex differentiation in J. curcas, which help in engineering high-yielding varieties of J. curcas.

Flower Development and Sex Determination between Male and Female Flowers in Vernicia fordii

Vernicia fordii is a monoecious and diclinous species with male and female flowers on the same inflorescence. Low female to male flower ratio is one of the main reasons for low yield in this species. However, little is known of its floral development and sex determination. Here, according to the results of scanning electron microscopy and histological analysis, the floral development of V. fordii was divided into 12 stages and the first morphological divergence between the male and female flowers was found to occur at stage 7. The male flowers are always unisexual, but the female flowers present bisexual characteristics, with sterile stamen (staminode) restricted to pre-meiosis of mother sporogenous cells and cell death occurring at later development stages. To further elucidate the molecular mechanism underling sex determination at the divergence stage for male and female flowers, comparative transcriptome analysis was performed. In total, 56,065 unigenes were generated and 608 genes were differentially expressed between male and female flowers, among which 310 and 298 DEGs (differentially expressed genes) showed high expression levels in males and females, respectively. The transcriptome data showed that the sexual dimorphism of female flowers was affected by jasmonic acid, transcription factors, and some genes related to the floral meristem activity. Ten candidate genes showed consistent expression in the qRT-PCR validation and DEGs data. In this study, we provide developmental characterization and transcriptomic information for better understanding of the development of unisexual flowers and the regulatory networks underlying the mechanism of sex determination in V. fordii, which would be helpful in the molecular breeding of V. fordii to improve the yield output.

Phytoplasma-conserved phyllogen proteins induce phyllody across the Plantae by degrading floral MADS domain proteins

ABCE-class MADS domain transcription factors (MTFs) are key regulators of floral organ development in angiosperms. Aberrant expression of these genes can result in abnormal floral traits such as phyllody. Phyllogen is a virulence factor conserved in phytoplasmas, plant pathogenic bacteria of the class Mollicutes. It triggers phyllody in Arabidopsis thaliana by inducing degradation of A- and E-class MTFs. However, it is still unknown whether phyllogen can induce phyllody in plants other than A. thaliana, although phytoplasma-associated phyllody symptoms are observed in a broad range of angiosperms. In this study, phyllogen was shown to cause phyllody phenotypes in several eudicot species belonging to three different families. Moreover, phyllogen can interact with MTFs of not only angiosperm species including eudicots and monocots but also gymnosperms and a fern, and induce their degradation. These results suggest that phyllogen induces phyllody in angiosperms and inhibits MTF function in diverse plant species.

Evolutionary Analysis of MIKCc-Type MADS-Box Genes in Gymnosperms and Angiosperms

MIKC^c-type MADS-box genes encode transcription factors that control floral organ morphogenesis and flowering time in flowering plants. Here, in order to determine when the subfamilies of MIKC^c originated and their early evolutionary trajectory, we sampled and analyzed the genomes and large-scale transcriptomes representing all the orders of gymnosperms and basal angiosperms. Through phylogenetic inference, the MIKC^c-type MADS-box genes were subdivided into 14 monophyletic clades. Among them, the gymnosperm orthologs of AGL6, SEP, AP1, GMADS, SOC1, AGL32, AP3/PI, SVP, AGL15, ANR1, and AG were identified. We identified and characterized the origin of a novel subfamily GMADS within gymnosperms but lost orthologs in monocots and Brassicaceae. ABCE model prototype genes were relatively conserved in terms of gene number in gymnosperms, but expanded in angiosperms, whereas SVP, SOC1, and GMADS had dramatic expansions in gymnosperms but conserved in angiosperms. Our results provided the most detailed evolutionary history of all MIKC^c gene clades in gymnosperms and angiosperms. We proposed that although the near complete set of MIKC^c genes had evolved in gymnosperms, the duplication and expressional transition of ABCE model MIKC^c genes in the ancestor of angiosperms triggered the first flower.

The Rosa chinensis cv. Viridiflora Phyllody Phenotype Is Associated with Misexpression of Flower Organ Identity Genes

Phyllody is a flower abnormality in which leaf-like structures replace flower organs in all whorls. Here, we investigated the origin and the molecular mechanism of phyllody phenotype in Rosa chinensis cv. Viridiflora, an ancient naturally occurring Chinese mutant cultivar. Reciprocal grafting experiments and microscopy analyses, demonstrated that the phyllody phenotype in Viridiflora is not associated with phytoplasmas infection. Transcriptome comparisons by the mean of RNA-Seq identified 672 up-regulated and 666 down-regulated genes in Viridiflora compared to its closely related genotype R. chinensis cv. Old Blush. A fraction of these genes are putative homologs of genes known to be involved in flower initiation and development. We show that in flower whorl 2 of Viridiflora, a down-regulation of the floral organ identity genes RcPISTILLATA (RcPI), RcAPETALA3 (RcAP3) and RcSEPALLATA3 (RcSEP3), together with an up-regulation of the putative homolog of the gene SUPPRESSOR of OVEREXPRESSION of CONSTANS1 (RcSOC1) are likely at the origin of the loss of petal identity and leaf-like structures formation. In whorl 3 of Viridiflora, ectopic expression of RcAPETALA2 (RcAP2) along with the down regulation of RcPI, RcAP3, and RcSEP3 is associated with loss of stamens identity and leaf-like structures formation. In whorl 4, the ectopic expression of RcAP2 associated with a down-regulation of RcSEP3 and of the C-class gene RcAGAMOUS correlate with loss of pistil identity. The latter also suggested the antagonist effect between the A and C class genes in the rose. Together, these data suggest that modified expression of the ABCE flower organ identity genes is associated with the phyllody phenotype in the rose Viridiflora and that these genes are important for normal flower organs development.

The wheat Phs-A1 pre-harvest sprouting resistance locus delays the rate of seed dormancy loss and maps 0.3 cM distal to the PM19 genes in UK germplasm

The precocious germination of cereal grains before harvest, also known as pre-harvest sprouting, is an important source of yield and quality loss in cereal production. Pre-harvest sprouting is a complex grain defect and is becoming an increasing challenge due to changing climate patterns. Resistance to sprouting is multi-genic, although a significant proportion of the sprouting variation in modern wheat cultivars is controlled by a few major quantitative trait loci, including Phs-A1 in chromosome arm 4AL. Despite its importance, little is known about the physiological basis and the gene(s) underlying this important locus. In this study, we characterized Phs-A1 and show that it confers resistance to sprouting damage by affecting the rate of dormancy loss during dry seed after-ripening. We show Phs-A1 to be effective even when seeds develop at low temperature (13 °C). Comparative analysis of syntenic Phs-A1 intervals in wheat and Brachypodium uncovered ten orthologous genes, including the Plasma Membrane 19 genes (PM19-A1 and PM19-A2) previously proposed as the main candidates for this locus. However, high-resolution fine-mapping in two bi-parental UK mapping populations delimited Phs-A1 to an interval 0.3 cM distal to the PM19 genes. This study suggests the possibility that more than one causal gene underlies this major pre-harvest sprouting locus. The information and resources reported in this study will help test this hypothesis across a wider set of germplasm and will be of importance for breeding more sprouting resilient wheat varieties.

Whole-transcriptome analysis of differentially expressed genes in the ray florets and disc florets of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most important global cut flower and pot plants, and has been cultivated worldwide. However, limited genomic resources are available and the molecular mechanisms involved in the two morphologically distinct floret developmental cycles in chrysanthemum remain unclear.

RESULTS

The transcriptomes of chrysanthemum ray florets, disc florets and leaves were sequenced using Illumina paired-end sequencing technology. In total, 16.9 G reads were assembled into 93,138 unigenes with an average length of 738 bp, of which 44,364 unigenes showed similarity to known proteins in the Swissprot or NCBI non-redundant protein databases. Additionally, 26,320, 22,304 and 13,949 unigenes were assigned to 54 gene ontology (GO) categories, 25 EuKaryotic Orthologous Groups (KOG) categories, and 280 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. A total of 1863 differentially expressed genes (DEGs) (1210 up-regulated and 653 down-regulated) were identified between ray florets and disc florets, including genes encoding transcription factors and protein kinases. GO and KEGG pathway enrichment analyses were performed on the DEGs to identify differences in the biological processes and pathways between ray florets and disc florets. The important regulatory genes controlling flower development and flower organ determination, as well as important functional genes in the anthocyanin biosynthetic pathway, were identified, of which two leucoanthocyanidin dioxygenase-encoding genes showed specific expression in ray florets. Lastly, reverse transcription quantitative PCR was conducted to validate the DEGs identified in our study.

CONCLUSIONS

Comparative transcriptome analysis revealed significant differences in patterns of gene expression and signaling pathways between ray florets and disc florets in Chrysanthemum morifolium. This study provided the first step to understanding the molecular mechanism of the differential development of ray florets and disc florets in chrysanthemum, and also provided valuable genomic resources for candidate genes applicable for the breeding of novel varieties in chrysanthemum.

Comprehensive transcriptome analysis of Crocus sativus for discovery and expression of genes involved in apocarotenoid biosynthesis

Background

Crocus sativus stigmas form rich source of apocarotenoids like crocin, picrocrocin and saffranal which besides imparting color, flavour and aroma to saffron spice also have tremendous pharmacological properties. Inspite of their importance, the biosynthetic pathway of Crocus apocarotenoids is not fully elucidated. Moreover, the mechanism of their stigma specific accumulation remains unknown. Therefore, deep transcriptome sequencing of Crocus stigma and rest of the flower tissue was done to identify the genes and transcriptional regulators involved in the biosynthesis of these compounds.

Results

Transcriptome of stigma and rest of the flower tissue was sequenced using Illumina Genome Analyzer IIx platform which generated 64,604,402 flower and 51,350,714 stigma reads. Sequences were assembled de novo using trinity resulting in 64,438 transcripts which were classified into 32,204 unigenes comprising of 9853 clusters and 22,351 singletons. A comprehensive functional annotation and gene ontology (GO) analysis was carried out. 58.5 % of the transcripts showed similarity to sequences present in public databases while rest could be specific to Crocus. 5789 transcripts showed similarity to transcription factors representing 76 families out of which Myb family was most abundant. Many genes involved in carotenoid/apocarotenoid pathway were identified for the first time in this study which includes zeta-carotene isomerase and desaturase, carotenoid isomerase and lycopene epsilon-cyclase. GO analysis showed that the predominant classes in biological process category include metabolic process followed by cellular process and primary metabolic process. KEGG mapping analysis indicated that pathways involved in ribosome, carbon and starch and sucrose metabolism were highly represented. Differential expression analysis indicated that key carotenoid/apocarotenoid pathway genes including phytoene synthase, phytoene desaturase and carotenoid cleavage dioxygenase 2 are enriched in stigma thereby providing molecular proof for stigma to be the site of apocarotenoid biosynthesis.

Conclusions

This data would provide a rich source for understanding the carotenoid/apocarotenoid metabolism in Crocus. The database would also help in investigating many questions related to saffron biology including flower development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1894-5) contains supplementary material, which is available to authorized users.

A genome-wide analysis of MADS-box genes in peach [Prunus persica (L.) Batsch]

BACKGROUND

MADS-box genes encode a family of eukaryotic transcription factors distinguished by the presence of a highly-conserved ~58 amino acid DNA-binding and dimerization domain (the MADS-box). The central role played by MADS-box genes in peach endodormancy regulation led us to examine this large gene family in more detail. We identified the locations and sequences of 79 MADS-box genes in peach, separated them into established subfamilies, and broadly surveyed their tissue-specific and dormancy-induced expression patterns using next-generation sequencing. We then focused on the dormancy-related SVP/AGL24 and FLC subfamilies, comparing their numbers and phylogenetic relationships with those of other sequenced woody perennial genomes.

RESULTS

We identified 79 MADS-box genes distributed across all eight peach chromosomes and frequently located in clusters of two or more genes. They encode proteins with a mean length of 248 ± 72 amino acids and include representatives from most of the thirteen Type II (MIKC) subfamilies, as well as members of the Type I Mα, Mβ, and Mγ subfamilies. Most Type I genes were present in species-specific monophyletic lineages, and their expression in the peach sporophyte was low or absent. Most Type II genes had Arabidopsis orthologs and were expressed at much higher levels throughout vegetative and fruit tissues. During short-day-induced growth cessation, seven Type II genes from the SVP/AGL24, AGL17, and SEP subfamilies showed significant changes in expression. Phylogenetic analyses indicated that multiple, independent expansions have taken place within the SVP/AGL24 and FLC lineages in woody perennial species.

CONCLUSIONS

Most Type I genes appear to have arisen through tandem duplications after the divergence of the Arabidopsis and peach lineages, whereas Type II genes appear to have increased following whole genome duplication events. An exception to the latter rule occurs in the FLC and SVP/AGL24 Type II subfamilies, in which species-specific tandem duplicates have been retained in a number of perennial species. These subfamilies comprise part of a genetic toolkit that regulates endodormancy transitions, but phylogenetic and expression data suggest that individual orthologs may not function identically across all species.

Characterization of TM8, a MADS-box gene expressed in tomato flowers

BACKGROUND

The identity of flower organs is specified by various MIKC MADS-box transcription factors which act in a combinatorial manner. TM8 is a MADS-box gene that was isolated from the floral meristem of a tomato mutant more than twenty years ago, but is still poorly known from a functional point of view in spite of being present in both Angiosperms and Gymnosperms, with some species harbouring more than one copy of the gene. This study reports a characterization of TM8 that was carried out in transgenic tomato plants with altered expression of the gene.

RESULTS

Tomato plants over-expressing either TM8 or a chimeric repressor form of the gene (TM8:SRDX) were prepared. In the TM8 up-regulated plants it was possible to observe anomalous stamens with poorly viable pollen and altered expression of several floral identity genes, among them B-, C- and E-function ones, while no apparent morphological modifications were visible in the other whorls. Oblong ovaries and fruits, that were also parthenocarpic, were obtained in the plants expressing the TM8:SRDX repressor gene. Such ovaries showed modified expression of various carpel-related genes. No apparent modifications could be seen in the other flower whorls. The latter plants had also epinastic leaves and malformed flower abscission zones. By using yeast two hybrid assays it was possible to show that TM8 was able to interact in yeast with MACROCALIX.

CONCLUSIONS

The impact of the ectopically altered TM8 expression on the reproductive structures suggests that this gene plays some role in the development of the tomato flower. MACROCALYX, a putative A-function MADS-box gene, was expressed in all the four whorls of fully developed flowers, and showed quantitative variations that were opposite to those of TM8 in the anomalous stamens and ovaries. Since the TM8 protein interacted in vitro only with the A-function MADS-box protein MACROCALYX, it seems that for the correct differentiation of the tomato reproductive structures possible interactions between TM8 and MACROCALYX proteins might be important.

The role of WOX genes in flower development

BACKGROUND

WOX (Wuschel-like homeobOX) genes form a family of plant-specific HOMEODOMAIN transcription factors, the members of which play important developmental roles in a diverse range of processes. WOX genes were first identified as determining cell fate during embryo development, as well as playing important roles in maintaining stem cell niches in the plant. In recent years, new roles have been identified in plant architecture and organ development, particularly at the flower level.

SCOPE

In this review, the role of WOX genes in flower development and flower architecture is highlighted, as evidenced from data obtained in the last few years. The roles played by WOX genes in different species and different flower organs are compared, and differential functional recruitment of WOX genes during flower evolution is considered.

CONCLUSIONS

This review compares available data concerning the role of WOX genes in flower and organ architecture among different species of angiosperms, including representatives of monocots and eudicots (rosids and asterids). These comparative data highlight the usefulness of the WOX gene family for evo-devo studies of floral development.

Two euAGAMOUS genes control C-function in Medicago truncatula

C-function MADS-box transcription factors belong to the AGAMOUS (AG) lineage and specify both stamen and carpel identity and floral meristem determinacy. In core eudicots, the AG lineage is further divided into two branches, the euAG and PLE lineages. Functional analyses across flowering plants strongly support the idea that duplicated AG lineage genes have different degrees of subfunctionalization of the C-function. The legume Medicago truncatula contains three C-lineage genes in its genome: two euAG genes (MtAGa and MtAGb) and one PLENA-like gene (MtSHP). This species is therefore a good experimental system to study the effects of gene duplication within the AG subfamily. We have studied the respective functions of each euAG genes in M. truncatula employing expression analyses and reverse genetic approaches. Our results show that the M. truncatula euAG- and PLENA-like genes are an example of subfunctionalization as a result of a change in expression pattern. MtAGa and MtAGb are the only genes showing a full C-function activity, concomitant with their ancestral expression profile, early in the floral meristem, and in the third and fourth floral whorls during floral development. In contrast, MtSHP expression appears late during floral development suggesting it does not contribute significantly to the C-function. Furthermore, the redundant MtAGa and MtAGb paralogs have been retained which provides the overall dosage required to specify the C-function in M. truncatula.

Genome-wide identification, characterisation and expression analysis of the MADS-box gene family in Prunus mume

MADS-box genes encode transcription factors that play crucial roles in plant development, especially in flower and fruit development. To gain insight into this gene family in Prunus mume, an important ornamental and fruit plant in East Asia, and to elucidate their roles in flower organ determination and fruit development, we performed a genome-wide identification, characterisation and expression analysis of MADS-box genes in this Rosaceae tree. In this study, 80 MADS-box genes were identified in P. mume and categorised into MIKC, Mα, Mβ, Mγ and Mδ groups based on gene structures and phylogenetic relationships. The MIKC group could be further classified into 12 subfamilies. The FLC subfamily was absent in P. mume and the six tandemly arranged DAM genes might experience a species-specific evolution process in P. mume. The MADS-box gene family might experience an evolution process from MIKC genes to Mδ genes to Mα, Mβ and Mγ genes. The expression analysis suggests that P. mume MADS-box genes have diverse functions in P. mume development and the functions of duplicated genes diverged after the duplication events. In addition to its involvement in the development of female gametophytes, type I genes also play roles in male gametophytes development. In conclusion, this study adds to our understanding of the roles that the MADS-box genes played in flower and fruit development and lays a foundation for selecting candidate genes for functional studies in P. mume and other species. Furthermore, this study also provides a basis to study the evolution of the MADS-box family.

MADS reloaded: evolution of the AGAMOUS subfamily genes

AGAMOUS subfamily proteins are encoded by MADS-box family genes. They have been shown to play key roles in the determination of reproductive floral organs such as stamens, carpels and ovules. However, they also play key roles in ensuring a fixed number of floral organs by controlling floral meristem determinacy. Recently, an enormous amount of sequence data for nonmodel species have become available together with functional data on AGAMOUS subfamily members in many species. Here, we give a detailed overview of the most important information about this interesting gene subfamily and provide new insights into its evolution.