Next Article in Journal
Enhancement of Morphological and Physiological Performance of Zea mays L. under Saline Stress Using ZnO Nanoparticles and 24-Epibrassinolide Seed Priming
Previous Article in Journal
Evolution of Maize Compost in a Mediterranean Agricultural Soil: Implications for Carbon Sequestration
Previous Article in Special Issue
Seed Phenotyping and Genetic Diversity Assessment of Cowpea (V. unguiculata) Germplasm Collection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Forage Quality in Alfalfa (Medicago sativa L.) via Selection for Increased Stem Fiber Digestibility

USDA, Agricultural Research Service, Plant Science Research Unit, 1991 Upper Buford Circle, St. Paul, MN 55108, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Agronomy 2023, 13(3), 770; https://doi.org/10.3390/agronomy13030770
Submission received: 3 February 2023 / Accepted: 6 March 2023 / Published: 7 March 2023
(This article belongs to the Special Issue Breeding Advances in Legume Diversification and Biofortification)

Abstract

:
The low digestibility of fiber in alfalfa (Medicago sativa L.) limits dry matter intake and energy availability in ruminant animal production systems. Previously, alfalfa plants were identified for low or high rapid (16 h) and low or high potential (96 h) in vitro neutral detergent fiber digestibility (IVNDFD) of plant stems. Here, two cycles of bidirectional selection for 16 h and 96 h IVNDFD were carried out. The resulting populations were evaluated for total herbage, percentage of stems to total biomass, IVNDFD, neutral detergent fiber (NDF), and acid detergent lignin as a proportion of NDF (ADL/NDF) at three maturity stages. Within these populations, 96 h IVNDFD was highly heritable (h2 = 0.71), while 16 h IVNDFD had lower heritability (h2 = 0.46). Selection for high IVNDFD reduced NDF and ADL/NDF in plant stems at the late flowering and green pod maturity stages and reduced seasonal variability in stem digestibility but did not alter the percentage of stems. Stability analyses across 12 harvest environments found that selection for high IVNDFD had little effect on environmental stability of the trait compared to the unselected population. Thus, selection for stem IVNDFD was a highly effective strategy for developing alfalfa populations with improved nutritional quality without changing the percentage of stems to total biomass.

1. Introduction

In the United States, alfalfa (Medicago sativa L.) is the third largest field crop produced, with an annual estimated value of over USD 10 billion [1]. Primarily used in animal feeds, it is also valued for its diverse environmental services [2] such as carbon sequestration and mitigation of nitrogen leaching from surface water and tile drainage. Alfalfa herbage is an excellent source of nutrients for milk production in dairy cattle and increasing muscle mass in beef cattle [3,4,5].
Alfalfa herbage is comprised of a protein-rich, highly digestible leaf fraction and a less digestible carbohydrate-rich stem fraction [6]. The stem fraction is higher in indigestible fiber concentration than the leaf fraction, and as the stems mature, the fiber concentration increases within stems while the fiber in leaves remains low and fairly constant [7]. Because stems make up approximately 46 to 60% of the dry matter in alfalfa herbage [8], increasing the digestibility of alfalfa stems would impact the dairy market substantially through increased animal productivity, reduced feed costs, and reduced animal waste [4,9]. Additionally, increasing the available energy of the stem at later maturity stages would also increase dry matter yield [10,11]. Thus, it has been recognized that alfalfa breeding should focus on improving stem digestibility to improve the overall quality of alfalfa forage [12].
In forage quality analysis, plant cell walls (CWs) are typically characterized by the detergent fiber fractions: neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL). Generally, forages that limit feed intake have greater NDF concentrations and more highly lignified CWs [13]. Because alfalfa stems contain more NDF than alfalfa leaves, morphological changes, such as the leaf-to-stem biomass ratio, usually affect NDF values and digestibility. After alfalfa matures beyond the vegetative stage, there is a rapid increase in NDF concentration and a concurrent rapid decrease in in vitro digestible dry matter [14]. Feeding studies have reported that improving NDF digestibility (NDFD) resulted in significantly greater dry matter intake and milk yield in dairy cows [5]. Therefore, focusing breeding efforts on in vitro neutral detergent fiber digestibility (IVNDFD) may increase alfalfa digestible biomass yield.
Progress in developing alfalfa cultivars with improved stem CW digestibility has been slow due to the quantitative nature of the trait and low heritability [15], even though extensive variability in CW components and IVNDFD occurs in alfalfa germplasm [16]. Further genetic improvement in alfalfa forage digestibility utilizing marker-assisted selection or recombinant DNA technologies will require the identification of biochemical traits and associated genes that impact the rate and/or the extent of CW digestibility.
Fiber digestibility in alfalfa stem CWs has been increased through conventional breeding [17]. Recurrent phenotypic selection was used to develop populations that were evaluated using enzyme-released glucose as a proxy trait for fiber digestibility. The results indicated that digestibility was increased significantly, and the selection method was effective [17]. This trait was correlated with IVNDFD and in vitro total digestibility (IVTD) and had the benefit of not depending on the cyclic and variable rumen composition of different animals. In another breeding program, alfalfa plants were identified with stems that had either low or high rapid (16 h) and low or high potential (96 h) IVNDFD. These plants also differed in stem NDF and ADL as a proportion of NDF (ADL/NDF) [18]. A stability analysis demonstrated that stem NDF and 96 h IVNDFD were more environmentally stable among the alfalfa clones in the field environments than stem 16 h IVNDFD or ADL/NDF, suggesting that stem 96 h IVNDFD, NDF, and ADL as a proportion of NDF would be the traits to use in a field selection program to develop alfalfa with improved nutrient value [10].
In this study, two cycles of bidirectional selection for 16 h and 96 h IVNDFD were carried out on alfalfa stems. This research was designed to evaluate the effectiveness of recurrent selection on stem fiber digestibility for improving the nutritional value of alfalfa. Specifically, the study aimed to measure the genetic gain, heritability, and environmental stability of IVNDFD and determine whether selection altered CW fiber fractions important in forage digestibility. In a subsequent paper, data will be presented on the changes in CW composition associated with the selection for IVNDFD (Heuschele et al., in preparation).

2. Materials and Methods

2.1. Development of Selected Populations

