Fonio Mali and Guinea accessions were phenotyped for 18 metabolomic, physiological and morphometric traits once the soil volumetric water content (VWC) achieved the target (at 21 DAE), during the drought period (21–41 DAE), and/or after recovery (61 DAE):
Metabolomics
Analysis of 30 metabolites in fonio leaf tissue known to contribute to drought tolerance in other plants showed that fonio leaves accumulate betaines, namely, alanine betaine, glycine betaine and proline betaine under drought stress (Fig. 2A-B, Table S1A-B). There were remarkable 8- and 60-fold increases (P ≤ 0.05) in alanine betaine in the drought-exposed (6%-VWC) plants of the Mali accession in trials 1 and 2, respectively, whereas it was 9 and 6-fold in the case of the Guinea accession (Fig. 2A-B). Similarly, glycine betaine increased 7–12 times across accessions, which was significantly greater (P ≤ 0.05) than in optimally watered (26%-VWC) plants (Fig. 2A-B, Table S1A-B). Proline betaine concentrations also increased 2-fold (P ≤ 0.10) Proline betaine concentrations also increased 2-fold(P ≤ 0.10)under drought, except in the Guinea accession in Trial 1. No significant difference were observed in the concentrations of the above three metabolites when comparing the two accessions under drought (Table S1A-B). A 2-fold increase in proline was noted in trial 1 for both genotypes, but this was not consistent in Trial 2 (Fig. 2A).
Short term above-ground effects
Physiological traits
SPAD at 10 day intervals was measured to determine the impact of drought on fonio leaf chlorophyll content. The Mali accession showed a general trend of superior SPAD values over Guinea under water limitation (Fig. 3A, Table 1A). Compared to optimal water, drought caused a general trend of decreased SPAD chlorophyll in both accessions, but it was only significant for the Mali accession in Trial 2. As Mali plants aged, they showed a trend of increasing SPAD values, regardless of the water treatment, whereas the Guinea accession was more static (Table 1A).
With respect to rates of photosynthesis, no significant differences were observed between plants grown at 26% vs 16%-VWC, but significant declines were observed for 6%-VWC plants starting after 10 days of drought exposure (31 DAE), with a progressive decline until the end date of the treatment. In the Mali accession, 6%-VWC plants showed a 47–51% decline in photosynthesis compared to 26%-VWC plants, whereas this decline was 28–34% for the Guinea accession at 41 DAE (both significant at p ≤ 0.05) (Fig. 3B, Table 1B). When comparing the accessions, no consistent significant differences in photosynthesis were observed, though the Mali accession had a significantly lower rate at the end of the drought treatment (41 DAE).
Next, we asked whether the fonio accessions could acclimate to drought by reducing transpiration. Consistent with photosynthesis, whereas optimally watered plants showed increases in transpiration as plants aged (21–41 DAE), notable declines in transpiration were observed in both accessions (52–57% for Mali; 37–45% for Guinea) after 20 days of drought exposure (Fig. 3C, Table 1C). However, within each accession, significant declines in photosynthesis between 6% and 26%-VWC were inconsistently observed starting after 10 days of drought exposure (31 DAE), becoming consistent only at 41 DAE.
As a potential contributor to restricted transpiration, stomatal conductance was assessed at 5 day intervals during the drought treatment. Similar to transpiration, stomatal conductance in optimally watered plants showed increases in transpiration as plants developed (21–41 DAE), whereas stomatal conductance declined in both accessions exposed to drought as the drought interval progressed (Fig. 3D, Table 1D). By 41 DAE, stomatal conductance declined 63–64% for the Mali accession and 53–54% for Guinea when comparing 26–6%-VWC plants. No significant differences were observed in stomatal conductance when comparing the two accessions.
Shoot morphology
In addition to altered stomatal conductance, we hypothesized that the reduced transpiration observed in white fonio under drought was partly due to shoot acclimation responses. Drought caused significant declines in shoot biomass in both accessions (Fig. 4A Table S3). The Mali accession at 6%-VWC had a 48% and 45% reduction in shoot weight per
Table 1
Experiment A: Effect of moisture treatments on physiological traits of white fonio. A. Leaf Chlorophyll (SPAD)
| 21 DAE | 31DAE | 41 DAE | 61 DAE (recovery) |
Mali | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 39.61a | 37.78ab | 40.65a | 41.57b | 42.36a | 41.71b | 44.46a | 48.88b |
16% VWC | 38.71a | 37.98ab | 42.20a | 46.96a | 44.96a | 47.54a | 45.83a | 54.06a |
26% VWC | 38.97a | 39.36a | 41.63a | 47.81a | 45.67a | 46.12a | 46.38a | 51.89ab |
Guinea | | | | | | | | |
6% VWC | 36.63a | 36.28ab | 34.13b | 36.85c | 35.15b | 37.16c | 38.44b | 38.60c |
16% VWC | 37.35a | 35.84b | 36.81b | 38.56c | 37.34b | 38.18c | 39.28b | 39.59c |
26% VWC | 37.22a | 37.69ab | 36.00b | 39.19bc | 38.34b | 39.01c | 39.03b | 39.71c |
B. Photosynthesis (µmolCO2m− 2s− 1) |
| 21 DAE | 26 DAE | 31DAE | 36 DAE | 41 DAE |
Mali | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 15.16ab | 12.67a | 15.94ab | 11.75a | 13.87c | 11.29b | 9.12d | 11.27b | 7.78d | 8.22c |
16% VWC | 15.13ab | 11.49a | 16.88a | 11.68a | 15.62ab | 14.94a | 15.33ab | 15.31a | 14.95ab | 14.79a |
26% VWC | 17.12a | 13.75a | 15.77ab | 13.10a | 16.69a | 15.40a | 16.12a | 15.34a | 15.77a | 15.36a |
Guinea | | | | | | | | | | |
6% VWC | 16.57a | 14.17a | 15.70ab | 11.86a | 13.78c | 12.56b | 11.24c | 10.61b | 10.26c | 10.25b |
16% VWC | 15.38ab | 11.5a | 15.49ab | 13.06a | 14.6bc | 15.30a | 13.87b | 15.05a | 13.13b | 14.57a |
26% VWC | 14.17b | 13.97a | 14.56b | 12.63a | 15.41abc | 15.58a | 14.09b | 15.82a | 14.16ab | 15.58a |
C. Transpiration (mmolH2Om− 2s− 1) |
| 21 DAE | 26 DAE | 31DAE | 36 DAE | 41 DAE | |
Mali | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 1.71a | 1.53a | 1.61ab | 1.58b | 1.35a | 1.24c | 1.20b | 1.17b | 0.74c | 0.73b |
16% VWC | 1.72a | 1.71a | 1.77a | 2.52a | 1.53a | 2.11a | 1.60a | 2.03a | 1.69a | 2.24a |
26% VWC | 1.79a | 1.72a | 1.70ab | 2.54a | 1.62a | 2.09a | 1.47ab | 2.18a | 1.70a | 2.29a |
Guinea | | | | | | | | | | |
6% VWC | 1.71a | 1.71a | 1.45b | 1.71b | 1.36a | 1.46bc | 1.30b | 1.27b | 1.07b | 0.94b |
16% VWC | 1.57a | 1.59a | 1.56ab | 2.27ab | 1.48a | 1.99ab | 1.66a | 1.90a | 1.59a | 1.85a |
26% VWC | 1.62a | 1.74a | 1.52ab | 2.25ab | 1.51a | 2.18a | 1.64a | 2.03a | 1.67a | 2.03a |
D. Stomatal Conductance (molH2Om− 2s− 1) |
Mali | 21 DAE | 26 DAE | 31DAE | 36 DAE | 41 DAE |
Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 0.083a | 0.075a | 0.088c | 0.140a | 0.077c | 0.102c | 0.065c | 0.081d | 0.051c | 0.052b |
16% VWC | 0.089a | 0.087a | 0.109ab | 0.124a | 0.121ab | 0.145a | 0.135a | 0.144a | 0.137ab | 0.140a |
26% VWC | 0.086a | 0.083a | 0.113a | 0.124a | 0.128a | 0.137a | 0.139a | 0.140ab | 0.139a | 0.146a |
Guinea | | | | | | | | | | |
6% VWC | 0.086a | 0.086a | 0.092bc | 0.092a | 0.087c | 0.112bc | 0.073c | 0.097cd | 0.063c | 0.066b |
16% VWC | 0.085a | 0.078a | 0.097abc | 0.124a | 0.116b | 0.129ab | 0.118b | 0.119bc | 0.122b | 0.139a |
26% VWC | 0.081a | 0.085a | 0.103abc | 0.115a | 0.117ab | 0.134ab | 0.119b | 0.128ab | 0.136ab | 0.142a |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence.
plant in trials 1 and 2, respectively, compared to 26%-VWC grown plants. The Guinea accession also had 51% and 32% lower shoot weight at 6%-VWC trial 1 and trial 2, respectively. The effect of drought on shoot growth was clearly visible by 41 DAE (Fig. 4B), though drought-treated plants showed substantial recovery after 20 days of restoration of optimal watering (Fig. 4C).
Morphologically, reduced shoot branching was visible after 10 days of drought initiation (31 DAE) for both accessions (P≤0.05) (Fig. 4D, Table 2A). For the Mali accession, in both trials, the number of branches was 60% lower at 6%-VWC compared to 26%-VWC, whereas for the Guinea accession, the reduction was 43%-53% (P≤0.05) by 41 DAE (Fig. 4D, Table 2A). There was no significant difference between the two accessions in the total number of branches across the moisture levels at 41 DAE. However, for the 6%-VWC treated plants, partial recovery in branching was observed after 20 days of optimal irrigation (61 DAE) in both accessions (Fig. 4D, Table 2A).
Following a similar trend as branching, plant height was also significantly reduced after 10 days of drought (31 DAE) for both accessions (P≤0.05) (Fig. 4E, Table 2B). By 41 DAE, for the Mali accession at 41 DAE, plants exposed to 6%-VWC were 29%-37% (16–19 cm) shorter than the plants exposed to 26%-VWC, and for the Guinea accessions, they were 27%-32% (16 cm) shorter (Fig. 4E, Table 2B). In both accessions, 20 days after restoration of normal watering, plant height of 6%-VWR treated plants remained significantly reduced compared to plants that had been optimally watered throughout the experiment (Fig. 4E, Table 2B). No significant difference in plant height was noted between the two accessions across different moisture levels.
The leaves of the Guinea accession had significantly greater width than the Mali accession across the moisture treatments for all time point measurements except for a few exceptions (Fig. 4F, Table 2C). Similar to previous traits, leaf width was significantly reduced (P≤0.05) after 10 days of drought exposure (at 31 DAE) for both accessions (6% versus 26% VWC). By 41 DAE, the Mali accession had 10%-12% (0.6–0.8 mm) reduced leaf width in plants exposed to 6%-VWC than the plants exposed to 26%-VWC, similar to the 9%-13% (0.6–1.1 mm) reductions observed for the Guinea accession. After normal watering (61 DAE), in both accessions, the width of new leaves from drought-exposed plants recovered, to match those from the other two moisture levels (16% and 26%-VWC) (P≤0.05) (Fig. 4F, Table 2C).
Whereas optimally watered plants (26%-VWC) showed progressively reduced leaf angle as the plants developed, leaf angle significantly increased (p ≤ 0.05) during drought (6%-VWC versus 26%-VWC) for both accessions (Fig. 4G, Table 2D), starting after 10 days of drought initiation. For the Mali accession at 41 DAE, a 50%-83% higher leaf angle was observed in 6%-VWC exposed plants, compared to 26%-VWC treated plants, whereas for the Guinea accession, it was 41%-62% higher. Numerically, the greatest leaf angles observed in the experiment belonged to the Mali accession exposed to drought during the drought exposure at 31–41 DAE (i.e. angles of 50–65˚). For both accessions, after 20 days of normal watering (61 DAE), plants previously exposed to drought had decreased leaf angles but remained the highest among the three moisture treatments (Fig. 4G, Table 2D).
Long term above-ground impacts
Next, we asked to what extent drought had long term impacts on fonio, approximately 8–10 weeks after the drought treatments had ended after optimal watering had resumed. At 96–110 DAE, the accessions were harvested, and shoot and grain traits were measured.
Table 2
Experiment A: Effect of moisture treatments on shoot morphological traits in white fonio.
Mali | 21DAE | 31DAE | 41DAE | 61DAE (Recovery) |
Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 5.90a | 5.60a | 11.0b | 7.80c | 15.60b | 14.50b | 48.00c | 42.60b |
16% VWC | 6.80a | 7.40a | 16.7a | 14.60ab | 37.50a | 32.80a | 67.70ab | 66.00a |
26% VWC | 6.30a | 7.40a | 17.9a | 16.40a | 38.90a | 36.40a | 72.70a | 68.90a |
Guinea | | | | | | | | |
6% VWC | 6.80a | 5.50a | 9.7b | 7.50c | 20.10b | 16.90b | 51.10c | 42.60b |
16% VWC | 6.10a | 6.80a | 17.8a | 13.30ab | 35.90a | 30.40a | 67.20ab | 61.80a |
26% VWC | 6.50a | 6.40a | 17.0a | 12.00b | 35.20a | 35.00a | 64.70b | 62.10a |
A. Number of branches |
B. Plant height (cm) |
Mali | 21 DAE | 31 DAE | 41DAE | 61 DAE (Recovery) |
Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 28.75a | 21.22b | 30.91c | 33.56b | 32.42b | 39.93c | 49.01c | 43.90d |
16% VWC | 30.45a | 24.43ab | 38.13a | 41.82a | 47.36a | 52.42b | 52.83bc | 58.53c |
26% VWC | 30.2a | 26.70a | 40.10a | 42.85a | 51.72a | 55.81ab | 58.96ab | 65.86bc |
Guinea | | | | | | | | |
6% VWC | 30.7a | 24.20ab | 32.80bc | 33.93b | 33.0b | 43.72c | 50.80c | 58.51c |
16% VWC | 29.95a | 27.23a | 37.12ab | 41.74a | 48.4a | 56.24ab | 53.55bc | 81.44a |
26% VWC | 29.10a | 24.45a | 38.80a | 43.91a | 48.5a | 59.73a | 61.61a | 74.41ab |
C. Leaf width (cm) |
Mali | 21DAE | 26DAE | 31DAE | 36DAE | 41DAE | 61DAE (Recovery) |
Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 0.56b | 0.5b | 0.57b | 0.56d | 0.53d | 0.57c | 0.53c | 0.58c | 0.55d | 0.59c | 0.61b | 0.69b |
16% VWC | 0.60b | 0.5b | 0.61b | 0.6cd | 0.65abc | 0.66b | 0.6b | 0.68b | 0.59cd | 0.68b | 0.59b | 0.7b |
26% VWC | 0.57b | 0.5b | 0.61b | 0.63bc | 0.64bc | 0.66b | 0.64ab | 0.67b | 0.61c | 0.67b | 0.60b | 0.71b |
Guinea | | | | | | | | | | | | |
6% VWC | 0.57b | 0.58a | 0.63b | 0.64abc | 0.61c | 0.66b | 0.62b | 0.66b | 0.63bc | 0.68b | 0.69a | 0.74ab |
16% VWC | 0.67a | 0.55ab | 0.70a | 0.66ab | 0.71a | 0.77a | 0.69a | 0.79a | 0.68ab | 0.80a | 0.68a | 0.79a |
26% VWC | 0.67a | 0.59a | 0.71a | 0.70a | 0.70ab | 0.74a | 0.70a | 0.77a | 0.69a | 0.79a | 0.70a | 0.8a |
D. Leaf Angle (º) |
Mali | 21 DAE | 31 DAE | 41DAE | 61 DAE (Recovery) |
Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% VWC | 42.0a | 47.5a | 50.0a | 54.5a | 49.5a | 65.0a | 39.0a | 41.5a |
16% VWC | 43.0a | 45.0a | 39.0b | 41.5bcd | 35.5b | 42.5bc | 34.5bc | 32.0bc |
26% VWC | 43.5a | 50.5a | 35.0b | 49.0ab | 33b | 35.5cd | 32.0c | 35.0b |
Guinea | | | | | | | | |
6% VWC | 42.0a | 47.5a | 46.5a | 48.5bc | 46.5a | 49.5b | 37.0ab | 35.5b |
16% VWC | 39.5a | 44.5a | 36.0b | 37.0d | 34b | 33.0d | 32.5c | 29.0c |
26% VWC | 44.0a | 43.5a | 34.0b | 40.5cd | 33b | 30.5d | 31.5c | 28.0c |
Table S1
A. Experiment A: Concentrations of selected metabolites in white fonio leaves from drought exposed (6%-VWC) and optimally watered (26%-VWC) plants at 41 DAE (end of drought treatment) in the Mali and Guinea accession (Trial 1).
Metabolites (nmol/per mg tissue) | Mali | Guinea | SEM (±) |
VWC6% | VWC26% | VWC6% | VWC26% |
Glucose | 21.9545a | 16.1543a | 20.6480a | 14.3926a | 1.9333 |
Fructose | 9.2732a | 6.1355a | 8.9343a | 6.1801a | 0.9906 |
Succinic acid | 0.1532a | 0.0897a | 0.1648a | 0.1251a | 0.0141 |
Alanine betaine * | 0.0462a | 0.0056ab | 0.023ab | 0.0025b | 0.0064 |
Proline betaine | 0.0002a | 0.0001a | 0.0003a | 0.0002a | 0.0001 |
Glycine betaine* | 3.7862a | 0.4992b | 3.2952a | 0.4927b | 0.3018 |
Ectoine | 0.0001a | 0.0002a | 0.0002a | 0.0002a | 0.0001 |
Spermidine | 0.0027a | 0.0017a | 0.0019a | 0.0021a | 0.0004 |
Spermine | 0.0013a | 0.0009a | 0.0010a | 0.0008a | 0.0001 |
4-OH Pro | 0.0025a | 0.0022a | 0.0025a | 0.0021a | 0.0003 |
GABA | 0.2356a | 0.2110a | 0.2840a | 0.2362a | 0.0221 |
Arg | 0.0416a | 0.0446a | 0.0433a | 0.0548a | 0.0049 |
Asn | 0.3177a | 0.695a | 0.4436a | 0.4656a | 0.0780 |
Asp | 0.8906a | 0.8555a | 0.8742a | 0.9594a | 0.0842 |
Cit | 0.0105a | 0.0107a | 0.0187a | 0.0093a | 0.0036 |
Gln | 0.3994a | 0.2904a | 0.4634a | 0.3132a | 0.0407 |
Glu | 1.3078a | 1.2060a | 1.4276a | 1.4155a | 0.1263 |
Gly | 0.0656a | 0.0672a | 0.0854a | 0.0725a | 0.0107 |
Ile | 0.0386a | 0.0498a | 0.0348a | 0.0525a | 0.0063 |
Leu | 0.1261a | 0.1097a | 0.1208a | 0.1138a | 0.0094 |
Orn | 0.0026a | 0.0021a | 0.0026a | 0.0019a | 0.0003 |
Phe | 0.0897a | 0.0805a | 0.0733a | 0.0757a | 0.0118 |
Pipecolic acid | 0.0659a | 0.0567a | 0.0506a | 0.0338a | 0.0063 |
Pro* | 0.3132a | 0.1330ab | 0.2315ab | 0.1023b | 0.0287 |
Putrescine | 0.0242a | 0.0217a | 0.0167a | 0.0224a | 0.0028 |
Ser | 0.1990a | 0.1850a | 0.2192a | 0.1970a | 0.0183 |
Thr | 0.2983a | 0.2522a | 0.3056a | 0.2958a | 0.0221 |
Trp | 0.0145a | 0.0270a | 0.0091a | 0.0241a | 0.0038 |
Tyr | 0.0551a | 0.0525a | 0.0460a | 0.0548a | 0.0051 |
Val | 0.1080a | 0.1250a | 0.0873a | 0.1254a | 0.0138 |
Different letters within each row (*) symbolize statistical differences at P ≤ 0.10. Abbreviation: VWC, volumetric water content. |
Table S2
Experiment A: Effect of moisture treatments on shoot biomass, root biomass and root/shoot ratio of white fonio at the end of the moisture treatment (at 41 DAE).
VWC | Accession | Shoot dry weight/plant 41 DAE (g) | Root dry weight/plant 41 DAE (g) | Root/shoot ratio 41DAE |
Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% | Mali | 2.98c | 3.02b | 0.42a | 0.56ab | 0.144a | 0.187a |
Guinea | 3.42c | 3.61b | 0.33b | 0.48b | 0.095b | 0.134b |
16% | Mali | 5.65b | 5.03a | 0.48a | 0.63a | 0.087b | 0.126b |
Guinea | 5.91ab | 5.40a | 0.48a | 0.64a | 0.084b | 0.121b |
26% | Mali | 5.77b | 5.44a | 0.50a | 0.60a | 0.087b | 0.109b |
Guinea | 7.04a | 5.37a | 0.53a | 0.66a | 0.077b | 0.124b |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence. |
Table S3
Experiment A: Effect of moisture treatments on harvest index, shoot yield and grain yield traits of white fonio at final harvest (96–110 DAE).
VWC | Accession | Harvest Index (%) | Shoot weight/plant (g) | Grain weight/plant (g) | No of panicles/plant | Panicle length (cm) | Grain weight/ panicle (mg) |
| | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 | Trial 1 | Trial 2 |
6% | Mali | 14.51a | 14.76a | 64.96c | 58.09d | 11.01b | 10.06bc | 153.50ab | 132.50bc | 8.58bc | 7.44b | 75.00a | 77.59a |
Guinea | 6.33c | 7.78c | 76.07c | 72.26c | 5.15c | 6.13d | 122.50b | 127.90c | 10.44a | 11.33a | 43.75b | 48.17b |
16% | Mali | 10.92b | 9.82b | 94.59b | 110.00b | 11.59a | 11.95ab | 146.30b | 149.80ab | 8.33bc | 8.33b | 80.74a | 80.06a |
Guinea | 6.05c | 7.55c | 105.45ab | 109.72b | 6.63c | 8.96c | 159.10ab | 160.20a | 10.2a | 10.93a | 42.12b | 56.15b |
26% | Mali | 11.88b | 11.34b | 100.43b | 104.14b | 13.53a | 13.33a | 192.70a | 170.40a | 8.08c | 7.63b | 73.67a | 78.64a |
Guinea | 7.78c | 7.77c | 114.21a | 122.34a | 9.61b | 10.33bc | 159.40ab | 168.60a | 9.74ab | 11.12a | 60.49ab | 61.83b |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence. |
Table S4
Experiment B: Effect of moisture treatments on crown root architecture (count, diameter, length and surface area) in the Mali accession of white fonio at different time points.
Trial 1 | VWC | 3 DAE | 6 DAE | 9 DAE | 12 DAE | 15 DAE |
No of CR | 6%VWC | 1.00a | 1.00a | 1.00b | 1.25b | 1.25b |
26%VWC | 1.00a | 1.25a | 2.75a | 3.00a | 3.50a |
CR diameter | 6%VWC | 0.18a | 0.19b | 0.19b+ | 0.19b | 0.19b |
26%VWC | 0.19a | 0.25a | 0.28a+ | 0.32a | 0.33a |
Length of the longest CR | 6%VWC | 5.00a | 8.355a | 11.54a | 15.98a | 16.90a |
26%VWC | 4.22a | 7.40a | 8.20a | 10.88b | 10.74b |
Average length of CRs | 6%VWC | 5.00a | 8.36a | 11.54a | 15.48a | 16.12a |
26%VWC | 4.18a | 7.08a | 6.18b | 8.14b | 7.46b |
Root area-diameter > 0.2 mm (cm2) | 6%VWC | 0.09a | 0.15a | 0.18a+ | 0.19b | 0.26b |
26%VWC | 0.12a | 0.24a | 0.48a+ | 0.78a | 1.23a |
Trial 2 | VWC | 3DAE | 6DAE | 9DAE | 12DAE | 15DAE | 18DAE | 21DAE | 24DAE |
No of CR | 6% | 1.00a | 1.00b | 1.00b | 1.00b | 1.25b | 1.25b | 1.50b | 1.50b |
26% | 1.25a | 1.70a | 2.50a | 2.25a | 3.00a | 3.00a | 3.00a | 3.50a |
CR diameter | 6% | 0.16a | 0.15a | 0.15a | 0.16a | 0.14a | 0.15b | 0.15b+ | 0.17b |
26% | 0.18a | 0.18a | 0.22a | 0.25a | 0.25a | 0.29a | 0.28a+ | 0.31a |
Length of the longest CR | 6% | 6.56a | 8.33a | 10.50a | 11.36a | 14.03a | 14.71a | 16.05a | 17.48a |
26% | 5.81a | 6.90a | 8.96a | 9.48a | 11.54a | 11.92a | 11.70b | 12.64b |
Average length of CRs | 6% | 4.59a | 5.45a | 9.62a | 11.50a | 12.85a | 13.03a | 13.85a | 14.75a |
26% | 3.70a | 4.11a | 6.69a | 6.81b | 8.65b | 8.56b | 9.50b | 9.20b |
Root area-diameter > 0.2 mm (cm2)) | 6% | 0.08a | 0.11a | 0.19b | 0.23b+ | 0.43b | 1.28b | 1.80b | 2.52b |
26% | 0.20a | 0.39a | 0.57a | 0.78a+ | 1.77a | 3.60a | 3.83a | 4.04a |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence. |
Table S5
Experiment B: Effect of moisture treatments on lateral root architecture (tips/branching, surface area) in the Mali accession of white fonio at different time points
Trial 1 | VWC | 3 DAE | 6 DAE | 9 DAE | 12 DAE | 15 DAE |
No of tips | 6% | 336.75a | 499.00a | 761.00a+ | 957.50a | 1143.50a |
26% | 346.75a | 437.00a | 582.75b+ | 669.50b | 723.75b |
Root area-diameter ≤ 0.2 mm (cm2) | 6% | 0.69a | 1.65a | 3.76a | 5.37a | 7.34a |
26% | 0.86a | 1.86a | 3.03a | 3.64b | 4.08b |
Trial 2 | VWC | 3DAE | 6DAE | 9DAE | 12DAE | 15DAE | 18DAE | 21DAE | 24DAE |
No of tips | 6% | 454.25a | 565.00a | 827.00a | 999.00a | 1157.00a | 1249.00a | 1349.00a | 1503.25a |
26% | 263.75b | 413.25a | 474.50b | 509.50b | 682.00b | 777.50b | 843.75b | 920.25b |
Root area-diameter ≤ 0.2 mm (cm2) | 6% | 0.71a | 2.24a | 4.55a | 5.47a | 7.45a | 9.64a | 10.75a | 11.24a |
26% | 1.10a | 1.56a | 2.74b | 2.34b | 3.50b | 5.43b | 6.68b | 7.87b |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence. |
Table S6
Experiment B: Effect of moisture treatments on the root hair architecture (length and density) of white fonio at different time points measured on three crown root segments (S1-S3). Root hair measurements on the longest crown root:
Trial 1 | VWC | 3DAE | 6 DAE | 9DAE | 12DAE | 15DAE |
Root hair length at S1 (x100 µm) | 6% | 3.64a | 3.82a | 4.40a+ | 3.50a | 4.56a+ |
26% | 3.80a | 4.05a | 3.05b+ | 3.35a | 3.47b+ |
Root hair length at S2 (x100 µm) | 6% | 4.18a | 0.41a | 4.23a | 4.18a | 5.90a+ |
26% | 3.13a | 3.90a | 3.65a | 3.75a | 3.78b+ |
Root hair length at S3 (x100 µm) | 6% | 2.73a | 2.22a | 3.76a | 4.00a | 4.78a |
26% | 1.50b | 2.03a | 1.69b | 2.71a | 2.12b |
Root hair density at S1 (count/300 µm) | 6% | 9.23a | 7.95a | 8.78a | 9.15a+ | 9.13a |
26% | 6.03a | 5.03a | 4.53a | 5.15b+ | 4.28b |
Root hair density at S2 (count/300 µm) | 6% | 3.72a | 5.69a | 5.82a | 7.00a | 7.54a+ |
26% | 2.47a | 2.44a | 4.12a | 2.75b | 2.94b+ |
Root hair density at S3 (count/300 µm) | 6% | 2.19a | 5.06a | 8.69a | 6.69a | 9.50a |
26% | 0.63a | 1.81a | 2.33b | 2.58b | 2.94b |
Trial 2 | VWC | 3DAE | 6DAE | 9DAE | 12DAE | 15DAE | 18DAE | 21DAE | 24DAE |
Root hair length at S1 (x100 µm) | 6% | 3.69a | 3.39a | 1.72a | 3.11a | 2.12a | 3.21a | 1.96a | 6.10a |
26% | 2.74a | 3.10a | 1.87a | 4.16a | 2.18a | 3.36a | 1.59a | 3.84a |
Root hair length at S2 (x100 µm) | 6% | 2.75a | 2.87a | 1.90a | 2.46a | 2.17a | 3.15a | 2.29a | 4.17a+ |
26% | 3.34a | 3.34a | 1.91a | 3.27a | 1.84a | 2.69a | 1.55b | 3.17b+ |
Root hair length at S3 (x100 µm) | 6% | 2.62a | 2.82a | 2.12a | 2.29b | 1.85a | 2.23 | 1.85a | 3.04a |
26% | 2.56a | 3.37a | 1.69a | 3.02a | 1.70a | NA | 1.40b | 2.00b |
Root hair density at S1 (count/300 µm) | 6% | 11.23a | 8.69a | 7.75a | 6.38a | 7.15a | 6.19a | 9.69a | 11.25a |
26% | 9.17a | 9.21a | 7.75a | 5.75a | 6.63a | 6.38a | 9.42a | 4.59b |
Root hair density at S2 (count/300 µm) | 6% | 8.58a | 3.23a | 7.86a | 2.87a | 4.88a | 7.39a+ | 7.79a | 11.85a |
26% | 6.75a | 6.01a | 4.85a | 3.67a | 2.38a | 3.72b+ | 4.82a | 3.19b |
Root hair density at S3 (count/300 µm) | 6% | 5.81a | 4.71a | 7.00a | 6.63a+ | 5.38a | 5.21a | 7.11a | 7.56a |
26% | 0.42b | 3.13a | 1.25b | 2.94b+ | 3.63a | 0.00b | 2.19b | 0.44b |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence. |
Table S7
Experiment B: Effect of moisture treatments on the root hair architecture (length and density) of white fonio at different time points measured on three crown root segments (S1-S3). Root hair measurement; average of ≤ 4 crown roots/plant:
Trial 1 | VWC | 3DAE | 6DAE | 9DAE | 12DAE | 15DAE |
Root hair length at S1 (x100 µm) | 6%VWC | 3.78a | 3.90a | 4.40a | 3.50a | 4.17a |
26%VWC | 3.65a | 3.75a | 3.45a | 3.03a | 3.43a |
Root hair length at S2 (x100 µm) | 6%VWC | 4.18a | 4.05a | 4.23a | 4.18a | 5.90a |
26%VWC | 3.13a | 3.90a | 3.45a | 3.58a | 3.73b |
Root hair length at S3 (x100 µm) | 6%VWC | 2.75a | 2.95a | 3.76a | 3.00a | 4.78a |
26%VWC | 1.50b | 2.23a | 2.40b | 2.33a | 2.10b |
Root hair density at S1 (count/300 µm) | 6%VWC | 9.19a | 7.94a | 8.75a+ | 9.13a | 9.13a |
26%VWC | 6.00a | 4.47a | 4.08b+ | 3.87b | 4.34b |
Root hair density at S2 (count/300 µm) | 6%VWC | 3.72a | 5.69a | 5.82a | 7.00a | 8.29a |
26%VWC | 2.47a | 2.19a | 3.74a | 2.70b | 3.39b |
Root hair density at S3 (count/300 µm) | 6%VWC | 2.19a | 5.06a | 8.69a | 6.69a | 9.50a |
26%VWC | 0.63a | 1.81a | 2.15b | 1.74b | 3.06b |
Trial 2 | VWC | 3DAE | 6DAE | 9DAE | 12DAE | 15DAE | 18DAE | 21DAE | 24DAE |
Root hair length at S1 (x100 µm) | 6%VWC | 2.92a | 3.23a | 1.77a | 3.11a | 2.12a | 3.21a | 1.93a | 5.65a |
26%VWC | 2.68a | 3.28a | 2.05a | 3.92a | 2.07a | 3.13a | 1.62a | 3.23b |
Root hair length at S2 (x100 µm) | 6%VWC | 2.56a | 2.87a | 1.97a | 2.46a | 2.17a | 3.15a | 2.70a | 4.36a |
26%VWC | 3.07a | 3.34a | 2.01a | 3.38a | 1.87a | 2.84a | 1.59b | 2.78b |
Root hair length at S3 (x100 µm) | 6%VWC | 2.28a | 2.63a | 1.69a | 2.29b | 1.85a | 2.23a | 1.86a | 3.76a+ |
26%VWC | 2.90a | 3.98a | 3.10a | 3.74a | 1.82a | 1.86b | 1.18b | 2.76b+ |
Root hair density at S1 (count/300 µm) | 6%VWC | 9.74a | 7.52a | 8.33a | 6.38a | 7.15a | 6.19a | 10.34a | 11.75a |
26%VWC | 5.69a | 7.71a | 7.17a | 3.56a | 5.71a | 5.92a | 8.26a | 5.68a |
Root hair density at S2 (count/300 µm) | 6%VWC | 8.32a | 4.42a | 7.18a | 3.82a | 4.88a | 7.39a | 8.50a | 10.19a |
26%VWC | 5.63a | 6.01a | 4.11a | 3.33a | 2.55a | 3.15b | 4.743a | 2.97b |
Root hair density at S3 (count/300 µm) | 6%VWC | 5.05a | 4.70a | 4.58a | 6.39a | 7.17a | 4.92a | 6.36a | 6.88a |
26%VWC | 0.67b | 3.33a | 1.87a | 2.02b | 2.95b | 0.14b | 0.68b | 1.37b |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence. |
Supplemental Figures |
Different letters symbolize statistical differences at P ≤ 0.05. Abbreviations: VWC, volumetric water content; DAE, days after emergence.
The harvest index (HI) for the Mali accession was consistently superior across the moisture treatments than the Guinea accession (P ≤ 0.05) in both trials (Fig. 5A, Table S3). Interestingly, for the Mali accession, plants exposed to 6%-VWC had the highest HI (14.5%-14.8%) among all the moisture treatments. By contrast, the Guinea accession had significantly lower HI (6.3%-7.8%) at 6%-VWC, which was not significantly different than other moisture levels. The higher HI for the Mali accession could be partially explained by it having a numerically lesser shoot biomass at harvest compared to the Guinea accession (though inconsistently significant), in combination with having significantly higher grain yield (p ≤ 0.05) (Fig. 5B-C, Table S3).
Drought exposed plants did not fully recover (Fig. 5B, Table S3). The Mali accession that had earlier been exposed to 6% VWC had a 35–44% reduction in shoot weight per plant, as compared to 26% VWC treated plants; for the Guinea accession, it was 33–41% (Fig. 5B, Table S3). Similarly, grain weight per plant for the Mali accession remained decreased by 19%-25% at harvest for plants previously exposed to drought (6%-VWC) compared to those that had been optimally watered (26% VWC) throughout; 46%-41% reductions were noted for the Guinea accession exposed to drought (Fig. 5C, Table S3).
Across all moisture levels, the average grain weight per plant for the Mali accession was significantly greater than the Guinea accession (Fig. 5C, Table S3); for plants that had been continuously grown under optimum moisture (26%-VWC), the Mali accession was 35% greater. To understand this difference, the components of grain yield were analyzed. The number of panicles per plant was not significantly different between the two genotypes across the moisture treatments (P ≤ 0.05) (Fig. 5D, Table S3). However, the earlier drought treatment showed a trend for reduced panicle number per plant. Upon visual examination, the Guinea panicles appeared to be longer than those from the Mali accession, but had lower grain fill (Fig. 5E). Subsequent measurements verified these observations: the panicle length was significantly greater for the Guinea accession than the Mali accession across the moisture treatments (P ≤ 0.05) (Fig. 5D-E, Table S3). Indeed, the Guinea accession had 21–46% longer panicles than the Mali accession under optimum moisture (26%-VWC) (Fig. 5F, Table S3). Soil moisture did not affect panicle length. Furthermore, as predicted, grain weight per panicle was consistently higher for the Mali accession than the Guinea accession across the moisture treatments with one exception (Fig. 5G, Table S3). In fact, grain weight per panicle for the Guinea accession was 18–21% lower than the Mali accession when continuously grown under optimum water (26%-VWC). Grain weight per panicle showed a strong positive correlation (r = 0.6–0.8, P < 0.0001) with grain weight per plant, whereas panicle length had a negative correlation (r = -0.5, P ≤ 0.0001) (Fig. 5H-I).
Principal Component Analysis (PCA)
Multivariate spatial PCA was undertaken to understand the relationships amongst physiological and morphological measurements at 41 DAE, and yield-related traits recorded at harvest (total = 17 traits). Principal components 1 and 2 explained a mean 36% and 30% variation of the data, respectively, and hence were presented (Fig. 7A-D). Several physio-morphological traits measured at 41 DAE clustered including: shoot weight, root weight, leaf width, plant height, photosynthesis, conductance, transpiration, and number of branches; these mostly decreased during drought stress. Along the PC1 axis, leaf angle and root/shoot ratio at 41 DAE at vector angles opposite to the above cluster, consistent with these trait values increasing during drought; along this axis, the harvest index clustering pattern differed between the Mali and Guinea accession (Fig. 7A-D). The traits that consistently clustered with grain weight included: root weight-41DAE, leaf width, shoot weight-41DAE, plant height and SPAD, suggestive of their potential positive correlation. Significant positive correlations (r = 0.7 and 0.8, P ≤ 0.001) between grain weight and SPAD chlorophyll were found in trials 1 and 2, respectively (Fig. 7E). In contrast, a negative correlation between grain weight and leaf width was revealed, although the correlation was not strong (Fig. 7F). Positive and significant but lower correlations with grain weight were found for plant height and root weight-41DAE (Fig. 7G-H).
Underground Effects
Root biomass
To evaluate the potential underground acclimation responses of fonio against drought, root biomass was calculated after the end of the drought treatment (41-DAE). The Guinea accession had a significant (P ≤ 0.05) reduction in root weight (27–37%) from plants exposed to 6%-VWC compared to 26% VWC. By contrast, the Mali accession did not significantly change root biomass after exposure to drought (Fig. 6A-B, Table S2). Furthermore, Mali accession plants exposed to drought (6%-VWC) showed dramatic increases in the root/shoot ratio, which were significantly greater (p ≤ 0.05) than all other treatments, and 40%-52% greater than the Guinea accession (Fig. 6C, Table S2).
Root system architecture
The above measurements suggested that the Mali accession root system was acclimating to help plants survive under drought by prioritizing root growth for water scavenging. When excavated Mali accession roots were visually examined, it appeared as if drought exposure caused roots to be more fibrous, suggestive of a change in root system architecture (RSA) (Fig. 6C). Therefore, a new experiment was designed to examine RSA in detail, focusing only on the Mali accession. Plants were exposed only to the extreme moisture treatments (6% and 26%-VWC) which were initiated at germination, and root system architecture was examined at 3-day intervals for finer temporal resolution, and at macro-, meso- and micro-scales for morphological resolution. In Trial 1, the experiment terminated at 15 DAE, which was extended to 24 DAE in Trial 2:
Macro scale
For the macro scale root system study, crown root number, diameter and length were evaluated (Fig. 8). The number of crown roots from plants exposed to 6%-VWC was significantly lower (2–3 fold) (p ≤ 0.05) than under 26%-VWC starting at 6–9 DAE and continuing until the final data point (Fig. 8A-B, S1, Table S4). A similar trend was observed for the crown root diameter in both trials, except it was significant only after 18 DAE in Trial 2; drought was associated with a 45% reduction in crown root diameter by 24 DAE in Trial 2 (Fig. 8A, C, Table S4). Contrastingly, the length of the longest crown root was greater for drought-exposed plants, which started to be significant at 12 DAE and 21 DAE in Trials 1 and 2, respectively; longest crown root length was 29% greater under drought by 24 DAE in Trial 2 (Fig. 8A, D, Table S4). Looking at the average crown root lengths, the difference between moisture treatments was more visible, starting at 9–12 DAE (Fig. 8A, E, Table S4). At 24 DAE, the average crown root length for drought-exposed plants was 60% greater than the plants watered optimally (Table S4). Total root surface area for roots with a diameter > 0.2 mm, corresponding to crown roots, was significantly greater (p ≤ 0.10) in plants grown at 26%-VWC compared to 6%-VWC starting at 9 DAE which was 38% lower at 24 DAE (Fig. 8F, Table S4).
Meso scale
At the meso scale, corresponding to lateral root branching, the earlier observation (Fig. 6C) that root systems under drought were more fibrous/highly branched than optimally-watered plants, was again observed (Fig. 9A). Upon counting, the number of root tips was significantly greater for drought-treated plants than optimally-watered plants starting at 9 DAE (P ≤ 0.10), becoming 63% greater by 24 DAE (Fig. 9B, Table S5). The total surface area of thin roots (diameter ≤ 0.2 mm) was significantly greater (P ≤ 0.05) in plants exposed to drought than optimal water, starting at 9–12 DAE; by 24 DAE, total thin root surface area was 43% greater under drought (Fig. 9C, Table S5).
Microscale
At the micro scale, light microscope images combined with manual tracing were used to measure root hair length and density on the longest crown roots and an average from the longest ≤ 4 crown roots per plant. For each crown root, root hairs from older (S1), intermediate (S2) and younger (S3) segments were phenotyped (Fig. 10A-C). In total, an estimated 27,000 root hairs were manually traced. Root hair length and density were generally significantly greater (p ≥ 0.05–0.10) under drought than optimal water at the latest time point in each trial, i.e. the longest exposure to drought (Fig. 10B-H). In general, in Trial 1, the trend of increased root hair length and density under drought was observed at earlier time points, than Trial 2 (Fig. 10D-H).
Numerically, at the end of trial 1 (15 DAE), root hair lengths on the longest crown root were, respectively, 31, 56 and 125% greater in S1, S2, and S3 for the plants exposed to drought, whereas at the end of trial 2 (24 DAE) they were 59, 32 and 52% greater (Fig. 10D-F, Table S6). Similar observations were made for root hairs averaged across crown roots (Fig. S2A-C, Table S7).
Similarly, in all crown root segments, there was a > 2-fold increase in root hair density by 15 DAE in trial 1 and 24 DAE in trial 2 for plants grown under drought (Fig. 10G-H); it became significantly greater starting at 9 DAE in S3 (Table S6). Root hair density, when averaged across crown roots, increased 2–5 fold in S2 and S3 by 24 DAE (Fig. S2D-F, Table S7).