Flower development in Coffea arabica L.: new insights into MADS-box genes

Hormonal crosstalk during the reproductive stage of Coffea arabica: interactions among gibberellin, abscisic acid, and ethylene

MAIN CONCLUSION

The application of gibberellin and abscisic acid in coffee plants resulted in increased floral bud formation and fruit production by regulating key genes involved in flowering and hormonal biosynthesis pathways. Despite ongoing efforts, understanding hormonal regulation in perennial and woody species with complex phenological cycles, such as Coffea arabica L., remains limited. Given the global importance of coffee, identifying the main regulators of reproductive development is crucial to guarantee production, especially in face of climate change. This study investigated the effects of gibberellin (GA) and abscisic acid (ABA) at different concentrations (5, 25 and 100 ppm) in the reproductive development of C. arabica. Phenological analyses, molecular identification of genes involved in GA and ABA biosynthesis, degradation, and signaling, as well as gene expression profiling in leaves and floral buds during floral induction and development, were conducted. Promoter analysis of CaFT, quantification of 1-aminocyclopropane-1-carboxylate (ACC), enzymatic activity of ACC oxidase (ACO), and ethylene content were also assessed. Results showed that GA irrespective of concentration and ABA at 25 ppm applied during the main period of floral induction (March) significantly increased the number of floral buds, with ABA also accelerating the development. Similarly, applying these regulators in plants with floral buds at more advanced stages (August) increased the number of floral buds and fruit production in the GA (5 and 100 ppm) and ABA (25 and 100 ppm) treatments. Phylogenetic and molecular analyses identified genes related to GA and ABA biosynthesis, degradation, and signaling in coffee plants. GA and ABA treatments affected the expression of genes related to floral induction and organ formation, such as CaDELLA in March, which may relate to the increased number of floral buds. Moreover, in August, plants treated with 5 and 100 ppm GA and 100 ppm ABA showed up-regulation of CaFT1 expression, likely due to the down-regulation of CaCO during this period. In addition to GA-ABA interactions, our results suggest that GA promotes ACC accumulation in leaves in August, which may act as a mobile signal transported to floral buds, where its conversion to ethylene could regulate anthesis, highlighting a GA-ACC-ethylene interaction in coffee flowering. However, no significant differences in ethylene biosynthesis were observed in March with the application of these hormones, underscoring the incipient role of ethylene during floral induction in coffee. These results suggest reciprocal regulation of floral development by GA-ABA pathways in a dose-dependent manner and interacting with other hormonal pathways such as the ethylene biosynthesis in leaves and floral buds. These findings provide new insights into the hormonal regulation of coffee flowering, guiding field practices and breeding programs to maximize coffee production.

Transcriptomic Analyses Reveal That Coffea arabica and Coffea canephora Have More Complex Responses under Combined Heat and Drought than under Individual Stressors

Increasing exposure to unfavorable temperatures and water deficit imposes major constraints on most crops worldwide. Despite several studies regarding coffee responses to abiotic stresses, transcriptome modulation due to simultaneous stresses remains poorly understood. This study unravels transcriptomic responses under the combined action of drought and temperature in leaves from the two most traded species: Coffea canephora cv. Conilon Clone 153 (CL153) and C. arabica cv. Icatu. Substantial transcriptomic changes were found, especially in response to the combination of stresses that cannot be explained by an additive effect. A large number of genes were involved in stress responses, with photosynthesis and other physiologically related genes usually being negatively affected. In both genotypes, genes encoding for protective proteins, such as dehydrins and heat shock proteins, were positively regulated. Transcription factors (TFs), including MADS-box genes, were down-regulated, although responses were genotype-dependent. In contrast to Icatu, only a few drought- and heat-responsive DEGs were recorded in CL153, which also reacted more significantly in terms of the number of DEGs and enriched GO terms, suggesting a high ability to cope with stresses. This research provides novel insights into the molecular mechanisms underlying leaf Coffea responses to drought and heat, revealing their influence on gene expression.

The floral development of the allotetraploid Coffea arabica L. correlates with a small RNA dynamic reprogramming

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression.

Comprehensive expression analyses of the ABCG subfamily reveal SvABCG17 as a potential transporter of lignin monomers in the model C4 grass Setaria viridis

