Skip to main content

Towards an ideotype for food-fuel dual-purpose wheat in Argentina with focus on biogas production



Wheat straw, one of the most abundant agricultural residues worldwide, can be used to produce biogas, which is considered one of the most efficiently produced renewable energies. Wheat grown with the dual-purpose of producing food and biogas should display simultaneously high grain and straw yields, low lodging susceptibility and high conversion efficiency of straw into biogas. The aims of this study were to determine the best food-fuel dual-purpose wheat candidates among 36 wheat genotypes—including French, CIMMYT and local (Criollo) germplasm—used in breeding programs in Argentina and to gain some insights into the relationships between key traits relevant for dual-purpose wheat genotypes.


High variability in individual key traits for dual-purpose wheat ideotype has been found. Genotypes of French origin displayed the highest grain yield, and those of CIMMYT origin, the lowest straw yield. Genotypes of Criollo origin showed the highest lodging susceptibility, and French ones, the lowest. Straw yield was positively correlated with grain yield in all genotypes, and negatively correlated with plant height in genotypes of Criollo origin. Straw conversion into biogas was measured in terms of the biogas potential production and kinetic parameters Bmax (maximum specific biogas production) and k (first-order kinetic constant) were analyzed. All key traits were analyzed together by a principal component analysis. Baguette 31 and SNR Nogal, two genotypes of French origin, showed high grain yield, high-to-very high straw yield, low lodging susceptibility, and moderate-to-high Bmax and k. Buck Guapo and Buck Baqueano, two genotypes of Criollo origin, displayed good values for grain yield, straw yield and Bmax and k. However, their high lodging susceptibility precludes their production in shallow soils or high-input systems. Lastly, some old genotypes (e.g. Klein Atlas) harbored a good combination of all key traits and could prove valuable to be included in future breeding programs for dual-purpose wheat.


While none of the genotypes excelled in every key trait, a few candidates showed potential for dual-purpose ideotype, particularly Baguette 31 and SNR Nogal. The challenge lays in gathering all attributes for food and 2G fuel in the same genotype.


The contribution of fossil fuels to global climate change and the growing demand for energy are the main drivers for worldwide biofuels development [1]. In this context, biomass emerges as one of the most important resources of energy for the sustainable development of societies, in a context of circular economy and environmental protection [2]. Lignocellulosic biomass, the most abundant organic carbon source in the world, is one of the most attractive alternatives for second-generation (2G) biofuel production due to its high and ubiquitous availability and low cost [3, 4]. Straw is a lignocellulosic material largely produced by agriculture, being wheat (Triticum aestivum L.) straw one of the most abundant agricultural residues on earth [5]. In Argentina, wheat production is projected to reach a record of 20.2 million tons by 2021 [6], anticipating an approximate amount of available straw of 26.3 million tons [7].

Wheat straw has multiple applications, including use as animal bedding, mushroom-production substrate, and feedstock for biomass-burning power stations [8]. Its supply greatly exceeds its demand, and a large amount of straw is incorporated into the soil after grain harvest [9]. Hence, a significant amount of wheat straw could be used as feedstock for bioenergy production without threatening other straw uses or soil fertility [10]. Different types of energy can be derived from wheat straw through thermochemical or biochemical processes: liquid fuels such as bioethanol, gaseous fuels such as biogas (composed mainly of methane and carbon dioxide), and electricity from direct combustion [11]. Among them, biogas production technology is considered one of the most efficient for producing renewable energy, since it can achieve high energy recovery, and economic and environmental benefits [2, 12]. Recently, an innovative concept of cereals as dual-purpose (DP) for food (wheat is the staple for 35% of the world’s population) and 2G-biofuel was introduced [13]. As mentioned by Townsend et al. [14], an ideal genotype (ideotype) for these purposes would be characterized by high grain and straw yields, good lodging resistance and high straw digestibility.

Effects of the environment, the genotype, and their interaction on wheat grain yield (GY) and its components have been largely studied [15,16,17,18]. Historically, wheat breeding has been focused on increasing GY per unit area and progress in breeding for GY potential is well documented in the literature [19,20,21,22]. Breeding for novel wheat genotypes with reduced plant height has allowed to increase genetic gains and significantly contributed to increase wheat productivity globally [23,24,25,26]. Some authors observed that this increase was achieved without a significant change in total aboveground biomass [27,28,29], while others claimed that GY increases came at the expense of straw yield (SY) [19]. Moreover, some agronomic practices designed to increase GY have resulted in weaker stems and surface roots, and, hence, greater lodging risk [30,31,32], which in turn impair further GY increase. However, it has been suggested that additional improvements in GY could be achieved by increasing photosynthetic capacity through the optimization of biomass production while maintaining lodging resistance [23].

To successfully manage wheat straw as a commodity and to help position it in the biofuel industry, it is crucial to understand the correlations between biomass, GY, and SY [33], since selecting cultivars with higher SY and without negative effects on GY may constitute a mean to increase the overall straw resource. Even though SY is known to be influenced by several environmental and management factors, including cultivar selection, it is rarely reported in the literature [14, 33,34,35]. One of the most frequently used approaches to assess SY is based on empirical models that infer SY from GY statistics [36]. However, the uncertainties related to varietal differences in biomass partitioning and its complex response to different environmental and management factors [37] are difficult to address and, therefore, it becomes necessary to record SY in dedicated field trials [36, 38].

Straw digestibility is another relevant trait for DP wheat and it refers to the proportion of fermentable sugars released from lignocellulosic biomass cell wall carbohydrate polymers during hydrolysis. It has been shown to be the rate-limiting step in biomass transformation into bioenergy [39, 40]. Anaerobic digestibility, usually measured in terms of its biogas or biomethane potential defines the quality of the biomass as a feedstock for biogas production [41]. Biogas potential production of wheat straw has been modeled by first-order kinetics [42, 43], where the maximum specific biogas production (Bmax) and the first-order kinetic constant (k) are used to characterize the suitability of a feedstock for biogas production. According to this model, Bmax represents the maximal biogas that can be produced by volatile solids present in the biomass while k estimates the biogas production rate. Kinetic parameters of biogas potential not only define the value of the feedstock for biogas production, but they also give important information for the design and operation of anaerobic digesters [44]. Straw digestibility is often restricted by the complex composition of lignocellulosic material, consisting of three main types of polymers: cellulose, hemicellulose, and lignin. Recalcitrance is determined not only by the wall polymer features but also by the interactions among its components [45,46,47,48]. However, no correlation between straw digestibility and GY has been found so far and, hence, it might be possible to improve straw digestibility through breeding without any critical negative effects on GY [49].

The aims of this study were to determine the best food-fuel wheat ideotype candidates among 36 wheat genotypes currently used in breeding programs in Argentina and to give some insights into the relationships between key traits relevant for DP wheat genotypes. To this end, we quantified straw and grain yields and characterized biomass both compositionally and morphometrically. Kinetic parameters of biogas potential production and methane potential yield from straw were determined for all genotypes. Furthermore, correlations among traits were analyzed and genotypes candidates for DP use selected in an effort to assist future breeding programs.


Thirty-six wheat genotypes carrying germplasm from three different origins (i.e. CIMMYT, Criollo and French) were selected and cultivated during seasons 2014 and 2017 (Table 1). Although most genotypes were modern, old genotypes not carrying dwarfing genes (i.e. genotypes 19, 20, 22, 23, 24, 32 and 36) were also included. Under the framework proposed by Townsend et al. [14], key traits of the wheat DP ideotype were analyzed: (i) grain yield and straw yield (“Biomass production and allocation” section); (ii) morphology and lodging susceptibility (“Wheat straw morphology and lodging susceptibility” section), and (iii) biomass composition and biogas potential production (“Straw biomass composition and biogas potential production” section). Finally, key traits were integrated in a principal component analysis (PCA, “Target traits for ideotype design” section) and the best candidates for DP wheat among these genotypes were identified.

Table 1 Agronomic characteristics of the wheat genotypes

Biomass production and allocation

Grain yield (GY), straw yield (SY), and harvest index (HI, the ratio between GY and aboveground biomass production) were determined for all genotypes during the 2014 season. Genotypes were split into two groups according to their GY (high yield: over 4000 kg/ha; low yield: less than 4000 kg/ha), and from the combination of the two yield groups and the three origins, six genotype groups were defined: CIMMYT-high yield (CI-HY); CIMMYT-low yield (CI-LY); Criollo-high yield (CR-HY); Criollo-low yield (CR-LY); French-high yield (FR-HY); French-low yield (FR-LY).

