Skip to main content

Genotype-environment interaction and stability of fiber properties and growth traits in triploid hybrid clones of Populus tomentosa



Clones provide a sensitive method for evaluating genotypic stability and detecting genotype-environment (G × E) interactions because of non-additive genetic effects among clones and there being no genetic effect among ramets of an ortet. With this study, we aimed to confirm and expand earlier findings, estimate stability parameters, and provide accurate estimates of clonal repeatabilities and genetic gains for a triploid breeding program of P. tomentosa Carr.


Six 5-year-old clonal trials established in Northern China were used to determine the clonal variation, clone × site interactions, and the stability parameters of fiber properties of wood and growth traits. 360 trees from ten hybrid clones were collected from six sites. The clonal and site effects had a highly significant effect (P < 0.001) for all studied traits. While the clone × site interactions had a highly significant effect (P < 0.001) on fiber length (FL), coarseness (C), and tree growth (tree height [H], diameter at breast height [DBH] and stem volume [SV]), and a moderate effect (P < 0.05) on fiber width (FW) and fiber length/width (FL/W). For FL and SV, most of the triploid hybrid clones had higher reaction norms to the improvement in growth conditions and higher phenotypic plasticity. The estimated clonal repeatability of FW (0.93) was slightly higher than for FL (0.89), FL/W (0.83), C (0.91), DBH (0.76), H (0.85), and SV (0.80). Three clonal testing sites were sufficient to estimate quantitative parameters of fiber properties, however, more clonal testing sites would help improve the accuracy of quantitative parameters of the growth traits.


Our results highlight that accurate estimation of quantitative parameters for growth traits in triploid hybrid clones of P. tomentosa requires more clonal testing sites than the fiber properties.


Chinese white poplar, P. tomentosa Carr. (section Populus, family Salicaceae, genus Populus), is a widely distributed tree species in northern China and generally grows in mixed-species forests with Robinia pseudoacacia L. Since the first triploid hybrid of P. tomentosa was artificially created in 1995, the triploid hybrid clones have exhibited superior growth in China [1]. The latest research progress in wood fiber has led to the use of short fibers and short fibers mixed with coniferous long fibers for the production of high-quality paper [2]. The Chinese forestry industry has shown a new focus on using hybrid poplar for manufacturing evenly short fibers. Thus, an obvious increase in Chinese white poplar-use is expected, which in turn increases interest in breeding selected triploid hybrid clones [3,4,5].

The triploid hybrid clones of P. tomentosa have been shown to display a high degree of variability in growth, fiber traits, wood basic density, and wood chemical properties [1, 3,4,5,6]. This variability indicates that intriguing gains may be obtained through selection. The potential advantages of utilizing clonal methods in poplar hybrid clones have been reviewed by several authors [7, 8]. Genotype × environment (G × E) interaction normally refers to the performance differences of genotypes between environments [9]. Testing and selection in tree improvement programs are complicated by this interaction, and it often leads to reduced total genetic gains. In general, the literature on G × E is plentiful for forest trees [10, 11], however, the number of studied clones is relatively limited [12]. It is expected that clones are more sensitive in different environments than either seed- or family-origin plants because there is no genetic effects among ramets of an ortet and non-additive genetic effects among clones are substantial. Therefore, clones can provide a sensitive method for evaluating genotypic stability and detecting G × E interactions [13].

To obtain a better understanding of the genetic control of growth traits, fiber properties, wood basic density and wood chemical properties, basic genetics such as G × E interactions, inter-trait and inter-site correlations, clonal genetic variation, and repeatabilities for triploid hybrid clones of P. tomentosa were recently estimated [4, 5, 14, 15]. Both site effects and clonal effects were found to be significant for all tested traits and G × E interactions were also significant for all studied traits except for fiber length and holocellulose [5, 14]. This basic information is crucial to estimate expected genetic gains and develop suitable strategies for clonal forestry. However, most of the studies were conducted on only three or four sites. For the purpose of tree improvement, accurate estimates of quantitative parameters require that different genotypes are tested at more locations, therefore, in this present study we study six sites.

The objectives of the present study were: (1) to confirm and expand the earlier findings; (2) to estimate stability parameters of these triploid clones; and (3) to provide accurate estimated values of genetic gains and clonal repeatabilities for improving the breeding program of Chinese white poplar.


Basic statistics and variation within and between sites

The mean values, standard error (SE) of the mean, ranges of maximum and minimum values, and coefficients of phenotypic variation of all traits in the six clonal trials are presented in Table 1. A fairly obvious difference in growth was found between the six locations, with trees at the Taiyuan site exhibiting the fastest growth. For fiber properties, the longest fiber length (FL) and the greatest fiber length/width ratio (FL/W) were found in the trees at the Gaotang site. Trees from the Pinggu site had the shortest FL and the lowest FL/W. The difference between the longest and shortest means of FL was 7.9% (Table 1). However, the thinnest fibers were observed at Xiangfen, and the thickest, at Pinggu. The combined analysis of all six trials displayed that all studied traits exhibited significant site effects (Table 2). FL, fiber width (FW), and FL/W showed small phenotypic variation (CV = 3.2–5.5%) and the values were much lower than those of coarseness (C), diameter at breast height (DBH; at 1.3 m above ground level), tree height (H), and stem volume (SV).

Table 1 Clone mean values, standard error (SE) of the mean, ranges of maximum and minimum values, and coefficients of phenotypic variation (CV p %) for fiber properties and growth traits at the six clonal trials
Table 2 Analysis of variance results, estimates of repeatability (± SE) at the clonal mean and individual-tree levels of all traits of the nine triploid hybrid clones at the six sites combined

Clonal variation and repeatability

The results of the analysis of variance for the combined six sites are presented in Table 2. Significant clonal effects for all tested traits were seen (Table 2). For all studied traits, the variance caused by error (that is, the differences between ramets within a clone within a site) accounted for most of the variation (48.4 to 75.4% of the overall variation). However, most of the variance in FW (35.7%) was caused by the clone. Hence, the highest estimated repeatability of the clone mean (0.93) and the estimated individual-tree clonal repeatability (0.36) were observed for the FW.

The estimated clonal repeatability and estimated individual-tree clonal repeatability of all the studied traits at each site are shown in Table 3. No significant differences were found between the clones for H, DBH, and SV at the Yanzhou site, FL/W at the Gaotong site, as well as FW and FL/W at the Pinggu site. Therefore, we did not estimate the clonal repeatability for the six traits. The estimated clonal repeatability varied from 0.73 to 0.90 for FL, from 0.68 to 0.95 for FW, from 0.54 to 0.90 for FL/W, and from 0.78 to 0.93 for C. The estimated clonal repeatability of the growth traits ranged from 0.58 to 0.90 (Table 3).

Table 3 The results of analysis of variance and estimates of repeatability of all studied traits of the nine triploid hybrid clones at individual site

Clone × site interaction, stability, and reaction norms of triploid hybrid clones

In this study, all fiber properties and growth traits exhibited a significant clone × site interaction (Table 2). The level of significance for FW (0.031) and FL/W (0.014) was moderate and a highly significant effect (P < 0.001) was observed for FL, C, and tree growth. A higher estimated ratio in H (57.9%) for the clone × site interaction variance component versus the total clone×site plus clone variance component was observed compared with the other traits.

The inter-site (B-type) genotypic correlations are shown in Table 4. Some inter-site genotypic correlations between the same traits were not estimated because no significant differences were observed between the clones for H, DBH, and SV at the Yanzhou site, FL/W at the Gaotong site, and FW and FL/W at the Pinggu site. Most of the inter-site genotypic correlations between the same traits at different sites were either moderate or strong, except for H at Zhengzhou-Taiyuan and FL at Pinggu-Taiyuan (Table 4). For FL, FL/W, and the growth traits, the correlations were lower when the Taiyuan site was included.

Table 4 The estimates of type B genetic correlations among all possible pairs of sites for fiber properties and growth traits (with approximate SE in brackets) of the nine triploid hybrid clones

The stability parameters for FL and SV of the individual triploid hybrid clones across the six clonal trials, including the Finlay-Wilkinson parameters and the Shukla stability variances that characterize the reaction norms of the clones are shown in Table 5. The stability parameters of all of the studied clones were estimated only for FL and SV. For FL, the Shukla’s stability of clones B330 and B331 was slightly higher than those of the other clones due to their lower stability variance (1.01 × 10 − 4 and 2.28 × 10 − 4, respectively). However, for SV, lower stability variances were observed in clones B302, B303, and B306. Therefore, clones B302, B303, and B306 had higher stability for SV than the other clones. For FL, using the Finlay-Wilkinson’s method, the regression coefficient (bi) of clones B302, B303, B304, and B306 varied from 1.214 to 1.429, indicating that the four clones had higher reaction norms. However, the clones, B304, B330, and B331 had higher reaction norms for SV because their regression coefficients (bi) were 1.315, 1.411, and 1.477, respectively.

Table 5 Stability parameters for FL and SV of the nine triploid hybrid clones across six sites

Inter-trait genetic correlations, correlated genetic response, and genetic gain

The genotypic and phenotypic correlations between all of the studied traits at individual sites are presented in Table 6. At all sites, a positive estimated correlation was found between FL and growth traits. A significant positive estimated genetic correlation between FW and C was also observed at each site. However, a negative estimated genotypic correlation between growth traits (except for H) and the C existed. These results suggest that selection for growth traits (except for H) might lead to a slight decrease in C. SV, H, and DBH had some strong positive autocorrelations since SV was derived from H and DBH.

Table 6 Phenotypic correlations (below the diagonal) and genotypic correlations (above the diagonal) among all traits at individual sites of the nine triploid hybrid clones

The predicted genetic gains of the hybrid clones at the six sites for different selections are presented in Table 7. The results are based on the assumption of the estimated repeatabilities and correlations from direct clonal selection and the related genetic responses in all the studied traits. The predicted genetic gains were identical in C for Pinggu and Zhengzhou (6.1%), while the Taiyuan site had the highest gain (8.4%). However, the gains in C were about twice the gains observed for FL. The selection for C led to equal gains in FW as direct selection for FW. The selection for DBH led to the gain in SV equal to that of the direct selection for SV but resulted in a slight decrease in C (Table 7). Among growth traits, selection for C resulted in the lowest gains in DBH and SV, however, it had a greater negative impact on SV.

Table 7 Expected response (\( \varDelta G/\overline{X}\times 100 \)) in fiber properties and growth traits at the six sites for different selection criteria of the nine triploid hybrid clones


Variations between sites

Site effects represent the response of trees to the combined influences of edaphic and regional climatic conditions [16]. Even though the current trials were not designed to separate these different effects, some conclusions can still be drawn [16,17,18]. The poor field performance in growth at Yanzhou may be owing to poor drainage (due to soil compression) and wet weather at the early stages of the experiment. Significant site effects were observed in the combined analysis of all traits (Table 2). Moreover, a positive correlation between SV and latitude, and negative correlations between SV and longitude, as well as rainfall and mean annual temperature, were observed (Fig. 1.), suggesting that edaphic and regional climatic conditions had significant effects on the growth traits. At all sites, the fiber properties of the triploid hybrid clones had less phenotypic plasticity than the growth traits, as evidenced by the lower CV values, indicating that the site impacts fiber properties less than the growth traits.

Fig. 1
figure 1

Correlations between stem volume and (A) latitude, (B) longitude, (C) rainfall, and (D) mean annual temperature

Clonal variation and clone × site interaction

Former studies recorded that the effect of clone × site interactions on the growth traits at three clonal trials was greater than the clonal effect [4, 5]. Such a high interaction would lead to an increase in the complexity of testing [19, 20]. If the correlations between test and planting environments were alike and the objectives of an improvement program were restricted enough to allow selection for certain types of locations, then these interactions could be utilized. The impact of clonal effects on the growth traits in this study was, on the contrary, greater than the clone × site interactions when the number of testing sites was six. This suggests that an increase in the number of clonal testing sites helps improve the accuracy of clonal effects, thus reducing the complexity of clonal testing.

In this study, no significant positive type B genotypic correlations between genotypic values at the Xiangfen-Taiyuan sites for all studied traits were observed (Table 4), which suggests that a true G × E interaction exists. The differences between these sites may contribute to the clone × site interaction. However, the type B genotypic correlations for FL and C between the sites were mostly positive, which indicates that the fiber properties were more stable than the growth traits [21].

Stability and reaction norms of clones

For forest tree breeding, it is highly important to predict the yield of a genotype in different environments. This concept of stability may be determined by the stability variances [22] and the regression coefficients [23]. In this study, the results indicated that each of the studied clones had various magnitudes of phenotypic plasticity and specific patterns of FL and SV at the six sites. The stability values for the FL of the triploid hybrid clones were different for the different methods. This could be partially explained by the different principles that are used to estimate the stability parameters [24,25,26,27]. High values of the Finlay-Wilkinson stability parameter, i.e., the regression coefficient of the clone values on location value, determine the clones that benefit from a productive location, while high intercept coefficients determine the clones that can grow well under limited resource conditions. As shown by the higher regression coefficient means of clone B304, B330, and B331 for SV, these clones have higher norms of reaction to the improvement in growth environments and higher phenotypic plasticity.

Genotypic correlation between traits and clonal repeatability

Generally, an increase in tree growth was related to a slightly higher FL. Several previous studies have shown that FL increases with increasing tree growth [28,29,30,31]. In the present study, the genetic correlations between growth traits and FL were also positive at the six sites, indicating that triploid breeding of Populus not only improves the FL but also increases the tree growth.

Clonal repeatability of the fiber properties and growth traits across all six sites varied from medium to high (Table 2). For fiber properties, the estimated clonal repeatabilities ranged from 0.83 to 0.93, which was consistent with those that have been recorded in the literature [4, 5, 32]. This suggests that three clonal testing sites were enough to estimate quantitative parameters for fiber properties. However, clonal repeatabilities were estimated to be 0.76 for H, 0.85 for DBH, and 0.80 for SV, which were slightly higher than those estimated for the three clonal trials [5]. This implies that increasing the number of testing sites would help improve the accuracy of clonal repeatabilities of growth traits.

The analyses of the present study were based on nine genotypes and a relatively small number of testing sites measured over only five years. Therefore, the results from this study should be further verified by subsequent larger and longer trials. The implication of the G × E interaction for estimating the gain of genetically improved triploid hybrid clones of P. tomentosa should be further studied.


The experiments and analyses lead to the conclusion that the clone and site have a highly significant effect (P < 0.001) for all studied traits. The clone × site interactions had a moderate effect on FW and FL/W, and a highly significant effect (P < 0.001) on FL, C, and tree growth. The stability parameter analysis showed that most of the triploid clones had higher reaction norms to the improvement in growth conditions and higher phenotypic plasticity. The estimated clonal repeatability of FW (0.93) was slightly higher than for C (0.91), FL (0.89), FL/W (0.83), H (0.85), DBH (0.76), and SV (0.80). Three clonal testing sites were good enough to estimate quantitative parameters of fiber properties. However, more than three clonal testing sites would help improve the accuracy of quantitative parameters of growth traits. The findings in this study were due to a certain number of triploid hybrid clones derived from several parents (n = 5), therefore, the implication of G × E interaction for calculating the gain of genetically improved triploid hybrid clones of P. tomentosa should be further investigated.


Experimental design

In the present study, the materials used were sampled from six triploid hybrid poplar clonal trials in northern China. These trials were established by Beijing Forestry University using cuttings of P. tomentosa on sandy loam with typical soil fertility. The clonal trials at Yanzhou, Gaotang, and Xiangfen were built in the spring of 2004 and the clonal trials at Pinggu, Taiyuan, and Zhengzhou were established in the spring of 2005. The characteristics of the six clonal experiments are shown in Table 8. The planted cuttings were 1.5–2.0 cm in diameter and 20–25 cm in length. The experiments were repeated three times in the nursery. The experiments included one diploid clone and nine triploid hybrid clones at each site (Table 9). The nine triploid hybrid clones were preselected among the triploids derived from five parent trees because they exhibited good growth rates [33]. The diploid clone (M1319) was an elite tree of P. tomentosa that grows faster. A randomized complete block design with three replicates was used in each clonal trial (240 trees per plot, resulting in a total of 720 trees per clone per trial). Ten clones were planted on rectangular plots at each location, and each plot contained 240 trees (4 × 60 trees) with 2 m × 3 m spacing. No thinning was conducted during the testing period.

Table 8 Location, climatic conditions, and description of six clonal trials
Table 9 Identity and origin of the hybrid clones

Measurement of fiber properties and growth traits

A total of 360 trees were randomly selected from the six trial sites. Trees were collected randomly, 6 ramets (the 6 trees) of each clone were sampled fully at random per location. In the triploid poplar breeding program of Beijing Forestry University, the overall H and DBH of all tree samples were measured by workers at the age of 5 years. The SV of each sample was estimated according to the volume function of DBH and H used by Chen [34]. All sampled trees were shipped to Beijing Forestry University. For laboratory measurements, a 10-cm-thick stem disk was harvested from the breast height of each sample. Matchstick-sized wood samples (each sample implying two pairs of annual rings) were cut from the stem disks to evaluate the fiber properties, and then immersed in a boiling 1:1 (v/v) mixed solution of hydrogen peroxide and acetic acid. Afterward, FL, FW, and C were measured according to the methods described by Zhang et al. [4] using the Lorentzen & Wettre (LW) Fiber Tester (AB Lorentzen & Wettre, Kista, Sweden).

Statistical analysis

When the main statistical analysis and the estimation of genetic statistics were conducted, the diploid clone (M1319) was not included. Thus, only the nine triploid clones were used in the main variance analysis and estimates of genetic parameters. Analyses of variance were performed according to the UNIVARIATE program of the SPSS software (SPSS for Windows, version 13, SPSS, Chicago, IL). Analysis of variance was used to analyze the variation between the ramets of the sample clones at a location based on the following linear model (Eq. (1)):

$$ {X}_{ik}=\mu +{C}_i+{\varepsilon}_{ik} $$

where Xik represents the observation of the kth ramet of the ith clone; μ represents the general mean; Ci represents the effect attributed to the ith clone and εik represents random error. The clone × replicate effects were not included in the linear model (Eq. (1)) as the six sampled trees per clone per location were randomly selected.

The clonal repeatabilities were calculated within a location as:

$$ {R}_c^2=\frac{{\hat{\sigma}}_c^2}{{\hat{\sigma}}_c^2+\frac{{\hat{\sigma}}_e^2}{k}} $$

and the clonal repeatability of the individual-tree (\( {R}_b^2 \)) was estimated according to the following equation:

$$ {R}_b^2={\hat{\sigma}}_c^2/\left({\hat{\sigma}}_c^2+{\hat{\sigma}}_e^2\right) $$

where k is the average value of tree samples per clone within a location; \( {\hat{\sigma}}_c^2 \) is the estimated variance of clone; \( {\hat{\sigma}}_e^2 \) is the variance among the ramets within the clones.

The six locations were jointly analyzed using the following linear model [17] (Eq. (4)):

$$ {X}_{ijk}=\mu +{C}_i+{L}_j+{C}_i{L}_j+{\varepsilon}_{ijk} $$

where Xijk represents the performance of the kth ramet of the ith clone within the jth location; μ represents the general mean; Ci represents the effect attributed to the ith clone; Lj represents the effect attributed to the jth location; CiLj represents the interaction between the ith clone and jth location, and εijk represents random error. Only the location was regarded as a fixed effect and all other terms were regarded as random effects.

The clonal repeatabilities were estimated using the equation:

$$ {R}_c^2=\frac{{\hat{\sigma}}_c^2}{\frac{k_2{\hat{\sigma}}_c^2}{k_2}+\frac{k_1{\hat{\sigma}}_{Lxc}^2}{k_2}+\frac{{\hat{\sigma}}_e^2}{k_2}} $$

and the clonal repeatability of individual-trees was calculated across locations according to the following equation:

$$ {R}_b^2={\hat{\sigma}}_c^2/\left({\hat{\sigma}}_c^2+{\hat{\sigma}}_{Lxc}^2+{\hat{\sigma}}_e^2\right) $$

where k1 is the coefficient related to the variance caused by the clone × site interaction item (\( {\hat{\sigma}}_{Lxc}^2 \)); k2 is the coefficient related to the variance caused by the clonal variation (\( {\hat{\sigma}}_c^2 \)). The standard errors (SE) for the estimates of clonal repeatability were estimated according to the following equation [35]:

$$ SE\left({R}_c^2\right)=\sqrt{\frac{2{\left(1-{R}_c^2\right)}^2{\left[1+\left({k}_2-1\right){R}_c^2\right]}^2}{k_2\left({k}_2-1\right)N-1}} $$

where N is the number of tested clones. The SE for the clonal repeatability estimates of individual-tree were also estimated according to eq. (7), where \( {R}_b^2 \) was used instead of \( {R}_c^2 \).

The genetic correlation coefficients (rA(XY)) between traits at a location were computed using the equation [36]:

$$ {r}_{A\left(X,Y\right)}=\frac{{\hat{\sigma}}_{c\left(x,y\right)}}{\sqrt{{\hat{\sigma}}_{c(x)}^2{\hat{\sigma}}_{c(y)}^2}} $$

where \( {\hat{\sigma}}_{c\left(x,y\right)} \) is the clonal component of covariance estimated among traits x and y; \( {\hat{\sigma}}_{c(x)}^2 \) is the clonal variance component estimated for trait x; \( {\hat{\sigma}}_{c(y)}^2 \) is the clonal variance component estimated for trait y. Similar to using the mean squares for estimating variance components, the cross-products were used to estimate the corresponding clone component of variance from the data collected from the same individual tree.

The approximate SE of genetic correlation estimates were estimated based on the following formula [9]:

$$ SE=\frac{1-{r}^2}{\sqrt{2}}\sqrt{\frac{{\hat{\sigma}}_{\left({R}_x^2\right)}{\hat{\sigma}}_{\left({R}_y^2\right)}}{R_x^2{R}_y^2}} $$

where r is the estimate of genetic correlation; \( {R}_x^2 \) is the estimate of the clonal repeatability of the trait x; \( {R}_y^2 \) is the estimate of the clonal repeatability of the trait y; \( {\hat{\sigma}}_{\left({R}_x^2\right)} \) and \( {\hat{\sigma}}_{\left({R}_y^2\right)} \) are the standard error for \( {R}_x^2 \) and \( {R}_y^2 \), respectively. Using Eq. (7) as the basis for the t-tests, the significance of the genetic correlations between the traits was tested.

The genetic correlation estimates among the same traits for pairs of locations were estimated to assess the contribution of each pair of treatments to the overall clone × environment. These type B genotypic correlations were computed according to the measurements of sample ramets at different locations from the same clones planted using the equation [37]:

$$ {r}_{B\left(X,Y\right)}=\frac{r_{P\left(x1,y2\right)}}{R_{C(x1)}{R}_{C(y2)}} $$

where rp(x1,y2) is the coefficient of phenotypic correlation between the clonal means estimated between x measured at location 1 and y measured at location 2; Rc(x1) and Rc(y2) are the square roots of the mean repeatability of the clones x and y estimated at location 1 and 2, respectively.

The SPSS PROC CORR software was used to calculate the phenotypic correlations of the clonal means for each pair of sites, and the significance test was performed using the F-test.

In the present study, it is assumed that the selection intensity is equal to 1.271, which corresponds to selecting four clones from 18 or three clones from 14 [16]. According to the relevant response of the clonal selection in trait x, the expected genetic gain in trait y ΔGy can be predicted using the formula [9]:

$$ \varDelta {G}_y={i}_x\sqrt{R_x^2}{\sigma}_y{r}_{xy} $$

where ix (ix = 1.271) is the selection intensity; \( {R}_x^2 \) is the clonal repeatabilities for trait x; σy is the standard deviation of clonal for trait y; rxy is the genetic correlation coefficients between trait x and y (rxy ≤ 1).

Estimation of stability parameters

To evaluate the stability of the individual clones across the six experimental locations, the stability variances were estimated using the Shukla model [22], where lower values indicate higher stability, and the significance test was conducted using the F-test. Subsequently, the stability of the individual clones across locations was also estimated using the Finlay and Wilkinson stability parameter [23]. The triploid clones were grouped according to the size of their regression coefficients (bi), less than, equal to, or > 1. Clones with bi > 1 are more adapted to favorable growing conditions, those with bi < 1 are adapted to unfavorable environmental conditions, and those with bi equal to 1 have an average adaptation to all environments [23].

Availability of data and materials

All data generated or analysed during this study are included in this published article.


CV p %:

coefficients of phenotypic variation


fiber length


fiber width


fiber length/width




tree height


diameter at breast height


stem volume

\( {R}_b^2 \) :

repeatability of individual-tree clonal

\( {R}_c^2 \) :

repeatability of clonal means


standard errors













F :

Fisher’s criteria

P :

probability levels

t :

Student’s criteria

b i :

regression coefficient


  1. Zhu ZT, Lin HB, Kang XY. Studies on allotriploid breeding of Populus tomentosa B301 clones. Sci Silvae Sin. 1995;31:499–505.

    Google Scholar 

  2. Ranua J. Why is the Finnish forest industry interested in aspen? Sorbifolia. 1996;4:178–80.

    Google Scholar 

  3. Chen TY, Wang B, Wu YY, Wen JL, Liu CF, Yuan TQ, et al. Structural variations of lignin macromolecule from different growth years of triploid of Populus tomentosa Carr. Int J Biol Macromol. 2017;101:747–57.

    CAS  Article  PubMed  Google Scholar 

  4. Zhang PD, Wu F, Kang XY. Genotypic variation in wood properties and growth traits of triploid hybrid clones of Populus tomentosa at three clonal trials. Tree Genet Genom. 2012;8(5):1041–50.

    Article  Google Scholar 

  5. Zhang PD, Wu F, Kang XY. Genetic control of fiber properties and growth in triploid hybrid clones of Populus tomentosa. Scand J For Res. 2013;28(7):621–30.

    Article  Google Scholar 

  6. Wu F, Zhang PD, Pei JC, Kang XY. Within-stem variation in wood properties of triploid hybrids of Populus tomentosa at three clonal trials. Scand J For Res. 2014;29(1):71–6.

    CAS  Article  Google Scholar 

  7. Michler CH, Bauer EO. High frequency somatic embryogenesis from leaf tissue of Populus spp. Plant Sci 1991;77:111–118. (91)90186-C, 1.

  8. Yu QB. Can physiological and anatomical characters be used for selecting high yielding hybrid aspen clone? Silva Fennica 2001;35:137–146. 10.14214/sf.591, 2.

  9. Falconer DS. Introduction to quantitative genetics. 3rd ed. London: Longman; 1981. p. 285–6.

    Google Scholar 

  10. Wu HB, Lei J, Li XY, Wang H, Duan AG, Zhang JG. Aggregation distributions across stand age in provenances of Cunninghamia lanceolata (lamb.) hook. Forest Ecol Manag. 2021;494:119317.

    Article  Google Scholar 

  11. Yuan CZ, Zhang Z, Jin GQ, Zheng Y, Zhou ZC, Sun LS, et al. Genetic parameters and genotype by environment interactions influencing growth and productivity in Masson pine in east and Central China. Forest Ecol Manag. 2021;487:118991.

    Article  Google Scholar 

  12. Hayatgheibi H, Haapanen M, Lundströmer J, Berlin M, Kärkkäinen K, Helmersson A. The impact of drought stress on the height growth of young Norway spruce full-sib and half-sib clonal trials in Sweden and Finland. Forests. 2021;12(4):498.

    Article  Google Scholar 

  13. Bentzer BG, Foster GS, Hellberg AR, Podzorski AG. Genotype × environment interaction in Norway spruce involving three levels of genetic control: seed source, clone mixture, and clone. Can J For Res. 1988;18(9):1172–81.

    Article  Google Scholar 

  14. Zhang PD, Wu F, Kang XY. Chemical properties of wood are under stronger genetic control than growth traits in Populus tomentosa Carr. Ann For Sci. 2015;72(1):89–97.

    CAS  Article  Google Scholar 

  15. Zhang PD, Wu F, Kang XY, Zhao CG, Li YJ. Genotypic variations of biomass feedstock properties for energy in triploid hybrid clones of Populus tomentosa. Bioenerg Res. 2015;8(4):1705–13.

    CAS  Article  Google Scholar 

  16. Pliura A, Zhang SY, Mackay J, Bousquet J. Genotypic variation in wood density and growth traits of poplar hybrids at four clonal trials. For Ecol Manag. 2007;238(1-3):92–106.

    Article  Google Scholar 

  17. Blazier MA, Hennessey T, Schimleck L, Abbey S, Holbrook R, Dahlen J. Long-term effects of stand density management and genotype on wood properties of loblolly pine (Pinus taeda L.) in the mid-South USA. Forest Ecol Manag. 2021;491:119176.

    Article  Google Scholar 

  18. Szala K, Ogonowska H, Lugowska B, Zmijewska B, Wyszynska R, Dmochowska BM, et al. Different sets of TaCKX genes affect yield-related traits in wheat plants grown in a controlled environment and in field conditions. BMC Plant Biol. 2020;20(1):496.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. Randall WK, Cooper DT. Predicted genotypic gains from cottonwood clonal tests. Silvae Genet. 1973;22:165–7.

    Google Scholar 

  20. Lindgren D. Prediction and optimization of genetic gain with regard to genotype × environment interactions. Studia Forestalia Suecica. 1984;166:15–23.

    Google Scholar 

  21. Yu Q, Pulkkinen P, Rautio M, Haapanen M, Alen R, Stener LG, et al. Genetic control of wood physiochemical properties, growth and phenology in hybrid aspen clones. Can J For Res. 2001;31(8):1348–56.

    Article  Google Scholar 

  22. Shukla GK. Some statistical aspects of partitioning genotype-environment components of variability. Hered. 1972;29(2):237–45.

    CAS  Article  Google Scholar 

  23. Finlay KW, Wilkinson GN. The analysis of adaptation in a plant breeding program. Aus J Agric Res. 1963;14(6):742–54.

    Article  Google Scholar 

  24. Swanckaert J, Akansake D, Adofo K, Acheremu K, De BB, Eyzaguirre R, et al. Variance component estimations and mega-environments for sweetpotato breeding in West Africa. Crop Sci. 2020;60(1):50–61.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Macholdt J, Styczen ME, Macdonald A, Piepho HP, Honermeier B. Long-term analysis from a cropping system perspective: yield stability, environmental adaptability, and production risk of winter barley. Eur J Agron. 2020;117:126056.

    CAS  Article  Google Scholar 

  26. Gouvêa LRL, Silva GAP, Verardi CK, Silva JQ, Scaloppi-Junior EJ, Gonçalves PDS. Temporal stability of vigor in rubber tree genotypes in the pre- and post-tapping phases using different methods. Euphytica. 2012;186(3):625–34.

    Article  Google Scholar 

  27. Mohammadi RE, Amri A. Comparison of parametric and non-parametric methods for selecting stable and adapted durum wheat genotypes in variable environments. Euphytica. 2008;159(3):419–32.

    Article  Google Scholar 

  28. Nazari N, Bahmani M, Kahyani S, Humar M, Koch G. Geographic variations of the wood density and fiber dimensions of the Persian oak wood. Forests. 2020;11(9):1003.

    Article  Google Scholar 

  29. Ahmed AKM, Fu ZX, Ding CJ, Jiang LP, Han XD, Yang AG, et al. Growth and wood properties of a 38-year-old Populus simonii × P. nigra plantation established with different densities in semi-arid areas of northeastern China. J Forestry Res. 2020;31(2):497–506.

    Article  Google Scholar 

  30. Yanchuk AD, Dancik BP, Micko MM. Variation and heritability of wood density and fiber length of trembling aspen in Alberta, Canada. Silvae Genet. 1984;33:11–6.

    Google Scholar 

  31. Xie HF, Yu ZK, Chen DL, Guo DC, Song FX. Selection and genetic analysis for wood properties in pulpwood varieties of black poplar. J Shandong For Sci. 1996;1:1–5.

    Google Scholar 

  32. Wu F, Zhang PD, Pei JC, Kang XY. Genotypic parameters of wood density and fiber traits in triploid hybrid clones of Populus tomentosa at five clonal trials. Ann For Sci. 2013;70(7):751–9.

    Article  Google Scholar 

  33. Zhu ZT. Genetic improvement of Populus tomentosa. Beijing: Chinese Forestry Press; 2006. p. 155–220.

    Google Scholar 

  34. Chen ZS. The establishment of binomial stock volume tables of poplar. J For Res. 1989;2:78–83.

    Google Scholar 

  35. Becker WA. Manual of quantitative genetics. Washington: Academic Enterprises, Pullman; 1992. p. 160.

  36. Becker WA. Manual of quantitative genetics. Washington:Academic Enterprises, Pullman; 1984. p. 190.

  37. Burdon RD. Genetic correlation as a concept for studying genotype × environment interaction in forest tree breeding. Silvae Genet. 1977;26:168–75.

    Google Scholar 

Download references


The authors would like to thank Dr. Jicheng Pei from Tianjin University of Science and Technology for his assistance with fiber properties measurements. We also thank Rebecca James in University of Otago for her language editing.


This study was mainly supported by the Fundamental Research Funds for Central Universities (2021JD01) and the National Key Research and Development Program of China (2016YFD0600403).

Author information




JW was responsible for generating the phenotypic data set, analyzing the results, and writing the paper. QZ and YS was responsible for tree growth measurements. XK was responsible for planning and guiding this research. PZ was responsible for planing and building the clonal trails, collecting wood discs sample, fiber properties measurements and revised the paper. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Pingdong Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

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

Wu, J., Zhou, Q., Sang, Y. et al. Genotype-environment interaction and stability of fiber properties and growth traits in triploid hybrid clones of Populus tomentosa. BMC Plant Biol 21, 405 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Genotype-environment interactions
  • Stability parameter
  • Clonal repeatability
  • fiber property
  • Growth trait
  • P. tomentosa