# Stage-specific metabolic divergence in flavonoid biosynthesis correlates with embryogenic capacity in rubber tree (Hevea brasiliensis)

**Authors:** Jia Miao, Xiao-Long Sun, Jin Liu, Ming-Chun Gui, Min Tang, Hai Tian, Wan-Yuan Shi, Ling Li

PMC · DOI: 10.3389/fpls.2026.1766162 · Frontiers in Plant Science · 2026-02-04

## TL;DR

The study finds that metabolic changes, especially in flavonoid production, differ between high and low embryogenic genotypes of rubber trees during specific stages of development.

## Contribution

The study identifies stage-specific metabolic divergence in flavonoid biosynthesis as a key factor influencing embryogenic capacity in rubber tree genotypes.

## Key findings

- Metabolic variation is primarily driven by developmental stage rather than genotype.
- Flavonoid biosynthesis is a key differentiating pathway during the callus-to-differentiation transition.
- High embryogenic genotypes show coordinated repression of flavonoid biosynthesis genes during differentiation.

## Abstract

Somatic embryogenesis (SE) is an essential propagation technology for Hevea brasiliensis, yet its application remains limited by the strong genotype dependence of embryogenic capacity.

To elucidate the metabolic basis of this variation, we conducted integrated metabolomic and transcriptomic analyses across four SE developmental stages in a high-embryogenic (HE) and a low-embryogenic (LE) genotype, including explants, induced callus, non-embryogenic / embryogenic callus, and cotyledonary embryos (HE-specific).

A total of 1,383 metabolites belonging to 11 major classes were identified, with flavonoids, phenolic acids, and amino acids being the predominant groups. PCA and hierarchical clustering revealed that metabolic variation was driven primarily by developmental stage rather than genotype. Differential metabolite profiling revealed strong stage specificity, with the callus-to-differentiation transition (LE-C vs. HE-EC) exhibiting the greatest metabolic divergence between genotypes. KEGG enrichment consistently highlighted flavonoid biosynthesis as a key differentiating pathway. Comparative analyses revealed a conserved-to-divergent pattern of metabolic regulation. During the explant-to-callus transition, both genotypes exhibited highly conserved flavonoid biosynthesis responses, with 67.5% of genes and 85.7% of metabolites showing concordant regulation (either both up-regulated or both down-regulated). In contrast, during the callus-to-differentiation transition, pronounced metabolic divergence emerged, with only 37.5% of genes and 6.7% of metabolites showing concordant regulation, and 11 flavonoid-related genes displaying opposite regulatory directions between genotypes. Notably, the HE genotype exhibited coordinated repression of CHS, CHI, F3H, UFGT, and anthocyanin biosynthesis, accompanied by decreased accumulation of naringenin and glycosylated flavonoids, along with an overall attenuation of dihydroflavonol accumulation. Conversely, the LE genotype maintained relatively active flavonoid biosynthesis and glycosylation, along with increased amino sugar and nucleotide sugar metabolism.

Our results provide comprehensive metabolomic evidence for stage-dependent metabolic reprogramming during SE in H. brasiliensis. The contrasting patterns of flavonoid metabolism between genotypes at the callus-to-differentiation transition—systematic downregulation in the HE genotype versus sustained activation in the LE genotype—are consistent with the hypothesis that a timely reallocation of metabolic flux from secondary to primary metabolism may favor somatic embryo development. This study identifies the callus-to-differentiation transition as a critical metabolic checkpoint and suggests flavonoid biosynthesis genes, particularly CHS and glycosyltransferases, as potential targets for improving SE efficiency in recalcitrant genotypes.

## Linked entities

- **Genes:** LYST (lysosomal trafficking regulator) [NCBI Gene 1130], Chi (Chip) [NCBI Gene 37837], F3H (flavanone 3-hydroxylase) [NCBI Gene 732548], UFGT (anthocyanidin 3-O-glucosyltransferase 2) [NCBI Gene 100233099]
- **Chemicals:** naringenin (PubChem CID 932), dihydroflavonol (PubChem CID 147806)
- **Species:** Hevea brasiliensis (taxon 3981)

## Full-text entities

- **Diseases:** CHS (MESH:D002249), SE (MESH:D013001), EC (MESH:D005955)
- **Chemicals:** luteolin (MESH:D047311), spermidine (MESH:D013095), tryptophan (MESH:D014364), epigallocatechin (MESH:C057580), 5-O-caffeoylshikimic acid (MESH:C036792), L-ascorbic acid (MESH:D001205), alkaloids (MESH:D000470), cellulose (MESH:D002482), anthocyanin (MESH:D000872), N-methyl-trans-4-hydroxy-L-proline (MESH:C000716210), flavonol (MESH:C041477), xanthine (MESH:D019820), hemicellulose (MESH:C007916), flavonoid (MESH:D005419), UDP-arabinose (MESH:C040817), 2,3-dihydroxy-3-methylpentanoic acid (MESH:C014491), Vomifoliol (MESH:C525026), auxin (MESH:D007210), ROS (MESH:D017382), HgCl2 (MESH:D008627), leucine (MESH:D007930), NAA (MESH:D009280), polyamine (MESH:D011073), apigenin-7-O-rutinoside (MESH:C111466), dihydromyricetin (MESH:C472036), water (MESH:D014867), ATP (MESH:D000255), glutathione (MESH:D005978), carotenoid (MESH:D002338), myricetin-3-O-rhamnoside (MESH:C529905), ursonic acid (MESH:C000657350), cytokinins (MESH:D003583), valine (MESH:D014633), zirconia (MESH:C028541), sucrose (MESH:D013395), phenolic acid (MESH:C017616), lipid (MESH:D008055), nucleotide (MESH:D009711), lithospermoside (MESH:C070237), quercetin-3-O-galactoside (MESH:C021304), terpenoids (MESH:D013729), UDP-xylose (MESH:D014540), spermine (MESH:D013096), apigenin (MESH:D047310), naringenin (MESH:C005273), malonyl-CoA (MESH:D008316), dihydrochrysin (MESH:C016063), 1,3,5-benzenetriol (MESH:D010696), genistein-8-C-glucoside (MESH:C415443), proanthocyanidins (MESH:D044945), nitrogen (MESH:D009584), isoleucine (MESH:D007532), quercetin (MESH:D011794), Amino acid (MESH:D000596), UDP-glucose (MESH:D014532), DA (MESH:C025953), acetonitrile (MESH:C032159), myricitrin (MESH:C008577), butin (MESH:C051437), CoA (MESH:D003065)
- **Species:** Picea (genus) [taxon 3328], Hevea brasiliensis (jebe, species) [taxon 3981], Vitis vinifera (wine grape, species) [taxon 29760], Larix kaempferi (karamatsu, species) [taxon 54800], Coffea arabica (arabica coffee, species) [taxon 13443], Pseudotsuga menziesii (Douglas-fir, species) [taxon 3357], Silybum marianum (blessed milkthistle, species) [taxon 92921], Paeonia ostii (species) [taxon 459177], Pinus koraiensis (channamu, species) [taxon 88728], Gossypium hirsutum (American cotton, species) [taxon 3635], Picea abies (Norway spruce, species) [taxon 3329], Quercus suber (cork oak, species) [taxon 58331], Arabidopsis thaliana (mouse-ear cress, species) [taxon 3702], H. brasiliensis [taxon 312095]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12913399/full.md

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12913399/full.md

## References

79 references — full list in the complete paper: https://tomesphere.com/paper/PMC12913399/full.md

---
Source: https://tomesphere.com/paper/PMC12913399