Mean GY, SY and HI obtained across all genotypes during 2014 season were 4300 ± 1485 kg/ha, 8211 ± 1776 kg/ha, and 0.34 ± 0.06, respectively. Significant effects of genotype on all traits were detected (p < 0.01) and variability, by means of CV% were 35% for GY, 22% for SY, and 17% for HI. Germplasm origin significantly affected GY (p = 0.047), with GY values for genotypes of French origin (4591 ± 218 kg/ha) higher than those of Criollo and CIMMYT ones (4117 ± 370 and 4213 ± 333 kg/ha, respectively). Both GY and origin independently affected SY (p < 0.01). Mean SY across low-yield genotypes was lower than that of high-yield ones (7029 ± 292 and 9443 ± 274 kg/ha, respectively). Although a positive correlation (r = 0.74, p < 0.0001) between GY and SY was found, high variability of SY still exists in each genotype group. For instance, when only high-yield genotypes were considered (CV% of GY = 7.3), SY still displayed high variability (CV% = 14.0). Genotypes with CIMMYT origin showed lower SY (7511 ± 356 kg/ha) than Criollo and French ones (8894 ± 498 and 8534 ± 359 kg/ha, respectively). By defining the 75th percentile as a threshold value for SY (9317 kg/ha), genotypes 13, 14, 15, 16 and 17 (CR-HY group), and 29 and 30 (FR-HY group) were identified as genotypes with both high GY (high-yield group) and high SY (above the threshold). On a 90th percentile-based threshold (10,442 kg/ha of SY), only 15, 16 and 30 could be regarded as very high SY genotypes.

Origin and GY interacted in their effects on HI (p < 0.01): low-yield genotypes had lower HI than high-yield genotypes when considering CIMMYT and Criollo origin, while HI did not change with GY for French genotypes. Figure 1 shows biomass allometric analyses for all genotypes according to their origin. Allometric constants (i.e. slopes of relationships in Fig. 1) for CIMMYT and Criollo genotypes were significantly higher than 1 for GY and lower than 1 for SY, indicating that allocation to grains increased alongside plant size at the expense of biomass allocation to straw (Fig. 1a, b). For plants of the same size, wheat genotypes from CIMMYT origin displayed higher allocation to grains (and thus, lower to straw) than Criollo genotypes, in agreement with HI values, while French genotypes showed an intermediate behavior. Slopes close to 1 indicated that French genotypes conserved their biomass allocation pattern regardless of plant size (Fig. 1c).

Fig. 1

Allometric analyses of wheat biomass. Allometric relations for genotypes with germplasm of different origins: CIMMYT (a), Criollo (b), and French (c). Continuous lines represent linear regression between Ln of grain or straw yield and the Ln of shoot biomass. Inset to figures contain linear equations with best fitting parameters for each data set

Wheat straw morphology and lodging susceptibility

Plant morphology and lodging susceptibility index (∆LS, Eq. 1) were measured for all genotypes during the 2014 season (see Additional file 1). Genotype effects were significant on all traits measured (p < 0.005). Table 2 shows morphometric properties of straw and ∆LS of genotype groups. Interactive effects of GY and germplasm origin were found on stem diameter (p = 0.028): low-yield genotypes from CIMMYT origin showed larger stem diameter than high-yield ones, but no differences were observed between yield groups of French or Criollo origin. Stem wall thickness was significantly and independently affected by GY and origin (p < 0.05): low-yield genotypes showed a thicker stem wall than high-yield ones (0.55 ± 0.01 vs. 0.49 ± 0.01 mm, respectively); CIMMYT genotypes had thicker walls than French ones (0.55 ± 0.02 vs. 0.49 ± 0.01 mm, respectively), and Criollo genotypes showed intermediate values (0.53 ± 0.01 mm). Interactive effects of GY and germplasm origin were found on plant height (p = 0.009): CIMMYT and Criollolow-yield genotypes were taller than high-yield ones, but no differences were observed between yield groups of French origin. Plant height was negatively correlated with SY in genotypes of Criollo origin (r = − 0.59, p < 0.05), while no correlation between plant height and SY was found in any other genotype. No effects of either origin or GY were observed on spike dry weight (p > 0.005).

Table 2 Morphometric characteristics and lodging susceptibility of wheat genotype groups

Germplasm origin was the only factor affecting ∆LS (p < 0.01). Genotypes with Criollo germplasm were the most susceptible to lodging (∆LS = 29.4 ± 1.9), French genotypes were the most resistant to lodging (∆LS = 16.5 ± 1.6), and CIMMYT ones displayed an intermediate susceptibility (∆LS = 23.7 ± 2.3). Four French genotypes (29, 30, 34 and 35) showed ∆LS below the 10th percentile (∆LS = 8.06) and were considered to have low lodging susceptibility.

For genotypes of CIMMYT origin, higher GY was related to stems with smaller diameter, thinner wall, and lower plant height. Criollo genotypes showed differences related to GY only in plant height and spike weight. Lastly, no differences were found between high- and low-yield genotypes of French origin for any of the morphometric traits neither for lodging susceptibility (see Additional file 1).

Straw biomass composition and biogas potential production

Cell wall main components were determined to characterize straw chemical composition of the 36 genotypes during the 2017 season. Cellulose content ranged between 0.36 and 0.47 mg/mg DW, hemicellulose between 0.27 and 0.33 mg/mg DW, and lignin between 0.12 and 0.21 mg/mg DW (see Additional file 2). Lignin content showed higher variability (CV = 14.4%) than cellulose (CV = 6.5%) and hemicellulose (CV = 5.1%). Table 3 shows mean values of cellulose, hemicellulose, and lignin for the genotype groups. No significant differences were found in cellulose and hemicellulose content among groups. Interactive effects of origin and GY were detected on lignin (p = 0.03): higher lignin content was observed in high-yield genotypes from CIMMYT origin, while the opposite effect was observed for French and Criollo genotypes.

Table 3 Concentration of main components of wheat straw

To assess the performance of the wheat genotypes as bioenergy feedstock, specific biogas production was measured over the 2014 and 2017 seasons. A first-order kinetic equation (Eq. 2) was fitted to experimental data and kinetic parameters Bmax and k were determined for each individual genotype as described in “Biogas potential production and kinetic analysis” section (see Additional file 3 for individual genotype information). Mean values for the best fitting kinetic parameters, methane content in biogas, and methane yield (MY) for the genotype groups are presented in Table 4. A significant effect of the season on Bmax (p = 0.002) was observed, with average values higher for season 2017 (433.1 ± 2.8) than for 2014 (420.7 ± 2.7). Additionally, there were interactive effects of genotype and season on Bmax and k (p < 0.0001). Although the best genotypes for each kinetic parameter were different in each season, some genotypes were consistent in both seasons: genotypes 16, 18, 25, 29 and 36 exceeded the 75th percentile for Bmax (436.0 Nml/g VS for season 2014 and 447.0 Nml/g VS for season 2017), while genotypes 6, 30, 34 and 35 had k values higher than the 75th percentile (0.132 day−1 for season 2014 and 0.140 day−1 for season 2017). Noteworthy, genotypes 6 and 18 showed values above the 75th percentile for both kinetic parameters during the 2017 season.

Table 4 Kinetic parameters of biogas production, methane (CH4) percentage and methane yield

No effects of GY nor origin were detected on the kinetic parameters in either season. Methane content did not differ among genotypes or seasons and showed a mean value of 51.1% (v/v). MY (Eq. 3), was strongly correlated to SY (r = 0.96, p < 0.001) and, consequently, was affected by GY (p < 0.001) and origin (p < 0.001). Genotypes from CIMMYT origin presented lower MY than French and Criollo ones, while high-yield genotypes presented higher MY than low-yield ones, regardless of their origin. Noteworthy, despite the high correlation between MY and SY, some high SY genotypes (13, 14 and 19) were not considered high MY genotypes, since their MY values were below the 75th percentile (1880 m3/ha). Genotype 16 stood out for its high MY as a consequence of its high SY and Bmax (see Additional file 3 for individual genotype information).

To evaluate potential relationships between biogas kinetic parameters (Bmax and k) and the main components of straw biomass (cellulose, hemicellulose, and lignin), regression analyses were performed (Fig. 2). Genotypes with higher cellulose and hemicellulose content showed a tendency to have higher values of Bmax, while lignin content effect on it was not significant (Fig. 2a–c). Cellulose and hemicellulose content were negatively related to k (Fig. 2d, e), while no relationship was observed between lignin content and k.

