# Agromorphological Characterization Revealed Three Phenotypic Groups in a Region-Wide Germplasm of Fonio (Digitaria exilis (Kippist) Stapf) from West Africa

Abdou R. Ibrahim Bio Yerima, E. G. Achigan-Dako, Aissata Mamadou, Emmanuel Sekloka, Claire Billot, Charlotte O. A. Adje, Adeline Barnaud, Yacoubou Bakasso

► To cite this version:

Abdou R. Ibrahim Bio Yerima, E. G. Achigan-Dako, Aissata Mamadou, Emmanuel Sekloka, Claire Billot, et al. Agromorphological Characterization Revealed Three Phenotypic Groups in a Region-Wide Germplasm of Fonio (Digitaria exilis (Kippist) Stapf) from West Africa. Agronomy, 2020, 10 (11), 10.3390/agronomy10111653. doi:10.3390/agronomy10111653

HAL Id: hal-03126724

https://hal.inrae.fr/hal-03126724

Submitted on 1 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Agromorphological Characterization Revealed Three Phenotypic Groups in a Region-Wide Germplasm of Fonio (Digitaria exilis (Kippist) Stapf) from West Africa

Abdou R. Ibrahim Bio Yerima<sup>1,2</sup>®, Enoch G. Achigan-Dako<sup>1,*</sup>®, Mamadou Aissata<sup>2</sup>, Emmanuel Sekloka<sup>3</sup>®, Claire Billot<sup>4,5</sup>, Charlotte O. A. Adje<sup>1</sup>, Adeline Barnaud<sup>6</sup>® and Yacoubou Bakasso<sup>7</sup>

¹ Laboratory of Genetics, Biotechnology and Seed Sciences (GBioS), Faculty of Agronomic Sciences (FSA), University of Abomey-Calavi (UAC), Cotonou 01 BP 526, Benin; abdourazakoui01@gmail.com (A.R.I.B.Y.); abikecharlotte@gmail.com (C.O.A.A.)

² Department of Rainfed Crop Production (DCP), National Institute of Agronomic Research of Niger (INRAN), Niamey BP 429, Niger; intaaoudi@gmail.com

³ Laboratory of Phytotechny, Plant Breeding and Plant Protection (LaPAPP), Department of Sciences and Techniques of Vegetal Production (STPV), Faculty of Agronomy, University of Parakou, Parakou BP 123, Benin; emmanuelsekloka@hotmail.com

⁴ Unité Mixte de Recherche Amélioration Génétique et Adaptation des Plantes (AGAP), Agricultural Research Centre for International Development (CIRAD), F-34398 Montpellier, France; claire.billot@cirad.fr

5 Amélioration Génétique et Adaptation des Plantes (AGAP), University of Montpellier, CIRAD, Institut National de la Recherche Agronomique (INRA), l'Institut Agro/Montpellier SupAgro, 2 Place Pierre Viala, 34060 Montpellier, France

6 Diversité, Adaptation et Développement des Plantes (DIADE), Institut de Recherche pour le Développement (IRD), University of Montpellier, 34060 Montpellier, France; adeline.barnaud@ird.fr

7 Faculty of Science and Techniques, University of Abdou Moumouni of Niamey, Niamey BP 10662, Niger; bakasso@yahoo.com

Received: 5 August 2020; Accepted: 17 September 2020; Published: 27 October 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Fonio is an ancient orphan cereal, cultivated by resource-poor farmers in arid and semi-arid regions of West Africa, who conserved and used the cereal for nutrition and income generation. However, the extent and patterns of phenotypic variation within the fonio germplasm remain scarcely evaluated to inform breeding decisions. In this study, we used alpha lattice design to assess the phenotypic variability of 180 fonio accessions using 20 agro morphological descriptors, including both qualitative and quantitative traits at Bengou research station in 2018 and 2019. Descriptive statistics, combined analysis and multivariate analyses, including principal components analysis (PCA), hierarchical clustering, and multiple factor analysis (MFA) were performed. As results, we found significant differences among fonio accessions and among agro ecological zones of provenance for most of the quantitative traits. Furthermore, highly significant positive correlations were found between grain yield and other yield-related attributes such as harvest index and thousand seeds weight, whereas significant negative correlations were found between grain yield and flowering and maturity times. Clusters analysis revealed three phenotypic groups. Group 1 was characterized by early maturing fonio accessions with higher grain yield. Groups 2 and 3 were characterized by late maturing accessions with intermediate to moderate grain yield. The accessions from Group 1 are candidate for yield improvement and development of fonio lines with enhanced lodging resistance. Accessions from Group 2 and 3 can be improved for yield through marker-assisted selection of best thousand seed weight. This study highlights how traits are correlated within the major phenotypic groups of fonio in West Africa, and we discussed how these groups could be further exploited for

<footnote>* Correspondence: e.adako@gmail.com; Tel.: +229-95-393283 or +227-96-401486</footnote>
improving traits of economic importance. Furthermore, this study exhibited agro morphological
descriptors that discriminate fonio accessions and provide useful information for parental selection
with economically important agronomic traits.

**Keywords:** fonio; fonio millet; white fonio; *Digitaria exilis*; agro morphological descriptors; phenotypic diversity; neglected and underutilized species (NUS); genetic improvement

**1. Introduction**

Fonio (*Digitaria exilis* (Kippist) Stapf), also referred to as white fonio or fonio millet, is one of the oldest cereal crops domesticated by farmers in West Africa. Its cultivation seems to have started about 7000 years ago [1]. Fonio plays a critical role for food and nutrition security for several millions of people in this region, especially during shortage period due to short time life cycle [2]. Besides, its grain is nutritiously valuable due to the presence of two human-vital amino acids methionine and cysteine, which are absent in other staple food cereals such as wheat (*Triticum aestivum* L.), rice (*Oryza sativa* L.), maize (*Zea mays* L.), or sorghum (*Sorghum bicolor* (L.) Moench) [3]. Additionally, some antithyroid properties have been reported due to high flavonoid content in the crude fonio grain [4]. The grains are poor in gluten and beneficial for people suffering from coeliac disease and overweight people, pregnant and breast feeding women, as well as diabetic people [5]. Furthermore, an efficient formulation of fonio flour enriched with local food resources for the complementary feeding of children under age two has been reported [6]. Due to its C4 metabolism, the crop resists to drought and is well adapted to local pedo-climatic conditions [7].

In traditional fonio farming systems of West Africa, all cultivated varieties are landraces [8]. Their productivity is limited due to the lack of improved varieties and good agronomic practices, lack of institutional support, lack of knowledge about some biotic and abiotic stresses, and tedious post-harvest activities [9]. The development of improved varieties well adapted to different agro ecological zones and preferred by farmers and consumers is of paramount importance in fonio cultivation systems [10]. This can only be achieved by implementing efficient genetic resources management, conservation and valorization, a prerequisite for any crop improvement. This requires a comprehensive knowledge of uses and of diversity patterns. Plant genetic diversity helps breeders in developing improved varieties with high quality, tolerant or resistant to biotic and abiotic stresses, and well adapted to different climate conditions [11].

Morphological markers based on phenotypic characterization of traits are among methods
used for assessing diversity within and between plant populations because they are simple, direct,
easy and cheap to characterize germplasm accessions. Besides, they are directly related to farmers'
perception. Even though these type of markers may present a low level of polymorphism, heritability,
and expression and are often susceptible to phenotypic plasticity [12], they were used by many authors
to study agro morphological diversity in other small grain cereals [13–16]. Various studies reported
the existence of a wide range of diversity in fonio germplasm accessions [17–21] for quantitative
traits (plant height, number of leaves per plant, internode length, days to 50% flowering, days to 50%
maturity, leaf width, panicle leaf length, stem girth, raceme length, dry biomass yield and grain yield,
seeds size, thousand seed weight) and for qualitative traits (collar color, green color of foliar limb,
anthocyanin coloration and distribution in different aerial organs, type of panicle and panicle exertion,
and grain color). However, previous works on fonio agromorphological characterization limited their
collection area to only one or two countries, or used few accessions, overlooking the wider distribution
range of the crop. Hence, this limits breeders' chance to get genes of interest, particularly in an orphan
crop with a narrow knowledge of the genetic base.

