- Research Paper
- Open access
- Published:

# Advantage of clonal deployment in Norway spruce (*Picea abies* (L.) H. Karst)

*Annals of Forest Science*
**volume 77**, Article number: 14 (2020)

## Abstract

###
*Key message*

*Key message*

**There is considerable genetic gain of tree volume from clonal deployment in Norway spruce (**
*Picea abies***(L.) H. Karst) and clonal deployment will have at least 50% more or double genetic gain than the seedling deployment.**

###
*Context*

*Context*

Genetic parameters and genetic gains for wood quality and growth traits were estimated in six large clonal progeny trials.

###
*Aims*

*Aims*

Develop the optimal clonal deployment strategy of Norway spruce in Sweden.

###
*Methods*

*Methods*

Wood quality and growth traits were measured in all clonal trials and additive and non-additive genetic variances are partitioned.

###
*Results*

*Results*

Additive and non-additive genetic variances were equally important for growth traits while non-additive variance was small or not significant for wood quality trait. The genetic gain predicted for clonal deployment was greater than any of the other four deployment strategies. Selecting the top 1% of tested clones (clonal forestry) would have 48.4% and 134.6% more gain than the gain predicted for the seedling deployment of selected full-sib families and half-sib family (family forestry), respectively, at the same selection intensity.

###
*Conclusion*

*Conclusion*

This study highlights that testing of 30–40 clones per family would maximize the realized genetic gain for different clonal selection scenarios, either selecting the best ten or 20 clones without any co-ancestry restrictions or selecting the best single clone from each of the best ten or 20 families (e.g., co-ancestry restriction). Clonal mean selection and vegetative deployment are the most effective.

## 1 Introduction

The total genetic variance in a biological trait can be partitioned into additive, dominance, and epistatic variances (Falconer and Mackay 1996). In forestry, however, quantitative geneticists have commonly ignored the interaction between alleles (dominance) at a locus and the interactions alleles between loci (epistasis) in the analysis of progeny trials (Foster and Shaw 1988). The main reasons are the lack of well-designed experiments with a full-sib family structure in a clonally replicated trial and limited number of families and clones within a family tested. Traditionally, two methods have been used to partition the genetic variance in forestry without using inbreeding methods. First, Foster and Shaw (1988) estimated the epistatic genetic variance from a clonally replicated experiment using the expected covariance among relatives. The model assumes that the epistatic genetic effect is mainly from numerous high-order interactions between loci. In contrast, another model proposed by Wu (1996) was based on the assumption that epistasis for a quantitative trait is limited to interactions between pairs of quantitative trait loci (QTL) only. Currently, many tree species based on genomic data has discovered significant epistatic genetic variance based on the interaction between QTL pairs, such as additive × additive, additive × dominance, and dominance × dominance interaction effects (Muñoz et al. 2014; Bouvet et al. 2016; Gamal El-Dien et al. 2016; Chen et al. 2019; Tan et al. 2018). Without clonally replicated experiments, however, it is difficult to explain how total genetic variance is distributed among the three components and make a sound judgment in optimal deployment strategy.

Mullin and Park (1992) defined five types of commonly used heritabilities applicable to full-sib clonal trials and its relevance for selection and deployment strategies in forestry. They are as follows: (1) narrow-sense heritability (*h*^{2}); (2) broad-sense heritability (*H*^{2}); (3) half-sib family mean heritability (\( {H}_{\overline{\mathrm{HS}}}^2 \)); (4) full-sib family mean heritability (\( {H}_{\overline{\mathrm{FS}}}^2 \)); and (5) clone mean heritability (\( {H}_{\overline{\mathrm{CL}}}^2 \)). For each of the five types of heritability, a corresponding selection and deployment strategy was proposed as: (1) forward mass selection of the best individuals and selected trees grafted into seed orchards for deployment of their seedlings (i.e., Seedling-based Open-Pollinated individual Forward selection (SOPF)); 2) forward mass selection of the best individuals and deployment with vegetative propagation (i.e., Clonal-based INdividual Forward selection (CINF)) using rooted cuttings (RC) or somatic embryogenesis (SE); (3) selecting the best parents based on their breeding values (backward general combining ability (GCA) selection) followed by repeated crossings in nature or in a grafted seed orchard in order to produce and deploy the progeny seedlings (i.e., Seedling-based Half-Sib family Backward selection (SHSB)); (4) selecting the best full-sib family (backward specific combining ability (SCA) selection) followed by repeated controlled crosses in a grafted seed orchard in order to produce and deploy progeny seedlings (i.e., Seedling-based Full-Sib family Backward selection (SFSB)); and 5) selecting the best clones, tested in a trial with replication, followed by deployment of these clones using vegetative propagation methods such as RC or SE (i.e., Clonal-based Replicated Clone Forward selection (CRCF)). In summary, *h*^{2}, *H*^{2}, and \( {H}_{\overline{\mathrm{CL}}}^2 \) were used to forward mass selection, while \( {H}_{\overline{\mathrm{HS}}}^2 \) and \( {H}_{\overline{\mathrm{FS}}}^2 \) are used to backward mass selection. The three selection and deployment strategies including SOPF, SHSB, and SFSB which are related to *h*^{2}, \( {H}_{\overline{\mathrm{HS}}}^2 \), and \( {H}_{\overline{\mathrm{FS}}}^2 \), respectively, are often referred to as *family forestry* strategies and are based on all seedlings produced or rooted cuttings from these seedlings in some conifer species such as radiata pine in Australia and Norway spruce in Sweden. CINF and CRCF are, on the other hand, referred to as clonal forestry because the best clones were selected for vegetative propagation and deployment. The only difference between CRCF and CINF is that the latter does not require the establishment of a clonally replicated test but rather deploy the best individuals using vegetative propagation.

Norway spruce (*Picea abies* (L.) H. Karst) is the most commercially important tree species in Swedish forestry (Hannrup et al. 2004). The Swedish breeding program for Norway spruce has been focused on growth, survival, branch quality, and wood density (Rosvall et al. 2011). Genetic parameters have been estimated for many full-sib and half-sib seedling progeny trials (Hannrup et al. 2004; Kroon et al. 2011; Chen et al. 2014) and also for clonal trials, but without family structure (Bentzer et al. 1989; Karlsson et al. 2001; Högberg and Dutkowski 2010; Isik et al. 2010). A most recent study of Norway spruce in clonal trials included family structure (Berlin et al. 2019), but considered just growth properties. There was a particular call to accurately estimate non-additive genetic variance for Norway spruce in a recent review of clonal forestry of Norway spruce (Wu 2018). In estimating genetic parameters using clone trials in conifer, another issue is that the number of families available has usually been in the low double-digit figures (Baltunis et al. 2007; Weng et al. 2008) or even less than ten (Isik et al. 2003; Isik and Kleinschmit 2005). The limited population size makes it difficult to accurately estimate the different sources of the genetic variance of growth traits. For the purpose of clonal forestry deployment, selection mainly focuses on the genetic variation within just a few elite families. Thus, it is important to capture enough of the variation (number clones per family) within these elite families. Therefore, one issue is to find the optimal number of clones per family in order to capture a large part of the within-family variation, achieve largest genetic gains, and maintain genetic diversity under constraints of the trial size and costs. This issue has not yet been sufficiently resolved for Norway spruce.

The aims of this study were to (1) partition the total genetic variance into additive and different non-additive sources of variation for growth and wood quality traits assessed at different ages in clonally replicated field trials of large family size; (2) estimate genetic parameters for growth and wood quality traits and compare them with previous studies in Norway spruce; and (3) evaluate the effects of family size and diversity on the genetic gain in a population. This study was done in four control-pollinated and two open-pollinated clonal trials of Norway spruce located in southern and central Sweden.

## 2 Materials and methods

### 2.1 Plant material and field trial

#### 2.1.1 Full-sib clonal trials

