# Genotype-Environment Interaction: Trade-Offs between the Agronomic Performance and Stability of Dual-Purpose Sorghum (Sorghum bicolor L. Moench) Genotypes in Senegal

Malick Ndiaye, Myriam Adam, Komla Kyky Ganyo, Aliou Guissé, Ndiaga Cissé, Bertrand Muller

► To cite this version:

Malick Ndiaye, Myriam Adam, Komla Kyky Ganyo, Aliou Guissé, Ndiaga Cissé, et al. Genotype-Environment Interaction: Trade-Offs between the Agronomic Performance and Stability of Dual-Purpose Sorghum (Sorghum bicolor L. Moench) Genotypes in Senegal. Agronomy, 2019, 9 (12), pp.1-16. 10.3390/agronomy9120867. [https://doi.org/10.3390/agronomy9120858](https://doi.org/10.3390/agronomy9120858)

HAL Id: hal-02620858

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

Submitted on 26 May 2020

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
Article

# Genotype-Environment Interaction: Trade-Offs between the Agronomic Performance and Stability of Dual-Purpose Sorghum (*Sorghum bicolor* L. Moench) Genotypes in Senegal

Malick Ndiaye ¹,²,*, Myriam Adam ³,⁴,⁵,⁶,*, Komla Kyky Ganyo ⁷, Aliou Guissé ⁸,⁹, Ndiaga Cissé ² and Bertrand Muller ²,³,¹⁰

¹ Institut Sénégalais de Recherches Agricoles (ISRA), Centre National de Recherches Agronomiques (CNRA) de Bambey, Bambey BP 53, Senegal

² Centre d'Etude Régional pour l'Amélioration de l'Adaptation à la Sécheresse (CERAAS), Route de Khombole, Thiès BP 3320, Senegal; ncisse@refer.sn (N.C.); bertrand.muller@cirad.fr (B.M.)

³ AGAP, Univ de Montpellier MUSE, CIRAD, INRA, 34090 Montpellier SupAgro, France

⁴ International Center for Agricultural Research for Development (CIRAD), UMR AGAP, Bobo-Dioulasso, Burkina Faso

5 Institut de l'Environnement et de la Recherche Agricole (INERA), Bobo Dioulasso 01 BP 910, Burkina Faso

⁶ International Crops Research Institute for Semi-Arid Tropics (ICRISAT), Bamako BP 320, Mali

⁷ Institut Togolais de Recherche Agronomique (ITRA), Lomé B.P. 1163, Togo; desireganyo@gmail.com

8 Département de Biologie Végétale, Université Cheikh Anta Diop de Dakar (UCAD), Dakar Fann BP 5005, Senegal; alguisse@orange.sn

9 Unité Mixte de Recherche Internationale Environnement, Santé, Sociétés (UMI 3189 ESS), Université Cheikh Anta Diop de Dakar (UCAD), Dakar Fann BP 5005, Senegal

¹⁰ International center for agricultural research for development (CIRAD), Ampandrianomby BP853, Antananarivo 101, Madagascar

* Correspondence: agromalick@yahoo.fr (M.N.); myriam.adam@cirad.fr (M.A.); Tel.: +221 77557 49 61 (M.N.)

Received: 27 August 2019; Accepted: 21 November 2019; Published: 10 December 2019

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

check for
updates

**Abstract:** Introducing sorghum (*Sorghum bicolor* L. Moench) genotypes into new environments is necessary for expanding the production of food and fuel, but these efforts are complicated by significant genotype × environment interactions that can reduce their effectiveness. This study set out to thoroughly analyze genotype × environment interactions and assess trade-offs between the agronomic performance and the stability of grain and biomass yields of ten contrasting genotypes under Sudano-Sahelian conditions. Experiments were carried out in a randomized complete block design with four replicates. They were conducted from 2013 to 2016 in Bambey, Sinthiou Malem and Nioro du Rip in Senegal. The joint analysis of variance revealed a highly significant effect (*p* < 0.0001) of genotypes (G), environments (E) and G × E interaction. Most genotypes showed specific adaptations. The best grain yields were obtained by the Nieleni and Fadda hybrids, while the improved varieties IS15401 and SK5912 were best for biomass production. An Additive Main effect and Multiplicative Interaction (AMMI) analysis showed that good grain yields were associated with environments having good soil fertility and good rainfall, while biomass yields were more influenced by the sowing date and rainfall. Similarly, we were able to confirm for our 10 sorghum genotypes that yield stability was generally associated with low performance, except for the Nieleni and Fadda hybrids, which performed well for grain and biomass production regardless of the environment. The Senegalese control genotype, 621B, showed particular susceptibility to growing conditions (soil), but remained very productive (more than 3 tons per hectare) under good agro-pedological conditions. These results lead us to recommend the Fadda and Nieleni hybrids for the entire study region, while 621B can also be recommended, but only for highly specific environments with good soils.
**Keywords:** G × E interaction; performance; stability; AMMI; hybrid genotype

# 1. Introduction

Sorghum (*Sorghum bicolor* (L.) Moench) is one of the main cereals grown in arid and semi-arid tropical regions [1]. Sorghum is well-adapted to warm regions and, given its plasticity, is able to grow in both temperate and tropical regions. With a global production of about 68.9 million tons in 2015, from around 49.9 million hectares, sorghum ranks fifth in cereal production after maize, wheat, rice and barley [2]. It is mainly used for animal feed in most developed countries, but in Africa and India it is a staple food for millions of people [3]. In addition, sorghum is one of the most important crops that can be used for bioethanol production [4]. In Senegal, after pearl millet and maize, sorghum is the third most important dryland cereal crop, with an estimated total area of more than 221,329 ha for a national production of 225,865 tons and a mean yield of 1,020 kg ha⁻¹ [5]. Sorghum production is essential for subsistence agriculture [6]. However, its production comes up against several constraints that lead to low yields, such as irregularities in rainfall distribution exacerbated by climate change, low soil fertility and sandy soils, and various crop diseases and pests [7].

Food security initiatives in Senegal include introducing new sorghum genotypes adapted to different soil and climate environments. However, when genotypes are evaluated for recommendation, a common problem arises: the high variability of their productivity from year to year and from environment to environment. Such variability creates difficulty in determining which genotypes can be recommended, so it deserves careful consideration. The different responses of a genotype in different environments are known as genotype × environment interaction (G × E). Understanding G × E interaction will help to (1) identify genotypes with a stable performance in fairly diverse growing conditions, and (2) match specific genotypes to specific environments [8].

Several statistical methods have been developed to characterize the effect of G × E interactions of genotypes and to predict phenotypic responses to environmental changes. However, statistical methods for characterizing stability are generally not able to provide an accurate and complete response model for this interaction [9], as the genotypic response to environmental variation is multivariate, while most stability indices have a univariate response [10]. Other methods have therefore been developed to explore G × E interaction models. Of these, the AMMI is a robust multivariate method for multi-environmental trials [11]. The additive main effect and multiplicative interaction (AMMI) method combines an analysis of variance (ANOVA) and a principal component analysis (PCA) in a unified approach that can be used to analyze multi-location trials [12–14]. The ANOVA studies the main effects of genotypes and environments and the principal component analysis (PCA) then focuses on the non-additive part of the model representing interaction (G × E). AMMI provides the G × E interaction sum of squares with a minimum number of degrees of freedom. In addition, AMMI concurrently quantifies the contribution of each genotype and environment to G × E interaction, and provides an easy graphical interpretation of the results using a biplot technique to classify genotypes and environments together [12,15]. This technique can therefore be used to identify productive genotypes with wide adaptability and mega-environments, and to delimit environments in which genotypes have specific adaptability [14–16].

The objective of this study was to: (1) analyze the genotype × environment interactions, adaptability and stability of 10 sorghum genotypes in several environments in Senegal using the AMMI method, and (2) identify genotypes that performed well in terms of grain and/or biomass yield (i.e., dual-purpose: food and feed).
## 2. Materials and Methods

### 2.1. Study Sites

The experiments were conducted during the rainy seasons in 2013, 2014, 2015 and 2016 at three locations in Senegal: the research stations in Sinthiou Salem (in 2013, 2014, 2015 and 2016), Bambey (in 2013) and Nioro du Rip (in 2015). The characteristics of the different sites are given in Table 1. Figure 1 shows the rainfall and temperature in the three locations over the trial period. Bambey is subject to a typical Sahelian climate characterized by a long 8 to 9-month dry season and a 3 to 4-month rainy season. Rainfall varies greatly from one year to another. The dominant soils are sandy with a very low water retention capacity of 80 to 100 mm m⁻¹ [17–19]. The Nioro du Rip and Sinthiou Salem stations are located at the interface between the Sahel and Sudanese zones. They benefit from a 4 to 5-month rainy season that is wetter than in Bambey. They are also characterized by strong inter-annual variability. The soils remain predominantly sandy but have slightly higher clay and silt contents, with a water retention capacity ranging from 90 to 120 mm.m⁻¹ [20–22].

### 2.2. Plant Material

The plant material consisted of ten genotypes from various regions of West and Central Africa, each of which was known to perform well in its area of origin. They were selected to make up a contrasting sample in terms of crop cycle duration (each adapted to its target region), morphology (height, stem diameter in particular), structural characteristics (lignin, cellulose), and grain and biomass production. The characteristics of these ten genotypes are presented in Table 2.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Rainfall and minimum and maximum temperatures at Bambey 2013 (A), Nioro du Rip 2015 (B) and Sinthiou Salem 2013 (C), 2014 (D), 2015 (E) and 2016 (F). Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December.</figcaption></figure>
<table><caption>Table 1. Main characteristics of the different trial sites.</caption><thead><tr><th>Environment</th><th>Zone</th><th>Code</th><th>Coordinates</th><th>Alt (m)</th><th>Soil Type *</th><th>SAN (%)</th><th>CS (%)</th><th>N (%)</th><th>OM (%)</th><th>Rain (mm)</th><th>R<sub>0-30</sub> (mm)</th><th>R<sub>30-60</sub> (mm)</th><th>R<sub>60-90</sub> (mm)</th><th>R<sub>90-120</sub> (mm)</th><th>Tmin (°C)</th><th>Tmax (°C)</th><th>Healthy **</th><th>Sowing Date</th><th>Previo-uscrop</th></tr></thead><tbody><tr><td>Sowing 1/2013</td><td>BBY</td><td>B13D1</td><td>14°42′N</td><td rowspan="2">20</td><td>Sandy</td><td>94.2</td><td>6.6</td><td>0.15</td><td>3.1</td><td>644</td><td>180</td><td>352</td><td>110</td><td>3</td><td>23</td><td>33.9</td><td>2</td><td>07/17/2013</td><td>Fallow</td><td></td></tr><tr><td>Sowing 2/2013</td><td>BBY</td><td>B13D2</td><td>16°29′W</td><td>Sandy</td><td>94.2</td><td>6.6</td><td>0.15</td><td>3.1</td><td>566</td><td>253</td><td>256</td><td>56</td><td>1</td><td>22.8</td><td>33.9</td><td>3</td><td>07/31/2013</td><td>Fallow</td><td></td></tr><tr><td>Sowing 1/2013</td><td>SIN</td><td>S13D1</td><td></td><td rowspan="6">23</td><td>Sandy-silty</td><td>89.4</td><td>11.6</td><td>0.21</td><td>4.3</td><td>575</td><td>146</td><td>365</td><td>59</td><td>6</td><td>21.4</td><td>35.3</td><td>5</td><td>07/25/2013</td><td>Fallow</td><td></td></tr><tr><td>Sowing 2/2013</td><td>SIN</td><td>S13D2</td><td></td><td>Sandy-silty</td><td>89.4</td><td>11.6</td><td>0.21</td><td>4.3</td><td>536</td><td>183</td><td>306</td><td>46</td><td>1</td><td>21.2</td><td>35.4</td><td>5</td><td>08/06/2013</td><td>Fallow</td><td></td></tr><tr><td>Sowing 1/2014</td><td>SIN</td><td>S14D1</td><td>13°49′N</td><td>Sandy</td><td>91.2</td><td>10.2</td><td>0.17</td><td>5.7</td><td>488</td><td>158</td><td>213</td><td>88</td><td>31</td><td>22.2</td><td>35.7</td><td>4</td><td>07/17/2014</td><td>Peanut</td><td></td></tr><tr><td>Sowing 2/2014</td><td>SIN</td><td>S14D2</td><td>13°55′W</td><td>Sandy</td><td>90.9</td><td>9.7</td><td>0.17</td><td>4.5</td><td>377</td><td>190</td><td>156</td><td>31</td><td>1</td><td>22.1</td><td>35.6</td><td>5</td><td>08/06/2014</td><td>Peanut</td><td></td></tr><tr><td>Sowing 1/2015</td><td>SIN</td><td>S15D1</td><td></td><td>Sandy</td><td>93.7</td><td>6.3</td><td>0.32</td><td>3.5</td><td>505</td><td>52</td><td>259</td><td>153</td><td>43</td><td>21.8</td><td>34.7</td><td>2</td><td>07/09/2015</td><td>Peanut</td><td></td></tr><tr><td>Sowing 2/2015</td><td>SIN</td><td>S15D2</td><td></td><td>Sandy</td><td>93.2</td><td>6.8</td><td>0.33</td><td>3.8</td><td>455</td><td>259</td><td>155</td><td>44</td><td>2</td><td>21.2</td><td>34.9</td><td>4</td><td>08/08/2015</td><td>Peanut</td><td></td></tr><tr><td>Sowing 1/2016</td><td>SIN</td><td>S16D1</td><td></td><td></td><td>Sandy-silty</td><td>84.1</td><td>15.9</td><td>0.55</td><td>10.6</td><td>447</td><td>230</td><td>155</td><td>24</td><td>38</td><td>22.5</td><td>35.6</td><td>5</td><td>07/25/2016</td><td>Fallow</td><td></td></tr><tr><td>Sowing 1/2015</td><td>NIO</td><td>N15D1</td><td>13°45′N</td><td rowspan="2">45</td><td>Sandy</td><td>92.4</td><td>7.6</td><td>0.31</td><td>3.5</td><td>943</td><td>196</td><td>361</td><td>261</td><td>126</td><td>20.6</td><td>33.8</td><td>3</td><td>07/16/2015</td><td>Cowpea</td><td></td></tr><tr><td>Sowing 2/2015</td><td>NIO</td><td>N15D2</td><td>15°45′W</td><td>Sandy-silty</td><td>87.0</td><td>13.0</td><td>0.43</td><td>6.1</td><td>747</td><td>329</td><td>273</td><td>145</td><td>0</td><td>19.7</td><td>33.8</td><td>4</td><td>08/13/2015</td><td>Fallow</td><td></td></tr></tbody></table>

BBY = Bambey, SIN = Sinthiou Malem, NIO=Nioro du Rip, Alt = Altitude, SAN = Sand, CS = Clay + Silt. * Classification according to the USDA method based on average data over the 0–30 cm horizon; R<sub>0-30</sub> = total rainfall between 0 and 30 days after sowing, R<sub>30-60</sub> total rainfall between 30 and 60 days after sowing, R<sub>60-90</sub> = total rainfall between 60 and 90 days after sowing, R<sub>90-120</sub> = total rainfall between 90 and 120 days after sowing, ** = score given to a given environment according to disease level: the favorable situation takes the score 5 (absence of disease) and the unfavorable situation the score 1 (strong presence of disease). Rain = total rainfall during the trial.

<table><caption>Table 2. Main characteristics of the ten genotypes studied.</caption><thead><tr><th>Genotype</th><th>Code</th><th>Type</th><th>Photoperiod-Sensitivity</th><th>Cycle Duration</th><th>Isohyets</th><th>Purpose</th><th>Plant Height</th><th>Yield Potential</th><th>Panicle Shape</th><th>Others</th><th>Origin</th></tr></thead><tbody><tr><td>Fadda</td><td>G1</td><td>Guinea (Hybride)</td><td>Moderate</td><td>110 days</td><td>700-1000 mm</td><td>Grain-biomass</td><td>2-3 m</td><td>4.5 t/ha</td><td>Semi-loose</td><td>Tolerant: mold, anthracnose</td><td>Mali, IER/ICRISAT selection, pedigree 02-SB-F5DT-12A xLata.</td></tr><tr><td>Nieleni</td><td>G2</td><td>Caudatum (Hybride)</td><td>Low</td><td>100 days</td><td>700-800 mm</td><td>Grain</td><td>3 m</td><td>4 t/ha</td><td>Semi-compact</td><td>Tolerant: mold, anthracnose</td><td>Mali, IER/ICRISAT selection</td></tr><tr><td>IS15401</td><td>G3</td><td>Guinea</td><td>High</td><td>120 days</td><td>900-1200 mm</td><td>Biomass</td><td>4-4.5 m</td><td>2 t/ha</td><td>Semi-compact</td><td>Resistant: mold, striga and midges</td><td>Cameroon, IER/ICRISAT selection</td></tr><tr><td>Pablo</td><td>G4</td><td>Guinea (Hybride)</td><td>Moderate</td><td>110 days</td><td>700-1000 mm</td><td>Biomass</td><td>4 m</td><td>4 t/ha</td><td>Loose</td><td>Tolerant: mold, anthracnose</td><td>Mali, IER/ICRISAT selection, pedigree FambeA x Lata.</td></tr><tr><td>CSM63E</td><td>G5</td><td>Guinea</td><td>Low</td><td>90 days</td><td>600-1000 mm</td><td>Grain</td><td>4 m</td><td>2 t/ha</td><td>Loose</td><td>Tolerant: diseases and insects</td><td>Mali, traditional variety</td></tr><tr><td>SK5912</td><td>G6</td><td>Caudatum</td><td>High</td><td>110 days</td><td>700-900 mm</td><td>Biomass</td><td>2 m</td><td>2.5-3.5 t/ha</td><td>Semi-compact</td><td>Tolerant: mold, anthracnose</td><td>Nigeria</td></tr><tr><td>Grinkan</td><td>G7</td><td>Caudatum</td><td>No</td><td>110 days</td><td>500-800 mm</td><td>Grain-biomass</td><td>1.2 m</td><td>4 t/ha</td><td>Semi-compact</td><td>Resistant: midges, insects</td><td>Mali, ICRISAT selection</td></tr><tr><td>Soumba</td><td>G8</td><td>Caudatum</td><td>Low</td><td>100 days</td><td>600-1000 mm</td><td>Grain-biomass</td><td>2.5 m</td><td>2.5 t/ha</td><td>Semi-compact</td><td>Tolerant: diseases and, insects, striga</td><td>Mali</td></tr><tr><td>621B</td><td>G9</td><td>Caudatum</td><td>No</td><td>90 days</td><td>600-900 mm</td><td>Grain</td><td>1.75 m</td><td>2.5-3 t/ha</td><td>Semi-compact</td><td>Mold resistant</td><td>Senegal, ISRA selection, pedigree CE 151-262 xSarvato-1</td></tr><tr><td>F2-20</td><td>G10</td><td>Caudatum</td><td>Low</td><td>110 days</td><td>600-900 mm</td><td>Grain</td><td>2.1m</td><td>3-5.3 t/ha</td><td>Semi-compact</td><td>Resistant: mold, striga</td><td>Senegal, ISRA selection, pedigree (MN1056 x 68-20) x 7410-195-1</td></tr></tbody></table>
## 2.3. Trial Management

Tillage at each site consisted of cross plowing with discs (depth about 25 cm) followed by harrowing. The seeds were treated with Granox (a combination of captafol-benomyl and carbofuran). Sowing was always carried out after a good rain event (Table 1). Crops were sown in hills with 5-6 seeds per hole, with 0.80 m spacing between rows and 0.20 m spacing between the hills along each row (i.e., planting density of 62,500 hills ha⁻¹). Around 15 days after emergence, the plots were thinned to one plant per hill. Mineral fertilizers were applied according to the research institute recommendation in Senegal: 150 kg ha⁻¹ of N-P-K (15-10-10) at sowing or emergence, and 100 kg ha⁻¹ urea applied twice, 50 kg ha⁻¹ just after thinning and 50 kg ha⁻¹ during vegetative growth. Weeding, pesticide and insecticide treatments (Decis and Dimethoate), and protection against birds, were provided as required to minimize the impact on crop growth and grain loss. Anti-erosion bunds were installed around the trials to limit runoff.

## 2.4. Experimental Design and Data Collection

All the trials were laid out in the same Randomized Complete Blocks design with four replicates (randomization was different from one trial to another), each containing the ten genotypes. This gave a total of 40 plots per trial. Each plot consisted of 7 rows of 40 hills, occupying an area of 44.8 m² (5.6 m × 8 m). At physiological maturity, plants were sampled from a well delimited 3.36 m² (3 rows × 7 hills) sub-plot to assess biomass (leaves and stems) and grain yields. Biomass dry weight was determined after air-drying in a greenhouse, followed by 48 h in an oven at 65 °C, and grain dry weight after panicle threshing. Grain and dry biomass yields were calculated in kg ha⁻¹.

## 2.5. Data Analysis

An initial analysis of variance was performed for each environment – defined in this study as an experimental situation, i.e., a site-year-seedling-date combination (11 in total), to verify the existence of differences between genotypes. Subsequently, a combined analysis of variance was conducted, considering the effect of the genotype and the environment as fixed, according to the following statistical model:

$$Y_{ijk} = \mu + G_i + E_j + B_k(E_j) + (GE)_{ij} + \varepsilon_{ijk} \quad (1)$$

where $Y_{ijk}$ represents the ith genotype in the jth environment and the kth block; $\mu$ is the overall mean; $B_k(E_j)$ corresponds to the block within the jth environment and in the kth block; $G_i$ is the effect of the ith genotype; $E_j$ is the effect of the jth environment; $(GE)_{ij}$ is the effect of interaction of the ith genotype with the jth environment; and $\varepsilon_{ijk}$ is the effect of experimental error. The homogeneity between residual variances was tested using Bartlett’s test [23].

Lastly, adaptability and phenotypic stability analyses were performed by the AMMI method as described in Zobel et al. [12] using the following statistical model:

$$Y_{ij} = \mu + g_i + e_j + \sum_{k=n}^{n} \lambda_k \alpha_{ik} y_{jk} + r_{ij} + \varepsilon_{ij} \quad (2)$$

where $Y_{ij}$ is the mean response of genotype i in environment j; $\mu$ is the overall mean; $g_i$ is the fixed effect of genotype i ($i = 1, 2, \dots, g$); $e_j$ is the fixed effect of environment j ($j = 1, 2, \dots, e$); $\varepsilon_{ij}$ is the average experimental error; $G \times E$ interaction is represented by the factors; $\lambda_k$ is a unique value of the kth interaction principal component axis (IPCA), ($k = 1, 2, \dots, p$, where p is the maximum number of estimable main components), $\alpha_{ik}$ is a singular value for the ith genotype in the kth IPCA, $y_{jk}$ is a unique value of the jth environment in the kth IPCA; $r_{ij}$ is the error for $G \times E$ interaction.

The sum of squares for $G \times E$ interaction was divided into n singular axes or main components of interaction (IPCA), which described the standard portion (ANOVA), with each axis corresponding to an AMMI model. Generally, when $G \times E$ interactions are significant, AMMI models with one or
two main axes (AMMI₁ and AMMI2 models respectively) are the most commonly used because of their simplicity in biplot graph representations. Biplot graph interpretation is based on the variation of the additive main effects (genotype and environment) and the multiplier effect of G × E interaction. According to Zobel et al. [12], for the AMMI₂ graph, genotypes that have low scores on IPCA₁ (first interaction principal component axis) or IPCA₂ (second interaction principal component axis) or both, contribute little to the interaction. This indicates a general adaptation. On the other hand, those with high scores, be it positive or negative, have strong interactions and are specifically adapted to the environment that has the same sign score.

To identify genotypes showing the best trade-offs between grain and biomass yield (dual-purpose potential), genotype performance was compared to the overall mean in a scatterplot via the IPO index described below:

$$I_{PO} = \frac{Y_{ij} - \bar{Y}}{\bar{Y}} \qquad (3)$$

where IPO = potential index of a given genotype i for grain (or biomass) yield for a given environment j; Yij = grain (or biomass) yield of a given genotype i for a given environment j; Y = Overall mean grain (or biomass) yield (all environments and genotypes included). Thus, for a given environment j, a positive IPO (IPO > 0) of a genotype i for grain or biomass yield indicates the good potential of this genotype i for this environment j. Conversely if the IPO of a genotype i is negative for grain or biomass yield, it indicates poor potential for grain or biomass, respectively. Positive (or negative) IPO values for a given genotype i for both grain and biomass yield will therefore indicate good (or poor) dual-purpose potential. Principal component analysis (PCA) associated with a Hierarchical clustering analysis were performed for the characterization of the study environments. All statistical analyses were performed using R version 3.2 software [24].

# 3. Results

## 3.1. Environment Characterization

The environments were characterized by quantitative and qualitative indicators of soil fertility, rainfall distribution during the growing cycle and the overall presence of diseases (chlorosis and plant necrosis causing heterogeneity in the field) (Table 1). The values of the indicators for each environment were summarized through a principal component analysis: this showed that 68.7 % of the initial information provided was returned (Figure 2). The PC1 axis tends to be correlated to soil fertility and the presence of diseases and the PC2 axis tends to be correlated to certain rainfall variables. The environments studied could be classed in six groups according to the two axes, PC1 and PC2:

Group 1: environments characterized by very good soil fertility, an absence of disease and low overall rainfall, but well distributed. Only S16D1 belonged to this group

Group 2: environments with relatively good soil fertility and almost no disease, and high rainfall. Only N15D2 belonged to this group

Group 3: characterized by very humid environments throughout the cycle, low soil fertility and the presence of diseases at a moderate level. Only N15D1 belonged to this group;

Group 4: environments characterized by good total rainfall, an early end of rainfall, relatively good soil fertility and an absence of disease. This group included environments S13D1, S13D2 and S14D1

Group 5: environments with many constraints: very low soil fertility, high disease occurrence and low rainfall at the end of the cycle. Environments B13D1, B13D2 and S15D1 belonged to this group

Group 6: this group was characterized by low soil fertility, low rainfall during the cycle, low rainfall accumulation at the end of the cycle, but a lower disease occurrence compared to group 5. This group included environments S14D2 and S15D2.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Principal component analysis (PCA) of the characteristics of the environments studied. The environments and indicators are explained in Table 1.</figcaption></figure>

## 3.2. Effects of Genotypes, Environments and Genotype × Environment Interactions

The results of the combined ANOVA and of the AMMI are presented in Table 3. The genotype, environment and G × E interaction effects were significant for grain and biomass yields (*p* < 0.001). The mean grain yield of the genotypes ranged from 2018 kg ha⁻¹ (Nieleni) to 807 kg ha⁻¹ (SK5912). The genotypes performed differently in all the environments, except Fadda and Nieleni, which performed relatively better in all environments (Figure 3). Three genotypes had a higher mean yield than the overall mean (1454 kg ha⁻¹): Nieleni, Fadda and Pablo, with yields of 2018 kg ha⁻¹, 1833 kg ha⁻¹, and 1615 kg ha⁻¹, respectively. The three genotypes with the poorest performance were F2-20, Grinkan and SK5912, with mean grain yields of 1333 kg ha⁻¹, 1281 kg ha⁻¹ and 807 kg ha⁻¹, respectively. Mean grain yields across the environments (Table 4) ranged from 530 kg ha⁻¹ (B13D2) to 2313 kg ha⁻¹ (S16D1). Six of the eleven environments exceeded the overall mean: S16D1 (2313 kg ha⁻¹), N15D2 (1,766 kg ha⁻¹), S13D1 (1714 kg ha⁻¹), S14D1 (1696 kg ha⁻¹), N15D1 (1610 kg ha⁻¹) and S13D2 (1570 kg ha⁻¹).

<table><caption>**Table 3.** Summary of the combined analysis of variance and decomposition of G × E interaction according to AMMI.</caption><thead><tr><th rowspan="2">Source of Variation</th><th colspan="3">Grain (kg ha<sup>-1</sup>)</th><th colspan="3">Biomass (kg ha<sup>-1</sup>)</th></tr><tr><th>DF</th><th>Mean Square</th><th>TSS Explained (%)</th><th>DF</th><th>Mean Square</th><th>TSS Explained (%)</th></tr></thead><tbody><tr><td>Genotype (G)</td><td>9</td><td>3,990,633 ***</td><td>17.9</td><td>9</td><td>178,164,830 ***</td><td>36.7</td></tr><tr><td>Environment (E)</td><td>10</td><td>7,936,033 ***</td><td>39.6</td><td>10</td><td>92,439,498 ***</td><td>21.2</td></tr><tr><td>Blocks (E)</td><td>33</td><td>523,880 ***</td><td>8.6</td><td>33</td><td>17,553,265 ***</td><td>13.3</td></tr><tr><td>Interaction (G × E)</td><td>89</td><td>759,922 ***</td><td>33.8</td><td>89</td><td>14,146,802 ***</td><td>28.8</td></tr><tr><td>IPCA<sub>1</sub></td><td>18</td><td>1,371,515 ***</td><td>36.6</td><td>18</td><td>32,487,030 ***</td><td>52</td></tr><tr><td>IPCA<sub>2</sub></td><td>16</td><td>1,117,060 ***</td><td>26.5</td><td>16</td><td>17,004,129 ***</td><td>24.2</td></tr><tr><td>IPCA<sub>3</sub></td><td>14</td><td>1,011,953 ***</td><td>21</td><td>14</td><td>10,009,335 *</td><td>12.5</td></tr><tr><td>IPCA<sub>4</sub></td><td>12</td><td>396,066</td><td>7</td><td>12</td><td>4,671,580</td><td>5</td></tr><tr><td>IPCA<sub>5</sub></td><td>10</td><td>386,334</td><td>5.7</td><td>10</td><td>3,398,772</td><td>3</td></tr><tr><td>Error</td><td>289</td><td>231,846</td><td></td><td>287</td><td>4,711,029</td><td></td></tr></tbody></table>

DF = degrees of freedom; **, * = significant at 0.1% and 5%, respectively; TSS = total sum of squares.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Mean grain yield (mean of replications) (kg ha⁻¹) of each of the ten genotypes in each of the eleven environments studied and values of the genotypic, environmental and overall means. (B13D1 = sowing date 1 Bambey, B13D2 = sowing date 2 Bambey, S13D1 = sowing date 1 Sinthiou Malem, S13D2 = sowing date 2 Sinthiou Malem, S14D1 = sowing date 1 Sinthiou Malem, S14D2 = sowing date 2 Sinthiou Malem, S14D1= sowing date 1 Sinthiou Malem, S15D2 = sowing date 2 Sinthiou Malem, S16D1= sowing date 1 Sinthiou Malem, N15D1= sowing date 1 Nioro du Rip, N15D2 = sowing date 2 Nioro Rip. Figures 13, 14, 15 and 16 correspond to the years 2013, 2014, 2015 and 2016. Gen. mean = genotypic mean, Env. mean = environmental mean).</figcaption></figure>

<table><caption>**Table 4.** Mean grain yield (kg ha⁻¹) of ten genotypes grown at eleven environments.</caption><thead><tr><th rowspan="2">Genotype</th><th colspan="10">Environment</th><th rowspan="2">Genotypic Mean</th></tr><tr><th>B13D1</th><th>B13D2</th><th>N15D1</th><th>N15D2</th><th>S13D1</th><th>S13D2</th><th>S14D1</th><th>S14D2</th><th>S15D1</th><th>S15D2</th><th>S16D1</th></tr></thead><tbody><tr><td>Fadda</td><td>1662</td><td>804</td><td>2417</td><td>1719</td><td>2329</td><td>1855</td><td>1604</td><td>2206</td><td>1634</td><td>2077</td><td>1857</td><td>1833</td></tr><tr><td>Nieleni</td><td>2011</td><td>972</td><td>1122</td><td>2824</td><td>2445</td><td>1742</td><td>2626</td><td>1326</td><td>2049</td><td>1871</td><td>2946</td><td>2018</td></tr><tr><td>IS15401</td><td>665</td><td>554</td><td>1431</td><td>1310</td><td>2182</td><td>1608</td><td>1297</td><td>1008</td><td>1524</td><td>1958</td><td>1883</td><td>1402</td></tr><tr><td>Pablo</td><td>1688</td><td>624</td><td>1432</td><td>1796</td><td>2111</td><td>1786</td><td>1592</td><td>1806</td><td>1435</td><td>1716</td><td>1780</td><td>1615</td></tr><tr><td>CSM63E</td><td>1423</td><td>346</td><td>1791</td><td>478</td><td>1628</td><td>1895</td><td>1939</td><td>1634</td><td>232</td><td>1247</td><td>2345</td><td>1360</td></tr><tr><td>SK5912</td><td>252</td><td>151</td><td>2050</td><td>1754</td><td>992</td><td>1071</td><td>459</td><td>358</td><td>477</td><td>503</td><td>-</td><td>807</td></tr><tr><td>Grinkan</td><td>888</td><td>502</td><td>1929</td><td>1707</td><td>1323</td><td>1475</td><td>1677</td><td>441</td><td>1171</td><td>881</td><td>1905</td><td>1281</td></tr><tr><td>Soumba</td><td>931</td><td>553</td><td>1381</td><td>2365</td><td>1016</td><td>1301</td><td>2307</td><td>900</td><td>1149</td><td>1339</td><td>2412</td><td>1443</td></tr><tr><td>621B</td><td>1367</td><td>491</td><td>1342</td><td>1503</td><td>1566</td><td>1665</td><td>1810</td><td>1064</td><td>533</td><td>572</td><td>3233</td><td>1392</td></tr><tr><td>F2-20</td><td>1223</td><td>409</td><td>1207</td><td>2205</td><td>1549</td><td>1302</td><td>1658</td><td>1110</td><td>1012</td><td>536</td><td>2453</td><td>1333</td></tr><tr><td>Mean</td><td>1211</td><td>530</td><td>1610</td><td>1766</td><td>1714</td><td>1570</td><td>1697</td><td>1192</td><td>1122</td><td>1270</td><td>2313</td><td>1454</td></tr></tbody></table>

With respect to biomass yield, mean yields per genotype ranged from 10,478 kg ha⁻¹ (IS15401) to 4384 kg ha⁻¹ (621B). The genotypes performed differently in all the environments, except IS15401 and SK5912, which performed relatively better in all the environments (Figure 4). Three out of the ten genotypes had a higher mean yield than the overall mean (6954 kg ha⁻¹): IS15401, SK5912 and Fadda, with respective values of 10,364 kg ha⁻¹, 10,115 kg ha⁻¹ and 7995 kg ha⁻¹. Mean yields across the environments ranged from 9536 kg ha⁻¹ (B13D2) to 4923 kg ha⁻¹ (S15D2). Five environments exceeded the overall mean: B13D2, S14D1, B13D1, S16D1 and N15D1, with respective yields of 9533 kg ha⁻¹, 9129 kg ha⁻¹, 8055 kg ha⁻¹, 7852 kg ha⁻¹ and 7660 kg ha⁻¹ (Table 5).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Mean biomass yield (mean of replications) (kg ha⁻¹) of each of the ten genotypes in each of the eleven environments studied and values of the genotypic, environmental and overall means.</figcaption></figure>

<table><caption>**Table 5.** Mean biomass yield (kg ha⁻¹) of ten genotypes grown at eleven environments.</caption><thead><tr><th rowspan="2">Genotype</th><th colspan="10">Environment</th><th rowspan="2">Genotypic Mean</th></tr><tr><th>B13D1</th><th>B13D2</th><th>N15D1</th><th>N15D2</th><th>S13D1</th><th>S13D2</th><th>S14D1</th><th>S14D2</th><th>S15D1</th><th>S15D2</th><th>S16D1</th></tr></thead><tbody><tr><td>Fadda</td><td>11,111</td><td>10,546</td><td>10,857</td><td>6496</td><td>7571</td><td>5409</td><td>10,508</td><td>4966</td><td>6056</td><td>5447</td><td>8972</td><td>7995</td></tr><tr><td>Nieleni</td><td>7667</td><td>13,322</td><td>3990</td><td>4004</td><td>5617</td><td>6473</td><td>8141</td><td>7277</td><td>5895</td><td>5360</td><td>8509</td><td>6784</td></tr><tr><td>IS15401</td><td>12,315</td><td>15,989</td><td>5655</td><td>10,803</td><td>10,611</td><td>7151</td><td>17,077</td><td>8211</td><td>5841</td><td>6276</td><td>14073</td><td>10364</td></tr><tr><td>Pablo</td><td>8137</td><td>7198</td><td>4712</td><td>5855</td><td>7109</td><td>5775</td><td>8728</td><td>5982</td><td>6822</td><td>5431</td><td>7460</td><td>6655</td></tr><tr><td>CSM63E</td><td>3643</td><td>3051</td><td>5376</td><td>3134</td><td>4926</td><td>4529</td><td>7650</td><td>3842</td><td>4129</td><td>3591</td><td>6460</td><td>4576</td></tr><tr><td>SK5912</td><td>14,806</td><td>13,623</td><td>11,591</td><td>8156</td><td>9827</td><td>7480</td><td>11,870</td><td>8783</td><td>6867</td><td>7332</td><td>-</td><td>10,115</td></tr><tr><td>Grinkan</td><td>9675</td><td>12,020</td><td>9638</td><td>3489</td><td>6870</td><td>4812</td><td>8295</td><td>5297</td><td>4882</td><td>4682</td><td>7094</td><td>6860</td></tr><tr><td>Soumba</td><td>4756</td><td>7497</td><td>6669</td><td>4483</td><td>4034</td><td>5109</td><td>7196</td><td>5718</td><td>5991</td><td>3290</td><td>5820</td><td>5459</td></tr><tr><td>621B</td><td>3581</td><td>5955</td><td>4969</td><td>3003</td><td>4300</td><td>4889</td><td>4258</td><td>6829</td><td>3308</td><td>2888</td><td>5413</td><td>4379</td></tr><tr><td>F2-20</td><td>4863</td><td>7219</td><td>13,140</td><td>4061</td><td>5261</td><td>5791</td><td>7568</td><td>7602</td><td>4833</td><td>4931</td><td>6870</td><td>6558</td></tr><tr><td>Mean</td><td>8055</td><td>9536</td><td>7660</td><td>5348</td><td>6613</td><td>5742</td><td>9129</td><td>6431</td><td>5426</td><td>4923</td><td>7852</td><td>6954</td></tr></tbody></table>

The AMMI analysis of variance of ten genotypes tested in eleven environments for grain yield showed that the main effect of genotypes and environments accounted for 17.9% and 39.6% of the variation respectively, and the G × E interaction effect amounted to 33.8%. For biomass yield, 36.7%, 21.2% and 28.8% of the total sum of squares were attributed to genotype, environment and G × E interaction effects, respectively. For the decomposition of the G × E interaction according to the AMMI model, the analysis showed that the first two main components of the interaction were significant (Table 4) for both yields. The first two main components explained 60.3% and 76.2% respectively of the sum of squares for grain and biomass yields (IPCA₁ and IPCA₂). These results indicated that genotype and environment scores on the first two main components of the interaction explained almost all of the interaction that occurred in the data matrix.

### 3.3. Which Genotype(s) for Which Environment(s)?

The AMMI2 biplot graph for grain yield shows that the S16D1, N15D1 and N15D2 environments best discriminated against the performance of the different genotypes evaluated because of their high score. They significantly contributed to interaction (Figure 5). However, the mean yields for these environments were among the highest, indicating that they were environments that were conducive to achieving high yields. The main reason for these high yields in these cited environments was the relatively high soil fertility (S16D1; group 1 in the characterization of environments; see Figure 2),
or good rainfall associated with relatively good soil fertility (cases of N15D1 and N15D2, belonging to groups 2 and 3 respectively). Similarly, environments S15D1, S15D2 and S14D2 were discriminating, but produced the lowest grain yields. They were characterized by low soil fertility, high disease occurrence (S15D1), and low rainfall at the end of the cycle (S14D2). However, B13D2 contributed significantly less to interaction and was the main factor contributing to the phenotypic stability of these genotypes (Figure 5). This environment had one of the lowest mean grain yields because of its high disease occurrence, low soil fertility and low rainfall at the end of the cycle.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** AMMI biplot of grain yield for the ten sorghum genotypes and eleven environments studied.</figcaption></figure>

Due to their position along both axes (scores close to zero), Grinkan and SK5912 were the most stable genotypes, but with the lowest yield. In contrast, genotypes 621B, Soumba, F2-20, IS15401 and CSM63E were very unstable due to their values far from the origin of the IPCA axis, with grain yields lower than the overall mean. Lastly, the Fadda, Nieleni and Pablo genotypes were also found to be unstable, but displayed high grain yields. Genotypes and environments close to each other in the biplot had positive associations, indicating specific adaptation. For instance, genotypes 621B and Nieleni had specific adaptations for the S16D1 and N15D2 environments, respectively. Likewise, S14D1 was found to be a suitable environment for Soumba and F2-20, S15D2 for genotypes IS15401 and Fadda, and lastly N15D1 for genotype CSM63E (Figure 5).

For biomass yield, environments B13D1, B13D2, S14D1 and N15D1, all with yields above the overall mean, contributed significantly to interaction, as indicated by values far from the origin of the IPCA axis (Figure 6). These environments were characterized by early sowing dates, with relatively good soil fertility (S14D1) and a very long and well-distributed rainy season (N15D1) conducive to biomass production (note that late cycle stress in B13D2 did not affect biomass production). Environments S13D2 and S15D1, very close together on the biplot (Figure 6), influenced genotypes in the same way, all with biomass yields lower than the overall mean. These environments were characterized by an early end to the season (S13D2) and low soil fertility, and the presence of disease (S15D1) affecting biomass production. In contrast, S13D1, S16D1 and N15D2 showed a smaller contribution to G × E interaction. These environments were the main contributors to the phenotypic stability of the genotypes (Figure 6). In addition, these environments recorded distinct levels of performance: above the overall mean for S16D1, close to the overall mean for S13D1, and below the overall mean for N15D2. In these environments, there was no occurrence of disease and relatively good soil fertility.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** AMMI biplot of biomass yields for the ten sorghum genotypes and eleven environments studied (B13D1 = Sowing date 1 Bambey, B13D2 = sowing date 2 Bambey, S13D1 = sowing date 1 Sinthiou Malem, S13D2 = sowing date 2 Sinthiou Malem, S14D1 = sowing date 1 Sinthiou Malem, S14D2 = sowing date 2 Sinthiou Malem, S14D1 = sowing date 1 Sinthiou Malem, S15D2 = sowing date 2 Sinthiou Malem, S16D1 = sowing date 1 Sinthiou Malem, N15D1 = sowing date 1 Nioro du Rip, N15D2 = sowing date 2 Nioro Rip).</figcaption></figure>

Due to their positions located near the origin of the biplot, some genotypes can be considered as stable (i.e., Soumba), or unstable with specific adaptations (i.e., IS15401, SK5912 and Fadda, F2-20, Nieleni and Pablo). Of these genotypes, IS15401, SK5912 and Fadda showed yields higher than the overall mean. SK5912 and Fadda displayed specific adaptation to environment B13D1. Likewise, N15D1 was found to be a suitable environment for F2-20, S14D1 for IS1540, and S15D1 and S13D2 for Nieleni and Pablo. On the other hand, IS15401 and SK5912 proved to be very poorly adapted to environments S15D1, S15D2, S13D2 and S14D2 (Figure 6). These environments were characterized by low fertility (except S14D2), high disease occurrence (especially S15D1), late sowing, and low rainfall at the end of the cycle (S14D2).

In general, the Fadda, Nieleni, IS15401, F2-20 and Pablo genotypes were very unstable for both grain and biomass yield. These genotypes showed specific adaptations to the Sinthiou Malem and Nioro du Rip environments for grain, and especially to the environments sown on date 1 for biomass yield. Of these genotypes, Fadda, Nieleni, IS15401 and Pablo were generally successful for both grain and biomass yield. The Soumba and Grinkan genotypes were stable across environments, but did not produce well.

### 3.4. Which Genotype(s) Showed Dual-Purpose Potential?

The dual-purpose potential indices for grain and biomass production in the eleven study environments for the ten genotypes are shown in Figure 7. Regardless of the environment (excluding B13D2), five genotypes had consistently higher indices for grain (Fadda, Nieleni and Pablo) or for biomass (IS15401 and SK5912). Supporting this result, these five genotypes also showed higher mean yields for grain and/or biomass production according to the AMMI analysis (Figures 3 and 4). Moreover, of these genotypes Fadda was the one that combined the best grain and biomass production in the most numerous environments. For dual-purpose potential, Fadda was therefore well positioned (in the upper right quadrant) in five environments, Nieleni and Pablo in three environments, and IS15401 in two environments. It should be noted that all the environments where these genotypes expressed dual-purpose potential were sown early (i.e., date 1).
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Potential index for dual production of the ten genotypes across the eleven study environments. (The dotted lines represent the lines of equations *x* = 0 and *y* = 0).</figcaption></figure>

For grain yield, Fadda ranked well (upper quadrants in Figure 7) in ten environments, Nieleni in eight, Pablo and CSM63E in six, and IS15401 in four. For biomass, SK5912 performed well (the right-hand side quadrants of the figures) in nine environments, IS15401 in seven environments, Fadda in six environments and Nieleni and Pablo in four environments.

Some other genotypes, such as CSM63E, Grinkan, Soumba, 621B and F2-20, all with mean grain and biomass yields lower than the overall mean, all showed poor dual-purpose potential (lower-left quadrant) in several environments. This was particularly true for Soumba and 621B in half of the environments. They were stable and inefficient for biomass production (Figure 4) and never demonstrated dual-purpose potential in any of the 11 environments. The environments in which they performed poorly for dual production were mostly late sowing (date 2) and belonged to groups of environments affected by stress to which these genotypes were susceptible: N15D1 for group 3, S13D2 for group 4, B13D1, B13D2 and S15D1 for group 5, and S14D2 and S15D2 for group 6 (Figure 2).

**4. Discussion**

The joint analysis of variance showed differences in sorghum performance due to environments (greater for grain), genotypes (greater for biomass) and G × E interactions. This indicates that these genotypes responded differently to environments, thereby confirming phenotypic diversity among the genotypes assessed. These results are in agreement with previous findings on sorghum [25–28]. In general, yields were higher for sowing date 1 than for sowing date 2, showing the potential importance of a longer cycle time and the existence of a photoperiod response of several of the sorghums studied, as already demonstrated in West Africa [29]. Further, it also highlights the importance of the choice of genotypes if the main objective is to obtain biomass, and of the environments if it is to obtain grain.

The low grain yield observed in Bambey in 2013 (B13D1 and B13D2), and Sinthiou Malem in 2014 (S14D2) could be explained by, among other things, the particularly sandy nature of the soil and by a high occurrence of diseases and/or deficiencies (signs observed but not clearly identified), constraints to which the 621B, Grinkan, Soumba and F2-20 genotypes were more susceptible than the others. It should be noted that these four genotypes are all *caudatum*, which are improved genotypes introduced into national sorghum breeding programs and which are known to be less hardy than *guinea* when edaphic conditions are not ideal [30]. For the Bambey trials, water stress at the end of the cycle (due to problems with the irrigation system) also occurred. This might explain the very low yield
of SK5912 (very late maturing variety), as this early end of the rainy season might have aggravated the effects of the constraints mentioned above.

## 4.1. Which Genotype(s) for Which Environment(s)?

Various multiparametric models for measuring the stability of genotype performance across environments are available in the literature. Currently, the most widely used model is AMMI [31–33], which involves both an ANOVA and a principal component analysis (PCA) to decompose $G \times E$ interaction. The ability to identify genotypes with a stable performance and genotypes showing specific adaptations to specific environments is a major advantage of the AMMI method over other commonly used methods [34].

In this study, two genotypes (Grinkan and SK5912) for grain and two genotypes for biomass (Soumba and Grinkan) were identified as being generally more stable according to this model, but they also showed yields below the mean across different environments. These results support those of Menad et al. [35], who stated that the stability of yields is independent of their values, and that high-yielding genotypes are generally relatively unstable. They also confirm Yan and Hunt's [36] conclusions that global stability is not necessarily a positive factor and is only desirable when it combines a high mean yield. In addition, the AMMI analysis also revealed for grain that genotypes 621B, F2-20, Fadda and CSM63E were close to environments S16D1, S14D1, S15D2 and N15D1, respectively. For biomass, genotypes IS15401 and SK5912 were close to environments S14D1 and B13D1, respectively, indicating specific adaptations. In contrast, IS15401 and SK5912 proved to be very poorly adapted to environments S15D1, S15D2, S13D2 and S14D2.

The specific adaptation of 621B to environment S16D1 might be explained by the good agro-pedological conditions that this environment benefited from. Thus, 621B was found to be particularly susceptible to growing conditions, particularly soil conditions. This result agrees with previous results [37]. The specific adaptation of IS15401 to environment S14D1 may have been due to both the rather early sowing date and the good agro-pedological conditions, which allowed this genotype to perform better despite low rainfall. The adaptation of SK5912 to B13D1 might be explained by the early sowing date and the fact that the diseases that affected the other genotypes (Figure 2) did not affect its biomass production. The poor adaptation of IS15401 and SK5912 to S15D1 and S14D2 for biomass production might have been due on the one hand to low fertility and the occurrence of diseases in S15D1, which attacked these genotypes, and on the other hand to late sowing and low rainfall at the end of the cycle observed in S14D2. In general, the differences in performance of our genotypes across our studied environments could be attributed to the type of soil (i.e., texture, with all genotypes performing better in clay-textured soils containing more organic matter), rainfall, their genetic nature and biotic constraints. These factors affected grain and biomass yields to different degrees. These results demonstrate that strategic choices must be made by breeders in the introduction of new sorghum genotypes. Breeders need to select specific lines according to their local environment and desired character.

## 4.2. Choice of Genotypes with Dual-Purpose Potential

The selection of genotypes showing a good trade-off between grain and biomass production, based on a comparison with the overall mean performance of the ten genotypes, revealed that the hybrids Fadda and Nieleni had the greatest dual-purpose potential followed by Pablo, then IS15401, which was relatively unstable for both grain and biomass production and therefore with specific adaptation. Fadda and Nieleni were identified as being particularly suitable for dual-purpose use in 5 and 3 environments, respectively, while Pablo and IS15401 were identified in 2 environments each (Figure 7). In addition, Nieleni appeared to be a poor producer of grains or biomass in only 3 and 6 environments, respectively. Fadda was found to be a poor producer of grains or biomass in only 1 and 5 environments, respectively. In contrast, Soumba and 621B showed poor grain and biomass production in half of the environments.
The higher dual-purpose potential of the three hybrids (Fadda, Nieleni and Pablo) definitely came from their genetic background, as already demonstrated for grain in Mali [38]. The added-value of this study was to study biomass production too, and to thus further recommend the two highest-yielding genotypes (Fadda and Nieleni) for both grain and biomass yield (dual-purpose) in our experimental zones (Sinthiou Malem and Nioro du Rip) with normal (early) sowing dates. However, it will also be important to investigate the good forage quality of these genotypes, with the hypothesis that Nieleni, a *caudatum* sorghum, will out-perform the Fadda genotype.

The results also showed that the dual-purpose potential of the genotypes was mainly expressed in environments with early sowing dates (date 1). Environments with late sowing dates tended to reduce biomass and grain production and were therefore not suitable for dual-purpose genotypes. These findings confirmed the merits of early sowing to improve the dual-purpose potential of genotypes. Similar results showing the beneficial effect of early sowing on sorghum performance in terms of biomass and/or grain production in different hot and dry growing environments in Asia, America and Africa were obtained by [39–42]. However, the last authors [42] did not find any significant effect of the sowing date on grain yield due to bird attacks that occurred in their experiments.

**5. Conclusions**

In this study, the AMMI model showed that grain and biomass yields were strongly influenced by genotypes, environments and genotype × environment interaction. The different environments resulted in different responses from the genotypes, with most of them displaying environmental adaptation. Soil fertility and rainfall during the experiment were major factors explaining the variation in genotype responses. Nieleni and Fadda had the highest grain yields, while IS15401 and SK5912 had the highest biomass yields. This study also showed that genotypes with good phenotypic stability performed poorly. The study found that out of the genotypes studied, those with the greatest dual-purpose potential were Nieleni and Fadda. An early sowing date was found to be beneficial for the expression of dual-purpose potential. In addition, the results indicated that Nieleni, Fadda, Pablo had the highest overall mean in terms of grain, whether the growing conditions were good or bad. This result is in line with those of various authors regarding the “superiority” of hybrids. Meanwhile, the Senegalese control genotype 621B appeared to have very good potential (being able to produce more than 3 tons per hectare under good agro-pedological conditions, such as those of S16D1), but it was particularly susceptible to growing conditions, particularly soil conditions. Hence, the AMMI statistical model showed its merits as a tool to help recommend sorghum genotypes: the Fadda and Nieleni varieties appeared to be the best genotypes for dual-purpose use in our study area. However, to gain a better understanding of these differences in genotypic performance depending on the environment (G × E interaction), further studies are needed to relate these results to the other phenology, morphology and growth characteristics of these genotypes, to explain in greater detail their phenotypic plasticity. In addition, future biochemistry and molecular biology analysis would be of great value to better understand the differences in genotypic performance depending on the environment; especially differences among grain type genotypes and biomass type genotypes.

**Author Contributions:** Conceptualization, M.N., M.A., N.C. and B.M.; methodology, M.N., M.A., and B.M.; software, M.N.; validation, M.A., B.M. and A.G.; formal analysis, M.N., K.K.G., M.A. and B.M.; investigation, M.N., K.K.G. and B.M.; resources, M.N., M.A., B.M. and N.C.; data curation, M.N., M.A. and B.M.; writing—original draft preparation, M.N.; writing—review and editing, M.A., B.M., K.K.G. and A.G.; visualization, M.A. and B.M.; supervision, M.A., B.M. and A.G.; project administration, B.M and M.A.; funding acquisition, B.M. and M.A.

**Funding:** This work constitutes part of doctoral research studies and was conducted with the financial support of the West Africa Agricultural Productivity Program for Senegal (WAAPP-Senegal).

**Acknowledgments:** The authors thank Mbaye Sarr Diop, Mamadou Sonko, Abdou Faye, Ibrahima Ndong, Massiré Badji, Bounama Mbengue, Yagouba Diao, Jean Gabriel Diatta and Thierno Ndiaye for their invaluable help and cooperation in the laboratory and field work.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Djè, Y.; Heuertz, M.; Ater, M.; Lefebvre, C.; Vekemans, X. Évaluation de la diversité morphologique des variétés traditionnelles de sorgho du Nord-ouest du Maroc. *Biotechnol. Agron. Soc. Environ.* **2007**, *11*, 30-40.

2. FAOSTAT 2015. *Statistiques des Données Année 2015*; United Nations Food and Agriculture Organization: Rome, Italy, 2015.

3. Agrama, H.A.; Tuinstra, M.R. La diversité phylogénétique et les relations entre le sorgho adhésions à l'aide SSR et RAPD. *Afr. J. Biotechnol.* **2003**, *2*, 334-340. [CrossRef]

4. Reddy, B.V.S.; Kumar, A.A.; Reddy, P.S.; Elangovan, M. Sorghum germplasm: Diversity and utilization. In *Sorghum Genetic Enhancement: Research Process, Dissemination and Impacts*; International Crops Research Institute for the Semi-Arid Tropics: Patancheru, Andhra Pradesh, India, 2008; pp. 153-169.

5. ANSD (Agence Nationale de la Statistique et de la Démographie). *Bulletin Mensuel des Statistiques Economiques de 2018*; Division des Statistiques Economiques Ministère de l'Economie, des Finances et du Plan: Dakar, Sénégal, 2018; p. 109.

6. Ba, K.; Tine, E.; Destain, J.; Cissé, N.; Thonart, P. Étude comparative des composés phénoliques, du pouvoir antioxydant de différentes variétés de sorgho sénégalais et des enzymes amylolytiques de leur malt. *Biotechnol. Agron. Soc. Environ.* **2010**, *14*, 131-139.

7. Seguin, B.; Soussana, J.F. Emissions de Gaz à Effet de Serre et Changement Climatique: Causes et Conséquences Observées Pour L'agriculture et L'élevage. *Courrierde L'environnement de l'INRA* **2008**, *55*, 79-91.

8. Cruz, C.D.; Regazzi, A.J. *Biometrical Models Applied to Plant Breeding*; Editora UFV: Viçosa, Brazil, 1997; p. 390.

9. Holhs, T. Analysis of genotype environment interactions. *S. Afr. J. Sci.* **1995**, *91*, 121-124.

10. Crossa, J. Statistical analyses of multilocation trials. *Adv. Agron.* **1990**, *44*, 55-85.

11. Romagosa, I.; Fox, P.N. Genotype x environment interaction and adaptation. In *Plant Breeding*; Springer: Dordrecht, The Netherlands, 1993; pp. 373-390. ISBN 9401046654.

12. Zobel, R.W.; Wright, M.W.; Gauch, H.G. Statistical analysis of a yield trial. *Agron. J.* **1988**, *80*, 388-393. [CrossRef]

13. Crossa, J.; Gauch, H.G.; Zobel, R.W. Additive main effect and multiplicative interaction analysis of two international maize cultivar trials. *Crop. Sci.* **1990**, *30*, 493-500. [CrossRef]

14. Gauch, H.G.; Zobel, R.W. AMMI analysis of yield trials. In *Genotype by Environment Interaction*; Kang, M.S., Gauch, H.G., Eds.; CRC Press: Boca Raton, FL, USA, 1996; pp. 85-122.

15. Kempton, R.A. The use of biplots in interpreting variety by environment interactions. *J. Agric. Sci.* **1984**, *103*, 123-135. [CrossRef]

16. Ferreira, D.F.; Demétrio, C.G.B.; Manly, B.F.J.; Machado, A.A.; Vencovsky, R. Statistical models in agriculture: Biometrical methods for evaluating phenotypic stability in plant breeding. *Cerne* **2006**, *12*, 373-388.

17. Vachaud, G.; Dancette, C.; Sonko, S.; Thony, J.L. Méthode de caractérisation hydrodynamique d'un sol non saturé: Application à deux types de sols du Sénégal. *Ann. Agron.* **1978**, *29*, 1-36.

18. Hamon, G. Mise en Œuvre et Critique des Méthodes de Caractérisation Hydrodynamique de la Zone non Saturée du sol. Application Aux Sols de Culture du Sénégal. Ph.D. Thesis, Institut de Mécanique, Grenoble, France, 1980.

19. Imberno, J. Variabilité Spatiale des Caractéristiques Hydrodynamiques d'un sol du Sénégal. Application au Calcul d'un Bilan Sous Culture. Ph.D. Thesis, USMG et Institut National Polytechnique de Grenoble, Grenoble, France, 1981.

20. Hamon, G. Caractérisation Hydrodynamique in situ d'un sol de Culture en Moyenne Casamance; Institut Sénégalais de Recherches Agricoles: Djibélor, Sénégal, 1978.

21. Baret, F. Caractérisation Hydrodynamique d'un Sol de la Région de Nioro; Institut Sénégalais de Recherches Agricoles (ISRA): Dakar, Sénégal, 1980.

22. Valet, S. Bilan Hydrique sous Cultures Dans les Essais Travail du Sol en Sols Sableux (Nioro) et en Sols Sablo-Argileux (Thysse); Institut Sénégalais de Recherches Agricoles: Kaolack, Sénégal, 1984; p. 16.

23. Bartlett, M.S. Properties of sufficiency and statistical tests. *Proc. R. Stat. Soc. Ser. A* **1937**, *160*, 268-282.

24. R Development Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2015.
25. Almeida Filho, J.E.; Tardin, F.D.; Daher, R.F.; Barbé, T.C.; Paula, C.M.; Cardoso, M.J.; Godinho, V.P.C. Stability and adaptability of grain sorghum hybrids in the off-season. *Genet. Mol. Res.* **2014**, *13*, 7626–7635. [CrossRef] [PubMed]

26. Ndiaye, M.; Adam, M.; Muller, B.; Guisse, A.; Cissé, N. Performances agronomiques et stabilité phénotypique de génotypes de Sorgho (*Sorghum bicolor* (L.) Moench) au Sénégal: Une étude des interactions génotypes-environnement. *Appl. Biosci.* **2018**, *125*, 12617–12629.

27. Ganyo, K.K.; Muller, B.; Gaglo, E.K.; Guisse, A.; Cissé, N. Optimisation du NPK et urée basée sur les informations climatiques pour accroitre la production du sorgho en zones soudano-sahéliennes du Sénégal. *Appl. Biosci.* **2018**, *131*, 13293–13307. [CrossRef]

28. Showemimo, F.A.; Echekwu, C.A.; Yeye, M.Y. Genotype x environment interaction in Sorghum trials and their implication for future variety evaluation in Sorghum growing areas of northern Nigeria. *Plant Sci.* **2000**, *1*, 24–31.

29. Kouressy, M.; Dingkuhn, M.; Vaksmann, M.; Heinemann, A.B. Adaptation to diverse semi-arid environments of sorghum genotypes having different plant type and sensitivity to photoperiod. *Agric. For. Meteorol.* **2008**, *148*, 357–371. [CrossRef]

30. Bazile, D.; Dembélé, S.; Soumaré, M.; Dembele, D. Utilisation de la diversité variétale du sorgho pour valoriser la diversité des sols au Mali. *Cah. Agric.* **2008**, *17*, 86–94. [CrossRef]

31. Raton, F.L.; Gauch, H.G. *Statistical Analysis of Regional Yield Trials. AMMI Analysis of Factorial Designs*; Elsevier: New York, NY, USA, 1992.

32. Raju, B.M.K. Study of AMMI model and its biplots. *J. Ind. Soc. Agric. Stat.* **2002**, *55*, 297–322.

33. Zali, H.; Farshadfar, E.; Sabaghpour, S.H.; Karimizadeh, R. Evaluation of genotype x environment interaction in chickpea using measures of stability from AMMI model. *Ann. Biol. Res.* **2012**, *3*, 3126–3136.

34. Silva Filho, J.L.; Morello, C.L.; Farias, F.J.C.; Lamas, F.M. Comparação de métodos para avaliar a adaptabilidade eestabilidade produtiva em algodoeiro. *Pesqui. Agropecuária Bras.* **2008**, *43*, 349–355. [CrossRef]

35. Menad, A.; Meziani, N.; Bouzerzour, H.; Benmahammad, A. Analyse de l'interaction génotype x milieux du rendement de l'orge (*Hordeum vulgare* L.): Application des modèles AMMI et la régression conjointe. *Nat. Biotechnol.* **2010**, *5*, 99–106.

36. Yan, W.; Hunt, L.A. Biplot analysis of multi-environment trial data. In *Kang MS: Quantitative Genetics, Genomics and Plant Breeding*; Louisiana State University: Baton Rouge, LA, USA, 2002; pp. 289–304.

37. Ganyo, K.K. Etude et Modélisation des Réponses de Variétés de Sorgho (*Sorghum bicolor* (L.) Moench) à des Stratégies Contrastées D'apports D'intrants. Ph.D. Thesis, Université Cheikh Anta Diop, Dakar, Senegal, 2019.

38. Rattunde, H.F.W.; Weltzien, E.; Diallo, B.; Diallo, A.G.; Sidibe, M.; Touré, A.O.; Rathore, A.; Das, R.R.; Leiser, W.L.; Touré, A. Yield of photoperiod-sensitive sorghum hybrids based on guinea-race germplasm under farmers' field conditions in Mali. *Crop. Sci.* **2013**, *53*, 1–8. [CrossRef]

39. Voigt, J.; Botha, P.R.; Gerber, H.S. The effect of planting date on the dry matter production of annual forage sorghum hybrids and hybrid millet cultivars. *Grassroots Newsl. Grassl. Soc. S. Afr.* **2008**, *8*, 18–24.

40. Erickson, J.E.; Helsel, Z.R.; Woodard, K.R.; Vendramini, J.M.B.; Wang, Y.; Sollenberger, L.E.; Gilbert, R.A. Planting date affects biomass and Brix of sweet sorghum grown for biofuel across Florida. *Agron. J.* **2011**, *103*, 1827–1833. [CrossRef]

41. Reddi, S.G.; Janawade, A.D.; Palled, Y.B. Influence of sowing dates on growth, grain and ethanol yield and economics. *Int. J. Agric. Sci. Vet. Med.* **2013**, *1*, 12–17.

42. Gutjahr, S.; Vaksmann, M.; Dingkuhn, M.; Thera, K.; Trouche, G.; Braconnier, S.; Luquet, D. Grain, sugar and biomass accumulation in tropical sorghums. *I. Trade-offs and effects of phenological plasticity*. *Funct. Plant Biol.* **2013**, *40*, 342–354.

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