Variations in growth traits and wood physicochemical properties among Pinus koraiensis families in Northeast China

Qinhui Zhang·Xiaona Pei·Xianbo Lu·Chunli Zhao·Guangzhi Dong·Wanling Shi·Liankui Wang·Yanlong Li·Xiyang Zhao·Mulualem Tigabu

Abstract This study aimed to explore and improve the different economic values of Pinus koraiensis (Siebold and Zucc.) by examining the variations in 6 growth traits and 9 physicochemical wood properties among 53 P.koraiensis half-sib families.Growth traits assessed included height，diameter at breast height，volume，degree of stem straightness，stem form，and branch number per node，while wood properties assessed included density，fiber length and width，fiber length to width ratio，and cellulose，hemicellulose，holocellulose，lignin，and ash contents.Except for degree of stem straightness and branch number per node，all other traits exhibited highly significant variations (P＜0.01)among families.The coefficients of variation ranged from 5.3 (stem form) to 66.7% (ash content)，whereas，the heritability ranged from 0.136 (degree of stem straightness) to 0.962 (ash content).Significant correlations were observed among growth traits and wood physicochemical properties.Principal component analysis identified four distinct groups representing growth traits，wood chemical and physical properties，and stem form traits.Multi-trait comprehensive evaluation identified three groups of elite families based on breeding objectives，including rapid growth，improved timber production for building and furniture materials，and pulpwood production.These specific families should be used to establish new plantations.

Keywords Growth and wood traits·Pinus koraiensis ·Principal component analysis·Family selection·Economic benefits

Introduction

Pinus koraiensis(Sieb.and Zucc.) is an economically valuable species in Northeast China (Wang et al.2015)，and an important protected plant in the country (Zhang et al.2015a，b).It is generally distributed in the mountains of Xiaoxing’an，Wanda，Zhangguangcai，Laoye，and Changbai(Xia et al.1991).P.koraiensiswood is of high quality (Son et al.2001) with strong corrosion resistance;the seeds have high nutritional value and are marketed as pine nuts (Imbs and Pham 1996).To improve these economic benefits，P.koraiensisshould be established in seed orchards as sources of quality regeneration material for future plantations (Xue et al.2009).

A seed orchard is important to improve forest productivity (Kjaer and Foster 1 996).Numerous coniferous plantations world-wide have originated from seed orchards (Lindgren and Prescher 2005).In some areas，the establishment of third-generation seed orchards has been achieved (Chen and Shen 2005).China began to establish the clonal seed orchard ofCedrus atlantica(Endl.) Manetti ex Carriere in 1964(Wang 2013).Today，over 20 tree species have established in seed orchards (Chen 2001).Similarly，large-scale breeding ofP.koraiensisin China began in the 1980s，and researchers mainly selected materials from naturalP.koraiensisforests(Liu 2017).However，most of these seed orchards are in the primary stages becauseP.koraiensisis characterized by slow growth (Liu et al.2015).

Currently，severalP.koraiensisresearch programs are primarily concerned with the study of chemical components of pine nuts (Kim et al.2012;Shpatov et al.2017;Lee et al.2018)，the identification of molecular markers (Kim et al.2005;Chen et al.2011)，asexual reproduction techniques(Liu 2004;Song et al.2007)，the establishment of seed orchards (Oh et al.2008;Feng et al.2010) and eco-physiological attributes (Makoto and Koike 2007;Sun et al.2016).There has been limited conventional breeding research onP.koraiensisalthough the evaluation of progenies from seed orchards can serve as materials for the selection of elite families and single plants.

In this study，variations in growth traits and wood physicochemical properties amongP.koraiensisfamilies were evaluated.Growth and wood traits are crucial reference standards in tree breeding，particularly as a reference for selecting high performance families and clones for the establishment of seed orchards (Zhang et al.2017).The economic value of wood is determined by its material properties，and this directly affects wood processing and utilization (Li et al.2018).Earlier studies focused on selection and application of growth trait variation (Zhen et al.2007)，but the genetic gain of a single pursuit of growth traits may lower wood quality (Li et al.2014).In addition，reports indicate that growth traits are weakly or negatively correlated with wood density (Ru et al.2001;Su 2001).Therefore，simultaneous selection of growth and wood traits is of considerable significance towards achieving greater economic benefits.P.koraiensiswood is used in building and furniture construction;however，reports onP.koraiensispulpwood are scarce.This paper investigates the economic value ofP.koraiensisto improve them.