The seed underlying the clonal material of this study originated from a sparse partial diallel mating design between selected and trees within southern Sweden’s seed orchard zone 8-9S. Forty-nine parents including 30 females and 22 males were used for making the crosses. In total, 32 full-sib crosses have been done. The produced seeds were sown in spring 2001 into pots. After 1-year growth, the seedling (donor) plants were transplanted to nursery beds and arranged in a row-wise manner with an average of approximately 60 donor plants per family (range from 11 to 109 surviving plants per family). In the early spring of 2004, clones for the propagation of cuttings were visually selected within each family based on good growth and vitality. Roughly 70% of the donor plants per family were selected as candidates for field testing (varying between 38 and 97% of clones within each family).

Finally, rooted cuttings from 1430 clones from 32 full-sib families were planted in four clonal trials in 2007, which were registered as S21S0721387 (S1387), S21S0721388 (S1388), S21S0721389 (S1389), and S21S0721390 (S1390). A randomized incomplete block design using the single-tree plot was used in all four trials. Finally, each family from four trials had an average of 44.7 clones (18–64 clones). Detailed descriptions about location, soil type, and climate condition of trials are given in Table 1.

#### 2.1.2 Half-sib clonal trials

Two clonal half-sib progeny trials of Norway spruce, S22S94202237 (S237) at Lugnet and S22S94202239 (S239) at Grangärde, were established in 1994 in central Sweden. The detailed description of the two trials was shown in Table 1. Rooted cuttings for 2244 clones from 138 half-sib families were planted with an average of 16.9 clones for each family. A randomized incomplete block design using a single-tree plot was used in both trials.

### 2.2 Measurement of Pilodyn penetration depth in standing trees

In order to indirectly assess wood density under the bark, the Pilodyn penetration depth (Pilo) was measured using a Pilodyn 6J Forest (PROCEQ, Zurich, Switzerland) with a 2.0-mm diameter pin, without removing the bark. The measurement was conducted at approximately 1.3 m above the ground for each tree and on the same side for all trees in the trial S239 (Table 2). In the trial S1387, Pilo was measured from two sides of the tree and the average was used for the analysis.

### 2.3 Measurement of the resistograph