Although several aspects of lignin metabolism have been extensively characterized, the mechanism(s) by which lignin monomers are transported across the plasma membrane remains largely unknown. Biochemical, proteomic, expression and co-expression analyses from several plant species support the involvement of active transporters, mainly those belonging to the ABC superfamily. Here, we report on the genome-wide characterization of the ABCG gene subfamily in the model C4 grass Setaria viridis and further identification of the members potentially involved in monolignol transport. A total of 48 genes encoding SvABCGs were found in the S. viridis genome, from which 21 SvABCGs were classified as full-size transporters and 27 as half-size transporters. Comprehensive analysis of the ABCG subfamily in S. viridis based on expression and co-expression analyses support a role for SvABCG17 in monolignol transport: (i) SvABCG17 is orthologous to AtABCG29, a monolignol transporter in Arabidopsis thaliana; (ii) SvABCG17 displays a similar expression profile to that of lignin biosynthetic genes in a set of different S. viridis tissues and along the elongating internode; (iii) SvABCG17 is highly co-expressed with lignin-related genes in a public transcriptomic database; (iv) SvABCG17displays particularly high expression in the top of the S. viridis elongating internode, a tissue undergoing active lignification; (v) SvABCG17 mRNA localization coincides with the histochemical pattern of lignin deposition; and (vi) the promoter of SvABCG17 is activated by secondary cell wall-associated transcription factors, especially by lignin-specific activators of the MYB family. Further studies might reveal further aspects of this potential monolignol transporter, including its real substrate specificity and whether it works redundantly with other ABC members during S. viridis lignification.

Expression of coffee florigen CaFT1 reveals a sustained floral induction window associated with asynchronous flowering in tropical perennials

The behavior of florigen(s) and environment-influenced regulatory pathways that control floral initiation in tropical perennials species with complex phenological cycles is poorly understood. Understanding the mechanisms underlying this process is important for food production in the face of climate change, thus, we used Coffea sp. L. (Rubiaceae) as a model to explore this issue. Homologs of FLOWERING LOCUS T (CaFT1) and environment-related regulators CONSTANS (CaCO), PHYTOCHROME INTERACTING FACTOR 4 (CaPIF4) and FLOWERING LOCUS C (CaFLC) were retrieved from coffee genomes and identified through phylogenetic analysis. Overexpression of CaFT1 in Arabidopsis caused early-flowering phenotype and yeast two hybrid studies indicated CaFT1 binding to bZIP floral regulator FD, which suggests that CaFT1 is a coffee florigen. Expression of CaFT1 and other floral regulators, together with carbohydrate analysis, were evaluated over one year using three contrasting genotypes, two C. arabica cultivars and C. canephora. All genotypes showed active and variable CaFT1 transcription from February until October, indicating the potential window for floral induction that reached a maximum in the cold period of June. CaCO expression, as expected, varied over a 24-hour day period and monthly with day length, whereas expression of temperature-responsive homologs, CaFLC and CaPIF4, did not correlate with temperature changes nor CaFT1 expression, suggesting alternative FT regulatory pathways in coffee. Based on our results, we suggest a continuum of floral induction that allows different starting points for floral activation, which explains developmental asynchronicity and prolonged anthesis events in tropical perennial species.

Multigenic regulation in the ethylene biosynthesis pathway during coffee flowering

Ethylene regulates different aspects of the plant's life cycle, such as flowering, and acts as a defense signal in response to environmental stresses. Changes induced by water deficit (WD) in gene expression of the main enzymes involved in ethylene biosynthesis, 1-aminocyclopropane-1-carboxylic acid synthase (ACS) and oxidase (ACO), are frequently reported in plants. In this study, coffee (Coffea arabica) ACS and ACO family genes were characterized and their expression profiles were analyzed in leaves, roots, flower buds, and open flowers from plants under well-watered (WW) and water deficit (WD) conditions. Three new ACS genes were identified. Water deficit did not affect ACS expression in roots, however soil drying strongly downregulated ACO expression, indicating a transcriptional constraint in the biosynthesis pathway during the drought that can suppress ethylene production in roots. In floral buds, ACO expression is water-independent, suggesting a higher mechanism of control in reproductive organs during the final flowering stages. Leaves may be the main sites for ethylene precursor (1-aminocyclopropane-1-carboxylic acid, ACC) production in the shoot under well-watered conditions, contributing to an increase in the ethylene levels required for anthesis. Given these results, we suggest a possible regulatory mechanism for the ethylene biosynthesis pathway associated with coffee flowering with gene regulation in leaves being a key point in ethylene production and ACO genes play a major regulatory role in roots and the shoots. This mechanism may constitute a regulatory model for flowering in other woody species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01235-y.