The cycle 0 (C0) parental population, UMN3097 (Table 1), was created by mixing seeds from six commercial alfalfa varieties (5312, Rushmore, Magnagraze, Wintergreen, Windstar, and WL 325HQ) that are adapted to the Minnesota environment [10,18]. As described previously, this population was established at the University of Minnesota Sand Plains Research Farm, Becker, MN [18]. Briefly, in April 1997, approximately 2400 seeds were hand sown to develop a spaced plant nursery with 15 cm within and between the rows. Herbage from individual plants was harvested when approximately 25% of the plants bloomed in spring (June) and early summer (July) in 1998 and 1999. After drying at 60 °C, stems and leaves were hand separated, and stems were ground to pass a 1 mm screen in a cyclone-type mill (SPEX, 8000M, SPEX SamplePrep, Mentuchen, NJ, USA) for further analysis. Each ground stem sample was scanned via near-infrared spectroscopy (NIRS) using a Foss Model 6500 (Foss North America Inc., Eden Prairie, MN, USA) to predict 16 h and 96 h IVNDFD using prediction equations developed previously [18].
Selected plants were divided into four populations from C0: high 16 h and 96 h digestibility (H16 × H96); high 16 h and low 96 h digestibility (H16 × L96); low 16 h and low 96 h digestibility (L16 × L96); and low 16 h and high 96 h digestibility (L16 × H96). The selected plants within each population were removed from the field and grown in the greenhouse for intermating. Plants were intercrossed by hand tripping flowers without emasculation. A Cycle 1 (C1) population for increased rate and extent of fiber digestion, UMN3355 (Table 1), was developed via random intercrossing of the plants identified with the highest 16 h IVNDFD and highest 96 h IVNDFD (H16 × H96) utilizing a total of 117 parents. Similarly, a divergent population, UMN3358, was developed via random intercrossing of plants with the lowest 16 h IVNDFD and 96 h IVNDFD (L16 × L96) using 33 parents. The number of intercrossing parent lines was determined based on availability and is listed in Table 1. To test gene effects, Population UMN3356 was developed by intermating plants with high 16 h IVNDFD and low 96 h IVNDFD (H16 × L96), and Population UMN3357 was developed by intermating plants with low 16 h IVNDFD and high 96 h IVNDFD (L16 × H96). For each population, equal amounts of seed were bulked from each plant for use in the next selection cycle. Seeds from the four C1 populations were established in a spaced plant nursery in 2001, as described above for the parental population. As with the first cycle, herbage was harvested three times a year over two years (2002 to 2003) from individual plants, was separated into leaf and stem fractions, and the stem fraction was analyzed via NIRS for 16 h and 96 h IVNDFD. Selected plants were removed from the field for intercrossing in the greenhouse as described for C1 to develop Cycle 2 (C2) populations UMN4016 (H16 × H96), UMN4017 (H16 × L96), UMN4018 (L16 × H96), and UMN4019 (L16 × L96) (Table 1).

2.2. Evaluation of Selected Populations

Experimental plots were established in 2009 at the Sand Plains Research Farm, Becker, MN, and the University of Minnesota Experimental Research Station, St. Paul, MN, with the parental population C0, and C1 and C2 populations. The plots at Becker were irrigated to meet plant moisture needs using the checkbook method [19], while plots in St. Paul were rainfed. The study layout was a randomized complete block design (RCBD) with four replicates. Plots were 1.2 × 1.2 m with plants seeded at 7.6 cm spacings resulting in 16 plants per 1.2 m2. A border of the alfalfa cultivar Agate was planted around each side of the plots to prevent edge effects. The experiment layout was a rectangle of 36 ranges by four rows in both Becker and St. Paul, Minnesota. Plots were managed with best practices, and pesticides were applied as needed. Fertilizers were used based on the recommended optimal level of P and K for high alfalfa yields [20].
Plots were sub-sampled in 2010 and 2011 to collect herbage at the early bud (EB) stage in which one or two nodes had visible buds with no flowers or seed pods present, at the late flowering (LF) stage when 25–30% of plants had open flowers, and at the green pod (GP) maturity stage in which a plant had more than four nodes with green seed pods. For each plot, the number of stems was counted, and the average number of stems per plant was calculated. All plants were cut at a 5 cm stubble height at each harvest, dried at 60 °C, and stems were hand-separated from the leaf material. Data on total herbage dry weight, leaf dry weight, stem dry weight, and percent stem were collected. Stems were then ground to pass a 1 mm screen in a cyclone-type mill and analyzed via NIRS for 16 h and 96 h IVNDFD, NDF, and ADL as a proportion of NDF [18].
The NIRS prediction equations were created using the software program Calibrate (NIRS 3 version 4.0, Infrasoft International) with the modified partial least squares regression option. A set of 138 samples, selected via stratified random sampling to provide an equal representation of each germplasm source, year, and harvest of the study, was used to calibrate the NIRS stem equations (Supplemental Table S1). The 138 alfalfa calibration stem samples were assayed for NDF, ADF, and ADL in sequential order [21] using the Ankom filter bag method and sulfuric acid for ADL determination [22]. The ADL was adjusted using NDF following the approach by Jung and Lamb [18]. The IVNDFD was determined as described previously by Jung and Lamb [22]. Rumen fluid was collected approximately 3 h post-feeding from a fistulated, lactating Holstein cow (Bos taurus L.). Samples were placed in Ankom filter bags (Ankom Technology Corp., Fairport, NY) and then incubated with a 20:80 (v/v) mixture of rumen fluid and McDougall’s buffer [23] at 39 °C for 16 h and 96 h in an Ankom Daisy oven. After the incubation, residues within the filter bags were extracted with the neutral detergent solution as described for dietary fibers [24] to calculate IVNDFD. All assays were conducted in duplicate. All filter bag assays contained an assay blank composed of glass fiber to correct for indigestible NDF entering the bags from the rumen fluid inoculum.

2.3. Statistical Analyses

A five-way analysis of variance (ANOVA) was conducted to estimate the interaction relationships between year (Y), location (L), genotype (G), maturity (M), and harvest (H) via function “aov” from the R package “stats” [25]. The R code for the model with all five main effects and all of the interactions was aov (digestibility ~ Y × L × M × H × G, data = digestibility.data). Fisher’s protected LSD test was applied to conduct the pairwise analysis via a function of LSD.test [26] from R package “agricolae” with alpha = 0.05.
Pearson correlation was used to estimate the relationship between the digestibility traits. Functions “cor” and “corrplot.mixed” from R package “corrplot” [27] were employed to estimate the correlation coefficients.
The function “slope.test” from the R package “smart” was used to test whether the regression of the genetic gain was significantly different from 0 using the ordinary least squares method [28]. The rate of genetic gain for all CW traits and IVNDFD was estimated using the best linear unbiased predictions (BLUPs) obtained from the linear mixed model as follows:
Genetic gain rate (GGR) = (BLUPs_C2 − BLUPs_C0)/BLUPs_C0 × 100
BLUPs_C2 are values of the variables from C2 and BLUPs_C0 are the BLUP values of the variables from C0.
Broad sense heritability was calculated for CW traits and IVNDFD as follows:
h2 = Vg/Vp = Vg/(Vg + Ve/r)
Vg stands for genotypic variance, Vp stands for phenotypic variance, Ve is the variance of the residuals, and r is the number of replicates for each measurement [29]. Broad sense heritability was calculated with data from both years and locations.
Linear regression was used to estimate whether genetic gain was significantly larger than 0. The significance of the slope was tested with function “slope.test” of an R package “smatr” [28].
slope.test(y = digestibility, x = cycle.number, alpha = 0.05, method = “OLS”, robust = True, test.value = 0)
Where alpha is the desired confidence level for the 100(1 − alpha)% confidence interval for the common slope; using the ordinary least square method, Robust = True shows that a robust regression was performed with all of the data and treats all of the data as true data compared with outliers. The “Test.value” = 0 is to test the null hypothesis of whether the slope is statistically equal to 0.
Digestibility variables were classified into three categories: Best: p-value of the slope <0.001 and genetic gain rate >0.5., Ok: p-value of the slope <0.05 and genetic gain rate <0.5, or No: p-value of the slope >0.1.
All figures were generated using the R package “ggplot2” [30] and “ggpubr” [31].

