Comparison of genetic impact on growth and wood traits between seedlings and clones from the same plus trees of Pinus koraiensis

To evaluate the relationships among clones and open pollinated families from the same plus trees and to select elite breeding materials, growth, and wood characteristics of 33-year-old Pinus koraiensis clones and families were measured and analyzed. The results show that growth and wood characters varied significantly. The variation due to clonal effects was higher than that of family effects. The ratio of genetic to phenotypic coefficient of variation of clones in growth and wood traits was above 90%, and the repeatability of these characteristics was more than 0.8, whereas the ratio of genetic to phenotypic coefficient of variation of families was above 90%. The broad-sense heritability of all characteristics exceeded 0.4, and the narrow-sense family heritability of growth traits was less than 0.3. Growth characteristics were positively correlated with each other, but most wood properties were weakly correlated in both clones and families. Fiber length and width were positively correlated between clones and families. Using the membership function method, eleven clones and four families were selected as superior material for improved diameter growth and wood production, and two families from clonal and open-pollinated trees showed consistently better performance. Generally, selection of the best clones is an effective alternative to deployment of families as the repeatability estimates from clonal trees were higher than narrow-sense heritability estimates from open pollinated families. The results provide valuable insight for improving P. koraiensis breeding programs and subsequent genetic improvement.


Introduction
The relationship between parents and offspring in numerous tree species has been studied and research shows that many growth and wood traits can be stably inherited, which provides important opportunities for early selection and genetic improvement Zhang et al. 2020Zhang et al. , 2021Zhang et al. , 2022. The evaluation and selection of elite families is therefore important for genetic improvement of Pinus koraiensis (Sieb. et Zucc). Consequently, the establishment of P. koraiensis seed orchards began since the 1980s (Zhang et al. 2015), and more than 20 seed orchards have been established according to the provenance (Liang 2021). Previous research mainly focused on the establishment of seed orchards (Sun 2006), molecular markers , the evaluation of superior material (Kaviriri et al. 2021), hybridization technology (Mo et al. 2017), seed nutrition (Wang et al. 2019a) and other aspects. Because of the slow growth of P. koraiensis and weak breeding technology, seed orchards of P. koraiensis are still in primary stages. In the evaluation of these orchards, most have parents for upgrading through the evaluation of progeny. However, the parents themselves are also superior individuals selected from the natural forest, and the evaluation and selection of parents are equally important.
P. koraiensis is an economically valuable species of the cold-temperate mixed forest (Ding et al. 1981) and belongs to the Pinaceae family (Zhou et al. 2020). It usually forms a mixed forest with broad-leaved species and is mainly distributed in China, Russia, the Korean Peninsula and Japan, among which the largest cover is in northeast China (Ma et al. 1992). P. koraiensis has ecological and economic values (Liang et al. 2016). In northeast China, it is an important native economic species. It is widely used in construction, shipbuilding, furniture, for planks and timber (Zhao et al. 2018;Liu et al. 2020) owing to its dense, straight trunk, beautiful texture, convenient processing, and is not easily damaged or deformed. Because of their rich fats, proteins, carbohydrates and other components, the seeds are extracted for oil or used as food . In addition, pine seed oil regulates blood lipids (Wang et al. 2019b). Finally, P. koraiensis forests positively influence the environment such as temperature, humidity, water regulation, and wind speed, and root systems stabilize the soil, conserve water, and control erosion (Wang et al. 2019c), so its ecological value is particularly important.
In the past several decades, P. koraiensis forests have been seriously degraded due to over-logging. With the expansion of development of the economy in China, P. koraiensis is often used for afforestation in the construction of the Three-North Shelterbelt program and in the conversion of abandoned farmland to forest. Due to its slow growth, its breeding cycle is lengthy, and therefore, it is important to plan an effective breeding strategy to achieve maximum gain. Forward and backward selection can obtain the highest genetic gain in the short term while reducing costs. Backward selection is a method of selecting superior parents according to their breeding value based on progeny measurement. At present, there have been reports on backward selection of good parents (Zhou et al. 2020). However, the genetic impact on growth and wood traits between seedlings and clones from the same plus trees is not well known and such a study would be important to select the best material for establishing seed orchards. Thus, this research was carried out with the following objectives: (1) to compare the effects of additive and non-additive effects on seedling traits by estimating the broad-sense heritability and narrow-sense heritability of each trait in open-pollinated families and the repeatability of each trait in the clones; (2) to compare the coefficients of variation and correlation between seedlings and grafted plants, and the genetic variation of growth and wood characters between parents and offspring; and, (3) to select the best clones and open-pollinated families in terms of growth and wood properties so as to provide the basis for breeding improvement. To achieve these objectives, the variations in growth traits and wood properties among 35 P. koraiensis parental clones and open-pollinated families were examined and elite clones and families selected.