So far, and for the last two decades, the following questions remain unresolved: (a) how is the
phenotypic variation distributed among West African fonio accessions? (b) Which agro morphological
traits contribute to distinctiveness in fonio germplasm? (c) How are they related? (d) How many phenotypic groups can be identified? The objectives of this study are to assess phenotypic variability among fonio accessions using agro morphological descriptors, to investigate association patterns among quantitative traits and to identify useful phenotypic groups for future improvement. We hypothesize that exploring a large set of fonio accessions collected from various origins and adaptation areas will enrich the existing germplasm. This will equally help breeders in identifying useful landraces with key agronomic traits for future crop improvement.

## 2. Materials and Methods

### 2.1. Plant Material

A total of 180 fonio accessions were used in this study. They were collected from five Western African countries including Benin (50 accessions), Burkina Faso (five accessions), Guinea (33 accessions), Mali (49 accessions), and Niger (43 accessions), maintained and available in the genebank of the Laboratory of Genetics, Biotechnology and Seed Science (GBioS) from University of Abomey-Calavi (Benin) in 2018. A total of 57 accessions were collected from the arid region, 93 from the semi-arid region, and 30 from the sub-humid region (Figure 1). The arid zone includes the Sahel with up to 750 mm of rainfall in a single short season followed by a long dry season. The vegetation is made up of grass, bush and thicket. In the semi-arid region, the annual rainfall varies from 750 mm to 1250 mm on average in one season followed by a dry season; the vegetation is a savanna grassland. The sub humid region is a deciduous forest and woodland savanna with an average annual rainfall varying between 1250 mm to 1500 mm in one season.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Distribution of collected fonio accessions according to countries and agro ecological zones.</figcaption></figure>

### 2.2. Description of Site

The field trial was carried out at Bengou research station (National Institute of Agronomic Research of Niger) during the rainy seasons 2018 and 2019. This site is in Northern Sudan Savanna
agroecological zone located at 11°58′44.4″ N, 3°32′51.6″ E, at 172 m above sea level and is at 260 km
from Niamey and at 17 km form Republic of Benin border. The soil was sandy loam with basic pH
medium. Annual rainfall was 1047.62 mm and 1211.2 mm in 2018 and 2019 respectively (Figure 2).
The temperature varied between 22.90 °C and 35.20 °C with an average of 29.05 °C. The District of
Bengou offers rich and varied vegetation where we can meet almost all the plant species of Niger.
The plant population can be subdivided according to the topography and the nature of the soils into
three groups: (i) the group of lateritic plateaus is made up of shrub from savannah dominated by
species such as *Combretum nigricans* [22], *Combretum micranthum* G.Don [23], *Acacia macrostachya* Rchb.
ex DC [24], etc.; (ii) the group of valleys mainly consisted of wooded savanna and classified forests
dominated by *Borassus aethiopium*, *Hyphaene thebaica* (L.) Mart [25], *Diospyros mespiliformis* Hochst. ex
A.DC, *Parinari macrophylla* (Sabine) Prance ex F.White [26], etc.; and (iii) the group of the old fixed dunes.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Cumulative rainfall at Bengou research station in 2018–2019.</figcaption></figure>

2.3. *Experimental Design and Management*

The experiment design used was alpha lattice (18 × 10) with two replicates and each accession had four rows per plot of 1 m length and 0.20 m apart. The plot area was 1 m² and distance between blocks and a replicate was 0.5 m and 1 m respectively. The observations were taken from randomly selected ten sampled plants from the four rows of a plot. Prior to planting, land was ploughed, harrowed and fonio seeds were sown (line sowing) at the rate of 2 g per plot on 19 July 2018 and 18 July 2019. Weeds were rouged out manually three weeks after planting and no fertilizer was applied. Weeding was subsequently done in order to keep the trials clean till harvest. No irrigation was applied during the experiment.

## 2.4. *Data Collection*

Data were collected on twenty traits including five qualitative (vigour at seedling. panicle exertion, panicle type, panicle leaf attitude and grain color) and 15 quantitative (plant height, number of internodes per plant, number of leaves per plant, leaf length, leaf width, culm diameter, days to 50% maturity, raceme length, number of racemes per panicle, number of grain par raceme, panicle weight, shoot dry matter yield, grain yield. harvest index and thousand seeds weight) (Table 1). These data were recorded following the Bioversity International rice descriptors list [27].
<table><caption>**Table 1.** Qualitative and quantitative descriptors for agro morphological characterization of fonio accessions.</caption><thead><tr><th colspan="3">Qualitative</th></tr><tr><th>Descriptors</th><th>Code</th><th>Description</th></tr></thead><tbody><tr><td>Vigour at seedling</td><td>1 = good establishment<br>2 = moderate establishment<br>3 = poor establishment</td><td>Recorded three weeks after planting</td></tr><tr><td>Panicle exertion</td><td>1 = well exerted<br>3 = moderately exerted<br>5 = lightly exerted<br>7 = partially exerted</td><td>Recorded near maturity</td></tr><tr><td>Panicle type</td><td>1 = compact<br>3 = intermediate<br>5 = open</td><td>Recorded near maturity</td></tr><tr><td>Panicle leaf attitude</td><td>1 = erect<br>3 = intermediate<br>5 = horizontal<br>7 = descending</td><td>Recorded after heading</td></tr><tr><td>Phenotypic grain colour</td><td>1 = brown<br>2 = greyed orange</td><td>Recorded after harvesting</td></tr><tr><th colspan="3">Quantitative</th></tr><tr><th>Descriptors</th><th>Code</th><th>Description</th></tr><tr><td>Plant height (cm)</td><td>PHT</td><td>Measured from the soil level to the tip of longest panicle</td></tr><tr><td>Number internode</td><td>NIN</td><td>Mean number of internodes on the main stem from 10 plants</td></tr><tr><td>Internode length (cm)</td><td>INL</td><td>Measured from the 2nd internode to last internode below the panicle</td></tr><tr><td>Number of leaves/plant</td><td>NFP</td><td>Mean number of leaves from 10 plants</td></tr><tr><td>Culm diameter (mm)</td><td>CD</td><td>Measured at the mid portion of the culm.</td></tr><tr><td>Days to 50% flowering (days)</td><td>FLO</td><td>Date on which 50% of the plants are flowering</td></tr><tr><td>Days to 50% maturity (days)</td><td>MAT</td><td>Date on which 50% of the plants are maturing</td></tr><tr><td>Number of racemes/panicle</td><td>NRP</td><td>Mean number of raceme from 10 panicles</td></tr><tr><td>Raceme length (cm)</td><td>RLT</td><td>Length of main axis measured from the raceme base to the tip</td></tr><tr><td>Number of grains/raceme</td><td>NGR</td><td>Mean number of grain from 10 racemes</td></tr><tr><td>Panicle length (cm)</td><td>PLT</td><td>Measured from panicle leaf insertion to the tip of longest raceme</td></tr><tr><td>Dry biomass yield (kg.ha<sup>-1</sup>)</td><td>DMY</td><td>Ratio of shoot dry weight (kg) to the plot area (ha)</td></tr><tr><td>Grain yield (kg.ha<sup>-1</sup>)</td><td>GRY</td><td>Ratio of grain weight (kg) to the plot area (ha)</td></tr><tr><td>Harvest index (%)</td><td>HI</td><td>Ratio of grain yield to total dry biomass × 100</td></tr><tr><td>Thousand seeds weight (g)</td><td>TSW</td><td>Counting and weighting of 1000 seeds</td></tr></tbody></table>
## 2.5. Data Analysis

Descriptive statistics were used to depict variations that emerge from qualitative and quantitative traits. When conditions are met, analysis of variance was performed for quantitative traits in order to firstly verify differences among accessions from various provenances and thereafter from agro ecological zones. Pearson's correlation analysis was carried out to investigate association patterns among quantitative traits. In order to describe the variation of the set of variables together principal components analysis (PCA) [28] was carried out using packages "FactoMineR" [29] and Factoextra [30] in R software version 3.5.1 [31]. Moreover, hierarchical clustering on principal component (HCPC) [32] analysis and linear discriminant analysis (LDA) [33] were used to define hierarchical typologies of fonio accessions and to describe phenotypic groups obtained. Multiple factor analysis (MFA) [34] was used for synthetic comparison of the clusters or phenotypic groups obtained.

# 3. Results

## 3.1. Variation of Qualitative Traits among Fonio Accessions

