Cyperus esculentus L., commonly known as yellow nutsedge, tiger nut, or chufa, is a perennial sedge in the family Cyperaceae[1]. This monocotyledon occurs both as a wild weedy species and as a deliberately cultivated food crop (Figure 1(A)). The species produces small, calorie-dense underground tubers as its primary storage organ (Figure 1(B)). These tubers have sustained human populations since antiquity and continue to support commercial food industries in Spain, China, Nigeria and Ghana [2]. Despite this long history of use and its documented nutritional value, tiger nut remains classified among the neglected and underutilized crop species of Africa. It has received comparatively little formal agronomic research investment relative to major staple crops, and its cultivation and value-chain development in East Africa in particular remain considerably less advanced than in West Africa or Spain.

Tuber size represents one of the most agronomically consequential traits in Cyperusesculentus L., since tuber dimensions relate directly to marketable yield, mechanical harvesting efficiency, and propagation performance [2]. Given ongoing efforts to combat malnutrition through nutritious locally sourced foods, understanding the tuber traits of this plant is valuable for potential crop improvement, domestication, and new adoptions [3].

Since yellow nutsedge propagates through tubers formed at the apical ends of rhizomes [4], rather than through seed, tuber-specific traits may be inherited or expressed somewhat independently of the broader vegetative morphology. Documenting tuber size and shape variation directly therefore serves as a necessary complement to, rather than a substitute for, whole-plant morphological characterization, and provides breeders and germplasm curators with trait-specific information relevant to propagation and yield-oriented selection decisions. This study therefore aimed to characterize and document the tuber trait diversity among Cyperus esculentus L. accessions collected from East Africa, Egypt and Spain.

Figure 1.Cyperus esculentus L. mature plant (A) and (B) C.esculentus L. tubers.

Table 2 shows accessions’ mean ± SD Tuber Dimensions (mm) for the four tuber dimensions in 27 Cyperus esculentus L. Accessions. Each mean is based on three pots, with one tuber measured per pot.

Table 2.Tuber dimensions (mm) for the selected Cyperus esculentus L. accessions.

The mean tuber length across the collection was 8.12 mm ± 2.40 mm, ranging more than three-fold from 5.00 mm in accession MAR to 15.00 mm in accession EGY. Width at the base averaged 5.10 mm ± 1.53 mm and ranged from 3.33 mm (MAR, OGEMBO) to 8.67 mm (BIO). Width at the mid-point, the broadest dimension recorded on average (mean = 6.26 mm), ranged from 4.33 mm (NKU-M) to a striking 11.67 mm in BIO, the single largest individual trait value recorded anywhere in the dataset. Width at the tip averaged 4.46 mm ± 1.32 mm and ranged from 3.00 mm (IDO, SOK) to 8.83 mm (BLACK). The four cultivated accessions had larger tubers than the wild accessions, with a mean length of 12.0 mm (range 8.83 to 15.00 mm) against 7.45 mm (range 5.00 to 10.17 mm), and a mean mid-point width of 9.42 mm against 5.71 mm.

Coefficients of variation (CV) remained consistently high and closely comparable across all four traits, ranging narrowly from 27.6% (width mid) to 30.0% (width base). This narrow spread in CV values across traits of otherwise quite different magnitude indicates that the four tuber dimensions scale together fairly proportionally across the collection. Accessions with unusually long tubers tend also to have unusually wide tubers, rather than length and width varying independently at markedly different rates.