Study site
The research was carried out in a P. koraiensis seed orchard in Lushuihe, Fusong County (127°29′ E, 42°29′ N), Baishan City, Jilin province. The seed orchard is located in the north temperate East Asian monsoon climate, with an average altitude of 780 m, average annual and maximum temperatures of 2.7 °C and 32.5 °C, respectively, a minimum temperature of − 44.4 °C, annual average precipitation of 871 mm, and annual average evaporation of 1166.44 mm. Annual sunshine hours are 2117 and the annual frost-free period is 115 days.

Plant materials
Thirty-five P. koraiensis clones and open-pollinated families were used as experimental materials. The plus trees were selected from the natural forest of Lushuihe Forestry Bureau of Jilin province in 1984, and seeds and scions collected in the autumn of the same year. The scions were grafted on local 4-year-old P. koraiensis rootstock for establishment of the seed orchard in the spring of 1985 in a randomized complete block design, with 10 blocks of 10 single plants of the same clone, each block containing one plant at 3 m × 3 m spacing. The seeds were germinated in the following year. The soil conditions used for seedlings and grafted materials were similar. In 1989, 4-year-old seedlings were used for the establishment of the progeny trail in a randomized complete block design with four blocks containing 24 single seedlings of the same family, each block containing 6 plants at 2 m × 2 m spacing.

Measurement of growth characteristics
In the autumn of 2017, the growth properties of 35 clones and 35 open-pollinated families were measured, including height (H) and diameter at breast height (DBH) for clones, and for open-pollinated families, ground diameter (GD), DBH, diameter at 3 m height (D 3 ), crown width (CR), stem straightness (SS), branch angle (BA), height to 1st branch whorl l from the ground (1BH), and height of the 6th branch whorl (6BH) from the ground for open pollinated families. Due to the high density of the seed orchard, height of individuals of open-pollinated families was not recorded. Height (H) was measured using the Vertex IV ultrasonic height range finder, and DBH, GD and D3 by a tape and volumes were calculated. Crown widths in the east-west and north-south directions were measured by tape. The 1BH and 6BH were measured by a sliding staff. The SS and BA were determined following the methods of Zhao (2010) and Liang et al. (2018a).

Measurement of wood properties
In the same year, wood cores from 35 clones in the first five blocks were taken at a height of 1.3 m in the north-south direction of each tree. Wood cores from the 35 open-pollinated families in the first block were also taken in the same pattern from each tree. The cores were marked, put in paper bags, and brought to the laboratory for determining wood traits. Wood density (WD) was determined by the drainage method (Yin et al. 2017) using five cores for each clone and six per family. Fiber lengths (FL) and widths (FW) were measured on five and six cores per clone and family, respectively, according to Yin (2017). Each clone and family consisted of 30 FLs and 30 FWs.