Proportion distributions of the qualitative traits (Table 2) revealed that three weeks after planting 30% of accession showed good seedling vigor, while 63% and 7% showed moderate and poor vigor, respectively. Three types of panicle were observed among accessions: compact type (50%), intermediate type (48%), and open type (2%). Panicle leaf attitude was erect (48%), intermediate (45%), or horizontal (7%). Moreover, three types of panicle exertion were noticed with 65% of accessions showing well exerted panicle, 27% moderately exerted, and 8% lightly exerted. Majority of fonio accessions' seed had 84% of brown color, while 16% had greyed orange color (Figure 3a,b).

<table><caption>**Table 2.** Qualitative parameters variations among fonio accessions.</caption><thead><tr><th>Descriptors</th><th>Modalities</th><th>Number of Accessions</th><th>Percentage (%)</th></tr></thead><tbody><tr><td rowspan="3">Vigour at seedling</td><td>Good</td><td>54</td><td>30</td></tr><tr><td>Moderate</td><td>114</td><td>63</td></tr><tr><td>Poor</td><td>12</td><td>7</td></tr><tr><td rowspan="3">Panicle leaf attitude</td><td>Erect</td><td>86</td><td>48</td></tr><tr><td>Intermediate</td><td>82</td><td>45</td></tr><tr><td>Horizontal</td><td>12</td><td>7</td></tr><tr><td rowspan="3">Panicle type</td><td>Compact</td><td>90</td><td>50</td></tr><tr><td>Intermediate</td><td>86</td><td>48</td></tr><tr><td>Open</td><td>4</td><td>2</td></tr><tr><td rowspan="3">Panicle exertion</td><td>Well exerted</td><td>117</td><td>65</td></tr><tr><td>Moderately exerted</td><td>48</td><td>27</td></tr><tr><td>Lightly exerted</td><td>15</td><td>8</td></tr><tr><td rowspan="2">Phenotypic grain color</td><td>Brown</td><td>152</td><td>84</td></tr><tr><td>Greyed orange</td><td>28</td><td>16</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Grain color in fonio accessions, (a): Greyed orange, (b): Brown.</figcaption></figure>

## 3.2. Fonio Performance Analysis Based on Country of Provenance

Table 3 presents the mean values of quantitative traits of fonio from the different countries. The accessions from Guinea produced significantly taller plant (81 ± 4.92 cm), higher number of internode (11 ± 2.04), and higher internode length (24.24 ± 11.6 cm). This was followed by accessions from Benin with 72 ± 4.62 cm plant height and 10 ± 1.76 internodes. Similar trend was observed for raceme length trait distribution. Accessions with longer racemes were frequent in Guinea (5.52 ± 0.59 cm) and Benin (5.10 ± 1.70 cm) whereas accessions with lower raceme were from Niger and Mali with 4.50 ± 0.60 cm and 4.41 ± 0.93 cm length respectively. However, shorter plants with reduced number of internodes were more predominant among accessions from Niger and Mali with respectively 70 ± 4.98 cm plant height and 7 ± 0.62 internodes, and 68 ± 6.27 cm plant height and 6 ± 0.78 internodes. Accessions from Burkina exhibited medium raceme length (4.80 ± 1.16 cm).

Flowering times (50% of plants flowering) in Guinean fonio accessions was 85 ± 7.75 days, and maturity happens after 105 ± 4.57 days. Flowering time in Malian accession is shorter with 69 ± 1.06 days only. The accessions from Benin exhibited higher number of grain per raceme (105 ± 12.03 grains) and higher dry biomass yield (5095 ± 145.3 kg.ha⁻¹) followed by Guinea with 94 ± 11.9 grains per raceme and 5075 ± 96.47 kg.ha⁻¹ of dry biomass yield, respectively. Accessions of Burkina Faso had 86 ± 22.3 grains per raceme and 4063± kg.ha⁻¹ of dry biomass. Accessions from Niger and Mali exhibited lower values with respectively 79 ± 9.18 and 76 ± 11.18 grains per raceme, and 4536 ± 125.8 kg.ha⁻¹ and 3398 ± 98.27 kg.ha⁻¹ for dry biomass yield. In contrast, accessions from Mali and Niger had higher number of leaves per plant with respectively 212 ± 70.02 leaves and 186 ± 64.9 leaves. In addition, accessions from Niger produced higher grain yield (1936 ± 62.25 kg.ha⁻¹) compared to those from Mali (1910 ± 68.75 kg.ha⁻¹), Burkina (1063 ± 70.53 kg.ha⁻¹), Benin (1057 ± 45.99 kg.ha⁻¹), and Guinea (924 ± 45.15 kg.ha⁻¹). Higher harvest index and larger seed size were observed in accessions from Mali and Niger (Table 3).
<table><caption>**Table 3.** Mean performances of fonio accessions in different countries.</caption><thead><tr><th rowspan="2">Quantitative Traits</th><th colspan="7">Countries</th></tr><tr><th>Benin</th><th>Burkina Faso</th><th>Guinea</th><th>Mali</th><th>Niger</th><th>Mean</th><th>SD</th><th>Pr (&gt;F)</th></tr></thead><tbody><tr><td>Plant height (cm)</td><td>72 ± 4.62b</td><td>71 ± 12.29bc</td><td>81 ± 4.92a</td><td>68 ± 6.27c</td><td>70 ± 4.98bc</td><td>72</td><td>5.03</td><td>***</td></tr><tr><td>Number of internodes</td><td>10 ± 0.85ab</td><td>10 ± 1.76b</td><td>11 ± 2.04a</td><td>6 ± 0.78c</td><td>7 ± 0.62c</td><td>9</td><td>2.17</td><td>***</td></tr><tr><td>Internode length (cm)</td><td>21.42 ± 6.1ab</td><td>22.62 ± 9.6ab</td><td>24.24 ± 11.6a</td><td>19.44 ± 9.3b</td><td>18.44 ± 3.9b</td><td>20.73</td><td>2.34</td><td>**</td></tr><tr><td>Number of leaves/plant</td><td>92 ± 10.8c</td><td>81 ± 10.7c</td><td>79 ± 8.1c</td><td>212 ± 70.2a</td><td>186 ± 64.9b</td><td>144</td><td>63.85</td><td>***</td></tr><tr><td>Culm diameter (mm)</td><td>1.48 ± 0.55a</td><td>1.40 ± 0.13ab</td><td>1.36 ± 0.20b</td><td>1.21 ± 0.11c</td><td>1.29 ± 0.09b</td><td>1.33</td><td>0.10</td><td>**</td></tr><tr><td>Days to flowering</td><td>81 ± 3.24b</td><td>80 ± 2.30b</td><td>85 ± 7.95a</td><td>69 ± 1.06d</td><td>74 ± 2.20c</td><td>77</td><td>6.30</td><td>***</td></tr><tr><td>Days to 50% Maturity</td><td>105 ± 3.75a</td><td>98 ± 3.87b</td><td>105 ± 4.57a</td><td>85 ± 5.07c</td><td>84 ± 1.94c</td><td>95</td><td>10.36</td><td>***</td></tr><tr><td>Number of racemes/panicle</td><td>4 ± 0.21a</td><td>4 ± 0.11a</td><td>3 ± 0.14b</td><td>3 ± 0.32b</td><td>3 ± 0.25b</td><td>3</td><td>0.55</td><td>***</td></tr><tr><td>Raceme length (cm)</td><td>5.10 ± 1.70b</td><td>4.80 ± 1.16bc</td><td>5.52 ± 0.59a</td><td>4.41 ± 0.93c</td><td>4.50 ± 0.60c</td><td>4.82</td><td>0.46</td><td>***</td></tr><tr><td>Number of grains/raceme</td><td>105 ± 12.03a</td><td>86 ± 22.3bc</td><td>94 ± 11.9b</td><td>76 ± 11.18c</td><td>79 ± 9.18c</td><td>88</td><td>11.77</td><td>***</td></tr><tr><td>Panicle length (cm)</td><td>15.1 ± 2.56c</td><td>15.9 ± 1.73bc</td><td>15.8 ± 1.07c</td><td>17.2 ± 3.2ab</td><td>17.4 ± 1.40a</td><td>16.4</td><td>0.98</td><td>**</td></tr><tr><td>Dry biomass yield (kg.ha<sup>-1</sup>)</td><td>5095 ± 145.3a</td><td>4603 ± 159.3ab</td><td>5075 ± 96.47a</td><td>3398 ± 98.27c</td><td>4536 ± 125.8b</td><td>4490</td><td>689.64</td><td>***</td></tr><tr><td>Grain yield (kg.ha<sup>-1</sup>)</td><td>1057 ± 45.99b</td><td>1063 ± 70.53b</td><td>924 ± 45.15b</td><td>1910 ± 68.75a</td><td>1936 ± 62.25a</td><td>1474</td><td>500.69</td><td>***</td></tr><tr><td>Harvest index (%)</td><td>26 ± 4.90c</td><td>27 ± 8.98c</td><td>25 ± 6.23c</td><td>59 ± 9.48a</td><td>54 ± 7.18b</td><td>41</td><td>16.81</td><td>***</td></tr><tr><td>Thousand seed weight (g)</td><td>0.59 ± 0.05d</td><td>0.63 ± 0.10c</td><td>0.68 ± 0.04b</td><td>0.82 ± 0.04a</td><td>0.69 ± 0.03b</td><td>0.70</td><td>0.09</td><td>***</td></tr></tbody></table>

