Transcriptome and carotenoid profiling of different varieties of Coffea arabica provides insights into fruit color formation

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.

Transcriptomic Analysis of Alternative Splicing Events during Different Fruit Ripening Stages of Coffea arabica L

To date, genomic and transcriptomic data on Coffea arabica L. in public databases are very limited, and there has been no comprehensive integrated investigation conducted on alternative splicing (AS). Previously, we have constructed and sequenced eighteen RNA-seq libraries of C. arabica at different ripening stages of fruit development. From this dataset, a total of 3824, 2445, 2564, 2990, and 3162 DSGs were identified in a comparison of different fruit ripening stages. The largest proportion of DSGs, approximately 65%, were of the skipped exon (SE) type. Biologically, 9 and 29 differentially expressed DSGs in the spliceosome pathway and carbon metabolism pathway, respectively, were identified. These DSGs exhibited significant variations, primarily in S1 vs. S2 and S5 vs. S6, and they involve many aspects of organ development, hormone transduction, and the synthesis of flavor components. Through the examination of research findings regarding the biological functions and biochemical pathways associated with DSGs and DEGs, it was observed that six DSGs significantly enriched in ABC transporters, namely, LOC113712394, LOC113726618, LOC113739972, LOC113725240, LOC113730214, and LOC113707447, were continually down-regulated at the fruit ripening stage. In contrast, a total of four genes, which were LOC113732777, LOC113727880, LOC113690566, and LOC113711936, including those enriched in the cysteine and methionine metabolism, were continually up-regulated. Collectively, our findings may contribute to the exploration of alternative splicing mechanisms for focused investigations of potential genes associated with the ripening of fruits in C. arabica.

A Systematic Analysis of the Correlation between Flavor Active Differential Metabolites and Multiple Bean Ripening Stages of Coffea arabica L

Coffee cherries contain a crucial flavor-precursor and chemical substances influencing roasted bean quality, yet limited knowledge exists on metabolite changes during cherry ripening. Our study identified 1078 metabolites, revealing 46 core differential metabolites using a KEGG pathway analysis. At the GF vs. ROF stage, amino acid synthesis dominated; ROF vs. BRF featured nucleotide catabolism; BRF vs. PRF exhibited glycoside and flavonoid synthesis; and PRF vs. PBF involved secondary metabolite synthesis and catabolism. The PRF stage emerged as the optimal cherry-harvesting period. A correlation analysis identified core differential metabolites strongly linked to taste indicators, suggesting their potential as taste markers. Notably, nucleotides and derivatives exhibited significant negative correlations with glycosides and flavonoids during ripening. This research systematically analyzed flavor and active substances in green coffee beans during cherry ripening, offering valuable insights into substance formation in Coffea arabica L.

Transcriptomics and Metabolomics Reveal the Critical Genes of Carotenoid Biosynthesis and Color Formation of Goji (Lycium barbarum L.) Fruit Ripening

Carotenoids in goji (Lycium barbarum L.) have excellent health benefits, but the underlying mechanism of carotenoid synthesis and color formation in goji fruit ripening is still unclear. The present study uses transcriptomics and metabolomics to investigate carotenoid biosynthesis and color formation differences in N1 (red fruit) and N1Y (yellow fruit) at three stages of ripening. Twenty-seven carotenoids were identified in N1 and N1Y fruits during the M1, M2, and M3 periods, with the M2 and M3 periods being critical for the difference in carotenoid and color between N1 and N1Y fruit. Weighted gene co-expression network analysis (WGCNA), gene trend analysis, and correlation analysis suggest that PSY1 and ZDS16 may be important players in the synthesis of carotenoids during goji fruit ripening. Meanwhile, 63 transcription factors (TFs) were identified related to goji fruit carotenoid biosynthesis. Among them, four TFs (CMB1-1, WRKY22-1, WRKY22-3, and RAP2-13-like) may have potential regulatory relationships with PSY1 and ZDS16. This work sheds light on the molecular network of carotenoid synthesis and explains the differences in carotenoid accumulation in different colored goji fruits.