A targeted evaluation method can be applied to select families and clones with superior genetic stability according to different breeding objectives (Matziris 1997;Wu et al.2018).Several evaluation methods are available (Luo et al.2010).In this study，the multi-character comprehensive evaluation method was used.Half-sib progenies ofP.koraiensisof the Jilin Provincial Forestry Bureau were used as research materials to:(1) evaluate the variations in growth traits and wood physicochemical properties among 53P.koraiensishalf-sib families，and search for links between growth and wood traits;(2) select high performance families based on different breeding goals;and，(3) provide a new approach for the improvement ofP.koraiensisto increase productivity and improve economic benefits.

Materials and methods

Experimental site

The experimental site was the Sanchazi seed orchard，126°56′ E，41°95′ N in Baishan City，Jilin Province at an altitude of 520 m.The area is characterized by a temperate continental monsoon climate with an average annual temperature of 2.5°C，annual precipitation of 725.5 mm，and annual sunshine of 2，300 h，with a frost-free period of 110 d.

Progeny plantation

For this study，seeds of 53 half-sib families [(family 1，10，49，65，77，82，83，84，85，87，88，90，91，92，93，94，96，97，98，100，101，102，103，105，106，108，109，110，111，112，113，114，115，117，118，120，122，125，129，130，132，133，134，135，136，137，138，140，143，144，145，148，163)，]each with 6 individuals，were acquired from the Sanchazi seed orchard，planted in the spring of 1990 in the nursery and cultivated for 4 years.The progeny plantation was established in 1995 in the spring as a completely randomized block design with 5 blocks and 6 replicates in each plot in rows at 2.0 m×3.0 m spacing.

Growth traits measurements

Growth traits included height (H)，diameter at breast height (DBH)，volume (V)，degree of stem straightness(SSD)，stem fullness (SF)，and branch number per node(BNN)，which were measured at 24 years.SSD was estimated according to Zhao et al.(2014) in a scale of 1-5.A tree with more than two obvious bends in the stem received a score of 1，those with more than two slight bends or with one obviously bend scored 2，those with two slight bends scored 3，those with one slight bend received a score of 4，and finally，those with a completely straight stem scored 5.The SSD values were squared before ANOVA.SF was estimated as the ratio of DBH to basal diameter.Volume was calculated following Zhang et al.(2014) as follows:

where，fdenotes the form factor ofP.koraiensisas 0.33;andD1.3denotes the DBH.

Wood traits measurements

Fifteen wood cores per family (one per tree and three trees per block) were taken from the north to south sides of the stem at 1.3 m (Luo et al.2004)，wrapped in straw and transported to the laboratory for analysis (Sun et al.2018).Wood density (WD)，fiber length (FL)，fiber width (FW)，fiber length to width ratio (LWR)，cellulose (CC)，hemicellulose (HEC)，holocellulose (HOC)，lignin (LC)，and ash (AC) contents were determined.WD was determined as described by Li et al.(2018) using five wood cores for each family.Thirty FLs and FWs were measured on five cores per each family that were mixed (Tang et al.2018)，and LWR was calculated.Dry matter production (DMP)，as the product of WD and V，and cellulose production(CP)，as the product of DMP and CC，were calculated for each family.

The remaining wood cores were ground in a ball mill and mixed to determine chemical properties using an ANKOM A2000i automatic fiber analyzer following Liu(2017).Neutral detergent fiber (NDF) was determined following the digestion of neutral washing solution，and the residues include hemicellulose，cellulose，and lignin.Acid detergent fiber (ADF) was determined following the digestion of the acid washing solution，and the residues include cellulose and lignin.Acid detergent lignin (ADL) treated with 72% sulfuric acid resulted in a lignin residue.Cellulose (CC)，HEC，HOC，and LC were determined using the following relationships:

whereP HEC，P CC，P HOCandP LCdenote percentage hemicellulose (HEC)，cellulose (CC)，holocellulose (HOC)，lignin (LC) contents;andP NDF，P ADF，P ADLandP ACdenote percentage neutral detergent fiber (NDF)，acid detergent fiber (ADF)，acid detergent lignin (ADL) and ash contents(AC).

Statistical analysis

The significant variations in growth traits among families and blocks were tested by ANOVA using the following linear model (Xia et al.2016):

whereXijkdenotes the observation value of thekth individual tree in theith family growing injth block，μthe family mean，αithe family effects，bjthe block effect，αbijthe effect of the interaction between familyiand blockj，eijkrepresent the random error.

Family variations in wood traits were analyzed using the following linear model:

whereYijdenotes the observation value ofjth individual tree inith family，the family mean，Eijthe random error.

Phenotypic coefficient of variation (PCV) was calculated using the following formula (Metougui et al.2017):

Heritability of growth (h 2) and heritability of wood traitswere calculated according to Eqs.9 and 10 (Li 2018):

where denotes variance component for the family effect，the variance component for family by block interaction effects，andthe residual error，b represents the number of blocks，n the number of plots，and r the number of values per family.

The phenotypic correlation between traits was computed according to Han et al.(2017):

whereCovP(x，y) denotes phenotypic covariance between trait x and y，andrepresent the phenotypic variance component for trait x and y，respectively.

The comprehensive evaluation values (Q i)for different families were computed following Zhao et al.(2016):

where ai is Xij/Xjmax;and Xij，Xjmax and n denote the average value of a trait，the maximum value of the trait and the number of traits，respectively.

The genetic gain (ΔG) was estimated following Silva et al.(2008)，and the real gain ( ΔGr) following Ye (2011) as follows:

whereh2is heritability; ΔSand S are the differential of selection;is the family mean.All statistical analysis were performed using Statistical Product and Service Solution(SPSS) software，version 26.0 (IBM Corp.，Armonk，NY，USA).

Results

Variations in growth and wood traits

Highly significant differences (P＜0.01) in growth traits among families and blocks were detected;except for SSD(degree of stem straightness) and BNN (branch number per node) but the interactions were significant for all growth traits(Table 1).Similarly，there were significant differences in wood traits among families (Table 2).There were moderate variations in growth traits，particularly in volume and high variation in ash content (Table 3).The PCV of traits ranged from 5.33% for stem form of 66.67% for ash content (Table 3).The PCV of volume，fiber length，fiber length to width ratio，and hemicellulose content were larger than 25%，whereas the PCV of stem fullness and wood density were less than 10%.Theh 2of growth traits was ＞ 0.7 except for height and degree of stem straightness，in which theh 2of SSD (0.136) was minimal(Table 3).Theh 2of each wood trait was ＞ 0.8.The highesth 2was for lignin content (0.960) and ash content (0.962).

Table 1 ANOVA results for examining variations in growth traits among families and blocks

Correlation analysis

There were highly significant (p＜0.01) positive correlations between H，DBH and V，with strong correlation (r=0.98)between DBH and V (Fig.1).SF (stem fullness) showed a significant (p＜0.01) correlation with H (0.49)，DBH (0.67)，and V (0.65).SSD (degree of stem straightness) and BNN(branch number per node) were significantly (p=0.04) correlated (0.30).For wood traits，fiber length was significantly(p＜0.01) correlated with fiber width (0.42) and fiber length to width ratio (p＜0.01;r=0.45)，but fiber width was significantly (p＜0.05) negatively correlated (-0.61) with fiber length to width ratio.In addition，holocellulose content was highly (p＜0.01) correlated with cellulose content(0.70) and hemicellulose content (0.69).Lignin content was significantly (p＜0.01) correlated with DBH (r=0.43) and volume (r=0.42).

Fig.1 Correlations between different growth and traits of P.koraiensis families (numbers are the correlation coefficients between different traits).Where H，DBH，V，SSD，SF，BNN，WD，F(xiàn)L，F(xiàn)W，LWR，CC，HEC，HOC.LC，and AC stand for height，diameter at breast height，volume，stem straightness degree，stem fullness，branch number per node，wood density，fiber length，fiber width，fiber length to width ratio，cellulose content，hemicellulose content，holocellulose content，lignin content and ash content，respectively.*Correlation is significant at 5% level;and **Correlation is significant at 1% level

Principal components analysis

Principal Component Analysis (PCA) summarizes the information in five significant components (Table 4).The eigenvalue for component I was 3.584 with a contribution of 23.9%;and the eigenvalues of H，DBH，V，and SF were higher，and therefore component 1 mainly represented the growth traits.The eigenvalue for component II was 2.263，with a contribution of 15.1%;and the eigenvalue of CC，HEC，and HOC were higher，and therefore this component represented the wood chemical properties.The eigenvalue for component III was 1.935，with a contribution of 12.9%;and the absolute eigenvalues of WD，F(xiàn)W，and LWR were higher，thus this component represented the physical properties of wood.The eigenvalue for component IV was 1.556，with a contribution of 10.4%;the eigenvalues of SSD and BNN were 0.643 and 0.659，respectively.Therefore，this component represented the stem form traits.The cumulative contribution of the five principal components was 71.1%，which provided most of the information on the growth and wood traits of the families.

Table 2 ANOVA results for examining variations in wood traits among families

Table 3 Descriptive statistics for different growth and wood traits together with phenotypic coefficient of variation (PCV)and heritability (h2 )

Family selection

To improve the growth ofP.koraiensis，H，DBH，V，and SF were applied as evaluation indices of growth traits through principal components analysis.The results of multiple-trait comprehensive evaluation analysis are shown in Table 5.At a 10% selection intensity，five families (115，117，120，133 and 138) were identified as families with rapid growth.The selected families hadQivalues of 1.997，1.915，1.863，1.870，and 1.917，and with average H，DBH，V，and SF of 0.48 m，2.66 cm，0.0304 m3，and 0.02，respectively whereas their genetic gains were 2.8%，12.1%，29.4%，and 1.9%，respectively.

High-quality materials for building and furniture construction could be obtained through the selection of families based on H，DBH，V，SF of growth traits and WD，F(xiàn)W，and LWR of physical properties of wood as evaluation indices.The results of the multiple-trait comprehensive evaluation analysis are shown in Table 6.At a 10% selection intensity，five families (112，115，117，135，and 138) were identified as elite families for production of high-quality timber.The selected families hadQivalues 2.160，2.145，2.219，2.140，and 2.126.Average H，DBH，V，SF，WD，F(xiàn)W and LWR were 8.0 m，18.9 cm，0.105 m3，0.8，488.4 kg m-3，28.4 μm and 66.6，whereas their genetic gains were 2.4%，11.4%，26.9%，1.8%，3.6%，-9.5% and 15.9%，respectively.

Families with superior pulp properties were selected via the multi-trait comprehensive evaluation of fiber length(FL)，length/width ratio (LWR) and cellulose production per plant (CP) as evaluation indices.TheQivalues are presentedin Table 7.At 10% selection intensity，five elite families(1，112，115，117，and 130) were chosen for the production of high-quality pulp.The selected families hadQivalues 1.603，1.593，1.637，1.619，and 1.582.Average FL，LWR，and CP were 1934.02 μm，70.52 and 16.01 kg，and their real gains were 9.1%，23.9%，and 43.2%，respectively.

Table 4 Principal component analysis of different traits resulting in four distinct groupings that represented growth traits (Component I)，chemical properties of wood (Component II)，physical properties of wood (Component III) and stem form traits (Component IV)

Table 5 Qi value of different families that are suitable for improving growth

Table 6 Qi value of different families for improving timber for building material and furniture

Table 7 Qi value of different families for improving pulpwood materials

Discussion

Variability in growth traits and wood physico-chemical properties

Genetic variation drives evolution，which in turn enables the adaptation and survival of organisms under different environmental conditions (Cardoso et al.2015).Therefore，to understand the genetic variation of forest growth，the growth and wood properties of 28-year-old half-sib families ofP.koraiensiswere evaluated.All wood traits and most of the growth traits were significantly different among families，which concurs with the findings of Yin(2017).Thus，there is room for improvement through the selection of elite families.

Assessment of inter-and intra-population variation is a prerequisite for developing effective forest genetic resource conservation strategies (Papi et al.2012).In this research，a lower average of growth and wood traits was reported than those reported by Wang et al.(2018a)，which may be ascribed to different tree ages.In addition，the HOC and LC of wood characteristics were lower，46.8% and 23.7%，respectively，which may be explained by the fact that the wood samples were too finely ground.Consequently，the exudation of the sample from the filter bag led to low experimental results.

The phenotypic coefficient of variation (PCV) is a widely adopted scaleless measure of variability (Wang et al.2018b).It reflects the potential genetic variation among various traits in each population.Our study revealed that the PCV of H，DBH，and V was 10.8%，16.3%，and 37.4%，respectively，lower than that ofEucalyptus urophyllaS.T.Blake (Lu et al.2010).The PCV of FL，F(xiàn)W，and LWR were 30.0%，23.5%，and 26.4%，respectively，higher than that ofEucalyptus camaldulensisDehnh.(Shang et al.2017.The PCV of CC，HEC and LC was 17.1%，11.1%，and 13.9%，respectively and higher than those ofTrachycarpus fortunei(Hook.) H.Wendl.reported by Hu and Zhao (2018).Results of the present study demonstrated abundant variations in growth and wood traits，which necessitates genetic improvement in different directions according to different breeding objectives.

Similarly，heritability is the measure of the degree to which parents pass on genetic characteristics to offspring(Abengmeneng et al.2015).The heritability of each wood trait in this study was over 0.8，significantly higher than that forPinus massonianaLamb.reported by Chen et al.(2015).These findings lay a foundation for the selection of families for pulpwood and building materials.Moreover，the heritability of H，DBH，and V in our study was 0.439，0.746，and 0.745，respectively.These values are higher than the results of Chen et al.(2016);higher heritability is beneficial for parents to pass on good characteristics to their offspring，making the selection of families more meaningful.

Association between growth traits and wood properties

There were highly significant positive correlations between height，diameter at breast height，and volume，which concur with the results of Liang et al.(2018) for different clones ofP.koraiensis.In most cases，height，diameter，and volume production are used to evaluate genetic variation in tree growth.Because these traits are relatively easy to assess，they have become a common indicator for the evaluation of forestry production (Rweyongeza et al.2005).In our study，there were strong correlations between height，diameter at breast height，and volume，which are beneficial to the evaluation and selection of fast-growingP.koraiensisfamilies.Among wood physical properties，fiber length was significantly positively correlated with fiber width (0.422)and fiber length/width ratio (0.452).Similar results were reported by Wang (2009) who studiedP.massoniana.The findings suggest the possibility of synchronized improvement of fiber length and width through selection of one of these traits.There was a highly significant correlation between holocellulose and cellulose (0.700) and hemicellulose contents (0.693)，but lignin content was not strongly correlated with cellulose and hemicellulose contents，which is different from the results of Yin et al.(2017).This may be attributed to the fine grinding of the sample，whereby the sample leaks out of the filter bag and affects the calculation of cellulose，hemicellulose，and lignin contents.On the other hand，lignin was strongly correlated with diameter at breast height (0.432) and volume (0.420);this is different and contrary to findings from previous studies.For instance，Zhang et al.(2015a，b) reported the lack of a significant correlation between lignin content and growth traits ofPopulus tomentosaCarriére，whereas Yu (2015) reported a significant negative correlation between lignin content and growth traits ofLarix gmeliniiLedebe ex Gordon.Therefore，the correlation of growth traits with wood properties is possibly related to different species and growing environments.The correlation analysis aimed to find the best linear combination between two sets of multivariate data that maximizes the correlation coefficient between them (Malacarne 2014).As a common statistical tool，correlation analysis can intuitively reflect the correlation between the datasets (Jendoubi and Strimmer 2019).

Principal components analysis (PCA) is a classic multivariate technique that reduces the dimensional complexity of a set of variables to a new set of unrelated variables，known as the principal component，which retains all the variable information (Fernández et al.2 015).It is applicable in the exploration of the relationship between different tree characters (Hoeber et al.2014).In this study，all traits were divided into five principal components with a cumulative contribution rate of 71.1%，lower than the results of a study by López de Heredia et al.(2009) on oaks.Following the analysis，principal component factors described growth traits，wood chemical properties，wood physical properties，and stem form.Each principal component objectively reflected the relationship among the analyzed traits and provided a basis for the evaluation and selection of elite families.

Selection of elite families

To select eliteP.koraiensisfamilies with different attributes，the families were screened according to three breeding objectives:fast-growth，production of high-quality construction timber，and production of pulpwood.

Fast growing families

Regarding best growth in height，diameter，volume，and stem fullness，five elite families were selected that had better growth for timber production and urban planting.These elite families showed 2.8%，12.1%，29.4%，and 1.9%genetic gain in height，diameter，volume，and stem fullness，respectively.The results were higher than that reported forP.koraiensisby Wang et al.(2018a)，and similar to a study on Douglas-fir (Clair et al.2004).The growth performance of these selected elite families was good;thus they can be used as source of planting material for the establishment of fast-growingP.koraiensisplantations.Elite families with higher genetic gains are beneficial to promote and apply in forestry production.

Families suitable for production of high quality building materials

To obtain high-quality materials for building and furniture，height，diameter，volume，stem fullness，density，f bier width，and fiber length/width ratio were used as indicators for comprehensive evaluation based on the results of principal component analysis.Principal component analysis found wood density and fiber length to width ratio as representing the physical properties of wood.Wood density is an important feature of trees that is central to many functional processes，and is considered a trait that powerfully integrates various wood characteristics (Woodall et al.2015).It directly affects wood quality.Fiber width and the ratio between fiber length and width are common wood fiber characteristics.Some studies have shown that fiber width is related to the wood structure，which is applied to judge the thickness of the wood(Huang and Chen 2014).A fine wood structure implies that the material is more compact，the surface is smoother，and is easier to cut and machine.Following the screening of 53P.koraiensisfamilies，5 elite families with these characteristics were identified.The genetic gains of each trait were 2.4%，11.4%，26.9%，1.8%，3.6%，-9.5%，and 15.9% for height，diameter，volume，stem fullness，density，f bier width，and fiber length/width ratio，respectively，higher than those reported by Wang et al.(2018a).The wood density of elite families was relatively high at 488.4 kg m-3.The selected elite families can be utilized for structural and decorative materials，including floor platforms and garden furniture.

Families suitable for pulpwood production

There are limited studies on the pulpwood production ofP.koraiensis.However，given its long fiber length(0.90-3.00 mm)，P.koraiensisis highly suitable to produce high-strength paper (Wang 2000).Therefore，the utilization value ofP.koraiensisfor pulpwood was evaluated.Fiber length，length/width ratio，and cellulose production were evaluation indices for pulpwood.Notably，the nature of wood determines its economic value and directly affects wood processing and utilization.Fiber length is the most important morphological parameter in papermaking.Wang et al.(2016) showed that longer fibers could achieve higher tear strength，tensile strength，and flexural strength.The length/width ratio of fiber is an important morphological index second only to fiber length.A larger length/width ratio is related to a larger bonding area of the fiber，a higher paper index，and higher paper strength (Chang et al.2007).In this study，5 elite families were identified using fiber length，length/width ratio，and cellulose production as evaluation indexes.The real gains were 9.1%，23.9%，and 43.2%，respectively.The genetic improvement is remarkable and can provide valuable material with fine quality fiber.

Conclusion

We examined variations in growth and wood properties of 53 half-sib families ofP.koraiensis.Substantial genetic variations in growth traits and wood physicochemical properties were found among families with good heritability，which in turn enabled selection of elite families.Consequently，five elite families that had better growth for timber production and urban planting;five families suitable for production of high-quality building materials;and five elite families suitable for pulpwood production were selected，with good genetic gain.Thus，these families should be used in future plantation establishment to achieve economic objectives.

Author contributionsQHZ，CLZ，MT and XYZ conceived and designed the project，and QHZ，XNP，XBL and CLZ performed the experiments and wrote the manuscript.GZD，WLS，LKW and YLL assisted in conducting the experiments and analyzing the data.MT revised and edited the manuscript.All the authors read and approved the final manuscript.Open AccessThis 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/.