A set of autotetraploid [polyploid Energo (PP-E); 2n = 4x = 76] genotypes was generated by colchicine treatment of axillary buds from Energo plants (Salix viminalis var. Energo diploid cultivar). The detailed characterization of these genotypes in comparison to the Energo plants was published earlier (Dudits et al., 2016). Using tetraploid plants, a crossbreeding program was initiated to produce triploid genotypes. As greenhouse crossings, the receptive female catkins of commercial energy willow cultivars (i.e., Inger and Tordis) were hand-pollinated with pollen grains collected from male flowers of the autotetraploid plants (PP-E7 and PP-E15). After crossing, hybrid seeds were germinated as in vitro cultures, and their genome size was determined by flow cytometry of nuclei isolated from their root tips. Vegetative propagation of the selected triploid hybrid lines was carried out in a breeding garden by using cuttings of woody stems.

In the spring of 2019, woody stems of the 2-year-old parental and triploid hybrid plants were cut back, and the daily growth rates of their newly developing shoots were determined over a 37-day-period from five plants by measuring three shoots per plant. Ten plants per genotype were grown in randomized plots of our breeding garden Szeged, Hungary [location I: global positioning system (GPS) coordinates: 46°15′10.8′′N, 20°8′53.66′′E]. Next spring in 2020, woody stems of the 1-year-old parental and triploid hybrid plants were cut back, and the daily growth rates of their newly developing shoots in 4 days were determined from 10 plants (measuring three shoots per plant). Twenty plants per genotype were grown in randomized plots of the breeding garden in Kiskunhalas, Hungary (location II. GPS coordinates: 46° 26′ 7.5876′′ N and 19° 29′ 0.3408′′ E). To characterize hybrid vigor, we calculated the mid-parent heterosis (MPH%) as the heterosis over mid-parent (MP%) = [(F1-MP)/MP × 100], where F1 is the numerical value trait measurement in the hybrid and MP values are the mean values of the parents (P1 + P2)/2. In addition, we present heterosis values relative to the cultivar parents (CPH%) = [(F1-CP)/CP × 100] for the selected traits.

Plants were annually harvested in winter during the dormant phase. We measured the woody biomass of 60 cm-long cuttings and collected woody stem samples for their chemical analysis of lignocellulose compositions and biogas fermentation potential.

Willow plants were grown in the soil of the complex stress diagnostic system (Cseri et al., 2020) under greenhouse conditions (at 21°C and illuminated with ∼400 μmol photons m^–2 s^–1). The analyzed plants were grown under an optimal water supply (60% soil water content) during the whole life cycle, and their water consumption data were stored automatically by the computer.

To measure the photosynthetic functioning, both stomatal conductance (g_sw; mol H₂O m^–2 s^–1) and CO₂ assimilation rate (A_N; μmol CO₂ m^–2 s^–1) were monitored in young, fully expanded leaves with a portable photosynthesis system (LI-6400, LI-COR, Inc., Lincoln, NE, United States) and an atmospheric CO₂ source, according to Poór et al. (2011). The CO₂ concentration was maintained at 400 ppm by a soda-lime reagent (Sigma-Aldrich ACS reagent). On each leaf, 2 cm² of its area was measured with the controlled CO₂ flow (300 μmol mol^–1) at 25°C under greenhouse conditions.

The carbohydrate (soluble sugars and starch) concentration was determined according to Poór et al. (2011). Briefly, 100 mg [fresh weight (FW)] of leaf material was ground in liquid nitrogen, from which soluble sugars were extracted with 1 ml of 80% ethanol at 80°C for 30 min. The homogenate was centrifuged at 2,600 g for 10 min. After a second extraction step, the combined supernatants were used to determine the soluble sugar content at 630 nm after its reaction with anthrone (Normapur, VWR Int., Leuven, Belgium) dissolved in 72% sulfuric acid, where glucose (Normapur, VWR Int., Leuven, Belgium) dissolved in 80% ethanol was then used as the standard. The remaining pellet was washed with 1 ml of deionized water, then hydrolyzed with 1 ml of 1.1% HCl at 100°C for 30 min and centrifuged for 10 min at 2,600 g. The starch concentration was also evaluated spectrophotometrically at 630 nm with anthrone reagent, and for this purpose starch (Normapur, VWR Int., Leuven, Belgium) dissolved in 1.1% HCl served as the standard.