2.4. Environmental Stability Analysis

Stability analysis evaluates the environmental stability of a genotype relative to the overall experimental mean of all genotypes [32]. The environmental index for 16 h and 96 h IVNDFD, detergent fiber components, and CW compositional traits was calculated by subtracting the mean value for each trait over each of the 12 harvest environments (two locations, three harvests per location from 2 years’ experiments) from the overall experimental mean. Linear regression analysis was used to compare the performance of an individual population across the 12 harvest environments to the calculated overall environmental index. The stability analysis was conducted with a function of genotype by environment regression “ge_reg” of R package “metan” for multi-environment trial analysis [32]. The slope of regression (b) was used to evaluate population environmental stability for each cycle of selection using the Type II stability definition of Lamb et al. [10]. A stability coefficient b = 1 indicated that the response of the population to different environments was parallel to the average response of all other populations and was stable across environments. If b > 1, the response of the population was greater than the average response to favorable environments, or conversely, the response was less than average in unfavorable environments. If b < 1, the response was greater than the average in unfavorable environments or less than average in response to favorable environments.

3. Results

3.1. Response to Selection for IVNDFD

The main effects of the year (Y), location (L), genotype (G), plant maturity (M), and harvest (H) were highly significant (p < 0.001) for stem 16 h IVNDFD (Table 2), indicating that all five factors affect stem digestibility. However, among the 15 two-to-five-way interaction terms for 16 h IVNDFD, only 3 of the 15 interactions, H:G, M:G, and L:M:G, had significant interactions (p-values < 0.01). There were no significant interactions between genotype by year or between genotype by location. There were no complex interactions among the four and five factors with genotypes. The most significant interaction was among maturity, harvest, and year (p = 2.80 × 10−208), with the second-largest interaction occurring between maturity and harvest (p = 2.50 × 10−162) (Supplemental Table S2), indicating that maturity and harvest time are the most critical factors for stem digestibility. While the interaction with maturity is expected, the interaction of harvest suggests that the time of sampling is a critical factor for the selection of plants for 16 h IVNDFD. The lack of interaction of genotype with location and year indicates the environmental stability of selected traits and will be addressed more fully later. Across locations, harvests, and maturities, stem 16 h IVNDFD was strongly negatively correlated with NDF (−0.94) (Figure 1).
In contrast to the limited interactions observed with 16 h IVNDFD, significant interactions were found for 96 h IVNDFD (Supplemental Table S2) between germplasm and year (p = 0.0059, p < 0.01, respectively) and germplasm and location (p = 0.02, p < 0.05, respectively) as well as between germplasm with year and location (p = 0.00012, p < 0.001, respectively). Interactions between year and location suggest that the environment plays a significant role in determining the extent of fiber digestibility when digestion time increases from 16 to 96 h. Stem 96 h IVNDFD was negatively correlated with NDF (−0.73), but the correlation strength decreased compared with 16 h IVNDFD from −0.94 to −0.73 for NDF (Figure 1).
Interestingly, the percentage of stems (Per_stem) was the most stable trait with no interactions among any of the five factors except for one significant interaction among location, maturity, and genotype with a p-value = 0.047. The observation of low interactions for the percentage of stems indicates that environmental factors did not affect the percentage of stems. Stem dry weight (SDW) was significantly affected by maturity and harvest, and stem dry weight significantly interacted with germplasm (Table 2). There were significant negative correlations of stem dry weight with 16 h IVNDFD (−0.47).

3.2. Genetic Gain by Cycle for IVNDFD and Fiber Digestibility

Stem 16 h and 96 h IVNDFD increased with each cycle of selection in the H16 × H96 populations and decreased in the L16 × L96 populations (Figure 2). A faster gain was seen for 96 h IVNDFD (b = 9.8) compared to 16 h IVNDFD (b = 4.0) from the H16 × H96 population (Figure 2). In contrast, the L16 × L96 populations had negative trait improvement gains, with negative improvement rates (b = −3.7 and −15.0 for 16 h IVNDFD and 96 h IVNDFD, respectively). NDF and ADL as a proportion of NDF increased with selection in the L × L populations and decreased in the H16 × H96 populations compared to the parental population. In the H16 × H96 populations, stem dry weight decreased with selection; in the L16 × L96 populations, stem dry weight increased with selection. Percent stem in the H × H and L × L populations did not change across the selection cycles (Supplemental Table S3).

3.3. Heritability and Genetic Gain of IVNDFD and Detergent Fiber Components

The values for broad sense mean heritability of fiber digestibility (16 h and 96 h IVNDFD) were significantly greater than the statistical threshold of 0.15 (Table 3). The high heritability of 96 h IVNDFD is likely linked to changes in ADL/NDF, which was highly negatively correlated with 96 h IVNDFD (ADL/NDF −0.82, p < 0.001) (Figure 1). ADL/NDF also had the highest heritability (0.75) and high rates of genetic gain (−3.98) (Table 3). Selection for IVNDFD did not result in significant changes in stem dry weight or the percentage of stems (p-values > 0.05).

3.4. Effects of Plant Maturity and Harvest on IVNDFD

The stem digestibility of the 16H × 96H populations was measured at three stages, early bud (EB), late flowering (LF), and green pod (GP), over three harvests each year to determine how selection affected digestibility and fiber traits during plant development. The mean 96 h IVNDFD increased with each selection cycle, and the increase was greater at the late flower and green pod stages (Table 4). There was no significant change in 96 h IVNDFD at the early bud stage among the selection cycles. At the late flowering and green pod stages, the second harvests showed the greatest differences among the selection cycles (bold fonts). Overall, the increase in 96 h IVNDFD was greater for C2 than for C1 across the three maturities and harvests (Table 4).

3.5. Additive Gene Effects May Control Alfalfa Stem Digestibility

The patterns and trends in 96 h IVNDFD for the four different populations suggest that additive genes may control this trait. In the H × H populations, 96 h IVNDFD increased with each cycle of selection at the three maturity stages, EB, LF, and GP, suggesting the accumulation of positive alleles from C0 to C2 (Figure 3A), and this trend was similar in the first, second, and third harvests (Figure 3B). In contrast, 96 h IVNDFD for the L × L populations decreased with each cycle of selection, suggesting the stacking of negative alleles. There were no unidirectional increases or decreases in the H × L and L × H populations.

3.6. Stability of IVNDFD with Cycles of Selection

Complex genotype × environment interactions for forage quality traits in alfalfa present a challenge to breeders in analyzing germplasm performance. Stem IVNDFD, ADL/NDF, and NDF were shown in this experiment and previous experiments [10] to be impacted by the environment. Environmental stability was measured for each of these traits across the 12 harvest environments for the H16 × H96 and L16 × L96 populations. Stability coefficients for 96 h IVNDFD decreased with selection cycles from 1.08 (C0) to 0.93 (C2) for the H × H populations (Figure 4). In contrast, stability coefficients increased in the L16 × L96 populations with selection cycles from 0.93 to 1.04 and from C0 to C2, respectively. The stability of 16 h IVNDFD followed similar patterns (Supplemental Figure S1). A stability coefficient of b > 1 indicates that the response of the population was better than the average response in favorable environments, while b < 1 suggests that the population performed better than average in unfavorable environments. Therefore, for 96 h IVNDFD, the C2 H16 × H96 population performed better in the less favorable environments, and the C2 L16 × L96 population performed better than average in the favorable environments. Overall, 16 h and 96 h IVNDFD were considered to be environmentally stable in H × H and L × L populations.