Statistical analysis
For each tree, volume (V) was calculated following Zhao et al. (2013) as follow: where V is volume, D is diameter at breast height, and H is tree height.
Analysis of variance (ANOVA) was used to examine the variations in growth and wood properties among clones and open-pollinated families  using the following model: (1) where μ is the overall mean, C i (F i ) the clone (family) effect, B j the block effect, FB ij the effect of family i within block j, and e ij is the random residual. Clones, families and interactions between blocks and families were treated as random factors and blocks as fixed factors. For each growth and wood trait, variance components were estimated for each source of variation based on the maximum likelihood method, and broad-sense heritability was estimated according to Xu (2006) as follows: where F value is from analysis of variance. The estimation of narrow-sense heritability (h 2 ) of a characteristic for the families was calculated according to Xu (2006) as follows: where 2 F and 2 e are the family and residual variances, and 2 FB the variance component of interaction between family and block. Since the wood cores in the same block were used for sampling, there was no block effect, the effect of environmental differences is small, and the results of narrow-sense heritability are unsatisfactory. Therefore, the narrow-sense heritability of wood traits was not analyzed. The repeatability (R) of clones of all investigated characteristics was calculated as follows (Zhao et al. 2015): where 2 d is the genetic variance component between clones and 2 e is the residual variance component. The phenotypic coefficient of variation (PCV) and genotypic coefficient of variation (GCV) were calculated using the following formula : where 2 p and 2 g are the phenotypic variance and genetic variance of a character, respectively, X is the phenotypic mean of a character.
The phenotype correlation r p12 of characters x and y was calculated as follows (Guerra et al. 2016): where Cov p12 denotes the phenotypic covariance between characters of x and y, and 2 p1 and 2 P2 denote the phenotypic variance of a character x and y, respectively.
The general combining ability was calculated as follows  where g is the general combining ability of the parent, x is the average value of the offspring of a certain mating combination of a parent in a certain trait, u is the total mean value of all combinations of this trait.
For selection of the best clones and families in growth and in wood production, the membership function of different clones (families) was calculated using DBH, WD and FL as selection indexes. As fiber length and fiber width were positively correlated, only fiber length was considered in the evaluation. The membership function was calculated as follows (Liu et al. 2017): When the trait was negatively correlated with wood physical properties: where Xij is the measured value for the j character of clone (family) i, X jmax is the maximum value of the character, X jmin is the minimum value of the character. All data were analyzed using SPSS 19.0 software (IBM Corp., Armonk, NY, USA).

Variations in growth traits
Highly significant variations were detected among parental clones and open-pollinated families ( Table 1). Estimates of variance show that the clonal effect was substantially higher than family effects (Table 2). They also show that family variation in growth traits was very low compared to the interaction between families and blocks as well as clonal effects. Descriptive statistics for each trait are shown in Table 3. Among clones, height ranged from 5.6 m to 13.2 m, DBH from 10.7 cm to 24.8 cm, and volume from 0.03 m 3 to 0.19 m 3 . Among open-pollinated families, GD (ground diameters) ranged from 15.4 cm to 24.6 cm, DBH 12.4 cm to 20.4 cm, the D3 (diameter at 3 m) from 11.2 cm to 18.6 cm, the BA (branch angle) from 130.8° to 157.1°, the CR (crown width) from 3.2 m to 5.2 m, SS (stem straightness) from 3.9 to 4.7, the 1BH (1st branch whorl) from 1.7 m to 1.9 m, and the 6BH (6th branch whorl) from 3.8 m to 4.4 m.
The phenotypic coefficient of variation (PCV) and the genotypic coefficient variation (GCV) were higher for clones than for open-pollinated families. Values of PCV and GCV among clones were closer as compared to values for open-pollinated families, suggesting that most of the variations in growth traits were genetically controlled. This is further evidenced from the high repeatability of each trait (R > 0.90) compared to the low narrow-sense heritability (h 2 < 0.3). However, family broad-sense heritability was high (H 2 ≥ 0.66) for all traits except stem straightness and height at the 1st branch whorl for which heritability was less than 50% (Table 3). Due to the high planting density of the materials in this study, the heritability and coefficient of variation of 1BH were not considered.