Fig. 2

Relationship between straw main components and biogas kinetic parameters for wheat genotypes. Linear regressions between biogas kinetic parameters: Bmax (ac) and k (df), and biomass components: cellulose (a, d), hemicellulose (b, e), and lignin (c, f) were performed for the 36 wheat genotypes. Continuous lines represent the fitting of the linear regressions and inset to figures contain equations with best fitting parameters for each data set

To determine the effects of cell wall component interactions on biogas kinetic parameters, correlation analyses were performed between cellulose-to-hemicellulose, cellulose-to-lignin, and hemicellulose-to-lignin ratios and Bmax or k. As no significant correlation was found (p > 0.05, data not shown), a cluster analysis of these correlations for genotype groups according to their lignin content was performed. Regression analyses between cellulose-to-hemicellulose ratio and both kinetic parameters were performed for each lignin content, either above or below 0.15 mg/mg DW. Although the determination coefficient value was too low to explain all data variation (R2 = 0.13), a positive relationship between cellulose-to-hemicellulose ratio and Bmax was observed, regardless of the lignin content (p = 0.0347). In opposition, a negative relationship between cellulose-to-hemicellulose ratio and k (R2 = 0.51; p < 0.0001) was found for genotypes with lignin content higher than 0.15 mg/mg DW (Fig. 3).

Fig. 3

Relationship between kinetic parameters and cellulose-to-hemicellulose ratio. Biogas kinetic parameters Bmax (a) and k (b) as a function of the cellulose-to-hemicellulose ratio. Continuous line represents the fitting of a linear regression to all data (a) or to data from genotypes with lignin content higher than 0.15 mg/mg DW (b)

Target traits for ideotype design

The ideotype for food and 2G-fuel DP wheat should bear simultaneously high GY, SY, Bmax, and k and low ∆LS. In the previous sections, key traits were separately analyzed. In the present section, a principal component analysis (PCA) was performed to analyze all key traits combined (Fig. 4). Grain and straw yields, Bmax, and ∆LS accounted for principal component 1 (PC1), explaining 38.7% of the variability of the data. Principal component 2 (PC2) represented 26.6% of the variability and was mainly composed by k, Bmax, and ∆LS. A large variation in both PCs was observed among the genotypes evaluated. Low-yield genotypes (empty symbols in Fig. 4) presented the lowest values of PC1, while high-yield genotypes (black symbols in Fig. 4) were located to the right side of the bi-plot and had the highest PC1 values. Genotypes of French and Criollo origins presented the highest values of PC1. Noteworthy, no significant relation with GY and SY was found for the kinetic parameters of biogas production or lodging susceptibility in this analysis.

Fig. 4

Principal Component Analysis for key traits for the wheat dual-purpose ideotype. A Principal Component Analysis (PCA) was performed based on GY, SY, lodging susceptibility, and biogas production kinetic parameters for the 36 wheat genotypes (see Table 1 for numeric reference). Correlation of traits with PC1 were 0.913, 0.868, 0.415, − 0.415 and 0.06 for GY, SY, Bmax, ∆LS, and k, respectively. Correlation of traits with PC2 were 0.843, − 0.599, − 0.504, − 0.080 and − 0.018 for k, ∆LS, Bmax, GY, and SY, respectively. Ellipse was empirically built

A closer analysis of the PCA in combination with the results mentioned in the previous sections showed the presence of genotypes with outstanding performance in some of the traits but not in all. Among the three genotypes with high GY and very high SY (“Biomass production and allocation” section), genotype 15-Buck Baqueano produced biogas with high k and genotype 16-Buck Guapo with very high Bmax (above 90th percentile), but both displayed high lodging susceptibility (∆LS above 75th percentile). Genotype 30-SNR Nogal showed very high SY, low ∆LS, high k value and moderate Bmax (above the 50th percentile). Among the genotypes with the lowest ∆LS (“Wheat straw morphology and lodging susceptibility” section), genotype 29-Baguette 31 showed high GY, SY, and Bmax and moderate k (above the 50th percentile), genotype 34-Baguette 10 had high k but low GY and SY, and genotype 35-Baguette 21 had high k but low GY, SY and low Bmax values. Finally, among the genotypes with high values for the kinetic parameters (“Straw biomass composition and biogas potential production” section), genotypes with high k values (7-ACA 907, 26-Baguette 9 and 35-Baguette 21) did not show any other suitable trait, while among those with high Bmax, genotype 16-Buck Guapo was a high GY and SY genotype that displayed low lodging resistance and low k.

The dotted circle in Fig. 4 defines an area where acceptable values for all traits can be expected. Among them, genotypes 6-Sursem LE 2331 and 28-Baguette 19 were found to have values above the 50th percentile for all traits (or below for ∆LS). Close to this cluster, two genotypes not carrying dwarfing genes (32-Klein Atlas and 33-Klein Centauro) showed relatively good traits for DP wheat.


In this work, we have addressed genotype effects on relevant traits for the production of wheat with the dual-purpose of food (grains) and 2G-biofuel (straw); more specifically, for biogas production. We studied 36 genotypes with contrasting GY that included old and modern cultivars bearing germplasm of 3 different origins (CIMMYT, Criollo and French). In agreement with previous reports, we found high variability in all agronomic traits (GY, SY, HI, and ∆LS) [35, 38, 50, 51]. When all genotypes were regarded together, we found a positive correlation between GY and SY, in agreement with other authors [52], which suggests that it is possible to increase SY without compromising GY, a central issue for developing DP cultivars [14]. However, despite this positive correlation, high variability of SY can still be found in groups of similar GY (Table 1, [35]). These discrepancies regarding SY and GY correlation point out the importance of gaining knowledge on biomass production and allocation in wheat plants [38].

Harvest index (HI) indirectly indicates the partitioning of biomass between grain and straw [53] and a number of studies have demonstrated that HI is affected by genetic variations [52, 54]. Both SY and biomass allocation pattern differed among groups of genotypes of the same origin (Table 1 and Fig. 1). To our knowledge, this is the first time that origin effects on wheat biomass allocation are reported. In genotypes of French origin, both grain and straw yields increased at the same rate with plant size, but for genotypes of CIMMYT and Criollo origins, higher GY was attained at the expense of SY, consistently with previous reports [19, 55]. According to this, breeding for bigger plants from CIMMYT or Criollo origin do not grant higher SY, while it could be effective to increase both GY and SY in French origin genotypes.

Plant height was significantly higher for low-yield genotypes of CIMMYT and Criollo origins, while no effect of GY on plant height was detected on genotypes of French origin. A common strategy to increase GY potential has been to increase the HI through plant height reduction [21, 56]. However, in the last 30 years, a raise in GY potential without changes in plant height has been achieved in countries such as the United Kingdom [57] and Argentina [51], and even a slight increase in plant height with GY progress in modern CIMMYT cultivars has been reported [58]. Plant height has been shown to be positively correlated [52, 59] and not correlated at all [34] with SY. Here, we found that the association between plant height and SY differed among genotypes of different germplasm origin: while a negative correlation exists in genotypes of Criollo origin, no correlation was found in genotypes of the other two origins.

Morphometric traits measured at maturity were combined in a simple model to evaluate cultivar lodging susceptibility in low lodging conditions [60]. Lodging tolerance has been directly associated to plant height [61] and the risk of lodging was notably reduced after the introgression of dwarfing genes [62,63,64]. However, we found a moderate correlation between plant height and lodging susceptibility (r = 0.60, p < 0.05). Genotypes carrying French germplasm showed the lowest lodging susceptibility while Criollo ones showed the highest (Table 2). When analyzed by genotype groups, we found that correlations between height and lodging were diverse, ranging from positive (e.g. for FR-HY genotypes) to negative (e.g. for CI-HY). This observation highlights the variability of plant structure and its response to genetic factors. It also points out to the importance of complementing morphology observations with biomass composition in order to understand complex traits such as lodging. For example, for CI-HY genotypes, where there is a negative correlation between plant height and lodging susceptibility, lignin content might play a key role in avoiding lodging in plants with heavier spikes. In this sense, breeding targeted to increase capsizing resistance has led to a decrease in the saccharification potential of straw due to changes in the stem components [65]. In contrast, lodging avoidance among CR-HY genotypes with increasing plant height could be explained by thicker stems with thicker walls rather than by changes in lignin content. Finally, for FR-HY genotypes lodging avoidance strategy might be in between CIMMYT and Criollo cultivars, presenting wider stems and thicker walls in combination with a moderate increase in lignin content.