4. Discussion

Recurrent selection for alfalfa stem IVNDFD was a successful strategy for reducing ADL as a proportion of NDF and NDF in stems without changing the proportion of leaves to stems in total herbage. Increased digestibility and decreased lignin occurred in the later maturity stages, late flowering and green pod, which would increase biomass yields while maintaining forage quality. The selection also reduced seasonal variability in stem digestibility.
Alfalfa leaves are more digestible than stems; therefore, increases in the proportion of leaves to stems result in significant changes in herbage quality [8]. A higher leaf-to-stem ratio is also positively correlated with forage intake [33]. When decreased total herbage lignin concentration was used as a selection trait, overall forage quality increased due to a larger fraction of leaves but the digestibility of the fiber fraction remained unchanged [33,34]. The current study focused on improving stem fiber digestibility to increase forage quality. Using this approach, we found that the percent stem did not change significantly with the increase in stem IVNDFD over each cycle of selection (Table 3). The regression slope for the percentage of stems was close to 0 (p = 0.26) from C0 to C2 (Table 3). The percentage of stems was the most stable trait with the least G × E interaction, which further indicated that the increase in digestibility was due to the improved fiber digestibility of the stems, rather than the proportion of leaves, which changes with plant maturity.
Alfalfa cultivars with improved total herbage digestibility have been developed using traditional breeding methods and genetic modification. The improved total herbage digestibility for conventionally bred Hi-Gest cultivars resulted from an increase in the leaf-to-stem ratio [35,36], whereas the improved forage nutritive value of genetically modified HarvXtra alfalfa was primarily a direct result of changes in ADL and IVNDFD within the stem fraction from the downregulation of caffeoyl coenzyme A-3-O-methyltransferase in the lignin biosynthetic pathway [37,38]. A detailed morphological evaluation of HarvXtra alfalfa in Minnesota reported that HarvXtra and conventional alfalfa cultivars had similar leaf-to-stem ratios that changed similarly with advancing maturity [38]. However, a multistate study found that HarvXtra had a higher leaf-to-stem ratio than the reference alfalfa [39]. The populations developed in the current study appear to be unique, although more research is needed to investigate the basis of increased stem IVNDFD. From observations of CW development and digestion in rumen fluid, Jung and Engles speculated that the rate of alfalfa stem CW digestion could be improved by lowering the syringyl to guaiacyl monolignol ratio in lignin, by reducing the amount of lignification, specifically in the secondary xylem, or by increasing pectin content [40]. Chemical analyses of alfalfa genotypes varying in CW digestibility indicated that the extent of CW digestion could be enhanced by decreasing total CW concentration and selecting for the reduced amount of xylem deposition with maturity [22].
Forage digestibility is measured by the rate and extent of CW breakdown, with rate measured by 16 h IVNDFD and the extent measured by 96 h IVNDFD. When selecting for plants with a high rate and/or extent of stem CW breakdown, the genetic gain and heritability of 16 h IVNDFD was lower compared to the genetic gain and heritability of 96 h IVNDFD. Therefore, breeding progress in increasing the rate of fiber digestibility will be slower than progress in increasing the extent of fiber digestion. The results from this study suggest that more genes may be involved in controlling the rate of digestion than in the extent of digestion. Additive gene control for 96 h IVNDFD was supported by the increase in digestibility in the H16 × H96 populations with each cycle of selection, likely from the accumulation of favorable alleles, while the decrease in digestibility in the L16 × L96 populations was likely from accumulated undesired alleles with each cycle. However, the rates of genetic gain were relatively low, from 4 to 10% in the H16 × H96 populations and −4% to −15% in the L16 × L96 populations. A low rate of genetic gain has been linked to strong genotype by environmental interactions and complex genetic architecture in alfalfa [41,42]. The populations in this study had limited genotype to environment interactions, suggesting that the slow rate of gain was due to either the genetic complexity of the trait or due to the fact that testing locations were not sufficiently diverse to detect interactions.
The heritability of ADL/NDF was high, 0.75, while NDF in the selected populations was moderate, 0.47. These traits are negatively associated with forage digestibility and negatively associated with forage intake [43]. Previously, divergent selection for ADL in total herbage indirectly affected NDF [33,34]. When populations were selected for high ADL, the NDF increased in total herbage and plants had shortened stem internodes compared to the low ADL populations. In the previous study, the improved digestibility occurred due to the dilution of stem NDF by the greater proportion of highly digestible leaves in total herbage. Thus, in this study, the selection for stem IVNDFD improved traits other than leaf digestibility traits.
Most alfalfa is harvested at the early bud stage to maximize forage quality. For ruminant animals, digestibility decreases as plants mature due to the deposition of lignin in vascular cells and the cross-linking of lignin with cellulose microfibrils, as well as the cross-linking of cellulose and hemicellulose. However, harvesting at the early bud stage sacrifices herbage biomass yield potential, increases stress on the plant, and incurs greater financial costs for harvesting since more harvests are required each season.
Harvest timing can also affect alfalfa forage quality. The first harvest, usually in late May in Minnesota, is typical of high quality due to the high leaf-to-stem ratio, but quality often decreases in the hotter and drier harvests obtained in July and August, particularly for later maturity stages. Selection for high IVNDFD resulted in populations that had increased digestibility and reduced ADL/NDF at the second and third harvests (Supplemental Table S4). Selection for low IVNDFD resulted in populations that had reduced digestibility and increased lignin at all three harvests. Thus, the selection methodology resulted in populations with the desired characteristics at later maturities and across harvests, reducing seasonal variation in forage quality.
Improving forage digestibility may impact forage yield. For example, alfalfa from divergent selections for total herbage and low lignin had lower DM yields than high-lignin lines [34]. This was attributed to the reduced main stem length of the low lignin lines. Additionally, genetically modified HarvXtra, with increased stem digestibility, had a lower biomass yield in a multistate study [39], although no biomass yield differences were found compared to conventional lines in a single-state study [37,38]. Despite a trend in yield differences between populations from selection C0 to C2, the only significant decrease in herbage dry weight (HDW) occurred at the green pod maturity stage for the C2 H16 × H96 population (Supplemental Table S5). The yield decreased because, during the second year, the number of plants per plot decreased later in the growing season—especially in the GP plots. This result indicated that an increase in alfalfa stem digestibility was not achieved at the expense of HDW at the early bud and late flowering maturity stages. Similarly, recurrent phenotypic selection for higher stem enzyme-released glucose concentration (i.e., higher digestibility) in mature alfalfa was shown to improve CW digestibility without affecting alfalfa DM yield [17]. Our study relied on hand sampling for determining HDW; thus, the future testing of herbage yields of improved populations in larger plots with mechanical harvesting is warranted.
The pattern of change in 96 h IVNDFD from the four populations shows apparent additive gene effects (Figure 3). The unidirectional increases in the H × H selection or decreases in the L × L selections are evidence of additive effects, while the lack of a pattern in the L × H and H × L populations indirectly confirms the hypothesis of additive effects. However, epistasis cannot be ruled out if the gene interactions are stable for specific germplasm groups as heterotic groups. In addition, stem digestibility is a complex trait, and genome-wide association studies (GWAS) confirmed the polygenic control of forage quality traits [44]. These studies found that many marker-trait associations differed depending on the environment, although a few were consistent across environments [41,45]. Thus, genomic selection may be a promising strategy for continuing to improve forage digestibility [44]. Genomic selection models also show that separating tissue types enhances accuracy, and prediction accuracy is higher for the stem NDFD than for leaf NDFD [44].
Genome-wide markers make it possible to quantify additive and non-additive gene effects. Parametric genomic selection (GS) models are known to capture additive genetic effects but are not efficient with epistatic effects due to the computational burden of high-order interactions [46,47]. Semi-parametric and non-parametric GS models capturing epistatic interactions have been developed and implemented, including high-order interactions in plant breeding [48,49]. Due to the high heritability for 96 h IVNDFD and the potential to capture epistasis and dominance, the use of genomic prediction for stem digestibility with the new non-parametric GS models might further improve breeding efforts to enhance alfalfa forage digestibility. The populations developed in this study will be valuable resources in future experiments to develop GS methods to accelerate breeding for fiber digestibility and to further dissect the morphological and/or compositional changes in alfalfa stems in response to selection.