Variations in wood properties
Significant variations in wood properties were observed among parental clones and open-pollinated families ( Table 1). Estimates of variance showed that clonal variations in wood properties was very high (78.6% − 86.2%) compared to family variations (20.0% − 53.0%; Table 2). The average FL ranged from 1855.47 µm to 2583.41 µm, the average FW from 29.47 µm to 41.93 µm, the average WD from 0.44 g cm −3 to 0.62 g cm −3 among parental clones (Table 3). Among open-pollinated families, the average FL ranged from 1314.86 µm to 2112.12 µm, the average FW from 34.76 µm to 50.43 µm, and the average WD from 0.47 g cm −3 to 0.57 g cm −3 .
The phenotypic coefficient of variation (PCV) and the genotypic coefficient of variation (GCV) were slightly higher for open-pollinated families than for clones, particularly for fiber length (FL) and fiber width (FW). The repeatability of FL and FW was slightly lower (R = 0.79 − 0.86) for clones than for family broad-sense heritability (H 2 = 0.96-0.97), but repeatability of WD was slightly higher than family broadsense heritability (R = 0.79, H 2 = 0.71; Table 3).

Correlation analysis
The correlation of height and diameter was highly significant (P < 0.01; 0.46 < r < 0.99) and positive for both clones and open-pollinated families (Table 4). For wood properties, there was a highly significant positive correlation between fiber length and fiber width in both clones and open-pollinated families (0.16 < r < 0.26). In open-pollinated families, stem form traits were significantly correlated with growth traits (-0.14 < r < 0.63). The correlation between growth traits and wood properties was insignificant in both clones and open-pollinated families. Most traits were not significantly correlated between clones and open-pollinated families. However, there were positive correlations among wood properties except for wood density (0.11 < r < 0.22).

General combining ability
The general combining ability of various traits of P. koraiensis open-pollinated families is shown in Table 5. The general combining ability of DBH of each family ranged from − 4.11 to 3.88, and it was higher for PK007, PK145, PK036 and PK429 than for the other families. This shows that using the parents of these families as hybrid parents, the genetic improvement of the progeny's DBH was greater. The general combining ability of CR (crown width) of each family ranged from − 1.03 to 1.00, being higher for PK150, PK025, PK429 and PK145 than for the other families. For SS (stem straightness) of each family, it varied from − 0.43 to 0.28; thus the clonal parents of PK424, PK066, PK064 and PK025 could be used as the hybridization parents, and the improvement of stem straightness in offspring will be stronger. The general combining ability of FL of each family ranged from − 512.99 to 284.27, and the clones PK036, PK004, PK014 and PK407 could be used as the hybridization parents. Note: H, tree height (m); DBH, diameter at breast height (cm); V, volume (m 3 ); FL, fiber length (µm); FW, fiber width (µm); WD, wood density (g cm −3 ); GD, ground diameter (cm); D 3 , diameter at 3 m height (cm); BA, branch angle (°); CR, crown width (m); SS, stem straightness; 1BH, height at 1 st branch whorl (m); 6BH, height at 6 th branch whorl (m)

Membership function comprehensive evaluation
The comprehensive evaluation index of the fuzzy mathematical membership function method was used to select superior clones and families with large DBH, high wood density and long fiber length (Table 6). According to the method, the standardized value of each trait is multiplied by the weight of each trait (0.4 for DBH and 0.3 for wood density and fiber length, according to Liu et al. 2017) in the evaluation system, and the membership function mean value of each trait for clones and families was obtained. The superior clones and families with average membership functions > 0.6 were selected. Consequently, eleven clones (PK124, PK051, PK042, PK003, PK415, PK036, PK162, PK303, PK004, PK025, PK423) and four families (PK003, PK066, PK124, PK401) were selected as superior breeding material for diameter growth and wood production in the establishment of Korean pine forests, as well as for the H tree height (m); DBH diameter at breast height (cm); V volume (m 3 ), FL fiber length (µm); FW fiber width (µm); WD wood density (g cm −3 ); GD ground diameter (cm); D 3 diameter at 3 m height (cm), BA branch angle (°); CR crown width (m); SS stem straightness; 1BH, height at 1st branch whorl (m); 6BH height at 6th branch whorl (m); SD standard deviation, PCV phenotypic coefficient of variation, GCV genotypic coefficient of variation; R repeatability; H 2 broad-sense heritability, h 2 narrow-sense heritability

