Productivity and Feed Quality Performance of Napier Grass (Cenchrus purpureus) Genotypes Growing under Different Soil Moisture Levels

In the semi-arid and arid environments of Sub-Sharan Africa, forage availability throughout the year is insufficient and highly limited during the dry seasons due to limited precipitation. Thus, the identification of drought stress-tolerant forage cultivars is one of the main activities in forage development programs. In this study, Napier grass (Cenchrus purpureus), an important forage crop in Eastern and Central Africa that is broadly adapted to produce across tropical environments, was evaluated for its water use efficiency and production performance under field drought stress conditions. Eighty-four Napier grass genotypes were evaluated for their drought stress tolerance from 2018 to 2020 using agro-morphological and feed quality traits under two soil moisture stress regimes during the dry season, i.e., moderate (MWS) and severe (SWS) water stress conditions, and under rainfed conditions in the wet season (wet). Overall, the results indicated the existence of genotype variation for the traits studied. In general, the growth and productivity of the genotypes declined under SWS compared to MWS conditions. High biomass-yielding genotypes with enhanced WUE were consistently observed across harvests in each soil moisture stress regime. In addition, the top biomass-yielding genotypes produced the highest annual crude protein yield, indicating the possibility of developing high-feed-quality Napier grass genotypes for drought stress environments.


Introduction
Global warming and climate change have been described as potential threats to agricultural production and productivity due to increasing temperatures and declining and more erratic rainfall in the semi-arid and arid environments of Sub-Saharan Africa (SSA) [1]. The changes in precipitation and temperature levels trigger environmental stresses such as drought and heat stress that influence normal plant growth and hence minimize yield and quality. Forage crops are one of the feed sources for livestock production in SSA, but the availability of forage throughout the year is insufficient and much more limited in the dry seasons of the year due to inadequate precipitation [2]. Forage performance and production are highly associated with available soil moisture, so forage species adapted to semi-arid and arid environments need to be drought-tolerant and should offer a higher yield and feed quality potential to maximize the availability of feed in seasons with inadequate soil moisture [3,4]. Thus, the identification of drought-and heat stress-tolerant forage cultivars is one of the key strategies in forage development programs [5,6]. The evaluation of genotype performance under field conditions is one of the strategies for identifying tolerant and relatively productive genotypes in support of the development of varieties resilient to drought stress [7]. The morphological and physiological traits of plants are very important for selection to improve drought tolerance due to their relation to the adaptation of future climate scenarios.
Napier grass, or Elephant grass (Cenchrus purpureus (Schumach.) Morrone syn. Pennisetum purpureum Schumach.), is an important forage grass in tropical and sub-tropical environments, growing mainly from sea level up to 2000 m.a.s.l., with the best growth at temperatures ranging from 25 to 40 • C [8,9]. Napier grass is also considered a shortterm drought-tolerant forage, an important characteristic in areas that frequently face drought stress conditions [10]. It is a popular forage crop in tropical environments, mainly due to its high biomass production per unit area when compared with other tropical forages such as Guinea grass and Rhodes grass [11]. The International Livestock Research Institute (ILRI) forage gene bank conserves a range of Napier grass genotypes that are comprised of collections from ILRI and the Brazilian Agricultural Research Corporation (EMBRAPA) [11,12]. Diversity studies conducted on these collections have revealed the existence of genetic variation with the potential to support the development of varieties with desirable traits suitable for forage production in SSA. Phenotypic and genotypic diversity studies conducted on these collections revealed the existence of genetic variation with the potential to support sustainable forage production in the region and subsequently contribute towards improved animal performance [12]. Likewise, the field characterization of those collections indicated the existence of genotypic variations for yield and feed quality traits under both wet and dry season growing conditions, further indicating the possibility of selecting genotypes that fit the trending climatic conditions [13]. Generally, the high biomass-yielding capability and adaptability to wider environments of Napier grass make it an alternative forage species for livestock feed in moisture stress environments. Many cultivars of Napier grass, and hybrids with pearl millet (Cenchrus purpureus (Schumach.) Morrone × Cenchrus americanus (L.) Morrone), have been developed in different parts of the world for their high yield potential, broad adaptability, resistance to diseases and feed quality traits [14]. Despite the information on the adaptability of Napier grass to short-term drought stress, there have been no genotypes identified that can tolerate and produce biomass to support the availability of feed for livestock in the longer dry season conditions of tropical environments. Thus, the present study was conducted to evaluate Napier grass genotypes for enhanced water use efficiency (WUE) under field drought stress conditions using different soil moisture regimes in the dry season, from 2018 to 2020. Overall, the results of the analysis of agro-morphological and feed quality traits indicated the existence of significant phenotypic diversity among the experimental genotypes. Consistently high biomass-yielding genotypes with enhanced WUE were observed across all harvests in each soil moisture regime in the dry season. Furthermore, the additive main effects and multiplicative interaction (AMMI) analysis identified productive and stable genotypes which were consistent with the identified productive genotypes.

Effects of Growing Season, Harvest Period and Genotype on Napier Grass Performance
To analyze the performance of Napier grass genotypes, 12 traits representing morphological, physiological and agronomic parameters (Table 1) were investigated from the plants grown under moderate water stress (MWS) and severe water stress (SWS) conditions in the dry season as well as in the main rainy season (wet). The principal component analysis, conducted using all trait values from all growing conditions and harvests, showed a clear separation of genotype performance between the dry and rainy seasons, with relatively little variation in agronomic or feed quality traits between the MWS and SWS conditions ( Figure 1). Furthermore, analysis of variance (ANOVA) indicated highly significant (<0.001) variation among genotypes and harvests in each of the dry season and wet season conditions for all the traits considered, while the difference between treatments (MWS and SWS) was only highly significant (<0.001) in the dry season harvests ( Table 1). Regardless of the production season, the interaction effects among different genotypes and harvests were highly significant (<0.001) for all traits; however, the genotype-by-treatment effect was only significant in the dry season. This finding indicates that the plants grown under the SWS condition recovered more quickly with the onset of the rainy season. This recovery potential is an important response to determine the overall growth and development performance of the genotypes. Table 1. Combined ANOVA for agro-morphological traits of 84 Napier grass genotypes tested under rainy (wet) and dry (dry) season conditions.  ments (MWS and SWS) was only highly significant (<0.001) in the dry season harvests ( Table 1). Regardless of the production season, the interaction effects among different genotypes and harvests were highly significant (<0.001) for all traits; however, the genotypeby-treatment effect was only significant in the dry season. This finding indicates that the plants grown under the SWS condition recovered more quickly with the onset of the rainy season. This recovery potential is an important response to determine the overall growth and development performance of the genotypes. In addition, a separate ANOVA conducted for the dry season revealed significant (<0.05) and highly significant (<0.001) differences among genotypes, treatments and harvests, as well as in the interaction effects, for all traits except for leaf width (LW) and chlorophyll fluorescence (Fv/Fm) (Table A2), but in the rainy season, only the genotype and harvest effects were statistically highly significant (<0.001) (Table A2). Overall, the results showed that the performance of Napier grass was primarily affected by the genotype and harvest round in each of the wet, MWS and SWS growing conditions. The observed genotypic variations are important to exploit the potential of the genotypes to maximize forage production under different water stress environments. Table 1. Combined ANOVA for agro-morphological traits of 84 Napier grass genotypes tested under rainy (wet) and dry (dry) season conditions. In addition, a separate ANOVA conducted for the dry season revealed significant (<0.05) and highly significant (<0.001) differences among genotypes, treatments and harvests, as well as in the interaction effects, for all traits except for leaf width (LW) and chlorophyll fluorescence (Fv/Fm) (Table A2), but in the rainy season, only the genotype and harvest effects were statistically highly significant (<0.001) (Table A2). Overall, the results showed that the performance of Napier grass was primarily affected by the genotype and harvest round in each of the wet, MWS and SWS growing conditions. The observed genotypic variations are important to exploit the potential of the genotypes to maximize forage production under different water stress environments.