5. Conclusions

Recurrent selection for alfalfa stem IVNDFD was a successful strategy for improving fiber digestibility and reducing lignin in stems without changing the proportion of leaves to stems in total herbage. Increased digestibility and decreased lignin occurred in later maturity stages, which would increase biomass yields while maintaining forage quality. The selection also reduced seasonal variability in stem digestibility. Successive selection cycles changed the environmental stability of IVNDFD and ADL/NDF, indicating that these traits should be evaluated in multiple environments when developing alfalfa with improved nutritional value. If the mode of action for increased IVNDFD is found to be different than current genetically modified reduced lignin lines, there is a potential for stacking improved stem digestibility traits.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/agronomy13030770/s1: Figure S1: Environmental stability of 16 h IVNDFD with selection cycles across the 12 harvest environments. (Left) H × H, intermating of plants with high 16 h IVDFD and high 96 h IVNDFD. (Right) L × L, intermating of plants with low 16 h IVDFD and low 96 h IVNDFD. Table S1: Calibration statistics for near-infrared reflectance spectroscopy (NIRS) prediction equations for in vitro neutral detergent fiber digestibility (IVNDFD), detergent fiber components, and CW traits of alfalfa stems; Table S2: p-values of the 5-way analysis of variance (ANOVA) for 16 h IVNDFD (NDFD16), 96 h IVNDFD (NDFD96), acid detergent lignin as a proportion of NDF (ADL), neutral detergent fiber (NDF), stem dry weight (SDW), and percent stems (Per_stem) among the following five variables: year (Y), location (L), germplasm (G), maturity (M), and harvest (H); Table S3: Summary of 16 h IVNDFD (NDFD16H), 96 h IVNDFD (NDFD96H), acid detergent lignin as a proportion of NDF (ADL), neutral detergent fiber (NDF), stem dry weight (SDW), and percent stems (Per_stem) of alfalfa populations across years, locations, harvests, and maturities. C0, unselected population; C1, Cycle 1 population; C2, Cycle 2 population; H × H, H16 × H96, intermating of plants with high 16 h IVDFD and high 96 h IVNDFD; L × L, L16 × L96, intermating of plants with low 16 h IVDFD and low 96 h IVNDFD. NDF is on a dry matter (DM) basis. ADL is on an NDF basis. Table S4: Change in ADL in the H16 × H96 populations with cycle of selection at the early bud (EB), late flowering (LF), and green pod (GP) maturity stages over three harvests. LSD is Fisher’s least significant difference test with different letters within a maturity stage indicating a significant difference; %/cycle is the percent change from C0 to C1 and C1 to C2. Table S5: Summary of stem dry weight (SDW) and herbage dry weight (HDW) changes among maturities, harvests, and selection cycles. LSD is Fisher’s least significant difference (LSD) test with different letters within a maturity stage indicating a significant difference; %/cycle is the percent change from C0 to C1 and C1 to C2.

Author Contributions

Conceptualization, J.F.S.L. and H.-J.G.J.; methodology, J.F.S.L. and H.-J.G.J.; formal analysis, Z.X., D.J.H. and D.A.S.; investigation, J.F.S.L. and H.-J.G.J.; resources, J.F.S.L. and H.-J.G.J.; data curation, Z.X.; writing—original draft preparation, Z.X., D.J.H. and D.A.S.; writing—review and editing, J.F.S.L. and H.-J.G.J.; visualization, Z.X.; supervision, D.A.S.; project administration, D.A.S.; funding acquisition, D.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the US Department of Agriculture, Agricultural Research Service, project 5062-12210-004-00D.

Data Availability Statement

All of the data and R codes associated with the analyses are available upon request from the authors.

Acknowledgments

This paper is a joint contribution from the Plant Science Research Unit, USDA-ARS, and the Minnesota Agricultural Experiment Station. We acknowledge Keith Henjum and Ted Jeo for excellent technical assistance in conducting these studies.

Conflicts of Interest

The authors declare no conflict of interest.

EOE Statement

The mention of any trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the US Department of Agriculture. USDA is an equal opportunity provider and employer.