Many efforts have been done to characterize the factors affecting total biomass of modern wheat cultivars [3, 16, 66, 67]. However, less attention has been paid to factors affecting wheat straw [52]. With the development of 2G-biofuel production, efforts on increasing the yield and quality of straw to maximize the overall value of wheat crop are expected to increase [35, 68]. An analysis run on GY and HI from many trials in the United States allowed to detect differences in SY among wheat classes [33] and, here, we report genotype effects on both SY and composition, in agreement with previous reports [35, 38, 69]. However, to assess the overall potential of biogas production, both SY and its digestibility must be considered. Even though SY has been shown to be the main driver of methane yield (MY) and thus, it has been considered the main factor in feedstock selection for bioenergy production [61], straw digestibility drives the efficiency of the transformation process and determines the inputs it requires, conditioning its overall economic revenue [50, 70,71,72,73]. Straw conversion into biogas was assessed using two kinetic parameters, Bmax and k. Bmax ranged between 368–459 cm3/g VS and 392–473 cm3/g VS for 2014 and 2017 seasons, respectively (Additional File 3: Table S3), in agreement with previously reported values [74]. Although MY was strongly correlated with SY, there were some genotypes that showed high SY but low MY, due to low Bmax values (e.g. 13-BIOINTA 1000, 14–Buck 75 Aniversario and 19-Barletta 77). Conversely, some genotypes attained high MY despite not having high SY, due to their high Bmax values (e.g. 2-Don Mario Arex and 27-Baguette 18). In the same sense, genotype 16-Buck Guapo, combining high SY and high Bmax, stood out as the best in terms of MY. For biogas production, increasing Bmax reduces biomass transport, handling, storage, and pretreatment costs, while increasing k reduces feedstock pretreatments needs, and lowers retention time during anaerobic digestion, allowing reactor downsizing and the possibility of performing the anaerobic digestion at lower temperatures [40].

Biomass digestibility is often restricted by the recalcitrant nature of lignocellulose [75] and is related to lignin content and composition, cellulose crystallinity, hemicellulose content and interaction between cell wall components [45, 76,77,78] Here, straw composition ranged between 0.36 and 0.46 mg/mg DW for cellulose, 0.27 and 0.33 mg/mg DW for hemicellulose, and 0.12 and 0.21 mg/mg DW for lignin (Additional File 2: Table S2), in agreement with previous reports [75, 79, 80]. We found that biogas production potential and rate were affected by the amount of cellulose and hemicellulose in the biomass, but not by lignin content (Fig. 2). Although lignin content has been associated to lower biogas yield [81,82,83], a low correlation of biogas potential and lignin has been reported [84]. The limited digestibility of lignin under anaerobic conditions could explain these results [85]. Additionally, the efficiency of cellulose and hemicellulose degradation into simple sugars has been reported to be dependent on the polysaccharide-to-lignin ratio, among other factors [46]. We found a negative relationship between cellulose-to-hemicellulose ratio and k for genotypes with lignin content higher than 15% (w/w, Fig. 3), suggesting a lignin threshold value at which relative content of structural carbohydrates have an effect on biogas production rate. Li et al. [77] showed that lignin caused a more severe inhibition on methane production from cellulose than from hemicellulose mainly because lignin-incrusted cellulose cannot be reached by cellulases [86]. Additionally, hemicellulose has been shown to reduce lignocellulose crystallinity and, therefore, to increase biomass digestibility [45]. Hence, predominance of hemicellulose in the biomass could counteract lignin’s inhibitory effect on cellulose degradability on wheat straw with lignin contents higher than the threshold value. At lignin contents lower than this threshold, differences in Bmax and k among genotypes could respond to cellulose’s inherent recalcitrance given by its degree of polymerization and crystallinity, and not by the interaction between cellulose and lignin [48]. Accordingly, Dandikas et al. [83] found that k was not correlated with lignin in biomass with less than 4% of lignin.

We found significant genetic variability in the four traits relevant for a DP wheat ideotype. Although this should facilitate the selection for improved individual traits, trade-offs between some of them require a more complex analysis when looking for the best cultivars currently available. For instance, increasing SY and digestibility could potentially have negative impacts on GY and lodging resistance, hence, reducing the feasibility of a single ideotype [68]. Using percentile-based thresholds, we identified genotypes bearing high GY and SY (e.g. genotype 15-Buck Baqueano, 16-Buck Guapo, 29-Baguette 31 and 30-SNR Nogal). However, among them, only genotypes 29-Baguette 31 and 30-SNR Nogal displayed high lodging resistance. The use of graphical analyses such as the PCA allowed highlighting the strengths and weaknesses of each genotype. For example, when only considering MY, genotype 16-Buck Guapo would stand out as the best candidate. However, low lodging resistance (∆LS = 29.5) would preclude its production in shallow soils or high-input level conditions. It would be interesting to explore whether it is possible to increase lodging resistance of genotype 16-Buck Guapo background without significantly affecting the remaining traits. When considering all traits together, genotypes 29-Baguette 31 and 30-SNR Nogal, both of French origin, emerged as the most promising among all genotypes explored. Interestingly, genotype 6-Sursem LE 2331 showed acceptable values for all traits and might represent an interesting genotype to be used as a starting point to improve CIMMYT origin genotypes for DP use.

By including low-yield genotypes in this study, it was possible to detect genotypes that, despite presenting low GY, SY and Bmax, had good values for attributes like k or ∆LS (e.g. genotype 7-ACA 907). Also included in this study were old genotypes not carrying dwarfing genes derived from Criollo or French origin. They showed both low GY and SY values with variable ∆LS, k and Bmax. Yet, genotypes like 32-Klein Atlas and 33-Klein Centauro could be of interest in future breeding programs, since they harbor a good combination of all key traits for DP wheat. Old wheat cultivars have been proposed to be particularly promising for climate change mitigation in organic and low input systems, since they could represent a valuable resources more suitable to adapt to future scenarios of sustainable agriculture in a climate change context [87, 88].

Engineering of wheat to produce the ideal DP appears as an attractive strategy to improve the overall economic revenue of the 2G-biofuel production process. In this context, we have analyzed already available wheat cultivars for their potential use as a feedstock for anaerobic digestion without compromising GY. Results here reported can contribute to guide breeding programs meant to improve wheat as a DP crop.


There is high variability in the key traits for DP wheat ideotype for food and 2G-biofuel. None of the genotypes excels in every key trait, but a few candidates presented acceptable values in many of them. Particularly, two genotypes of French origin (29-Baguette 31 and 30-SNR Nogal) showed high GY, high-to-very high SY, low ∆LS, and moderate-to-high Bmax and k. The challenge lays in attaining a genotype with excelling values in all key traits. Breeding efforts in this direction will increase the resources available for lignocellulosic biofuel production without compromising existing food production.


Plant material

Thirty-six genotypes of wheat (Triticum aestivum L.) were grown during the 2014 and 2017 seasons. The selected genotypes are a sample of genetic materials used in the past and/or present by breeding programs in Argentina and included French, CIMMYT and local (Criollo) germplasm (French, CIMMYT and Criollo origin, respectively). Genotypes were chosen based on available information on GY and HI in order to span a wide range of variability for these traits [51, 89]. Old cultivars that do not have dwarfing genes (i.e. 19-Barletta 77; 20-Buck Naposta; 22-Klein Impacto; 23-Klein Rendidor; 24-Oleata Artillero; 32-Klein Atlas and 36-Sinvalocho) were also included to explore biomass attributes that could have been lost in more modern high-yield cultivars.

Site and crop details

Field experiments were carried out at Balcarce Agricultural Experimental Station of Instituto Nacional de Tecnología Agropecuaria (INTA), Argentina, (37° 84′ 50″ S; 55° 81′ 80″ W; 130 m.a.s.l.). Soil was a complex of fine mixed Typic Argiudoll and fine thermic Petrocalcic Paleudoll (petrocalcic horizon at 140 cm). Experiments were conducted under conventional tillage, with genotypes arranged in randomized complete blocks with two replications. Plots comprised seven rows of 5.5 m in length with a distance between rows of 0.2 m. Genotypes were sown on recommended dates and at recommended densities to achieve optimal yield according to plant lifecycle duration. Soil was fertilized with 25 kg/ha of phosphorous before sowing, and with 300 kg/ha of nitrogen, distributed in two moments, at tillering and at the beginning of the stem elongation period. Experiments were conducted under no nutrient or water limitations. Weeds, pests, and fungal diseases were chemically controlled.