Partitioning Quantitative Genetic Variation
The phenotypic coefficient of variation (PCV), genotypic coefficient of variation (GCV) and broad sense heritability (H 2 ) were calculated to assess the contribution of these factors to the respective traits (Table 2). Generally, the PCV values are greater than the GCV values for all growing conditions, a finding that indicates that both genotypes and environmental factors contributed to the observed variation in the traits under investigation. The PCV and GCV values under the wet condition were higher than those under the dry season treatments (MWS/SWS) for the corresponding traits. In the wet condition, the highest PCV and GCV values were recorded for LSR (115% and 69%, respectively), followed by TFW (89% and 55%, respectively) and TDW (82% and 59%, respectively). In the wet condition, the PCV values for LSR, PH and IL were approximately two-fold higher than the GCV values, indicating that the environment plays a significant role in the variation observed among genotypes for these traits. Table 2. Variations and heritability in agro-morphological traits of Napier grass genotypes grown under wet and dry season conditions-either exposed to moderate water stress (MWS) or severe water stress (SWS)-for two years.

Traits
Growing The PCV and GCV values were similar to each other for the corresponding traits in the MWS and SWS conditions, which suggests that the effects of genotype and environmental factors on trait expression in both dry season treatments were similar. The traits with the highest PCV values, when growing under the MWS and SWS conditions, were WUE (80% and 89%, respectively), TFW (69% and 71%, respectively) and TDW (68% and 70%, respectively), and the corresponding GCV values for these traits when growing under these conditions were (39% and 41%, respectively), (41% and 42%, respectively) and (56% and 43%, respectively). The broad sense heritability estimates were higher for LW, TN, TFW and TDW when growing under the wet condition than those under both the MWS and SWS conditions, while the heritability for PH and LL was lower under the wet condition than that under both the MWS and SWS conditions. The heritability estimates of the corresponding traits were similar for the MWS and SWS conditions ( Table 2).

The Effect of Soil Moisture Stress Levels on Napier Grass Performance
The current study revealed that soil water stress conditions had a different effect on the physiological parameters, such as water use efficiency (WUE), as the genotype mean values for SWS were higher than the values for the corresponding genotypes grown under MWS. In addition, this observation showed the capability of different genotypes to show enhanced water use efficiency when exposed to different soil moisture stress conditions. Under MWS, the highest WUE was observed for genotypes 16819 (3.94), CNPGL 92-66-3 (3.73) and BAGCE 30 (3.72), and under SWS, the highest WUE was observed for genotypes 16819 (4.12), 16802 (4.1) and CNPGL 92-66-3 (3.98).
The mean trait values for wet, MWS and SWS conditions for 12 agro-morphological and physiological traits are presented in Supplementary Table S1. The genotype values of all traits decreased under the dry season treatments compared to the rainy season conditions, except for tiller number (TN). The physiological parameters that indicate photosynthetic efficiency, Fv/Fm and PI, were lower in the dry season (MWS and SWS) for corresponding genotypes, indicating that the imposed field water stress treatments impaired the normal physiological function of the plants. The lower Fv/Fm and PI mean values for genotypes growing under SWS compared to the corresponding genotypes growing under MWS support the assumption that the water stress severity under SWS was higher than that under MWS. This difference in severity level was shown by the lower mean values for genotypes grown under SWS, except for LW and LL, than the corresponding genotypes growing under MWS. However, the degree of decline was different for the various genotypes and at different harvests, as revealed by the genotype-by-treatment interaction effects (Supplementary Table S2).
A stress tolerance index (STI), one of the most commonly used drought stress tolerance indices [15], was used to estimate the tolerance levels of the genotypes when growing under both MWS and SWS conditions. Generally, the STI value was lower for genotypes growing under SWS conditions compared to the corresponding genotype growing under MWS conditions (Supplementary Table S2), a finding that reflects that the more severe the treatment conditions, the greater the influence on genotype performance. Under MWS, the highest STI was observed for genotypes 16819 (0.53), 16791 (0.46) and BAGCE 30 (0.45), and under SWS, the highest STI was observed for genotypes 16819 (0.49), 16791 (0.47) and CNPGL 93-37-5 (0.41). In addition, this observation showed the capability of different genotypes to show enhanced water use efficiency when exposed to different soil moisture stress conditions. Under MWS, the highest WUE was observed for genotypes 16819 (3.94), CNPGL 92-66-3 (3.73) and BAGCE 30 (3.72), and under SWS, the highest WUE was observed for genotypes 16819 (4.12), 16802 (4.1) and CNPGL 92-66-3 (3.98).
The WUE and STI profiles for Napier grass genotypes growing under both MWS and SWS conditions are presented in Figure 2. Regardless of the treatment under dry season conditions, genotypes showed a similar performance trend for both WUE and STI, indicating the potential importance of these traits for the screening and identification of genotypes for drought stress tolerance.
The WUE and STI profiles for Napier grass genotypes growing under both MWS and SWS conditions are presented in Figure 2. Regardless of the treatment under dry season conditions, genotypes showed a similar performance trend for both WUE and STI, indicating the potential importance of these traits for the screening and identification of genotypes for drought stress tolerance.