Woody stem cuttings (each 10 cm in length) of triploid hybrid and parental genotypes were placed into holes of a plastic foam floating in a pool filled with the Knop-solution. The dormant buds were activated, and the outgrowing shoots grew in a thermostat room at 21°C under artificial light. The level of illumination was ∼200 μmol photons m^–2s^–1. To determine hormone contents (pmol g^–1 FW) in the shoot tip meristems, tissues samples (30–80 mg fresh weight) were purified and analyzed by following the methodology of Dobrev and Kamínek (2002) and Dobrev and Vankova (2012) as described in our previous publication (Dudits et al., 2016). The shoot tip tissues were collected from two replicates of plants whose mean value is provided.

For the determination of fiber composition on a dry weight basis, the wood samples were dried at 55°C in an air dryer for 12–16 h. After drying, each sample was ground to pass through a 1-mm forage mill. Analysis of fibers was then carried out following the protocols described by Van Soest and Robertson (1985). The acid detergent fiber (ADF) was determined gravimetrically, as the amount of residue remaining upon ignition after receiving the 72% H₂SO₄ treatment. Hemicellulose% was estimated as neutral detergent fiber (%)—ADF (%). Cellulose (%) was estimated as ADF (%)—lignin (%). Three biological duplicates of wood samples per genotype were used in each fiber analysis.

The batch experiments were carried out in 160-ml glass reactors in triplicates. The milled wood particles were first inoculated with a fresh sample from an industrial scale mesophilic biogas plant, then fed with pig slurry and maize silage mix (Zöldforrás Biogas Plant, Szeged, Hungary), and finally filtered through a 2-mm mesh. The CH₄ content was measured by an Agilent 6,890 N gas chromatograph (Agilent Technologies, Santa Clara, United States). Triplicate batch digesters were assembled for each combination of substrates having a substrate–inoculum ratio (S/I) of 0.5, for which the reaction temperature was set at 37°C. Full details of the fermentation protocol are provided in our recent article (Kakuk et al., 2021).

Prior to statistical analyses, statistically significant outliers (if there were any) were identified by the extreme studentized deviate (ESO) method using the graph pad online outlier calculator. The significant outliers were removed from the dataset if the normal distribution was confirmed without the outliers detected by the ESO test. Statistical analyses were implemented by SigmaPlot 12.0 software (Systat Software Inc.), the significant differences among genotypes were evaluated by using Duncan’s multiple range test after one-way ANOVA if all the preconditions for the test were met. If any of the preconditions were not met, then the Student–Newman–Keuls test (SNK test) or Dunn’s multiple comparisons test was used after ANOVA on ranks as indicated in the figure captions. Differences were considered significant if p ≤ 0.05.

In the greenhouse, the receptive female catkins of the Inger and Tordis cultivars were hand-pollinated with pollens collected from male flowers of Energo tetraploid plants (PP-E 7 and PP-E15). Embryos were excised and further cultured in vitro, and outgrowing plantlets were propagated in agar cultures. The hybrid and parental plants were first planted into the soil of pots and grown in the greenhouse. Cuttings of woody stems were rooted in water and used for nuclei isolation from their root tips and characterized by flow cytometry. Based on these DNA content indicators, we were able to establish a collection of triploid hybrid (TH) lines, all of which were propagated by the cutting of woody stems in an outside nursery during the study years.

