Introduction
Maize production and the production of dairy milk represent two of the main economic activities in México (Reta et al., 2015; Espinoza et al., 2007). According to SAGARPA (2016), there is a population of 2.3 million dairy cows, of which 85 % are located on small-scale farms, contributing to approximately 70 % of the national milk supply peryear, with a reported annual production per cow of 5190 L (Posadas et al., 2016). In addition, the national production of dairy milk for the second quarter of 2017 reached 5670 million liters (SIAP, 2017). In 2016, Mexico imported 209,803 t of milk powder to cover national supply needs, a number which is expected to increase (Brio Agropecuario, 2016).
In 2013, Mexico was identified as one of the countries most affected by climate change (SEMARNAT, 2014). Furthermore, maize production is the main farming activity in Mexico. Nearly 2 million producers participate in this activity, and 85 % have less than 5 ha of land. In Mexico, maize represents the main use crop to produce dairy milk and human consumption (Jiménez-Leyva et al., 2016), and an undetermined amount is allocated as straw, green fodder, and to a lesser extent for the preparation of silages for cattle feed (Celis Álvarez et al., 2016; Jiménez Leyva et al., 2016). According to the SIAP (2016) reports, in Mexico in 2015, an area of 445,775 ha was planted in the rainy season and 161,623 ha in irrigation for fodder maize, with yields of 19.29 and 47.55 t/ha of dry matter (DM) and green matter (GM), respectively. Given the great heterogeneity of agroclimatic conditions that present negative impacts on agricultural-animal production, this results in a disparity of yields per hectare and per animal, which is why there is a need to optimize the use of forage.
The StAnD (sustainable animal diets) method (FAO, 2014) is a tool that integrates several dimensions of sustainability, including the three P (people, planet, and profitability) dimensions, and gives an overall picture of the current state of a production system. The indicators corresponding to each dimension allow for the detection of specific problems or limitations that may be addressed to improve the sustainability of the system (FAO, 2014). One indicator of the StAnD method is “do not use cereals in animal diets and improve the use of native resources” (Planet dimension). This study used the StAnD method to evaluate the sustainability of native and hybrid silages in Mexico and can help to guide agricultural practices and policies in accordance with the economic and environmental performance of different maize production systems. The objective of the present study was to perform a literature search on the yield and chemical composition of maize silage samples produced in tropical and temperate regions of Mexico and to characterize maize silage according to chemical composition, forage yield (t/ha), as well as predicted milk production.
Materials and methods
Data collection.
The selection process limited the results to studies published from 2001 to 2016. For inclusion in the final database, the studies had to include agronomic and chemical variables such as: yield of dry matter (t/ha), density of plants (number of plants/ha), dry matter content (DM), crude protein (CP), neutral detergent fiber (NDF), non-fibrous carbohydrates soluble (NFCS), ash or organic matter (OM), dry matter digestibility (DMD), and neutral detergent fiber digestibility (NDFD). Thus, in this study, 142 maize silage samples were included (135 hybrids and seven native varieties); 41 maize varieties from temperate regions and 101 maize varieties from tropical origin.
Calculations
Studies with missing values for starch and fat were calculated according to the National Research Council (2001). The net energy for lactation (NEL, Mcal/kg DM), total digestible nutrients (TDN), kilograms of milk per t of dry matter (kg of milk/t DM) and kilograms of milk per hectare (kg milk/ha) were determined using the MILK2006®spreadsheet. The model utilizes the concentration of NEL, estimated from adaptation of the equations for NEL provided by the National Research Council (NRC, 2001), in combination with the NDF in vitro digestibility (NDFD), to predict the milk yield with a fat content of 35 g/kg. The estimated fat-corrected milk yield (FCM) is a means of adjusting the milk yield for the amount of fat in the milk to reflect the relative energy concentration in the milk, thus it reflects the energy required to produce the given amount of milk. The feed value of the fiber affects the supply of NEL, both through altering the concentration of NEL and through its effect on forage intake. In the model, milk production is estimated from a cow with a live weight of 612 kg, fed a diet with 300 g/kg dietary NDF concentration, with maize silage as the only forage, and a production of 35 kg of milk (kg per day) (Shaver, 2006). The missing values for NDFD were calculated using the regression equation derived from both NDF and NDFD data available in the papers used in the present study as:
A descriptive analysis of the chemical composition, forage yield (t/ha), as well as the potential milk production of the maize silage was carried out using the SAS statistical software (Statistical Analysis System [SAS], 2004). To identify differences between varieties of maize silage, Kolmogorov-Smirnov tests were applied to determine if the resulting scores varied significantly with respect to a normal distribution. The data were normally distributed and analyzed with a completely randomized design model, using the Tukey’s test (P < 0.05).
Statistical analysis
The analysis of the variables described above was carried out in two independent stages: 1) With the information of the sixteen variables, an analysis was made of the interrelationships between the variables (multivariate analysis) and their estimated contribution to the total variance, applying a Principal Components Analysis (PCA) (PROC PRINCOMP, SAS Institute INC., 2007). 2) The grouping of information sources into homogeneous levels was done through the hierarchical clustering analysis (PROC CLUSTER, WARD method) from SAS Institute INC., (2007) using the selected variables. One hundred and forty-two maize silages were evaluated by taking the mean of treatments of each, to determine which groups of variables were correlated with each other. Six clusters were obtained. According to the six integrated groups, an analysis of variance was carried out to determine if there were significant statistical differences between the six levels (P < 0.05). From the maize silage samples included in this study, 101 varieties were from tropical zones and 41 varieties from temperate zones.
Finally, maize groups that had similar productive conditions were determined from the variables considered for this study. The effects were considered significant if they were lower than P < 0.05. To determine if there were differences between the effects, the Tukey test was used for multiple comparisons of means (Steel & Torrie, 1997, pp. 179-180).
Results
Descriptive statistics of chemical composition, forage yield, and milk production potential are shown in Table 1. The PCA results are presented in Table 2. Five factors (F) were obtained and named according to the variables included in each of them. Factor 1 (F1) refers to kg milk/t DM integrated by the variables NDF, NDFD, NFCS, starch, TDN, and kg milk/t DM, and F1 had the highest contribution (38.48 %) to the total variance. Factor 2 (F2) refers to DM yield/t and kg milk/ha, factor 3 (F3) to fat content, factor 4 (F4) to CP and OM, and factor 5 (F5) to DM and DMD. F2- F5 contributed 15.24 %, 11.49 %, 8.09 %, and 7.57 %, respectively, of the total variance.
Chemical composition is expressed in g/100g Dry Matter. t/ha = yield in tons per hectare of dry matter; Plant density/ha = highest plant density; DM = dry matter content; DMD = dry matter digestibility; CP = crude protein; NDF = neutral detergent fiber; NDFD = neutral detergent fiber digestibility; NFCS = non fibrous carbohydrates soluble; OM = organic matter content; TDN = total digestible nutrients; NEL = net energy for lactation (Mcal/kg DM).
CL1_DMD, CP, NDFD, NFCS,TDN, NEL, and kg milk/t DM, C2_highest CP content, C3_highest plant density (number plats/ha), NDF and the lowest NFCS, TDN, NEL, kg milk/t DM and kg milk/ha, C4_highest NDF and OM content, C5_higher forage yield (t DM/ha), DM % content and kg milk/ha C6_higher DM digestibility and the lowest OM content.
In the PCA, 16 groups were obtained, however only the first five have an eigenvalue greater than one (Table 2), and they explained 81 % of the overall variance. Subsequently, the variables were used to perform a cluster analysis (cluster, CL) (Table 3) and resulted in a hierarchical dendrogram program (Figure 1), forming six cluster groups (CL) of maize silage. The first group (CL1) (22 varieties) was called “Digestibility”, the variables associated to it were: DMD, CP, NDFD, NEL, and kg of milk/t DM. The second cluster (CL2) (40 varieties) was characterized by maize silage varieties with the highest number of plants per hectare and NDF content. The third (CL3) (21 varieties) called “Ash” had varieties with the highest ash content; the fourth (CL4) (24 varieties) was called “Average” because it was made up of the average of all variables. The fifth (CL5) (7 varieties) was called “milk yield per ha”, which had the highest forage yield (t DM/ha), and kg of milk/ha, and the sixth cluster (CL6) (7 varieties) was called “kg Milk/t” which had the greater DM, NDF, and NEL digestibility (Mcal/kg DM), and kg milk/DM t.
DM = dry matter content; DMD = dry matter digestibility-, CP = crude protein; NDF = neutral detergent fiber; NDFD = neutral detergent fiber digestibility; NFC = nonfibrous carbohydrates- OM = organic matter- tDn = total digestible nutrients- NEL = net energy for lactation (MJ/kg DM); kg milk/t DM = kilograms of milk per t of dry matter; Kg milk/ha = kilograms of milk per hectare; CL1_DMD, CP, NDFD, NFCS, TDN, NEL, and kg milk/t DM; CL2_highest CP content; CL3_ highest plant density (number plats/ha), NDF, and the lowest NFCS, TDN, NeL, kg milk/t DM and kg milk/ha; CL4_highest NdF % and OM % content; CL5_higher forage yield (t DM/ha), DM % content, and kg milk/ha; CL6_higher DM digestibility and the lowest OM content.
Discussion
Maize silages were characterized by the variability in their chemical composition due to different factors such as the use of local varieties (natives), hybrids, the agronomic management carried out with the crop, and the content of dry matter (Darby & Lauer, 2002; Komainda et al., 2016). In Mexico, different maize genotypes are cultivated, however, most of them are planted for their yield (t DM/ha) rather than for their nutritional quality, which is very common in most countries producing maize (Núñez et al., 2003). Similarly, for the different types of silage, development and characterization should be considered as a function of their chemical composition, maize varieties, and harvesting management (fertilization rate, irrigation, crop labors) (Darby & Lauer, 2002; Lynch et al., 2012; Marchesini et al., 2019). These aspects play a key role in the development of a good silage (silage processing after cutting, i.e. particle size, pressing, bailing, fermentation rate) with a potential for improving milk production.
Productivity of maize varieties in Mexico has been reported to range from 10 to 22 t/ha, whereas plant density varies between 58,000 to 96,000 plants/ ha. In addition, DM digestibility of those cultivars ranged from 63 to 79 % (Núñez et al., 2001; Elizondo & Boschini, 2002; Núñez et al., 2003; SAGARPA, 2016). In Mexico, maize silage is considered a secondary product since the main use is the production of grain for human consumption.
In this study, the PCA help to characterize maize silage according to its chemical composition, forage yield and milk production potential, resulting in five factors, which could be used for the selection of high and low-quality maize silages. F1 grouped most of the variables such as kg milk/t DM as well as chemical composition variables, for example the NDFD, which is useful for the characterization of forages quality (Gallo et al., 2013).
Six cluster groups were obtained from the chemical and productive characteristics of the different maize silages. One of the main factors determining the nutritional quality of the silage is the NDF digestibility, which has variable degradability and influences the energy value and animals’ intake (Spanghero et al., 2008; Khan et al., 2015). In a study on forage yield and quality of maize silage, Lynch et al. (2012), reported an increase in NDF in maize harvested at 22 days of maturity, which indicates that the determination of maize harvesting time is a key factor for silage quality according to its chemical characteristics and its potential for milk production.
In CL2, maize silage varieties were grouped according to characteristics for greater density, which is consistent with the tropical maize varieties from this study. The forage maize of tropical origin, known for having smaller stems and fewer leaves, decreased DM production, in addition to requiring more plants to reach an optimal yield. Some authors reported that this increase is due to the lower number of leaves and the fibrous and highly lignified stems (Núñez et al., 2001; Elizondo & Boschini, 2002). A disadvantage of this group is the high NDF concentration in the forage, which limits total DM intake, due to a filling effect at the ruminal level, as well as a decrease in the degradation of the forage (Spanghero et al., 2008).
The amount of ash is the inorganic content of the forage, which can be affected by maturity stage. It has been reported that maize harvested in the last phenological stages has higher contents of ash than in early growth stages (Darby & Lauer, 2002; Komainda et al., 2016). In this sense, it has been indicated that maize silages can be harvested around 115 days post sowing when 95 % of OM is present (Filya, 2004). Marchesini et al. (2019) reported how the effects of yield potential and maturity at harvest vary depending on whether early or late hybrids are considered resulting in a progressive reduction in the fermentation quality at higher maturity stages by increasing their ash and NDF contents. However, it has been reported that neither the density of plants by area, nor the use of higher densities of maize plants allow an increase in crop yield per unit area (Khan et al., 2015).
In Mexico, the selection of maize silage for livestock feed does not consider varieties based on nutritional quality. Instead of selecting maize varieties for grain production, determining the fermentation quality index (QFI) could be a better option, which is an index with values ranging from 1 to 100 and is calculated based on the concentrations of lactic acid, ammonia, ethanol, butyric, acetic acids (Kung et al., 2018), and pH. However, this information is not always available for comparison or to determine with greater certainty the yield, quality, and production potential of silage maize varieties.
In this sense, it has been shown that when the percentage of dry matter of maize silage increases from 25 % to 30 %, the intake and production of milk are higher (Khan et al., 2015). This increase is related to the percentage of DM, the intake and nutritional value of the maize, due to a higher content of grains, water structural carbohydrates, and energy. Similarly, higher concentrations of NEL and CP stimulate the production of microbial crude protein in the rumen, increasing the concentration of protein in the milk (Ferraretto & Shaver, 2015; Lascano et al., 2016), positively influencing milk yield production (kg of milk/t DM).
The results of the CL6 may be due to DMD and NDFD contents resulting from the degree of maturity obtained at the time of cutting and by the variety of maize, as well as the factors affecting forage yield and milk production per t of dry matter (kg of milk/t of DM) (Lynch et al., 2012; Khan et al., 2015).
Conclusion
Cluster five (CL5) was characterized for having a higher forage yield (t DM/ha), DM content, and kg milk/ha. CL5 also showed the highest DM and energy to produce milk (kg of milk per t of DM). Our results suggest that the most suitable option was the maize silage that had a higher forage yield (t DM per ha) and more than 35 % DM, since it produced more kg of milk/ha. In the same way, the characteristics of CL6 were of interest, considering its DM and NDF digestibility. Results from this review could help farmers when decisions are needed for silage conservation.