Agronomic characterization

Plants were harvested at physiological maturity and dried at 60 °C for 48 h. To determine grain and straw yields, the five central rows of the plot cultivated during the 2014 season were mechanically harvested. After shoot biomass determination, spikes and grains were separated from straw and weighed. Straw dry weight was calculated as the difference between shoot biomass and grain weight. All weights were expressed on an area basis. HI was calculated as the quotient between grain weight and shoot biomass weight.

Lodging susceptibility

Lodging susceptibility was estimated as in Mirabella et al. [60]. Plant height was determined as the distance from the ground to the tip of the spike and calculated as the average of two measurements per plot. Stem diameter and stem wall thickness were measured using a digital caliper at 1 cm below the first internode. Lodging susceptibility was calculated according to the following equation:

$$\Delta {\text{LS}} = 0.87\Delta H + 27.5\Delta {\text{SW}} - 3.62\Delta J_{x} ,$$

where ∆LS represents deviation of the lodging susceptibility of the ith genotype from the mean value across all genotypes (∆LS = LSi – LS), ∆H, ∆SW and ∆Jx are the deviations of plant height (H), spike weight and the inertia moment of the area of the stem base values of individual genotypes from the mean traits, respectively. To transform ∆LS to positive values, values were corrected by adding the minimum ∆LS value obtained.

Biomass allocation

Allometric relationships were calculated as the linear regression between the natural logarithms of GY or SY and the natural logarithm of total plant biomass [90]. The slope of this relationship is called allometric constant (m) and it is the ratio of the logarithmic growth rates of the two components under study. Any value of m other than unity implies a discrepancy between these two rates [91].

Biomass analysis

Wheat straw samples were obtained by cutting ten consecutive plants above soil surface and removing the spikelets. Plants were dried at 65 °C until constant weight and then milled and sieved through a 2-mm mesh. Total solids, ash content, volatile solids (VS), and lignin content were determined according to NREL Biomass Analytical Methods [92,93,94,95]. Acid detergent fiber (ADF), neutral detergent fiber (NDF), and acid detergent lignin (ADL) were determined using the Ankom filter bag system (Ankom Technology, Fairport, NY, USA) [96]. Cellulose was calculated as the difference between ADF and ADL and hemicellulose as the difference between NDF and ADF.

Biogas potential production and kinetic analysis

Biogas potential production was assessed by performing batch anaerobic digestion experiments in 120 cm3 capacity bottles. The inoculum used consisted of equal parts (in VS) of chicken manure, pig slurry and active sludge from a potato-industry biodigester, stabilized as described by Steinmetz et al. [97]. Samples and inoculum were mixed at 1:2 ratio on a VS basis. Basal medium containing all the necessary micro- and macro-nutrients for an optimum anaerobic microbial growth was added to all bottles [98]. To maintain anaerobic conditions, the bottles were flushed with CO2 and then sealed with rubber stoppers and screwcaps. Anaerobic digestion was run at 37 °C in a water bath with daily manual stirring. Each sample was digested in duplicate. As a negative control, biogas production by the inoculum without any substrate added was measured and later subtracted from the biogas values obtained in each bottle. As a positive control, microcrystalline cellulose (Biopack) was used as a standard substrate.

The biogas volume was estimated by measuring each bottle´s inner pressure with an electronic barometer (Digital high-pressure manometer, Testo 512), and then normalized to 0 °C and 101.3 kPa. Methane content was determined by transferring 5 cm3 of the bottle’s headspace gas (v1) into an Einhorn’s saccharimeter containing 1 N NaOH with a syringe. After CO2 and H2S scrubbing in NaOH, the gas in the saccharimeter (v2) was considered to be methane. The ratio v2/v1 provided the proportion of methane (CH4%) in the biogas occupying the bottle’s headspace [98, 99]. Determinations were done after cumulative biogas reached a plateau.

A first-order kinetics model was chosen to describe biogas production over time [100], fitting the data to the exponential equation:

$$B = B_{{{\text{max}}}} \times \left( {1 - e^{ - kt} } \right),$$

where B represents cumulative biogas, Bmax is the maximum specific biogas production, k is the first-order kinetic constant (which allows to estimate biogas production rate) and t the digestion time.

Methane yield (MY) was calculated according to the following equation:

$${\text{MY}} = {\text{SY}} \times {\text{VS}} \times B_{\max } \times {\text{CH}}_{4} \% ,$$

where SY represents straw yield, VS are the volatile solids, Bmax is the maximum specific biogas production, and CH4% is the percentage of methane present in the biogas.

Data analysis

Genotypes were split into two groups according to their GY (HY: higher than 4000 kg/ha; LY: lower than 4000 kg/ha). From the combination of two yield groups and three origins (CIMMYT, Criollo and French), six genotype groups were defined as indicated in Table 1. Data were analyzed using a multifactorial ANOVA analysis [101]. Factors were germoplasm origin (CIMMYT, Criollo and French) and yield group (high and low yield). Differences between treatments were evaluated by the Fisher´s test (p < 0.05). Linear and first-order kinetic equations were fitted by non-linear simple regressions using the software Sigmaplot 10.0 [102]. Principal Component Analysis was performed using InfoStat [101].

Availability of data and materials

All data supporting the conclusions of this article are included within the article and in Additional files.



Second generation


Acid detergent fiber


Acid detergent lignin

B max :

Maximum specific biogas production


CIMMYT-high yield


CIMMYT-low yield


Criollo-High yield


Criollo-Low yield


Coefficient of variation




French-high yield


French-low yield


Grain yield


Plant height


Harvest index


Methane yield

k :

First-order kinetic constant


Neutral detergent fiber


Principal component one


Principal component two


Principal component analysis


Straw yield


Volatile solids