SD = standard deviation. Means followed by the same letter are not significantly different at *p* = 0.05, *** = significant at 0.001, ** = significant at 0.01.
### 3.3. Fonio Performance Analysis Based on Ecological Zones

Fonio accessions from arid zone produced significantly higher grain yield (2937 ± 327.81 kg.ha⁻¹), higher harvest index (60%), and higher thousand seed weight (0.53 ± 0.07 g) (Figure 4m-o). Accessions from sub humid zones revealed lower value with 1364 ± 246.58 kg.ha⁻¹ grain yield, 30% harvest index, and 0.42 ± 0.05 g thousand seed weight. Furthermore, 50% flowering and 50% maturity times were shorter in accessions from arid zone (72 ± 2.95 days and 93 ± 4.03 days, respectively), whereas accessions from sub humid areas exhibited longer 50% flowering (86 ± 4.55 days) and 50% maturity (108 ± 5.19 days) times (Figure 4f,g). Accessions from the semi-arid zone showed relatively an intermediate 50% flowering and 50% maturity times. Highest number of internode (8 ± 1.19) and the higher internode length (25.62 ± 5.56 cm) were observed in accessions from the sub humid zone. A similar trend was observed for raceme length (9.37 ± 3.97 cm), the number of grain per raceme (88 ± 13.60 grains), panicle length (18.3 ± 1.16 cm), and dry biomass yield (5818 ± 543.87 kg.ha⁻¹). However, accessions with wider culm diameter (1.32 ± 0.26 mm) were abundant in the semi-arid zone compared to those from sub humid and arid zones (Figure 4e).

### 3.4. Relationship among Quantitative Traits

We observed highly positive correlations between grain yield (GRY) and harvest index (HI) (r = 0.89), between maturity time (MAT) and flowering time (FLO) (r = 0.80), and between maturity (MAT) and raceme length (RLT) (r = 0.81). Likewise, strong significant positive correlations were observed between number of internode (NIN) and maturity (MAT) (r = 0.77) and between grain yield (GRY) and number of leaves per plant (NLP) (r = 0.68), and thousand seeds weight (TSW) (r = 0.60). However, negative correlations (highly significant and moderately significant) were observed between grain yield (GRY) and flowering time (FLO) (r = -0.79), maturity time (MAT) (r = -0.77), number of internode (NIN) (r = -0.76), and plant height (PHT) (r = -0.35). Moreover, strong significant negative correlations were found between harvest index (HI) and flowering (FLO) (r = -0.84), and maturity (MAT) (r = -0.79) (Table 4).

<table><caption>**Table 4. Pearson's correlation matrix on fifteen quantitative traits used in fonio genotypes description.**</caption><thead><tr><th></th><th>PHT</th><th>NIN</th><th>INL</th><th>NLP</th><th>CD</th><th>MAT</th><th>FLO</th><th>PLT</th><th>RLT</th><th>NRP</th><th>NGR</th><th>DBY</th><th>GRY</th><th>HI</th><th>TSW</th></tr></thead><tbody><tr><td>PHT</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NIN</td><td>0.53</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>INL</td><td>0.55</td><td>0.77</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NLP</td><td>-0.38</td><td>-0.75 *</td><td>-0.70</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CD</td><td>0.35</td><td>0.53</td><td>0.51</td><td>-0.51</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MAT</td><td>0.36</td><td>0.77 **</td><td>0.70</td><td>-0.71</td><td>0.51</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>FLO</td><td>0.35</td><td>0.77</td><td>0.64</td><td>-0.71</td><td>0.56</td><td>0.80 **</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PLT</td><td>0.46</td><td>0.25</td><td>0.34</td><td>-0.19</td><td>0.02</td><td>0.24</td><td>0.05</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>RLT</td><td>0.57</td><td>0.74</td><td>0.73</td><td>-0.64</td><td>0.42</td><td>0.81 **</td><td>0.66</td><td>0.49</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NRP</td><td>-0.02</td><td>0.21</td><td>0.16</td><td>-0.27</td><td>0.31</td><td>0.27</td><td>0.23</td><td>-0.21</td><td>0.13</td><td>1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NGR</td><td>0.43</td><td>0.69</td><td>0.64</td><td>-0.60</td><td>0.50</td><td>0.64</td><td>0.67</td><td>0.15</td><td>0.61</td><td>0.30</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>DBY</td><td>0.48</td><td>0.65</td><td>0.52</td><td>-0.50</td><td>0.43</td><td>0.51</td><td>0.66</td><td>0.17</td><td>0.53</td><td>0.11</td><td>0.59</td><td>1</td><td></td><td></td><td></td></tr><tr><td>GRY</td><td>-0.35 **</td><td>-0.76 **</td><td>-0.70</td><td>0.68 *</td><td>-0.42</td><td>-0.77 *</td><td>-0.79 **</td><td>-0.18</td><td>-0.69</td><td>-0.17</td><td>-0.56</td><td>-0.43</td><td>1</td><td></td><td></td></tr><tr><td>HI</td><td>-0.37</td><td>-0.76</td><td>-0.74</td><td>0.69</td><td>-0.48</td><td>-0.79 **</td><td>-0.84 **</td><td>-0.21</td><td>-0.69</td><td>-0.16</td><td>-0.61</td><td>-0.62</td><td>0.89 **</td><td>1</td><td></td></tr><tr><td>TSW</td><td>-0.24</td><td>-0.63</td><td>-0.45</td><td>0.61</td><td>-0.53</td><td>-0.64</td><td>-0.79</td><td>0.08</td><td>-0.49</td><td>-0.34</td><td>-0.62</td><td>-0.60</td><td>0.60 **</td><td>0.66</td><td>1</td></tr></tbody></table>

PHT = plant height. NIN = number of internodes, INL = internode length, NFP = number of leaves per plant, CD = culm diameter, FLO = 50% flowering, MAT = 50% maturity, NRP = number of racemes per panicle, PLT = panicle length, RLT = raceme length, NGR = number of grain raceme, DBY = dry biomass yield, GRY = grain yield, HI = harvest index, TSW = thousand seeds weight, * = significance at 0.05, ** = significance at 0.01.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Boxplot showing fonio accessions performance based on ecological zones. Variation of: (a) plant height, (b) number of internodes, (c) internodes length, (d) number of leaves, (e) culm diameter, (f) 50% flowering, (g) 50% maturity, (h) number of racemes per panicle, (i) raceme length, (j) number of grains per raceme, (k) panicle length, (l) dry biomass yield, (m) grain yield, (n) harvest index, (o) 1000 seeds weight.</figcaption></figure>
### 3.5. Clustering of Accessions Based on Country of Provenance