Genome-Wide Analyses of MADS-Box Genes in Humulus lupulus L. Reveal Potential Participation in Plant Development, Floral Architecture, and Lupulin Gland Metabolism

MADS-box transcription factors (TFs) are involved in multiple plant development processes and are most known during the reproductive transition and floral organ development. Very few genes have been characterized in the genome of Humulus lupulus L. (Cannabaceae), an important crop for the pharmaceutical and beverage industries. The MADS-box family has not been studied in this species yet. We identified 65 MADS-box genes in the hop genome, of which 29 encode type-II TFs (27 of subgroup MIKC^C and 2 MIKC*) and 36 type-I proteins (26 α, 9 β, and 1 γ). Type-II MADS-box genes evolved more complex architectures than type-I genes. Interestingly, we did not find FLOWERING LOCUS C (FLC) homologs, a transcription factor that acts as a floral repressor and is negatively regulated by cold. This result provides a molecular explanation for a previous work showing that vernalization is not a requirement for hop flowering, which has implications for its cultivation in the tropics. Analysis of gene ontology and expression profiling revealed genes potentially involved in the development of male and female floral structures based on the differential expression of ABC homeotic genes in each whorl of the flower. We identified a gene exclusively expressed in lupulin glands, suggesting a role in specialized metabolism in these structures. In toto, this work contributes to understanding the evolutionary history of MADS-box genes in hop, and provides perspectives on functional genetic studies, biotechnology, and crop breeding.

Crosstalk Between Ethylene and Abscisic Acid During Changes in Soil Water Content Reveals a New Role for 1-Aminocyclopropane-1- Carboxylate in Coffee Anthesis Regulation

Coffee (Coffea arabica L.) presents an asynchronous flowering regulated by an endogenous and environmental stimulus, and anthesis occurs once plants are rehydrated after a period of water deficit. We evaluated the evolution of Abscisic Acid (ABA), ethylene, 1-aminocyclopropane-1-carboxylate (ACC) content, ACC oxidase (ACO) activity, and expression analysis of the Lysine Histidine Transporter 1 (LHT1) transporter, in the roots, leaves, and flower buds from three coffee genotypes (C. arabica L. cv Oeiras, Acauã, and Semperflorens) cultivated under field conditions with two experiments. In a third field experiment, the effect of the exogenous supply of ACC in coffee anthesis was evaluated. We found an increased ACC level, low ACO activity, decreased level of ethylene, and a decreased level of ABA in all tissues from the three coffee genotypes in the re-watering period just before anthesis, and a high expression of the LHT1 in flower buds and leaves. The ethylene content and ACO activity decreased from rainy to dry period whereas the ABA content increased. A higher number of opened and G6 stage flower buds were observed in the treatment with exogenous ACC. The results showed that the interaction of ABA-ACO-ethylene and intercellular ACC transport among the leaves, buds, and roots in coffee favors an increased level of ACC that is most likely, involved as a modulator in coffee anthesis. This study provides evidence that ACC can play an important role independently of ethylene in the anthesis process in a perennial crop.

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

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

Pitfall Flower Development and Organ Identity of Ceropegia sandersonii (Apocynaceae-Asclepiadoideae)

Deceptive Ceropegia pitfall flowers are an outstanding example of synorganized morphological complexity. Floral organs functionally synergise to trap fly-pollinators inside the fused corolla. Successful pollination requires precise positioning of flies headfirst into cavities at the gynostegium. These cavities are formed by the corona, a specialized organ of corolline and/or staminal origin. The interplay of floral organs to achieve pollination is well studied but their evolutionary origin is still unclear. We aimed to obtain more insight in the homology of the corona and therefore investigated floral anatomy, ontogeny, vascularization, and differential MADS-box gene expression in Ceropegia sandersonii using X-ray microtomography, Light and Scanning Electronic Microscopy, and RT-PCR. During 10 defined developmental phases, the corona appears in phase 7 at the base of the stamens and was not found to be vascularized. A floral reference transcriptome was generated and 14 MADS-box gene homologs, representing all major MADS-box gene classes, were identified. B- and C-class gene expression was found in mature coronas. Our results indicate staminal origin of the corona, and we propose a first ABCDE-model for floral organ identity in Ceropegia to lay the foundation for a better understanding of the molecular background of pitfall flower evolution in Apocynaceae.

Homeologous regulation of Frigida-like genes provides insights on reproductive development and somatic embryogenesis in the allotetraploid Coffea arabica