Genotype Diversity and Trait Selection under Water Stress Conditions
To evaluate the diversity of genotypes and select traits for both MWS and SWS growing conditions, the morphological, physiological and agronomic traits were subjected to principal component analysis (PCA). In both stress conditions, the first two principal components (PCs) explained approximately 80% of the genotype variation ( Figure 3). Under MWS, PC1 explained about 62% of the total variation, and this was principally associated with the traits STI, TDW, TFW, WUE and LL. PC2 under MWS explained about 18% of the total variation and was mostly associated with the traits PI, TN and Fv/Fm ( Figure 3A and Table A3). Similarly, under SWS, PC1 explained about 61% of the total variation and was associated with STI, TDW, TFW, WUE and LL, and PC2 explained about 19% of the   Figure 2. Water use efficiency (WUE-Bar graph) and stress tolerance index (STI-Line graph) profile of Napier grass genotypes grown under moderate water stress (light brown) or severe water stress (blue black) conditions. The order of genotypes from left to right is based on WUE rank (1 to 84) under severe water stress.

Genotype Diversity and Trait Selection under Water Stress Conditions
To evaluate the diversity of genotypes and select traits for both MWS and SWS growing conditions, the morphological, physiological and agronomic traits were subjected to principal component analysis (PCA). In both stress conditions, the first two principal components (PCs) explained approximately 80% of the genotype variation ( Figure 3). Under MWS, PC1 explained about 62% of the total variation, and this was principally associated with the traits STI, TDW, TFW, WUE and LL. PC2 under MWS explained about 18% of the total variation and was mostly associated with the traits PI, TN and Fv/Fm ( Figure 3A and Table A3). Similarly, under SWS, PC1 explained about 61% of the total variation and was associated with STI, TDW, TFW, WUE and LL, and PC2 explained about 19% of the total variation, mostly associated with PI, Fv/Fm and TN ( Figure 3B). Overall, the first principal component is mainly related to forage biomass and architecture traits, while the second principal component is related to physiological parameters that are indicative of photosynthetic efficiency.

Genotype and Trait Clusters under Water Stress Treatments
Biplot and hierarchical cluster analyses were conducted separately for both MWS and SWS conditions using TDW, WUE and STI to partition genotypes based on their drought stress response (Figures 4 and 5). Under both MWS and SWS conditions, the first two principal components (PC1 and 2) explained more than 90% of the genotype variation ( Figures 4A and 5A). The PCA categorized genotypes into highly drought stress-tolerant

Genotype and Trait Clusters under Water Stress Treatments
Biplot and hierarchical cluster analyses were conducted separately for both MWS and SWS conditions using TDW, WUE and STI to partition genotypes based on their drought stress response (Figures 4 and 5). Under both MWS and SWS conditions, the first two principal components (PC1 and 2) explained more than 90% of the genotype variation ( Figures 4A and 5A). The PCA categorized genotypes into highly drought stress-tolerant (high-yielding and water use-efficient), moderately tolerant and susceptible genotypes. Genotypes that were on the positive side of PC1 were high in WUE, TDW and STI, while genotypes on the negative side were low for the respective traits. Similarly, the hierarchical clustering categorized genotypes into highly tolerant, moderately tolerant and susceptible clusters in both soil moisture stress conditions. Furthermore, the drought stress response grouping of genotypes in both moisture stress conditions revealed that the highly tolerant and susceptible genotype groups were similar and showed the consistent performance of these genotypes regardless of the drought stress conditions ( Figures 4B and 5B). The performance of the moderately tolerant genotypes varied depending on the soil moisture stress levels.   In both stress levels, genotypes 16819, 16791, BAGCE 30, CNPGL 92-66-3 and BAGCE 93 were considered highly tolerant to drought stress, with consistently high TDW and WUE performances in both conditions, while genotypes 16834, 16805, 16621, 16790 and 16797 were the most susceptible.

Biomass Productivity of Napier Grass Genotypes under Water Stress and Optimum Soil Moisture Conditions
To assess the biomass productivity potential and maximize the availability of feed across seasons, genotypes grown under dry (MWS and SWS) and wet season conditions were compared by assessing TDW production. The top biomass-yielding genotypes identified in the wet season were 16791, 16819, BAGCE 30, CNPGL 93-37-5, similarly 16819, 16803, 16839 and BAGCE 30 under MWS and 16819, CNPGL 93-37-5, CNPGL 92-66-3 and 16839 under SWS (Supplementary Table S1). Furthermore, a general linear model regression between rainy and dry (MWS and SWS) season conditions for TDW revealed a positive correlation between rainy and MWS (R 2 = 0.53) conditions and between rainy and SWS (R 2 = 0.66) conditions ( Figure A1A,B). This result suggests that the high potential yield under optimal conditions could also result in improved yield under water stress conditions.

Biomass Yield Stability across Harvests
To improve the efficiency of selecting the best genotypes for wet, MWS and SWS conditions, Additive Multiplicative Interaction (AMMI) analysis was employed for TDW, where the sum of squares for the yield was partitioned into genotype (G), environment (E) and GE interaction (GEI), where the GEI was further partitioned by the interaction of principal components (IPCs) ( Table A4). The results of the AMMI analysis indicated a significant GEI effect, and further partitioning of the total sum of squares revealed that the genotype effect accounted for the largest proportion of the sum of squares for TDW, followed by E and then GEI, under both MWS and SWS conditions, while in the wet condition, the environment factor was the most significant, followed by genotype and GEI. The first two IPCs contributed the highest GEI partitioning for all growing conditions. The AMMI stability value (ASV) ranked the genotypes based on their score, with the lowest scores representing the most stable genotypes (Table A5). Under the wet condition, the most stable genotype for TDW yield was CNPGL 96-27-3, followed by 16789 and 16821, which had ASV rankings of 1, 2 and 3, respectively, while the least stable genotype was 16819, followed by 16791 and 16621, ranking 84, 83 and 82, respectively. Under MWS, the most stable genotypes for TDW yield were 15743, followed by 16806 and 16803, while the least stable genotype was 16839, followed by BAGCE 30 and CNPGL 93-37-5. Under SWS, the most stable genotype for TDW yield was 16792, followed by 15357 and 14389, while the least stable genotype was 16839, followed by BAGCE 93 and BAGCE 30. Thus, high-yielding genotypes were identified as less stable and vice versa based on ASV analysis. To include productivity as a selection criterion for drought stress tolerance, yield stability index (YSI) was used, as it encompasses the rank from both ASV and overall yield. The genotypes that had low YSI scores under the wet condition include CNPGL 92-66-3, 16839 and 14983, while for the MWS condition, genotypes 16803, 16806 and 14984 had the lowest scores; similarly, under SWS, genotypes 16792, CNPGL 92-66-3 and 14389 had the lowest YSI score.

Feed Quality Trait Variation among Genotypes, Soil Moisture Levels and Harvest
The combined ANOVA under wet, MWS and SWS conditions revealed highly significant genotypic and harvest round differences for feed quality traits (NDF, ADF, ADL, OM, CP, IVOMD and Me) ( Table 3). The observed differences between the genotypes, harvest round and genotypes × harvest indicated that the studied feed quality traits were greatly influenced by both genotype and harvest round (Table A6). Table 3. Combined ANOVA for seven feed quality traits of Napier grass genotypes grown under rainy season and dry season conditions across six harvests.

NDF
The PCV, GCV and broad-sense heritability (H 2 ) for feed quality traits are presented in Table 4. The mean values for the fiber components NDF, ADF and ADL were reduced in moisture stress conditions (MWS and SWS) compared to the wet condition; conversely, the mean value for CP, IVOMD and Me increased under moisture stress conditions (Supplementary Table S2). Similarly, a radar plot shows that almost all genotypes studied had reduced levels of fiber components under the stress conditions ( Figure A1A-C), while the levels of CP, IVOMD and Me increased ( Figure A1D-F). The PCV values were higher than the respective GCV values in all soil moisture conditions, indicating that the environment has a considerable effect on the feed quality traits (Table 4). Generally, PCV values were highest under the SWS condition, followed by the MWS and wet conditions, except for NDF and ADL, suggesting that the environmental impact on the expression of feed quality traits increased as soil moisture levels decreased. The maximum estimated value for PCV and GCV was for CP, followed by ADL and ADF, whereas the minimum estimated values were for OM and NDF. Heritability is a measure of the degree of influence of genotype and environment on the expression of the parameters. In our study, the estimated heritability values were medium to low; the heritability estimates for feed quality traits were lower in the SWS condition, except for NDF and ADL. Overall, the highest heritability was observed for OM and NDF for the studied genotypes. Table 4. Variations and heritability of the feed quality traits of 84 Napier grass genotypes grown under rainy season (wet) and dry season conditions-either exposed to moderate water stress (MWS) or severe water stress (SWS)-for two years.

Association of Feed Quality Traits and Total Dry Weight under Moisture Stress Conditions
The principal component analysis (PCA) showed that the first two principal components (PC) for MWS and SWS conditions accounted for 68.3% (44.4 and 23.9%) and 71.5% (49.3 and 22.2%) of the variation, respectively ( Figure A2). PC1 was positively correlated with ADF, ADL and NDF and negatively correlated with CP, IVOMD and Me. PC2 was mainly positively correlated with OM and NDF. The PCA further revealed the association between different feed quality traits and genotypes, as demonstrated by the PC biplots ( Figure A2A,B) for both water stress conditions. A smaller angle between different feed quality traits in the same direction indicated a high association between corresponding traits for classifying genotypes. Under both soil moisture stress conditions, the variations in IVOMD, CP, Me and ADF were the biggest contributors to PC1, while the variances for OM and NDF were the biggest contributors to PC2. Along PC1, genotypes located on the right side of the origin had high values for IVOMD, CP and Me, while the genotypes located on the left side had low values for these traits. Similarly, along PC2, genotypes at the top of the origin had high values for OM and NDF, while genotypes at the bottom revealed low values for these traits. Genotypes superior for a particular parameter were located close together and along the same direction of the vector line.

Annual Total Dry Weight and Crude Protein Yield Performance of Genotypes
To assess the forage value of genotypes, cumulative annual biomass and crude protein yield were evaluated. The range of annual TDW for the MWS condition was between 2.65-68.05 t/ha/year, with a mean value of 34.14 t/ha/year ( Table 5). The TDW range for the SWS condition was between 2.47-67.35 t/ha/year, with a mean value of 32.95 t/ha/year. Genotypes 16791, 16819, BAGCE 30, CNPGL 93-37-5 and 16802 produced the most annual TDW in both conditions. The genotypes that produced the highest annual TDW also produced the highest annual crude protein yield (CPY).

Variation between Leaf and Stem Tissue Samples for Feed Quality Traits in the Wet Season
The combined ANOVA of feed quality traits for leaf and stem tissue samples under the wet condition is presented in Table A7. All the studied feed quality traits showed highly significant genotypic and harvest round variation, while prior moisture stress treatments (MWS and SWS) in the dry season did not appear to affect them.
Leaf and stem tissue samples differed in their feed quality traits; generally, leaf samples had higher CP levels, while stem samples had higher fiber components (NDF and ADF) ( Figure 6 and Table A8). Interestingly, the study revealed genotypes that had stem samples with a high IVOMD and leaf samples with high ADL (Figure 6A,B).     Table S3).

Discussion
Drought stress is one of the main challenges that influence forage production in the arid and semi-arid tropics, particularly in the dry season. Developing forage cultivars that are high biomass-yielding and drought-tolerant is a priority in the tropical regions where droughts are increasing. Generally, biomass yield is determined by overall plant growth and developmental processes; thus, plant growth and development attributes are used to screen for drought adaptability. Traditionally, forage breeders utilize morphological and physiological indicators to screen genotypes for drought tolerance, often coupled with feed quality analyses. Phenotypic assessments based on adaptive morphophysiological-including yield-related and feed quality-traits were the key indicators to identify drought-tolerant, productive and higher-quality forage genotypes [16,17]. Napier grass is one of the best adapted tropical forage species that can withstand drought stress [18,19]. In this study, genotypes from the ILRI and EMBRAPA collections were used to identify genotypes better adapted to produce under soil moisture stress conditions. The present study examined the performances of 84 Napier grass genotypes for morphophysiological, agronomic and feed quality traits under moderate and severe water stress conditions in the dry season, as well as under rainfed conditions in the wet season. In our study, the field drought stresses (MWS and SWS) in the dry season altered the performance of morpho-agronomic and feed quality traits in the crop compared to plants grown during the rainy season. Both MWS and SWS conditions negatively influenced the plants' morphological development by decreasing PH, LL and LW in all of the studied Napier grass genotypes, since optimum soil moisture is important for nutrient uptake and plant growth and development [20]. However, the number of tillers did not decrease under drought stress treatments. This may be attributed to the adaptive nature of tiller density, which indicates increased resource use efficiency under stress conditions, enhancing the chances of survival and fast recovery [21,22]. Such performance changes have also been noted in previous studies that reported on the morphological and physiological adjustments of Napier grass plants when exposed to drought stress conditions [19]. The reduction in morphological development could be associated with a decline in cell growth and expansion that limits the overall plant architecture when growing under drought stress. Similarly, drought stress affected the photosynthetic efficiency (Fv/Fm and PI) and forage biomass yield (TFW and TDW) compared to growth in the rainy season. Thus, the reduction in photosynthesis may decrease assimilate accumulation, which ultimately determines biomass yield. Taken together, the observed changes in growth performance are related to differences in soil moisture levels between the rainy and dry seasons [13]. This underscores the fact that soil moisture is one of the driving factors in determining the overall growth and development of Napier grass.
Generally, highly significant genotypic differences were observed for all the studied morpho-agronomic and feed quality traits in both rainy and dry season conditions, indicating that the tested genotypes had inherent genetic differences for the studied traits. Genotype × harvest interaction was also significant for TDW, indicating that TDW was quantitatively inherited and differentially expressed in response to different harvest rounds.

Trait Expression under Rainfed, Moderate and Severe Water Stress Conditions
In the dry season, there was significant variation between drought stress treatments (MWS and SWS) for most agro-morphological and feed quality traits (Table 1). Genotypeby-drought-stress-treatment interactions were also significant for most traits, indicating that genotypes performed differently under MWS and SWS conditions. Thus, independent selection for MWS or SWS conditions can be used as a strategy to identify the best performing genotypes, to exploit their potential and to subsequently maximize forage production in the dry season with minimal irrigation support.
In the rainy season, genotypes showed significant genotypic variation for morphological, agronomic and feed quality traits (Tables 1 and 4). These phenotypic variations between genotypes are indicative of the opportunity to select genotypes for increased performance under optimum soil moisture conditions. However, under the wet season condition, there was no statistically significant difference for all morpho-agronomic traits between the blocks that were treated with either MWS or SWS in the dry season. For example, in the wet season, we detected no significant difference in total dry weight production between genotypes exposed to MWS and SWS conditions in the dry season. However, the reduction in biomass yield was higher under SWS than it was under MWS, showing that severe water stress had reduced the biomass yield. Thus, the similar annual above-ground biomass production could be attributed to an enhanced growth rate stimulated by the increased assimilation rate of genotypes grown under SWS. The re-allocation of stored carbohydrates in the root may have contributed to the stronger compensatory growth of plants exposed to severe water stress [23]. This indicates that genotypes which were grown under the SWS condition in the dry season had a greater recuperative ability during the rainy season when the soil moisture became optimum. Hence, such recovery potential is an important response to determine the overall productivity performance of Napier grass genotypes for tropical environments with a long dry season and a short rainy season.
Furthermore, the feed quality traits in the rainy season were not affected by the preceding drought stress treatments (MWS and SWS) in the dry season (Table 3). Since feed quality traits are strongly associated with plant development, the similarities displayed in the morphological and agronomic performance of the plants between treatments might have also resulted in similarities in feed quality status. The results also indicated that significant variation existed between harvests for the quality traits considered, which could predominantly be influenced by the soil moisture conditions due to the differences in the received precipitation during growth periods in the wet season (Table A6).
The present study revealed different levels of variation and heritability among genotypes for all traits measured. Higher levels of variation and heritability were observed during the rainy season compared to the dry season for TN, TFW and TDW traits. These traits recorded the highest values for PCV and GCV in all conditions, indicating the presence of high genetic variation. In the drought stress treatments, TDW, WUE and TN had the highest PCV and GCV values, indicating the presence of high genetic variation, while Fv/Fm and LW had low PCV and GCV values, indicating low genetic variation under both MWS and SWS conditions. Thus, the adaptive traits with high amounts of variation under water stress treatments are worth considering for the effective screening of genotypes for performance under drought stress. The observation of higher PCV than the respective GCV values indicates the significant contribution of environmental effects on trait expression. On the other hand, the high heritability estimates for LL, TN, LW and TDW indicate that the trait expression was predominantly governed by additive gene actions, so direct selection might be effective to improve these traits. Our findings concur with previous studies that reported considerable variation in plant growth traits in perennial forage grasses [24]. Subsequently highly heritable traits, such as LW and LL, can be efficiently used to support improved yield in Napier grass [25]. Under drought stress treatments, strong and positive associations were observed between biomass yield and plant growth and development traits (PH, LL, LW and TN). These morphological traits can play an important role in adaptation under soil moisture stress conditions; thus, the indirect selection for biomass productivity based on these morphological traits could be effective [26]. Plant species usually implement multiple drought resistance strategies, varying morphological and physiological traits in response to drought stress [27]. The present study demonstrated that, under both levels of soil moisture stress, the most productive genotypes generally have superior plant height, leaf length and leaf width and increased tillering compared to less productive genotypes. Therefore, as an adaptation strategy, we suggest that the drought-tolerant genotypes have the ability to enhance post-stress regrowth [28].
In contrast to the agro-morphological traits, the level of variation and heritability was generally low for feed quality traits. In all conditions, relatively high amounts of variation and heritability were observed for the traits CP, ADL and ADF. Generally, PCV values were more than two-fold the respective GCV values, indicating that environmental effects strongly influenced trait expression.

Drought Stress-Responsive Agro-Morphological Traits
A PCA analysis was conducted to explain and describe important indicators of stress resistance in the plants that can be utilized for the selection of drought-tolerant genotypes. The PCA indicated that, under drought stress treatments, STI, WUE, PH, LL and LW significantly influenced biomass yield (Table A3). In both drought stress treatments, the PCA biplot revealed that most of the variation was contributed by PC1, where WUE and STI accounted for the highest positive contribution, whereas Fv/Fm and PI contributed negatively to the variation. Thus, the major positive contributing traits can act as a promising indicator for screening Napier grass genotypes growing under moisture stress conditions. Previous studies revealed that WUE is an important characteristic of a plant's response to drought stress, where a high WUE indicates adaptability to a drought environment [27]. The biplot created between PC1 and PC2 showed the grouping of genotypes along the vector line. Genotypes 16839 and 16801 were located at a considerable distance along the line from the origin and hence can be considered drought-tolerant genotypes. The PCA biplot also showed that genotypes 16839, 16791 and 16819 were located close to the vector line that is associated with drought-tolerant indicators (WUE and STI); thus, these traits could be employed for identifying groups of Napier grass genotypes with drought stress tolerance and susceptibility. These findings indicate that the maintenance of high WUE and STI makes the greatest contribution to biomass yield compared with other morphological traits.

Drought Stress Effects on Genotype Performance
The evaluation of genotypes growing at different drought stress severity levels is useful to support the selection of high-yielding and stable genotypes for production in drought-prone environments. The current study indicated that the SWS condition resulted in a greater reduction in the growth and development of most genotypes compared to the MWS treatment. Plants utilize various adaptation mechanisms to cope with different drought stresses and ensure growth and development. The present study revealed that Napier grass genotypes are able to survive under reduced soil moisture conditions, as no genotype perished under the SWS condition, although growth and development were negatively affected. The observed genotypic variation for WUE and STI under both MWS and SWS conditions highlights the potential to screen Napier grass genotypes for enhanced drought stress tolerance. The genotypic variation in WUE could be linked to a reduction in stomatal density and leaf width and length so that the plant maintains its internal water balance [29] or increased water status due to the increased root depth and biomass under drought stress [3,4]. These suggestions may be the subject of future research.
A PCA identified that enhanced WUE, STI and TDW were the key indicators for the analysis of the response of genotypes to drought stress. Accordingly, these traits were able to distinguish drought-tolerant and -susceptible genotypes maintained under both MWS and SWS conditions. In both stress conditions, genotypes were distinctly categorized into highly tolerant and susceptible groups, while some genotypes that were grouped in the highly tolerant group under MWS were grouped in the moderately tolerant group under SWS. Among the genotypes, 16819, 16791, BAGCE 30 and 16839 were considered to be highly drought-tolerant. On the other hand, genotypes such as 16834, 16790, 16797 and 18662 were found to be susceptible. Our findings confirmed a previous report that indicated that genotype 16790 was a low-performing genotype, while genotype 16791 was highly productive [30].

Biomass Yield of Genotypes
The evaluation of the total dry weight production of genotypes across the wet, MWS and SWS conditions revealed that the highly productive genotypes identified under the rainy season were also high yielders under both MWS and SWS conditions. A similar observation was also reported by [13], where the top yielding genotypes on an annual basis also performed best under dry season production conditions. Furthermore, the general linear model regression for TDW between rainy and MWS conditions and between rainy and SWS conditions also substantiated the observation of consistent, productive genotypes across the seasons and drought stress treatments. The highly productive genotypes during the wet season include 16791, 16819, BAGCE 30, CNPGL 93-37-5 and 16802. Our findings align with a previous study that reported that genotypes 16791 and 16819 were high biomass yielders, and these accessions have already been released as commercial cultivars in Ethiopia [21]. Under MWS, genotypes 16819, 16803, 16839, BAGCE 30 and 16811 produced the highest forage yield. Similarly, under SWS, genotypes 16819, CNPGL-93-37-5, 16839, BAGCE 100 and BAGCE 30 were the top forage producers. Therefore, we suggest that the selection of Napier grass genotypes under optimal and water stress conditions is a suitable approach to obtain genotypes that can reliably produce high biomass yields in the dry season and maximum yields in optimum soil moisture conditions. Regarding the annual biomass production, the highest dry biomass yield was obtained from highly productive genotypes identified under wet season and drought stress conditions (Table 5). These top biomass-yielding genotypes were tall [31] with low LSR values, indicating that they had a high stem biomass, although this could, at least in part, be associated with their accelerated development due to their high production. A separate study that agreed with our findings showed that WUE was higher in vigorous genotypes than it was in low-yielding genotypes [30]. The top biomass-yielding genotypes could: (a) avoid drought stress by regulating water loss or improving water uptake, or (b) store water in the stem [32], thereby enhancing WUE [4]. Other reports also suggest that productive genotypes have the capacity to compensate for the reduced soil moisture by the enhanced recovery potential in the rainy season, whereas less productive genotypes are slow-growing in addition to being sensitive to the stress that further limits plant development [33].
The present study revealed significant genotype-by-harvest-round interaction in each of the three growing conditions, as indicated by the crossover performances for genotypes on different harvest rounds using AMMI analysis (Table A4). This might lead to variations in the mean ranks of the genotypes at different harvests. Thus, AMMI stability value analysis was performed to identify genotypes whose forage biomass yield fluctuated less over the seasons. Here, genotypes with low ASV values are scored as highly stable, whilst those with high ASV values are scored as less stable. However, ASV stability alone for biomass yield might not indicate productivity because consistently low-yielding genotypes, such as 16790 and 16621, have been shown to be highly stable. Therefore, the yield performance index (YSI), the sum rank of the mean yield across the time of harvests with the rank of the ASV of genotypes, was employed to determine both stable and productive genotypes. Genotypes with a low YSI are considered for selection, as they combine a high yield and stable performance. The results of the present study showed that high-yielding genotypes such as 16791, 16819, BAGCE 30 and BAGCE 100 had high ASV values and were categorized as moderately to lowly stable genotypes using the yield stability index. A similar observation has been reported by a previous study [21], where high-yielding genotypes, such as 16791, had low stability [34].

Feed Nutrient Quality Performance
The feed quality of forage is an essential component to consider during selection. Significant genotypic variation for feed quality traits was observed in each of the wet, MWS and SWS growing conditions, suggesting the potential to select for these traits. This study revealed the existence of genotypic variation for feed quality traits, similar to previous reports [35,36]. In addition, the growing conditions influenced the feed quality; for example, the mean values of NDF, ADF and ADL were higher during the rainy season, while CP, IVOMD and ME were high in the dry season [37]. The decrease in CP and IVOMD during the rainy season may be associated with a dilution effect due to increased phenological development. The observation of reduced NDF in the dry season was also reported by a previous study [30]. However, in contrast to these studies, another study [21] reported increased NDF in the dry season, although these conflicting observations could be associated with differences in the maturity level at harvest between the two trials [33]. Our findings have revealed that genotypes that produced the highest CP content were consistent across the three growing conditions, indicating the inherent characteristics of the genotypes. As expected, the fiber components (NDF, ADF and ADL) were negatively correlated with CP, IVOMD and Me [37].
The combined biplot analysis of Napier grass genotypes across the two drought stress treatments identified clusters of genotypes with high fiber components (NDF, ADF and ADL) such as 16839, BAGCE 100 and 16819 [34] or high CP and IVOMD, such as 16811, BAGCE 81 and BAGCE 17 ( Figure A2). These results indicate that there is considerable opportunity for the improvement of Napier grass in terms of different forage quality traits. Genotypes that produced a high biomass yield under MWS and SWS also had a higher fiber content, presumably to support increased growth, while their CP content was low. A similar result was also reported in a previous study [36], where high dry matter-producing Napier grass cultivars exhibited low nutritional quality. However, the present study showed that the CP content of the productive genotypes was more than the 7% that is the minimum required level for rumen microbial activity [38]. It is intriguing to note that genotypes that combined both a good feed quality and a high biomass yield, such as 16811 and 14983, were identified.
Leaf-to-stem ratio (LSR) has been shown to be a key component in determining feed quality. In the rainy season, the highest LSR was observed in genotypes 18448, BAGCE 17, 16902 and 16787-genotypes that also had high CP and IVOMD values. The high biomass-yielding genotypes 16791, 16819 and BAGCE 100 had low LSR values, while their fiber components were high. Nevertheless, considering the annual CP yield, these high biomass-yielding genotypes produced the highest CP yield, suggesting that annual biomass production coupled with annual feed nutrition values such as CP yield are important traits as far as forage value is concerned.
Partitioning the forage value into leaf and stem fractions in terms of nutritional composition is a key breeding strategy to develop higher-quality forages. Generally, the leaf is more digestible, higher in CP and lower in fiber components than the stems [39]. Such differences in feed nutrition between the leaf and the stem might be attributed to the photosynthesis/metabolic role of leaves and the structural role of stems [40]. Consequently, improving feed quality relies, at least partially, on increasing the proportion of the leaf when harvested. Consistent with previous reports, our findings show that the leaves generally have higher CP and IVOMD values and lower fiber components than the stems [36,41]. However, it is interesting to note that the stem tissues of some Napier grass genotypes have high CP and IVOMD values compared to the leaves. This may provide a biochemical basis for the selection of Napier grass with increased feed value.