The micro-drill Resistograph IML-RESI PD300 (Instrumenta Mechanic Labor, Germany) was used, as an alternative method in order to determine the wood density of standing trees in the trial S237. The resistograph was used for drilling trees bark to pith to bark at a tree height of ca. 1.3 m and special attention was paid to avoiding drilling through knots or visible stem damages. Each profile was checked immediately after drilling on the tool’s screen and the measurement was repeated when necessary. Custom software is available from a web URL (https://forestquality.shinyapps.io/EucalyptResiProcessor/) (Downes et al. 2018). Finally, the mean resistance (resistograph (Resi)) values were provided from the data profile.

### 2.4 Measurement of acoustic velocity

The Hitman ST300 tool (Fiber-gen, Christchurch; New Zealand) was used to determine the acoustic velocity (AV) at the tree side with fewer branches in two half-sib trials (Table 2). All obstructing branches below the height of approximate 2 m were removed from the trees before the measurement. To avoid knots as much as possible, the upper probe was usually inserted just below the higher branch whorl and the lower probe was placed just above the lower whorl.

### 2.5 Estimating the dynamic modulus of elasticity

Based on the observed high correlations between Pilo and wood density, between AV and microfibril angle (MFA) (Chen et al. 2015), and between resistograph (Resi) and wood density (Fundova et al. 2018), either Pilo or Resi was used as a surrogate for estimation of wood density of a standing tree. The dynamic modulus of elasticity (MOE) was thus calculated using the following two models, respectively:

### 2.6 Measurement of growth traits

The tree height (cm) was measured by the tangent (clinometer) method and the diameter at breast height, ca. 1.3 m (DBH, mm), was measured by digital calipers.

### 2.7 Estimating the tree volume

Volume (VOL) of each tree (dm^{3}) was estimated with a function by Brandel (1990):

where DBH is the diameter at breast height (1.3 m) (cm) and HT is the tree height (m).

### 2.8 Statistical analysis

#### 2.8.1 Full-sib trials

Across-site analyses for traits measured in the four southern Swedish full-sib clonal trials were conducted with the following single-tree model

where *y* is the vector of the observations for a specific trait in the four full-sib clonal trials, *b* is the vector of fixed effects (i.e., overall mean, trial, and complete blocks within a trial), *a*, *f*, *c*, al, fl, and cl are the vectors of random additive genetic effects of individual genotypes, random special combining ability (SCA) effects, clones within full-sib family effects, additive genetic effect by location (trial) interaction effects, SCA-by-location effects, and clones within full-sib family-by-location effects, respectively. ** X**,

*Z*_{1},

*Z*_{2}

*Z*_{3},

*Z*_{4},

*Z*_{5}, and

*Z*_{6}are the known incidence matrices related to the

*b*,

*a*,

*f*,

*c*, al, fl, and cl, respectively. The variances associated with the random effects:

*a*,

*f*,

*c*, al, fl, and cl, are referred to as A\( {\hat{\sigma}}_A^2 \), I\( {\hat{\sigma}}_{\mathrm{SCA}}^2 \), I\( {\hat{\sigma}}_C^2 \), A\( {\hat{\sigma}}_{AxL}^2 \), I\( {\hat{\sigma}}_{\mathrm{SCA} xL}^2 \), and I\( {\hat{\sigma}}_{CxL}^2 \), respectively, where

*A*is the pedigree relationship and

*I*is the incidence matrix. The vector

*e*represents residual effects ~

*N*(0,

*R*), with heterogeneity of variance among

*j*sites \( R={\sigma}_{e_1}^2{I}_{n_1}\bigoplus {\sigma}_{e_2}^2{I}_{n_2}\bigoplus \dots \bigoplus {\sigma}_{e_j}^2{I}_{n_j} \), where ⨁ is the direct sum The average of the four residual variances from the four full-sib trials \( {\sigma}_e^2 \) was also estimated. Traits of HT, VOL, and Pilo at age 12 were only measured in the S1387 full-sib trial (Table 2); thus, the following single-site and single-tree model was used:

where all effects are the same as Eq. (1), except that *b* only includes the overall mean and block effects.

#### 2.8.2 Half-sib trials

The across-site analysis for all traits in the two half-sib trials used the following single-tree model

where *y* is the vector of the observations for trait in two half-sib clonal trials, *b* is the vector of fixed effects (i.e., overall mean and trial), and *a* and *al* are defined as above, except for *b*_{0}, *n*, and nl, which are the vectors of incomplete block within a trial, clones within half-sib family (i.e., non-additive effects), and clones within half-sib family-by-location effects, respectively. ** X**,

**Z**

_{0},

*Z*_{1},

*Z*_{3},

*Z*_{4}, and

*Z*_{6}are the known incidence matrices related to the

*b*,

*b*

_{0},

*a*,

*n*, al, and nl, respectively. The variances associated with the random effects:

*b*

_{0,}

*a*,

*n*, al, and nl, were referred to \( {\hat{\sigma}}_b^2 \), \( {\hat{\sigma}}_A^2 \), \( {\hat{\sigma}}_N^2 \), \( {\hat{\sigma}}_{AxL}^2 \), and \( {\hat{\sigma}}_{NxL}^2 \), respectively. In the half-sib trials, to avoid confounding of incomplete block and genetic effects, the block effects were here assumed to be random. The vector

*e*represents residual effects ~

*N*(0,

*R*), where heterogeneity of variance between the two sites is modeled as \( R={\sigma}_{e_1}^2{I}_{n_1}\bigoplus {\sigma}_{e_2}^2{I}_{n_2} \), and the average of the two residual variances from the two half-sib trials \( {\sigma}_e^2 \) was also estimated. All traits measured in only a single half-sib family clonal trial (see Table 2) were analyzed with the single-site and single-tree model as follows

where all the definitions of the respective effects are the same as in Eq. (3).

The variance associated with

is the estimate of additive genetic variance in both full-sib and half-sib trials. *V*_{A}, *V*_{AA}, and *V*_{AAA} pertain to the additive genetic variance and epistatic genetic variance due to interactions of additive effects at two and three loci, respectively.

is the estimate of dominance genetic variance in full-sib clonal trials. *V*_{D}, *V*_{AA}, *V*_{AD}, and *V*_{DD} correspond to dominance genetic variance and epistatic genetic variance due to interactions of additive-by-additive, additive-by-dominance, and dominance-by dominance effects at two loci.

is the estimate of epistatic genetic variance in full-sib clonal trials and

is the non-additive genetic (dominance and epistatic) variance in the half-sib clonal trials. To estimate the total genetic variance (\( {\hat{V}}_{\mathrm{G}} \)), the equations:

and

were used to full-sib and half-sib clonal trials, respectively. Furthermore, the total phenotypic variance was estimated as

and

for full-sib and half-sib clonal trial analyses, respectively. Also, we adjusted the equation (Weng et al. 2008) to calculate the phenotypic variances of half-sib family means as

and

for full-sib and half-sib clonal trials, respectively. The coefficients of denominator degrees of freedom are *k*_{1} – *k*_{6} and *k*_{7} – *k*_{10} for the model terms of full-sib and half-sib clonal trials, respectively, and were necessary for the correct calculation of family mean variances. To obtain these coefficients, additional analyses were made where the respective individual linear models (Eqs. (1–4)) were replaced with equivalent parental family models, but with homogenous residual variances for different trials. From the results of those analyses, the coefficients were extracted by using the DDF function in ASReml 4.1. The principal difference between individual and family models was that parental general combining ability (GCA) term replaced the additive genotype effect term. This implies a variance relationship between GCA variance and additive genetic variance as \( {\hat{\sigma}}_{\mathrm{GCA}}^2=\frac{1}{4}{\hat{\sigma}}_A^2 \) which was accounted for when applying the coefficients for downstream family mean variance calculations. In addition to \( {\hat{V}}_{P_{\overline{\mathrm{HS}}}} \), we also adjusted the equation (Weng et al. 2008) to calculate the phenotypic variance of full-sib family means in full-sib clonal trials as:

Finally, phenotypic variances of clonal means were calculated as:

and

for analyses of full-sib and half-sib clonal trials, respectively. When those estimated variances are related to single-site analyses, genetic-by-location interaction effects were excluded from the four estimated phenotypic variances above.

### 2.9 Heritability estimates

With respect to heritability, \( {h}^2={\hat{V}}_{\mathrm{A}}/{\hat{V}}_{\mathrm{P}} \) is the estimate of individual-tree narrow-sense heritability; \( {H}^2={\hat{V}}_{\mathrm{G}}/{\hat{V}}_{\mathrm{P}} \) is the estimate of individual-tree broad-sense heritability; \( {H}_{\overline{\mathrm{HS}}}^2=\frac{1}{4}{\hat{V}}_{\mathrm{A}}/{\hat{V}}_{P_{\overline{\mathrm{HS}}}} \) is the half-sib family mean heritability; \( {H}_{\overline{\mathrm{FS}}}^2=\left(\frac{1}{2}{\hat{V}}_{\mathrm{A}}+{\hat{\sigma}}_{\mathrm{SCA}}^2\right)/{\hat{V}}_{P_{\overline{\mathrm{FS}}}} \) is the full-sib family mean heritability; and \( {H}_{\overline{\mathrm{CL}}}^2={\hat{V}}_{\mathrm{G}}/{\hat{V}}_{P_{\overline{\mathrm{CL}}}} \) is the clonal mean heritability. Approximate standard errors for variance components and heritabilities were estimated using the Taylor series expansion method in ASReml 4.1 (Gilmour et al. 2015).

### 2.10 Expected genetic gain

The expected genetic gain (*∆G*) as a percentage of the overall trait mean (\( \overline{y} \)) was estimated for five types of selection and deployment strategies (Mullin and Park 1994) under the selection intensity of *i* = 2.67 (i.e., 1%). The selection and deployment strategies included (1) individual forward mass selection and grafted seed orchard (SOPF); (2) individual forward mass selection and cloning (CINF); (3) backward selection of parents (GCA) and polycross (SHSB); (4) backward selection of specific parent pairs and repeat crossing (SFSB); and (5) forward selection based on clonal means followed by further vegetative propagation (CRCF). Based on each selection and deployment strategy, the following formulae for each of expected genetic gain (White et al. 2007) were used:

is the genetic gain associated with SOPF;.

is the genetic gain associated with CINF;

is the genetic gain associated with SHSB;

is the genetic gain associated with SFSB;

is the genetic gain associated with CRCF.

### 2.11 The optimal family size for clonal testing

In order to estimate the optimal family size (number of clones within a full-sib family) for both maximized genetic gain and then maintain recommended diversity (10–20 clones based on simulation (Ingvarsson and Dahlberg 2018)) in clonal forestry (CRCF), random subsets of clones at the range 5, 10, 15, 20, 25, 30, 35, and 40 clones without replacement from each of the 32 full-sib families present in the full-sib clonal trials were sampled to estimate the trend of the genetic gain. To account for the long-term risks associated with the deployment of a few and possibly interrelated clones in Swedish forest stands (Rosvall et al. 2019), four selection strategies were evaluated: (1) selecting the best ten clones without any restrictions on the relationships among the selected clones and (2) selecting the best single clone within each of the ten best families, and (3) selecting the best 20 clones without relationship restrictions and (4) selecting the best single clone within each of the 20 best families. Thus, 32 combinations of family-size and selection strategy were evaluated and the random subset sampling procedure was repeated ten times for each scenario. The genetic gain (%) was calculated as

where *p*_{i} is the fractional contribution to the clonal mixture by the *i*th clone and GV_{i} is the *i*th clonal genetic value (GV = *Z*_{1}*a**+ Z*_{2}*f**+ Z*_{3}*c*, from Eq. (1), which is related to the CRCF strategy with genetic gains predicted by \( \boldsymbol{\Delta }{G}_{H_{\overline{\mathrm{CL}}}^2} \). The mean of *∆G* across the ten repetitions was calculated for each strategy/family-size combination. ASReml-R4 was used to conduct the simulations (Butler et al. 2017).

## 3 Results

### 3.1 Distribution of genetic variance components

Estimates of variance components and genetic parameters for growth and wood quality traits, obtained from four full-sib and two half-sib clonal trials, are presented in Tables 3 and 4. In the full-sib clonal trials, additive genetic variance (\( {\hat{V}}_{\mathrm{A}} \)) accounted for 45.1 to 67.8% of the total genetic variance for the different traits. Dominance variance estimates (\( {\hat{V}}_{\mathrm{D}} \)) were zero for most traits but DBH_{12} and Pilo_{12} which had \( {\hat{V}}_{\mathrm{D}} \) estimates accounting for 10.6% and 54.5% of the total genetic variation, respectively, and their standard errors were larger than their variance components. Epistatic variance estimates of HT_{6}, DBH_{12}, and VOL_{12} varied from 32.2 to 55.2% of the total genetic variance for the different traits while HT_{12} and Pilo_{12} showed no epistatic genetic variation at all. For traits measured at multiple sites, additive-by-location interaction effects (\( {\upsigma}_{AxL}^2 \)), SCA-by-location interaction effects (\( {\upsigma}_{\mathrm{SCA} xL}^2 \)), and clone within family-by-location interaction effects (\( {\upsigma}_{CxL}^2\Big) \) explained only small amounts of the phenotypic variances (0.0–1.8%).

In the two half-sib clonal trials, all growth traits (HT_{3}, HT_{6}, HT_{13}, DBH_{13}, DBH_{25}, and VOL_{13}) showed non-additive genetic variance estimates (\( {\hat{V}}_{\mathrm{N}} \)) that are higher than the corresponding additive genetic variances. For example, \( {\hat{V}}_{\mathrm{N}} \) for HT_{6} accounted for 64.2% of the total genetic variance, which was close to twice the additive variance (35.8%). For wood quality traits, \( {\hat{V}}_{\mathrm{N}} \) was estimated at zero, except for Resi_{25} where \( {\hat{V}}_{\mathrm{N}} \) accounted for 29.0% of the total genetic variance. For growth traits, percentages of the additive-by-site variance \( {\upsigma}_{AxL}^2 \) (4.7–7.4%) were higher than the corresponding percentages for the full-sib trials (0.6–1.8%). The residual variance component was the largest source of variation in both full-sib and half-sib clonal trials for all traits, accounting for 78.6–87.8% of the phenotypic variance for growth traits and 47.4–61.6% of the phenotypic variance for wood quality traits.

### 3.2 Heritability

Five types of heritabilities for different traits are shown in Tables 5 and 6 for full-sib and half-sib clonal trials, respectively. In the full-sib and half-sib clonal trials, clonal mean heritability (\( {H}_{\overline{\mathrm{CL}}}^2 \)) estimates were always larger than individual-tree narrow-sense (*h*^{2}) and broad-sense (*H*^{2}) heritabilities for all traits. \( {H}_{\overline{\mathrm{CL}}}^2 \) varied from 0.22 for DBH_{25} to 0.66 for Resi_{25} (Tables 5 and 6) and *h*^{2} varied from 0.03 for DBH_{25} to 0.43 for Pilo_{25}, while *H*^{2} varied from 0.09 for DBH_{25} to 0.49 for Resi_{25}.The *h*^{2} estimates of most wood quality traits including AV, Resi, Pilo, and MOE (0.21–0.43) were systematically higher than those of growth traits (0.03–0.09).

In full-sib clonal trials, both half-sib (\( {H}_{\overline{\mathrm{HS}}}^2\Big) \) and full-sib (\( {H}_{\overline{\mathrm{FS}}}^2 \)) family mean heritabilities were usually higher than \( {H}_{\overline{\mathrm{CL}}}^2 \) estimates, except for Pilo_{12} where the \( {H}_{\overline{\mathrm{CL}}}^2 \) estimate (0.63) was higher than \( {H}_{\overline{\mathrm{HS}}}^2 \) (0.58). For growth traits measured only in a single full-sib trial (HT_{12} and VOL_{12}), estimates of \( {H}_{\overline{\mathrm{HS}}}^2 \) and \( {H}_{\overline{\mathrm{FS}}}^2 \) were more than twice as high as those of \( {H}_{\overline{\mathrm{CL}}}^2 \). In half-sib clonal trials, however, \( {H}_{\overline{\mathrm{HS}}}^2 \) estimates for growth traits (0.19–0.49) were lower than the corresponding \( {H}_{\overline{\mathrm{CL}}}^2 \) estimates (0.22–0.53), except for HT_{13} measured at a single site. In contrast, several wood quality traits, such as AV_{25}, MOE_{25}, and Resi_{25}, showed that \( {H}_{\overline{\mathrm{HS}}}^2 \) estimates (0.56–0.74) are systematically higher than the corresponding \( {H}_{\overline{\mathrm{CL}}}^2 \) estimates (0.52–0.66).

### 3.3 Genetic gain

Given a common selection intensity of 1%, the expected genetic gains for different traits under five different selection and deployment strategies are shown in Table 7 and b for full-sib and half-sib clonal trials. In full-sib clonal trials, the expected genetic gain (EGG) achieved by mass selection and grafted orchard (SOPF, \( \Delta {G}_{h^2} \)) ranged from 3.7% for HT_{12} to 12.4% for VOL_{12}. EGGs achieved by mass selection and cloning (CINF, \( \Delta {G}_{H^2} \)) ranged from 7.3% for HT_{12} to 18.7% for VOL_{12}. EGGs by clonal mean selection and cloning (CRCF, \( \Delta {G}_{H_{\overline{CL}}^2}\Big) \) varied from 11.3% for HT_{12} to 27.7% for VOL_{12} and were consistently higher than the corresponding EGGs for backward parental (GCA) selection (SHSB, \( \Delta {G}_{H_{\overline{\mathrm{HS}}}^2} \), 5.3 to 19.8%) and EGGs for backward selection of specific controlled crosses (SFSB, \( \Delta {G}_{H_{\overline{\mathrm{FS}}}^2} \), 8.6 to 27.5%).

For the half-sib clonal data analyses, \( \Delta {G}_{h^2} \) of SOPF varied from 2.1% for DBH_{25} to 19.3% for MOE_{25}. For the CINF strategy, \( \Delta {G}_{H^2} \) varied from 4.5% for AV_{25} to 29.5% for VOL_{13}. Finally EGGs for the CRCF strategy (\( \Delta {G}_{H_{\overline{\mathrm{CL}}}^2} \)) varied from 6.1% for AV_{25} to 48.5% for VOL_{13} and were systematically higher than the EGGs for all other strategies available for half-sib clonal trials (SOPF, CINF, and SHSB).

### 3.4 The optimal family size to capture the most clonal gain

To detect the optimal family size considering both genetic gain and the maintenance of genetic diversity at deployment, additional predictions of the genetic gain for DBH_{12} were made based on four putative selection scenarios. The results of this subsampling and selection response predictions are shown in Fig. 1. The realized genetic gain (RGG) achievable decreased from scenario 1 to scenario 4 irrespective of family size but the scenarios exhibited similar trends where the genetic gain increased by an increasing number of available clones per family. The RGG using scenario 1 increased quickly from 15.1% when five clones were available per family to 19.7% when 25 clones/family were available and then a slight further increase to 20.1% when 40 clones/family were available. To consider diversity for deployment using scenarios 3 and 4 compared with scenarios 1 and 2, the RGG decreased, especially for scenario 4 which imposed restrictions on the co-ancestry between selected individuals.

## 4 Discussion

### 4.1 Partitioning of genetic variances and genetic parameters

Traditionally, partitioning epistasis requires either inbred lines or a vegetatively propagated (clonal) population (Foster and Shaw 1988). In forestry, inbred lines have seldom been used as it is usually prohibitively time-consuming to conduct the required multiple generations of inbreeding in tree species (Russell and Ferguson 2008) due to slower generational turnover rate. Moreover, most tree species are known to exhibit severe inbreeding depression (Eriksson et al. 1973; Wu et al. 2016). Thus, clonal field trials with family structure are a more common method to estimate non-additive, in particular, epistatic genetic variance. Currently, genomic data with numerous markers distributed genome-wide using exome capture (Thistlethwaite et al. 2017; Chen et al. 2019), Genotyping-by-sequencing (GBS) (Ratcliffe et al. 2015) or SNP-chips (Tan et al. 2018) are becoming gradually available and genomic relationship matrices for additive, dominance, and epistatic effects can be calculated to estimate the additive, dominance, and epistasis genetic variances. Predicted genomic breeding values can also be compared with the traditional pedigree-based breeding values (Muñoz et al. 2014). Given a multitude of markers, all genetic effects could be estimated and dissected even using a simple half-sib progeny trial of large enough sample size. However, so far, experiments were usually not large in genomic studies with a limited number of individuals within family sampled; the estimates are not accurate enough for studying various deployment strategies.

#### 4.1.1 Growth traits

In this study, large clonal trials with half-sib and full-sib family structures were used to partition the total genetic variance into its basic components (*V*_{A}, *V*_{D}, and *V*_{I}). We found that *V*_{A} accounted for a large proportion of the total genetic variance (45.1–67.8% in full-sib clonal trials and 33.5–48.1% in half-sib trials), indicating that the use of breeding values can capture considerable additive genetic gain. For growth traits in full-sib clonal trials, only DBH showed a small amount of *V*_{D} (10.6%) in relation to *V*_{G}, but epistatic effects accounted for a considerable amount of the total genetic variance (32.2–55.2%, with large standard errors), indicating that the use of clonal deployment strategies (CRCF and maybe also CINF) could greatly improve genetic gain compared with the deployment of seedlings where *V*_{A} (SOPF and SHSB) plus, at best, a fraction of *V*_{D}(SFSB) can be exploited. In contrast, Chen et al. (2019) used both pedigree-based and genomic-based relationship matrices to estimate all genetic effects in Norway spruce and could not find any epistatic variance at all for tree height but instead obtained a considerable *V*_{D} estimate. The result of our study suggests, however, that *V*_{I} for growth traits in Norway spruce is much more important than *V*_{D} in line with other studies such as in white spruce (*Picea glauca* (Moench) Voss) by Weng et al. (2008) and black spruce (*Picea mariana* (Mill.) B.S.P) by Mullin and Park (1994). In a radiata pine (*Pinus radiata* L.) study, DBH at age 5 years showed a large amount of epistatic variance (35.7%) in relation to *V*_{G} (Baltunis et al. 2009). In a loblolly pine (*Pinus taeda* D. Don) study, however, the tree height at age 2 in six trials exhibited a negative estimate of epistasis (Baltunis et al. 2007). In a disconnected partial diallels of coastal Douglas fir (*Pseudotsuga menziesii var. menziesii* DOUGL.) progeny study, Yanchuk (1996) reported that *V*_{D}/*V*_{A} varied from 75 to 260% for three growth traits. The relative comparison of the *V*_{D} and *V*_{I} estimates of our study should be interpreted with caution because the full-sib trial design offered only a limited representation of crosses for each investigated parent (on average 1.3 crosses per parent) and very high standard error thereby making an accurate estimation of *V*_{D} difficult. In conclusion, the non-additive effects seem to be an important part for growth traits and comparable with additive effects in Norway spruce, but the relative amount of variance for dominance and epistatic effects may change with the population, and their size, and the number of full-sib within a test population, and species.

#### 4.1.2 Wood quality traits

For the wood quality traits on three sites, Pilo in full-sib clonal trials and Resi in half-sib clonal trials showed some appreciable non-additive genetic effects (Table 8). Pilo in the full-sib clonal trials showed a *V*_{D} estimate that accounted for 54.5% of the total genetic variance but in the two half-sib clonal trials, the non-additive genetic variance for Pilo was estimated at zero. In *Eucalyptus*, Costa e Silva et al. (2004) reported that epistatic effects explained 22.7% of the total genetic variation for Pilo while dominance genetic effects explained no variation at all. In this study, we did not observe any non-additive effects for acoustic velocity (AV). In contrast, Chen et al. (2019) found small but appreciable dominance and epistatic effects for AV in two full-sib progeny trials of Norway spruce in northern Sweden. It seems there was less non-additive genetic variance in wood quality traits than in growth traits in Norway spruce as similarly observed in other pines (Wu et al. 2008), which may partly be a result of less field environmental error or nursery treatment effects (e.g., sizes of cutting or positions of ortet influence, one of C effects (Burdon and Shelbourne 1974) in forestry) affecting the wood properties compared with growth traits.

### 4.2 Heritability

Both narrow-sense and broad-sense heritabilities are trait-specific and population-specific and are greatly influenced by the homogeneity of the environment containing the genetic test (White et al. 2007). In this study, the individual-tree heritability (*h*^{2} and *H*^{2}) estimates were low for growth traits (0.03–0.09 and 0.09–0.20, respectively), but the large number of clones and families in the full-sib (1430 and 32, respectively) and half-sib (2244 and 138, respectively) trials with low standard errors, nonetheless, indicated that estimates of *h*^{2}and *H*^{2} were accurate (Perron et al. 2013). For wood quality traits, moderate estimates of *h*^{2} and *H*^{2} were in the ranges 0.21–0.43 and 0.29–0.49, respectively, which is in line with a previous study for Norway spruce (Chen et al. 2015). The additive by environmental effect (GEI) was non-significant or low (0–1.8%) from the four full-sib trials, resulting in insignificant impact on the estimation of heritability. The GEI were higher (4.7–9.3%) in two half-sib trials than the corresponding percentages for the full-sib trials. Such GEI may decrease heritability estimates in the joint-site analyses.

#### 4.2.1 Heritability of growth traits

In this study, we used joint-site models to estimate *h*^{2} and *H*^{2} for most the growth traits. Moreover, G × E interaction was accounted for by excluding G × E sources of variance in numerator in the equation of heritability estimates in contrast to many other studies for Norway spruce (Hannrup et al. 2004; Kroon et al. 2011). Estimates of *h*^{2} for tree height were lower (0.04–0.07) than the values of other published results with seedlings using joint-site models for Norway spruce (Chen et al. 2015, 2017; Skrøppa et al. 2015). Estimates of *H*^{2} for tree height were in the range 0.12–0.16 which is comparable with other studies, such as 0.13 for tree height at age 6 (excluding G × E) (Karlsson and Högberg 1998), 0.17–0.29 for tree height at age 9 (including G × E) (Högberg and Karlsson 1998), 0.25–0.29 for tree height at ages 3 to 10 years (Bentzer et al. 1989), and 0.09–0.43 for tree height at ages 9 to 18 years (including G × E) (Högberg and Dutkowski 2010), and 0.10–0.22 for height at age 17 (including G × E). Compared with those published studies with high prior intense selection in the nursery, the difference here may be that only a slight prior selection was performed and that may produce less bias in the estimation of *h*^{2} and *H*^{2}. The *H*^{2} estimates for DBH (0.09–0.12 from ages 12 to 25) and VOL (0.12–0.19 from ages 12 to 13) are also comparable with previously published results (Bentzer et al. 1989) where *H*^{2} was 0.12 for both DBH and VOL at the age of 10.

The total genetic variance and *H*^{2} are expected to be higher than the corresponding additive variance and *h*^{2} if there is a non-additive genetic variance. In the present study, the ratios of *V*_{A}/*V*_{G} ranged from 0.34 to 0.68 (*V*_{N}/*V*_{A} = (*V*_{D} + *V*_{I})/*V*_{A}, 0.50–1.94) for growth traits (tree height, DBH, and VOL) leaving considerable amounts of genetic variation to be explained by non-additive sources. Kroon et al. (2011) reported slightly higher but still relatively small ratios of *V*_{A}/*V*_{G} (0.60–0.84) using three full-sib clonal trials in Norway spruce. These results indicate that non-additive effects are important for the Norway spruce breeding program.

In this study, we found that \( {H}_{\overline{\mathrm{CL}}}^2 \) estimates were lower than numerous corresponding estimates in other conifer species (Mullin and Park 1994; Baltunis et al. 2007, 2009; Weng et al. 2008). The most probable explanation for this discrepancy is that Norway spruce clonal progeny test series in Sweden often only plant approximately 14 ramets per clone, compared with several times this number for test series of other species. Consequently, growth traits that were measured only at a single site, thus encompassing only approximately three ramets per clone. Traits were measured in several trials, such as HT_{6} in full-sib clonal trials. \( {H}_{\overline{\mathrm{CL}}}^2 \) estimates for tree height (0.38–0.60 for tree aged 3 to 6) were comparable with the other results (Lepistö 1993 0.75–0.78, Bentzer et al. 1988 0.63–0.95 for tree aged 1, 0.45–0.73 for tree aged 5), which all excluded G × E effects in the estimate of clonal mean heritability). In our half-sib clonal trials, the \( {H}_{\overline{\mathrm{CL}}}^2 \) estimate for tree height at age 6 years (0.38) was much lower than that at age 3 (0.55). Compared with the higher \( {H}_{\overline{\mathrm{CL}}}^2 \) estimate for HT_{3}, we found that higher estimates of G × E interactions (\( {\upsigma}_{AxL}^2 \) and \( {\upsigma}_{NxL}^2 \)) for HT_{6} as the most probable reason for the lower \( {H}_{\overline{\mathrm{CL}}}^2 \) estimate of this trait. C effects in the nursery might increase \( {H}_{\overline{\mathrm{CL}}}^2 \) at a young age as well. For HT_{6} that \( {\upsigma}_{AxL}^2 \)/\( {\upsigma}_A^2 \)>1 and \( {\upsigma}_{CxL}^2 \)/\( {\upsigma}_C^2 \)=0.42, it suggests that the deployment of such clones in multiple environments, such as in two breeding zones with different requirements of phenology (Berlin et al. 2019), should properly consider G × E interactions.

#### 4.2.2 Heritability of wood quality traits

In the present study, *h*^{2} estimates for wood traits (0.21–0.45) were comparable with those of other published results (0.15–0.50, e.g., AV, Pilo, and MOE) (Hallingbäck et al. 2010; Chen et al. 2015). The ratios of *V*_{A}/*V*_{G} ranged from 0.45 to 1.00 (*V*_{N}/*V*_{A} = (*V*_{D} + *V*_{I})/*V*_{A} at 0.00–0.90 for wood properties and the values were usually lower than those for growth traits, indicating that non-additive effects would contribute less to the prospective increased gains of wood quality traits in a clonal deployment situation than what would be possible for growth traits.

### 4.3 Response for different selection scenarios

To maximize genetic gain is one of the most important objectives for tree improvement programs. However, accurate predictions of genetic gain require reliable estimates of all parameters (Mullin and Park 1992). In this study, we used 1% of the available number of clones as the selection intensity implying a selection of 14 and 24 clones from the full-sib and half-sib clonal populations respectively (CRCF strategy). These numbers have been suggested suitable for clonal forestry in order to maintain genetic diversity in Sweden at the stand level (Rosvall et al. 2019). The estimated genetic gains (EGGs) from such a selection were 17.3 to 20.8% for tree height (if multi-site analyses are considered). These estimates were close to previous estimates for Norway spruce, e.g., 15.0% and 25.0% (Karlsson 1993), 11.5–14.8% for tree height at age 5 years (Bentzer et al. 1988) in Sweden, and 13.1–19.1% in Finland in the same species (Lepistö 1993). Based on previously published documents, Wu (2018) summarizes that in conifer species, an extra genetic gain of 5–25% would be possible from clone testing and deployment, effectively doubling that achievable from backward family-based selection strategies (e.g., SHSB and SFSB). Similarly, a study on radiata pine reported 24% genetic gain, a relative improvement of more than 100% over family forestry (SFSB) (Baltunis et al. 2009). The differences with respect to estimated extra gains were mainly dependent on the size of the population, the number of tested clones within-family, the accuracy of the progeny testing, and the ratio of additive to non-additive genetic variances (Wu 2018).

If the top 1% of replicated clones were selected for vegetative propagation (CRCF), genetic gains for DBH at age 12 years were estimated to reach 19.0%, a relative improvement of 48.4% over full-sib family deployment (SFSB) and 134.6% over half-sib family deployment (SHSB) at the same selection intensity. A 1% selection of replicated clones would also be doable in practice (14 and 24 clones from full-sib and half-sib-trials, respectively). However, with respect to the backward family selection strategies (SFSB, SHSB), it should be noted that in the Swedish breeding program, each breeding population comprise only 50 parents and 50 crosses on average. Thus, a 1% selection intensity would imply a single family being selected and deployed as family forestry which would not be practicable from the point of genetic diversity unless several breeding zones were combined as a single seed orchard. However, for deployment, each seed orchard is usually established with ca. 16–20 parents/clones, thus producing 16–20, or even more, half-sib or full-sib families for deployment (Lindgren and Prescher 2005). For SFSB and SHSB strategies, a realistic selection intensity would be ca. 32–40% (*i* in the range 1.12–1.40) rather than 1% (*i* = 2.67). Thus, from a practical point of view, the relative superiority of CRCF in comparison with SFSB and SHSB would likely be greater in terms of genetic gain than the apparent theoretical estimates shown in Table 7.

In order to estimate the optimal family size for the clonal deployment of the best-selected clones (CRCF strategy), genetic gains were estimated by simulations and clones were sampled into prospective subsets. Genetic gains predicted from this procedure are therefore better comparable with \( \boldsymbol{\Delta }{G}_{H_{\overline{\mathrm{CL}}}^2} \) rather than to \( \Delta {G}_{H^2} \). We found that the realized genetic gain from clonal propagation will be maximized given the availability of 30–40 clones per family, which is also the current Swedish operational design in clonal tests of Norway spruce, both with respect to selecting the best ten/20 clones without any restrictions on the relationships among the selected clones and the best single clone from the best ten/20 families. The maximum realized genetic gain was then found to be 20.2/18.9% and 19.3/16.2% for each selection scenario, respectively. To consider the genetic diversity, a deep simulation based on the current data could be done in the future as other studies (Rosvall et al. 1998; Weng et al. 2010). The study of genomic selection (GS) has been conducted currently in many commercial tree species, such as Norway spruce (Chen et al. 2018), white spruce, loblolly pine, radiate pine, Maritime pine (Bartholomé et al. 2016; Isik et al. 2016), and *Eucalyptus* (Resende et al. 2012; Denis and Bouvet 2013; Tan et al. 2017), but focusing on the seedling progeny trials. If a prospective deployment strategy would entail the establishment of a conventional Norway spruce clonal field test in combination with SE technology (Egertsdotter 2019), the whole process would require ca. 14 years in Sweden, including a nursery stage (2 years), a field test (10 years), and an additional nursery stage (2 years) with implying an annual genetic gain of 1.4% using replicated clonal selection (CRCF). However, it should be noted that if GS methods were applied for clonal forestry, then the process would only require 2 years in the nursery to make rooted cuttings with the maximum two additional months in the laboratory for DNA extraction and genotyping. In such a GS scenario, the annual genetic gain would instead be about 9% if we consider that GS has the same accuracy as a traditional pedigree-based selection. Furthermore, if somatic embryogenesis using cryostorage could be combined with GS, then the time required by the process could be further reduced from 2 years to 6 months in the nursery stage.

## 5 Conclusions

The aim of clonal deployment in forestry is to utilize non-additive genetic effects plus the fact that individual clones can be replicated and therefore better tested in the field. The estimates of genetic parameters and genetic gain for different traits will assist the design of clonal forestry and guide the Norway spruce clonal deployment strategy. Here we conclude that:

- 1)
Both additive and non-additive genetic variances seem to be almost equally important for Norway spruce growth traits.

- 2)
Both individual-tree heritabilities

*h*^{2}and*H*^{2}were higher for wood quality traits (Resi, Pilo, AV, and MOE) than for the corresponding heritabilities for growth traits. \( {H}_{\overline{\mathrm{CL}}}^2 \)–estimates based on clonal means were always the largest in magnitude in comparison with*h*^{2}and*H*^{2}for all studied traits. - 3)
The generation of 30–40 clones per full-sib family will assure the maximum realized genetic at future clonal propagation following a selection scenario of either selecting the best ten/20 clone from the population (no restrictions on co-ancestry) or selecting the single best clone from the best ten/20 families.

- 4)
In most cases, the predicted genetic gains using clonal deployment based on clonal mean selection (CRCF) yielded the largest genetic gains, followed by the deployment of specific full-sib families (SFSB) and seedling deployments based on open-pollinated selected parents (SHSB) or selected individuals (SOPF). If the top 1% of replicated clones were selected for vegetative propagation (CRCF), genetic gains for DBH at age 12 years were estimated to reach 19.0%, a relative improvement of 48.4% over full-sib family deployment (SHSB), and 134.6% over half-sib family deployment (SFSB) at the same selection intensity.

## Data availability

The data sets generated and/or analyzed during the current study are available from the corresponding author and Skogforsk on reasonable request.

## References

Baltunis BS, Huber DA, White TL, Goldfarb B, Stelzer HE (2007) Genetic gain from selection for rooting ability and early growth in vegetatively propagated clones of loblolly pine. Tree Genet Genomes 3:227–238. https://doi.org/10.1007/s11295-006-0058-9

Baltunis BS, Wu HX, Dungey HS, Mullin TJT, Brawner JT (2009) Comparisons of genetic parameters and clonal value predictions from clonal trials and seedling base population trials of radiata pine. Tree Genet Genomes 5:269–278. https://doi.org/10.1007/s11295-008-0172-y

Bartholomé J, Van Heerwaarden J, Isik F, Boury C, Vidal M, Plomion C, Bouffier L (2016) Performance of genomic prediction within and across generations in maritime pine. BMC Genomics 17:604. https://doi.org/10.1186/s12864-016-2879-8

Bentzer BG, Foster GS, Hellberg AR, Podzorski AC (1988) Genotype × environment interaction in Norway spruce involving three levels of genetic control: seed source, clone mixture, and clone. Can J For Res 18:1172–1181. https://doi.org/10.1139/x88-180

Bentzer BG, Foster GS, Hellbersg AR, Podzorski AC (1989) Trends in genetic and environmental parameters, genetic correlations, and response to indirect selection for 10-year volume in a Norway spruce clonal experiment. Can J For Res 19:897–903

Berlin M, Jansson G, Högberg K-A, Helmersson A (2019) Analysis of non-additive genetic effects in Norway spruce. Tree Genet Genomes 15:42–17. https://doi.org/10.1007/s11295-019-1350-9

Bouvet J-M, Makouanzi G, Cros D, Vigneron P (2016) Modeling additive and non-additive effects in a hybrid population using genome-wide genotyping: prediction accuracy implications. Heredity 116:146–157

Brandel G (1990) Volume functions for individual trees; Scots pine (

*Pinus sylvestris*), Norway spruce (*Picea abies*) and birch (*Betula pendula*&*Betula pubescens*)Burdon RD, Shelbourne CJA (1974) The use of vegetative propagules for obtaining genetic information. N Z J For Sci 4:418–425

Butler DG, Cullis BR, Gilmour AR, Gogel BJ, Thompson R (2017) ASReml-R Reference Manual Version 4. VSN International Ltd, Hemel Hempstead, HP1 1ES, UK

Chen Z-Q, García-Gil MR, Karlsson B, Lundqvist S-O, Olsson L, Wu HX (2014) Inheritance of growth and solid wood quality traits in a large Norway spruce population tested at two locations in southern Sweden. Tree Genet Genomes 10:1291–1303. https://doi.org/10.1007/s11295-014-0761-x

Chen Z-Q, Karlsson B, Lundqvist S-O, García-Gil MR, Olsson L, Wu HX (2015) Estimating solid wood properties using Pilodyn and acoustic velocity on standing trees of Norway spruce. Ann For Sci 72:499–508. https://doi.org/10.1007/s13595-015-0458-9

Chen Z-Q, Karlsson B, Wu HX (2017) Patterns of additive genotype-by-environment interaction in tree height of Norway spruce in southern and central Sweden. Tree Genet Genomes 13:25–14. https://doi.org/10.1007/s11295-017-1103-6

Chen Z-Q, Baison J, Pan J, Karlsson B, Andersson B, Westin J, García-Gil MR, Wu HX (2018) Accuracy of genomic selection for growth and wood quality traits in two control-pollinated progeny trials using exome capture as the genotyping platform in Norway spruce. BMC Genomics 19:946. https://doi.org/10.1186/s12864-018-5256-y

Chen Z-Q, Baison J, Pan J, Westin J, García Gil MR, Wu HX. 2019. Increased prediction ability in Norway spruce trials using a marker x environment interaction and non-additive genomic selection model. J. of Heredity, esz061, https://doi.org/10.1093/jhered/esz061

Costa e Silva J, Borralho NG, Potts B (2004) Additive and non-additive genetic parameters from clonally replicated and seedling progenies of

*Eucalyptus globulus*. Theor Appl Genet 108:1113–1119. https://doi.org/10.1007/s00122-003-1524-5Denis M, Bouvet J-M (2013) Efficiency of genomic selection with models including dominance effect in the context of

*Eucalyptus*breeding. Tree Genet Genomes 9:37–51. https://doi.org/10.1007/s11295-012-0528-1Downes GM, Lausberg M, Potts B, Pilbeam D, Bird M, Bradshaw B (2018) Application of the IML Resistograph to the infield assessment of basic density in plantation eucalypts. Aust For 81:177–185

Egertsdotter U (2019) Plant physiological and genetical aspects of the somatic embryogenesis process in conifers. Scand J For Res 34:360–369. https://doi.org/10.1080/02827581.2018.1441433

Eriksson G, Schelander B, Åkebrand V (1973) Inbreeding depression in an old experimental plantation of Picea abies. Hereditas 73:185–193

Falconer D, Mackay T (1996) Introduction to quantitative genetics. Longman, New York

Foster GS, Shaw DV (1988) Using clonal replicates to explore genetic variation in a perennial plant species. Theor Appl Genet 76:788–794. https://doi.org/10.1007/bf00303527

Fundova I, Funda T, Wu HX (2018) Non-destructive wood density assessment of Scots pine (

*Pinus sylvestris*L.) using Resistograph and Pilodyn. PLoS One 13:e0204518Gamal El-Dien O, Ratcliffe B, Klápště J, Porth I, Chen C, El-Kassaby YA (2016) Implementation of the realized genomic relationship matrix to open-pollinated white spruce family testing for disentangling additive from nonadditive genetic effects. G3: Genes, Genom, Genet 6:743–753

Gilmour AR, Gogel BJ, Cullis BR, Welham SJ, Thompson R (2015) ASReml user guide release 4.1. VSN International Ltd, Hemel Hempstead, UK

Hallingbäck HR, Jansson G, Hannrup B (2010) Genetic correlations between spiral grain and growth and quality traits in

*Picea abies*. Can J For Res 40:173–183. https://doi.org/10.1139/x09-173Hannrup B, Cahalan C, Chantre G, Grabner M, Karlsson B, Le Bayon I, Jones GL, Muller U, Pereira H, Rodrigues JC, Rosner S, Rozenberg P, Wilhelmsson L, Wimmer R (2004) Genetic parameters of growth and wood quality traits in

*Picea abies*. Scand J For Res 19:14–29. https://doi.org/10.1080/02827580310019536Högberg K-A, Dutkowski GW (2010) Genetic correlations among field trials of Norway spruce clones representing different propagation cycles. Silvae Genet 59:182–189. https://doi.org/10.1515/sg-2010-0022

Högberg K-A, Karlsson B (1998) Nursery selection of

*Picea abies*clones and effects in field trials. Scand J For Res 13:12–20. https://doi.org/10.1080/02827589809382957Ingvarsson PK, Dahlberg H (2018) The effects of clonal forestry on genetic diversity in wild and domesticated stands of forest trees. Scand J For Res 1–10

Isik K, Kleinschmit J (2005) Similarities and effectiveness of test environments in selecting and deploying desirable genotypes. Theor Appl Genet 110:311–322. https://doi.org/10.1007/s00122-004-1840-4

Isik F, Li B, Frampton J (2003) Estimates of additive, dominance and epistatic genetic variances from a clonally replicated test of loblolly pine. For Sci 49:77–88

Isik K, Kleinschmit J, Steiner W (2010) Age-age correlations and early selection for height in a clonal genetic test of Norway Spruce. For Sci 56:212–221

Isik F, Bartholomé J, Farjat A, Chancerel E, Raffin A, Sanchez L, Plomion C, Bouffier L (2016) Genomic selection in maritime pine. Plant Sci 242:108–119. https://doi.org/10.1016/j.plantsci.2015.08.006

Karlsson B (1993) Twenty years of clonal forestry in Sweden. In Norway spruce provenances and breeding. Proceedings of the IUFRO S22-11 symposium in Latvia Rone V, editor Riga: Lativian Forestry Research Institute; p 208–212

Karlsson B, Högberg K (1998) Genotypic parameters and clone x site interaction in clone tests of Norway spruce (

*Picea abies*(L.) Karst.). For Genet 5:21–30Karlsson B, Wellendorf H, Roulund H, Werner M (2001) Genotype× trial interaction and stability across sites in 11 combined provenance and clone experiments with

*Picea abies*in Denmark and Sweden. Can J For Res 31:1826–1836Kroon J, Ericsson T, Jansson G, Andersson B (2011) Patterns of genetic parameters for height in field genetic tests of

*Picea abies*and*Pinus sylvestris*in Sweden. Tree Genet Genomes 7:1099–1111. https://doi.org/10.1007/s11295-011-0398-yLepistö M (1993) Genetic variation, heritability and expected gain of height in

*Picea abies*in 7 to 9-year-old clonal tests. Scand J For Res 8:480–488Lindgren D, Prescher F (2005) Optimal clone number for seed orchards with tested clones. Silvae Genet 54:80–92

Mullin TJ, Park YS (1992) Estimating genetic gains from alternative breeding strategies for clonal forestry. Can J For Res 22:14–23. https://doi.org/10.1139/x92-003

Mullin TJ, Park YS (1994) Genetic parameters and age–age correlations in a clonally replicated test of black spruce after 10 years. Can J For Res 24:2330–2341. https://doi.org/10.1139/x94-301

Muñoz PR, Resende MFR, Gezan SA, Resende MDV, de los Campos G, Kirst M, Huber D, Peter GF (2014) Unraveling additive from nonadditive effects using genomic relationship matrices. Genetics 198:1759–1768. https://doi.org/10.1534/genetics.114.171322

Perron M, DeBlois J, Desponts M (2013) Use of resampling to assess optimal subgroup composition for estimating genetic parameters from progeny trials. Tree Genet Genomes 9:129–143. https://doi.org/10.1007/s11295-012-0540-5

Ratcliffe B, El-Dien OG, Klápště J, Porth I, Chen C, Jaquish B, El-Kassaby Y (2015) A comparison of genomic selection models across time in interior spruce (

*Picea engelmannii*×*glauca*) using unordered SNP imputation methods. Heredity 115:547–555Resende MDV, Resende MFR, Sansaloni CP, Petroli CD, Missiaggia AA, Aguiar AM, Abad JM, Takahashi EK, Rosado AM, Faria DA, Pappas GJ, Kilian A, Grattapaglia D (2012) Genomic selection for growth and wood quality in

*Eucalyptus*: caspturing the missing heritability and accelerating breeding for complex traits in forest trees. New Phytol 194:116–128. https://doi.org/10.1111/j.1469-8137.2011.04038.xRosvall O, Lindgren D, Mullin TJ (1998) Sustainability robustness and efficiency of a multi-generation breeding strategy based on within-family clonal selection. Silvae Genet 47:307–320

Rosvall O, Ståhl P, Almqvist C, Anderson B, Berlin M, Ericsson T, Eriksson M, Gregorsson B, Hajek J, Hallander J (2011) Review of the Swedish tree breeding programme

Rosvall O, Bradshaw RHW, Egertsdotter U, Ingvarsson PK, Mullin TJ, Wu H (2019) Using Norway spruce clones in Swedish forestry: implications of clones for management. Scand J For Res 1−48. https://doi.org/10.1080/02827581.2019.1590631

Russell JH, Ferguson DC (2008) Preliminary results from five generations of a western redcedar (Thuja plicata) selection study with self-mating. Tree Genet Genomes 4:509–518. https://doi.org/10.1007/s11295-007-0127-8

Skrøppa T, Solheim H, Steffenrem A (2015) Genetic variation, inheritance patterns and parent-offspring relationships after artificial inoculations with Heterobasidion parviporum and Ceratosystis polonica in Norway spruce seed orchards and progeny tests. Silva Fenn 49

Tan B, Grattapaglia D, Martins GS, Ferreira KZ, Sundberg B, Ingvarsson PK (2017) Evaluating the accuracy of genomic prediction of growth and wood traits in two Eucalyptus species and their F1 hybrids. BMC Plant Biol 17:110. https://doi.org/10.1186/s12870-017-1059-6

Tan B, Grattapaglia D, Wu HX, Ingvarsson PK (2018) Genomic relationships reveal significant dominance effects for growth in hybrid Eucalyptus. Plant Sci 267:84–93. https://doi.org/10.1016/j.plantsci.2017.11.011

Thistlethwaite FR, Ratcliffe B, Klápště J, Porth I, Chen C, Stoehr MU, El-Kassaby YA (2017) Genomic prediction accuracies in space and time for height and wood density of Douglas-fir using exome capture as the genotyping platform. BMC Genomics 18:930. https://doi.org/10.1186/s12864-017-4258-5

Weng YH, Park YS, Krasowski MJ, Tosh KJ, Adams G (2008) Partitioning of genetic variance and selection efficiency for alternative vegetative deployment strategies for white spruce in Eastern Canada. Tree Genet Genomes 4:809–819. https://doi.org/10.1007/s11295-008-0154-0

Weng YH, Park YS, Krasowski MJ (2010) Managing genetic gain and diversity in clonal deployment of white spruce in New Brunswick, Canada. Tree Genet Genomes 6:367–376. https://doi.org/10.1007/s11295-009-0255-4

White TL, Adams WT, Neale DB (2007) Forest genetics. CABI publishing, Cambridge

Wu RL (1996) Detecting epistatic genetic variance with a clonally replicated design: models for low- vs high-order nonallelic interaction. Theor Appl Genet 93:102–109. https://doi.org/10.1007/bf00225734

Wu HX (2018) Benefits and risks of using clones in forestry – a review. Scand J For Res 1−8. https://doi.org/10.1080/02827581.2018.1487579

Wu HX, Ivković M, Gapare WJ, Matheson AC, Baltunis BS, Powell MB, McRae TA (2008) Breeding for wood quality and profit in radiata pine: a review of genetic parameters. N Z J For Sci 38:56–87

Wu HX, Hallingbäck HR, Sánchez L (2016) Performance of seven tree breeding strategies under conditions of inbreeding depression. G3: Genes, Genom, Genet 6:529–540. https://doi.org/10.1534/g3.115.025767

Yanchuk A (1996) General and specific combining ability from disconnected partial diallels of coastal Douglas-fir. Silvae Genet 45:37–45

## Acknowledgments

We thank Dr. Ru Wang, Iker Pardo, and Amara Santiesteban for the help of field assistance. We also thank Ola Rosvall for the comments on an earlier version of the manuscript.

## Funding

Open access funding provided by Swedish University of Agricultural Sciences. Financial support was received from Föreningen Skogsträdsförädling (grant number 230-2014-427) and the Swedish Foundation for Strategic Research (SSF, grant number RBP14-0040).

## Author information

### Authors and Affiliations

### Corresponding author

## Ethics declarations

### Conflict of interest

The authors declare that they have no conflict of interest.

## Additional information

Handling Editor: Ricardo Alia

### Publisher’s note

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

### Contribution of the co-authors

HXW, MB, ZQC, AH, and ML initiated the project. ZQC, AF, HNTH, AH, ML, and HRH conducted the data collection. ZQC conducted the data analysis and wrote the manuscript. HXW, HRH, AF, MB, ML, AH, and HNTH contributed to the writing of the manuscript. All the authors read and approved the final manuscript.

## 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 http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Chen, ZQ., Hai, H.N.T., Helmersson, A. *et al.* Advantage of clonal deployment in Norway spruce (*Picea abies* (L.) H. Karst).
*Annals of Forest Science* **77**, 14 (2020). https://doi.org/10.1007/s13595-020-0920-1

Received:

Accepted:

Published:

DOI: https://doi.org/10.1007/s13595-020-0920-1