The results of PCA showed that the first two dimensions accounted for 50.90% of the total variation among fonio accessions. The first dimension (Dimension 1) was highly correlated with flowering time, number of internodes, maturity time, internode length, raceme length, number of grain per raceme, dry biomass yield, thousand seed weight, number of leaves per plant, grain yield and harvest index while the second dimension (Dimension 2) was moderately correlated to number of raceme per panicle, panicle length, phenotypic grain color, culm diameter, and plant height (Figure 5). Furthermore, the results of principal component analysis (PCA) showed that fonio accessions from Benin, Burkina Faso, and Guinea were characterized by late flowering and maturing accessions, with taller plants, long raceme length and higher number of grains per raceme compared to the group of accessions from Mali and Niger that were characterized by early flowering and early maturing accessions, with higher grain yield, higher thousand seeds weight and higher number of leaves per plant. Moreover, the PCA results showed that Dim1 differentiates two groups of fonio accessions: one group composed essentially of accessions from Mali and Niger, while the second group was composed of accessions from Benin, Burkina, and Guinea. Likewise, Dimension 2 differentiates two groups of accessions, the first group was composed of accessions from Guinea, Benin, and Burkina Faso, while the second group was again a mixed of accessions from these three countries (Figure 5).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Principal components analysis biplot showing variables and individuals on Dimensions 1 and 2.</figcaption></figure>

### 3.6. Partition of Accessions into Phenotypic Groups

The hierarchical clustering on principal components analysis grouped fonio accessions into three different clusters or phenotypic groups (Figure 6). Cluster 1 (50.56% of total accessions), essentially gathered accessions from Mali and Niger, characterized by early flowering and early maturing genotypes (respectively means of 69 ± 3.06 days and 86 ± 3.99 days), and short plants (PHT = 69 ± 5.76 cm), and high number of leaves per plant (NLP = 202 ± 69.74 leaves per plant). Accessions of Cluster 1 produced higher average grain yield (GRY = 1280 ± 660.10 kg.ha⁻¹), with higher harvest index (HI = 50%) and the higher thousand seeds weight (TSW = 0.76 ± 0.07 g). Cluster 1
was predominantly composed of brown color accessions with well exerted panicle. Even though, accessions from this phenotypic group showed high performance in terms of yield, they produced less vigorous plant (CD = 0.77 ± 0.11 mm), and recorded lower values of number of grains per raceme (NGR = 77 ± 10.39 grains), with intermediate raceme length (8 ± 0.79 cm) and panicle length (14 ± 1.63 cm) (Table 5). Cluster 2 (23.89% of total accessions) and Cluster 3 (25.55% of total accessions) were composed of late flowering (FLO = 84 ± 5.21 days for Cluster 2 and 86 ± 3.38 days for Cluster 3) and late maturing accessions (MAT = 104 ± 5.50 days for Cluster 2 and 105 ± 3.38 days for Cluster 3). These phenotypic groups (Cluster 2 and 3) were predominantly composed of accessions with long to intermediate plant height (PHT = 79 ± 5.42 cm for Cluster 2 and PHT = 72 ± 5.33 cm for Cluster 3), with more vigorous plants (CD = 1.14 ± 0.38 mm for Cluster 3 and CD = 0.87 ± 0.09 mm for Cluster 2), and with long racemes (RLT = 9.90 ± 0.61 cm for Cluster 2 and RLT = 9.40 ± 0.58 cm for Cluster 3) and long panicles (PLT = 17 ± 3.54 cm for Cluster 2, PLT = 14 ± 1.05 cm for Cluster 3). Accessions of these two Cluster recorded intermediate number of leaves per plant and intermediate grain yield (GRY = 2048 ± 450.40 kg.ha⁻¹ for Cluster 3 and GRY = 2019 ± 550.80 kg.ha⁻¹ for Cluster 2). However, accessions from Cluster 3 recorded higher number of grains per raceme (NGR = 104 ± 14.27 grains) than accessions from Cluster 2 (NGR = 96 ± 13.22 grains). The same trend was observed for number of racemes per panicle (NRP = 4 ± 0.21 racemes for Cluster 3 and NRP = 3 ± 0.23 racemes for Cluster 2).

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Hierarchical clustering on principal components analysis (HCPC) showing number of clusters and individuals within clusters in fonio accessions.</figcaption></figure>
<table><caption>**Table 5.** Description of phenotypic groups in fonio accessions.</caption><thead><tr><th rowspan="2">Variables</th><th>Cluster 1, N = 91</th><th>Cluster 2, N = 43</th><th>Cluster 3, N = 46</th><th rowspan="2">Pr (&gt;F)</th></tr><tr><th>ML = 48 and NG = 43</th><th>BN = 6, BF = 2,<br>GN = 33, ML = 2</th><th>BN = 43 and BF = 3</th></tr><tr><th></th><th>Mean ± SD</th><th>Mean ± SD</th><th>Mean ± SD</th><th></th></tr></thead><tbody><tr><td>Plant height (cm)</td><td>69 ± 5.76</td><td>79 ± 5.42</td><td>72 ± 5.33</td><td>***</td></tr><tr><td>Number of internodes/plant</td><td>7 ± 0.77</td><td>11 ± 1.20</td><td>10 ± 1.02</td><td>***</td></tr><tr><td>Internode length (cm)</td><td>19 ± 7.69</td><td>39 ± 6.16</td><td>37 ± 6.68</td><td>**</td></tr><tr><td>Number of leaves/plant</td><td>202 ± 69.74</td><td>90 ± 23.84</td><td>83 ± 17.75</td><td>***</td></tr><tr><td>Culm diameter (mm)</td><td>0.77 ± 0.11</td><td>0.87 ± 0.09</td><td>1.14 ± 0.38</td><td>***</td></tr><tr><td>Flowering (days)</td><td>69 ± 3.06</td><td>84 ± 5.21</td><td>86 ± 3.38</td><td>***</td></tr><tr><td>Maturity (days)</td><td>86 ± 3.99</td><td>104 ± 5.50</td><td>105 ± 3.38</td><td>**</td></tr><tr><td>Panicle exertion</td><td>1 ± 0.68</td><td>3 ± 1.42</td><td>3 ± 1.48</td><td>***</td></tr><tr><td>Raceme length (cm)</td><td>8.1 ± 0.79</td><td>9.90 ± 0.61</td><td>9.4 ± 0.58</td><td>***</td></tr><tr><td>Panicle length (cm)</td><td>14 ± 1.63</td><td>17 ± 3.54</td><td>14 ± 1.05</td><td>***</td></tr><tr><td>Number of racemes/panicle</td><td>3 ± 0.29</td><td>3 ± 0.23</td><td>4 ± 0.21</td><td>***</td></tr><tr><td>Number of grains/raceme</td><td>77 ± 10.39</td><td>96 ± 13.22</td><td>104 ± 14.27</td><td>**</td></tr><tr><td>Grain yield (kg/ha)</td><td>3840 ± 660.10</td><td>2019 ± 550.80</td><td>2048 ± 450.40</td><td>**</td></tr><tr><td>Harvest index (%)</td><td>57 ± 8.78</td><td>26 ± 7.68</td><td>26 ± 5.43</td><td>**</td></tr><tr><td>Thousand seeds Weight (g)</td><td>0.76 ± 0.07</td><td>0.68 ± 0.06</td><td>0.59 ± 0.04</td><td>ns</td></tr></tbody></table>

SD = standard deviation, BN = Benin, BF = Burkina Faso, GN = Guinea, ML = Mali, NG = Niger. *** = significance at 0.001, ** = significance at 0.01, ns = not significant.

The HCPC analysis was supported by linear discriminant analysis (LDA) (Figure 7). The LDA results showed that linear discriminators LD1 and LD2 clearly separate the three phenotypic groups without admixture, while LD1 and LD3, LD1 and LD4 are fairly good at separating two to three Clusters. The result of multiple factor analysis (MFA) showed a synthetic comparison of the three phenotypic groups. This method groups variables in different classes of morphology, phenology and yield before making comparison between phenotypic groups. Dimension 1 clearly differentiates Cluster 1 from Clusters 2 and 3, while Dimension 2 differentiates Cluster 1 from Cluster 2.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Linear discriminant analysis (LDA) in fonio.</figcaption></figure>
# 4. Discussion

## 4.1. Qualitative Traits Variation in Fonio

