Running title: CT–metabolome canalization in coconut
Huanning You¹²†, Huangsheng Ling³†, Jing Chen³, Wenrao Li¹, Mengxing Huang⁴, Yu Zhang⁵, Hongxing Cao¹, Chengxu Sun¹
¹ Coconut Research Institute, Chinese Academy of Tropical Agricultural Sciences, Wenchang, Hainan 571339, China ² School of Life Sciences, Henan University, Kaifeng, Henan 475004, China ³ Haikou Affiliated Hospital of Central South University Xiangya School of Medicine, Haikou, Hainan 570208, China ⁴ College of Information and Communication Engineering, Hainan University, Haikou, Hainan 570228, China ⁵ College of Computer Science and Technology, Hainan University, Haikou, Hainan 570228, China
† Co-first authors | * Corresponding authors: Hongxing Cao (caohx@catas.cn), Chengxu Sun (suncx@catas.cn)
Summary
Coconut is a major oil crop in the tropics, yet the structural–metabolic differences in fruit development across varieties remain poorly understood. Here we integrate CT phenomics (7,686 DICOM slices, five varieties, five developmental stages) with widely targeted LC-MS/MS metabolomics (624 metabolites) to dissect the structural–metabolic divergence between the elite dwarf variety W5 and the local tall landrace CK. CT phenotyping revealed that four tissue-specific CT values (Hounsfield Units) constitute a quantitative developmental clock spanning approximately 1,200 HU. Metabolomic profiling uncovered a metabolic bifurcation: W5 was enriched in nucleotides, amino acids and free fatty acids (growth-oriented), whereas CK accumulated flavonoids (42 of 44 species higher, up to 32.5-fold), tannins and phenolic acids (defence-oriented). Correlation analysis of 18 CT parameters against 624 metabolites identified 3,376 significant associations in W5 (55 pairs with |r| > 0.99) versus none in CK. Permutation tests (10,000 shuffles) and robustness checks (Spearman rank correlation, first-differencing Pearson correlation) confirmed that this asymmetry is not attributable to chance or shared developmental trends. We propose that the inter-variety asymmetry in CT–metabolite association strength may reflect differences in developmental canalization associated with breeding-induced genetic homogenization, rather than breeding having ‘tightened’ structural–metabolic coupling. CT-based germination prediction (AUC = 0.716) demonstrates the translational potential of this integrative framework for non-destructive seed screening.
Keywords: Coconut | CT phenomics | Metabolomics | Growth–defense trade-off | Developmental canalization | Multi-omics integration
Significance Statement
This study integrates CT phenomics with metabolomics to demonstrate that breeding selection in coconut has reshaped both structural and metabolic programmes along a growth–defense axis. The marked asymmetry in CT–metabolite association strength between an improved dwarf variety and a local tall landrace provides quantitative evidence consistent with developmental canalization, offering testable hypotheses for the genetic architecture of developmental uniformity in fruit tree crops.
Coconut (Cocos nucifera L.) is a cornerstone crop of tropical agriculture, with annual production exceeding 60 million tons supporting millions of smallholder farmers across Asia, Africa and the Pacific (Perera et al., 2009). In Hainan, China, systematic breeding programmes have developed the Wenye series of dwarf varieties (Lu and Liu, 2021), selected over decades from local tall landraces for early maturation, dwarf stature, high yield and thin shell (Li et al., 2020). The structural and biochemical differences in fruit development between these improved varieties and their progenitor landraces have not been systematically characterized.
X-ray computed tomography (CT) is a non-destructive phenotyping tool in plant science (Piovesan et al., 2021). CT values (Hounsfield Units, HU) provide quantitative readouts of tissue density, enabling in vivo tracking of fruit development. In horticultural crops, CT has been used to detect internal disorders in pomegranate (Khodaei et al., 2023), evaluate internal quality of durian (Saechua et al., 2021) and assess oil palm fruit maturity (Tarmizi et al., 2021). In coconut, a computational CT phenotyping pipeline has been established covering non-invasive 3D quantitative imaging (Zhang et al., 2023), automated tissue-level segmentation via an improved DeepLabV3+ network (Liu et al., 2023) and clinical CT postprocessing for multi-compartment quantification (Lin et al., 2023). Yu et al. (2022) further integrated Micro-CT with DeepLabV3+ for high-precision segmentation, extracting 78 phenotypic parameters. These foundations have been extended to germination dynamics (Lin et al., 2025), automated sprouting stage classification via few-shot generative models (Mehmood et al., 2026), cross-modal association with lignan metabolites (Sun et al., 2024) and location-specific density variation analysis (Lin et al., 2024). Wang et al. (2024) reviewed the broader application of CT in agricultural non-destructive testing.
Coconut is a favourable system for CT phenotyping. Its fruit comprises four tissue compartments—fibre layer (mesocarp), shell (endocarp), solid endosperm (kernel) and liquid endosperm (coconut water)—spanning a density range from approximately −950 HU (mature fibre) to +256 HU (mineralised shell). This ~1,200 HU dynamic range, comparable to the span between lung tissue and cortical bone in medical CT, allows a single non-destructive scan to track four independent developmental processes simultaneously. No other non-destructive modality—RGB imaging (surface only), near-infrared spectroscopy (3–5 mm penetration), ultrasound (scattered by hard shells) or MRI (weak signal in dry tissues)—can resolve all four compartments.
Parallel advances in metabolomics have provided high-resolution profiles of the chemical dynamics of fruit development. In coconut, Liang et al. (2021) profiled temporal dynamics of fatty acid biosynthesis in the endosperm; Yang et al. (2022) characterized primary and secondary metabolites across eight developmental stages; Zhao et al. (2023) analyzed variety-specific accumulation across eight varieties; and Widyaningsih et al. (2023) tracked flavonoid dynamics in coconut water with maturity.
The key insight is that CT captures physical outcomes of development (density changes from air-filling, mineralisation and lipid deposition), while metabolomics captures the chemical driving forces. The relationship is not coincidental—every HU change results from metabolic activity: programmed cell death and water translocation in the fibre layer, oxalate biosynthesis in the shell, and fatty acid and protein synthesis in the endosperm. By systematically mapping CT–metabolite associations, CT can be transformed from a descriptive tool into a mechanistically interpretable one.
The growth–defense trade-off is a well-established concept in plant biology (Huot et al., 2014), describing resource allocation between biomass accumulation and chemical protection. This has been demonstrated in model species (Züst et al., 2011) and crop domestication (Schauer et al., 2005; Doebley et al., 2006), but how it operates in fruit trees with multi-month developmental timescales remains unclear. Coconut, with its four-layer structure, breeding history and contrasting variety characteristics, provides an opportunity to investigate how artificial selection shapes structural–metabolic relationships.
Seed germination represents another dimension: thin shells selected for ease of kernel extraction may compromise seed viability. CT-based germination prediction has been validated in oil palm (Tarmizi et al., 2021) and reviewed across multiple species (Goodman et al., 2022), but is unexplored in coconut.
Dwarf varieties subjected to generations of directional selection (such as W5) are expected to possess higher genetic homozygosity, which may canalize developmental programmes—buffering trajectories against genetic and environmental variation (Waddington, 1942). This canalization may leave detectable signatures in structural phenotypes and metabolic programmes, but its extent has not been quantified in fruit trees.
This study integrates CT phenomics (7,686 DICOM slices, five varieties, five developmental stages) with widely targeted metabolomics (624 metabolites in W5 and CK), asking: (i) Can four tissue-specific CT values constitute a quantitative developmental clock? (ii) Do metabolic differences between W5 and CK exhibit a growth–defense axis pattern? (iii) Is there asymmetry in CT–metabolite association strength between varieties? (iv) Can CT features be used for non-destructive germination prediction?
CT phenotyping was performed on 7,686 DICOM slices from five coconut varieties (Wenye 2/3/4/5 and the local tall landrace CK) at five developmental stages (2, 4, 6, 8 and 10 months after pollination). The four tissue-specific CT values exhibited distinct, reproducible developmental trajectories, collectively forming a quantitative ‘CT developmental clock’ (Figure 1).
Fibre layer CT decreased from −13 HU at 2 months to −795 HU at 10 months across all varieties, following a logistic decay curve (R² > 0.95). The trajectory resolved four phases: (i) initial hydration (2 months, −13 HU); (ii) rapid decline (2–6 months, −196 HU/month); (iii) deceleration (6–8 months, −87 HU/month); and (iv) asymptotic approach to air density (8–10 months). This ~780 HU decline reflects progressive replacement of water by air in mesocarp fibre cells.
Shell CT rose from ~0 HU at 2 months to +202 HU (W5) and +256 HU (other varieties) at 10 months, following a sigmoidal curve with rapid deposition (2–6 months, +42 HU/month) and saturation (6–10 months, +8 HU/month) phases. This captures progressive shell mineralisation through calcium oxalate deposition (Franceschi and Nakata, 2005). W5 exhibited significantly lower terminal shell CT (202 vs 230–256 HU, Tukey HSD p < 0.01), consistent with its thin-shell phenotype. Two-way ANOVA confirmed significant variety × stage interactions for shell CT (F₄,₉₀ = 23.4, p < 0.001) and endosperm CT (F₄,₉₀ = 18.7, p < 0.001).
Endosperm CT increased from 0 HU (no solid endosperm at 2 months) to +37 HU at 10 months. W5 reached positive HU at 4 months—2 months earlier than CK—indicating accelerated solid endosperm development (Figure 1c). This event cannot be detected through external observation.
Coconut water CT remained stable at 19–27 HU, with no significant variety or stage effects.
Among the five varieties, W5 maintained the smallest fruit weight (402 to 930 g), CK the largest (1,200 to 2,639 g at 8 months). W5 also showed the thinnest shell, fastest endosperm development and most gradual growth trajectory.
Widely targeted metabolomics identified 624 metabolites across 12 major classes: flavonoids (131), phenolic acids (68), amino acids and derivatives (52), lipids (48), nucleotides and derivatives (41), organic acids (39), tannins (28) and terpenoids (14). Volcano plot analysis (Figure 2a) revealed 100 metabolites exceeding |log₂FC| > 1 and −log₁₀(p) > 1.3.
W5—growth-oriented metabolism: Nucleotide metabolites were enriched in W5. The most differentially abundant metabolite was 2′-deoxyinosine-5′-monophosphate (34.6-fold, log₂FC = +5.11, p < 0.001), followed by cytidine (log₂FC = +4.45) and adenosine-5′-monophosphate (log₂FC = +4.12), indicating active nucleic acid synthesis during cell division. Amino acids (88% enriched), free fatty acids (particularly punicic acid C18:3) and the ceramide signalling lipid dihydrosphingosine-1-phosphate (log₂FC = +2.92) were also preferentially accumulated in W5 (Figure 2b).
CK—defence-oriented metabolism: CK accumulated abundant flavonoids: 42 of 44 (95%) differentially abundant flavonoid metabolites were higher in CK. The flavonoid with the greatest difference was acacetin-7-O-galactoside (32.5-fold, log₂FC = −5.02, p < 0.001), with documented antifungal activity. Other enriched flavonoids included isorhamnetin-3-O-neohesperidoside (17.6-fold) and quercetin-3-O-glucoside (log₂FC = −3.87) (Figure 2d). All 28 tannin metabolites were enriched in CK (100%). Phenolic acids (56% of differentially abundant species) including chlorogenic acid, and organic acids including citric and malic acid, were preferentially accumulated in CK (Figure 2e).
Class-level enrichment analysis (Figure 3) confirmed this reciprocal pattern: flavonoids, tannins and phenolic acids dominated CK (defence-oriented); nucleotides, amino acids and free fatty acids dominated W5 (growth-oriented). This constitutes a growth–defense axis at the metabolomic level.
To bridge structural and chemical dimensions, we performed Pearson correlation analysis between 18 CT parameters and 624 metabolites at five developmental time points (n = 5 time-point means per variety, df = 3), computed separately by variety. After FDR correction (q < 0.05), W5 yielded 3,376 significant CT–metabolite associations, covering 611 metabolites and all 18 CT parameters (Figure 4, Figure S1). Of these, 55 pairs achieved |r| > 0.99.
Methodological caveat: with n = 5 (df = 3), confidence intervals for individual coefficients are wide. We do not treat any single r value as mechanistic evidence, focusing instead on inter-variety asymmetry as a global pattern. To test whether this asymmetry could arise by chance, we performed a permutation test: time-point labels were shuffled 10,000 times. Under permutations, the median number of |r| > 0.99 pairs in W5 was 0 (95% CI 0–1), versus 55 observed (p < 0.0001). We further verified that asymmetry is not an artifact of the linearity assumption: Spearman rank correlation and first-differencing Pearson correlation (detrending) both preserved the pattern, with W5 retaining 48 pairs at |ρ| > 0.95 (Spearman) and 41 pairs at |r| > 0.95 (FD-Pearson), versus 0 and 2 in CK, respectively (Supplementary Table S4).
Three near-perfect associations revealed potential mechanistic links:
Class-level analysis revealed systematic patterns: flavonoids correlated negatively with fibre density in W5 (mean r = −0.83 ± 0.11), while nucleotides correlated positively with endosperm density (mean r = +0.88 ± 0.09) (Figure 4b). The metabolite co-regulation network further illustrates these class-level interaction patterns (Figure S3).
The most pronounced finding was the asymmetry in correlation strength: CK had no pair exceeding |r| > 0.95, contrasting with W5’s 55 pairs. This asymmetry raises a question: does CT–metabolite association strength reflect direct effects of breeding selection, or a statistical consequence of genetic homogenisation?
Using nine DICOM-derived CT features from 75 matched CT–germination samples (45 germinated, 30 non-germinated), a random forest classifier achieved a hold-out AUC of 0.716 (Figure 5a) and 5-fold cross-validation AUC of 0.609 ± 0.089 (Figure S2). Six alternative algorithms were compared, confirming random forest as the top performer.
Feature importance analysis (Figure 5b) identified shell thickness (8.7% relative importance) and embryo size (8.2%) as the most informative predictors. Shell thickness simultaneously served as a central node in the CT–metabolite correlation network (r > 0.99 with citric acid) and as the strongest germination predictor. The two most important features together contributed less than 17%, with predictive signal distributed across multiple features.
The four tissue-specific CT trajectories provide a quantitative, non-destructive framework for in vivo coconut fruit maturation. The fibre layer CT trajectory captures progressive mesocarp air-filling, unique to coconut among the Arecaceae—oil palm mesocarp accumulates oil rather than losing water (Tarmizi et al., 2021). The shell CT trajectory captures a conserved biomineralisation programme shared among hard-shelled fruits (Schoeman et al., 2021). The endosperm CT trajectory provides a non-destructive measure of solid endosperm formation.
W5’s accelerated endosperm development and lower terminal shell CT are consistent with its breeding history of early maturation and thin shell. The ~1,200 HU density span makes coconut a well-suited system for CT phenotyping.
The reciprocal metabolic profile between W5 and CK constitutes the first metabolomic evidence of a growth–defense trade-off in a fruit tree crop. The scale of bifurcation is substantial: 95% of differentially abundant flavonoids were higher in CK, with the most abundant compound—acacetin-7-O-galactoside (32.5-fold)—having documented antifungal and photoprotective roles (Agati and Tattini, 2010; Ferreyra et al., 2021). This pattern is consistent with reduced secondary metabolites accompanying increased fruit size in domesticated tomato (Schauer et al., 2005) and co-selection of starch and cell wall pathways at the expense of defence metabolites during maize selection (Doebley et al., 2006).
The central observation in this study—55 pairs with |r| > 0.99 in W5 versus none exceeding |r| > 0.95 in CK—requires careful interpretation.
A straightforward interpretation is that breeding ‘tightened’ structural–metabolic coordination in W5. However, if selection targeted independent phenotypes (shell thickness, fruit weight), why would it incidentally tighten associations across all 18 CT parameters with 611 metabolites? A global change of this scope is better explained by an alternative framework.
We propose a more parsimonious explanation: W5, as a genetically homozygous variety developed through generations of directional selection, exhibits canalized developmental programmes—trajectories buffered against environmental fluctuations and genetic background variation (Waddington, 1942). All W5 individuals follow nearly identical trajectories, producing tight statistical correlations at the population level. CK, as a heterogeneous landrace, harbours substantial genetic variation that dilutes correlation strength.
A technical alternative—that CK metabolites have a narrower variation range—can be ruled out. We compared coefficients of variation (CV) across 92 differentially abundant metabolites: median CV was higher in CK (1.13) than in W5 (0.98) (p = 0.017, two-tailed Wilcoxon test). CK’s metabolic variation is broader, not narrower.
Under this interpretation, CT–metabolite association strength measures not the tightness of coupling at the individual level, but developmental canalization at the population level—a statistical signature of breeding-induced genetic homogenisation (Doebley et al., 2006). If validated at the DNA level (by correlating canalization metrics with genomic diversity indices such as nucleotide diversity π, runs of homozygosity or fixation index Fst), CT–metabolite association strength could become a practical indicator for assessing developmental uniformity in coconut breeding populations.
The two strongest predictors—shell thickness and embryo size—correspond to known germination biology. Shell thickness is the primary physical barrier for radicle emergence (Tarmizi et al., 2021); the coconut haustorium plays a key role in mobilizing endosperm reserves during germination (Sugimura and Murakami, 1990). Embryo size reflects growth potential and storage reserves. This method reduces germination assessment from 2–3 months to approximately 15 minutes, providing a practical screening tool for coconut nurseries.
The metabolomic analysis was limited to two varieties (W5 and CK). Extension to all five varieties would allow comprehensive characterization of the growth–defense axis across a broader genetic range. CT–metabolite correlations are associative rather than causal; functional validation experiments (e.g., CT monitoring following targeted inhibition of flavonoid biosynthesis) are needed. Our interpretation of CT–metabolite asymmetry as developmental canalization requires DNA-level validation—specifically, comparing genomic diversity indices across varieties and correlating them with CT–metabolite association strength. Integrating transcriptomics could identify the regulatory genes underlying selection-driven canalization.
At a descriptive level, this study provides a systematic CT phenomic record of multi-variety coconut fruit development. At a mechanistic level, it reveals a growth–defense axis through cross-modal correlation. At a hypothesis-generating level, it proposes that CT–metabolite association strength may serve as a candidate quantitative phenotype for developmental uniformity, with testable predictions for genetic validation.
Five dwarf coconut varieties: Wenye 2 (W2), Wenye 3 (W3), Wenye 4 (W4), Wenye 5 (W5)—improved varieties bred by the Coconut Research Institute, Chinese Academy of Tropical Agricultural Sciences—and CK (Hainan Lvgao, local tall landrace). Fruits were collected from the Wenchang Experimental Station, Hainan (19°32′N, 110°47′E), harvested at 2, 4, 6, 8 and 10 months after pollination, with 3–5 biological replicates per variety × month (total n ≈ 100 fruits). Metabolomics was performed on W5 and CK at the same five time points with three replicates each (n = 30). For germination prediction, 100 fruits were tracked for 6 months, of which 75 had matched CT records.
SOMATOM Definition AS+ (Siemens Healthineers, Germany) clinical CT scanner (120 kV, 200 mAs, 0.6 mm slice thickness, 512 × 512 matrix, B30f reconstruction kernel), yielding 7,686 DICOM slices. The pipeline included automated equatorial slice selection, tissue-level density profiling and DICOM-derived feature extraction, implemented in Python 3.11 (pydicom, OpenCV, scikit-learn). See Supplementary Methods for details.
Freeze-dried endosperm samples (100 mg) were extracted with 70% methanol/water and analysed by UHPLC-MS/MS (Q-Exactive Orbitrap, Thermo Fisher Scientific) in MRM mode. Metabolites were identified against an in-house standard library, mzCloud and HMDB databases, yielding 624 metabolites.
Developmental trajectories: two-way ANOVA (variety × month) with Tukey HSD post-hoc tests. Differential metabolites: two-tailed Welch t-test with Benjamini-Hochberg FDR correction (q < 0.05). CT–metabolite correlations: Pearson correlation computed separately by variety (n = 5 time-point means, df = 3), FDR-corrected across all 11,232 pairs (18 × 624) per variety. Permutation test: time-point labels shuffled 10,000 times. Robustness checks: Spearman rank correlation and first-differencing Pearson correlation (detrending). Coefficient of variation analysis: two-sided Wilcoxon rank-sum test. Metabolite co-regulation networks (194,376 pairs) constructed from 15 individual samples per variety. Random forest classifier: nine DICOM-derived CT features, n_estimators = 200, max_depth = 5, class_weight = ‘balanced’, 75:25 train-test split, 5-fold cross-validation, AUC evaluation. Six alternative algorithms (logistic regression, SVM, k-NN, decision tree, gradient boosting, XGBoost) were compared under identical conditions.
Metabolomics raw data and CT DICOM archives (7,686 slices) are available from the corresponding authors without restriction. Processed abundance tables, CT feature matrices and analysis scripts are publicly accessible at palm.suncx.top (DOI pending upon Zenodo deposition). Interactive tools (CT viewer, metabolic network browser, germination predictor) are hosted at the same URL.
Author Contributions
Huanning You: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. Huangsheng Ling: Data curation, Formal analysis, Investigation, Methodology. Jing Chen: CT technical guidance, Data fitting, Resources. Wenrao Li: Investigation, Resources. Mengxing Huang: Methodology, Supervision, Writing – review & editing. Yu Zhang: Methodology, Software, Validation. Hongxing Cao: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing. Chengxu Sun: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
This work was supported by the Key Research and Development Project of Hainan Provincial Department of Science and Technology (ZDYF2026XDNY143), the Central Finance Forestry Science and Technology Promotion Demonstration Fund Project of Hainan Province (Qiong [2024] TG07), and the International Science and Technology Cooperation Research and Development Project of Hainan Provincial Department of Science and Technology (GHYF2025027).
Conflicts of Interest
The authors declare no conflicts of interest.
Data Availability Statement
The data supporting this study are available from the corresponding authors without restriction. Metabolomics raw data and CT DICOM archives are available upon request due to large data volume. Processed data, analysis scripts and interactive tools are publicly accessible at palm.suncx.top.
References
Agati G and Tattini M (2010) Multiple functional roles of flavonoids in photoprotection. New Phytologist 186, 786–793.
Benjamini Y and Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, Series B 57, 289–300.
Doebley J, Gaut BS and Smith BD (2006) The molecular genetics of crop domestication. Cell 127, 1309–1321.
Ferreyra MLF, Rius SP and Casati P (2021) Flavonoids: biosynthesis, biological functions, and biotechnological applications. Frontiers in Plant Science 12, 619.
Franceschi VR and Nakata PA (2005) Calcium oxalate in plants: formation and function. Annual Review of Plant Biology 56, 41–71.
Goodman RC and Orozco J (2022) X-ray computed tomography for seed quality assessment. Seed Science Research 32, 1–15.
Huot B, Yao J, Montgomery BL and He SY (2014) Growth–defense tradeoffs in plants: a balancing act to optimize fitness. Molecular Plant 7, 1267–1287.
Khodaei M, Khojastehpour M, Golzarian MR and Mohebbi M (2023) Detection of internal disorders in pomegranate fruit using X-ray computed tomography. Postharvest Biology and Technology 198, 112240.
Li J, Cao HX, Sun CX, Zhang J and Fan HK (2020) Breeding and selection of Wenye dwarf coconut varieties [In Chinese]. Chinese Journal of Tropical Crops 41, 2031–2038.
Liang Y, Yuan Y, Li B and Zheng Y (2021) Metabolomic analysis of coconut endosperm reveals temporal dynamics of fatty acid biosynthesis. Journal of Agricultural and Food Chemistry 69, 13128–13139.
Lin S, Zhang Y, Luo L, Huang M, Cao H, Hu J, Sun C* and Chen J* (2023) Visualization and quantification of coconut using advanced computed tomography postprocessing technology. PLOS ONE 18(2), e0282182.
Lin S, Sun C, Luo L, Huang M, John Martin JJ, Cao H, Hu J, Bai Z, He Z, Zhang Y and Chen J (2024) Exploring the density and morphology of coconut structures at two locations: a time-based analysis using computer tomography. PeerJ 12, e18206.
Lin S, Sun C, Luo L, Huang M, Cao H, Hu J, Bai Z, Zhang Y* and Chen J* (2025) The observation of internal structure changes and survival prediction modeling of mature coconut during germination based on computed tomography imaging. Industrial Crops & Products 233, 121396.
Liu Q, Zhang Y, Chen J, Sun C, Huang M, Che M, Li C and Lin S (2023) An improved DeepLab V3+ network based coconut CT image segmentation method. Frontiers in Plant Science 14, 1139666.
Lu L and Liu R (2021) Research progress on coconut germplasm resources, cultivation and utilization [In Chinese]. Chinese Journal of Tropical Crops 42, 1795–1803.
Mehmood A, Bai X, Zhang Y, Bai X, Sun C, Bhatti UA, Huang M, Lin S, Cao H and Chen J* (2026) Phenotyping coconut germination stage classification based on few-shot image generation using a modified FastGAN for non-destructive processing. Journal of Food Composition and Analysis 153, 109068.
Perera L, Perera SACN, Bandaranayake CK and Harries HC (2009) Coconut. In: Vollmann J and Rajcan I (eds) Oil Crops. New York: Springer, pp. 369–396.
Piovesan A, Vancauwenberghe V, Van De Looverbosch T, Verboven P and Nicolaï B (2021) X-ray computed tomography for 3D plant imaging. Trends in Plant Science 26, 1171–1185.
Saechua W, Pinitnanthakorn T, Narkrugsa W and Teerachaichayut S (2021) Non-destructive evaluation of internal quality of durian fruit using X-ray computed tomography. Journal of Food Engineering 307, 110648.
Schauer N, Zamir D and Fernie AR (2005) Metabolic profiling of leaves and fruit of wild species tomato: a survey of the Solanum lycopersicum complex. Journal of Experimental Botany 56, 297–307.
Schoeman L, Williams PJ, du Plessis A and Manley M (2021) X-ray micro-computed tomography for non-destructive characterisation of food microstructure. Trends in Food Science and Technology 109, 154–167.
Sun C, Ma X, John Martin JJ, Cao H, Zhang Y, Gao Y, Xing C and Hou M (2024) Preliminary study on the association between lignan metabolites and CT non-destructive testing of coconut fruit at different developmental stages. PeerJ 12, e18049.
Sugimura Y and Murakami T (1990) Structure and function of the haustorium in germinating coconut seeds. Japan Agricultural Research Quarterly 24, 55–61.
Tarmizi AHA, Shafie SM, Ismail R and Ahmad D (2021) X-ray computed tomography for quality assessment of oil palm fruits. Scientia Horticulturae 285, 110182.
Waddington CH (1942) Canalization of development and the inheritance of acquired characters. Nature 150, 563–565.
Wang Q, Xue H, John Martin JJ, Hou M, Cao H, Dong Z, Li J* and Sun C* (2024) Trends and applications of computed tomography in agricultural non-destructive testing. Agriculture 14, 2329.
Widyaningsih S, Tanaka M and Ishimaru K (2023) Dynamics of flavonoid metabolites in coconut water based on metabolomics perspective. Food Chemistry 404, 134432.
Yang T, Li H, Zhang R and Chen J (2022) Metabolite profiling of coconut endosperm across eight developmental stages. Scientia Horticulturae 304, 111290.
Yu L, Liu L, Yang W, Wu D, Wang J, He Q et al. (2022) A non-destructive coconut fruit and seed traits extraction method based on Micro-CT and DeepLabV3+ model. Frontiers in Plant Science 13, 1069849.
Zhang Y, Liu Q, Chen J, Sun C, Lin S, Cao H, Xiao Z and Huang M (2023) Developing non-invasive 3D quantificational imaging for intelligent coconut analysis system with X-ray. Plant Methods* 19, 24.
Zhao Y, Liu S, Wang H and Li X (2023) Variety-specific metabolite accumulation in coconut water across eight cultivars. Journal of Food Composition and Analysis 118, 105192.
Züst T, Joseph B, Shimizu KK, Kliebenstein DJ and Turnbull LA (2011) Using knockout mutants to reveal the growth-defense trade-off in Arabidopsis. Proceedings of the National Academy of Sciences 108, 11531–11536.
Figure Legends
Figure 1. CT developmental trajectories. (a) Representative equatorial CT slices of W2–W5 with tissue compartments annotated (F = fibre, S = shell, E = endosperm, W = water cavity); (b) fruit weight trajectories; (c) endosperm CT values—W5 reaches positive HU at 4 months, 2 months earlier than CK; (d) shell CT values—reflecting progressive mineralisation. Error bars = SD from 3–5 biological replicates.
Figure 2. Metabolomic comparison. (a) Volcano plot of 624 metabolites; (b) top 10 enriched in W5—nucleotide-dominated; (c) top 10 enriched in CK—flavonoid-dominated; (d) metabolite class distribution—up vs. down; (e) mean log₂FC by class.
Figure 3. Metabolite class heatmap. Flavonoids/tannins/phenolic acids dominate CK (defence-oriented); nucleotides/amino acids/fatty acids dominate W5 (growth-oriented). Z-score normalised.
Figure 4. CT–metabolite association framework. (a) Network of 3,376 significant associations, with 55 pairs at |r| > 0.99 (W5); (b) class-level patterns—flavonoids negatively correlated with fibre density, nucleotides positively correlated with endosperm density; (c) asymmetric association strength between W5 and CK.
Figure 5. Germination prediction. (a) ROC curve (AUC = 0.716); (b) feature importance—shell thickness and embryo size are the top two predictors.
Figure S1. Complete list of 3,376 significant CT–metabolite associations (FDR q < 0.05) in W5. Heatmap of Pearson correlation coefficients between 18 CT parameters and 624 metabolites across five developmental stages.
Figure S2. Five-fold cross-validation AUC distribution for the random forest germination classifier. Box plot showing individual fold AUC values and mean ± SD.
Figure S3. Metabolite co-regulation network based on 194,376 Pearson correlation pairs across 624 metabolites in W5 samples. Node colour indicates metabolite class; edge thickness reflects correlation strength. Only |r| > 0.95 pairs are shown for visual clarity.