Lodging susceptibility


  1. 1.

    Karagöz P, Rocha IV, Özkan M, Angelidaki I. Alkaline peroxide pretreatment of rapeseed straw for enhancing bioethanol production by Same Vessel Saccharification and Co-Fermentation. Bioresour TechnoL. 2012;104:349–57.

    PubMed  Article  CAS  Google Scholar 

  2. 2.

    Chandra R, Takeuchi H, Hasegawa T. Methane production from lignocellulosic agricultural crop wastes: a review in context to second generation of biofuel production. Renew Sust Energ Revs. 2012;16:1462–76.

    CAS  Article  Google Scholar 

  3. 3.

    Sanchez-Garcia M, Royo C, Aparicio N, Martin-Sanchez JA, Alvaro F. Genetic improvement of bread wheat yield and associated traits in Spain during the 20th century. J Agric Sci. 2013;151:105–18.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Badshah M, Lam DM, Liu J, Mattiasson B. Use of an Automatic Methane Potential Test System for evaluating the biomethane potential of sugarcane bagasse after different treatments. Bioresour Technol. 2012;114:262–9.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Wang G, Gavala HN, Skiadas IV, Ahring BK. Wet explosion of wheat straw and codigestion with swine manure: effect on the methane productivity. Waste Manage. 2009;29:2830–5.

    CAS  Article  Google Scholar 

  6. 6.

    Grain and Feed Annual Country: Argentina, AR2020-0009. In: United States Department of Agriculture, Foreign Agricultural Service.

  7. 7.

    Talebnia F, Karakashev D, Angelidaki I. Production of bioethanol from wheat straw: an overview on pretreatment, hydrolysis and fermentation. Bioresour Technol. 2010;101:4744–53.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Copeland J, Turley D. National and regional supply/demand balance for agricultural straw in Great Britain. York: National Non-Food Crops Centre; 2008.

    Google Scholar 

  9. 9.

    Glithero NJ, Wilson P, Ramsden SJ. Straw use and availability for second generation biofuels in England. Biomass Bioenergy. 2013;55:311–21.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Blanco-Canqui H, Lal R. Crop residue removal impacts on soil productivity and environmental quality. CRC Crit Rev Plant Sci C. 2009;28:139–63.

    CAS  Article  Google Scholar 

  11. 11.

    Menon V, Rao M. Trends in bioconversion of lignocellulose: Biofuels, platform chemicals & biorefinery concept. Prog Energy Combust Sci. 2012;38:522–50.

    CAS  Article  Google Scholar 

  12. 12.

    Ferreira LC, Nilsen PJ, Fdz-Polanco F, Pérez-Elvira SI. Biomethane potential of wheat straw: influence of particle size, water impregnation and thermal hydrolysis. Chem Eng. 2014;242:254–9.

    CAS  Article  Google Scholar 

  13. 13.

    Townsend TJ, Ramsden SJ, Wilson P. Towards sustainable intensification of cropping systems: analysing reduced tillage practices within a bio-economic modelling framework; 2015.

  14. 14.

    Townsend TJ, Roy J, Wilson P, Tucker GA, Sparkes DL. Food and bioenergy: exploring ideotype traits of a dual-purpose wheat cultivar. Field Crops Res. 2017;201:210–21.

    Article  Google Scholar 

  15. 15.

    Slafer GA, Savin R, Sadras VO. Coarse and fine regulation of wheat yield components in response to genotype and environment. Field Crops Res. 2014;157:71–83.

    Article  Google Scholar 

  16. 16.

    Bustos DV, Hasan AK, Reynolds MP, Calderini DF. Combining high grain number and weight through a DH-population to improve grain yield potential of wheat in high-yielding environments. Field Crops Res. 2013;145:106–15.

    Article  Google Scholar 

  17. 17.

    Liu H, Searle IR, Mather DE, Able AJ, Able JA. Morphological, physiological and yield responses of durum wheat to pre-anthesis water-deficit stress are genotype-dependent. Crop Pasture Sci. 2015;66:1024–38.

    CAS  Article  Google Scholar 

  18. 18.

    Tshikunde NM, Mashilo J, Shimelis H, Odindo A. Agronomic and physiological traits, and associated quantitative trait loci (QTL) affecting yield response in wheat (Triticum aestivum L.): a review. Front Plant Sci. 2019;10:1428.

    PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Austin RB, Bingham J, Blackwell RD, Evans LT, Ford MA, Morgan CL, et al. Genetic improvements in winter wheat yields since 1900 and associated physiological changes. J Agric Sci. 1980;94:675–89.

    Article  Google Scholar 

  20. 20.

    Slafer GA, Andrade FH. Changes in physiological attributes of the dry matter economy of bread wheat (Triticum aestivum) through genetic improvement of grain yield potential at different regions of the world. Euphytica. 1991;58:37–49.

    Article  Google Scholar 

  21. 21.

    Brancourt-Hulmel M, Doussinault G, Lecomte C, Bérard P, Le Buanec B, Trottet M. Genetic improvement of agronomic traits of winter wheat cultivars released in France from 1946 to 1992. Crop Sci. 2003;43:37–45.

    Article  Google Scholar 

  22. 22.

    Crespo-Herrera LA, Crossa J, Huerta-Espino J, Vargas M, Mondal S, Velu G, et al. Genetic gains for grain yield in CIMMYT’s semi-arid wheat yield trials grown in suboptimal environments. Crop Sci. 2018;58:1890–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Beche E, Benin G, da Silva CL, Munaro LB, Marchese JA. Genetic gain in yield and changes associated with physiological traits in Brazilian wheat during the 20th century. Eur J Agron. 2014;61:49–59.

    Article  Google Scholar 

  24. 24.

    Gummadov N, Keser M, Akin B, Cakmak M, Mert Z, Taner S, et al. Genetic gains in wheat in Turkey: winter wheat for irrigated conditions. Crop J. 2015;3:507–16.

    Article  Google Scholar 

  25. 25.

    Würschum T, Langer SM, Longin CFH. Genetic control of plant height in European winter wheat cultivars. Theor Appl Genet. 2015;128:865–74.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  26. 26.

    Zhang Y, Xu W, Wang H, Dong H, Qi X, Zhao M, et al. Progress in genetic improvement of grain yield and related physiological traits of Chinese wheat in Henan Province. Field Crops Res. 2016;199:117–28.

    Article  Google Scholar 

  27. 27.

    De Vita P, Di Paolo E, Fecondo G, Di Fonzo N, Pisante M. No-tillage and conventional tillage effects on durum wheat yield, grain quality and soil moisture content in southern Italy. Soil Tillage Res. 2007;92:69–78.

    Article  Google Scholar 

  28. 28.

    Royo C, Álvaro F, Martos V, Ramdani A, Isidro J, Villegas D, et al. Genetic changes in durum wheat yield components and associated traits in Italian and Spanish varieties during the 20th century. Euphytica. 2007;155:259–70.

    Article  Google Scholar 

  29. 29.

    Sun Y, Wang X, Wang N, Chen Y, Zhang S. Changes in the yield and associated photosynthetic traits of dry-land winter wheat (Triticum aestivum L.) from the 1940s to the 2010s in Shaanxi Province of China. Field Crops Res. 2014;167:1–10.

    Article  Google Scholar 

  30. 30.

    Sparkes DL, King M. Disentangling the effects of PAR and R:FR on lodging-associated characters of wheat (Triticum aestivum). Ann App Biol. 2008;152:1–9.

    Article  Google Scholar 

  31. 31.

    Sparkes DL, Berry P, King M. Effects of shade on root characters associated with lodging in wheat (Triticum aestivum). Ann App Biol. 2008;152:389–95.

    Article  Google Scholar 

  32. 32.

    Peake AS, Bell KL, Carberry PS, Poole N, Raine SR. Vegetative nitrogen stress decreases lodging risk and increases yield of irrigated spring wheat in the subtropics. Crop Pasture Sci. 2016;67:907–20.

    Article  Google Scholar 

  33. 33.

    Dai J, Wang Z, Li M, He G, Li Q, Cao H, et al. Winter wheat grain yield and summer nitrate leaching: long-term effects of nitrogen and phosphorus rates on the Loess Plateau of China. Field Crops Res. 2016;196:180–90.

    Article  Google Scholar 

  34. 34.

    Donaldson E, Schillinger WF, Dofing SM. Straw production and grain yield relationships in winter wheat. Crop Sci. 2001;41:100–6.

    Article  Google Scholar 

  35. 35.

    Jørgensen H, van Hecke J, Zhang H, Malik PL, Felby C, Schjoerring JK. Wheat as a dual crop for biorefining: straw quality parameters and their interactions with nitrogen supply in modern elite cultivars. GCB Bioenergy. 2019;11:400–15.

    Article  CAS  Google Scholar 

  36. 36.

    García-Condado S, López-Lozano R, Panarello L, Cerrani I, Nisini L, Zucchini A, et al. Assessing lignocellulosic biomass production from crop residues in the European Union: modelling, analysis of the current scenario and drivers of interannual variability. GCB Bioenergy. 2019;11:809–31.

    Article  CAS  Google Scholar 

  37. 37.

    Unkovich M, Baldock J, Forbes M. Variability in harvest index of grain crops and potential significance for carbon accounting: examples from Australian agriculture. Adv Agron. 2010;10:173–219.

    Article  Google Scholar 

  38. 38.

    Blümmel M, Updahyay SR, Gautam N, Barma NCD, Hakim MA, Hussain M, et al. Comparative assessment of food-fodder traits in a wide range of wheat germplasm for diverse biophysical target domains in South Asia. Field Crops Res. 2019;236:68–74.

    Article  Google Scholar 

  39. 39.

    Wu Z, Zhang M, Wang L, Tu Y, Zhang J, Xie G, et al. Biomass digestibility is predominantly affected by three factors of wall polymer features distinctive in wheat accessions and rice mutants. Biotechnol Biofuels. 2013;6:183.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Pavlostathis SG, Giraldo-Gomez E. Kinetics of anaerobic treatment: a critical review. Critical reviews in environmental control, vol. 21. Routledge: Taylor & Francis; 1991. p. 411–90.

    Google Scholar 

  41. 41.

    Esposito G, Frunzo L, Liotta F, Panico A, Pirozzi F. Bio-methane potential tests to measure the biogas production from the digestion and co-digestion of complex organic substrates. Open Environ Eng J. 2012;5:1.

    CAS  Article  Google Scholar 

  42. 42.

    Dumas C, Silva Ghizzi Damasceno G, Barakat A, Carrère H, Steyer JP, Rouau X. Effects of grinding processes on anaerobic digestion of wheat straw. Ind Crops Prod. 2015;74:450–6.

    CAS  Article  Google Scholar 

  43. 43.

    Mancini G, Papirio S, Lens PNL, Esposito G. Increased biogas production from wheat straw by chemical pretreatments. Renew Energ. 2018;119:608–14.

    CAS  Article  Google Scholar 

  44. 44.

    Buffiere P, Loisel D, Bernet N, Delgenes J-P. Towards new indicators for the prediction of solid waste anaerobic digestion properties. Water Sci Technol. 2006;53:233–41.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

    Xu N, Zhang W, Ren S, Liu F, Zhao C, Liao H, et al. Hemicelluloses negatively affect lignocellulose crystallinity for high biomass digestibility under NaOH and H2SO4 pretreatments in Miscanthus. Biotechnol Biofuels. 2012;5:58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Klimiuk E, Pokój T, Budzyński W, Dubis B. Theoretical and observed biogas production from plant biomass of different fibre contents. Bioresour Technol. 2010;101:9527–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    Ma S, Wang H, Li J, Fu Y, Zhu W. Methane production performances of different compositions in lignocellulosic biomass through anaerobic digestion. Energy. 2019;189:116190.

    CAS  Article  Google Scholar 

  48. 48.

    Xu N, Liu S, Xin F, Zhou J, Jia H, Xu J, et al. Biomethane production from lignocellulose: biomass recalcitrance and its impacts on anaerobic digestion. Front Bioeng Biotechnol. 2019;7:191.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Jensen S, Oestdal H, Clausen MR, Andersen ML, Skibsted LH. Oxidative stability of whole wheat bread during storage. LWT Food Sci Technol. 2011;44:637–42.

    CAS  Article  Google Scholar 

  50. 50.

    Bellucci A, Torp AM, Bruun S, Magid J, Andersen SB, Rasmussen SK. Association mapping in scandinavian winter wheat for yield, plant height, and traits important for second-generation bioethanol production. Front Plant Sci. 2015;6:1046.

    PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Mirabella NE. Identificación de atributos asociados con la tolerancia al vuelco en trigo pan para su utilización en el mejoramiento genético. Tesis Ingeniería Agronómica, Facultad de Ciencias Agrarias, Universidad Nacional de Mar del Plata, Balcarce; 2012.

  52. 52.

    Larsen SU, Bruun S, Lindedam J. Straw yield and saccharification potential for ethanol in cereal species and wheat cultivars. Biomass Bioenergy. 2012;45:239–50.

    CAS  Article  Google Scholar 

  53. 53.

    Donald CM, Hamblin J. The biological yield and harvest index of cereals as agronomic and plant breeding criteria. Adv Agron. 1976;28:361–405.

    Article  Google Scholar 

  54. 54.

    Engel RE, Long DS, Carlson GR. Predicting straw yield of hard red spring wheat. Agron J. 2003;95:1454–60.

    Article  Google Scholar 

  55. 55.

    Shearman VJ, Sylvester-Bradley R, Scott RK, Foulkes MJ. Physiological processes associated with wheat yield progress in the UK. Crop Sci. 2005;45:175–85.

    Google Scholar 

  56. 56.

    Canevara MG, Romani M, Corbellini M, Perenzin M, Borghi B. Evolutionary trends in morphological, physiological, agronomical and qualitative traits of Triticum aestivum L. cultivars bred in Italy since 1900. Eur J Agron. 1994;3:175–85.

    Article  Google Scholar 

  57. 57.

    Berry PM, Kendall S, Rutterford Z, Orford S, Griffiths S. Historical analysis of the effects of breeding on the height of winter wheat (Triticum aestivum) and consequences for lodging. Euphytica. 2015;203:375–83.

    Article  Google Scholar 

  58. 58.

    Aisawi KAB, Reynolds MP, Singh RP, Foulkes MJ. The physiological basis of the genetic progress in yield potential of CIMMYT spring wheat cultivars from 1966 to 2009. Crop Sci. 2015;55:1749–64.

    Article  Google Scholar 

  59. 59.

    Long DS, McCallum JD. Mapping straw yield using on-combine light detection and ranging (lidar). Int J Remote Sens. 2013;34:6121–34.

    Article  Google Scholar 

  60. 60.

    Mirabella NE, Abbate PE, Alonso MP, Panelo JS, Pontaroli AC. Identifying traits at crop maturity and models for estimation of lodging susceptibility in bread wheat. Crop Pasture Sci. 2019;70:95–106.

    Article  Google Scholar 

  61. 61.

    Pinthus MJ. Lodging in wheat, barley, and oats: the phenomenon, its causes, and preventive measures. In: Advances in agronomy. New York: Academic Press; 1974. p. 209–63.

    Google Scholar 

  62. 62.

    Reitz LP, Salmon SC. Origin, history, and use of Norin 10 wheat. Crop Sci. 1968;8:686–9.

    Article  Google Scholar 

  63. 63.

    Worland T, Snape JW, Bonjean AP, Angus WJ. The world wheat book: a history of wheat breeding. Paris: Lavoisier; 2001.

    Google Scholar 

  64. 64.

    Berry PM, Sterling M, Spink JH, Baker CJ, Sylvester-Bradley R, Mooney SJ, et al. Understanding and reducing lodging in cereals. Adv Agron. 2004;84:215–69.

    Google Scholar 

  65. 65.

    Travis AJ, Murison SD, Hirst DJ, Walker KC, Chesson A. Comparison of the anatomy and degradability of straw from varieties of wheat and barley that differ in susceptibility to lodging. J Agric Sci. 1996;127:1–10.

    Article  Google Scholar 

  66. 66.

    Pedro A, Savin R, Habash DZ, Slafer GA. Physiological attributes associated with yield and stability in selected lines of a durum wheat population. Euphytica. 2011;180:195–208.

    Article  Google Scholar 

  67. 67.

    Ferrante A, Savin R, Slafer GA. Differences in yield physiology between modern, well adapted durum wheat cultivars grown under contrasting conditions. Field Crops Res. 2012;136:52–64.

    Article  Google Scholar 

  68. 68.

    Townsend TJ. Agronomic, economic and environmental analyses of dual-purpose wheat cultivars for bioenergy. Thesis . University of Nottingham; 2014.

  69. 69.

    Townsend TJ, Sparkes DL, Wilson P. Food and bioenergy: reviewing the potential of dual-purpose wheat crops. GCB Bioenergy. 2017;9:525–40.

    Article  Google Scholar 

  70. 70.

    Ercolano MR, Gomez LD, Andolfi A, Simister R, Troise C, Angelino G, et al. Residual biomass saccharification in processing tomato is affected by cultivar and nitrogen fertilization. Biomass Bioenergy. 2015;72:242–50.

    CAS  Article  Google Scholar 

  71. 71.

    van der Weijde T, Kiesel A, Iqbal Y, Muylle H, Dolstra O, Visser RGF, et al. Evaluation of Miscanthus sinensis biomass quality as feedstock for conversion into different bioenergy products. GCB Bioenergy. 2017;9:176–90.

    Article  CAS  Google Scholar 

  72. 72.

    Torres AF, Noordam-Boot CM, Dolstra O, van der Weijde T, Combes E, Dufour P, et al. Cell wall diversity in forage maize: genetic complexity and bioenergy potential. Bioenergy Res. 2015;8:187–202.

    Article  Google Scholar 

  73. 73.

    Mahmood A, Ullah H, Ijaz M, Javaid MM, Shahzad AN, Honermeier B. Evaluation of sorghum hybrids for biomass and biogas production. Aust J Crop Sci. 2013;7:1456.

    CAS  Google Scholar 

  74. 74.

    Rajput AA, Zeshan, Visvanathan C. Effect of thermal pretreatment on chemical composition, physical structure and biogas production kinetics of wheat straw. J Environ Manage. 2018;221:45–52.

    CAS  PubMed  Article  Google Scholar 

  75. 75.

    Chen Y, Zhao Z, Zou H, Yang H, Sun T, Li M, et al. Digestive performance of sludge with different crop straws in mesophilic anaerobic digestion. Bioresour Technol. 2019;289:121595.

    CAS  PubMed  Article  Google Scholar 

  76. 76.

    Yoshida M, Liu Y, Uchida S, Kawarada K, Ukagami Y, Ichinose H, et al. Effects of cellulose crystallinity, hemicellulose, and lignin on the enzymatic hydrolysis of miscanthus sinensis to monosaccharides. Biosci Biotechnol Biochem. 2008;72:805–10.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  77. 77.

    Li W, Khalid H, Zhu Z, Zhang R, Liu G, Chen C, et al. Methane production through anaerobic digestion: participation and digestion characteristics of cellulose, hemicellulose and lignin. Appl Energy. 2018;226:1219–28.

    CAS  Article  Google Scholar 

  78. 78.

    Méchin V, Argillier O, Rocher F, Hébert Y, Mila I, Pollet B, et al. In search of a maize ideotype for cell wall enzymatic degradability using histological and biochemical lignin characterization. J Agric Food Chem. 2005;53:5872–81.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  79. 79.

    Ali SS, Khan M, Mullins E, Doohan F. The effect of wheat genotype on ethanol production from straw and the implications for multifunctional crop breeding. Biomass Bioenergy. 2012;42:1–9.

    CAS  Article  Google Scholar 

  80. 80.

    Fitria, Ruan H, Fransen SC, Carter AH, Tao H, Yang B. Selecting winter wheat straw for cellulosic ethanol production in the Pacific Northwest, USA. Biomass Bioenergy. 2019;123:59–69.

    CAS  Article  Google Scholar 

  81. 81.

    Triolo JM, Sommer SG, Møller HB, Weisbjerg MR, Jiang XY. A new algorithm to characterize biodegradability of biomass during anaerobic digestion: Influence of lignin concentration on methane production potential. Bioresour Technol. 2011;102:9395–402.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  82. 82.

    Monlau F, Sambusiti C, Barakat A, Guo XM, Latrille E, Trably E, et al. Predictive models of biohydrogen and biomethane production based on the compositional and structural features of lignocellulosic materials. Environ Sci Technol. 2012;46:12217–25.

    CAS  PubMed  Article  Google Scholar 

  83. 83.

    Dandikas V, Heuwinkel H, Lichti F, Eckl T, Drewes JE, Koch K. Correlation between hydrolysis rate constant and chemical composition of energy crops. Renew Energ. 2018;118:34–42.

    CAS  Article  Google Scholar 

  84. 84.

    Thomas HL, Arnoult S, Brancourt-Hulmel M, Carrère H. Methane production variability according to miscanthus genotype and alkaline pretreatments at high solid content. Bioenerg Res. 2019;12:325–37.

    CAS  Article  Google Scholar 

  85. 85.

    Koch K, Lübken M, Gehring T, Wichern M, Horn H. Biogas from grass silage—measurements and modeling with ADM1. Bioresour Technol. 2010;101:8158–65.

    CAS  PubMed  Article  Google Scholar 

  86. 86.

    Jördening H-J, Winter J. Environmental biotechnology: concepts and applications. Hoboken: John Wiley & Sons; 2005.

    Google Scholar 

  87. 87.

    Ceccarelli S, Grando S, Tutwiler R, Baha J, Martini AM, Salahieh H, et al. A methodological study on participatory barley breeding I. Selection phase Euphytica. 2000;111:91–104.

    Article  Google Scholar 

  88. 88.

    Carranza-Gallego G, Guzmán GI, García-Ruíz R, González de Molina M, Aguilera E. Contribution of old wheat varieties to climate change mitigation under contrasting managements and rainfed Mediterranean conditions. J Clean Prod Jr. 2018;195:111–21.

    Article  Google Scholar 

  89. 89.

    Abbate PE, Pontaroli AC, Lorenzo M, Retamar JO, Castaño MC, Toledo JI, et al. Red de evaluación de cultivares de trigo pan (RET): Resultados obtenidos en INTA Balcarce sin y con funguicida durante la campaña 2010/11. INTA, EEA Balcarce, AR; 2011.

  90. 90.

    Coleman JS, McConnaughay KD, Ackerly DD. Interpreting phenotypic variation in plants. Ecol Evol. 1994;9:187–91.

    CAS  Article  Google Scholar 

  91. 91.

    Hunt R. Plant growth analysis studies on biology, vol. 96. London: Edward Arnold Publishers Ltd.; 1978.

    Google Scholar 

  92. 92.

    Hames B, Ruiz R, Scarlata C, Sluiter A, Sluiter J, Templeton D. Preparation of samples for compositional analysis. Lab Anal Procedure (LAP). 2008;1617:65–71.

    Google Scholar 

  93. 93.

    Sluiter A, Hames B, Hyman D, Payne C, Ruiz R, Scarlata C, et al. Determination of total solids in biomass and total dissolved solids in liquid process samples, vol. 9. Citeseer: National Renewable Energy Laboratory; 2008.

    Google Scholar 

  94. 94.

    Sluiter A, Hames B, Ruiz R, Scarlata C, Sluiter J, Templeton D. Determination o f ash in biomass. Technical Report NREL/TP–510-42622. National Renewabl e Energy Laboratory; 2008.

  95. 95.

    Sluiter A, Hames B, Ruiz R, Scarlata C, Sluiter J, Templeton D, et al. Determination of structural carbohydrates and lignin in biomass. Laboratory analytical procedure, vol. 1617. Golden: National Renewable Energy Laboratory; 2008. p. 1–16.

    Google Scholar 

  96. 96.

    Dumas JBA. Procedes de I’analyse organique. Ann Chim Phys. 1831;47:198–205.

    Google Scholar 

  97. 97.

    Steinmetz RLR, Mezzari MP, da Silva MLB, Kunz A, do Amaral AC, Tapparo DC, et al. Enrichment and acclimation of an anaerobic mesophilic microorganism’s inoculum for standardization of BMP assays. Bioresour Technol. 2016;219:21–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  98. 98.

    Angelidaki I, Alves M, Bolzonella D, Borzacconi L, Campos JL, Guwy AJ, et al. Defining the biomethane potential (BMP) of solid organic wastes and energy crops: a proposed protocol for batch assays. Water Sci Technol. 2009;59:927–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  99. 99.

    Ergüder TH, Güven E, Demirer GN. Anaerobic treatment of olive mill wastes in batch reactors. Process Biochem. 2000;36:243–8.

    Article  Google Scholar 

  100. 100.

    Tong X, Smith LH, McCarty PL. Methane fermentation of selected lignocellulosic materials. Biomass. 1990;21:239–55.

    CAS  Article  Google Scholar 

  101. 101.

    Di Rienzo JA, Casanoves F, Balzarini MG, González L, Tablada M, Robledo y C. InfoStat versión 2011. Grupo InfoStat, FCA, Universidad Nacional de Córdoba, Argentina; 2011; 8:195–9.

  102. 102.

    Inc SS. SigmaPlot for windows (development and testing procedures): Version 11.0. Bochum: Triestram & Partner GmbH Bochum; 2008.

    Google Scholar 