The observed dimensions fall within the range reported for the species, with the wild accessions toward its lower end. The cultivated Spanish and Egyptian accessions (8.8 to 15.0 mm long and 7.3 to 11.7 mm wide at the mid-point) resemble the morphotypes described from Burkina Faso (9.8 to 13.1 mm long and 9.0 to 11.9 mm wide) [5] and the site means reported from Russia (10.9 to 12.7 mm long) [6]. Most wild East African accessions produced tubers 5 to 9 mm long and 4 to 6 mm wide, in line with descriptions of smaller tubers in weedy than in cultivated forms of the species [7]. The CVs of 27.6 to 30.0% exceed the moderate CVs (15.6 to 19.1%) reported for traits including tuber dimensions among 42 accessions grown in China [8], which likely reflects the inclusion of both wild and cultivated material in the present collection. Length-to-mid-width ratios ranged from 0.89 (BIO) to 2.05 (EGY), spanning shapes comparable to the round, oval, and long classes described by Yang et al. [8].

EGY, BLACK, BIO, KHA, and KARI-E were the five longest-tubered accessions, while MAR, KAM, EMBU G, SOK, and OGEMBO formed the five shortest-tubered accessions. One-way ANOVA was conducted separately for each of the four tuber traits to test whether accession identity explained a statistically significant proportion of the observed variation. Accession identity explained a highly statistically significant proportion of variation in all four traits (p < 0.001 in every case), with F-statistics ranging from 9.18 (width at the mid) to 17.00 (length) (Table 3). The comparatively larger F-statistic obtained for length indicates that this trait discriminated most sharply among accessions relative to the within-accession replicate variability, consistent with the wide three-fold range and high CV. The smaller, though still highly significant, F-statistic for width at the mid reflects the somewhat greater within-accession replicate variability observed for this trait, particularly evident in accessions such as NGA (SD = 2.65 mm) and KHA (SD = 2.02 mm), where individual tuber measurements varied considerably around the genotype mean. Because accessions were propagated clonally and grown in one environment, repeatability here combines genetic effects, carry-over effects from the mother tuber, and accession-specific responses to the greenhouse. It therefore represents an upper limit for broad-sense heritability rather than an estimate of it [9].

Table 3. One-way ANOVA results for tuber traits across 27 accessions.

Note. Significance assessed at α= 0.05 (df = 26, 54).

Pearson correlation coefficients were calculated among the four tuber traits using genotype means and presented in Table 4.

Table 4. Pearson correlation matrix for tuber traits.

All pairwise correlations were positive and, with the exception of the length-width mid pairing, uniformly strong (r > 0.82), consistent with the expectation that linear dimensions of the same physical organ tend to increase together as overall tuber size increases. The strongest single correlation was recorded between width at the mid and width at the tip (r = 0.912), indicating that these two width measurements track one another very closely across the collection and may carry substantially overlapping information for classification purposes. Width at the base likewise correlated strongly with both width at the tip (r = 0.851) and length (r = 0.841), positioning it as a dimension broadly representative of overall tuber girth and size.

The comparatively weaker correlation between length and width at the mid (r = 0.682), while still statistically meaningful, represents the weakest pairwise relationship in the matrix and is directly responsible for the emergence of a distinct shape axis in the subsequent PCA. In practical terms, this weaker correlation means that knowing an accession’s tuber length alone provides a less reliable prediction of its mid-tuber girth than knowing its base width or tip width would predict the other width measurements. This is precisely the pattern underlying the length-to-mid-width shape ratio where accessions such as EGY (long, comparatively narrow) and BIO (comparatively short, very wide at the mid-point) diverge most sharply from a simple proportional size relationship.

PCA was performed on standardized (z-score) genotype means for the four tuber traits to summarize the correlation structure into a smaller number of interpretable axes. PCA of the standardized accession means showed the first two principal components together explained 95.6% of total variation among accessions (Table 5), an unusually high cumulative proportion that follows directly from the strong positive intercorrelation among all four traits and the small number of variables analysed. With only four highly correlated linear measurements available, most of the meaningful variation in the dataset is necessarily concentrated into very few dimensions.

Table 5.Percentage variance explained and trait loadings for PC 1 and PC 2.

PC1, accounting for 86.3% of total variation on its own, carried positive loadings of broadly similar magnitude from all four traits (0.467 to 0.519), confirming its interpretation as a general tuber size axis rather than a contrast between specific dimensions. Accessions positioned toward the positive end of PC1, most prominently BIO, BLACK, and EGY, followed by KHA, KARI-E, DOM, BUB, and MICRO, are distinguished by consistently larger measurements across all four traits simultaneously. Accessions positioned toward the negative end of PC1, including MAR, KAM, SOK, EMBU G, OGEMBO, NKU-M, and IDO, are correspondingly smaller across all four traits. This size gradient along PC1 aligns closely with the raw ranking of accessions by mean tuber length reported earlier, confirming that PC1 essentially recovers and formalizes the size ordering already evident from the descriptive statistics.

PC2, accounting for a further 9.3% of variation, presented a markedly different loading pattern: length loaded strongly and positively (0.793), while width at the mid and width at the tip loaded moderately and negatively (−0.465 and −0.378, respectively), with width at the base contributing only a small positive loading (0.108). This pattern identifies PC2 as a tuber elongation or shape axis, orthogonal to the general size axis captured by PC1. Accessions scoring highly on PC2, namely EGY, SUK and CHU 18, combine substantial length with comparatively narrow mid- and tip-width measurements, producing visibly elongated tuber shapes. BIO, despite its high PC1 score reflecting large overall size, scored the lowest of any accession on PC2, reflecting its unusually compact, rounded tuber shape already noted in the length-to-mid-width ratio analysis. This separation between size (PC1) and shape (PC2) is visualized directly in the PCA biplot (Figure 2Figure 2), where BIO and EGY occupy near-opposite vertical positions despite both scoring toward the positive end of the horizontal PC1 axis. The orthogonality of PC1 and PC2 does not show that tuber size and shape are biologically independent. This is because principal components are uncorrelated by construction, regardless of the underlying biology [10].

Figure 2. PCA biplot of 27 Cyperus esculentus L. accessions based on four standardized tuber morphometric traits.

The remaining two components, PC3 and PC4, together accounted for only 4.4% of total variation and are not interpreted further, since the near-complete summarization achieved by PC1 and PC2 indicates that a two-dimensional size-and-shape representation captures all systematic morphological variation present in this four-trait tuber dataset essentially.

Ward’s minimum-variance hierarchical clustering was applied to Euclidean distances calculated from the standardized genotype means. Bootstrap support values were generated by resampling the four available traits with replacement across 1,000 iterations and recomputing the dendrogram topology on each resampled dataset. The resulting dendrogram (Figure 3Figure 3) divided the collection into two principal branches at a Ward distance of approximately 11, though bootstrap support for this primary split was weak (28%), reflecting the limited discriminating power available from only four correlated size-related traits. The smaller of the two principal branches, shown in green, grouped eight accessions distinguished by consistently larger tuber dimensions: BIO, BLACK, EGY, KHA, MICRO, BUB, DOM, and KARI-E. This grouping corresponds closely, though not perfectly, to the accessions scoring most positively along PC1 in the preceding PCA, reinforcing that both analytical approaches are detecting the same underlying large-tuber size gradient rather than providing independent structural evidence.

The larger principal branch, shown in orange, grouped the remaining 19 small-to-medium tubered accessions and subdivided further into several smaller sub-clusters at generally low-to-moderate bootstrap support, including a grouping of MAR, EMBU G, and OGEMBO; a grouping of NKU-M, IDO, KAM, and SOK; a grouping built around CHU 18 and SUK; and a grouping built around NGA, MAZ, and BUWA. Several individual accession pairs within these sub-clusters returned comparatively strong bootstrap support despite the limited trait set available, most notably BIO-BLACK (86%), DOM-KARI-E (84%), MAR-EMBU G (74%), and CHU 18-SUK (52%). These higher-support pairings identify the accession relationships that are most consistently recovered regardless of which subset of the four traits is emphasized in a given bootstrap resample, and thus represent the most reliable tuber-based groupings identified in this analysis.