Qualitative traits analysis revealed that most of the fonio accessions had well exerted panicle and brown grain color. Panicle exertion is important in cereal research though this was rarely documented in fonio for which a persistent interest for higher yielding cultivars is nowadays expected. Panicle exertion was reported as a crucial morphological marker that could be used for identification of hybrid plants derived from crosses between cultivated rice and African wild rice (Oryza longistaminata A. Chev. & Roehrich) [35]. The effect of genetic and abiotic factors on the panicle exertion in fonio is a topic of interest for future investigations. Seed color is one of the most important traits that guide breeders, farmers, processors, as well as consumers for cultivar preference [36,37]. In fonio, seed color was often associated with other traits including crop cycle, yield, tolerance to biotic and abiotic stresses, harvesting and processing abilities [3], though this needs to be fully documented. The brown seeded accessions, commonly known as white fonio (reference to *D. exilis*), is the most diverse and widely cultivated cultivar in the region [38]. Except in Niger where the majority of fonio accessions had greyed orange color [21], the predominance of brown seeded color in West Africa could be explained by its well adaptability to various pedo-climatic conditions and its association with other desirable agronomic traits. Lule et al. [39] reported that farmer's preferences of grain color (purple black) in finger millet (Eleusine coracana (L.) Gaertn.) was associated with adaptability to environment factors, yielding ability, major use of the grain, or cultural role and ease of production in northern Ethiopia. Therefore, future research should investigate the genetic basis of the seed color and its association of functional agronomic traits.

## 4.2. Performance of Fonio Accessions Based on Provenance

Our findings revealed that grain yield in fonio from various countries and ecological zones were higher than that reported previously [18]. Similar trend was observed for number of internode and internode length which were also higher compared to results obtained by Saidou et al. on Niger accessions only [21]. The mean thousand seeds weight recorded by Nyam [20] was also lower than that obtained in this study. However, higher values for plant height, number of internode, number of raceme per panicle, raceme length, and panicle length were reported by Sekloka et al. [19] in the characterization at the University of Parakou (Rep. of Benin) of twenty fonio accessions collected from Boukoumbe. This is an illustration of the morphological diversity among fonio accessions which could be due to the expression of genes controlling complex traits, or environmental factors including agricultural practices [10,18,40]. The significant differences observed among fonio accessions for most of the quantitative traits could also be attributed to the variability in rainfall patterns because the cumulative rainfall recorded in 2019 (1211.20 mm) was much higher than that of 2018 (1047.62 mm). This suggests that the amount and distribution of rainfall impact the performance of fonio accessions during the two consecutive years of experiment. Dehgahi et al. [41] reported also that grain yield was positively correlated with the average annual rainfall in wheat. That variability of rainfall can trigger the disruption in crop production and food insecurity [42].

Differences observed among accessions of various provenances indicated that the variation of most quantitative traits may be likely due to adaptation to specific environment and sowing date. Differences for flowering and maturity times were conspicuous among accessions. The early maturing accessions were more frequent in relatively hot and dry areas of arid zone from Mali and Niger whereas the late maturing accessions occurred abundantly in semi-arid and sub-humid areas from Benin, Burkina Faso and Guinea. This suggests that fonio accessions from the arid zone might have lost their photoperiod sensitivity in trying to adapt to prolonged dry spells and intensified heat waves compare to those from sub-humid zone which might be photoperiod sensitive with less drought spells and heat waves. Earlier works reported highly significant differences for number of days to flowering as affected by photoperiods within *D. exilis* [43]. In pearl millet (*Pennisetum glaucum* (L.) R. Br.), one of major crops
cultivated in the dryland areas of West Africa, similar results were reported by Haussmann et al. [44] that long-cycle millets showed a stronger sensitivity to photoperiod than shorter-cycle cultivars in West and Central African pearl millet landraces.

In this study, fonio accessions were sown around mid-July with the consequence of exposing them to flower late in long days. Ordinary, farmers from semi-arid and sub humid zones sow fonio earlier in May. This shift in sowing date might have led to the lower grain yield observed in fonio accessions from sub-humid zone compared to those from arid zones. Gueye et al. [45] also reported that the delay in sowing dates reduced grain yield by 50% to 87% in fonio. Our result was in accordance with those of Wolabu and Tadege [46], who reported that regardless of the day length, non-photoperiod sorghum genotypes flower earlier than photoperiod sensitive genotypes. Future research will certainly decipher the effect of specific environment on the yield and yield components through a multiple environmental trial in order to reveal the ability of each ecological zone to contribute to yield accumulation in Fonio. This will clearly highlight the vulnerability and adaptability of fonio along the ecological gradient of west Africa as this was investigated for wheat in Pakistan [47].

### 4.3. *Quantitative Traits Association and Partition of Phenotypic Groups in Fonio*

Correlation analysis revealed highly significant positive relationships between grain yield and harvest index, thousand seed weight, number of leaves. Traits with high significant positive correlations would contribute to grain yield improvement in fonio. This implies that those traits can be used as proxy for yield in indirect selection and harvest index could be recommended as a selection criterion for increasing yield in fonio. In finger millet, similar results were reported that grain yield was positively associated with harvest index [48]. However, number of internode and plant height negatively impacted grain yield. Likewise, long flowering and maturing times negatively influenced grain yield. Similar results were found on wheat (*Triticum aestivum*) [49]. However, some findings showed that grain yield had highly significant positive correlation with maturity and flowering times in finger millet [50]. In addition, grain yield exhibited moderate negative correlation with plant height indicating that higher grain yield is associated with shorter fonio plants. This is an important feature in tackling lodging problem which is the main cause of yield losses and quality reduction in fonio production in West Africa. So, breeding for shorter plant could enhance lodging resistance in fonio varieties. Similar results were reported in tef (*Eragrotis tef* (Zucc.) Trotter) [51].

Yield, yield components, and phenological traits were the contributors to the variation observed in this study. Clusters analysis and linear discriminant analysis (LDA) showed three clusters of fonio accessions. A previous report indicated four clusters while characterizing 20 fonio accessions only from Benin [19]. Another report revealed two groups in the evaluation of 30 fonio accessions from Nigeria [20]. Our findings cover accessions from five countries (Benin, Burkina Faso, Guinea, Mali, and Togo) and three phytogeographical areas. However, accessions from Nigeria should be added for a broader perspective on the morphological variation in fonio in relation to different agro ecological zones in order to enrich the cultivated germplasm from West Africa. With the clear separation of phenotypic groups, it becomes critical to revisit the level of genetic diversity within fonio germplasm [18] and the association between genetic and phenotypic features for proper improvement of the crop. We found that fonio accessions clustered in phenotypic groups based on their provenance; though accessions from Benin and Mali that were distributed throughout Cluster 1 and Cluster 2. Amplified fragment length polymorphisms (AFLPs) analysis of 122 fonio accessions [18] revealed similar trends where fonio accessions from Benin, Burkina Faso, Guinea, and Togo clustered in identical group based on their origin except accessions from Guinea and Mali. In addition to site specificity, genes flow is a powerful mechanism that contribute to homogenization of variability distribution. For instance, grouping of accessions regardless of boundaries reflect the importance of farmers’ seed exchange in fonio seed systems. As informal seed systems, seed exchange is commonly practiced by fonio farmers as reported by previous works [2,10,52].
Group 1 was composed essentially of early flowering and maturing accessions, and recorded higher values for grain yield, harvest index and thousand seeds weight. The accessions of this group are adequate candidates for developing early maturing fonio lines for arid and semi-arid regions where drought stress occurs frequently. Furthermore, this group is composed of fonio accessions with shorter plants that could be used as source for lodging resistance since short stature plants have been reported in cereals to be preferred in varietal development for that trait. Phenotypic Groups 2 and 3 were composed predominantly of late maturing accessions with long to intermediate plant height, with higher values of culm diameter. Moreover, these groups were composed of accessions with long racemes and long panicles. The difference observed in the groups could be likely due to the genetic make-up of each genotype and prevailing environmental factors. However, despite their phenotypical similitude, accessions from Group 3 recorded higher number of grains per raceme and number raceme per panicle than accessions from Group 2 revealing another agronomic trait of importance for the improvement of the crop. Marker assisted index selection of fonio lines with long raceme, and long panicle alongside with high number of grains per raceme and high number of racemes per panicle can be explored as crop a improvement avenue. Moreover, genome-wide association studies (GWAS) are needed to understand the genetic basis of key important agronomic traits in fonio.

**5. Conclusions**