References

  1. Putnam, D.; Meccage, E. Profitable Alfalfa Production Sustains the Environment. Available online: https://calhaysymposium.com/wp-content/uploads/2022/11/PROFITABLE-ALFALFA-PRODUCTION-SUSTAINS-THE-ENVIRONMENT-Dan-Putnam.pdf (accessed on 2 March 2023).
  2. Fernandez, A.; Sheaffer, C.; Tautges, N.; Putnam, D.; Hunter, M. Alfalfa, Wildlife & the Environment; National Alfalfa and Forage Alliance: St. Paul, MN, USA, 2019. [Google Scholar]
  3. De Ondarza, M.; Tricarico, J. Advantages and limitations of dairy efficiency measures and the effects of nutrition and feeding management interventions. Prof. Anim. Sci. 2017, 33, 393–400. [Google Scholar] [CrossRef]
  4. Jung, H.G.; Allen, M.S. Characteristics of plant cell walls affecting intake and digestibility of forages by ruminants. J. Anim. Sci. 1995, 73, 2774–2790. [Google Scholar] [CrossRef] [PubMed]
  5. Oba, M.; Allen, M. Evaluation of the importance of the digestibility of neutral detergent fiber from forage: Effects on dry matter intake and milk yield of dairy cows. J. Dairy Sci. 1999, 82, 589–596. [Google Scholar] [CrossRef] [PubMed]
  6. Wilman, D.; Altimimi, M.A. The in-vitro digestibility and chemical composition of plant parts in white clover, red clover and lucerne during primary growth. J. Sci. Food Agric. 1984, 35, 133–138. [Google Scholar] [CrossRef]
  7. Engels, F.; Jung, H. Alfalfa stem tissues: Cell-wall development and lignification. Ann. Bot. 1998, 82, 561–568. [Google Scholar] [CrossRef] [Green Version]
  8. Sheaffer, C.C.; Martin, N.P.; Lamb, J.F.; Cuomo, G.R.; Jewett, J.G.; Quering, S.R. Leaf and stem properties of alfalfa entries. Agron. J. 2000, 92, 733–739. [Google Scholar] [CrossRef]
  9. Hatfield, R.; Ralph, J.; Grabber, J. Cell wall structural foundations: Molecular basis for improving forage digestibilities. Crop Sci. 1999, 39, 27–37. [Google Scholar] [CrossRef]
  10. Lamb, J.F.S.; Jung, H.J.G.; Samac, D.A. Environmental Variability and/or Stability of Stem Fiber Content and Digestibility in Alfalfa. Crop Sci. 2014, 54, 2854–2863. [Google Scholar] [CrossRef]
  11. Buxton, D.R.; Redfearn, D.D. Plant limitations to fiber digestion and utilization. J. Nutr. 1997, 127, 814S–818S. [Google Scholar] [CrossRef] [Green Version]
  12. Buxton, D.; Marten, G.; Hornstein, J. Genetic variation for forage quality of alfalfa stems. Can. J. Plant Sci. 1987, 67, 1057–1067. [Google Scholar] [CrossRef]
  13. Buxton, D.R. Quality-related characteristics of forages as influenced by plant environment and agronomic factors. Anim. Feed Sci. Technol. 1996, 59, 37–49. [Google Scholar] [CrossRef]
  14. Albrecht, K.A.; Wedin, W.F.; Buxton, D.R. Cell-wall composition and digestibility of alfalfa stems and leaves. Crop Sci. 1987, 27, 735–741. [Google Scholar] [CrossRef]
  15. Wang, Z.; Qiang, H.; Zhao, H.; Xu, R.; Zhang, Z.; Gao, H.; Wang, X.; Liu, G.; Zhang, Y. Association Mapping for Fiber-Related Traits and Digestibility in Alfalfa (Medicago sativa). Front. Plant Sci. 2016, 7, 331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Jung, H.; Sheaffer, C.; Barnes, D.; Halgerson, J. Forage quality variation in the US alfalfa core collection. Crop Sci. 1997, 37, 1361–1366. [Google Scholar] [CrossRef]
  17. Bertrand, A.; Claessens, A.; Thivierge, M.-N.; Rocher, S.; Lajeunesse, J.; Castonguay, Y.; Seguin, P. Field Assessment of Alfalfa Populations Recurrently Selected for Stem Cell Wall Digestibility. Crop Sci. 2018, 58, 1632–1643. [Google Scholar] [CrossRef] [Green Version]
  18. Jung, H.G.; Lamb, J.F.S. Stem Morphological and Cell Wall Traits Associated with Divergent In Vitro Neutral Detergent Fiber Digestibility in Alfalfa Clones. Crop Sci. 2006, 46, 2054–2061. [Google Scholar] [CrossRef] [Green Version]
  19. Laboski, C.; Lamb, J.; Dowdy, R.; Baker, J.; Wright, J. Irrigation scheduling for a sandy soil using mobile frequency domain reflectometry with a checkbook method. J. Soil Water Conserv. 2001, 56, 97–100. [Google Scholar]
  20. Kaiser, D.E.; Lamb, J.A.; Eliason, R. Fertilizer Guidelines for Agronomic Crops in Minnesota. Available online: https://conservancy.umn.edu/bitstream/handle/11299/198924/Fertilizer%20Guidelines%20for%20Agronomic%20Crops%20in%20Minnesota.pdf?sequence=1 (accessed on 2 September 2022).
  21. Van Soest, P.v.; Robertson, J.B.; Lewis, B.A. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 1991, 74, 3583–3597. [Google Scholar] [CrossRef]
  22. Jung, H.-J.G.; Lamb, J.F. Identification of lucerne stem cell wall traits related to in vitro neutral detergent fibre digestibility. Anim. Feed Sci. Technol. 2003, 110, 17–29. [Google Scholar] [CrossRef]
  23. McDougall, E. Studies on ruminant saliva. 1. The composition and output of sheep’s saliva. Biochem. J. 1948, 43, 99. [Google Scholar] [CrossRef] [Green Version]
  24. Vogel, K.P.; Pedersen, J.F.; Masterson, S.D.; Toy, J.J. Evaluation of a filter bag system for NDF, ADF, and IVDMD forage analysis. Crop Sci. 1999, 39, 276–279. [Google Scholar] [CrossRef]
  25. Team, R.C. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  26. Mendiburu, F.D.; Yaseen, M. Agricolae: Statistical procedures for agricultural research. R package version 1.4. 0. 2020. Available online: https://cran.r-project.org/web/packages/agricolae/index.html (accessed on 2 March 2023).
  27. Wei, T.; Simko, V.; Levy, M. Package “Corrplot”: Visualization of a Correlation Matrix. 2017. Version 0.84 2021., R package, online resources. Available online: https://cran.microsoft.com/snapshot/2016-08-01/web/packages/corrplot/vignettes/corrplot-intro.html#:~:text=The%20corrplot%20package%20is%20a%20graphical%20display%20of,choosing%20color%2C%20text%20labels%2C%20color%20labels%2C%20layout%2C%20etc. (accessed on 2 March 2023).
  28. Warton, D.I.; Duursma, R.A.; Falster, D.S.; Taskinen, S. smatr 3– an R package for estimation and inference about allometric lines. Methods Ecol. Evol. 2012, 3, 257–259. [Google Scholar] [CrossRef]
  29. Falconer, D.S. Introduction to Quantitative Genetics, 4th ed.; Longman: Harlow, UK, 1996. [Google Scholar]
  30. Villanueva, R.A.M.; Chen, Z.J. ggplot2: Elegant Graphics for Data Analysis (2nd ed.). Meas. Interdiscip. Res. Perspect. 2019, 17, 160–167. [Google Scholar] [CrossRef]
  31. Kassambara, A.; Kassambara, M.A. Package ‘Ggpubr’. Available online: https://CRAN.R-project.org/package=ggpubr (accessed on 31 January 2023).
  32. Olivoto, T.; Lúcio, A.D.C. metan: An R package for multi-environment trial analysis. Methods Ecol. Evol. 2020, 11, 783–789. [Google Scholar] [CrossRef] [Green Version]
  33. Kephart, K.D.; Buxton, D.R.; Hill, R.R., Jr. Digestibility and cell-wall components of alfalfa following selection for divergent herbage lignin concentration. Crop Sci. 1990, 30, 207–212. [Google Scholar] [CrossRef]
  34. Kephart, K.D.; Buxton, D.; Hill, R., Jr. Morphology of alfalfa divergently selected for herbage lignin concentration. Crop Sci. 1989, 29, 778–782. [Google Scholar] [CrossRef]
  35. Damiran, D.; Biligetu, B.; Pearce, L.; Lardner, H. PSXI-15 Evaluation of low-lignin alfalfa ‘Hi-Gest® 360′ on the Canadian prairies: Productivity, nutrient profile, and rumen degradation kinetics. J. Anim. Sci. 2021, 99, 348. [Google Scholar] [CrossRef]
  36. Jungers, J.; Cherney, J.; Martinson, K.; Jaqueth, A.; Sheaffer, C. Forage nutritive value of modern alfalfa cultivars. Crop Forage Turfgrass Manag. 2020, 6, e20076. [Google Scholar] [CrossRef]
  37. Grev, A.M.; Wells, M.S.; Samac, D.A.; Martinson, K.L.; Sheaffer, C.C. Forage Accumulation and Nutritive Value of Reduced Lignin and Reference Alfalfa Cultivars. Agron. J. 2017, 109, 2749–2761. [Google Scholar] [CrossRef] [Green Version]
  38. Grev, A.M.; Wells, M.S.; Catalano, D.N.; Martinson, K.L.; Jungers, J.M.; Sheaffer, C.C. Stem and leaf forage nutritive value and morphology of reduced lignin alfalfa. Agron. J. 2020, 112, 406–417. [Google Scholar] [CrossRef]
  39. Arnold, A.M.; Cassida, K.A.; Albrecht, K.A.; Hall, M.H.; Min, D.; Xu, X.; Orloff, S.; Undersander, D.J.; Santen, E.; Sulc, R.M. Multistate Evaluation of Reduced-Lignin Alfalfa Harvested at Different Intervals. Crop Sci. 2019, 59, 1799–1807. [Google Scholar] [CrossRef] [Green Version]
  40. Jung, H.; Engels, F. Alfalfa stem tissues: Cell wall deposition, composition, and degradability. Crop Sci. 2002, 42, 524–534. [Google Scholar] [CrossRef]
  41. Lin, S.; Medina, C.A.; Boge, B.; Hu, J.; Fransen, S.; Norberg, S.; Yu, L.-X. Identification of genetic loci associated with forage quality in response to water deficit in autotetraploid alfalfa (Medicago sativa L.). BMC Plant Biol. 2020, 20, 1–18. [Google Scholar] [CrossRef]
  42. Annicchiarico, P.; Nazzicari, N.; Li, X.; Wei, Y.; Pecetti, L.; Brummer, E.C. Accuracy of genomic selection for alfalfa biomass yield in different reference populations. BMC Genom. 2015, 16, 1–13. [Google Scholar] [CrossRef] [Green Version]
  43. Marten, G.; Buxton, D.; Barnes, R. Feeding value (forage quality). Alfalfa Alfalfa Improv. 1988, 29, 463–491. [Google Scholar]
  44. Biazzi, E.; Nazzicari, N.; Pecetti, L.; Brummer, E.C.; Palmonari, A.; Tava, A.; Annicchiarico, P. Genome-wide association mapping and genomic selection for alfalfa (Medicago sativa) forage quality traits. PLoS ONE 2017, 12, e0169234. [Google Scholar] [CrossRef] [Green Version]
  45. Jia, C.; Zhao, F.; Wang, X.; Han, J.; Zhao, H.; Liu, G.; Wang, Z. Genomic prediction for 25 agronomic and quality traits in alfalfa (Medicago sativa). Front. Plant Sci. 2018, 9, 1220. [Google Scholar] [CrossRef] [Green Version]
  46. Moore, J.H.; Williams, S.M. Epistasis and its implications for personal genetics. Am. J. Hum. Genet. 2009, 85, 309–320. [Google Scholar] [CrossRef] [Green Version]
  47. Howard, R.; Carriquiry, A.L.; Beavis, W.D. Parametric and nonparametric statistical methods for genomic selection of traits with additive and epistatic genetic architectures. G3 Genes Genomes Genet. 2014, 4, 1027–1046. [Google Scholar] [CrossRef] [Green Version]
  48. Gianola, D.; de Los Campos, G. Inferring genetic values for quantitative traits non-parametrically. Genet. Res. 2008, 90, 525–540. [Google Scholar] [CrossRef]
  49. De los Campos, G.; Gianola, D.; Rosa, G.J.; Weigel, K.A.; Crossa, J. Semi-parametric genomic-enabled prediction of genetic values using reproducing kernel Hilbert spaces methods. Genet. Res. 2010, 92, 295–308. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Pearson correlation among the digestibility and CW traits. Numeric values in the bottom lower left section are the correlation coefficients. Names along the diagonal are the variables. Size and color of the circles on the top right section show how strongly the variables are correlated.