The close correspondence between the PCA and HCA results in this analysis is an expected consequence of the small number of highly intercorrelated traits examined: with only four variables, most of which correlate at r > 0.68, both ordination and clustering methods necessarily converge on describing the same dominant size gradient, offering limited scope for genuinely independent structural insight.

Figure 3.Hierarchical cluster analysis of 27 C. esculentus L. accessions based on four tuber morphometric traits.

This study concludes that the mean tuber length of the 27 Cyperus esculentus L. accessions collected from wild and cultivated sources across East Africa, Egypt, and Spain ranged more than three-fold, from 5.00 mm (MAR) to 15.00 mm (EGY). All four traits showed consistently high coefficients of variation (27.6% - 30.0%). One-way ANOVA confirmed that accession identity explained a highly significant proportion of variation in every trait (p < 0.001), and Pearson correlation showed all four dimensions positively intercorrelated, most strongly between width at the mid and tip (r = 0.912), most weakly between length and width at the mid (r = 0.682). Principal component analysis (PCA) reduced this correlation structure to two axes explaining 95.6% of total variation: a dominant size axis (PC1, 86.3%) and a secondary elongation axis (PC2, 9.3%). The orthogonality of these axes is a property of PCA and does not demonstrate independent control of tuber size and shape. The four cultivated accessions from Spain and Egypt produced the largest tubers, while KHA, KARI-E, DOM, and BUB were the largest-tubered wild East African accessions. Hierarchical cluster analysis (HCA) grouped eight large-tubered accessions (BIO, BLACK, EGY, KHA, MICRO, BUB, DOM, KARI-E) apart from the remaining nineteen, though with weak bootstrap support (28%). However, because the data come from one environment and from a mixture of wild and imported cultivated material, they establish phenotypic differences under these conditions but not the heritable basis of the diversity or the crop potential of East African accessions.

It is recommended that future work on this Cyperus esculentus germplasm collection should: 1) evaluate the accessions in replicated multi-location field trials over more than one season, with several tubers measured per plant, to estimate genotype-by-environment interaction and broad-sense heritability, 2) genotype the accessions to test whether tuber-size groupings reflect genetic relationships, 3) prioritize the largest-tubered wild East African accessions (KHA, KARI-E, DOM, and BUB), alongside the cultivated references BIO, BLACK, and EGY, for nutritional analysis to test whether tuber size predicts oil, starch, protein and mineral content, 4) compare elongated and rounded morphotypes for oil extraction and processing quality, 5) test congruence among tuber, vegetative and nutritional datasets with Mantel tests and combined ordination, and 6) extend collection and characterization to additional East African collection sites to identify populations best suited to advancing food and nutritional security in the region.

Accession-level summaries (mean ± SD of three pots per accession) are reported in Table 2. The individual pot-level measurements (81 records: 27 accessions × 3 pots, one tuber per pot, four traits) can be provided by the corresponding author if needed. The analyzed data generated are included in this published article.

Conceptualization: Rosbella Jemurgor, Elizabeth W. Njenga, Emily J. Too, Beatrice A. Were, and Anders S. Carlsson; Methodology: Rosbella Jemurgor, Elizabeth W. Njenga, Emily J. Too, Beatrice A. Were, and Anders S. Carlsson; Formal analysis: Rosbella Jemurgor, Elizabeth W. Njenga, Emily J. Too, Beatrice A. Were, and Anders S. Carlsson; Investigation: Rosbella Jemurgor; Resources: Anders S. Carlsson; Data curation: Rosbella Jemurgor; Writing of original draft preparation: Rosbella Jemurgor; Review and editing: Elizabeth W. Njenga and Emily J. Too; Supervision: Elizabeth W. Njenga, Emily J. Too, Beatrice A. Were, and Anders S. Carlsson. All authors have read and contributed to the published version of the manuscript.