Integrated metabolomics and transcriptomics analysis highlight key pathways involved in the somatic embryogenesis of Darjeeling tea

Darjeeling tea is a globally renowned beverage, which faces numerous obstacles in sexual reproduction, such as self-incompatibility, poor seed germination, and viability, as well as issues with vegetative propagation. Somatic embryogenesis (SE) is a valuable method for rapid clonal propagation of Darjeeling tea. However, the metabolic regulatory mechanisms underlying SE in Darjeeling tea remain largely unknown. To address this, we conducted an integrated metabolomics and transcriptomics analysis of embryogenic callus (EC), globular embryo (GE), and heart-shaped embryo (HE).

Integrated metabolome and transcriptome analysis revealed enriched metabolic pathways, including auxin biosynthesis and signal transduction, brassinosteroid, gibberellin, phenylpropanoid biosynthesis, amino acids metabolism, and transcription factors (TFs) during SE in Darjeeling tea. Notably, EC displayed lower endogenous IAA levels, conducive to maintaining differentiation, while higher IAA concentration in GE and HE was crucial for preserving embryo identity. Additionally, a negative correlation between bioactive gibberellin A4 (GA4) and IAA was observed, impacting callus growth in EC. The high accumulation of chlorogenic acid, a phenolic compound, might contribute to the low success rate in GE and HE formation in Darjeeling tea. TFs such as BBM1, LEC1, FUS3, LEA, WOX3, and WOX11 appeared to regulate gene expression, influencing SE in Darjeeling tea.

The integrated analyses showed that various genes and metabolites involved in the phenylpropanoid pathway, auxin biosynthesis pathway, gibberellin, brassinosteroid and amino acids biosynthesis pathways were differentially enriched in EC, GE, and HE. Our results revealed that despite highly up-regulated auxin biosynthesis genes YUC1, TAR1 and AAO1 in EC, endogenous indole-3-acetic acid (IAA) was significantly lower in EC than GE and HE. However, bioactive Gibberellin A4 displayed higher accumulation in EC. We also found higher BABY BOOM (BBM) and Leafy cotyledon1 (LEC1) gene expression in GE along with high accumulation of castasterone, a brassinosteroid. Total flavonoids and phenolics levels were elevated in GE and HE compared to EC, especially the phenolic compound chlorogenic acid was highly accumulated in GE. Conclusions: Integrated metabolome and transcriptome analysis revealed enriched metabolic pathways, including auxin biosynthesis and signal transduction, brassinosteroid, gibberellin, phenylpropanoid biosynthesis, amino acids metabolism, and transcription factors (TFs) during SE in Darjeeling tea. Notably, EC displayed lower endogenous IAA levels, conducive to maintaining differentiation, while higher IAA concentration in GE and HE was crucial for preserving embryo identity. Additionally, a negative correlation between bioactive gibberellin A4 (GA4) and IAA was observed, impacting callus growth in EC. The high accumulation of chlorogenic acid, a phenolic compound, might contribute to the low success rate in GE and HE formation in Darjeeling tea. TFs such as BBM1, LEC1, FUS3, LEA, WOX3, and WOX11 appeared to regulate gene expression, influencing SE in Darjeeling tea.