Figure 1. Pearson correlation among the digestibility and CW traits. Numeric values in the bottom lower left section are the correlation coefficients. Names along the diagonal are the variables. Size and color of the circles on the top right section show how strongly the variables are correlated.
Agronomy 13 00770 g001
Figure 2. Responses in IVNDFD and CW components over two cycles of divergent selection. Populations selected for high 16 h and high 96 h IVNDFD are in red, and populations selected for low 16 h IVNDFD and low 96 h IVNDFD are in blue. Gray shading around the regression line is the 95% confidence interval around the mean across years, locations, harvests, and maturities. Regression equations and R2 over the two cycles are shown above or below each line.
Figure 2. Responses in IVNDFD and CW components over two cycles of divergent selection. Populations selected for high 16 h and high 96 h IVNDFD are in red, and populations selected for low 16 h IVNDFD and low 96 h IVNDFD are in blue. Gray shading around the regression line is the 95% confidence interval around the mean across years, locations, harvests, and maturities. Regression equations and R2 over the two cycles are shown above or below each line.
Agronomy 13 00770 g002
Figure 3. Effect of maturity and harvest on 96 h IVNDFD in selected populations. (A) Ninety-six h IVNDFD at early bud (EB), late flower (LF), and green pod (GP) maturity stages across two years, two locations, and three harvests. (B) Ninety-six h IVNDFD for combined maturity stages at Harvest 1, 2, and 3 across two years, two locations, and three maturities. The small dots are the individual values, and the large black dots are outliers. The red diamond inside each box is the mean value. C0, unselected population; C1, cycle 1 population; C2, cycle 2 population; H × H, intermating of plants with high 16 h IVDFD and high 96 h IVNDFD; L × L, intermating of plants with low 16 h IVDFD and low 96 h IVNDFD.
Figure 3. Effect of maturity and harvest on 96 h IVNDFD in selected populations. (A) Ninety-six h IVNDFD at early bud (EB), late flower (LF), and green pod (GP) maturity stages across two years, two locations, and three harvests. (B) Ninety-six h IVNDFD for combined maturity stages at Harvest 1, 2, and 3 across two years, two locations, and three maturities. The small dots are the individual values, and the large black dots are outliers. The red diamond inside each box is the mean value. C0, unselected population; C1, cycle 1 population; C2, cycle 2 population; H × H, intermating of plants with high 16 h IVDFD and high 96 h IVNDFD; L × L, intermating of plants with low 16 h IVDFD and low 96 h IVNDFD.
Agronomy 13 00770 g003
Figure 4. Environmental stability of 96 h IVNDFD with selection cycles across the 12 harvest environments. (Left) H × H, intermating of plants with high 16 h IVDFD and high 96 h IVNDFD. (Right) L × L, intermating of plants with low 16 h IVDFD and low 96 h IVNDFD.
Figure 4. Environmental stability of 96 h IVNDFD with selection cycles across the 12 harvest environments. (Left) H × H, intermating of plants with high 16 h IVDFD and high 96 h IVNDFD. (Right) L × L, intermating of plants with low 16 h IVDFD and low 96 h IVNDFD.
Agronomy 13 00770 g004
Table 1. Summary of populations used in the study. H16, high 16 h IVNDFD; H96, high 96 h IVNDFD; L16, low 16 h IVNDFD; L96, low 96 h IVNDFD. C0, unselected population; C1, cycle 1 population; C2, cycle 2 population.
Table 1. Summary of populations used in the study. H16, high 16 h IVNDFD; H96, high 96 h IVNDFD; L16, low 16 h IVNDFD; L96, low 96 h IVNDFD. C0, unselected population; C1, cycle 1 population; C2, cycle 2 population.
Population
Names
Population
Combinations
Cycle
Number
Number of Plants
Intermated
UMN3097ParentalC0
UMN3355H16 × H96C1117
UMN3356H16 × L96C128
UMN3357L16 × H96C126
UMN3358L16 × L96C133
UMN4016H16 × H96C260
UMN4017H16 × L96C2~30
UMN4018L16 × H96C2~30
UMN4019L16 × L96C2~30
Table 2. Mean squares summary from the analysis of variance (ANOVA) for 16 h IVNDFD (NDFD16), 96 h IVNDFD (NDFD96), acid detergent lignin (ADL) as a proportion of NDF, neutral detergent fiber (NDF), stem dry weight (SDW), and percent stems (Per_stem) among the following five variables: year (Y), location (L), germplasm (G), maturity (M), and harvest (H).
Table 2. Mean squares summary from the analysis of variance (ANOVA) for 16 h IVNDFD (NDFD16), 96 h IVNDFD (NDFD96), acid detergent lignin (ADL) as a proportion of NDF, neutral detergent fiber (NDF), stem dry weight (SDW), and percent stems (Per_stem) among the following five variables: year (Y), location (L), germplasm (G), maturity (M), and harvest (H).
TermNDFD16
g kg−1
NDFD96
g kg−1
ADL g kg−1 NDFNDF
g kg−1 DM
SDW
g plant−1
Per_stem
%
Y 61,131 ***37,685 ***3868 ***79,412 ***851 *1471 ***
M 136,873 ***292,141 ***20,875 ***789,443 ***71,707 ***21,300 ***
L 11,393 ***191,234 ***49,474 ***12,798 ***3510 ***996 ***
H 20,023 ***39,845 ***10,591 ***62,982 ***106,104 ***33,054 ***
G 5023 ***51,724 ***5257 ***27,206 ***2587 ***201 ***
H:G 173 ***846 ***31919 ***317 *24
L:G 89768 *100 ***689 *464 *15
M:G 136 **1609 ***75 **823 ***714 ***33
Y:G 99893**6159911818
L:H:G 663122432010925
L:M:G 137 **756 **50478914 ***38 *
M:H:G 702561937212825
Y:H:G 89274323835517
Y:L:G 142 *1258 ***128 ***50516415
Y:M:G 383375327023618
L:M:H:G 472782426827322
Y:L:H:G 50184202707118
Y:L:M:G 664303738619334
Y:M:H:G 602821721911424
Y:L:M:H:G 401201520111014
* Significant at the 0.05 level of probability. ** Significant at the 0.01 level of probability. *** Significant at the 0.001 level of probability.
Table 3. Summary of heritability, genetic gain (GG), genetic regression slope, and p-values of the regression.
Table 3. Summary of heritability, genetic gain (GG), genetic regression slope, and p-values of the regression.
TraitHeritabilityGGSlopep-Value
ADL (g kg−1 NDF)0.75−3.98−3.421.17 × 10−9
NDFD96H (g kg−1)0.704.5410.132.50 × 10−9
NDF (g kg−1 DM)0.47−3.15−9.794.60 × 10−6
NDFD16H (g kg−1)0.464.274.181.84 × 10−5
SDW_gm (g plant−1)0.15−16.1−1.318.34 × 10−2
Per_Stem (%)0.11−2.02−0.012.62 × 10−1
Table 4. Responses in 16 h and 96 h NDFD in the H16 × H96 populations with cycles of selection at the early bud (EB), late flowering (LF), and green pod (GP) maturity stages over three harvests. LSD is Fisher’s least significant difference test with different letters within a maturity stage indicating a significant difference; %/cycle is the percent change from C0 to C1 and C1 to C2.
Table 4. Responses in 16 h and 96 h NDFD in the H16 × H96 populations with cycles of selection at the early bud (EB), late flowering (LF), and green pod (GP) maturity stages over three harvests. LSD is Fisher’s least significant difference test with different letters within a maturity stage indicating a significant difference; %/cycle is the percent change from C0 to C1 and C1 to C2.
MaturityHarvestCycleNDFD
16H (g kg−1)
LSD
16H
%/cycle
16H
NDFD
96H (g kg−1)
LSD
96H
%/Cycle
96H
EB1C2218a2.25485a2.03
EB1C1213a3.02476a0.25
EB1C0207a 475a
EB2C2210a1.39468a1.27
EB2C1207ab0.89462a1.55
EB2C0205b 455a
EB3C2192a3.44440a3.84
EB3C1186a−0.41424a−1.28
EB3C0187a 430a
LF1C2182a0.54442a2.47
LF1C1181a3.81431a4.11
LF1C0174b 414b
LF2C2188a1.50456a3.86
LF2C1186ab3.05439b2.86
LF2C0180b 427c
LF3C2197a4.27445a5.33
LF3C1189ab1.43422b1.21
LF3C0186b 417b
GP1C2178a1.66417a3.22
GP1C1175ab1.29404b1.42
GP1C0173b 398b
GP2C2191a4.75450a4.08
GP2C1182b2.81433b3.65
GP2C0177b 417c
GP3C2199a4.11452a5.09
GP3C1191b1.91430b2.41
GP3C0188b 420b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, Z.; Heuschele, D.J.; Lamb, J.F.S.; Jung, H.-J.G.; Samac, D.A. Improved Forage Quality in Alfalfa (Medicago sativa L.) via Selection for Increased Stem Fiber Digestibility. Agronomy 2023, 13, 770. https://doi.org/10.3390/agronomy13030770

AMA Style

Xu Z, Heuschele DJ, Lamb JFS, Jung H-JG, Samac DA. Improved Forage Quality in Alfalfa (Medicago sativa L.) via Selection for Increased Stem Fiber Digestibility. Agronomy. 2023; 13(3):770. https://doi.org/10.3390/agronomy13030770

Chicago/Turabian Style

Xu, Zhanyou, Deborah J. Heuschele, JoAnn F. S. Lamb, Hans-Joachim G. Jung, and Deborah A. Samac. 2023. "Improved Forage Quality in Alfalfa (Medicago sativa L.) via Selection for Increased Stem Fiber Digestibility" Agronomy 13, no. 3: 770. https://doi.org/10.3390/agronomy13030770

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop