Estimation of Genetic Parameters and Wood Yield Selection Index in a Clonal Trial of Korean Pine ( Pinus koraiensis ) in Northeastern China

: In order to determine suitable traits for selecting high-wood-yield Korean pine materials, eleven morphological characteristics (tree height, basal diameter, diameter at breast height, diameter at 3 meter height, stem straightness degree, crown breadth, crown height, branch angle, branch number per node, bark thickness, and stem volume) were investigated in a 38-year-old Korean pine clonal trial at Naozhi orchard. A statistical approach combining variance and regression analysis was used to extract appropriate traits for selecting elite clones. Results of variance analysis showed signiﬁcant difference in variance sources in most of the traits, except for the stem straightness degree, which had a p -value of 0.94. Moderate to high coefﬁcients of variation and clonal repeatability ranged from 10.73% to 35.45% and from 0.06% to 0.78%, respectively. Strong signiﬁcant correlations on the phenotypic and genotypic levels were observed between the straightness traits and tree volume, but crown breadth was weakly correlated to the volume. Four principal components retaining up to 80% of the total variation were extracted, and stem volume, basal diameter, diameter at breast height, diameter at 3 meter height, tree height, and crown height displayed high correlation to these components (r ranged from 0.76 to 0.98). Based on the Type III sum of squares, tree height, diameter at breast height, and branch number showed signiﬁcant information to explain the clonal variability based on stem volume. Using the extracted characteristics as the selection index, six clones (PK105, PK59, PK104, PK36, PK28, and K101) displayed the highest Qi values, with a selection rate of 5% corresponding to the genetic gain of 42.96% in stem volume. This study provides beneﬁcial information for the selection of multiple traits for genetically improved genotypes of Korean pine.


Introduction
The Korean pine breeding program remains among the largest coniferous tree improvement programs in Northeastern China that are carried out [1]. Pinus koraiensis was primarily selected to increase its wood yield. A few years later, however, seed yield was taken into account, responding to users' demand for the seeds and nuts of Korean pine and following the official decree in 2016 against the cutting of natural forests. Korean pine offers many economic benefits. Its wood is more appreciated in construction, and its seeds are variously used in diet and medicine [2,3].
P. koraiensis is much desired, being preferred in community reforestation projects following the services it provided in Northeastern China [4]. However, the planting of stock with desirable genetic properties is still lower than the planting demand [5]. Because of the limitation of asexual propagation techniques due to the difficulty of producing large amounts of planting materials by cuttings, tissue culture, somatic embryogenesis, and other asexual methods, seed orchards took an important place in improving Korean pine [6].
The improvement of Korean pine dates back more than a half century [7]. To date, significant gains have been realized, and several studies reported some improvements in growth rate, and seed and wood quality [8][9][10][11]. Facing multiple needs and a permanent global change, multi-objective approaches have been suggested for the sustainable management of forest resources [12]. In tree improvement, selection for multiple characteristics has been introduced to meet multiple purposes.
From this perspective, the establishment of selection criteria in a tree improvement program remains a critical step for choosing appropriate individual trees on which future breeding objectives will focus [13]. Several selection methods, for a particular trait or multiple traits, can be used for the purpose of improvement. Among these methods, the comprehensive evaluation by Qi value, for instance, assigns importance to the weighting of correlated characteristics when simultaneously enhancing various traits [14,15]. On another side, Leakey and Page [16] referred to the breeding ideotypes approach, which consists in the combination of biomass-related traits to define a hypothetical model of suitable plants producing an economic yield.
Multiple traits selection has been applied to select superior Korean pine materials in different phases for its breeding, including clonal, provenance, and progeny tests [17,18]. The selection of multiple traits presents some advantages over tree selection for a single particular characteristic [19], especially when the characteristics are positively correlated and display high heritability [20]. However, correlation coefficients, although important for indirect selection, are not sufficient to designate more specific traits to be used as the selection index because of positive and negative effects between linked traits, due principally to uncontrollable and complex trait interaction [21].
Pathway analysis, consisting of measurement of the weight of variance of dependent variables and independent ones [22], involves the conversion of correlation coefficients into direct or indirect effects through a multiple regression equation that calculates the impact of correlated variables with precision. Pathway analysis has revealed associations between cause and effect among yield and yield-related traits of genotypes in agronomic plants [23], as well as different growing environments or habitats [24,25]. By reducing a large set of variables into underlying correlated factors [26], principal component analysis (PCA) appears to be a useful tool to highlight trait combinations that greatly contribute to the variation among tree populations under characterization [27,28].
Combining correlation, pathway analysis, and PCA in a multiple-trait comprehensive evaluation could be a powerful statistical method for selecting suitable values in multiple traits to achieve improvement goals.
Given that first generation Korean pine was selected for timber production, and that stem volume represents the major part of the above-ground tree biomass, knowledge of any effects of other growth characteristics on volume would allow the prediction of various gains from the orchard. As there has been little research previously on these aspects, especially on a multiple-usage tree species, this study tries to determine the effects of other growth traits on the volume of Korean pine. It specifically aims to (1) investigate the growth performances of different clones; (2) establish correlations between growth characteristics and their effects on wood yield; (3) highlight the optimum number of traits proper for a comprehensive evaluation index to select superior clones; and (4) select elite clones with high clonal values for future breeding purposes. This work is written from the perspective of setting up a basic theory for the efficient selection of trees for genetic improvement.

Data Collection and Trial Design
The morphological characteristics were measured in Naozhi forestry seed orchard (41 • 05 N, 126 • 06 E), located on the western hillside of Changbai Mountain in Linjiang town, Northeastern China ( Figure 1). This region belongs to the northern cold-temperate continental monsoon climate; summer is warm and rainy, and winter is cold and long. The annual average temperature varies between 4 and 6 • C and the annual average rainfall range between 750 and 1000 mm. The frost-free period is 135 days. The soil type corresponds to the Albi-Boric Argosols, according to the American soil taxonomy [29], and it is dominated by dark brown soil more than 40 cm thick with important textural constituents of sand (15.13%), silt (63.31%), and clay (21.56%) [30].
The annual average temperature varies between 4 and 6 °C and the annual average rainfall range between 750 and 1000 mm. The frost-free period is 135 days. The soil type corresponds to the Albi-Boric Argosols, according to the American soil taxonomy [29], and it is dominated by dark brown soil more than 40 cm thick with important textural constituents of sand (15.13%), silt (63.31%), and clay (21.56%) [30].
One hundred and ten clones (named from PK01 to PK110; see Table 11) were used as experimental materials. These clones were chosen, based on their growth performances, from the natural distribution area of P. koraiensis in Linjiang town in Jilin Province in 1978. The materials were grafted and planted with 4-year-old seedlings in the spring of 1980. The experimental design consisted of a completely randomized block, with 10 blocks containing 110 different clones, and one ramet of each clone planted with inter-plant and row spacing of 7.0 m × 7.0 m.

Growth Trait Measurement
The following traits were measured in living ramets of each clone in July 2018: tree height (Ht), basal diameter (BD), diameter at breast height (DBH), diameter at 3 meter height (DIAM3), stem straightness degree (SSD), crown breadth (CB), crown height (CH), branch number per node (BNN), branch angle (BA), and bark thickness (BT). Ht and CH One hundred and ten clones (named from PK01 to PK110; see Table 11) were used as experimental materials. These clones were chosen, based on their growth performances, from the natural distribution area of P. koraiensis in Linjiang town in Jilin Province in 1978. The materials were grafted and planted with 4-year-old seedlings in the spring of 1980. The experimental design consisted of a completely randomized block, with 10 blocks containing 110 different clones, and one ramet of each clone planted with inter-plant and row spacing of 7.0 m × 7.0 m.

Growth Trait Measurement
The following traits were measured in living ramets of each clone in July 2018: tree height (Ht), basal diameter (BD), diameter at breast height (DBH), diameter at 3 m height (DIAM3), stem straightness degree (SSD), crown breadth (CB), crown height (CH), branch number per node (BNN), branch angle (BA), and bark thickness (BT). Ht and CH were measured with a meter ruler; BD, DBH, and DIAM3 were measured by a caliper; SSD was estimated and squared before calculation according to Zhao et al. [31]. CB was calculated by averaging the values of north-south and east-west crown width, and CH was obtained by subtracting the value of tree height at the first branch from the total tree height. The average number of BA was measured with a protractor at 3 m height. The number of branches (BNN) was counted on the first five nodes. Pieces of bark were removed at the base of the trunk and at a 1.3 m height of the stem, and were measured by caliper. Their average gave the BT. The stem volume (Vol.) of each tree was calculated according to Louppe's method [32], using the coefficient of 0.41 for coniferous trees in the formula proposed by the following equation:

Statistical Analysis
All the data were analyzed using XLSTAT 2019 [33], and IBM SPSS version 20 [34,35]. Firstly, the traits were tested for normality by Shapiro-Wilk's test, and standardized by subtracting the mean value and dividing by the standard deviation to meet the assumption of the normal distribution of residuals. Average values, ranges, and standard error (SE) were calculated. To estimate clonal/genotype and random/environmental effects, the mixed linear model (Equation (2)) was employed to join together the fixed and random effects, and F-tests of one-way analysis of variance were used to test the significance in the variance sources [36].
where P ij is the phenotypic performance of an individual tree j of clone i, µ is the overall mean, a ij is the random clonal effects, b i is the fixed effect of the block, and ε ij is the residual value of the trait. Because each block contained a unique ramet from each clone, there were not computable interaction effects of clones and blocks from the mixed model. Clonal values were calculated through the relationship in Equation (3), and the variance component was estimated using the relationship in the Equation (4), considering the micro-environmental conditions in the orchard to be equal.
where σ 2 Yij , σ 2 ai , and σ 2 ε symbolize the total phenotypic variance, genotypic variance, and variance of within-plot error, respectively [37]. A discriminatory factor analysis using PCA was done to determine traits that retained most of the variation in investigated clones [38]. The PCA analysis aimed to summarize the information contained in the eleven original traits into a smaller set of new composite variants (principal components) with a minimum loss of information, and to discover the relationship between new factors with single traits [39]. Standardized data were used prior to PCA regarding different scales of units of variables. The repeatability of clones was estimated using the clonal variance component divided by the total phenotypic variance that was equal to the subtraction of the inverse of F (Fisher-Yates) coefficients in ANOVA (1/F) from the unit value (1) as indicated by Zheng et al. [40].
To determine correlation of traits and determine whether multiple-trait selection supported the genetic breeding purpose, Pearson's correlation coefficient r A (xy) was calculated using phenotypic values, and broad-sense genetic correlation was computed using genotypic variance of clones for all phenotypic traits according to Porth et al. [41] and Pan et al. [42], as per following equations: r g (xy) = COV g(x,y) where COV P(x,y) and COV g(x,y) are covariance between phenotypic and genetic values of traits x and y; and σ p (x) and σ p (y) are variance components of traits x and y, respectively.
To determine the extent of phenotypic variation, the coefficient of variation (PCV) was calculated by dividing the standard deviation by mean values from each trait, following Hai et al. [43]: where X and SD were the mean value and standard deviation, respectively, for a given trait.
To determine the effects of growth characteristics on the wood yield of clones, a regression analysis was performed taking tree volume as the independent variable in the linear modal (Equation (9)). Trait contributions were then deduced from regression coefficients according to Couto et al. [44].
where Y represents the dependent variable, the stem volume; X 1 , X 2 . . . , X n were the explanatory independent variables representing the other growth traits; and p 1 , p 2 . . . , p n were regression coefficients. The pathway diagram showing trait effects was constructed in Amos Software, IBM SPSS Amos 23 (http://www-03.ibm.com/spss-amos, accessed on 21 December 2020) [35]. Multicollinearity statistics of variables were done to detect variables with identical effects. Twelve statistical parameters-including R 2 , adjusted R 2 , mean square error (MSE), root mean square error (RMSE), mean absolute percentage error (MAPE), Durbin Watson statistic (DW), Mallows' Cp, Akaike information criterion (AIC), Bayesian information criterion (SBC), predicted residual error sum of squares (PRESS), and quality index (Q)-were computed to verify the goodness of the fit of statistics of the regression model. A comprehensive evaluation (Equation (10)) was made to rank the clones for selection using clonal phenotypic performances in the Qi value formula according to Zhao et al. [31]: where Q i is the comprehensive evaluation value, a i = X ij /X jmax , X ij is the mean value of a given trait, and X jmax is the maximum mean value of the trait. To estimate the genetic gain from the selected clones, the formula developed by Francisco et al. [45] was applied (Equation (11)). Considering P. koraiensis logs, the expected economic gain was deduced by extrapolating the unitary economic value per cubic meter of wood multiplied by the logged production and realized gain in stem volume (Equation (12)).
∆Eco. = Price × volume unit × realized gain in stem vol.

ANOVA and Variation Parameters
To characterize the phenotypic variation of different clones of Korean pine in the Naozhi orchard, 11 morphological traits were measured. Significant differences were found in all the investigated traits, providing a selection basis for breeding purposes (Table 1). ANOVA results showed that the clonal phenotypic average values were significantly different in most of the traits (p ≤ 0.001), except for SSD, for which clone averages did not differ significantly (p = 0.94). Moreover, the growth was significantly different for blocks with certain characteristics, showing the environmental effects in the orchard (Table 1). Average values, ranges, and variation parameters of different traits are shown in Table 2. The average Ht was 11.53 m, with a maximum value of 14.22 m, which was higher than the minimum by over 53%. For diameters at different levels, the highest mean value was observed in BD, which differed from the DBH and DIAM3 by 1.82 cm and 4.36 cm, respectively. Between the maximum average value of volume (0.709 m 3 ) and the minimum, a 40.39% gap was recorded. For stem BT, the maximum average value (9.00 mm) was higher than the minimum one (5.00 mm) by 44%. The SSD showed a gap of less than one unit scale among clones (30% of gap). The maximum average values in crown-related traits were observed in BNN, at over 85%, and the lowest were observed in CB, at 48%. Apart from SSD and CH, the traits with high differences (over 50%) among clonal averages were useful for selecting clones of great genetic value. To delineate the extent of variation of each trait, phenotypic coefficients of variation (PCV) were calculated. The PCV values ranged from 10.73% to 35.45% ( Table 2). The highest PCV (35.45%) was observed in Vol., followed by 29.18% in BNN, and 24.17% in BT. The lowest PCV value (10.76%) was observed in SSD. Characteristics showing greater PCV values, higher than 15%, efficiently supported the selection of the 110 Korean pine clones.
To determine the proportion of the variation that was due to the genotypes of clones, clonal repeatability (R) was computed to separate genotypic and environmental effects on clone growth. The proportion of variation related to the genotypes of clones was thus calculated for each trait. The repeatability values ranged from 0.060 to 0.78 ( Table 2). The highest R (0.78) was found in Vol., followed by diameters, DBH, BD, and DIAM3, with 0.767, 0.766, and 0.724, respectively. Ht (0.641), CH (0.601), and BA (0.523) indicated important genotypic influence on the growth of clones in these traits. The lowest R value was observed in SSD (0.060).

Correlation Analysis
The relationship between traits at the phenotypic and genotypic level was tested by Pearson correlation coefficient. The phenotypic relationships are shown in Table 3. There was strong significant positive correlation between Ht, BD, DBH, DIAM3, CH, and Vol., and between CH and Ht. The highest correlation coefficient was 0.928, between DBH and BD; followed by 0.904 (between the Vol. and DBH and Ht with CH); 0.897 (between the DBH and DIAM3); 0.876 (between Vol. and H); 0.829 (between the BD and DIAM3); 0.801 (between Vol. and DIAM3); 0.708 (between the Vol. and Ht); and 0.671 (between the V and CH). Traits related to the crown (BT and SSD, BA, BNN, CB and CH) were weakly correlated to other traits and to each other. Genotypic correlation coefficients were close to phenotypic correlation coefficients (Figure 2), indicating balanced genotypic and environmental effects on the clones' growth.

Principle Components Analysis
Traits displaying high variation in clones were extracted by principal component analysis (PCA) to describe the correlation between growth characteristics and the intended variation concentrated in the principal components. Four principal components were retained from the PCA analysis with a cumulative variance of 82.11% on the total explained variance ( Table 4). The four principal components comprised 45.36%, 15.11%, 11.23%, and 10.40%, respectively, of the total variation. In the first principal component (F1), the traits Ht, DB, DBH, DIAM3, Vol., and CH exhibited positive eigenvalues, and loading values ranged from 0.37 to 0.44 and from 0.68 to 0.92 (Tables 5 and 6). This shows that the F1 primarily represented stem traits, and that these traits would be useful for selecting stem growth, as there is a high positive correlation between these traits and F1 (Table 7). In the second principal component (F2), the absolute values of Ht, CH, BT, and SSD in the eigenvectors and loading vectors were higher but negative, indicating that F2 represented the stem characteristics that would negatively affect the stem selection. In the third principal component (F3), the eigenvalue and loading values of the SSD, CB, and BNN were higher, and the value of SSD was inferior to that shown in F2, which indicated that F3 represented a part of the crown characteristics. In the fourth PC (F4), BNN and BA were higher that might principally describe the crown characteristics. Traits with high eigenvalues and loading values also showed high correlation in the corresponding components (Table 7). 1.000 ** indicates strong significant correlation at the 0.01 level (2-tailed); * indicates significant correlation at the 0.05 level (2tailed).

Principle Components Analysis
Traits displaying high variation in clones were extracted by principal component analysis (PCA) to describe the correlation between growth characteristics and the intended variation concentrated in the principal components. Four principal components were retained from the PCA analysis with a cumulative variance of 82.11% on the total explained variance ( Table 4). The four principal components comprised 45.36%, 15.11%, 11.23%, and 10.40%, respectively, of the total variation. In the first principal component (F1), the traits Ht, DB, DBH, DIAM3, Vol., and CH exhibited positive eigenvalues, and loading values ranged from 0.37 to 0.44 and from 0.68 to 0.92 (Tables 5 and 6). This shows that the F1 primarily represented stem traits, and that these traits would be useful for selecting stem growth, as there is a high positive correlation between these traits and F1 (Table 7). In the second principal component (F2), the absolute values of Ht, CH, BT, and SSD in the eigenvectors and loading vectors were higher but negative, indicating that F2 represented the stem characteristics that would negatively affect the stem selection. In the third principal component (F3), the eigenvalue and loading values of the SSD, CB, and BNN were higher, and the value of SSD was inferior to that shown in F2, which indicated that F3 represented a part of the crown characteristics. In the fourth PC (F4), BNN and BA were higher that might principally describe the crown characteristics. Traits with high eigenvalues and loading values also showed high correlation in the corresponding components (

Regression Analysis
The best regression model for all 11 morphological traits was selected using the statistical parameters of goodness of fit statistics, including MSE, R 2 , adjusted R 2 , Mallows' Cp, Akaike's AIC, Schwarz's SBC, and Amemiya's PC ( Table 8). The analysis of variance on the dependent variable (Vol.) showed a significant difference (Table 9), indicating that the regression analysis was feasible when taking the stem volume as the dependent variable. Out of 11 independent variables, only BDH and Ht showed high positive regression coefficients on stem volume (0.71 and 0.39, respectively), t-values of 39.93 and 20.64, Pr < 0.0001 (Table 10), and high positive effects on stem volume (Figure 3). The standardized regression coefficients ranged from −0.97 to 1.64. Based on the Type III sum of squares, Ht, DBH, BTH, BA, CB, and SSD showed significant information to explain the variability of stem volume. Among the explanatory variables, DBH was the most influential.

Discussion
Different ramets of cloned genotypes may have various phenotypic performances in different microenvironmental conditions [46], which may lead to a certain amount of variation among trees within genotypes. Variation analysis helped to evaluate the magnitude of difference, and also to see if the difference was significant at the individual level, or at the level of the genotype, the blocks, and/or the interaction of blocks and genotypes [47]. In this study, 11 growth characteristics of stem and crown were analyzed on 38-year-old Korean pine clones. Significant differences were found in variance sources (clone and blocks, p ≤ 0.001) for most of the investigated phenotypic characteristics, indicating the feasible selection of clones. These results were almost similar to previous research on the growth, seed, and wood properties of Korean pine [48], although the level of significance in this study seems slightly high, owing probably to different environmental conditions and different genotypes and ages of the studied materials.
When different clones are planted in a macro-site with approximately similar environmental conditions, the clonal average allows the selection of excellent genotypes by minimizing environmental effects at the individual level [49]. Average performances of different genotypes helped to establish a clonal ranking to select the elite ones [42]. In this study, average values of 14.22 m, 39.40 cm, 37.60 cm, 0.71 m 3 , and 11.40 m in Ht, BD, DBH, DIAM3, Vol., and CH, respectively, and their annual increment (Table 2), were greater than those reported previously for Norway spruce and Casuarina junghuhniana [50,51], reflecting the good growing conditions for Pinus koraiensis at the Naozhi seed orchard.
Variation parameters are mostly considered as the main factors to select trees for breeding purposes [52,53]. In this study, the difference in variation of different traits as given by PCV values varied from 10.73% to 35.45%, which were higher (considering the stand age) in comparison to the study of Lee et al. [54]. This demonstrates that the selection made on this population was well founded, although, compared to those values obtained with Betula platyphylla [43] and Eucalyptus sp. [55], the investigated traits here seem to be less varied. However, these results were close to those reported previously on Korean pine by Liang et al. [11] and Wang et al. [17].
BA, Ht, BD, DBH, DIAM3, CH, and Vol. were highly inherited (0.523 to 0.778) ( Table 2), meaning that these traits will be recognized by different clones than CB, BNN, and SSD, which showed lower repeatability of clones (R, 0.060 to 0.392), even for the unique macro site. The high repeatability in this study was in accordance with Yin et al. [56], which reported repeatability of 0.849, 0.687, and 0.362 in Ht, DBH, and SSD, respectively, in Larix olgensis clones.
The interdependence of traits was revealed by correlation coefficients [57]. This relationship is an important parameter in tree improvement when selecting for multiple characteristics [58]. High genotypic correlation is a high degree of pleiotropy, signifying that many common gene loci affect several traits [24]. In the present study, moderate to high phenotypic and genotypic correlation was found between Ht, BD, DBH, DIAM3, CH, and Vol., indicating that these traits would be improved simultaneously when selecting elite clones using tree volume as an index [59], whereas SSD, CB, or AB would be slightly improved or would need to be selected independently, as they showed weak and negative correlations with stem volume [59]. The correlation results were in agreement with previous studies by Zhao et al. [41] and Hui et al. [60], suggesting indirect selection of Ht, BD, and DBH by volume in Betula platyphylla and Larix olgensis families in Northeastern China.
Although correlation coefficients explained the inference for traits' indirect selection, regression coefficients provided more precision to appreciate the influences of trait interactions and traits' contributions to the related wood yield [61]. Korean pine was firstly bred to improve timber production, and growth traits such as volume were primarily taken into account in the selection process [4]. In this study, tree height, diameter at breast height, and branch number positively influenced the stem volume ( Figure S1). DBH was the most influential trait, followed by Ht. In fact, in artificial plantations, growth in tree height and diameter at breast height were reported to be greatly influenced by interline and tree spacing. Wide spacing, for example, promotes growth in diameter, and reduced spacing promotes growth in height [62,63]. Considering individual influence diagnostics, from three to nine genotypes-including PK105, PK104, PK87, PK86, PK84, PK83, PK74, PK69, PK59, and PK57-showed outlier values, taking observation values of the cloned genotypes as a predictor, stem volume in the studentized deleted residuals with threshold was equal to 3, Cook's distances had a threshold of 0.039, DFFits (Std.) a threshold of 0.603, and DFBetas (Std.) a threshold of 0.191 ( Figure S2). Most of these clones were selected for their high growth performances ( Figure S3). Regarding the goodness of fit statistics parameters, the fitted regression model significantly explained the instructive variables' contribution to the wood yield based on stem volume (Table S1). The findings in this study were in agreement with Wang et al. [63], who found that diameter at breast height and tree height were the traits that most affected the stem volume in Korean pine.
PCA is a valuable technique in the selection process because it reduces a large number of indicators to new principal factors (components) that effectively summarize much of the variation in a complex dataset [64]. The method is very important in the selection of multiple traits, determining suitable traits that contribute most to the variation of the materials to be selected [65]. PCA helped to determine the trait that displayed great contributions to the variations in growth, wood, and seed traits [66]. In this study, variance contributions of Vol., DBH, BD, DIAM3, Ht, and CH were high in the PCI, demonstrating that these characteristics are the key traits accounting for most of the variations in these cloned genotypes, and further, they contribute more to wood yield production in the characterized clone population. The accumulating contribution of the first three principal axes accounted for over 82% of the total variation, which was higher compared to Singh et al. [67] and Du et al. [68]. This finding was likewise high compared to the cumulative variation of 78% found in the study of Wang et al. [17], who used PCA to select suitable traits in growth and wood properties of Korean pine. However, the high cumulative amount of variation in PCA in this study may be due to high multicollinearity among traits (Table S2).
Multiple-trait comprehensive evaluation has served as an accurate practice in multipletrait selection referring to Qi value [41]. Since multiple-trait selection may reduce the genetic gain for some uncorrelated or negatively correlated traits, it is crucial to select appropriate traits to be used as the selection index [69]. Based on correlation coefficients and PCA results, six growth traits-including Ht, BD, DBH, DIAM3, CH, and Vol.-were retained and used to select superior clones in this study. The six elite clones (PK109, PK101, PK14, PK47, PK107, and PK27) displayed high Qi values using a selection rate of 5%. The generated genetic gain reached 42.96% in stem volume (Table 11), which was higher compared to the study of Liang et al. [11], possibly due to high variation (great PCV value) in stem volume in these materials.
Based on growth performance, the economic outcome of the Naozhi orchard was projected. In fact, estimating economic outcome from improved materials is of great importance to know whether outlay-related investments in tree improvement will be covered [70]. The current price reported for P. koraiensis logs was approximately 145.31 USD per cubic meter in Jilin Province [18]. Considering the average volume of clones (0.418 m 3 ), and the genetic gain of 42.96% in volume for elite clones, the price of a single tree will be increased by 26.093 USD (i.e., 15% of the initial price) and the value of the stand volume per hectare will increase by 5325.10 USD with a plant spacing of 7.0 m × 7.0 m, or a density of 204.081 trees/hectare. This achievement was superior to that reported on improved Scots pine in northern Finland after clear-cutting [71]. However, the engendered economic value was significantly higher than that reported on sib-families Korean pine [18], which would be due to the difference in planting age and stand density per hectare, which has been mentioned as having an effect on wood yield per hectare in seed orchards [72].

Conclusions
Determining the selection criterion of trees for a breeding program is of great importance in achieving the goal of tree improvement. Out of 11 investigated traits in this study, six were highlighted based on correlation, regression coefficients, and principal component analysis. High coefficient of variance and heritability were observed, indicating that elite clone selection was viable. Strong positive correlation between traits suggested that those growth traits together promote the clone selection. Using the chosen growth traits as evaluation indicators, six P. koraiensis clones were selected as elite clones using a selection rate of 5%. The estimated economic gain was significantly substantial compared to other tree species. The combination of correlation analysis, regression coefficients, and principal components analysis, in addition to multiple-trait comprehensive evaluation, could provide a basis for highlighting useful growth traits and selecting improved varieties of other tree species. A selection integrating logs, wood properties, and seed production would have to be envisioned in the future in order to achieve a performance that takes into account the full economic value of Korean pine.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/su13084167/s1, Figure S1: Regression analyses on 11 morphological traits of Korean pine grown at the Naozhi orchard. Tree volume was taken as the independent variable and 10 trunk and branch traits as explanatory variables. The analysis was done in XLSTAT software version 2019-1, (https://www.xlstat.com/en/news/xlstat-version-2019-1, accessed on 27 December 2020); Figure  S2: Influence diagnostics of 11 morphological traits of stem and branch in 110 P. koraensis clones at the Naozhi orchard. The stem volume was considered as independent variable. The analysis was done in XLSTAT software version 2019-1 (https://www.xlstat.com/en/news/xlstat-version-2019-1 accessed on 27 December 2020); Figure S3: Selected parent clones of P. koraensis grown in the Naozhi orchard. The figure shows the normal distribution of Qi values on 110 cloned parents. At the right side, the vertical red line highlights the clones in the top 10% of performances. The graph was made using XLSTAT software version 2019-1 (https://www.xlstat.com/en/news/xlstat-version-2019-1, accessed on 27 December 2020); Table S1: Goodness of fit statistics (Vol.); Table S2: Multicollinearity statistics of 11 morphological traits in Korean pine grown at the Naozhi orchard.