This study assessed phenotypic variability among 180 fonio accessions from five West African countries. It revealed the existence of wide range of agro morphological variability. The results revealed interesting correlations between yield and yield components. Furthermore, clusters analysis revealed three phenotypic groups. Group 1 was composed of early flowering and maturing accessions with higher grain yield, higher harvest index, higher thousand seeds weight, and higher number of leaves per plant. Accessions from Group 2 and 3 were characterized by late flowering and maturing times, with intermediate grain yield, thousand seeds weight, and number of leaves per plant, with taller plants, higher values of number of grains per raceme, and number of racemes per panicle. Accessions from Group 1 are candidates for yield improvement and enhancing lodging resistance, whereas accessions from Group 2 and Group 3 are candidate also for yield improvement via indirect selection and good parent for hybrids development. The study exhibited agro morphological descriptors that discriminate fonio accessions and provide also useful information for parental selection with key economically important agronomic traits.

**Author Contributions:** A.R.I.B.Y. and E.G.A.-D. conceptualized the study. A.R.I.B.Y. and M.A. conducted the experiments. A.R.I.B.Y. and E.G.A.-D. wrote the manuscript. A.R.I.B.Y. and E.G.A.-D. edited the manuscript. E.G.A.-D., M.A., E.S., C.O.A.A., C.B., A.B., and Y.B. critically revised the draft and updated the manuscript for publication. All authors have read and agreed to the published version of the manuscript.

**Funding:** Ibrahim Bio Yerima, from the National Institute of Agronomic Research of Niger (INRAN), is a scholar of the "Intra-Africa Academic Mobility Scheme" under the project grant number 2016-2988 on "Enhancing training and research mobility for novel crops breeding in Africa (MoBreed)" funded by the Education, Audiovisual and Culture Executive Agency (EACEA) of the European Commission. The project provided a scholarship for academic training and research mobility and a research grant to the first author to complete a PhD degree at the University of Abomey-Calavi (Republic of Benin).

**Acknowledgments:** The authors wish to thank MoBreed program for the financial support. We are thankful to all farmers and agriculture extension services in Benin, Burkina Faso, Guinea, Mali, and Niger. We also thank Dèdéou Tchopkonhoue, Djamilou Ibrahim, Doumbia Zan Ibrahima, Diarra Cheik, and Habou Kano for assistance during the survey, field activities and data analysis.

**Conflicts of Interest:** Authors declare no conflict of interest.

**References**

1. Garí, J.A. Review of the African millet diversity. In *International Workshop on Fonio, Food Security and Livelihood among the Rural Poor in West Africa*; IPGRI/IFAD: Bamako, Mali, 2002.
2. Adoukonou-Sagbadja, H.; Dansi, A.; Vodouhè, R.; Akpagana, K. Indigenous knowledge and traditional conservation of fonio millet (*Digitaria exilis*, *Digitaria iburua*) in Togo. *Biodivers. Conserv.* **2006**, *15*, 2379–2395. [CrossRef]

3. Ballogou, V.Y.; Soumanou Mohamed, M.; Toukourou, F.; Hounhouigan Joseph, D. Indigenous knowledge on landraces and fonio-based food in Benin. *Ecol. Food Nutr.* **2014**, *53*, 390–409. [CrossRef]

4. Sartelet, H.; Serghat, S.; Lobstein, A.; Ingenbleek, Y.; Anton, R.; Petitfrere, E.; Aguie-Aguie, G.; Martiny, L.; Haye, B. Flavonoids extracted from fonio millet (*Digitaria exilis*) reveal potent antithyroid properties. *Nutrition* **1996**, *12*, 100–106. [CrossRef]

5. Vodouhè, S.; Achigan Dako, E.G. *Digitaria exilis* (Kippist) Stapf; Plant Resources of Tropical Africa 1. Cereals and Pulses; Brink, M., Belay, G., Eds.; PROTA Foundations: Wageningen, The Netherlands; Backhuys Publishers: Leiden, The Netherlands, 2006; Volume 1.

6. Fogny, N.F.; Madode, E.Y.; Laleye, F.F.; Amoussou-Lokossou, Y.; Kayode, A.P. Formulation de farine de fonio enrichie en ressources alimentaires locales pour l'alimentation complémentaire des jeunes enfants au Bénin. *Int. J. Biol. Chem. Sci.* **2017**, *11*, 2745–2755. [CrossRef]

7. Cruz, J.; Beavogui, F.; Dramé, D. *Fonio, an African Cereal* [*Le Fonio, une Céréale Africaine*]; CTA: Wageningen, The Netherlands, 2011.

8. Sekloka, E.; Adoukonou-Sagbadja, H.; Paraïso Armand, A.; Yoa Brigitte, K.; Bachabi, F.-X.; Zoumarou-Wallis, N. Evolution de la diversité des cultivars de fonio pratiqués dans la commune de Boukoumbé au Nord-Ouest du Bénin. *Int. J. Biol. Chem. Sci.* **2015**, *9*, 2446–2458. [CrossRef]

9. Ayenan, M.A.T.; Sodedji, K.A.F.; Nwankwo, C.I.; Olodo, K.F.; Alladassi, M.E.B. Harnessing genetic resources and progress in plant genomics for fonio (*Digitaria* spp.) improvement. *Genet. Resour. Crop Evol.* **2017**, *65*, 373–386. [CrossRef]

10. Dansi, A.; Adoukonou-Sagbadja, H.; Vodouhè, R. Diversity, conservation and related wild species of Fonio millet (*Digitaria* spp.) in the northwest of Benin. *Genet. Resour. Crop Evol.* **2010**, *57*, 827–839. [CrossRef]

11. Bhanu, A.N. Assessment of Genetic Diversity in Crop Plants—An Overview. *Adv. Plants Agric. Res.* **2017**, *7*, 279–286. [CrossRef]

12. Mondini, L.; Noorani, A.; Pagnotta, M. Assessing Plant Genetic Diversity by Molecular Tools. *Diversity* **2009**, *1*, 19–35. [CrossRef]

13. Odjo, T.C.; Dossou-Aminon, I.; Dansi, A.; Bonou-Gbo, Z.; Kombaté, K. Agro-Morphological Characterization and Assessment of Variability within a Germplasm of Benin Rice (*Oryza sativa* L.) Varieties. *Int. J. Curr. Res. Biosci. Plant Biol.* **2017**, *4*, 1–16.

14. Babu, B.K.; Sood, S.; Agrawal, P.; Chandrashekara, C.; Kumar, A.; Kumar, A. Molecular and Phenotypic Characterization of 149 Finger Millet Accessions Using Microsatellite and Agro-Morphological Markers. *Proc. Natl. Acad. Sci. India Sect. B Biol. Sci.* **2016**, *87*, 1217–1228. [CrossRef]

15. Sattler, F.; Sanogo, M.; Kassari, I.; Angarawai, I.; Gwadi, K.; Dodo, H.; Haussmann, B. Characterization of West and Central African accessions from a pearl millet reference collection for agro-morphological traits and Striga resistance. *Plant Genet. Resour. Charact. Util.* **2017**, *16*, 260–272. [CrossRef]

16. Naoura, G.; Sawadogo, N.; Atchozou, E.A.; Emendack, Y.; Hassan, M.A.; Reoungal, D.; Amos, D.N.; Djirabaye, N.; Tabo, R.; Laza, H. Assessment of agro-morphological variability of dry-season sorghum cultivars in Chad as novel sources of drought tolerance. *Sci. Rep.* **2019**, *9*, 19581. [CrossRef]

17. Clottey, V.; Agyare, W.; Bayorbor, T.; Abanga, J.; Kombiok, J. Genetic relatedness of fonio (*Digitaria* spp.) landraces assembled in Ghana. *Plant Genet. Resour. Newsl.* (*Bioversity Int. FAO*) **2006**, *147*, 6–11.

18. Adoukonou-Sagbadja, H.; Wagner, C.; Dansi, A.; Ahlemeyer, J.; Dainou, O.; Akpagana, K.; Ordon, F.; Friedt, W. Genetic diversity and population differentiation of traditional fonio millet (*Digitaria* spp.) landraces from different agro-ecological zones of West Africa. *Theor. Appl. Genet.* **2007**, *115*, 917–931. [CrossRef]

19. Sekloka, E.; Kanlindogbe, C.; Biaou Samadori, S.H.; Adoukonou-Sagbadja, H.; Kora, A.; Motouama, F.T.; Seidou, M.; Zinsou Valérien, A.; Afouda, L.; Baba-Moussa, L. Agro-morphological characterization of fonio millet accessions (*Digitaria exilis* Stapf.) collected from Boukoumb, Northwest of Benin. *J. Plant Breed. Crop Sci.* **2016**, *8*, 211–222.