Coffea arabica is an allotetraploid of high economic importance. C. arabica transcriptome is a combination of the transcripts of two parental genomes (C. eugenioides and C. canephora) that gave rise to the homeologous genes of the species. Previous studies have reported the transcriptional dynamics of C. arabica. In these reports, the ancestry of homeologous genes was identified and the overall regulation of homeologous differential expression (HDE) was explored. One of these genes is part of the FRIGIDA-like family (FRL), which includes the Arabidopsis thaliana flowering-time regulation protein, FRIGIDA (FRI). As nonfunctional FRI proteins give rise to rapid-cycling summer annual ecotypes instead of vernalization-responsive winter-annuals, allelic variation in FRI can modulate flowering time in A. thaliana. Using bioinformatics, genomic analysis, and the evaluation of gene expression of homeologs, we characterized the FRL gene family in C. arabica. Our findings indicate that C. arabica expresses 10 FRL homeologs, and that, throughout flower and fruit development, these genes are differentially transcribed. Strikingly, in addition to confirming the expression of FRL genes during zygotic embryogenesis, we detected FRL expression during direct somatic embryogenesis, a novel finding regarding the FRL gene family. The HDE profile of FRL genes suggests an intertwined homeologous gene regulation. Furthermore, we observed that FLC gene of C. arabica has an expression profile similar to that of CaFRL genes.

Cytological analysis of ginseng carpel development

Panax ginseng Meyer, commonly known as ginseng, is considered one of the most important herbs with pharmaceutical values due to the presence of ginsenosides and is cultivated for its highly valued root for medicinal purposes. Recently, it has been recognized that ginseng fruit contains high contents of triterpene such as ginsenoside Re as pharmaceutical compounds. However, it is unclear how carpel, the female reproductive tissue of flowers, is formed during the three-year-old growth before fruit is formed in ginseng plants. Here, we report P. ginseng carpel development at the cytological level, starting from the initial stage of ovule development to seed development. The carpel of P. ginseng is composed of two free stigmas, two free styles, and one epigynous bilocular ovary containing one ovule in each locule. Based on our cytological study, we propose that the female reproductive development in P. ginseng can be classified into seven stages: early phase of ovule development, megasporogenesis, megagametogenesis, pre-fertilization, fertilization, post-fertilization, and seed development. We also describe the correlation of the female and male gametophyte development and compare morphological differences in carpel development between ginseng and other higher plants. One unique feature for ginseng seed development is that it takes 40 days for the embryo to develop to the early torpedo stage and that the embryo is small relative to the seed size, which could be a feature of taxonomic importance. This study will provide an integral tool for the study of the reproductive development and breeding of P. ginseng.

Ectopic expression of Jatropha curcas APETALA1 (JcAP1) caused early flowering in Arabidopsis, but not in Jatropha

Jatropha curcas is a promising feedstock for biofuel production because Jatropha oil is highly suitable for the production of biodiesel and bio-jet fuels. However, Jatropha exhibits a low seed yield as a result of unreliable and poor flowering. APETALA1 (AP1) is a floral meristem and organ identity gene in higher plants. The flower meristem identity genes of Jatropha have not yet been identified or characterized. To better understand the genetic control of flowering in Jatropha, an AP1 homolog (JcAP1) was isolated from Jatropha. An amino acid sequence analysis of JcAP1 revealed a high similarity to the AP1 proteins of other perennial plants. JcAP1 was expressed in inflorescence buds, flower buds, sepals and petals. The highest expression level was observed during the early developmental stage of the flower buds. The overexpression of JcAP1 using the cauliflower mosaic virus (CaMV) 35S promoter resulted in extremely early flowering and abnormal flowers in transgenic Arabidopsis plants. Several flowering genes downstream of AP1 were up-regulated in the JcAP1-overexpressing transgenic plant lines. Furthermore, JcAP1 overexpression rescued the phenotype caused by the Arabidopsis AP1 loss-of-function mutant ap1-11. Therefore, JcAP1 is an ortholog of AtAP1, which plays a similar role in the regulation of flowering in Arabidopsis. However, the overexpression of JcAP1 in Jatropha using the same promoter resulted in little variation in the flowering time and floral organs, indicating that JcAP1 may be insufficient to regulate flowering by itself in Jatropha. This study helps to elucidate the function of JcAP1 and contributes to the understanding of the molecular mechanisms of flower development in Jatropha.