Discussion
Genetic variation is the driving factor of tree breeding (Mwase et al. 2008). In this study, analysis of variance showed highly significant variations among clones (P < 0.01) and open-pollinated families (P < 0.01) in growth and wood properties. This indicates that there was a great variation in the clones and families that can be used for further selection of superior breeding material. The phenotypic and genotypic coefficients of variation (PCV and GCV, respectively) reflect the degree of variation of each trait, and a high value indicates the suitability for selection of superior clones and families (Metougui et al. 2017). In the present study, the PCVs of clones ranged from 7.3 to 37.2%, and the GCVs from 6.5 to 36.3%; volume had the largest PCV and GCV. These results are similar to those of Xu et al. (2012) for Olga Bay larch (Larix olgensis A. Henry), where the improvement potential was high Table 5 General combining ability of phenotypic traits in P. koraiensis clones Note: H, tree height (m); DBH, diameter at breast height (cm); V, volume (m 3 ); FL, fiber length (µm); FW, fiber width (µm); WD, wood density (g cm −3 ); GD, ground diameter (cm); D 3 , diameter at 3 m height (cm); BA, branch angle (°); CR, crown width (m); SS, stem straightness; 1BH, height at 1 st branch whorl (m); 6BH, height at 6 th branch whorl (m) when volume was used as an index. The PCV of family traits ranged from 3.2 to 12.4%, the GCVs from 0.5 to 10.2%; among them, crown width had the largest PCV, which is similar to the results of Chen (2016) for Eucalyptus grandis W. Hill. This indicates that there is great potential for improvement using crown width as an index to evaluate and select. The ratio of genetic and phenotypic coefficients of variation of growth and wood traits among clones was more than 90%. The same results were obtained for families for wood traits, which are higher than those reported by Ji et al. (2007) for Chinese red pine (Pinus massoniana Lamb.). Overall, these results show that the growth and wood property traits were more genetically controlled. This study shows that the repeatability of clones was higher than 75.0%, higher than the results of Jiang et al. (2019) for 38-year-old P. koraiensis. The broad-sense heritability of open-pollinated families was > 50% except for stem straightness, which is similar to a previous study on 24-yearold P. koraiensis (Wang et al. 2019d). These results show that all traits of different clones and families could be stably inherited, as repeatability is an index to measure trait stability. The higher the repeatability, the less the trait is affected by the external environment (Wu and Mao 2012). Heritability is also an index that informs about the ability of a parent to pass on a trait to their offspring. Higher heritability means that the trait is more stable to be passed on (Maniee et al. 1998). In this study, for growth traits, the repeatability of clones were higher than the family broad-sense as well as narrow-sense heritability. For wood traits, particularly for fiber length and width, repeatability of clones was slightly lower than the family broad-sense heritability. Therefore, the selection of clonal material is a valuable alternative for the establishment or upgrading of seed orchards. The broadsense heritability of families for wood properties was higher than that for growth characteristics, indicating that wood characters are strongly genetically controlled. The narrowsense heritability of families for growth traits is lower than 0.3, implying that non-additive effects dominate variations in growth traits.
The correlation between traits is important in the process of tree improvement, and can provide reference for comprehensive evaluation of multiple traits. In this study, the correlation coefficient between DBH and volume was greater than 0.9, indicating that they were not inherited independently. Our result is consistent with Liang et al. (2018b) on 32-year-old P. koraiensis clones. In openpollinated families, most of the growth traits were highly positive correlated and stem straightness was negatively correlated with the ground diameter, diameter at breast height, the diameter at the 3-m height, and the crown width. The greater the stem straightness, the smaller all these were, which may be caused by the high stand density and the longitudinal growth of the trees. Correlations between growth traits and wood properties were not significant in both clones and open-pollinated families, indicating that these are not synchronized, which are similar to results by Zhou et al. (1991) and Liang et al. (2016). Some results indicate that the growth rate may not significantly affect wood properties (Zhou et al. 1991). The correlation of growth traits between clones and open-pollinated families was insignificant, and the heritability of DBH was lower than repeatability. The phenotypic variation was mainly influenced by heredity and environment. In this study, growth traits of families were more influenced by environment; studies have shown that phenotypic traits are more easily affected by environmental factors (Liu et al. 2021). Wood physical properties, except for density, showed positive correlations between clones and open-pollinated families, which is similar to Wang (2020) for L. olgensis. The heritability of wood traits was higher than repeatability, indicating that wood properties may be strongly genetically controlled. Thus, superior breeding material can be selected by the combined analysis of clones and families for subsequent research. General combining ability (GCA) refers to the relative ability of parents to pass elite traits on to their offspring (Zhou et al. 2004;Moterle et al. 2012). Superior parents can be selected through progeny testing, which is widely used in crop breeding. In this study, 35 open-pollinated families were analyzed for GCA, and the parents were evaluated and selected with a 10% selection rate. When considering wood production as the main breeding objective, DBH used as an evaluation index, the general combining ability of PK7, PK36, PK145 and PK429 families was higher. For landscaping, crown width and stem straightness are used as evaluation indexes; PK25 performed well. When highquality wood is considered as a breeding objective, fiber length can be used as evaluation index; families such as PK4, PK14, PK36 and PK407 had higher GCA, indicating that the genetic improvement of fiber length was stronger in offspring. The selected superior materials can provide references for subsequent research.
Because the traits with high general combining ability (GCA) were different among families, it is difficult to combine the traits to screen. Therefore, the membership function comprehensive evaluation method was used to evaluate and select clones and open-pollinated families separately. With DBH, fiber length, and wood density as evaluation indexes, superior clones and families were selected. The superior clones were PK124, PK051, PK042, PK003, PK415, PK036, PK162, PK303, PK004, PK025, and PK423. The superior families were PK003, PK066, PK124 and PK401. PK124 and PK003. These clones and families selected in this study can be used in the establishment of plantations to produce furniture and high-quality wood, and provide superior materials for genetic improvement of P. koraiensis. However, since the repeatability estimates from the clonal trees were higher than h 2 estimates, selection of the best clones would be an effective alternative.
Recently, with further advancement of forest genetics, the establishment of high-quality seed orchards has great significance. In this study, the clones and families of PK124 and PK003 performed well. The superior families provide the basis for the reconstruction and upgrading of the seed orchard. Based on breeding objectives, superior clones can be propagated asexually to obtain higher ecological and economic value, and provide materials for the popularization and establishment of high-quality P. koraiensis plantations. These results also provide reference for backward and forward selection for P. koraiensis propagation.

Conclusions
The results demonstrate that the genetic impact on growth traits of grafted plants is higher than seedlings of openpollinated families as shown by their higher repeatability. Correlation in growth traits between clones and open-pollinated families were insignificant but significant in wood traits. Considering diameter at breast height, wood density, and fiber length as indexes, a comprehensive evaluation of clones and families was carried out by the membership function method to select superior families. Consequently, two families, PK003 and PK124, performed well in both clones and open-pollinated families. Thus, these selected clones and families can provide the foundation for establishing or upgrading P. koraiensis seed orchards to provide material for the establishment of high-quality P. koraiensis plantations in the future.
Author contributions Xiaoting Liu, Zhimin Lu and Jingui Wu contributed equally to this work.
Funding Open access funding provided by Swedish University of Agricultural Sciences. The work was supported by the Jilin Province Science and Technology Department Key Research and Development Project (20200402094 N C).
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.