In the third year of the willow plantation, woody stems were cut back in spring. Their newly developed shoots were used to characterize the genotype-dependent early growth and stem parameters of willow plants, and this was done between May 5 and June 11, 2020 (Table 1, location I). Depending on the crossing combinations, significant differences in shoot length and growth rate could be recorded between the parental plants and THs. Both Swedish cultivars (Inger, Tordis) and the two tetraploid lines (PP-E7 and PP-E15) showed differences in shoot length and daily growth rate in June 2020. Crossing the fast-growing Inger with the tetraploid PP-E7 plants led to reduced growth in THs (TH21/2), exhibiting increased growth intensity. In this crossing combination, in June 2020, we detected an MPH value of 37.88% for shoot length and one of 42.96% for growth rate. A similar hybridization effect was recorded in the TH17/17 hybrids with heterosis for shoot length (43.46% MPH) and growth rate (39.57% MPH) values in the June measurements. Combining the two parental genotypes that have low growth rates (cultivar Tordis and PP-E7) generated considerable hybrid vigor in shoot length (54.73% MPH) and growth rate (60.55% MPH) in the TH3/12 plants sampled in June 2020. Comparing the MPH percentages determined at the two measurement time points, the heterosis of shoot length evidently was decreased in the TH17/17 genotype, while MPH% values increased in plants of the TH3/12 and the TH21/2 lines. This pattern suggested the possibility of increasing the expression of heterosis in the later period vegetative growth for these two genotypes. The diameter values in the growing phase of the green stems differed significantly between the two cultivars, and the autotetraploid genomes (PP-E7 and PP-E15) had developed wider stems (Table 1). For this trait, triploid heterosis was expressed in all three hybrid genotypes as follows: TH3/12: MPH: 38.98%; TH17/17: MPH: 42.34%, and TH21/2: MPH:29.71%.

In 2021, the same growth characteristics were recorded in late April by examining shoots developed in the second year of the willow plantation (Table 1, location II). Considerable hybrid vigor in plant height and growth rate was evident in all three crossing combinations under these distinct conditions. Crossing the Tordis plants with PP-E7 autotetraploid resulted in the highest MPH for shoot length. This was found in TH3/12 plants (MPH:68.85%), and likewise for daily growth rate (MPH: 97.18%). In general, these second-year data collected under better field conditions indicate that triploid hybrid vigor in terms of growth rate was higher than in the preceding year, except plants in TH21/2.

Early growth vigor of triploid hybrid willow plants could be recorded under field conditions. As shown in Figure 1, the TH3/12 plants produced substantially more biomass than did the parental plants (Tordis, PP-E7). Although the TH21/2 triploid hybrid plants accrued a slightly greater amount of the green biomass than did the leading Swedish cultivar (Inger), they were nonetheless considerably more productive than their tetraploid parental plants (PP-E7).

On searching for the physiological basis of detected hybrid vigor, leaf functions are considered suitable primary targets. Differences in leaf sizes can also be expressed by weight data (Figure 2). Doubling the genome size resulted in the largest leaves, which could have contributed to the lack of MPH%, except in the TH21/2 hybrid. When the THs were compared with the commercial cultivars (Inger, Tordis), the CPH% indicated significant levels of hybrid vigor for all THs but especially for the TH21/2 plants.

The length and width of leaves differed significantly between the parental and TH21/2 plants (Figure 3). The two parents produced leaves with different morphologies, with shortened but wider leaves characteristic of tetraploid plants (PP-E 7). By contrast, in the triploid hybrid TH21/2 plants this is clearly a feature of the combination of long and wide leaf characteristics.

It was shown previously that a positive correlation between photosynthetic rate and biomass production could only be detected in a restricted set of willow genotypes (Andralojc et al., 2014). Still, willow breeding programs could consider photosynthetic capacity as a promising target to increase biomass production. Hence, we tested the expression of hybrid vigor in traits related to that (Table 2).

In agreement with our study, it was found that tetraploid plants were capable of higher photosynthetic efficiency (Dudits et al., 2016), and in the present greenhouse study, both PP-E7 and PP-E15 plants were capable of high net photosynthetic CO₂ uptake. Based on this a low MPH was recorded for the TH3/12 and TH21/2 hybrid plants. The stomatal conductance values indicated a higher evaporation rate than in their parental plants. Considerable hybrid vigor could be seen in these crossing combinations. Table 2 also provides pertinent information on water requirements, as inferred from stomatal conductance values of willow plants with different genetic backgrounds. The two diploid cultivars exhibited relatively reduced levels of transpiration. Significant differences between the THs raise the possibility for selecting drought-tolerant energy willow genotypes. The highest positive MPH value for stomatal conductance was measured in plants of the TH3/12 hybrid. The TH21/2 hybrids expressed positive heterosis in this trait, yet low heterosis was found in water-use efficiency (MPH:0.56%; CPH:-1.79%). As presented in Table 2, the autotetraploid energy willow parental plants (PP-E15) exhibited the highest stomatal conductance and water use efficiency values. These parameters have contributed to the negative heterosis of stomatal conductance in the TH17/17 hybrid plants.

The primary products of photosynthesis are soluble carbohydrates which form starch, which is an insoluble and non-structural carbohydrate. The present set of energy willow genotypes let us analyze hybrid vigor in terms of the production of these compounds (Table 2). Sucrose concentration showed positive heterosis in the leaves of TH21/2 hybrids (MPH:25.31%; CPH:46.40%). Limited heterosis could be detected for starch accumulation, namely in leaves from the TH3/12 and TH21/2 hybrid plants.

In short-rotation energy plants, the growth rate and later biomass formation are the key parameters that are governed by the hormonal status of the shoots. In this study, stem cuttings from the THs and their parents were cultured in the Knop-solution, and the outgrowing shoots were used to quantify a large set of plant growth regulators. In Tables 3A–C the concentrations of hormones that indicated genotype-dependent, variable MPH (MPH%) in shoots of different triploid hybrid plants when compared to their parents.

From the view of hybrid vigor, positive MPH could be detected for jasmonic acid (MPH: 91.91, 67.08, 29.73%)

in all three crossing combinations. Regarding the hormone composition in young shoots, the TH21/2 triploid hybrid plants differed significantly from the two other hybrids (TH17/17 and TH3/12) in having a high concentration of auxin-type hormones (Table 3C). In the TH21/2 shoots, levels of significant positive MPH were found in indole-3-acetic acid (MPH: 207.49%), phenylacetic acid (MPH:223.51%), salicylic acid (27.72%), and benzoic acid (85.75%). By contrast, expression of MPH for cytokinin-related compounds reached higher values in the shoots of TH3/12 and TH17/17 hybrid plants, especially in TH3/12 in comparison with the TH21/2 plants; for the TH3/12 shoots: isopentenyladenosine (43.95%), cis-zeatin riboside (72.42%), trans-zeatin riboside-O-glucoside (86.97%), and dihydrozeatin riboside (158.29%) and in TH17/17 shoots: isopentenyladenosine (20.90%) and trans-zeatin riboside monophosphate (500.00%, Tables 3A,B). Despite the complexity of young shoots’ hormonal status, the presented snapshot pointed to some key regulators involved in MPH.

The efficiency of anaerobic digestion (AD) of lignocellulose for biomethane production is directly determined by the compositional properties of wood biomass (recently reviewed by Xu et al., 2019). Accordingly, we also analyzed the possible involvement of hybrid vigor in determining wood tissue weight (g/cm) and the ratios of the most significant components of the woody cell wall (Table 4). Relatively low MPH was expressed for wood weight after crossing Inger and the PP-E7 parents (MPH:8.98%). Lignin is the most effective factor limiting the biodegradability of lignocellulose in the three-dimensional network inside the cell wall (Thomsen et al., 2016). Of the three THs, lower lignin contents were present in wood tissues from stems of the TH21/2 and TH3/12 genotypes, with negative MPH values (–14.36 and –8.01%, respectively). In stark contrast, the analyzed data revealed a positive MPH value (15.59%) for acid detergent lignin content from the wood of the TH17/17 hybrid. In the AD process, the amounts and polymerization degree of the cellulose component also play a central role. As evinced by Table 4, portions of cellulose in wood samples from the TH hybrids were intermediate between the cultivars and tetraploid parents. This indicated a significant increase in the hemicellulose content of wood samples from the TH3/12 THs. Conceivably, this could influence methane fermentation efficiency (refer below).

Recent work by Kakuk et al. (2021) and Nyári et al. (2021) reported that the efficiency of using both woody and green biomass of willow as raw material for the anaerobic digestion (AD) and biomethane production depended on the genome size of willow plants. With that in mind, here, we evaluated whether the expression of triploid hybrid vigor in willow biomass-related traits could influence the biogas yield. As seen in Figure 4, methane yields when using woody tissues from the Swedish cultivars were the same, whereas using wood samples from the tetraploid genotypes led to moderate differences in the fermentation responses. For methane yield as a trait, the MPH varied from –27.40 to 27.80%, and a significantly positive value was detected in the fermentation of wood tissues from the TH3/12 THs. Biogas yield from the TH21/2 triploid wood also surpassed that of its parental plants but with lower heterosis. Hybrid vigor in terms of biomethane production from TH3/12 wood samples was detected during the whole fermentation cycle as shown by its cumulative values (Figure 5). The negative heterosis for methane yield from the wood samples of the TH17/17 hybrids originated from specific wood structural attributes, such as its higher lignin content (Table 4).

Generation of hybrid vigor by crossing plants coupled with multiplying their genome size via polyploidy offers an efficient approach to optimize traits of forest trees, especially those that are grown in short-rotation systems. The wider use of Salix species for efficient bioenergy production depends on sufficient biomass productivity and promising techniques in fermentation. Biologically, there is a need to increase the genetic capacities for both biomass yield and the optimal composition of raw material. Together, these efforts could contribute to addressing global issues, such as climate change or environmental pollution. Previously, we reported on several traits of autotetraploid willow genotypes as their larger leaf and root systems, greater net photosynthetic CO₂ uptake, improved photosynthetic functions, slower primary growth, and increased shoot diameter could jointly determine biomass yield (Dudits et al., 2016). In the present study, these tetraploid plants were crossed with two leading diploid Swedish cultivars to exploit positive effects arising from both hybrid vigor and polyploidy. This experimental approach demonstrates a novel way forward characterizing the interaction between hybrid vigor and polyploidy, which has been extensively studied for alloploid plants derived from hybridization between different species followed by chromosome doubling events (refer to review by Chen, 2010). Study by Carlson and Smart (2021) reported greater total above-ground dry biomass for a suite of interspecific Salix THs in comparison with either diploid or tetraploid parents.

The present crossings with different parents resulted in a variable degree of MPH% in biomass-related traits in the THs. In all crossing combinations, examined here, considerable MPH values were recorded for early shoot length and daily growth rate under field or laboratory conditions (Tables 1, 3A–C). Despite the differing MPH% values between genotypes characterized in two growing seasons and at two locations, the ranking hybrid vigor of various crossings was similar across the different experimental conditions. In the Knop-solution as an optimal condition, tetraploid parents (PP-E7 and PP-E15) grew taller shoots than did triploid hybrid plants with reduced MPH% (Tables 3A–C). Our findings are in agreement with a study that found that environmental factors are capable of modulating the manifestation of hybrid vigor (Li et al., 2018). When evaluating the relationship between daily growth rate and actual plant height, stem diameter is a useful variable as noted for other woody species (Gleason et al., 2018). As Table 1 shows, plants of the Inger cultivar with the highest shoot growth rate attained maximum plant height and stem diameter. At the same time, the MPH% in its hybrid TH21/2 was relatively low. By contrast, the plants of the TH3/12 hybrids generated by crossing Tordis with PP-E7 willow were found to express significant hybrid vigor relative to the parental values.

Larger leaf size can be considered a significant contributing component to heterosis (Liu et al., 2020). In the present study, MPH for leaf size was detectable only in plants from the TH21/2 triploid hybrid (MPH:33.99%). However, when leaf weight values from all the three hybrids were compared with their cultivar parents, the CPH values provided clear evidence of heterosis (TH3/12:12.62%; TH17/17:36.00%; TH21/2:67.20%). This greater leaf area enabled a greater production of photosynthates. Further, the triploid willow hybrids featured a more efficient CO₂ assimilation rate, especially relative to the Tordis parents (Table 2). Metabolome and proteome studies on maize hybrids and their inbred parents have shown positive MPH in the photosynthetic pathway (Li et al., 2020). Our study of willow clearly shows heterotic responses of different degrees in the CO₂ assimilation rate, as well as in sucrose concentration and of starch accumulation in leaves of various triploid hybrid genotypes. In parallel, our data did not uncover direct relationships between these parameters. This finding might be due to a technical limitation caused by the different timing of CO₂ assimilation rate measurements vs. the collection of leaf samples. In rice, heterotic responses were found to be related to the circadian rhythm pathway (Shen et al., 2015).

Water shortage and drought are the major concerns for plants generated by ongoing climate changes. Therefore, the water-use parameters of crop plants, including those trees used in agroforestry, can determine their end-use capacity for bioenergy production. Stomatal conductance values indicated positive MPH for the TH3/12 (72.04%) and TH21/2 (43.27%) hybrids. These higher evaporation rates could increase the drought sensitivity of these THs that might be compensated in part by an enlarged root system (unpublished data). For example, water use efficiency displayed a positive MPH in plants of TH3/12 genotypes, whereas a negative MPH was detected in TH21/2 plants. These water-use patterns differed markedly from the TH17/17 genotype, these hybrid plants evaporated water to the same extent as the Tordis parental plants, and this water loss was lower than that from the tetraploid PP-E 15 parental plants. These factors resulted in a negative MPH% for the trait of water use efficiency. This detected variability in water utilization provides the basis for further breeding, which is essential for overcoming these limitations to enable the wide use of short-rotation systems even under a changing climate. We should mention that higher water usage could be a positive factor in optimizing wastewater treatment efficiency of agroforestry (refer to review by Dimitriou et al., 2009 and study by Amiot et al., 2020).

In searching for a hormonal basis to triploid hybrid vigor, our results show that the expression of heterosis may be linked with increased concentrations of several hormones in the shoot meristems of THs arising from various crossing combinations. As a novel finding, all three THs expressed MPH for jasmonic acid accumulation (TH3/12: MPH:91.91; TH17/17: MPH:67.08; TH21/2: MPH:29.73%). In Table 3, the MPH value (102.19%) was significant for the jasmonate isoleucine level in the TH3/12 plants. Presently, we lack an explanation for the detected jasmonic acid response in shoots of triploid willow hybrids, especially if we consider results from proteomic and transcriptomic analysis of seedling leaf tissues of the maize hybrid, B73xMo17, and its inbred parents (Birdseye et al., 2021). These authors reported repression of jasmonic acid biosynthetic enzyme levels in that maize hybrid. Our data may therefore suggest a potential link to biotic stress tolerance in willow. The hormone results here highlight hybrid heterosis in the accumulation of another stress-related hormone, abscisic acid (ABA), and in addition to its physiologically inactive metabolites such as dihydrophaseic acid, neophaseic acid, and ABA–glucose ester. Intensive ABA metabolism was detected only in plants of the TH17/17 willow hybrids, these showing a 33% MPH for shoot height under in vitro growth conditions (Table 3B). Detection of ABA as positive factors in heterosis is an unexpected finding that might partly explain the negative MPH% in stomatal conductance and lower CO₂ assimilation rate (refer to Table 2). The essential role of ABA in the movement of stomata is well-established (reviewed by Vishwakarma et al., 2017). In contrast, in stems of intraspecific hybrids of balsam poplar (Populus balsamifera L.), the ABA content was negatively correlated with heterosis in stem volume as well as stem dry biomass (Hu and Thomas, 2019). Ectopic expression of an ABA inactivation gene (ZmABA8ox1b) enhanced plant vigor in the early stages of maize seed germination in its hybrid B73/Mo17 relative to its parental inbred lines (Li et al., 2016).

Triploid hybrid plants (TH21/2) produced by the crossing of Inger plants and PP-E7 tetraploid willows share a unique hormone composition (Table 3C). In these plants, considerable hybrid vigor was expressed in natural auxins (indole-3-acetic acid: MPH207.49%; phenylacetic acid: MPH-223.51%) that was not detectable in plants of the other two willow hybrids. The indole-3-acetic acid content of stem of intraspecific hybrids of balsam poplar (P. balsamifera) showed a strong positive correlation with heterosis in stem volume and stem dry biomass (Hu and Thomas, 2019). Higher auxin concentrations were detected in the seedling of canola (Brassica napus L.) hybrids (Zhu et al., 2020). Considering the potential roles of benzoic acid and salicylic acid in heterosis, we obtained contradictory results, in that we detected MPH for concentrations of these compounds (salicylic acid: 27.72%; benzoic acid:85.75%) as part of the heterosis response in willow. By contrast, in Arabidopsis F1 hybrids, an increasing salicylic acid concentration diminished their heterosis (Groszmann et al., 2015). Considering the multiple roles of cytokinins in plant developmental processes, including regulation of cell division and differentiation (by Kieber and Schaller, 2018), hybrid vigor in the concentrations of various cytokinin metabolites (Table 3) could have been responsible for the physiological responses that we detected in the willow THs. Apart from cytokinin glucosides, high levels of cytokinin ribosides were recorded in the three triploid hybrid combinations (TH3/12: 24.4, TH17/17: 21.02, and TH21/2: 8.67% + 14.26% + 26.92%). In the fast-growing young shoots of TH17/17 plants, a generally very high level of precursors for cytokinin phosphates was detected (MPH 186.09%), including that of the most physiologically active cytokinin trans-zeatin, the trans-zeatin riboside monophosphate: 500%. This result is consistent with the differential accumulation of trans-zeatin riboside in two maize hybrids that expressed heterosis in the development of the ear inflorescence (Shi et al., 2019).

Short-rotation coppice willow is a lignocellulose-rich energy crop whose green and woody biomass components may be considered as substrates for the anaerobic digestion (AD) process of biomethane production (Kakuk et al., 2021; Nyári et al., 2021). Our study demonstrates that breeding triploid hybrid willow genotypes can open a way forward for applying hybrid vigor to improve the bioenergy production system. As shown in Figure 4, positive or negative MPH in the methane yield (ml methane/1 g FW organic material) was obtainable when using the woody biomass from different THs. In the fermentation of woody tissues from the TH3/12 hybrids, the positive hybrid vigor (MPH:23.25%) in methane production was correlated with the lowest acid detergent lignin (ADL) content (13.89%). The greater efficiency in methane production was evident throughout the fermentation process (Figure 3). The highest ADL ratio (19.17%) was detected in wood samples from the Inger plants. Yet, this ADL level was reduced to 14.64% in the case of TH21/2 hybrids featuring negative heterosis in methane yield (MPH: –14.33%). Concerning their chemical composition, wood tissues from the TH17/17 hybrids contain more lignin than do tissues of its parents (MPH: 15.59%); hence, fermentation of this type of woods resulted in a reduced methane yield (MPH:- 27.80%). These trends are consistent with the previous findings that methane production decreases as the lignin content increases (Herrmann et al., 2016). In this respect, a negative correlation (r = −0.82) was obtained when using green or woody willow biomass for methane generation (Kakuk et al., 2021). Our experimental findings focus our attention upon fermentation efficiency as a special and important breeding target. Triploid heterosis offers a novel genetic basis for the selective improvement of valuable traits.