Download references


Maria Mercedes Echarte, Paola Mónica Talia, Ignacio Durruty, Ana Clara Pontaroli acknowledge CONICET as career research members. Nadia Gabbanelli, Elisa Erbetta and Maria Eugenia Sanz Smachetti acknowledge the CONICET fellowships.


This work was supported by Instituto Nacional de Tecnología Agropecuaria (INTA, PD I149), the Research Council of Argentina (CONICET), the Research Council of Buenos Aires Province, the National University of Mar del Plata and Agencia Nacional de Promoción Científica y Tecnológica (ANPCyT) Proyectos de Investigación Científica y Tecnológica (PICT) 2018 No. 4149. Nadia Gabbanelli, Elisa Erbetta and Maria Eugenia Sanz Smachetti hold scholarships from CONICET. This work is part of the thesis by Nadia Gabbanelli (Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Mar del Plata, Argentina) and Elisa Erbetta (Facultad Ciencias Agrarias, Universidad Nacional de Mar del Plata, Argentina) in their partial fulfillment of their requirements for the Doctor’s degree.

Author information




ML, IR, MV performed experiments at the field (crop sowing, management and harvest) and agronomical measurements. NG, EE, MV, ID and MME carried out chemical analysis and biogas production assays. NG, EE, MESS PMT, ID, ACP and MME performed data analysis, discussion of results and paper elaboration. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Maria Mercedes Echarte.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors consent for publications.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Morphometric characteristics and lodging susceptibility of the wheat genotypes.

Additional file 2.

Concentration of main components of wheat straw.

Additional file 3.

Kinetic parameters of biogas production and methane (CH4) percentage for wheat genotypes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gabbanelli, N., Erbetta, E., Sanz Smachetti, M.E. et al. Towards an ideotype for food-fuel dual-purpose wheat in Argentina with focus on biogas production. Biotechnol Biofuels 14, 85 (2021).

Download citation


  • 2G biofuels
  • Wheat straw
  • Anaerobic digestion
  • Straw quality