20. Nyam, D.; Kwon-Ndung, E.; Wuyep, A. Genetic affinity and breeding potential of phenologic traits of acha (fonio) in Nigeria. *J. Sci. Eng. Res.* **2017**, *4*, 91–101.
21. Saidou, S.I.; Bakasso, Y.; Inoussa, M.M.; Zaman-Allah, M.; Atta, S.; Barnaud, A.; Billot, C.; Saadou, M. Diversité agro-morphologique des accessions de fonio [*Digitaria exilis* (Kippist.) Stapf.] au Niger. *Int. J. Biol. Chem. Sci.* **2014**, *8*, 1710. [CrossRef]

22. Peyre, D. *Catalogue des Plantes Vasculaires du Niger*; Institut d'Elevage et de Médecine Vétérinaire des Pays Tropicaux: Paris, France, 1976; Volume 1, pp. 405-419.

23. Peyre de Fabregues, B.; Lebrun, J. *Catalogue des Plantes Vasculaires du Niger (Etude Botanique n 3.)* IEMVT; Laboratoire National d'Elevage: Maisons-Alfort, France; Niamey, Niger, 1976; p. 433.

24. White, F. *The Vegetation of Africa*; Natural Resources Research; UNESCO: Paris, France, 1983; Volume 20, p. 356.

25. Dobignard, A.; Chatelain, C. *Index Synonymique de la Flore d'Afrique du Nord*; Éditions des Conservatoire et Jardin Botaniques de la Ville de Genève: Geneve, Switzerland, 2010; Volume 1, p. 455.

26. Prance, G.; Sothers, C. *Chrysobalanaceae 1 & 2, Species Plantarum: Flora of the World, Parts 9 &10*; Australian Biological Resources Study: Canberra, Australia, 2003.

27. Bioversity International, I.W. Descriptors for Wild and Cultivated Rice (Oryza spp.); Bioversity International: Rome, Italy, 2007.

28. Kassambara, A. *Practical Guide to Principal Component Methods in R (Multivariate Analysis)*; CreateSpace Independent Publishing Platform: Scotts Valley, CA, USA, 2017; Volume 2.

29. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. *J. Stat. Softw.* **2008**, *25*, 1-18.
[CrossRef]

30. Kassambara, A.; Mundt, F. Package 'Factoextra'. Extract and Visualize the Results of Multivariate Data Analyses; R package. 2017. Available online: http://www.sthda.com/english/rpkgs/factoextra (accessed on 5 August 2020).

31. Team, R.C. Version 3.5. 1. R Foundation for Statistical Computing, Vienna, Austria. Available online: https://www.r-project.org (accessed on 3 November 2018).

32. Argüelles, M.; Benavides, C.; Fernández, I. A new approach to the identification of regional clusters: Hierarchical clustering on principal components. *Appl. Econ.* **2014**, *46*, 2511-2519. [CrossRef]

33. Rencher, A. Discriminant analysis: Description of group separation. In *Methods of Multivariate Analysis*, 2nd ed.; John Wiles and Sons: New York, NY, USA, 2002; pp. 270-298.

34. Pagès, J. *Multiple Factor Analysis by Example Using R*; CRC Press: Boca Raton, FL, USA, 2014.

35. Kanya, J.I.; Hauser, T.P.; Kinyamario, J.I.; Amugune, N.O. Hybridization potential between cultivated rice *Oryza sativa* and African wild rice *Oryza longistaminata*. *Int. J. Agric. Res.* **2012**, *7*, 291-302. [CrossRef]

36. Govindaraj, M.; Rao, A.S.; Shivade, H.; Rai, K. Effect of grain colour on iron and zinc density in pearl millet.
*Indian J. Genet. Plant Breed.* **2018**, *78*, 247. [CrossRef]

37. Garg, M.; Chawla, M.; Chunduri, V.; Kumar, R.; Sharma, S.; Sharma, N.K.; Kaur, N.; Kumar, A.; Mundey, J.K.; Saini, M.K.; et al. Transfer of grain colors to elite wheat cultivars and their characterization. *J. Cereal Sci.* **2016**, *71*, 138-144. [CrossRef]

38. Blench, R.M. Vernacular names for African millets and other minor cereals and their significance for agricultural history. *Archaeol. Anthropol. Sci.* **2016**, *8*, 1-8. [CrossRef]

39. Lule, D.; Tesfaye, K.; Fetene, M. Qualitative traits diversity and eco-geographical distribution in finger millet (*Eleusine coracana*) landraces from eastern and south eastern Africa: An implication for germplasm collection and conservation. *Afr. J. Plant Sci.* **2012**, *6*, 346-354.

40. Vigouroux, Y.; Barnaud, A.; Scarcelli, N.; Thuillet, A.-C. Biodiversity, evolution and adaptation of cultivated crops. *Comptes Rendus Biol.* **2011**, *334*, 450-457. [CrossRef]

41. Dehgahi, R.; Joniyas, A.; Latip, S.N.H.B.M. Rainfall Distribution and Temperature Effects on Wheat Yield in Torbate Heydarie. *Int. J. Sci. Res. Knowl. ISSN 2014*, 2, 2322-4541.

42. Kyei-Mensah, C.; Kyerematen, R.; Adu-Acheampong, S. Impact of Rainfall Variability on Crop Production within the Worobong Ecological Area of Fanteakwa District, Ghana. *Adv. Agric.* **2019**, *2019*. [CrossRef]

43. Aliero, A.; Morakinyo, J. Photoperiodism in *Digitaria exilis* (Kipp) Stapf accessions. *Afr. J. Biotechnol.* **2005**, *4*, 241-243.

44. Haussmann, B.I.; Fred Rattunde, H.; Weltzien-Rattunde, E.; Traoré, P.; Vom Brocke, K.; Parzies, H.K. Breeding strategies for adaptation of pearl millet and sorghum to climate variability and change in West Africa. *J. Agron. Crop Sci.* **2012**, *198*, 327-339. [CrossRef]
45. Gueye, M.; Kanfany, G.; Fofana, A.; Noba, K.; Grove, J. Effect of planting date on growth and grain yield of fonio millet (*Digitaria exilis Stapf*) in the Southeast of Senegal. *Int. J. Biol. Chem. Sci.* **2015**, *9*, 581–592.
[CrossRef]

46. Wolabu, T.W.; Tadege, M. Photoperiod response and floral transition in sorghum. *Plant Signal. Behave.* **2016**, *11*, e1261232. [CrossRef]

47. Sultana, H.; Ali, N.; Iqbal, M.M.; Khan, A.M. Vulnerability and adaptability of wheat production in different climatic zones of Pakistan under climate change scenarios. *Clim. Chang.* **2009**, *94*, 123–142. [CrossRef]

48. Jyothsna, S.; Patro, T.; Ashok, S.; Rani, Y.S.; Neeraja, B. Studies on genetic parameters, character association and path analysis of yield and its components in finger millet (*Eluesine coracana* L. Gaertn). *Int. J. Theor. Appl. Sci.* **2016**, *8*, 25.

49. Assefa, E. Correlation and Path Coefficient Studies of Yield and Yield Associated Traits in Bread Wheat (*Triticum aestivum* L.) Genotypes. *Adv. Plants Agric. Res.* **2017**, *6*, 1–10. [CrossRef]

50. Chavan, B.; Jawale, L.N.; Shinde, A.V. Correlation and path analysis studies in finger millet for yield and yield contributing traits (*Eleusine coracana* L. Gaertn). *Int. J. Chem. Stud.* **2020**, *8*, 2911–2914.

51. Jifar, H.; Dagne, K.; Tesfaye, K.; Assefa, K.; Tadele, Z. Agro-Morphological Traits Diversity in Tef [Eragrostis Tef (Zucc.) Trotter] Genotypes from Various Sources. *Ethiop. J. Agric. Sci.* **2018**, *28*, 131–148.

52. Sidibé, A.; Meldrum, G.; Coulibaly, H.; Padulosi, S.; Traore, I.; Diawara, G.; Sangaré, A.R.; Mbosso, C. Revitalizing cultivation and strengthening the seed systems of fonio and Bambara groundnut in Mali through a community biodiversity management approach. *Plant Genet. Resour.* **2020**, *19*, 31–48.

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).