Description of the Experimental Area and Planting Material
The study was conducted in Bishoftu, Ethiopia (008 • 4702000 N and 038 • 5901500 E) from 2018 to 2020. The altitude of Bishoftu is 1890 m.a.s.l., with an Alfisol soil type. A panel of 84 Napier grass genotypes sourced from the ILRI forage gene bank were used for the study (Table A1).

Field Trial Set Up
The trial was set up using a partially replicated (P-rep) design in four blocks. The field planting and growing conditions are as described in a previous report [42]. Land preparation, planting, weeding, harvesting and related management practices were uniformly applied across all plots. Inorganic fertilizer, urea and di ammonium phosphate (DAP) (50:50) were applied at a rate of 6.2 g/plant in each of the dry and wet seasons. The stress treatments, defined as moderate water stress (MWS) with 20% soil volumetric water content and severe water stress (SWS) with 10% soil volumetric water content, were imposed in the dry seasons (November to May), as indicated in our previous report [42] (Figure A3). As indicated in a report by [4], genotypes with drought stress tolerance characteristics also need to possess a high yield potential in non-stress environments to maximize annual yield; hence, genotypes were also evaluated during the rainy season (wet), which runs from June to September. The soil moisture content of the field plots was monitored using a Delta soil moisture probe (HD, UK) ( Figure A5). In addition, a weather station (Speck ware technologies) was installed to monitor the daily climatic variables [4]. A composite soil sample from each replication was collected at the start and end of the experiment using an auger from a depth of up to 40 cm. The physical and chemical characteristics of the analyzed soil in the blocks are presented in Table 6. Organic-carbon-to-nitrogen ratio (C:N), parts per million (ppm).

Data Collection
After the genotypes were well established, the plants were harvested at about 5 cm above ground level after every eight weeks of regrowth, resulting in a total of 12 harvests, conducted between June 2018 and May 2020. In each harvest, morphological, physiological, agronomic and feed quality traits were collected. Three randomly selected plants were used for all trait measurements described below.
Morphological traits: Plant height (PH-in cm), leaf length (LL-in cm), leaf width (LW-in mm), stem thickness (ST-in mm), tiller number (TN) and internode length (IL-in cm). PH was measured from the base to the tip of a randomly selected plant tiller, the same tiller was used to measure ST and IL, LL and LW were measured from the third leaf from the top of the plant.
Physiological traits: Water use efficiency (WUE) was calculated by dividing the total dry weight per plant by the total volume of irrigated water applied per plant during the dry seasons [19]. Chlorophyll fluorescence (Fv/Fm) and performance index (PI) measurements were conducted on fully expanded leaves using a plant efficiency analyzer (Handy PEA; Hansatech Instrument Ltd., Lynn, UK) in the morning hours. The measured leaves were dark-adapted for 20 min before the measurement. The measurement consisted of a single strong 1 s light pulse (3000 L/mol photons m −2 s −1 ).
Feed quality traits: In each harvest, the oven-dried whole plant, leaf and stem biomass fractions were ground separately for feed quality analysis. The feed quality traits for each component were measured as indicated by [13]; the traits measured were: neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), organic matter (OM), crude protein (CP), in vitro organic matter digestibility (IVOMD) and metabolizable energy (Me). The feed quality traits are corrected and expressed on a percentage of dry matter (DM) basis, except for Me, which is expressed as MJ/Kg DM. The annual CP yield (CPY) was derived from the annual CP content and TDW yield by multiplying the annual CP by the annual TDW.

Data Analysis
An averaged data value per trait per genotype was generated in each of the three soil water regime conditions (MWS, SWS and Wet). There were six harvests for the MWS and SWS treatments, each with two replicates; however, in the wet season, as the soil moisture was similar across all blocks ( Figure A3), the data from all four blocks were used for the six harvests. Hence, twelve harvests were used in this case. The averaged values for all traits were corrected for spatial variation using spatial analysis, assessed as described in a previous study [11,13].
The fitted data were used for the analysis described below. The values of the trait data were checked for normality using the bartlett test for heterogeneity. Statistical analysis was conducted using analysis of variance (ANOVA) in the GenStat software (version 19, VSN international Ltd., Hemel Hempstead, UK) [43] to determine the significance of the main effects and the interactions using the model: where Y ijk is the response during each season, µ = overall mean, G i = effect of the ith genotype, T j = effect of the jth treatment in the dry season (MWS/SWS), H k = effect of the kth harvest period, G i * T ij = the interaction of ith genotype and jth treatment, G i * H ik = the interaction of ith genotype and kth harvest, T j * H jk = the interaction of the jth treatment and kth harvest, G i * T j * H k = three-factor interaction among genotypes, treatments and harvests and ε ijk = the residual error. The least significant difference (LSD), for the comparison of the mean values of the traits, was employed to compare genotypes for traits with significant differences. The genetic parameters, genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV) were estimated using the formulae [44]:

of 34
where GCV = genotypic coefficient of variation, PCV = phenotypic coefficient of variation, σ 2 g = genotypic variance, σ 2 p = phenotypic variance and X = grand mean. The broad-sense heritability (H 2 ) for the traits was captured using the equation [25]: where σ 2 g and σ 2 e are the variance components for the genotype effect and the residual error, respectively. To categorize genotypes in clusters in both the MWS and SWS of the dry season, a hierarchical co-cluster algorithm was employed using the R package 'cluster' version 2.1.2 [45]. The visualization of the cluster analysis was prepared using DeltaGen 3_1 [46].
Principle component analysis (PCA) based on the correlation matrix was performed using the packages ggplot2, factoextra and FactoMiner [47] to identify influential traits for selection. The Eigen values, latent vectors and PCA biplot were extracted from the PCA. PCA biplots were plotted separately for MWS and SWS using DeltaGen 3_1 [46] to show the relationships among the studied genotypes based on recorded traits. Principle component analysis (PCA) based on the correlation matrix was performed using the packages ggplot2, factoextra and FactoMiner [48] to identify influential traits for selection. The Eigen values, latent vectors and PCA biplot were extracted from the PCA. PCA biplots were plotted separately for MWS and SWS using DeltaGen 3_1 [46] to show the relationships among the studied genotypes based on recorded traits.
For the graphical presentation of the mean values of traits, bar and line graphs were used. The relative magnitude of change in the trait values due to the two stress conditions, MWS and SWS, was presented as a radar plot. The graphical presentations of the bar graph, line graph and radar plots were prepared using Microsoft office software.
To select high-yielding genotypes growing under stressed conditions, a stress tolerance index (STI) was calculated using the following formula [47]: where Ys = total dry weight (TDW) of test genotypes growing under MWS or SWS conditions; Yp = total dry weight (TDW) and Xp = mean TDW of genotypes growing under the wet condition. Principal component analysis (PCA) and hierarchical clustering were generated using Manhattan distance to classify the tolerance ranking via the average linkage method, using STI, WUE and TDW for both MWS and SWS conditions separately.
The main additive effects and multiplicative interaction (AMMI) analysis model was applied for MWS, SWS and wet conditions using the AMMI function in the R package Agricole [49] separately for 84 Napier grass genotypes (G) and 6 harvest time points (Environments = E) as additive effects and genotype environment interaction (GEI) as a multiplicative term. The AMMI analysis first fits additive effects for host genotypes and environments by the usual additive ANOVA procedure and then fits multiplicative effects for G × E (genotype × environment) by PCA. (6) where Y ij is the TDW of the ith genotype in the jth environment, µ is the grand mean, G i and E j are the ith genotypic effect and jth environment effect, respectively, λk is the square root of the eigenvalue of the PCA axis k and α i k and γ j k are the principal component scores for the PCA axis k of the ith genotype and the jth environment. The AMMI stability value (ASV) was computed as described in [50] using GenStat V19. Smaller ASV scores indicate a more stable genotype characteristic across environments.
The yield stability index was also calculated using the sum of the ranking based on yield and ranking based on the ASV value.
where RASV is the rank of the genotypes based on the ASV value; RY is the rank of the genotypes based on TDW across environments/harvests (RY). YSI incorporates both the mean yield and stability in a single criterion. Low values of both parameters show desirable genotypes with a high mean yield and stability.
To describe the magnitude of the relationships among agro-morphological and feed quality traits, Pearson's correlation coefficients (r) were calculated separately for the dry and wet season conditions using GenStat (version 19, VSN international Ltd., Hemel Hempstead, UK) [43].
To evaluate the genotypes based on annual TDW and CP, cumulative TDW and CP yields across harvests were computed. The cumulative annual TDW was obtained by adding the mean TDW of the wet and dry season harvests for each genotype. The cumulative CP yield was obtained by adding the mean CP of wet and dry season harvests and multiplying by the respective cumulative TDW.

Conclusions
The field evaluation of 84 Napier grass genotypes sourced from the ILRI gene bank using moderate and severe soil moisture stress treatments elucidated significant variation in terms of agro-morphological, physiological and forage quality traits. The observed variations in phenotypic performance and the significant difference among genotypes in terms of measured traits indicated the possibility of developing varieties with improved water use efficiency and adaptability to drought stress environments. In general, the growth and development of most genotypes declined under SWS compared to MWS conditions. However, the genotypes that performed well in terms of WUE and TDW showed some level of consistency in both MWS and SWS conditions. For example, under MWS, genotypes 16819, 16803, 16839, BAGCE 30 and 16811 produced the highest forage yield; similarly, under SWS, the most productive genotypes were 16819, CNPGL-37-5, 16839, BAGCE 100 and BAGCE 30. In terms of forage quality traits such as crude protein content, high biomass-producing genotypes had low CP per kg of dry matter, but the total CP yield (CPY) of the high biomass-producing genotypes was also high, indicating the possibility of developing high-feed-quality Napier grass varieties under drought stress environments.
Supplementary Materials: The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/plants11192549/s1, Table S1: Mean Napier grass genotype values for agromorphological traits across harvests for wet, MWS and SWS conditions, Table S2: Mean Napier grass genotype values for feed quality traits across harvests for wet, MWS and SWS conditions, Table S3: Mean Napier grass genotype values for leaf and stem tissue samples for feed quality traits under the wet growing condition.

Acknowledgments:
The authors would like to thank Asebe Abdena, Tesfaye Tadesse, Yirsaw Wubete and Fetene Argaw for their contribution to the field data collection. The authors would like to thank Yonas Asmare and Prasad KVSV for the NIRS analysis and data prediction. The authors would also like to thank EMBRAPA for making their germplasm and breeding lines available for the study.

Conflicts of Interest:
The authors declare no conflict of interest.   Table A2. ANOVA summary of agro-morphological traits of 84 Napier grass genotypes in wet, moderate water stress (MWS) and severe water stress (SWS) conditions.          Neutral detergent fiber (NDF); acid detergent fiber (ADF); acidic detergent lignin (ADL); organic matter (OM); crude protein (CP); in vitro organic matter digestibility (IVOMD); metabolizable energy (ME); coefficient of variation (CV). Values indicate the significance probability level; no significance (NS). Neutral detergent fiber (NDF); acid detergent fiber (ADF); acidic detergent lignin (ADL); organic matter (OM); crude protein (CP); in vitro organic matter digestibility (IVOMD); metabolizable energy (ME); coefficient of variation (CV). Values indicate the significance probability level; no significance (NS). Figure A2. PCA biplots of feed quality traits and total dry weight (TDW) of Napier grass genotypes growing under: (A) moderate water stress conditions and (B) severe water stress conditions. The color of genotypes indicates different cluster groups using pattern analysis. The angles between vectors indicate the degree of correlation between traits. Neutral detergent fiber (NDF); acid detergent fiber (ADF); acid detergent lignin (ADL); organic matter (OM); crude protein (CP); in vitro organic matter digestibility (IVOMD); metabolizable energy (ME).

Traits/Sources of Variation
Tissue Samples  Table A7. ANOVA summary for seven feed quality traits of Napier grass genotypes for leaf and stem tissue in the wet season.